Reanalysis of RNA sequencing data ends diagnostic odyssey and expands the phenotypic spectrum of congenital titinopathy.

Although next-generation sequencing has enabled diagnoses for many patients with Mendelian disorders, the majority remain undiagnosed. Here, we present a sibling pair who were clinically diagnosed with Escobar syndrome, however targeted gene testing was negative. Exome sequencing (ES), and later genome sequencing (GS), revealed compound heterozygous TTN variants in both siblings, a maternally inherited frameshift variant [(NM_133378.4):c.36812del; p.(Asp12271Valfs*10)], and a paternally inherited missense variant [(NM_133378.4):c.12322G > A; p.(Asp4108Asn)]. This result was considered nondiagnostic due to poor clinical fit and limited pathogenicity evidence for the missense variant of uncertain significance (VUS). Following initial nondiagnostic RNA sequencing (RNAseq) on muscle and further pursuit of other variants detected on the ES/GS, a reanalysis of noncanonical splice sites in the muscle transcriptome identified an out-of-frame exon retraction in TTN, near the known VUS. Interim literature included reports of patients with similar TTN variants who had phenotypic concordance with the siblings, and a diagnosis of a congenital titinopathy was given 4 years after the TTN variants had been initially reported. This report highlights the value of reanalysis of RNAseq with a different approach, expands the phenotypic spectrum of congenital titinopathy and also illustrates how a perceived phenotypic mismatch, and failure to consider known variants, can result in a prolongation of the diagnostic journey.

Full-length Isoform Sequencing for Resolving the Molecular Basis of Charcot-Marie-Tooth 2A.

Objectives: Transcript sequencing of patient-derived samples has been shown to improve the diagnostic yield for solving cases of suspected Mendelian conditions, yet the added benefit of full-length long-read transcript sequencing is largely unexplored. Methods: We applied short-read and full-length transcript sequencing and mitochondrial functional studies to a patient-derived fibroblast cell line from an individual with neuropathy that previously lacked a molecular diagnosis. Results: We identified an intronic homozygous MFN2 c.600-31T>G variant that disrupts the branch point critical for intron 6 splicing. Full-length long-read isoform complementary DNA (cDNA) sequencing after treatment with a nonsense-mediated mRNA decay (NMD) inhibitor revealed that this variant creates 5 distinct altered splicing transcripts. All 5 altered splicing transcripts have disrupted open reading frames and are subject to NMD. Furthermore, a patient-derived fibroblast line demonstrated abnormal lipid droplet formation, consistent with MFN2 dysfunction. Although correctly spliced full-length MFN2 transcripts are still produced, this branch point variant results in deficient MFN2 levels and autosomal recessive Charcot-Marie-Tooth disease, axonal, type 2A (CMT2A). Discussion: This case highlights the utility of full-length isoform sequencing for characterizing the molecular mechanism of undiagnosed rare diseases and expands our understanding of the genetic basis for CMT2A.

Whole-Exome Sequencing Identifies Homozygote Nonsense Variants in LMOD2 Gene Causing Infantile Dilated Cardiomyopathy.

As an essential component of the sarcomere, actin thin filament stems from the Z-disk extend toward the middle of the sarcomere and overlaps with myosin thick filaments. Elongation of the cardiac thin filament is essential for normal sarcomere maturation and heart function. This process is regulated by the actin-binding proteins Leiomodins (LMODs), among which LMOD2 has recently been identified as a key regulator of thin filament elongation to reach a mature length. Few reports have implicated homozygous loss of function variants of LMOD2 in neonatal dilated cardiomyopathy (DCM) associated with thin filament shortening. We present the fifth case of DCM due to biallelic variants in the LMOD2 gene and the second case with the c.1193G>A (p.W398*) nonsense variant identified by whole-exome sequencing. The proband is a 4-month male infant of Hispanic descent with advanced heart failure. Consistent with previous reports, a myocardial biopsy exhibited remarkably short thin filaments. However, compared to other cases of identical or similar biallelic variants, the patient presented here has an unusually late onset of cardiomyopathy during infancy. Herein, we present the phenotypic and histological features of this variant, confirm the pathogenic impact on protein expression and sarcomere structure, and discuss the current knowledge of LMOD2-related cardiomyopathy.

Full-length isoform sequencing for resolving the molecular basis of Charcot-Marie-Tooth 2A.

Objectives: Transcript sequencing of patient derived samples has been shown to improve the diagnostic yield for solving cases of likely Mendelian disorders, yet the added benefit of full-length long-read transcript sequencing is largely unexplored. Methods: We applied short-read and full-length isoform cDNA sequencing and mitochondrial functional studies to a patient-derived fibroblast cell line from an individual with neuropathy that previously lacked a molecular diagnosis. Results: We identified an intronic homozygous MFN2 c.600-31T>G variant that disrupts a branch point critical for intron 6 spicing. Full-length long-read isoform cDNA sequencing after treatment with a nonsense-mediated mRNA decay (NMD) inhibitor revealed that this variant creates five distinct altered splicing transcripts. All five altered splicing transcripts have disrupted open reading frames and are subject to NMD. Furthermore, a patient-derived fibroblast line demonstrated abnormal lipid droplet formation, consistent with MFN2 dysfunction. Although correctly spliced full-length MFN2 transcripts are still produced, this branch point variant results in deficient MFN2 protein levels and autosomal recessive Charcot-Marie-Tooth disease, axonal, type 2A (CMT2A). Discussion: This case highlights the utility of full-length isoform sequencing for characterizing the molecular mechanism of undiagnosed rare diseases and expands our understanding of the genetic basis for CMT2A.

Genetic and functional analysis of a Pacific hagfish opioid system.

The actions of endogenous opioids and nociceptin/orphanin FQ are mediated by four homologous G protein-coupled receptors that constitute the opioid receptor family. However, little is known about opioid systems in cyclostomes (living jawless fish) and how opioid systems might have evolved from invertebrates. Here, we leveraged de novo transcriptome and low-coverage whole-genome assembly in the Pacific hagfish (Eptatretus stoutii) to identify and characterize the first full-length coding sequence for a functional opioid receptor in a cyclostome. Additionally, we define two novel endogenous opioid precursors in this species that predict several novel opioid peptides. Bioinformatic analysis shows no closely related opioid receptor genes in invertebrates with regard either to the genomic organization or to conserved opioid receptor-specific sequences that are common in all vertebrates. Furthermore, no proteins analogous to vertebrate opioid precursors could be identified by genomic searches despite previous claims of protein or RNA-derived sequences in several invertebrate species. The presence of an expressed orthologous receptor and opioid precursors in the Pacific hagfish confirms that a functional opioid system was likely present in the common ancestor of all extant vertebrates some 550 million years ago, earlier than all previous authenticated accounts. We discuss the premise that the cyclostome and vertebrate opioid systems evolved from invertebrate systems concerned with antimicrobial defense and speculate that the high concentrations of opioid precursors in tissues such as the testes, gut, and activated immune cells are key remnants of this evolutionary role.

Mosaic de novo SNRPN gene variant associated with Prader-Willi syndrome.

BACKGROUND: Prader-Willi syndrome (PWS) is an imprinting disorder caused by the absence of paternal expressed genes in the Prader-Willi critical region (PWCR) on chromosome 15q11.2-q13. Three molecular mechanisms have been known to cause PWS, including a deletion in the PWCR, uniparental disomy 15 and imprinting defects. RESULTS: We report the first case of PWS associated with a single-nucleotide SNRPN variant in a 10-year-old girl presenting with clinical features consistent with PWS, including infantile hypotonia and feeding difficulty, developmental delay with cognitive impairment, excessive eating with central obesity, sleep disturbances, skin picking and related behaviour issues. Whole-exome sequencing revealed a de novo mosaic nonsense variant of the SNRPN gene (c.73C>T, p.R25X) in 10% of DNA isolated from buccal cells and 19% of DNA from patient-derived lymphoblast cells. DNA methylation study did not detect an abnormal methylation pattern in the SNRPN locus. Parental origin studies showed a paternal source of an intronic single-nucleotide polymorphism within the locus in proximity to the SNRPN variant. CONCLUSIONS: This is the first report that provides evidence of a de novo point mutation of paternal origin in SNRPN as a new disease-causing mechanism for PWS. This finding suggests that gene sequencing should be considered as part of the diagnostic workup in patients with clinical suspicion of PWS.

Commonalities across computational workflows for uncovering explanatory variants in undiagnosed cases.

PURPOSE: Genomic sequencing has become an increasingly powerful and relevant tool to be leveraged for the discovery of genetic aberrations underlying rare, Mendelian conditions. Although the computational tools incorporated into diagnostic workflows for this task are continually evolving and improving, we nevertheless sought to investigate commonalities across sequencing processing workflows to reveal consensus and standard practice tools and highlight exploratory analyses where technical and theoretical method improvements would be most impactful. METHODS: We collected details regarding the computational approaches used by a genetic testing laboratory and 11 clinical research sites in the United States participating in the Undiagnosed Diseases Network via meetings with bioinformaticians, online survey forms, and analyses of internal protocols. RESULTS: We found that tools for processing genomic sequencing data can be grouped into four distinct categories. Whereas well-established practices exist for initial variant calling and quality control steps, there is substantial divergence across sites in later stages for variant prioritization and multimodal data integration, demonstrating a diversity of approaches for solving the most mysterious undiagnosed cases. CONCLUSION: The largest differences across diagnostic workflows suggest that advances in structural variant detection, noncoding variant interpretation, and integration of additional biomedical data may be especially promising for solving chronically undiagnosed cases.

Synaptic processes and immune-related pathways implicated in Tourette syndrome.

Tourette syndrome (TS) is a neuropsychiatric disorder of complex genetic architecture involving multiple interacting genes. Here, we sought to elucidate the pathways that underlie the neurobiology of the disorder through genome-wide analysis. We analyzed genome-wide genotypic data of 3581 individuals with TS and 7682 ancestry-matched controls and investigated associations of TS with sets of genes that are expressed in particular cell types and operate in specific neuronal and glial functions. We employed a self-contained, set-based association method (SBA) as well as a competitive gene set method (MAGMA) using individual-level genotype data to perform a comprehensive investigation of the biological background of TS. Our SBA analysis identified three significant gene sets after Bonferroni correction, implicating ligand-gated ion channel signaling, lymphocytic, and cell adhesion and transsynaptic signaling processes. MAGMA analysis further supported the involvement of the cell adhesion and trans-synaptic signaling gene set. The lymphocytic gene set was driven by variants in FLT3, raising an intriguing hypothesis for the involvement of a neuroinflammatory element in TS pathogenesis. The indications of involvement of ligand-gated ion channel signaling reinforce the role of GABA in TS, while the association of cell adhesion and trans-synaptic signaling gene set provides additional support for the role of adhesion molecules in neuropsychiatric disorders. This study reinforces previous findings but also provides new insights into the neurobiology of TS.

Disseminated Coccidioidomycosis Treated with Interferon-γ and Dupilumab.

We describe a case of life-threatening disseminated coccidioidomycosis in a previously healthy child. Like most patients with disseminated coccidioidomycosis, this child had no genomic evidence of any known, rare immune disease. However, comprehensive immunologic testing showed exaggerated production of interleukin-4 and reduced production of interferon-γ. Supplementation of antifungal agents with interferon-γ treatment slowed disease progression, and the addition of interleukin-4 and interleukin-13 blockade with dupilumab resulted in rapid resolution of the patient's clinical symptoms. This report shows that blocking of type 2 immune responses can treat infection. This immunomodulatory approach could be used to enhance immune clearance of refractory fungal, mycobacterial, and viral infections. (Supported by the Jeffrey Modell Foundation and the National Institutes of Health.).

Contribution of common and rare variants to bipolar disorder susceptibility in extended pedigrees from population isolates.

Current evidence from case/control studies indicates that genetic risk for psychiatric disorders derives primarily from numerous common variants, each with a small phenotypic impact. The literature describing apparent segregation of bipolar disorder (BP) in numerous multigenerational pedigrees suggests that, in such families, large-effect inherited variants might play a greater role. To identify roles of rare and common variants on BP, we conducted genetic analyses in 26 Colombia and Costa Rica pedigrees ascertained for bipolar disorder 1 (BP1), the most severe and heritable form of BP. In these pedigrees, we performed microarray SNP genotyping of 838 individuals and high-coverage whole-genome sequencing of 449 individuals. We compared polygenic risk scores (PRS), estimated using the latest BP1 genome-wide association study (GWAS) summary statistics, between BP1 individuals and related controls. We also evaluated whether BP1 individuals had a higher burden of rare deleterious single-nucleotide variants (SNVs) and rare copy number variants (CNVs) in a set of genes related to BP1. We found that compared with unaffected relatives, BP1 individuals had higher PRS estimated from BP1 GWAS statistics (P = 0.001 ~ 0.007) and displayed modest increase in burdens of rare deleterious SNVs (P = 0.047) and rare CNVs (P = 0.002 ~ 0.033) in genes related to BP1. We did not observe rare variants segregating in the pedigrees. These results suggest that small-to-moderate effect rare and common variants are more likely to contribute to BP1 risk in these extended pedigrees than a few large-effect rare variants.

Correction: GATAD2B-associated neurodevelopmental disorder (GAND): clinical and molecular insights into a NuRD-related disorder.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

GATAD2B-associated neurodevelopmental disorder (GAND): clinical and molecular insights into a NuRD-related disorder.

PURPOSE: Determination of genotypic/phenotypic features of GATAD2B-associated neurodevelopmental disorder (GAND). METHODS: Fifty GAND subjects were evaluated to determine consistent genotypic/phenotypic features. Immunoprecipitation assays utilizing in vitro transcription-translation products were used to evaluate GATAD2B missense variants' ability to interact with binding partners within the nucleosome remodeling and deacetylase (NuRD) complex. RESULTS: Subjects had clinical findings that included macrocephaly, hypotonia, intellectual disability, neonatal feeding issues, polyhydramnios, apraxia of speech, epilepsy, and bicuspid aortic valves. Forty-one novelGATAD2B variants were identified with multiple variant types (nonsense, truncating frameshift, splice-site variants, deletions, and missense). Seven subjects were identified with missense variants that localized within two conserved region domains (CR1 or CR2) of the GATAD2B protein. Immunoprecipitation assays revealed several of these missense variants disrupted GATAD2B interactions with its NuRD complex binding partners. CONCLUSIONS: A consistent GAND phenotype was caused by a range of genetic variants in GATAD2B that include loss-of-function and missense subtypes. Missense variants were present in conserved region domains that disrupted assembly of NuRD complex proteins. GAND's clinical phenotype had substantial clinical overlap with other disorders associated with the NuRD complex that involve CHD3 and CHD4, with clinical features of hypotonia, intellectual disability, cardiac defects, childhood apraxia of speech, and macrocephaly.

Genetic screening of a large series of North American sporadic and familial frontotemporal dementia cases.

INTRODUCTION: The Advancing Research and Treatment for Frontotemporal Lobar Degeneration (ARTFL) and Longitudinal Evaluation of Familial Frontotemporal Dementia Subjects (LEFFTDS) consortia are two closely connected studies, involving multiple North American centers that evaluate both sporadic and familial frontotemporal dementia (FTD) participants and study longitudinal changes. METHODS: We screened the major dementia-associated genes in 302 sporadic and 390 familial (symptomatic or at-risk) participants enrolled in these studies. RESULTS: Among the sporadic patients, 16 (5.3%) carried chromosome 9 open reading frame 72 (C9orf72), microtubule-associated protein tau (MAPT), and progranulin (GRN) pathogenic variants, whereas in the familial series we identified 207 carriers from 146 families. Of interest, one patient was found to carry a homozygous C9orf72 expansion, while another carried both a C9orf72 expansion and a GRN pathogenic variant. We also identified likely pathogenic variants in the TAR DNA binding protein (TARDBP), presenilin 1 (PSEN1), and valosin containing protein (VCP) genes, and a subset of variants of unknown significance in other rare FTD genes. DISCUSSION: Our study reports the genetic characterization of a large FTD series and supports an unbiased sequencing screen, irrespective of clinical presentation or family history.

Diagnostic utility of transcriptome sequencing for rare Mendelian diseases.

PURPOSE: We investigated the value of transcriptome sequencing (RNAseq) in ascertaining the consequence of DNA variants on RNA transcripts to improve the diagnostic rate from exome or genome sequencing for undiagnosed Mendelian diseases spanning a wide spectrum of clinical indications. METHODS: From 234 subjects referred to the Undiagnosed Diseases Network, University of California-Los Angeles clinical site between July 2014 and August 2018, 113 were enrolled for high likelihood of having rare undiagnosed, suspected genetic conditions despite thorough prior clinical evaluation. Exome or genome sequencing and RNAseq were performed, and RNAseq data was integrated with genome sequencing data for DNA variant interpretation genome-wide. RESULTS: The molecular diagnostic rate by exome or genome sequencing was 31%. Integration of RNAseq with genome sequencing resulted in an additional seven cases with clear diagnosis of a known genetic disease. Thus, the overall molecular diagnostic rate was 38%, and 18% of all genetic diagnoses returned required RNAseq to determine variant causality. CONCLUSION: In this rare disease cohort with a wide spectrum of undiagnosed, suspected genetic conditions, RNAseq analysis increased the molecular diagnostic rate above that possible with genome sequencing analysis alone even without availability of the most appropriate tissue type to assess.

Genome-wide survey of copy number variants finds MAPT duplications in progressive supranuclear palsy.

BACKGROUND: Progressive supranuclear palsy is a neurodegenerative tauopathy manifesting clinically as a progressive akinetic-rigid syndrome. In this study, we sought to identify genetic variants influencing PSP susceptibility through a genome-wide association analysis of a cohort of well-characterized patients who had participated in the Neuroprotection and Natural History in Parkinson Plus Syndromes and Blood Brain Barrier in Parkinson Plus Syndromes studies. METHODS: We genotyped single-nucleotide polymorphisms in 283 PSP cases from the United Kingdom, Germany, and France and compared these with genotypes from 4472 controls. Copy number variants were identified from genotyping data. RESULTS: We observed associations on chromosome 17 within or close to the MAPT gene and explored the genetic architecture at this locus. We confirmed the previously reported association of rs1768208 in the MOBP gene (P = 3.29 × 10-13 ) and rs1411478 in STX6 (P = 3.45 × 10-10 ). The population-attributable risk from the MAPT, MOBP, and STX6 single-nucleotide polymorphisms was found to be 0.37, 0.26, and 0.08, respectively. In addition, we found 2 instances of copy number variants spanning the MAPT gene in patients with PSP. These copy number variants include tau but few other genes within the chromosome 17 haplotype region, providing additional support for the direct pathogenicity of MAPT in PSP. CONCLUSIONS: Clinicians should also be aware of MAPT duplication as a possible genetic cause of PSP, especially in patients presenting with young age at onset. © 2019 International Parkinson and Movement Disorder Society.

Interrogating the Genetic Determinants of Tourette's Syndrome and Other Tic Disorders Through Genome-Wide Association Studies.

OBJECTIVE: Tourette's syndrome is polygenic and highly heritable. Genome-wide association study (GWAS) approaches are useful for interrogating the genetic architecture and determinants of Tourette's syndrome and other tic disorders. The authors conducted a GWAS meta-analysis and probed aggregated Tourette's syndrome polygenic risk to test whether Tourette's and related tic disorders have an underlying shared genetic etiology and whether Tourette's polygenic risk scores correlate with worst-ever tic severity and may represent a potential predictor of disease severity. METHODS: GWAS meta-analysis, gene-based association, and genetic enrichment analyses were conducted in 4,819 Tourette's syndrome case subjects and 9,488 control subjects. Replication of top loci was conducted in an independent population-based sample (706 case subjects, 6,068 control subjects). Relationships between Tourette's polygenic risk scores (PRSs), other tic disorders, ascertainment, and tic severity were examined. RESULTS: GWAS and gene-based analyses identified one genome-wide significant locus within FLT3 on chromosome 13, rs2504235, although this association was not replicated in the population-based sample. Genetic variants spanning evolutionarily conserved regions significantly explained 92.4% of Tourette's syndrome heritability. Tourette's-associated genes were significantly preferentially expressed in dorsolateral prefrontal cortex. Tourette's PRS significantly predicted both Tourette's syndrome and tic spectrum disorders status in the population-based sample. Tourette's PRS also significantly correlated with worst-ever tic severity and was higher in case subjects with a family history of tics than in simplex case subjects. CONCLUSIONS: Modulation of gene expression through noncoding variants, particularly within cortico-striatal circuits, is implicated as a fundamental mechanism in Tourette's syndrome pathogenesis. At a genetic level, tic disorders represent a continuous spectrum of disease, supporting the unification of Tourette's syndrome and other tic disorders in future diagnostic schemata. Tourette's PRSs derived from sufficiently large samples may be useful in the future for predicting conversion of transient tics to chronic tic disorders, as well as tic persistence and lifetime tic severity.

Genetic screen in a large series of patients with primary progressive aphasia.

INTRODUCTION: Primary progressive aphasia (PPA) is a neurological syndrome, associated with both frontotemporal dementia and Alzheimer's disease, in which progressive language impairment emerges as the most salient clinical feature during the initial stages of disease. METHODS: We screened the main genes associated with Alzheimer's disease and frontotemporal dementia for pathogenic and risk variants in a cohort of 403 PPA cases. RESULTS: In this case series study, 14 (3.5%) cases carried (likely) pathogenic variants: four C9orf72 expansions, nine GRN, and one TARDBP mutation. Rare risk variants, TREM2 R47H and MAPT A152T, were associated with a three- to seven-fold increase in risk for PPA. DISCUSSION: Our results show that while pathogenic variants within the most common dementia genes were rarely associated with PPA, these were found almost exclusively in GRN and C9orf72, suggesting that PPA is more TDP43- than tau-related in our series. This is consistent with the finding that PPA frequency in dominantly inherited dementias is the highest in kindreds with GRN variants.

Frontotemporal dementia spectrum: first genetic screen in a Greek cohort.

Frontotemporal dementia (FTD) is a heterogeneous group of neurodegenerative syndromes associated with several causative and susceptibility genes. Herein, we aimed to determine the incidence of the most common causative dementia genes in a cohort of 118 unrelated Greek FTD spectrum patients. We also screened for novel possible disease-associated variants in additional 21 genes associated with FTD or amyotrophic lateral sclerosis. Pathogenic or likely pathogenic variants were identified in 16 cases (13.6%). These included repeat expansions in C9orf72 and loss-of-function GRN variants, and likely pathogenic variants in TARDBP, MAPT, and PSEN1. We also identified 14 variants of unknown significance in other rarer FTD or amyotrophic lateral sclerosis genes that require further segregation and functional analysis. Our genetic screen revealed a high genetic burden in familial Greek FTD cases (30.4%), whereas only two of the sporadic cases (3.5%) carried a likely pathogenic variant. A substantial number of familial cases still remain without an obvious causal variant, suggesting the existence of other FTD genetic causes besides those currently screened in clinical routine.

De Novo Sequence and Copy Number Variants Are Strongly Associated with Tourette Disorder and Implicate Cell Polarity in Pathogenesis.

We previously established the contribution of de novo damaging sequence variants to Tourette disorder (TD) through whole-exome sequencing of 511 trios. Here, we sequence an additional 291 TD trios and analyze the combined set of 802 trios. We observe an overrepresentation of de novo damaging variants in simplex, but not multiplex, families; we identify a high-confidence TD risk gene, CELSR3 (cadherin EGF LAG seven-pass G-type receptor 3); we find that the genes mutated in TD patients are enriched for those related to cell polarity, suggesting a common pathway underlying pathobiology; and we confirm a statistically significant excess of de novo copy number variants in TD. Finally, we identify significant overlap of de novo sequence variants between TD and obsessive-compulsive disorder and de novo copy number variants between TD and autism spectrum disorder, consistent with shared genetic risk.