Genome and RNA sequencing were essential to reveal cryptic intronic variants associated to defective ATP6AP1 mRNA processing.

The diagnosis of Mendelian disorders has notably advanced with integration of whole exome and genome sequencing (WES and WGS) in clinical practice. However, challenges in variant interpretation and uncovered variants by WES still leave a substantial percentage of patients undiagnosed. In this context, integrating RNA sequencing (RNA-seq) improves diagnostic workflows, particularly for WES inconclusive cases. Additionally, functional studies are often necessary to elucidate the impact of prioritized variants on gene expression and protein function. Our study focused on three unrelated male patients (P1-P3) with ATP6AP1-CDG (congenital disorder of glycosylation), presenting with intellectual disability and varying degrees of hepatopathy, glycosylation defects, and an initially inconclusive diagnosis through WES. Subsequent RNA-seq was pivotal in identifying the underlying genetic causes in P1 and P2, detecting ATP6AP1 underexpression and aberrant splicing. Molecular studies in fibroblasts confirmed these findings and identified the rare intronic variants c.289-233C > T and c.289-289G > A in P1 and P2, respectively. Trio-WGS also revealed the variant c.289-289G > A in P3, which was a de novo change in both patients. Functional assays expressing the mutant alleles in HAP1 cells demonstrated the pathogenic impact of these variants by reproducing the splicing alterations observed in patients. Our study underscores the role of RNA-seq and WGS in enhancing diagnostic rates for genetic diseases such as CDG, providing new insights into ATP6AP1-CDG molecular bases by identifying the first two deep intronic variants in this X-linked gene. Additionally, our study highlights the need to integrate RNA-seq and WGS, followed by functional validation, in routine diagnostics for a comprehensive evaluation of patients with an unidentified molecular etiology.

CRISPR/Cas9-based functional genomics strategy to decipher the pathogenicity of genetic variants in inherited metabolic disorders.

The determination of the functional impact of variants of uncertain significance (VUS) is one of the major bottlenecks in the diagnostic workflow of inherited genetic diseases. To face this problem, we set up a CRISPR/Cas9-based strategy for knock-in cellular model generation, focusing on inherited metabolic disorders (IMDs). We selected variants in seven IMD-associated genes, including seven reported disease-causing variants and four benign/likely benign variants. Overall, 11 knock-in cell models were generated via homology-directed repair in HAP1 haploid cells using CRISPR/Cas9. The functional impact of the variants was determined by analyzing the characteristic biochemical alterations of each disorder. Functional studies performed in knock-in cell models showed that our approach accurately distinguished the functional effect of pathogenic from non-pathogenic variants in a reliable manner in a wide range of IMDs. Our study provides a generic approach to assess the functional impact of genetic variants to improve IMD diagnosis and this tool could emerge as a promising alternative to invasive tests, such as muscular or skin biopsies. Although the study has been performed only in IMDs, this strategy is generic and could be applied to other genetic disorders.

Functional Evidence of CCDC186 as a New Disease-Associated Gene with Endocrine and Central Nervous System Alterations.

CCDC186 protein is involved in the maturation of dense-core vesicles (DCVs) in the trans-Golgi network in neurons and endocrine cells. Mutations in genes involved in DCV regulation, other than CCDC186, have been described in patients with neurodevelopmental disorders. To date, only one patient, within a large sequencing study of 1000 cases, and a single case report with variants in CCDC186, had previously been described. However, no functional studies in any of these two cases had been performed. We identified three patients from two gypsy families, unrelated to each other, with mutations in the CCDC186 gene. Clinically, all patients presented with seizures, frontotemporal atrophy, hypomyelination, recurrent infections, and endocrine disturbances such as severe non-ketotic hypoglycemia. Low levels of cortisol, insulin, or growth hormone could only be verified in one patient. All of them had a neonatal onset and died between 7 months and 4 years of age. Whole exome sequencing identified a homozygous variant in the CCDC186 gene (c.2215C>T, p.Arg739Ter) in the index patients of both families. Protein expression studies demonstrated that CCDC186 was almost undetectable in fibroblasts and muscle tissue. These observations correlated with the transcriptomic analysis performed in fibroblasts in one of the patients, which showed a significant reduction of CCDC186 mRNA levels. Our study provides functional evidence that mutations in this gene have a pathogenic effect on the protein and reinforces CCDC186 as a new disease-associated gene. In addition, mutations in CCDC186 could explain the combined endocrine and neurologic alterations detected in our patients.

Leigh syndrome is the main clinical characteristic of PTCD3 deficiency.

Mitochondrial translation defects are a continuously growing group of disorders showing a large variety of clinical symptoms including a wide range of neurological abnormalities. To date, mutations in PTCD3, encoding a component of the mitochondrial ribosome, have only been reported in a single individual with clinical evidence of Leigh syndrome. Here, we describe three additional PTCD3 individuals from two unrelated families, broadening the genetic and phenotypic spectrum of this disorder, and provide definitive evidence that PTCD3 deficiency is associated with Leigh syndrome. The patients presented in the first months of life with psychomotor delay, respiratory insufficiency and feeding difficulties. The neurologic phenotype included dystonia, optic atrophy, nystagmus and tonic-clonic seizures. Brain MRI showed optic nerve atrophy and thalamic changes, consistent with Leigh syndrome. WES and RNA-seq identified compound heterozygous variants in PTCD3 in both families: c.[1453-1G>C];[1918C>G] and c.[710del];[902C>T]. The functional consequences of the identified variants were determined by a comprehensive characterization of the mitochondrial function. PTCD3 protein levels were significantly reduced in patient fibroblasts and, consistent with a mitochondrial translation defect, a severe reduction in the steady state levels of complexes I and IV subunits was detected. Accordingly, the activity of these complexes was also low, and high-resolution respirometry showed a significant decrease in the mitochondrial respiratory capacity. Functional complementation studies demonstrated the pathogenic effect of the identified variants since the expression of wild-type PTCD3 in immortalized fibroblasts restored the steady-state levels of complexes I and IV subunits as well as the mitochondrial respiratory capacity. Additionally, minigene assays demonstrated that three of the identified variants were pathogenic by altering PTCD3 mRNA processing. The fourth variant was a frameshift leading to a truncated protein. In summary, we provide evidence of PTCD3 involvement in human disease confirming that PTCD3 deficiency is definitively associated with Leigh syndrome.

Diagnostic Odyssey in an Adult Patient with Ophthalmologic Abnormalities and Hearing Loss: Contribution of RNA-Seq to the Diagnosis of a PEX1 Deficiency.

Peroxisomal biogenesis disorders (PBDs) are a heterogeneous group of genetic diseases. Multiple peroxisomal pathways are impaired, and very long chain fatty acids (VLCFA) are the first line biomarkers for the diagnosis. The clinical presentation of PBDs may range from severe, lethal multisystemic disorders to milder, late-onset disease. The vast majority of PBDs belong to Zellweger Spectrum Disordes (ZSDs) and represents a continuum of overlapping clinical symptoms, with Zellweger syndrome being the most severe and Heimler syndrome the less severe disease. Mild clinical conditions frequently present normal or slight biochemical alterations, making the diagnosis of these patients challenging. In the present study we used a combined WES and RNA-seq strategy to diagnose a patient presenting with retinal dystrophy as the main clinical symptom. Results showed the patient was compound heterozygous for mutations in PEX1. VLCFA were normal, but retrospective analysis of lysosphosphatidylcholines (LPC) containing C22:0-C26:0 species was altered. This simple test could avoid the diagnostic odyssey of patients with mild phenotype, such as the individual described here, who was diagnosed very late in adult life. We provide functional data in cell line models that may explain the mild phenotype of the patient by demonstrating the hypomorphic nature of a deep intronic variant altering PEX1 mRNA processing.

Clinical implementation of RNA sequencing for Mendelian disease diagnostics.

BACKGROUND: Lack of functional evidence hampers variant interpretation, leaving a large proportion of individuals with a suspected Mendelian disorder without genetic diagnosis after whole genome or whole exome sequencing (WES). Research studies advocate to further sequence transcriptomes to directly and systematically probe gene expression defects. However, collection of additional biopsies and establishment of lab workflows, analytical pipelines, and defined concepts in clinical interpretation of aberrant gene expression are still needed for adopting RNA sequencing (RNA-seq) in routine diagnostics. METHODS: We implemented an automated RNA-seq protocol and a computational workflow with which we analyzed skin fibroblasts of 303 individuals with a suspected mitochondrial disease that previously underwent WES. We also assessed through simulations how aberrant expression and mono-allelic expression tests depend on RNA-seq coverage. RESULTS: We detected on average 12,500 genes per sample including around 60% of all disease genes-a coverage substantially higher than with whole blood, supporting the use of skin biopsies. We prioritized genes demonstrating aberrant expression, aberrant splicing, or mono-allelic expression. The pipeline required less than 1 week from sample preparation to result reporting and provided a median of eight disease-associated genes per patient for inspection. A genetic diagnosis was established for 16% of the 205 WES-inconclusive cases. Detection of aberrant expression was a major contributor to diagnosis including instances of 50% reduction, which, together with mono-allelic expression, allowed for the diagnosis of dominant disorders caused by haploinsufficiency. Moreover, calling aberrant splicing and variants from RNA-seq data enabled detecting and validating splice-disrupting variants, of which the majority fell outside WES-covered regions. CONCLUSION: Together, these results show that streamlined experimental and computational processes can accelerate the implementation of RNA-seq in routine diagnostics.