Identifying Genes Whose Mutant Transcripts Cause Dominant Disease Traits by Potential Gain-of-Function Alleles

Genotype-phenotype relationship in RDH12 retinopathy: a perspective from a pediatric age group

INTRODUCTION

Monocentric retrospective case series to describe clinical and molecular peculiarities in a series of pediatric patients in attempting a possible genotype-phenotype correlation.

METHODS

We included 13 pediatric patients from 7 unrelated families (ages 1-18) with biallelic pathogenic and likely pathogenic variants in RDH12 gene. For all our patients segregation analyses were performed in their parents and affected siblings. According to their cooperation, patients underwent a complete ophtalmic examination and imaging with full field standard electroretinography (ffERG), spectral domain optical coherence tomography (SD-OCT) and fundus autofluorescence (FAF).

RESULTS

According to previous studies, we did not observe a conclusive genotype-phenotype correlation in our series, nevertheless we reported 2 new RDH12 likely pathogenic variants. Also, we report clinical data on pediatric patients, with the fundus imaging in the youngest child (2 yo) described in the literature in whom retinal dystrophic changes are already present.

DISCUSSION

This study includes a collection of genotypic and phenotypic data from children with RDH12-associated IRD. These findings will help further characterize RDH12-related retinopathy. Determining the time window of onset of dystrophic changes is critical to research the correct timing for administering possible therapies. More extensive and functional studies are needed in view of the opportunity of gene replacement therapy for RDH12 associated IRD.

A series of reviews in familial cancer: genetic cancer risk in context variants of uncertain significance in MMR genes: which procedures should be followed?

Interpreting variants of uncertain significance (VUS) in mismatch repair (MMR) genes remains a major challenge in managing Lynch syndrome and other hereditary cancer syndromes. This review outlines recommended VUS classification procedures, encompassing foundational and specialized methodologies tailored for MMR genes by expert organizations, including InSiGHT and ClinGen's Hereditary Colorectal Cancer/Polyposis Variant Curation Expert Panel (VCEP). Key approaches include: (1) functional data, encompassing direct assays measuring MMR proficiency such as in vitro MMR assays, deep mutational scanning, and MMR cell-based assays, as well as techniques like methylation-tolerant assays, proteomic-based approaches, and RNA sequencing, all of which provide critical functional evidence supporting variant pathogenicity; (2) computational data/tools, including in silico meta-predictors and models, which contribute to robust VUS classification when integrated with experimental evidence; and (3) enhanced variant detection to identify the actual causal variant through whole-genome sequencing and long-read sequencing to detect pathogenic variants missed by traditional methods. These strategies improve diagnostic precision, support clinical decision-making for Lynch syndrome, and establish a flexible framework that can be applied to other OMIM-listed genes.

De novo variants in RYBP are associated with a severe neurodevelopmental disorder and congenital anomalies

PURPOSE

Polycomb group proteins are key epigenetic transcriptional regulators. Multiple neurodevelopmental disorders are associated with pathogenic variants of the genes encoding Polycomb group proteins. RYBP is a core component of the noncanonical Polycomb Repressor Complex 1; however, its role in disease is unclear.

METHODS

Functional consequences of RYBP variants were assessed using in vitro cellular and in vivo Drosophila melanogaster studies.

RESULTS

We described 7 individuals with heterozygous de novo variants of RYBP and their clinical findings, including severe developmental delay, dysmorphisms, and multiple congenital anomalies. We showed that all single-nucleotide variants in RYBP localize to the N-terminal domain of the gene, which encodes the zinc-finger domain and ubiquitin-binding moiety. In vitro studies have demonstrated that the RYBP c.132C>G p.(Cys44Trp) variant causes reduced protein expression but does not affect the binding of YY1, RING1B, or ubiquitin. In vivo overexpression studies in Drosophila melanogaster showed a dramatic functional difference between human RYBP and its variant forms, affecting the C44 amino acid residue. DNA methylation studies suggested a possible episignature associated with RYBP-related disorder.

CONCLUSION

Heterozygous de novo variants in RYBP are associated with an identifiable syndromic neurodevelopmental disorder with multiple congenital anomalies.

Heterozygous variants in AP4S1 are not associated with a neurological phenotype

Biallelic loss-of-function variants in AP4S1 cause childhood-onset hereditary spastic paraplegia. A recent report suggested that heterozygous AP4S1 variants lead to a syndrome of lower limb spasticity and dysregulation of sphincter function. We critically evaluate this claim against clinical observations in 28 heterozygous carriers of the same AP4S1 variant (NM_007077.3: c.289C>T, p.Arg97Ter). In these 14 males and 14 females (mean age: 37.6 ± 4.9 years [SD], range: 30-50 years), we ascertain no increased prevalence of neurological manifestations. Alternative causes should be considered when evaluating patients with heterozygous AP4S1 variants and neurological symptoms, as misattribution of pathogenicity can impact clinical care and genetic counseling.

Biallelic variation in the choline and ethanolamine transporter FLVCR1 underlies a severe developmental disorder spectrum

PURPOSE

FLVCR1 encodes a solute carrier protein implicated in heme, choline, and ethanolamine transport. Although Flvcr1-/- mice exhibit skeletal malformations and defective erythropoiesis reminiscent of Diamond-Blackfan anemia (DBA), biallelic FLVCR1 variants in humans have previously only been linked to childhood or adult-onset ataxia, sensory neuropathy, and retinitis pigmentosa.

METHODS

We identified individuals with undiagnosed neurodevelopmental disorders and biallelic FLVCR1 variants through international data sharing and characterized the functional consequences of their FLVCR1 variants.

RESULTS

We ascertained 30 patients from 23 unrelated families with biallelic FLVCR1 variants and characterized a novel FLVCR1-related phenotype: severe developmental disorders with profound developmental delay, microcephaly (z-score -2.5 to -10.5), brain malformations, epilepsy, spasticity, and premature death. Brain malformations ranged from mild brain volume reduction to hydranencephaly. Severely affected patients share traits, including macrocytic anemia and skeletal malformations, with Flvcr1-/- mice and DBA. FLVCR1 variants significantly reduce choline and ethanolamine transport and/or disrupt mRNA splicing.

CONCLUSION

These data demonstrate a broad FLVCR1-related phenotypic spectrum ranging from severe multiorgan developmental disorders resembling DBA to adult-onset neurodegeneration. Our study expands our understanding of Mendelian choline and ethanolamine disorders and illustrates the importance of anticipating a wide phenotypic spectrum for known disease genes and incorporating model organism data into genome analysis to maximize genetic testing yield.

Novel Splice-Altering Variants in the CHM and CACNA1F Genes Causative of X-Linked Choroideremia and Cone Dystrophy

BACKGROUND

An estimated 10-15% of all genetic diseases are attributable to variants in noncanonical splice sites, auxiliary splice sites and deep-intronic variants. Most of these unstudied variants are classified as variants of uncertain significance (VUS), which are not clinically actionable. This study investigated two novel splice-altering variants, CHM NM_000390.4:c.941-11T>G and CACNA1F NM_005183.4:c.2576+4_2576+5del implicated in choroideremia and cone dystrophy (COD), respectively, resulting in significant visual loss.

METHODS

Next-generation sequencing was employed to identify the candidate variants in CHM and CACNA1F, which were confirmed using Sanger sequencing. Cascade analysis was undertaken when additional family members were available. Functional analysis was conducted by cloning genomic regions of interest into gateway expression vectors, creating variant and wildtype midigenes, which were subsequently transfected into HEK293 cells. RNA was harvested and amplified by RT-PCR to investigate the splicing profile for each variant compared to the wildtype. Novel variants were reclassified according to ACMG/AMP and ClinGen SVI guidelines.

RESULTS

Midigene functional analysis confirmed that both variants disrupted splicing. The CHM NM_000390.4:c.941-11T>G variant caused exon 8 skipping, leading to a frameshift and the CACNA1F NM_005183.4:c.2576+4_2576+5del variant caused a multimodal splice defect leading to an in-frame insertion of seven amino acids and a frameshift. With this evidence, the former was upgraded to likely pathogenic and the latter to a hot VUS.

CONCLUSIONS

This study adds to the mutational spectrum of splicing defects implicated in retinal degenerations by identifying and characterising two novel variants in CHM and CACNA1F. Our results highlight the importance of conducting functional analysis to investigate the consequences of intronic splice-altering variants and the significance of reclassifying VUS to confirm a genetic diagnosis.

GREGoR: Accelerating Genomics for Rare Diseases

Rare diseases are collectively common, affecting approximately one in twenty individuals worldwide. In recent years, rapid progress has been made in rare disease diagnostics due to advances in DNA sequencing, development of new computational and experimental approaches to prioritize genes and genetic variants, and increased global exchange of clinical and genetic data. However, more than half of individuals suspected to have a rare disease lack a genetic diagnosis. The Genomics Research to Elucidate the Genetics of Rare Diseases (GREGoR) Consortium was initiated to study thousands of challenging rare disease cases and families and apply, standardize, and evaluate emerging genomics technologies and analytics to accelerate their adoption in clinical practice. Further, all data generated, currently representing ~7500 individuals from ~3000 families, is rapidly made available to researchers worldwide via the Genomic Data Science Analysis, Visualization, and Informatics Lab-space (AnVIL) to catalyze global efforts to develop approaches for genetic diagnoses in rare diseases (https://gregorconsortium.org/data). The majority of these families have undergone prior clinical genetic testing but remained unsolved, with most being exome-negative. Here, we describe the collaborative research framework, datasets, and discoveries comprising GREGoR that will provide foundational resources and substrates for the future of rare disease genomics.

Clinical and molecular spectrum along with genotype-phenotype correlation of 25 patients diagnosed with 3 M syndrome: a study from Turkey

3 M syndrome is a well-known autosomal recessive skeletal genetic disorder caused by biallelic pathogenic variants in the CUL7, OBSL1, and CCDC8 genes. Affected individuals exhibit profound pre- and postnatal growth retardation, distinctive facial features with normal intelligence. This study aims to provide insight into the comprehensive evaluation of clinical, laboratory, and radiological findings, expand the mutational spectrum of the disease, and establish a genotype-phenotype correlation in the present cases. A total of 25 patients from 19 unrelated families were included in the study. Genetic etiology was determined in probands through the utilization of Sanger sequencing and/or targeted gene panel analysis. The clinical, laboratory, and genetic features of all patients at admission and during follow-up were documented. Genotype-phenotype correlation was carried out in the CUL7 and OBSL1 groups. The genetic etiology was established in all patients (n = 25/25, 100%). We identified 15 distinct variants in CUL7, OBSL1, and CCDC8 genes, with eleven being novel. CUL7 variants were present in 13 patients (n = 13/25, 52%), while OBSL1 variants were found in 11 patients (n = 11/25, 44%). No notable distinctions were found in mean birth weight, height, and standard deviation scores between the CUL7 and OBSL1 mutation groups (p > 0.05). Patients with CUL7 variants exhibited significantly lower height standard deviation scores both at admission and at the last examination, as well as lower weight standard deviation scores at the last examination, compared to those with OBSL1 variants (p < 0.05).

CONCLUSION

To date, genotype-phenotype correlations have been identified in a limited number of studies. Further research involving larger cohorts is necessary to solidify these correlations.

WHAT IS KNOWN

• 3M syndrome is a well-known skeletal dysplasia caused by biallelic pathogenic variants in CUL7, OBSL1, and CCDC8 genes. • Despite genetic heterogeneity, clinical, and radiologic features show homogeneity in affected individuals.

WHAT IS NEW

• Genotype-phenotype correlations have been established in limited studies. • The CUL7 group exhibited significantly lower height SDS at both admission and the final evaluation and lower weight SDS at the final examination compared to the OBSL1 group. • The frequency of variants in the OBSL1 gene among Turkish patients exceeds the rates reported in the literature. • Gradenigo syndrome is being reported for the first time in a patient with 3M syndrome.

Using multiplexed functional data to reduce variant classification inequities in underrepresented populations

BACKGROUND

Multiplexed Assays of Variant Effects (MAVEs) can test all possible single variants in a gene of interest. The resulting saturation-style functional data may help resolve variant classification disparities between populations, especially for Variants of Uncertain Significance (VUS).

METHODS

We analyzed clinical significance classifications in 213,663 individuals of European-like genetic ancestry versus 206,975 individuals of non-European-like genetic ancestry from All of Us and the Genome Aggregation Database. Then, we incorporated clinically calibrated MAVE data into the Clinical Genome Resource's Variant Curation Expert Panel rules to automate VUS reclassification for BRCA1, TP53, and PTEN.

RESULTS

Using two orthogonal statistical approaches, we show a higher prevalence (p ≤ 5.95e - 06) of VUS in individuals of non-European-like genetic ancestry across all medical specialties assessed in all three databases. Further, in the non-European-like genetic ancestry group, higher rates of Benign or Likely Benign and variants with no clinical designation (p ≤ 2.5e - 05) were found across many medical specialties, whereas Pathogenic or Likely Pathogenic assignments were increased in individuals of European-like genetic ancestry (p ≤ 2.5e - 05). Using MAVE data, we reclassified VUS in individuals of non-European-like genetic ancestry at a significantly higher rate in comparison to reclassified VUS from European-like genetic ancestry (p = 9.1e - 03) effectively compensating for the VUS disparity. Further, essential code analysis showed equitable impact of MAVE evidence codes but inequitable impact of allele frequency (p = 7.47e - 06) and computational predictor (p = 6.92e - 05) evidence codes for individuals of non-European-like genetic ancestry.

CONCLUSIONS

Generation of saturation-style MAVE data should be a priority to reduce VUS disparities and produce equitable training data for future computational predictors.

Clinical and Molecular Findings in Patients with Knobloch Syndrome 1: Case Series Report

BACKGROUND/OBJECTIVES

Knobloch syndrome 1 (KS) is an autosomal recessive inherited ocular syndrome characterized by a combination of high myopia, vitreoretinal degeneration, and occipital encephalocele. KS is caused by biallelic pathogenic variants in the COL18A1 gene. Diagnosing KS can be challenging due to its clinical heterogeneity and the rarity of the syndrome.

METHODS

We conducted comprehensive clinical and instrumental ophthalmological examinations, whole-exome sequencing, Sanger sequencing, and segregation analysis to evaluate affected families.

RESULTS

Two patients presenting with high myopia, low visual acuity, chorioretinal atrophy, and occipital skin/skull defects were diagnosed with Knobloch syndrome 1 (KS). In Case 1, a 14-year-old boy, the COL18A1 variants identified were c.2673dup and c.3523_3524del in a compound heterozygous state. Case 2 involved a 3-year-old girl, the c.1637_1638dup and c.3523_3524del variants were identified in a compound heterozygous state. In Case 3, a retrospectively observed boy of 3 y.o. with KS, the variants c.929-2A>G and c.3523_3524del were defined earlier.

CONCLUSIONS

We confirmed KS molecularly in two novel families. Additionally, in Case 3 of a retrospectively analyzed third family and in both novel cases, one of the biallelic causative variants was the same known 2bp deletion in exon 40 of the collagen XVIII gene. Cases 1 and 3 were characterized by connective tissue dysplasia features and a pathognomonic Knobloch triad. No neurological manifestations and no trends in the genotype-phenotype relationship were found. The heterogeneity of phenotype in the case series is likely to be the result of further factors and/or genetic background.

Considerations for reporting variants in novel candidate genes identified during clinical genomic testing

Since the first novel gene discovery for a Mendelian condition was made via exome sequencing, the rapid increase in the number of genes known to underlie Mendelian conditions coupled with the adoption of exome (and more recently, genome) sequencing by diagnostic testing labs has changed the landscape of genomic testing for rare diseases. Specifically, many individuals suspected to have a Mendelian condition are now routinely offered clinical ES. This commonly results in a precise genetic diagnosis but frequently overlooks the identification of novel candidate genes. Such candidates are also less likely to be identified in the absence of large-scale gene discovery research programs. Accordingly, clinical laboratories have both the opportunity, and some might argue a responsibility, to contribute to novel gene discovery, which should, in turn, increase the diagnostic yield for many conditions. However, clinical diagnostic laboratories must necessarily balance priorities for throughput, turnaround time, cost efficiency, clinician preferences, and regulatory constraints and often do not have the infrastructure or resources to effectively participate in either clinical translational or basic genome science research efforts. For these and other reasons, many laboratories have historically refrained from broadly sharing potentially pathogenic variants in novel genes via networks such as Matchmaker Exchange, much less reporting such results to ordering providers. Efforts to report such results are further complicated by a lack of guidelines for clinical reporting and interpretation of variants in novel candidate genes. Nevertheless, there are myriad benefits for many stakeholders, including patients/families, clinicians, and researchers, if clinical laboratories systematically and routinely identify, share, and report novel candidate genes. To facilitate this change in practice, we developed criteria for triaging, sharing, and reporting novel candidate genes that are most likely to be promptly validated as underlying a Mendelian condition and translated to use in clinical settings.

Preclinical alternative drug discovery programs for monogenic rare diseases. Should small molecules or gene therapy be used? The case of hereditary spastic paraplegias

Patients diagnosed with rare diseases and their and families search desperately to organize drug discovery campaigns. Alternative models that differ from default paradigms offer real opportunities. There are, however, no clear guidelines for the development of such models, which reduces success rates and raises costs. We address the main challenges in making the discovery of new preclinical treatments more accessible, using rare hereditary paraplegia as a paradigmatic case. First, we discuss the necessary expertise, and the patients' clinical and genetic data. Then, we revisit gene therapy, de novo drug development, and drug repurposing, discussing their applicability. Moreover, we explore a pool of recommended in silico tools for pathogenic variant and protein structure prediction, virtual screening, and experimental validation methods, discussing their strengths and weaknesses. Finally, we focus on successful case applications.

Targeted resequencing reveals high-level mosaicism for a novel frameshift variant in WDR45 associated with beta-propeller protein-associated neurodegeneration

OBJECTIVES

Beta-propeller protein-associated neurodegeneration (BPAN) is a rare X-linked dominant neurodegenerative disease, which is characterized by iron accumulation in the basal ganglia. BPAN is associated with pathogenic variation in WDR45, which has been reported almost exclusively in females most probably due to male lethality in the hemizygous state.

METHODS

Whole exome sequencing (WES) and targeted deep sequencing were performed for a male with a clinical diagnosis of BPAN at the age of 37.

RESULTS

The novel frameshift variant in WDR45 detected by WES was further analyzed with targeted resequencing to detect a mosaic variant with a level of 85.5% in the blood sample of the proband.

DISCUSSION

Although the main role of WDR45 remains elusive, recent studies show that WDR45 may contribute to neurodegeneration through defects in autophagy, iron storage and ferritin metabolism, mitochondria organization, and endoplasmic reticulum homeostasis. The extend of spatiotemporal haploinsufficiency of WDR45 frameshifting variants caused by mosaicism in males may lead to variable clinical severity, which may be hard to elaborate clinically. Promising genetic analysis strategies using targeted deep sequencing may help determine the clinical outcome of somatic mosaicism in neurological disorders including BPAN. Additionally, we suggest that deep sequencing should be conducted in cerebrospinal fluid samples to provide more reliable results in terms of reflecting the mosaicism level in the brain for future studies.

Biallelic TYR and TKFC variants in Egyptian patients with OCA1 and new expanded TKFC features

BACKGROUND

Oculocutaneous albinism type1 (OCA1) is caused by the TYR gene's homozygous and compound heterozygous variants. TKFC gene variants cause triokinase & FMN cyclase deficiency syndrome with variable multisystemic disorders.

OBJECTIVES

To determine the potential disease-causing variants in two deceased patients presenting atypical OCA1 features by demonstrating three generations for a single family. The two deceased neonates had severe skeletal abnormalities and fatal hypertrophic cardiomyopathy. We also explored the potential mechanisms for the causative relationship between TKFC and multisystem disorders.

PATIENTS AND METHODS

Due to the new emerging symptoms that weren't reported before with the TYR gene, the following methods were performed: Sanger sequencing for the TYR gene, followed by whole exome sequencing, co-segregation, and computational analyses.

RESULTS

Extensive parental consanguinity was found, and consequently an autosomal recessive mode of inheritance was prioritized. Upon performing sequencing and segregation data, the following has been confirmed: positive co-segregation of nonsense homozygous NM_000372.5:c.346C > T p.(Arg116*) variant in TYR gene and multisystem disease-missense homozygous NM_015533.4:c.598G > A p.(Val200Ile) variant in TKFC gene in the two affected index patients who deceased due to hypertrophic cardiomyopathy. Using computational analysis, we found that c.598G > A p.(Val200Ile) pathogenicity has led to the failure of L2-K1 active site closure due to the potential differential fluctuation between valine and isoleucine residues. Subsequently, disruption of endogenous DHA phosphorylation was found. Two potential mechanisms exploring the causative relationship between TKFC gene and multisystem disorders have been suggested.

CONCLUSIONS

This study presented a first family with the co-existence of biallelic variants in TYR and TKFC genes associating severe skeletal abnormalities and lethal hypertrophic cardiomyopathy. Neither of these genes would have been pursued in the standard genetic counseling. Such discovery is paving the way for more efficient genetic counseling. Comparing TKFC results with literature data showed that our relevant expanded TKFC variant is the 3rd worldwide.

Maple syrup urine disease diagnosis in Brazilian patients by massive parallel sequencing

Biallelic pathogenic variants cause maple syrup urine disease (MSUD) in one of the branched-chain α-keto acid dehydrogenase (BCKDH) complex genes (BCKDHA, BCKDHB, DBT, DLD, and PPM1K) leading to the accumulation of leucine, isoleucine, and valine. This study aimed to perform a molecular diagnosis of Brazilian patients with MSUD using gene panels and massive parallel sequencing. Eighteen Brazilian patients with a biochemical diagnosis of MSUD were analyzed by massive parallel sequencing in the Ion PGM Torrent Server using a gene panel with the BCKDHA, BCKDHB, and DBT genes. The American College of Medical Genetics and Genomics guidelines were used to determine variant pathogenicity. Thirteen patients had both variants found by massive parallel sequencing, whereas 3 patients had only one variant found. In 2 patients, the variants were not found by this analysis. These 5 patients required additional Sanger sequencing to confirm their genotype. Twenty-five pathogenic variants were identified in the 3 MSUD-related genes (BCKDHA, BCKDHB, and DBT). Most variants were present in the BCKDHB gene, and no common variants were found. Nine novel variants were observed: c.922 A > G, c.964C > A, and c.1237 T > C in the BCKDHA gene; and c.80_90dup, c.384delA, c.478 A > T, c.528C > G, c.977 T > C, and c.1039-2 A > G in the BCKDHB gene. All novel variants were classified as pathogenic. Molecular modeling of the novel variants indicated that the binding of monomers was affected in the BCKDH complex tetramer, which could lead to a change in the stability and activity of the enzyme. Massive parallel sequencing with targeted gene panels seems to be a cost-effective method that can provide a molecular diagnosis of MSUD.

Glycogen myophosphorylase loss causes increased dependence on glucose in iPSC-derived retinal pigment epithelium

Loss of glycogen myophosphorylase (PYGM) expression results in an inability to break down muscle glycogen, leading to McArdle disease-an autosomal recessive metabolic disorder characterized by exercise intolerance and muscle cramps. While previously considered relatively benign, this condition has recently been associated with pattern dystrophy in the retina, accompanied by variable sight impairment, secondary to retinal pigment epithelial (RPE) cell involvement. However, the pathomechanism of this condition remains unclear. In this study, we generated a PYGM-null induced pluripotent stem cell line and differentiated it into mature RPE to examine structural and functional defects, along with metabolite release into apical and basal media. Mutant RPE exhibited normal photoreceptor outer segment phagocytosis but displayed elevated glycogen levels, reduced transepithelial resistance, and increased cytokine secretion across the epithelial layer compared to isogenic WT controls. Additionally, decreased expression of the visual cycle component, RDH11, encoding 11-cis-retinol dehydrogenase, was observed in PYGM-null RPE. While glycolytic flux and oxidative phosphorylation levels in PYGM-null RPE were near normal, the basal oxygen consumption rate was increased. Oxygen consumption rate in response to physiological levels of lactate was significantly greater in WT than PYGM-null RPE. Inefficient lactate utilization by mutant RPE resulted in higher glucose dependence and increased glucose uptake from the apical medium in the presence of lactate, suggesting a reduced capacity to spare glucose for photoreceptor use. Metabolic tracing confirmed slower 13C-lactate utilization by PYGM-null RPE. These findings have key implications for retinal health since they likely underlie the vision impairment in individuals with McArdle disease.

RNA Sequencing in Disease Diagnosis

RNA sequencing (RNA-seq) enables the accurate measurement of multiple transcriptomic phenotypes for modeling the impacts of disease variants. Advances in technologies, experimental protocols, and analysis strategies are rapidly expanding the application of RNA-seq to identify disease biomarkers, tissue- and cell-type-specific impacts, and the spatial localization of disease-associated mechanisms. Ongoing international efforts to construct biobank-scale transcriptomic repositories with matched genomic data across diverse population groups are further increasing the utility of RNA-seq approaches by providing large-scale normative reference resources. The availability of these resources, combined with improved computational analysis pipelines, has enabled the detection of aberrant transcriptomic phenotypes underlying rare diseases. Further expansion of these resources, across both somatic and developmental tissues, is expected to soon provide unprecedented insights to resolve disease origin, mechanism of action, and causal gene contributions, suggesting the continued high utility of RNA-seq in disease diagnosis.

Low-level brain somatic mutations in exonic regions are collectively implicated in autism with germline mutations in autism risk genes

Low-level somatic mutations in the human brain are implicated in various neurological disorders. The contribution of low-level brain somatic mutations to autism spectrum disorder (ASD), however, remains poorly understood. Here, we performed high-depth exome sequencing with an average read depth of 559.3x in 181 cortical, cerebellar, and peripheral tissue samples to identify brain somatic single nucleotide variants (SNVs) in 24 ASD subjects and 31 controls. We detected ~2.4 brain somatic SNVs per exome per single brain region, with a variant allele frequency (VAF) as low as 0.3%. The mutational profiles, including the number, signature, and type, were not significantly different between the ASD patients and controls. Intriguingly, when considering genes with low-level brain somatic SNVs and ASD risk genes with damaging germline SNVs together, the merged set of genes carrying either somatic or germline SNVs in ASD patients was significantly involved in ASD-associated pathophysiology, including dendrite spine morphogenesis (p = 0.025), mental retardation (p = 0.012), and intrauterine growth retardation (p = 0.012). Additionally, the merged gene set showed ASD-associated spatiotemporal expression in the early and mid-fetal cortex, striatum, and thalamus (all p < 0.05). Patients with damaging mutations in the merged gene set had a greater ASD risk than did controls (odds ratio = 3.92, p = 0.025, 95% confidence interval = 1.12-14.79). The findings of this study suggest that brain somatic SNVs and germline SNVs may collectively contribute to ASD-associated pathophysiology.

Mechanisms and Delivery of tRNA Therapeutics

Transfer ribonucleic acid (tRNA) therapeutics will provide personalized and mutation specific medicines to treat human genetic diseases for which no cures currently exist. The tRNAs are a family of adaptor molecules that interpret the nucleic acid sequences in our genes into the amino acid sequences of proteins that dictate cell function. Humans encode more than 600 tRNA genes. Interestingly, even healthy individuals contain some mutant tRNAs that make mistakes. Missense suppressor tRNAs insert the wrong amino acid in proteins, and nonsense suppressor tRNAs read through premature stop signals to generate full length proteins. Mutations that underlie many human diseases, including neurodegenerative diseases, cancers, and diverse rare genetic disorders, result from missense or nonsense mutations. Thus, specific tRNA variants can be strategically deployed as therapeutic agents to correct genetic defects. We review the mechanisms of tRNA therapeutic activity, the nature of the therapeutic window for nonsense and missense suppression as well as wild-type tRNA supplementation. We discuss the challenges and promises of delivering tRNAs as synthetic RNAs or as gene therapies. Together, tRNA medicines will provide novel treatments for common and rare genetic diseases in humans.

Considerations for reporting variants in novel candidate genes identified during clinical genomic testing

Since the first novel gene discovery for a Mendelian condition was made via exome sequencing (ES), the rapid increase in the number of genes known to underlie Mendelian conditions coupled with the adoption of exome (and more recently, genome) sequencing by diagnostic testing labs has changed the landscape of genomic testing for rare disease. Specifically, many individuals suspected to have a Mendelian condition are now routinely offered clinical ES. This commonly results in a precise genetic diagnosis but frequently overlooks the identification of novel candidate genes. Such candidates are also less likely to be identified in the absence of large-scale gene discovery research programs. Accordingly, clinical laboratories have both the opportunity, and some might argue a responsibility, to contribute to novel gene discovery which should in turn increase the diagnostic yield for many conditions. However, clinical diagnostic laboratories must necessarily balance priorities for throughput, turnaround time, cost efficiency, clinician preferences, and regulatory constraints, and often do not have the infrastructure or resources to effectively participate in either clinical translational or basic genome science research efforts. For these and other reasons, many laboratories have historically refrained from broadly sharing potentially pathogenic variants in novel genes via networks like Matchmaker Exchange, much less reporting such results to ordering providers. Efforts to report such results are further complicated by a lack of guidelines for clinical reporting and interpretation of variants in novel candidate genes. Nevertheless, there are myriad benefits for many stakeholders, including patients/families, clinicians, researchers, if clinical laboratories systematically and routinely identify, share, and report novel candidate genes. To facilitate this change in practice, we developed criteria for triaging, sharing, and reporting novel candidate genes that are most likely to be promptly validated as underlying a Mendelian condition and translated to use in clinical settings.

Expanding the phenotype of PPP1R21-related neurodevelopmental disorder

PPP1R21 encodes for a conserved protein that is involved in endosomal maturation. Biallelic pathogenic variants in PPP1R21 have been associated with a syndromic neurodevelopmental disorder from studying 13 affected individuals. In this report, we present 11 additional individuals from nine unrelated families and their clinical, radiological, and molecular findings. We identified eight different variants in PPP1R21, of which six were novel variants. Global developmental delay and hypotonia are neurological features that were observed in all individuals. There is also a similar pattern of dysmorphic features with coarse faces as a gestalt observed in several individuals. Common findings in 75% of individuals with available brain imaging include delays in myelination, wavy outline of the bodies of the lateral ventricles, and slight prominence of the bodies of the lateral ventricles. PPP1R21-related neurodevelopmental disorder is associated with a consistent phenotype and should be considered in highly consanguineous individuals presenting with developmental delay/intellectual disability along with coarse facial features.

Clinical and molecular characterization of patients with YWHAG-related epilepsy

OBJECTIVE

YWHAG variant alleles have been associated with a rare disease trait whose clinical synopsis includes an early onset epileptic encephalopathy with predominantly myoclonic seizures, developmental delay/intellectual disability, and facial dysmorphisms. Through description of a large cohort, which doubles the number of reported patients, we further delineate the spectrum of YWHAG-related epilepsy.

METHODS

We included in this study 24 patients, 21 new and three previously described, with pathogenic/likely pathogenic variants in YWHAG. We extended the analysis of clinical, electroencephalographic, brain magnetic resonance imaging, and molecular genetic information to 24 previously published patients.

RESULTS

The phenotypic spectrum of YWHAG-related disorders ranges from mild developmental delay to developmental and epileptic encephalopathy (DEE). Epilepsy onset is in the first 2 years of life. Seizure freedom can be achieved in half of the patients (13/24, 54%). Intellectual disability (23/24, 96%), behavioral disorders (18/24, 75%), neurological signs (13/24, 54%), and dysmorphisms (6/24, 25%) are common. A genotype-phenotype correlation emerged, as DEE is more represented in patients with missense variants located in the ligand-binding domain than in those with truncating or missense variants in other domains (90% vs. 19%, p < .001).

SIGNIFICANCE

This study suggests that pathogenic YWHAG variants cause a wide range of clinical presentations with variable severity, ranging from mild developmental delay to DEE. In this allelic series, a genotype-phenotype correlation begins to emerge, potentially providing prognostic information for clinical management and genetic counseling.

Clustering of predicted loss-of-function variants in genes linked with monogenic disease can explain incomplete penetrance

BACKGROUND

Genetic variants that severely alter protein products (e.g. nonsense, frameshift) are often associated with disease. For some genes, these predicted loss-of-function variants (pLoFs) are observed throughout the gene, whilst in others, they occur only at specific locations. We hypothesised that, for genes linked with monogenic diseases that display incomplete penetrance, pLoF variants present in apparently unaffected individuals may be limited to regions where pLoFs are tolerated. To test this, we investigated whether pLoF location could explain instances of incomplete penetrance of variants expected to be pathogenic for Mendelian conditions.

METHODS

We used exome sequence data in 454,773 individuals in the UK Biobank (UKB) to investigate the locations of pLoFs in a population cohort. We counted numbers of unique pLoF, missense, and synonymous variants in UKB in each quintile of the coding sequence (CDS) of all protein-coding genes and clustered the variants using Gaussian mixture models. We limited the analyses to genes with ≥ 5 variants of each type (16,473 genes). We compared the locations of pLoFs in UKB with all theoretically possible pLoFs in a transcript, and pathogenic pLoFs from ClinVar, and performed simulations to estimate the false-positive rate of non-uniformly distributed variants.

RESULTS

For most genes, all variant classes fell into clusters representing broadly uniform variant distributions, but genes in which haploinsufficiency causes developmental disorders were less likely to have uniform pLoF distribution than other genes (P < 2.2 × 10-6). We identified a number of genes, including ARID1B and GATA6, where pLoF variants in the first quarter of the CDS were rescued by the presence of an alternative translation start site and should not be reported as pathogenic. For other genes, such as ODC1, pLoFs were located approximately uniformly across the gene, but pathogenic pLoFs were clustered only at the end, consistent with a gain-of-function disease mechanism.

CONCLUSIONS

Our results suggest the potential benefits of localised constraint metrics and that the location of pLoF variants should be considered when interpreting variants.

Exome-wide analysis implicates rare protein-altering variants in human handedness

Handedness is a manifestation of brain hemispheric specialization. Left-handedness occurs at increased rates in neurodevelopmental disorders. Genome-wide association studies have identified common genetic effects on handedness or brain asymmetry, which mostly involve variants outside protein-coding regions and may affect gene expression. Implicated genes include several that encode tubulins (microtubule components) or microtubule-associated proteins. Here we examine whether left-handedness is also influenced by rare coding variants (frequencies ≤ 1%), using exome data from 38,043 left-handed and 313,271 right-handed individuals from the UK Biobank. The beta-tubulin gene TUBB4B shows exome-wide significant association, with a rate of rare coding variants 2.7 times higher in left-handers than right-handers. The TUBB4B variants are mostly heterozygous missense changes, but include two frameshifts found only in left-handers. Other TUBB4B variants have been linked to sensorineural and/or ciliopathic disorders, but not the variants found here. Among genes previously implicated in autism or schizophrenia by exome screening, DSCAM and FOXP1 show evidence for rare coding variant association with left-handedness. The exome-wide heritability of left-handedness due to rare coding variants was 0.91%. This study reveals a role for rare, protein-altering variants in left-handedness, providing further evidence for the involvement of microtubules and disorder-relevant genes.

Increasing histone acetylation improves sociability and restores learning and memory in KAT6B-haploinsufficient mice

Mutations in genes encoding chromatin modifiers are enriched among mutations causing intellectual disability. The continuing development of the brain postnatally, coupled with the inherent reversibility of chromatin modifications, may afford an opportunity for therapeutic intervention following a genetic diagnosis. Development of treatments requires an understanding of protein function and models of the disease. Here, we provide a mouse model of Say-Barber-Biesecker-Young-Simpson syndrome (SBBYSS) (OMIM 603736) and demonstrate proof-of-principle efficacy of postnatal treatment. SBBYSS results from heterozygous mutations in the KAT6B (MYST4/MORF/QFK) gene and is characterized by intellectual disability and autism-like behaviors. Using human cells carrying SBBYSS-specific KAT6B mutations and Kat6b heterozygous mice (Kat6b+/-), we showed that KAT6B deficiency caused a reduction in histone H3 lysine 9 acetylation. Kat6b+/- mice displayed learning, memory, and social deficits, mirroring SBBYSS individuals. Treatment with a histone deacetylase inhibitor, valproic acid, or an acetyl donor, acetyl-carnitine (ALCAR), elevated histone acetylation levels in the human cells with SBBYSS mutations and in brain and blood cells of Kat6b+/- mice and partially reversed gene expression changes in Kat6b+/- cortical neurons. Both compounds improved sociability in Kat6b+/- mice, and ALCAR treatment restored learning and memory. These data suggest that a subset of SBBYSS individuals may benefit from postnatal therapeutic interventions.

Genome Sequencing of Idiopathic Speech Delay

Genetic investigations of people with speech and language disorders can provide windows into key aspects of human biology. Most genomic research into impaired speech development has so far focused on childhood apraxia of speech (CAS), a rare neurodevelopmental disorder characterized by difficulties with coordinating rapid fine motor sequences that underlie proficient speech. In 2001, pathogenic variants of FOXP2 provided the first molecular genetic accounts of CAS aetiology. Since then, disruptions in several other genes have been implicated in CAS, with a substantial proportion of cases being explained by high-penetrance variants. However, the genetic architecture underlying other speech-related disorders remains less well understood. Thus, in the present study, we used systematic DNA sequencing methods to investigate idiopathic speech delay, as characterized by delayed speech development in the absence of a motor speech diagnosis (such as CAS), a language/reading disorder, or intellectual disability. We performed genome sequencing in a cohort of 23 children with a rigorous diagnosis of idiopathic speech delay. For roughly half of the sample (ten probands), sufficient DNA was also available for genome sequencing in both parents, allowing discovery of de novo variants. In the thirteen singleton probands, we focused on identifying loss-of-function and likely damaging missense variants in genes intolerant to such mutations. We found that one speech delay proband carried a pathogenic frameshift deletion in SETD1A, a gene previously implicated in a broader variable monogenic syndrome characterized by global developmental problems including delayed speech and/or language development, mild intellectual disability, facial dysmorphisms, and behavioural and psychiatric symptoms. Of note, pathogenic SETD1A variants have been independently reported in children with CAS in two separate studies. In other probands in our speech delay cohort, likely pathogenic missense variants were identified affecting highly conserved amino acids in key functional domains of SPTBN1 and ARF3. Overall, this study expands the phenotype spectrum associated with pathogenic SETD1A variants, to also include idiopathic speech delay without CAS or intellectual disability, and suggests additional novel potential candidate genes that may harbour high-penetrance variants that can disrupt speech development.

Deep phenotyping of PROM1-associated retinal degeneration

BACKGROUND/AIMS

The purpose of this study was to investigate retinal structure in detail of subjects with autosomal-dominant (AD) and autosomal-recessive (AR) PROM1-associated retinal degeneration (PROM1-RD), study design: institutional, cross-sectional study.

METHODS

Four eyes from four subjects (three with AD and one with AR) PROM1-RD were investigated by ophthalmic examination including best-corrected visual acuity (BCVA) and multimodal retinal imaging: fundus autofluorescence (FAF), spectral-domain optical coherence tomography (SD-OCT) and adaptive optics scanning light ophthalmoscopy. Quantitative assessment of atrophic lesions determined by FAF, thickness of individual retinal layers and cone photoreceptor quantification was performed.

RESULTS

BCVA ranged from 20/16 to 20/200. Initial pathological changes included the presence of hyperautofluorescent spots on FAF imaging, while later stages demonstrated discrete areas of atrophy. In all patients, thinning of the outer retinal layers on SD-OCT with varying degrees of atrophy could be detected depending on disease-causing variants and age. Cone density was quantified both in central and/or at different eccentricities from the fovea. Longitudinal assessments were possible in two patients.

CONCLUSIONS

PROM1-RD comprises a wide range of clinical phenotypes. Depending on the stage of disease, the cone mosaic in PROM1-RD is relatively preserved and can potentially be targeted by cone-directed interventions.

Deciphering complexity: TULP1 variants linked to an atypical retinal dystrophy phenotype

Introduction: TULP1 exemplifies the remarkable clinical and genetic heterogeneity observed in inherited retinal dystrophies. Our research describes the clinical and molecular characteristics of a patient manifesting an atypical retinal dystrophy pattern, marked by the identification of both a previously unreported and a rarely encountered TULP1 variant. Methods: Whole-exome sequencing was performed to identify potential causative variants. The pathogenicity of the identified TULP1 variants was evaluated through in silico predictors and a minigene splice assay, specifically designed to assess the effect of the unreported TULP1 variant. Results: We identified two TULP1 gene variants in a patient exhibiting unusual and symmetrical alterations in both retinas, characterized by an increase in autofluorescence along the distribution of retinal vessels. These variants included a known rare missense variant, c.1376T>C, and a novel splice site variant, c.822G>T. For the latter variant (c.822G>T), we conducted a minigene splice assay that demonstrated the incorporation of a premature stop codon. This finding suggests a likely activation of the nonsense-mediated mRNA decay mechanism, ultimately resulting in the absence of protein production from this allele. Segregation analysis confirmed that these variants were in trans. Discussion: Our data support that individuals with biallelic TULP1 variants may present with a unique pattern of macular degeneration and periarteriolar vascular pigmentation. This study highlights the importance of further clinical and molecular characterization of TULP1 variants to elucidate genotype-phenotype correlations in the context of inherited retinal dystrophies.

The impact of the Turkish population variome on the genomic architecture of rare disease traits

Purpose

The variome of the Turkish (TK) population, a population with a considerable history of admixture and consanguinity, has not been deeply investigated for insights on the genomic architecture of disease.

Methods

We generated and analyzed a database of variants derived from exome sequencing data of 773 TK unrelated, clinically affected individuals with various suspected Mendelian disease traits and 643 unaffected relatives.

Results

Using uniform manifold approximation and projection, we showed that the TK genomes are more similar to those of Europeans and consist of 2 main subpopulations: clusters 1 and 2 (N = 235 and 1181, respectively), which differ in admixture proportion and variome (https://turkishvariomedb.shinyapps.io/tvdb/). Furthermore, the higher inbreeding coefficient values observed in the TK affected compared with unaffected individuals correlated with a larger median span of long-sized (>2.64 Mb) runs of homozygosity (ROH) regions (P value = 2.09e-18). We show that long-sized ROHs are more likely to be formed on recently configured haplotypes enriched for rare homozygous deleterious variants in the TK affected compared with TK unaffected individuals (P value = 3.35e-11). Analysis of genotype-phenotype correlations reveals that genes with rare homozygous deleterious variants in long-sized ROHs provide the most comprehensive set of molecular diagnoses for the observed disease traits with a systematic quantitative analysis of Human Phenotype Ontology terms.

Conclusion

Our findings support the notion that novel rare variants on newly configured haplotypes arising within the recent past generations of a family or clan contribute significantly to recessive disease traits in the TK population.

Oberg, Clark RD, Wenman C, Loughlin S, Saad R, Ashraf T, Male A, Tadros S, Boostani R, Abdel-Salam GM, Zaki M, Abdalla E, Manzini MC, Pehlivan D, Posey JE, Gibbs RA, Houlden H, Alkuraya FS, Bujakowska K, Maroofian R, Lupski JR, Nguyen LN. Biallelic variation in the choline and ethanolamine transporter FLVCR1 underlies a pleiotropic disease spectrum from adult neurodegeneration to severe developmental disorders

FLVCR1 encodes Feline leukemia virus subgroup C receptor 1 (FLVCR1), a solute carrier (SLC) transporter within the Major Facilitator Superfamily. FLVCR1 is a widely expressed transmembrane protein with plasma membrane and mitochondrial isoforms implicated in heme, choline, and ethanolamine transport. While Flvcr1 knockout mice die in utero with skeletal malformations and defective erythropoiesis reminiscent of Diamond-Blackfan anemia, rare biallelic pathogenic FLVCR1 variants are linked to childhood or adult-onset neurodegeneration of the retina, spinal cord, and peripheral nervous system. We ascertained from research and clinical exome sequencing 27 individuals from 20 unrelated families with biallelic ultra-rare missense and predicted loss-of-function (pLoF) FLVCR1 variant alleles. We characterize an expansive FLVCR1 phenotypic spectrum ranging from adult-onset retinitis pigmentosa to severe developmental disorders with microcephaly, reduced brain volume, epilepsy, spasticity, and premature death. The most severely affected individuals, including three individuals with homozygous pLoF variants, share traits with Flvcr1 knockout mice and Diamond-Blackfan anemia including macrocytic anemia and congenital skeletal malformations. Pathogenic FLVCR1 missense variants primarily lie within transmembrane domains and reduce choline and ethanolamine transport activity compared with wild-type FLVCR1 with minimal impact on FLVCR1 stability or subcellular localization. Several variants disrupt splicing in a mini-gene assay which may contribute to genotype-phenotype correlations. Taken together, these data support an allele-specific gene dosage model in which phenotypic severity reflects residual FLVCR1 activity. This study expands our understanding of Mendelian disorders of choline and ethanolamine transport and demonstrates the importance of choline and ethanolamine in neurodevelopment and neuronal homeostasis.

Exome sequencing improves the molecular diagnostics of paediatric unexplained neurodevelopmental disorders

BACKGROUND

Neurodevelopmental disorders (NDDs) and/or associated multiple congenital abnormalities (MCAs) represent a genetically heterogeneous group of conditions with an adverse prognosis for the quality of intellectual and social abilities and common daily functioning. The rapid development of exome sequencing (ES) techniques, together with trio-based analysis, nowadays leads to up to 50% diagnostic yield. Therefore, it is considered as the state-of-the-art approach in these diagnoses.

RESULTS

In our study, we present the results of ES in a cohort of 85 families with 90 children with severe NDDs and MCAs. The interconnection of the in-house bioinformatic pipeline and a unique algorithm for variant prioritization resulted in a diagnostic yield of up to 48.9% (44/90), including rare and novel causative variants (41/90) and intragenic copy-number variations (CNVs) (3/90). Of the total number of 47 causative variants, 53.2% (25/47) were novel, highlighting the clinical benefit of ES for unexplained NDDs. Moreover, trio-based ES was verified as a reliable tool for the detection of rare CNVs, ranging from intragenic exon deletions (GRIN2A, ZC4H2 genes) to a 6-Mb duplication. The functional analysis using PANTHER Gene Ontology confirmed the involvement of genes with causative variants in a wide spectrum of developmental processes and molecular pathways, which form essential structural and functional components of the central nervous system.

CONCLUSION

Taken together, we present one of the first ES studies of this scale from the central European region. Based on the high diagnostic yield for paediatric NDDs in this study, 48.9%, we confirm trio-based ES as an effective and reliable first-tier diagnostic test in the genetic evaluation of children with NDDs.

Expanding the genetic spectrum of mitochondrial diseases in Tunisia: novel variants revealed by whole-exome sequencing

Introduction: Inherited mitochondrial diseases are the most common group of metabolic disorders caused by a defect in oxidative phosphorylation. They are characterized by a wide clinical and genetic spectrum and can manifest at any age. In this study, we established novel phenotype-genotype correlations between the clinical and molecular features of a cohort of Tunisian patients with mitochondrial diseases. Materials and methods: Whole-exome sequencing was performed on five Tunisian patients with suspected mitochondrial diseases. Then, a combination of filtering and bioinformatics prediction tools was utilized to assess the pathogenicity of genetic variations. Sanger sequencing was subsequently performed to confirm the presence of potential deleterious variants in the patients and verify their segregation within families. Structural modeling was conducted to study the effect of novel variants on the protein structure. Results: We identified two novel homozygous variants in NDUFAF5 (c.827G>C; p.Arg276Pro) and FASTKD2 (c.496_497del; p.Leu166GlufsTer2) associated with a severe clinical form of Leigh and Leigh-like syndromes, respectively. Our results further disclosed two variants unreported in North Africa, in GFM2 (c.569G>A; p.Arg190Gln) and FOXRED1 (c.1261G>A; p.Val421Met) genes, and we described the first case of fumaric aciduria in a Tunisian patient harboring the c.1358T>C; p.Leu453Pro FH variant. Conclusion: Our study expands the mutational and phenotypic spectrum of mitochondrial diseases in Tunisia and highlights the importance of next-generation sequencing to decipher the pathomolecular mechanisms responsible for these disorders in an admixed population.

Systematic analysis of variants escaping nonsense-mediated decay uncovers candidate Mendelian diseases

Protein-truncating variants (PTVs) near the 3' end of genes may escape nonsense-mediated decay (NMD). PTVs in the NMD-escape region (PTVescs) can cause Mendelian disease but are difficult to interpret given their varying impact on protein function. Previously, PTVesc burden was assessed in an epilepsy cohort, but no large-scale analysis has systematically evaluated these variants in rare disease. We performed a retrospective analysis of 29,031 neurodevelopmental disorder (NDD) parent-offspring trios referred for clinical exome sequencing to identify PTVesc de novo mutations (DNMs). We identified 1,376 PTVesc DNMs and 133 genes that were significantly enriched (binomial p < 0.001). The PTVesc-enriched genes included those with PTVescs previously described to cause dominant Mendelian disease (e.g., SEMA6B, PPM1D, and DAGLA). We annotated ClinVar variants for PTVescs and identified 948 genes with at least one high-confidence pathogenic variant. Twenty-two known Mendelian PTVesc-enriched genes had no prior evidence of PTVesc-associated disease. We found 22 additional PTVesc-enriched genes that are not well established to be associated with Mendelian disease, several of which showed phenotypic similarity between individuals harboring PTVesc variants in the same gene. Four individuals with PTVesc mutations in RAB1A had similar phenotypes including NDD and spasticity. PTVesc mutations in IRF2BP1 were found in two individuals who each had severe immunodeficiency manifesting in NDD. Three individuals with PTVesc mutations in LDB1 all had NDD and multiple congenital anomalies. Using a large-scale, systematic analysis of DNMs, we extend the mutation spectrum for known Mendelian disease-associated genes and identify potentially novel disease-associated genes.

The genetic and molecular basis of haploinsufficiency in flowering plants

In diploid organisms, haploinsufficiency can be defined as the requirement for more than one fully functional copy of a gene. In contrast to most genes, whose loss-of-function alleles are recessive, loss-of-function alleles of haploinsufficient genes are dominant. However, forward and reverse genetic screens are biased toward obtaining recessive, loss-of-function mutations, and therefore, dominant mutations of all types are underrepresented in mutant collections. Despite this underrepresentation, haploinsufficient loci have intriguing implications for studies of genome evolution, gene dosage, stability of protein complexes, genetic redundancy, and gene expression. Here we review examples of haploinsufficiency in flowering plants and describe the underlying molecular mechanisms and evolutionary forces driving haploinsufficiency. Finally, we discuss the masking of haploinsufficiency by genetic redundancy, a widespread phenomenon among angiosperms.

Combination of RNAseq and RADseq to Identify Physiological and Adaptive Responses to Acidification in the Eastern Oyster (Crassostrea virginica)

Ocean acidification (OA) is a major stressor threatening marine calcifiers, including the eastern oyster (Crassostrea virginica). In this paper, we provide insight into the molecular mechanisms associated with resilience to OA, with the dual intentions of probing both acclimation and adaptation potential in this species. C. virginica were spawned, and larvae were reared in control or acidified conditions immediately after fertilization. RNA samples were collected from larvae and juveniles, and DNA samples were collected from juveniles after undergoing OA-induced mortality and used to contrast gene expression (RNAseq) and SNP (ddRADseq) profiles from animals reared under both conditions. Results showed convergence of evidence from both approaches, particularly in genes involved in biomineralization that displayed significant changes in variant frequencies and gene expression levels among juveniles that survived acidification as compared to controls. Downregulated genes were related to immune processes, supporting previous studies demonstrating a reduction in immunity from exposure to OA. Acclimation to OA via regulation of gene expression might confer short-term resilience to immediate threats; however, the costs may not be sustainable, underscoring the importance of selection of resilient genotypes. Here, we identified SNPs associated with survival under OA conditions, suggesting that this commercially and ecologically important species might have the genetic variation needed for adaptation to future acidification. The identification of genetic features associated with OA resilience is a highly-needed step for the development of marker-assisted selection of oyster stocks for aquaculture and restoration activities.

GBA1 variants in Brazilian Gaucher disease patients

Gaucher disease (GD) is an autosomal recessive lysosomal disorder caused by pathogenic variants in GBA1 which result in the deficient activity of glucocerebrosidase (GCase). There are few data on the genetic characterization of Brazilian GD patients. This study aimed at characterizing the genotype of 72 unrelated Brazilian GD patients (type I = 63, type II = 4, type III = 5; male = 31). Forty patients were from South Brazil (SB), and 32 were from other regions of Brazil (Others). The exons and exon/intron junctions of GBA1 were analyzed by Sanger sequencing in 8 patients, or by massive parallel sequencing followed by Sanger of exons 9 and 10 in 64 patients. In total, 31 pathogenic variants were identified. The most frequent allele found was N370S (p.(Asn409Ser)) (41.0%), and the most frequent genotype was N370S/RecNciI p.[Asn409Ser];[Leu483Pro;Ala495Pro;Val499=](23.6%). Three variants (N370S - in exon 9, and RecNciI and L444P (p.(Leu483Pro), in exon 10) correspond to 76.3% of total alleles in SB and 59.4% in Others. Two novel variants were described: c.326del(p.(Gln109Argfs*9)) and c.690G>A (p.(?)). Although sequencing all the exons of GBA1 is the gold-standard method for the genetic analysis of GD patients, a step analysis can be proposed for Brazilian patients, starting with analysis of exons 9 and 10. The N370S allele is the most frequently associated with GD in Brazil.

Mining cancer genomes for change-of-metabolic-function mutations

Enzymes with novel functions are needed to enable new organic synthesis techniques. Drawing inspiration from gain-of-function cancer mutations that functionally alter proteins and affect cellular metabolism, we developed METIS (Mutated Enzymes from Tumors In silico Screen). METIS identifies metabolism-altering cancer mutations using mutation recurrence rates and protein structure. We used METIS to screen 298,517 cancer mutations and identify 48 candidate mutations, including those previously identified to alter enzymatic function. Unbiased metabolomic profiling of cells exogenously expressing a candidate mutant (OGDHLp.A400T) supports an altered phenotype that boosts in vitro production of xanthosine, a pharmacologically useful chemical that is currently produced using unsustainable, water-intensive methods. We then applied METIS to 49 million cancer mutations, yielding a refined set of candidates that may impart novel enzymatic functions or contribute to tumor progression. Thus, METIS can be used to identify and catalog potentially-useful cancer mutations for green chemistry and therapeutic applications.

Transcriptomics and chromatin accessibility in multiple African population samples

Mapping the functional human genome and impact of genetic variants is often limited to European-descendent population samples. To aid in overcoming this limitation, we measured gene expression using RNA sequencing in lymphoblastoid cell lines (LCLs) from 599 individuals from six African populations to identify novel transcripts including those not represented in the hg38 reference genome. We used whole genomes from the 1000 Genomes Project and 164 Maasai individuals to identify 8,881 expression and 6,949 splicing quantitative trait loci (eQTLs/sQTLs), and 2,611 structural variants associated with gene expression (SV-eQTLs). We further profiled chromatin accessibility using ATAC-Seq in a subset of 100 representative individuals, to identity chromatin accessibility quantitative trait loci (caQTLs) and allele-specific chromatin accessibility, and provide predictions for the functional effect of 78.9 million variants on chromatin accessibility. Using this map of eQTLs and caQTLs we fine-mapped GWAS signals for a range of complex diseases. Combined, this work expands global functional genomic data to identify novel transcripts, functional elements and variants, understand population genetic history of molecular quantitative trait loci, and further resolve the genetic basis of multiple human traits and disease.

RNA-Based Strategies for Cancer Therapy: In Silico Design and Evaluation of ASOs for Targeted Exon Skipping

Precision medicine in oncology has made significant progress in recent years by approving drugs that target specific genetic mutations. However, many cancer driver genes remain challenging to pharmacologically target ("undruggable"). To tackle this issue, RNA-based methods like antisense oligonucleotides (ASOs) that induce targeted exon skipping (ES) could provide a promising alternative. In this work, a comprehensive computational procedure is presented, focused on the development of ES-based cancer treatments. The procedure aims to produce specific protein variants, including inactive oncogenes and partially restored tumor suppressors. This novel computational procedure encompasses target-exon selection, in silico prediction of ES products, and identification of the best candidate ASOs for further experimental validation. The method was effectively employed on extensively mutated cancer genes, prioritized according to their suitability for ES-based interventions. Notable genes, such as NRAS and VHL, exhibited potential for this therapeutic approach, as specific target exons were identified and optimal ASO sequences were devised to induce their skipping. To the best of our knowledge, this is the first computational procedure that encompasses all necessary steps for designing ASO sequences tailored for targeted ES, contributing with a versatile and innovative approach to addressing the challenges posed by undruggable cancer driver genes and beyond.

De novo PHF5A variants are associated with craniofacial abnormalities, developmental delay, and hypospadias

PURPOSE

The SF3B splicing complex is composed of SF3B1-6 and PHF5A. We report a developmental disorder caused by de novo variants in PHF5A.

METHODS

Clinical, genomic, and functional studies using subject-derived fibroblasts and a heterologous cellular system were performed.

RESULTS

We studied 9 subjects with congenital malformations, including preauricular tags and hypospadias, growth abnormalities, and developmental delay who had de novo heterozygous PHF5A variants, including 4 loss-of-function (LOF), 3 missense, 1 splice, and 1 start-loss variant. In subject-derived fibroblasts with PHF5A LOF variants, wild-type and variant PHF5A mRNAs had a 1:1 ratio, and PHF5A mRNA levels were normal. Transcriptome sequencing revealed alternative promoter use and downregulated genes involved in cell-cycle regulation. Subject and control fibroblasts had similar amounts of PHF5A with the predicted wild-type molecular weight and of SF3B1-3 and SF3B6. SF3B complex formation was unaffected in 2 subject cell lines.

CONCLUSION

Our data suggest the existence of feedback mechanisms in fibroblasts with PHF5A LOF variants to maintain normal levels of SF3B components. These compensatory mechanisms in subject fibroblasts with PHF5A or SF3B4 LOF variants suggest disturbed autoregulation of mutated splicing factor genes in specific cell types, that is, neural crest cells, during embryonic development rather than haploinsufficiency as pathomechanism.

Rare variant associations with plasma protein levels in the UK Biobank

Integrating human genomics and proteomics can help elucidate disease mechanisms, identify clinical biomarkers and discover drug targets1-4. Because previous proteogenomic studies have focused on common variation via genome-wide association studies, the contribution of rare variants to the plasma proteome remains largely unknown. Here we identify associations between rare protein-coding variants and 2,923 plasma protein abundances measured in 49,736 UK Biobank individuals. Our variant-level exome-wide association study identified 5,433 rare genotype-protein associations, of which 81% were undetected in a previous genome-wide association study of the same cohort5. We then looked at aggregate signals using gene-level collapsing analysis, which revealed 1,962 gene-protein associations. Of the 691 gene-level signals from protein-truncating variants, 99.4% were associated with decreased protein levels. STAB1 and STAB2, encoding scavenger receptors involved in plasma protein clearance, emerged as pleiotropic loci, with 77 and 41 protein associations, respectively. We demonstrate the utility of our publicly accessible resource through several applications. These include detailing an allelic series in NLRC4, identifying potential biomarkers for a fatty liver disease-associated variant in HSD17B13 and bolstering phenome-wide association studies by integrating protein quantitative trait loci with protein-truncating variants in collapsing analyses. Finally, we uncover distinct proteomic consequences of clonal haematopoiesis (CH), including an association between TET2-CH and increased FLT3 levels. Our results highlight a considerable role for rare variation in plasma protein abundance and the value of proteogenomics in therapeutic discovery.

aenmd: annotating escape from nonsense-mediated decay for transcripts with protein-truncating variants

SUMMARY

DNA changes that cause premature termination codons (PTCs) represent a large fraction of clinically relevant pathogenic genomic variation. Typically, PTCs induce transcript degradation by nonsense-mediated mRNA decay (NMD) and render such changes loss-of-function alleles. However, certain PTC-containing transcripts escape NMD and can exert dominant-negative or gain-of-function (DN/GOF) effects. Therefore, systematic identification of human PTC-causing variants and their susceptibility to NMD contributes to the investigation of the role of DN/GOF alleles in human disease. Here we present aenmd, a software for annotating PTC-containing transcript-variant pairs for predicted escape from NMD. aenmd is user-friendly and self-contained. It offers functionality not currently available in other methods and is based on established and experimentally validated rules for NMD escape; the software is designed to work at scale, and to integrate seamlessly with existing analysis workflows. We applied aenmd to variants in the gnomAD, Clinvar, and GWAS catalog databases and report the prevalence of human PTC-causing variants in these databases, and the subset of these variants that could exert DN/GOF effects via NMD escape.

AVAILABILITY AND IMPLEMENTATION

aenmd is implemented in the R programming language. Code is available on GitHub as an R-package (github.com/kostkalab/aenmd.git), and as a containerized command-line interface (github.com/kostkalab/aenmd_cli.git).

Monoallelic variation in DHX9, the gene encoding the DExH-box helicase DHX9, underlies neurodevelopment disorders and Charcot-Marie-Tooth disease

DExD/H-box RNA helicases (DDX/DHX) are encoded by a large paralogous gene family; in a subset of these human helicase genes, pathogenic variation causes neurodevelopmental disorder (NDD) traits and cancer. DHX9 encodes a BRCA1-interacting nuclear helicase regulating transcription, R-loops, and homologous recombination and exhibits the highest mutational constraint of all DDX/DHX paralogs but remains unassociated with disease traits in OMIM. Using exome sequencing and family-based rare-variant analyses, we identified 20 individuals with de novo, ultra-rare, heterozygous missense or loss-of-function (LoF) DHX9 variant alleles. Phenotypes ranged from NDDs to the distal symmetric polyneuropathy axonal Charcot-Marie-Tooth disease (CMT2). Quantitative Human Phenotype Ontology (HPO) analysis demonstrated genotype-phenotype correlations with LoF variants causing mild NDD phenotypes and nuclear localization signal (NLS) missense variants causing severe NDD. We investigated DHX9 variant-associated cellular phenotypes in human cell lines. Whereas wild-type DHX9 was restricted to the nucleus, NLS missense variants abnormally accumulated in the cytoplasm. Fibroblasts from an individual with an NLS variant also showed abnormal cytoplasmic DHX9 accumulation. CMT2-associated missense variants caused aberrant nucleolar DHX9 accumulation, a phenomenon previously associated with cellular stress. Two NDD-associated variants, p.Gly411Glu and p.Arg761Gln, altered DHX9 ATPase activity. The severe NDD-associated variant p.Arg141Gln did not affect DHX9 localization but instead increased R-loop levels and double-stranded DNA breaks. Dhx9-/- mice exhibited hypoactivity in novel environments, tremor, and sensorineural hearing loss. All together, these results establish DHX9 as a critical regulator of mammalian neurodevelopment and neuronal homeostasis.

HNRNPC haploinsufficiency affects alternative splicing of intellectual disability-associated genes and causes a neurodevelopmental disorder

Heterogeneous nuclear ribonucleoprotein C (HNRNPC) is an essential, ubiquitously abundant protein involved in mRNA processing. Genetic variants in other members of the HNRNP family have been associated with neurodevelopmental disorders. Here, we describe 13 individuals with global developmental delay, intellectual disability, behavioral abnormalities, and subtle facial dysmorphology with heterozygous HNRNPC germline variants. Five of them bear an identical in-frame deletion of nine amino acids in the extreme C terminus. To study the effect of this recurrent variant as well as HNRNPC haploinsufficiency, we used induced pluripotent stem cells (iPSCs) and fibroblasts obtained from affected individuals. While protein localization and oligomerization were unaffected by the recurrent C-terminal deletion variant, total HNRNPC levels were decreased. Previously, reduced HNRNPC levels have been associated with changes in alternative splicing. Therefore, we performed a meta-analysis on published RNA-seq datasets of three different cell lines to identify a ubiquitous HNRNPC-dependent signature of alternative spliced exons. The identified signature was not only confirmed in fibroblasts obtained from an affected individual but also showed a significant enrichment for genes associated with intellectual disability. Hence, we assessed the effect of decreased and increased levels of HNRNPC on neuronal arborization and neuronal migration and found that either condition affects neuronal function. Taken together, our data indicate that HNRNPC haploinsufficiency affects alternative splicing of multiple intellectual disability-associated genes and that the developing brain is sensitive to aberrant levels of HNRNPC. Hence, our data strongly support the inclusion of HNRNPC to the family of HNRNP-related neurodevelopmental disorders.

Variants in BSN gene associated with epilepsy with favourable outcome

BACKGROUND

BSN gene encodes Bassoon, an essential protein to assemble the cytomatrix at the active zone of neurotransmitter release. This study aims to explore the relationship between BSN variants and epilepsy.

METHODS

Whole-exome sequencing was performed in a cohort of 313 cases (trios) with epilepsies of unknown causes. Additional cases with BSN variants were collected from China Epilepsy Gene V.1.0 Matching Platform. The Clinical Validity Framework of ClinGen was used to evaluate the relationship between BSN variants and epilepsy.

RESULTS

Four pairs of compound heterozygous variants and one cosegregating heterozygous missense variant in BSN were identified in five unrelated families. These variants presented statistically higher frequency in the case cohort than in controls. Additional two de novo heterozygous nonsense variants and one cosegregating heterozygous missense variant were identified in three unrelated cases from the gene matching platform, which were not present in the Genome Aggregation Database. The missense variants tended to be located in C-terminus, including the two monoallelic missense variants. Protein modelling showed that at least one missense variant in each pair of compound heterozygous variants had hydrogen bond alterations. Clinically, two cases were diagnosed as idiopathic generalised epilepsy, two as focal epilepsy and the remaining four as epilepsy with febrile seizures plus. Seven out of eight probands showed infancy or childhood-onset epilepsy. Eight out of 10 affected individuals had a history of febrile convulsions. All the cases were seizure-free. The cases with monoallelic variants achieved seizure-free without treatment or under monotherapy, while cases with biallelic missense variants mostly required combined therapy. The evidence from ClinGen Framework suggested an association between BSN variants and epilepsy.

CONCLUSION

The BSN gene was potentially a novel candidate gene for epilepsy. The phenotypical severity was associated with the genotypes and the molecular subregional effects of the variants.

PE-STOP: A versatile tool for installing nonsense substitutions amenable for precise reversion

The rapid advances in genome editing technologies have revolutionized the study of gene functions in cell or animal models. The recent generation of double-stranded DNA cleavage-independent base editors has been suitably adapted for interrogation of protein-coding genes on the basis of introducing premature stop codons or disabling the start codons. However, such versions of stop/start codon-oriented genetic tools still present limitations on their versatility, base-level precision, and target specificity. Here, we exploit a newly developed prime editor (PE) that differs from base editors by its adoption of a reverse transcriptase activity, which enables incorporation of various types of precise edits templated by a specialized prime editing guide RNA. Based on such a versatile platform, we established a prime editing-empowered method (PE-STOP) for installation of nonsense substitutions, providing a complementary approach to the present gene-targeting tools. PE-STOP is bioinformatically predicted to feature substantially expanded coverage in the genome space. In practice, PE-STOP introduces stop codons with good efficiencies in human embryonic kidney 293T and N2a cells (with medians of 29% [ten sites] and 25% [four sites] editing efficiencies, respectively), while exhibiting minimal off-target effects and high on-target precision. Furthermore, given the fact that PE installs prime editing guide RNA-templated mutations, we introduce a unique strategy for precise genetic rescue of PE-STOP-dependent nonsense mutation via the same PE platform. Altogether, the present work demonstrates a versatile and specific tool for gene inactivation and for functional interrogation of nonsense mutations.

Deep intronic mutation in CRTAP results in unstable isoforms of the protein to induce type I collagen aggregation in a lethal type of osteogenesis imperfecta type VII

Genetic mutations are involved in Mendelian disorders. Unbuffered intronic mutations in gene variants can generate aberrant splice sites in mutant transcripts, resulting in mutant isoforms of proteins with modulated expression, stability, and function in diseased cells. Here, we identify a deep intronic variant, c.794_1403A>G, in CRTAP by genome sequencing of a male fetus with osteogenesis imperfecta (OI) type VII. The mutation introduces cryptic splice sites in intron-3 of CRTAP, resulting in two mature mutant transcripts with cryptic exons. While transcript-1 translates to a truncated isoform (277 amino acids) with thirteen C-terminal non-wild-type amino acids, transcript-2 translates to a wild-type protein sequence, except that this isoform contains an in-frame fusion of non-wild-type twenty-five amino acids in a tetratricopeptide repeat sequence. Both mutant isoforms of CRTAP are unstable due to the presence of a unique 'GWxxI' degron, which finally leads to loss of proline hydroxylation and aggregation of type I collagen. Although type I collagen aggregates undergo autophagy, the overall proteotoxicity resulted in death of the proband cells by senescence. In summary, we present a genetic disease pathomechanism by linking a novel deep intronic mutation in CRTAP to unstable mutant isoforms of the protein in lethal OI type VII.

Jansen-de Vries syndrome: Expansion of the PPM1D clinical and phenotypic spectrum in 34 families

Jansen-de Vries syndrome (JdVS) is a neurodevelopmental condition attributed to pathogenic variants in Exons 5 and 6 of PPM1D. As the full phenotypic spectrum and natural history remain to be defined, we describe a large cohort of children and adults with JdVS. This is a retrospective cohort study of 37 individuals from 34 families with disease-causing variants in PPM1D leading to JdVS. Clinical data were provided by treating physicians and/or families. Of the 37 individuals, 27 were male and 10 female, with median age 8.75 years (range 8 months to 62 years). Four families document autosomal dominant transmission, and 32/34 probands were diagnosed via exome sequencing. The facial gestalt, including a broad forehead and broad mouth with a thin and tented upper lip, was most recognizable between 18 and 48 months of age. Common manifestations included global developmental delay (35/36, 97%), hypotonia (25/34, 74%), short stature (14/33, 42%), constipation (22/31, 71%), and cyclic vomiting (6/35, 17%). Distinctive personality traits include a hypersocial affect (21/31, 68%) and moderate-to-severe anxiety (18/28, 64%). In conclusion, JdVS is a clinically recognizable neurodevelopmental syndrome with a characteristic personality and distinctive facial features. The association of pathogenic variants in PPM1D with cyclic vomiting bears not only medical attention but also further pathogenic and mechanistic evaluation.

A recurrent de novo variant in NUSAP1 escapes nonsense-mediated decay and leads to microcephaly, epilepsy, and developmental delay

NUSAP1 encodes a cell cycle-dependent protein with key roles in mitotic progression, spindle formation, and microtubule stability. Both over- and under-expression of NUSAP1 lead to dysregulation of mitosis and impaired cell proliferation. Through exome sequencing and Matchmaker Exchange, we identified two unrelated individuals with the same recurrent, de novo heterozygous variant (NM_016359.5 c.1209C > A; p.(Tyr403Ter)) in NUSAP1. Both individuals had microcephaly, severe developmental delay, brain abnormalities, and seizures. The gene is predicted to be tolerant of heterozygous loss-of-function mutations, and we show that the mutant transcript escapes nonsense mediated decay, suggesting that the mechanism is likely dominant-negative or toxic gain of function. Single-cell RNA-sequencing of an affected individual's post-mortem brain tissue indicated that the NUSAP1 mutant brain contains all main cell lineages, and that the microcephaly could not be attributed to loss of a specific cell type. We hypothesize that pathogenic variants in NUSAP1 lead to microcephaly possibly by an underlying defect in neural progenitor cells.