Increased Diagnostic Yield by Reanalysis of Whole Exome Sequencing Data in Mitochondrial Disease.

Background: The genetic diagnosis of mitochondrial disorders is complicated by its genetic and phenotypic complexity. Next generation sequencing techniques have much improved the diagnostic yield for these conditions. A cohort of individuals with multiple respiratory chain deficiencies, reported in the literature 10 years ago, had a diagnostic rate of 60% by whole exome sequencing (WES) but 40% remained undiagnosed. Objective: We aimed to identify a genetic diagnosis by reanalysis of the WES data for the undiagnosed arm of this 10-year-old cohort of patients with suspected mitochondrial disorders. Methods: The WES data was transferred and processed by the RD-Connect Genome-Phenome Analysis Platform (GPAP) using their standardized pipeline. Variant prioritisation was carried out on the RD-Connect GPAP. Results: Singleton WES data from 14 individuals was reanalysed. We identified a possible or likely genetic diagnosis in 8 patients (8/14, 57%). The variants identified were in a combination of mitochondrial DNA (n = 1, MT-TN), nuclear encoded mitochondrial genes (n = 2, PDHA1, and SUCLA2) and nuclear genes associated with nonmitochondrial disorders (n = 5, PNPLA2, CDC40, NBAS and SLC7A7). Variants in both the NBAS and CDC40 genes were established as disease causing after the original cohort was published. We increased the diagnostic yield for the original cohort by 15% without generating any further genomic data. CONCLUSIONS: In the era of multiomics we highlight that reanalysis of existing WES data is a valid tool for generating additional diagnosis in patients with suspected mitochondrial disease, particularly when more time has passed to allow for new bioinformatic pipelines to emerge, for the development of new tools in variant interpretation aiding in reclassification of variants and the expansion of scientific knowledge on additional genes.

Unravelling undiagnosed rare disease cases by HiFi long-read genome sequencing.

Solve-RD is a pan-European rare disease (RD) research program that aims to identify disease-causing genetic variants in previously undiagnosed RD families. We utilised 10-fold coverage HiFi long-read sequencing (LRS) for detecting causative structural variants (SVs), single nucleotide variants (SNVs), insertion-deletions (InDels), and short tandem repeat (STR) expansions in extensively studied RD families without clear molecular diagnoses. Our cohort includes 293 individuals from 114 genetically undiagnosed RD families selected by European Rare Disease Network (ERN) experts. Of these, 21 families were affected by so-called 'unsolvable' syndromes for which genetic causes remain unknown, and 93 families with at least one individual affected by a rare neurological, neuromuscular, or epilepsy disorder without genetic diagnosis despite extensive prior testing. Clinical interpretation and orthogonal validation of variants in known disease genes yielded thirteen novel genetic diagnoses due to de novo and rare inherited SNVs, InDels, SVs, and STR expansions. In an additional four families, we identified a candidate disease-causing SV affecting several genes including an MCF2 / FGF13 fusion and PSMA3 deletion. However, no common genetic cause was identified in any of the 'unsolvable' syndromes. Taken together, we found (likely) disease-causing genetic variants in 13.0% of previously unsolved families and additional candidate disease-causing SVs in another 4.3% of these families. In conclusion, our results demonstrate the added value of HiFi long-read genome sequencing in undiagnosed rare diseases.

Partial loss of desmin expression due to a leaky splice site variant in the human DES gene is associated with neuromuscular transmission defects.

Recessive desminopathies are rare and often present as severe early-onset myopathy. Here we report a milder phenotype in three unrelated patients from southern India (2 M, 1F) aged 16, 21, and 22 years, who presented with childhood-onset, gradually progressive, fatigable limb-girdle weakness, ptosis, speech and swallowing difficulties, without cardiac involvement. Serum creatine kinase was elevated, and repetitive nerve stimulation showed decrement in all. Clinical improvement was noted with pyridostigmine and salbutamol in two patients. All three patients had a homozygous substitution in intron 5: DES(NM_001927.4):c.1023+5G>A, predicted to cause a donor splice site defect. Muscle biopsy with ultrastructural analysis suggested myopathy with myofibrillar disarray, and immunohistochemistry showed partial loss of desmin with some residual staining, while western blot analysis showed reduced desmin. RT-PCR of patient muscle RNA revealed two transcripts: a reduced normal desmin transcript and a larger abnormal transcript suggesting leaky splicing at the intron 5 donor site. Sequencing of the PCR products confirmed the inclusion of intron 5 in the longer transcript, predicted to cause a premature stop codon. Thus, we provide evidence for a leaky splice site causing partial loss of desmin associated with a unique phenotypic presentation of a milder form of desmin-related recessive myopathy overlapping with congenital myasthenic syndrome.

Improving Diagnostic Precision: Phenotype-Driven Analysis Uncovers a Maternal Mosaicism in an Individual with RYR1-Congenital Myopathy.

Congenital myopathies (CMs) are rare genetic disorders for which the diagnostic yield does not typically exceed 60% . We performed deep phenotyping, histopathological studies, clinical exome and trio genome sequencing and a phenotype-driven analysis of the genomic data, that led to the molecular diagnosis in a child with CM. We identified a heterozygous variant in RYR1 in the affected child, inherited from her asymptomatic mother. Given the alignment of the clinical and histopathological phenotype with RYR1-CM, we considered the potential existence of a missing second variant in trans in the proband, but also hypothesized that the variant might be mosaic in the mother, as subsequently demonstrated. Our study is an example of how heterozygous variants inherited from asymptomatic parents are frequently dismissed. When the genotype-phenotype correlation is strong, it is recommended to consider a parental mosaicism.

Clinical and genetic characterisation of a large Indian congenital myasthenic syndrome cohort.

Congenital myasthenic syndromes (CMS) are a rare group of inherited disorders caused by gene defects associated with the neuromuscular junction and potentially treatable with commonly available medications such as acetylcholinesterase inhibitors and ß2 adrenergic receptor agonists. In this study, we identified and genetically characterized the largest cohort of CMS patients from India to date. Genetic testing of clinically suspected patients evaluated in a South Indian hospital during the period 2014-19 was carried out by standard diagnostic gene panel testing or using a two-step method that included hotspot screening followed by whole-exome sequencing. In total, 156 genetically diagnosed patients (141 families) were characterized and the mutational spectrum and genotype-phenotype correlation described. Overall, 87 males and 69 females were evaluated, with the age of onset ranging from congenital to fourth decade (mean 6.6 ± 9.8 years). The mean age at diagnosis was 19 ± 12.8 (1-56 years), with a mean diagnostic delay of 12.5 ± 9.9 (0-49 years). Disease-causing variants in 17 CMS-associated genes were identified in 132 families (93.6%), while in nine families (6.4%), variants in genes not associated with CMS were found. Overall, postsynaptic defects were most common (62.4%), followed by glycosylation defects (21.3%), synaptic basal lamina genes (4.3%) and presynaptic defects (2.8%). Other genes found to cause neuromuscular junction defects (DES, TEFM) in our cohort accounted for 2.8%. Among the individual CMS genes, the most commonly affected gene was CHRNE (39.4%), followed by DOK7 (14.4%), DPAGT1 (9.8%), GFPT1 (7.6%), MUSK (6.1%), GMPPB (5.3%) and COLQ (4.5%). We identified 22 recurrent variants in this study, out of which eight were found to be geographically specific to the Indian subcontinent. Apart from the known common CHRNE variants p.E443Kfs*64 (11.4%) and DOK7 p.A378Sfs*30 (9.3%), we identified seven novel recurrent variants specific to this cohort, including DPAGT1 p.T380I and DES c.1023+5G>A, for which founder haplotypes are suspected. This study highlights the geographic differences in the frequencies of various causative CMS genes and underlines the increasing significance of glycosylation genes (DPAGT1, GFPT1 and GMPPB) as a cause of neuromuscular junction defects. Myopathy and muscular dystrophy genes such as GMPPB and DES, presenting as gradually progressive limb girdle CMS, expand the phenotypic spectrum. The novel genes MACF1 and TEFM identified in this cohort add to the expanding list of genes with new mechanisms causing neuromuscular junction defects.

Phenotypic similarity-based approach for variant prioritization for unsolved rare disease: a preliminary methodological report.

Rare diseases (RD) have a prevalence of not more than 1/2000 persons in the European population, and are characterised by the difficulty experienced in obtaining a correct and timely diagnosis. According to Orphanet, 72.5% of RD have a genetic origin although 35% of them do not yet have an identified causative gene. A significant proportion of patients suspected to have a genetic RD receive an inconclusive exome/genome sequencing. Working towards the International Rare Diseases Research Consortium (IRDiRC)'s goal for 2027 to ensure that all people living with a RD receive a diagnosis within one year of coming to medical attention, the Solve-RD project aims to identify the molecular causes underlying undiagnosed RD. As part of this strategy, we developed a phenotypic similarity-based variant prioritization methodology comparing submitted cases with other submitted cases and with known RD in Orphanet. Three complementary approaches based on phenotypic similarity calculations using the Human Phenotype Ontology (HPO), the Orphanet Rare Diseases Ontology (ORDO) and the HPO-ORDO Ontological Module (HOOM) were developed; genomic data reanalysis was performed by the RD-Connect Genome-Phenome Analysis Platform (GPAP). The methodology was tested in 4 exemplary cases discussed with experts from European Reference Networks. Variants of interest (pathogenic or likely pathogenic) were detected in 8.8% of the 725 cases clustered by similarity calculations. Diagnostic hypotheses were validated in 42.1% of them and needed further exploration in another 10.9%. Based on the promising results, we are devising an automated standardized phenotypic-based re-analysis pipeline to be applied to the entire unsolved cases cohort.

Innovative Computerized Dystrophin Quantification Method Based on Spectral Confocal Microscopy.

Several clinical trials are working on drug development for Duchenne and Becker muscular dystrophy (DMD and BMD) treatment, and, since the expected increase in dystrophin is relatively subtle, high-sensitivity quantification methods are necessary. There is also a need to quantify dystrophin to reach a definitive diagnosis in individuals with mild BMD, and in female carriers. We developed a method for the quantification of dystrophin in DMD and BMD patients using spectral confocal microscopy. It offers the possibility to capture the whole emission spectrum for any antibody, ensuring the selection of the emission peak and allowing the detection of fluorescent emissions of very low intensities. Fluorescence was evaluated first on manually selected regions of interest (ROIs), proving the usefulness of the methodology. Later, ROI selection was automated to make it operator-independent. The proposed methodology correctly classified patients according to their diagnosis, detected even minimal traces of dystrophin, and the results obtained automatically were statistically comparable to the manual ones. Thus, spectral imaging could be implemented to measure dystrophin expression and it could pave the way for detailed analysis of how its expression relates to the clinical course. Studies could be further expanded to better understand the expression of dystrophin-associated protein complexes (DAPCs).

Remote visualization of large-scale genomic alignments for collaborative clinical research and diagnosis of rare diseases.

The Solve-RD project objectives include solving undiagnosed rare diseases (RD) through collaborative research on shared genome-phenome datasets. The RD-Connect Genome-Phenome Analysis Platform (GPAP), for data collation and analysis, and the European Genome-Phenome Archive (EGA), for file storage, are two key components of the Solve-RD infrastructure. Clinical researchers can identify candidate genetic variants within the RD-Connect GPAP and, thanks to the developments presented here as part of joint ELIXIR activities, are able to remotely visualize the corresponding alignments stored at the EGA. The Global Alliance for Genomics and Health (GA4GH) htsget streaming application programming interface (API) is used to retrieve alignment slices, which are rendered by an integrated genome viewer (IGV) instance embedded in the GPAP. As a result, it is no longer necessary for over 11,000 datasets to download large alignment files to visualize them locally. This work highlights the advantages, from both the user and infrastructure perspectives, of implementing interoperability standards for establishing federated genomics data networks.

A Solve-RD ClinVar-based reanalysis of 1522 index cases from ERN-ITHACA reveals common pitfalls and misinterpretations in exome sequencing.

PURPOSE: Within the Solve-RD project (https://solve-rd.eu/), the European Reference Network for Intellectual disability, TeleHealth, Autism and Congenital Anomalies aimed to investigate whether a reanalysis of exomes from unsolved cases based on ClinVar annotations could establish additional diagnoses. We present the results of the "ClinVar low-hanging fruit" reanalysis, reasons for the failure of previous analyses, and lessons learned. METHODS: Data from the first 3576 exomes (1522 probands and 2054 relatives) collected from European Reference Network for Intellectual disability, TeleHealth, Autism and Congenital Anomalies was reanalyzed by the Solve-RD consortium by evaluating for the presence of single-nucleotide variant, and small insertions and deletions already reported as (likely) pathogenic in ClinVar. Variants were filtered according to frequency, genotype, and mode of inheritance and reinterpreted. RESULTS: We identified causal variants in 59 cases (3.9%), 50 of them also raised by other approaches and 9 leading to new diagnoses, highlighting interpretation challenges: variants in genes not known to be involved in human disease at the time of the first analysis, misleading genotypes, or variants undetected by local pipelines (variants in off-target regions, low quality filters, low allelic balance, or high frequency). CONCLUSION: The "ClinVar low-hanging fruit" analysis represents an effective, fast, and easy approach to recover causal variants from exome sequencing data, herewith contributing to the reduction of the diagnostic deadlock.

BiP inactivation due to loss of the deAMPylation function of FICD causes a motor neuron disease.

PURPOSE: The chaperone protein BiP is the master regulator of the unfolded protein response in the endoplasmic reticulum. BiP chaperone activity is regulated by the post-translational modification AMPylation, exclusively provided by FICD. We investigated whether FICD variants identified in patients with motor neuron disease could interfere with BiP activity regulation. METHODS: Exome sequencing was performed to identify causative pathogenic variants associated with motor neuron diseases. Functional studies were conducted on fibroblasts from patients to explore the molecular mechanism of the disease. RESULTS: We identified biallelic variants in FICD causing a neurodegenerative disease of upper and lower motor neurons. Affected individuals harbor a specific missense variant, Arg374His, positioned in the catalytic motif of the enzyme and important for adenosine triphosphate binding. The mutated residue abolishes intramolecular interaction with the regulatory residue Glu234, essential to inhibit AMPylation and to promote de-AMPylation by FICD. Consequently, fibroblasts from patients with FICD variants have abnormally increased levels of AMPylated and thus inactivated BiP. CONCLUSION: Loss of BiP chaperone activity in patients likely results in a chronic impairment of the protein quality control system in the endoplasmic reticulum. These findings will guide the development of therapeutic strategies for motoneuron and related diseases linked to proteotoxic stress.

Systematic Collaborative Reanalysis of Genomic Data Improves Diagnostic Yield in Neurologic Rare Diseases.

Many patients experiencing a rare disease remain undiagnosed even after genomic testing. Reanalysis of existing genomic data has shown to increase diagnostic yield, although there are few systematic and comprehensive reanalysis efforts that enable collaborative interpretation and future reinterpretation. The Undiagnosed Rare Disease Program of Catalonia project collated previously inconclusive good quality genomic data (panels, exomes, and genomes) and standardized phenotypic profiles from 323 families (543 individuals) with a neurologic rare disease. The data were reanalyzed systematically to identify relatedness, runs of homozygosity, consanguinity, single-nucleotide variants, insertions and deletions, and copy number variants. Data were shared and collaboratively interpreted within the consortium through a customized Genome-Phenome Analysis Platform, which also enables future data reinterpretation. Reanalysis of existing genomic data provided a diagnosis for 20.7% of the patients, including 1.8% diagnosed after the generation of additional genomic data to identify a second pathogenic heterozygous variant. Diagnostic rate was significantly higher for family-based exome/genome reanalysis compared with singleton panels. Most new diagnoses were attributable to recent gene-disease associations (50.8%), additional or improved bioinformatic analysis (19.7%), and standardized phenotyping data integrated within the Undiagnosed Rare Disease Program of Catalonia Genome-Phenome Analysis Platform functionalities (18%).

The RD-Connect Genome-Phenome Analysis Platform: Accelerating diagnosis, research, and gene discovery for rare diseases.

Rare disease patients are more likely to receive a rapid molecular diagnosis nowadays thanks to the wide adoption of next-generation sequencing. However, many cases remain undiagnosed even after exome or genome analysis, because the methods used missed the molecular cause in a known gene, or a novel causative gene could not be identified and/or confirmed. To address these challenges, the RD-Connect Genome-Phenome Analysis Platform (GPAP) facilitates the collation, discovery, sharing, and analysis of standardized genome-phenome data within a collaborative environment. Authorized clinicians and researchers submit pseudonymised phenotypic profiles encoded using the Human Phenotype Ontology, and raw genomic data which is processed through a standardized pipeline. After an optional embargo period, the data are shared with other platform users, with the objective that similar cases in the system and queries from peers may help diagnose the case. Additionally, the platform enables bidirectional discovery of similar cases in other databases from the Matchmaker Exchange network. To facilitate genome-phenome analysis and interpretation by clinical researchers, the RD-Connect GPAP provides a powerful user-friendly interface and leverages tens of information sources. As a result, the resource has already helped diagnose hundreds of rare disease patients and discover new disease causing genes.

High diagnostic rate of trio exome sequencing in consanguineous families with neurogenetic diseases.

Consanguineous marriages have a prevalence rate of 24% in Turkey. These carry an increased risk of autosomal recessive genetic conditions, leading to severe disability or premature death, with a significant health and economic burden. A definitive molecular diagnosis could not be achieved in these children previously, as infrastructures and access to sophisticated diagnostic options were limited. We studied the cause of neurogenetic disease in 246 children from 190 consanguineous families recruited in three Turkish hospitals between 2016 and 2020. All patients underwent deep phenotyping and trio whole exome sequencing, and data were integrated in advanced international bioinformatics platforms. We detected causative variants in 119 known disease genes in 72% of families. Due to overlapping phenotypes 52% of the confirmed genetic diagnoses would have been missed on targeted diagnostic gene panels. Likely pathogenic variants in 27 novel genes in 14% of the families increased the diagnostic yield to 86%. Eighty-two per cent of causative variants (141/172) were homozygous, 11 of which were detected in genes previously only associated with autosomal dominant inheritance. Eight families carried two pathogenic variants in different disease genes. De novo (9.3%), X-linked recessive (5.2%) and compound heterozygous (3.5%) variants were less frequent compared to non-consanguineous populations. This cohort provided a unique opportunity to better understand the genetic characteristics of neurogenetic diseases in a consanguineous population. Contrary to what may be expected, causative variants were often not on the longest run of homozygosity and the diagnostic yield was lower in families with the highest degree of consanguinity, due to the high number of homozygous variants in these patients. Pathway analysis highlighted that protein synthesis/degradation defects and metabolic diseases are the most common pathways underlying paediatric neurogenetic disease. In our cohort 164 families (86%) received a diagnosis, enabling prevention of transmission and targeted treatments in 24 patients (10%). We generated an important body of genomic data with lasting impacts on the health and wellbeing of consanguineous families and economic benefit for the healthcare system in Turkey and elsewhere. We demonstrate that an untargeted next generation sequencing approach is far superior to a more targeted gene panel approach, and can be performed without specialized bioinformatics knowledge by clinicians using established pipelines in populations with high rates of consanguinity.

Genome-wide variant calling in reanalysis of exome sequencing data uncovered a pathogenic TUBB3 variant.

Almost half of all individuals affected by intellectual disability (ID) remain undiagnosed. In the Solve-RD project, exome sequencing (ES) datasets from unresolved individuals with (syndromic) ID (n = 1,472 probands) are systematically reanalyzed, starting from raw sequencing files, followed by genome-wide variant calling and new data interpretation. This strategy led to the identification of a disease-causing de novo missense variant in TUBB3 in a girl with severe developmental delay, secondary microcephaly, brain imaging abnormalities, high hypermetropia, strabismus and short stature. Interestingly, the TUBB3 variant could only be identified through reanalysis of ES data using a genome-wide variant calling approach, despite being located in protein coding sequence. More detailed analysis revealed that the position of the variant within exon 5 of TUBB3 was not targeted by the enrichment kit, although consistent high-quality coverage was obtained at this position, resulting from nearby targets that provide off-target coverage. In the initial analysis, variant calling was restricted to the exon targets ± 200 bases, allowing the variant to escape detection by the variant calling algorithm. This phenomenon may potentially occur more often, as we determined that 36 established ID genes have robust off-target coverage in coding sequence. Moreover, within these regions, for 17 genes (likely) pathogenic variants have been identified before. Therefore, this clinical report highlights that, although compute-intensive, performing genome-wide variant calling instead of target-based calling may lead to the detection of diagnostically relevant variants that would otherwise remain unnoticed.

A mosaic PIK3CA variant in a young adult with diffuse gastric cancer: case report.

Hereditary diffuse gastric cancer (HDGC) is associated with germline deleterious variants in CDH1 and CTNNA1. The majority of HDGC-suspected patients are still genetically unresolved, raising the need for identification of novel HDGC predisposing genes. Under the collaborative environment of the SOLVE-RD consortium, re-analysis of whole-exome sequencing data from unresolved gastric cancer cases (n = 83) identified a mosaic missense variant in PIK3CA in a 25-year-old female with diffuse gastric cancer (DGC) without familial history for cancer. The variant, c.3140A>G p.(His1047Arg), a known cancer-related somatic hotspot, was present at a low variant allele frequency (18%) in leukocyte-derived DNA. Somatic variants in PIK3CA are usually associated with overgrowth, a phenotype that was not observed in this patient. This report highlights mosaicism as a potential, and understudied, mechanism in the etiology of DGC.