Search

Abstract Background Rett syndrome (RTT) is a neurodevelopmental disorder mainly caused by mutations in the methyl-CpG-binding protein 2 gene (MECP2). MeCP2 is a multi-functional protein involved in many cellular processes, but the mechanisms by which its dysfunction causes disease are not fully understood. The duplication of the MECP2 gene causes a distinct disorder called MECP2 duplication syndrome (MDS), highlighting the importance of tightly regulating its dosage for proper cellular function. Additionally, some patients with mutations in genes other than MECP2 exhibit phenotypic similarities with RTT, indicating that these genes may also play a role in similar cellular functions. The purpose of this study was to characterise the molecular alterations in patients with RTT in order to identify potential biomarkers or therapeutic targets for this disorder. Methods We used a combination of transcriptomics (RNAseq) and proteomics (TMT mass spectrometry) to characterise the expression patterns in fibroblast cell lines from 22 patients with RTT and detected mutation in MECP2, 15 patients with MDS, 12 patients with RTT-like phenotypes and 13 healthy controls. Transcriptomics and proteomics data were used to identify differentially expressed genes at both RNA and protein levels, which were further inspected via enrichment and upstream regulator analyses and compared to find shared features in patients with RTT. Results We identified molecular alterations in cellular functions and pathways that may contribute to the disease phenotype in patients with RTT, such as deregulated cytoskeletal components, vesicular transport elements, ribosomal subunits and mRNA processing machinery. We also compared RTT expression profiles with those of MDS seeking changes in opposite directions that could lead to the identification of MeCP2 direct targets. Some of the deregulated transcripts and proteins were consistently affected in patients with RTT-like phenotypes, revealing potentially relevant molecular processes in patients with overlapping traits and different genetic aetiology. Conclusions The integration of data in a multi-omics analysis has helped to interpret the molecular consequences of MECP2 dysfunction, contributing to the characterisation of the molecular landscape in patients with RTT. The comparison with MDS provides knowledge of MeCP2 direct targets, whilst the correlation with RTT-like phenotypes highlights processes potentially contributing to the pathomechanism leading these disorders.

Abstract Background Molecular measurements of the genome, the transcriptome, and the epigenome, often termed multi-omics data, provide an in-depth view on biological systems and their integration is crucial for gaining insights in complex regulatory processes. These data can be used to explain disease related genetic variants by linking them to intermediate molecular traits (quantitative trait loci, QTL). Molecular networks regulating cellular processes leave footprints in QTL results as so-called trans-QTL hotspots. Reconstructing these networks is a complex endeavor and use of biological prior information can improve network inference. However, previous efforts were limited in the types of priors used or have only been applied to model systems. In this study, we reconstruct the regulatory networks underlying trans-QTL hotspots using human cohort data and data-driven prior information. Methods We devised a new strategy to integrate QTL with human population scale multi-omics data. State-of-the art network inference methods including BDgraph and glasso were applied to these data. Comprehensive prior information to guide network inference was manually curated from large-scale biological databases. The inference approach was extensively benchmarked using simulated data and cross-cohort replication analyses. Best performing methods were subsequently applied to real-world human cohort data. Results Our benchmarks showed that prior-based strategies outperform methods without prior information in simulated data and show better replication across datasets. Application of our approach to human cohort data highlighted two novel regulatory networks related to schizophrenia and lean body mass for which we generated novel functional hypotheses. Conclusions We demonstrate that existing biological knowledge can improve the integrative analysis of networks underlying trans associations and generate novel hypotheses about regulatory mechanisms.

Bi-allelic variants in ASCC1 cause the ultrarare bone fragility disorder “spinal muscular atrophy with congenital bone fractures-2” (SMABF2). However, the mechanism by which ASCC1 dysfunction leads to this musculoskeletal condition and the nature of the associated bone defect are poorly understood. By exome sequencing, we identified a novel homozygous deletion in ASCC1 in a female infant. She was born with severe muscular hypotonia, inability to breathe and swallow, and virtual absence of spontaneous movements; showed progressive brain atrophy, gracile long bones, very slender ribs, and a femur fracture; and died from respiratory failure aged 3 months. A transiliac bone sample taken postmortem revealed a distinct microstructural bone phenotype with low trabecular bone volume, low bone remodeling, disordered collagen organization, and an abnormally high bone marrow adiposity. Proteomics, RNA sequencing, and qPCR in patient-derived skin fibroblasts confirmed that ASCC1 was hardly expressed on protein and RNA levels compared with healthy controls. Furthermore, we demonstrate that mutated ASCC1 is associated with a downregulation of RUNX2, the master regulator of osteoblastogenesis, and SERPINF1, which is involved in osteoblast and adipocyte differentiation. It also exerts an inhibitory effect on TGF-β/SMAD signaling, which is important for bone development. Additionally, knockdown of ASCC1 in human mesenchymal stromal cells (hMSCs) suppressed their differentiation capacity into osteoblasts while increasing their differentiation into adipocytes. This resulted in reduced mineralization and elevated formation of lipid droplets. These findings shed light onto the pathophysiologic mechanisms underlying SMABF2 and assign a new biological role to ASCC1 acting as an important pro-osteoblastogenic and anti-adipogenic regulator.

Newborn screening (NBS) for hepatorenal tyrosinemia type I (HT1) based on a determination of succinylacetone is performed in countries worldwide. Recently, biallelic pathogenic variants in GSTZ1 underlying maleylacetoacetate isomerase (MAAI) deficiency have been described as a differential diagnosis in individuals with slightly elevated succinylacetone detected by NBS. We report the experience with NBS for HT1 over 53 months in a large German NBS center and the identification and characterization of additional cases with MAAI deficiency, including one individual with a natural history over 32 years. A total of 516,803 children underwent NBS for HT1 at the NBS center in Heidelberg between August 2016 and December 2020. Of 42 children with elevated succinylacetone, HT1 was confirmed in two cases (1 in 258.401). MAAI deficiency was suspected in two cases and genetically confirmed in one who showed traces of succinylacetone in urine. A previously unreported pathogenic GSTZ1 variant was found in the index in a biallelic state. Segregation analysis revealed monoallelic carriership in the index case‘s mother and homozygosity in his father. The 32-year-old father had no medical concerns up to that point and the laboratory work-up was unremarkable. MAAI has to be considered a rare differential diagnosis in NBS for HT1 in cases with slight elevations of succinylacetone to allow for correct counselling and treatment decisions. Our observation of natural history over 32 years adds evidence for a benign clinical course of MAAI deficiency without specific treatment.

Abstract 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.

ABSTRACT Although whole‐exome sequencing and whole‐genome sequencing has tremendously improved our understanding of the genetic etiology of human disorders, about half of the patients still do not receive a molecular diagnosis. The high fraction of variants with uncertain significance and the challenges of interpretation of noncoding variants have urged scientists to implement RNA sequencing (RNA‐seq) in the diagnostic approach as a high throughput assay to complement genomic data with functional evidence. RNA‐seq data can be used to identify aberrantly spliced genes, detect allele‐specific expression, and identify gene expression outliers. Amongst eight studies utilizing RNA‐seq, a mean diagnostic uplift of 15% has been reported. Here, we provide an overview of how RNA‐seq has been implemented to aid in identifying the causal variants of Mendelian disorders.

Serum urate concentration can be studied in large datasets to find genetic and epigenetic loci that may be related to cardiometabolic traits. Here the authors identify and replicate 100 urate-associated CpGs, which provide insights into urate GWAS loci and shared CpGs of urate and cardiometabolic traits.

Many genetic loci have been identified to be associated with kidney disease, but the molecular mechanisms are not well understood. Here, the authors perform epigenome-wide association studies on kidney function measures to identify epigenetic marks and pathways involved in kidney function.

Abstract Mitochondrial disorders are clinically and genetically diverse, with isolated complex III (CIII) deficiency being relatively rare. Here, we describe two affected cousins, presenting with recurrent episodes of severe lactic acidosis, hyperammonaemia, hypoglycaemia and encephalopathy. Genetic investigations in both cases identified a homozygous deletion of exons 2 and 3 of UQCRH, which encodes a structural complex III (CIII) subunit. We generated a mouse model with the equivalent homozygous Uqcrh deletion (Uqcrh−/−), which also presented with lactic acidosis and hyperammonaemia, but had a more severe, non‐episodic phenotype, resulting in failure to thrive and early death. The biochemical phenotypes observed in patient and Uqcrh−/− mouse tissues were remarkably similar, displaying impaired CIII activity, decreased molecular weight of fully assembled holoenzyme and an increase of an unexpected large supercomplex (SXL), comprising mostly of one complex I (CI) dimer and one CIII dimer. This phenotypic similarity along with lentiviral rescue experiments in patient fibroblasts verifies the pathogenicity of the shared genetic defect, demonstrating that the Uqcrh−/− mouse is a valuable model for future studies of human CIII deficiency.

Stavsky et al. examined the effects of deleting the mitochondrial sodium/lithium/calcium exchanger, NCLX, on mitochondrial and synaptic calcium homeostasis, synaptic activity, and plasticity in mice. Having identified a human mutation that impairs NCLX activity and is associated with mental retardation, they show that NCLX is crucial for defining synaptic strength and plasticity, which are pivotal elements of learning and memory.

Aberrant splicing is a major contributor to rare disease, but detection accuracy using current methods is limited. Here, the authors develop an algorithm that detects aberrant splicing and intron retention events from RNA-seq data and apply it to diagnosis in mitochondrial disease.

Abstract Objective Pantothenate kinase 2‐associated neurodegeneration (PKAN) is a rare neurodegenerative disease caused by mutations in the pantothenate kinase 2 (PANK2) gene. PKAN is associated with iron deposition in the basal ganglia and, occasionally, with the occurrence of misshaped erythrocytes (acanthocytes). The aim of this study was to assess residual activity of PANK2 in erythrocytes of PKAN patients and to correlate these data with the type of PANK2 mutations and the progression of neurodegeneration. Methods Residual PANK2 activities in erythrocytes of 14 PKAN patients and 14 related carriers were assessed by a radiometric assay. Clinical data on neurodegeneration included the Barry–Albright Dystonia Scale (BAD‐Scale) besides further general patient features. A molecular visualization and analysis program was used to rationalize the influence of the PKAN causing mutations on a molecular level. Results Erythrocytes of PKAN patients had markedly reduced or no PANK2 activity. However, patients with at least one allele of the c.1583C > T (T528M) or the c.833G > T (R278L) variant exhibited 12–56% of residual PANK2 activity. In line, molecular modeling indicated only minor effects on enzyme structure for these point mutations. On average, these patients with c.1583C > T or c.833G > T variant had lower BAD scores corresponding to lower symptom severity than patients with other PANK2 point mutations. Interpretation Residual erythrocyte PANK2 activity could be a predictor for the progression of neurodegeneration in PKAN patients. Erythrocytes are an interesting patient‐derived cell system with still underestimated diagnostic potential.

Osteogenesis imperfecta (OI) is an inherited genetic disorder characterized by frequent bone fractures and reduced bone mass. Most cases of OI are caused by dominantly inherited heterozygous mutations in one of the two genes encoding type I collagen, COL1A1 and COL1A2. Here we describe a five-year-old boy with typical clinical, radiological and bone ultrastructural features of OI type I. Establishing the molecular genetic cause of his condition proved difficult since clinical exome and whole exome analysis was repeatedly reported negative. Finally, manual analysis of exome data revealed a silent COL1A2 variant c.3597 T > A (NM_000089.4), which we demonstrate activates a cryptic splice site. The newly generated splice acceptor in exon 50 is much more accessible than the wild-type splice-site between the junction of exon 49 and 50, and results in an in-frame deletion of 24 amino acids of the C-terminal propeptide. In vitro collagen expression studies confirmed cellular accumulation and decreased COL1A2 secretion to 45%. This is the first report of a cryptic splice site within the coding region of COL1A2. which results in abnormal splicing causing OI. The experience from this case demonstrates that routine diagnostic approaches may miss cryptic splicing mutations in causative genes due to the lack of universally applicable algorithms for splice-site prediction. In exome-negative cases, in-depth analysis of common causative genes should be conducted and trio-exome analysis is recommended.

ObjectiveThe cytochrome c oxidase assembly factor 7 (COA7) gene encodes a protein localized to mitochondria that is involved in the assembly of mitochondrial respiratory chain complex IV. Here, we report the clinical, genetic and biochemical analysis of a female patient with suspected mitochondrial disorder and novel variants in COA7, that presented with a considerably different phenotype and age of onset than the five COA7 patients reported to date.MethodsWe performed trio-exome sequencing in the affected patient and both parents. To verify the pathogenicity of the detected variants in COA7, mitochondrial enzyme activities and oxygen consumption rate were investigated in fibroblasts of the patient and her parents.ResultsA Chinese girl was referred at 9 months of age with a history of developmental delay and regression since 3 months of age. In the following months, she lost previously acquired skills and developed progressive spasticity of the lower extremities. Trio-exome sequencing revealed compound heterzygous variants in COA7 (c.511G > A/p.Ala171Thr and c.566A > G/p.Asn189Ser). Functional validation experiments revealed isolated complex IV deficiency and a significantly reduced mitochondrial respiration rate in patient-derived fibroblasts.InterpretationHitherto, characteristic features of COA7 patients were described as slowly progressing neuropathy and spinocerebellar ataxia, starting at the toddler age and progressing into adulthood. In contrast, our patient was reported to show developmental delay from 3 months of age, which was found to be due to a rapidly progressive encephalopathy and brain atrophy seen at 9 months of age. Unexpectedly, the genetic investigation revealed a COA7-associated mitochondrial disease, which was confirmed functionally. Thus, this report broadens the genetic and clinical spectrum of this heterogeneous mitochondriopathy and highlights the value of the presented unbiased approach.

Rare diseases, although individually rare, collectively affect approximately 350 million people worldwide. Currently, nearly 6,000 distinct rare disorders with a known molecular basis have been described, yet establishing a specific diagnosis based on the clinical phenotype is challenging. Increasing integration of whole exome sequencing into routine diagnostics of rare diseases is improving diagnostic rates. Nevertheless, about half of the patients do not receive a genetic diagnosis due to the challenges of variant detection and interpretation. During the last years, RNA sequencing is increasingly used as a complementary diagnostic tool providing functional data. Initially, arbitrary thresholds have been applied to call aberrant expression, aberrant splicing, and mono-allelic expression. With the application of RNA sequencing to search for the molecular diagnosis, the implementation of robust statistical models on normalized read counts allowed for the detection of significant outliers corrected for multiple testing. More recently, machine learning methods have been developed to improve the normalization of RNA sequencing read count data by taking confounders into account. Together the methods have increased the power and sensitivity of detection and interpretation of pathogenic variants, leading to diagnostic rates of 10–35% in rare diseases. In this review, we provide an overview of the methods used for RNA sequencing and illustrate how these can improve the diagnostic yield of rare diseases.

Abstract Background Complex I (CI or NADH:ubiquinone oxidoreductase) deficiency is the most frequent cause of mitochondrial respiratory chain defect. Successful attempts to rescue CI function by introducing an exogenous NADH dehydrogenase, such as the NDI1 from Saccharomyces cerevisiae (ScNDI1), have been reported although with drawbacks related to competition with CI. In contrast to ScNDI1, which is permanently active in yeast naturally devoid of CI, plant alternative NADH dehydrogenases (NDH-2) support the oxidation of NADH only when the CI is metabolically inactive and conceivably when the concentration of matrix NADH exceeds a certain threshold. We therefore explored the feasibility of CI rescue by NDH-2 from Arabidopsis thaliana (At) in human CI defective fibroblasts. Results We showed that, other than ScNDI1, two different NDH-2 (AtNDA2 and AtNDB4) targeted to the mitochondria were able to rescue CI deficiency and decrease oxidative stress as indicated by a normalization of SOD activity in human CI-defective fibroblasts. We further demonstrated that when expressed in human control fibroblasts, AtNDA2 shows an affinity for NADH oxidation similar to that of CI, thus competing with CI for the oxidation of NADH as opposed to our initial hypothesis. This competition reduced the amount of ATP produced per oxygen atom reduced to water by half in control cells. Conclusions In conclusion, despite their promising potential to rescue CI defects, due to a possible competition with remaining CI activity, plant NDH-2 should be regarded with caution as potential therapeutic tools for human mitochondrial diseases.

Abstract Despite their importance in determining protein abundance, a comprehensive catalogue of sequence features controlling protein‐to‐mRNA (PTR) ratios and a quantification of their effects are still lacking. Here, we quantified PTR ratios for 11,575 proteins across 29 human tissues using matched transcriptomes and proteomes. We estimated by regression the contribution of known sequence determinants of protein synthesis and degradation in addition to 45 mRNA and 3 protein sequence motifs that we found by association testing. While PTR ratios span more than 2 orders of magnitude, our integrative model predicts PTR ratios at a median precision of 3.2‐fold. A reporter assay provided functional support for two novel UTR motifs, and an immobilized mRNA affinity competition‐binding assay identified motif‐specific bound proteins for one motif. Moreover, our integrative model led to a new metric of codon optimality that captures the effects of codon frequency on protein synthesis and degradation. Altogether, this study shows that a large fraction of PTR ratio variation in human tissues can be predicted from sequence, and it identifies many new candidate post‐transcriptional regulatory elements.

Abstract Background Most research into myocardial infarctions (MIs) have focused on preventative efforts. For survivors, the occurrence of an MI represents a major clinical event that can have long-lasting consequences. There has been little to no research into the molecular changes that can occur as a result of an incident MI. Here, we use three cohorts to identify epigenetic changes that are indicative of an incident MI and their association with gene expression and metabolomics. Results Using paired samples from the KORA cohort, we screened for DNA methylation loci (CpGs) whose change in methylation is potentially indicative of the occurrence of an incident MI between the baseline and follow-up exams. We used paired samples from the NAS cohort to identify 11 CpGs which were predictive in an independent cohort. After removing two CpGs associated with medication usage, we were left with an “epigenetic fingerprint” of MI composed of nine CpGs. We tested this fingerprint in the InCHIANTI cohort where it moderately discriminated incident MI occurrence (AUC = 0.61, P = 6.5 × 10−3). Returning to KORA, we associated the epigenetic fingerprint loci with cis-gene expression and integrated it into a gene expression-metabolomic network, which revealed links between the epigenetic fingerprint CpGs and branched chain amino acid (BCAA) metabolism. Conclusions There are significant changes in DNA methylation after an incident MI. Nine of these CpGs show consistent changes in multiple cohorts, significantly discriminate MI in independent cohorts, and were independent of medication usage. Integration with gene expression and metabolomics data indicates a link between MI-associated epigenetic changes and BCAA metabolism.

Mitochondrial encephalomyopathy, lactic acidosis, and stroke-like episodes (MELAS) is a maternally inherited mitochondrial disease. Most cases of MELAS are caused by the m.3243A > G variant in the MT-TL1 gene encoding tRNALeu(UUR). However, the genetic cause in 10% of patients with MELAS is unknown. We investigated the pathogenicity of the novel mtDNA variant m.9396G > A/MT-CO3 (p.E64K), which affects an extremely conserved amino acid in the CO3 subunit of mitochondrial respiratory chain (MRC) complex IV (CIV) in a patient with MELAS. Biochemical assays of a muscle biopsy confirmed remarkable CIV deficiency, and pathological examination showed ragged red fibers and generalized COX non-reactive muscle fibers. Transfer of the mutant mtDNA into cybrids impaired CIV assembly, followed by remarkable mitochondrial dysfunction and ROS production. Our findings highlight the pathogenicity of a novel m.9396G > A variant and extend the spectrum of pathogenic mtDNA variants.

We present two unrelated Chinese patients with CAD deficiency manifesting with a triad of infantile-onset psychomotor developmental delay with regression, drug-refractory epilepsy, and anaemia with anisopoikilocytosis. Timely translation into uridine supplementation, within 2-months of disease onset, allowed us to stop conventional anti-epileptic drugs and led to dramatic improvement in the clinical symptoms, with prompt cessation of seizures, resolution of anaemia, developmental progress, and prevention of development of severe and non-reversible manifestations. The remarkable recovery and prevention of advanced disease with prompt treatment, highlights the need to act immediately upon genetic diagnosis of a treatable disease. This further reinforces CAD deficiency as a treatable neurometabolic disorder and emphasises the need for a biomarker or genetic new born screening for early identification.

AARS1 deficiency belongs to the group of disorders affecting aminoacyl-tRNA synthetases. To date, AARS1 deficiency has only been linked to neurologic disorders. We report a 6-year-old girl with microcephaly and developmental delay who presented with repeated episodes of acute liver failure. Whole-exome sequencing revealed compound heterozygosity for two missense variants within the AARS1 gene, p.[Leu298Gln];[Arg751Gly]), whose functional relevance was demonstrated by decreased enzymatic activity in fibroblasts. This is the first report that shows that AARS1 variants may be associated with recurrent acute liver failure.

The concept of a mitochondrial disorder was initially described in 1962, in a patient with altered energy metabolism. Over time, mitochondrial energy metabolism has been discovered to be influenced by a vast number of proteins with a multitude of functional roles. Amongst these, defective oxidative phosphorylation arose as the hallmark of mitochondrial disorders. In the premolecular era, the diagnosis of mitochondrial disease was dependent on biochemical criteria, with inherent limitations such as tissue availability and specificity, preanalytical and analytical artifacts, and secondary effects. With the identification of the first mitochondrial disease-causing mutations, the genetic complexity of mitochondrial disorders began to unravel. Mitochondrial dysfunctions can be caused by pathogenic variants in genes encoded by the mitochondrial DNA or the nuclear DNA, and can display heterogenous phenotypic manifestations. The application of next generation sequencing methodologies in diagnostics is proving to be pivotal in finding the molecular diagnosis and has been instrumental in the discovery of a growing list of novel mitochondrial disease genes. In the molecular era, the diagnosis of a mitochondrial disorder, suspected on clinical grounds, is increasingly based on variant detection and associated statistical support, while invasive biopsies and biochemical assays are conducted to an ever-decreasing extent. At present, there is no uniform biochemical or molecular definition for the designation of a disease as a “mitochondrial disorder”. Such designation is currently dependent on the criteria applied, which may encompass clinical, genetic, biochemical, functional, and/or mitochondrial protein localization criteria. Given this variation, numerous gene lists emerge, ranging from 270 to over 400 proposed mitochondrial disease genes. Herein we provide an overview of the mitochondrial disease associated genes and their accompanying challenges.

Summary: The proteasome is the main proteolytic system for targeted protein degradation in the cell and is fine-tuned according to cellular needs. Here, we demonstrate that mitochondrial dysfunction and concomitant metabolic reprogramming of the tricarboxylic acid (TCA) cycle reduce the assembly and activity of the 26S proteasome. Both mitochondrial mutations in respiratory complex I and treatment with the anti-diabetic drug metformin impair 26S proteasome activity. Defective 26S assembly is reversible and can be overcome by supplementation of aspartate or pyruvate. This metabolic regulation of 26S activity involves specific regulation of proteasome assembly factors via the mTORC1 pathway. Of note, reducing 26S activity by metformin confers increased resistance toward the proteasome inhibitor bortezomib, which is reversible upon pyruvate supplementation. Our study uncovers unexpected consequences of defective mitochondrial metabolism for proteasomal protein degradation in the cell, which has important pathophysiological and therapeutic implications.

Mitochondrial diseases are amongst the most genetically and phenotypically diverse groups of inherited diseases. The vast phenotypic overlap with other disease entities together with the absence of reliable biomarkers act as driving forces for the integration of unbiased methodologies early in the diagnostic algorithm, such as whole exome sequencing (WES) and whole genome sequencing (WGS). Such approaches are used in variant discovery and in combination with high-throughput functional assays such as transcriptomics in simultaneous variant discovery and validation. By capturing all genes, they not only increase the diagnostic rate in heterogenous mitochondrial disease patients, but accelerate novel disease gene discovery, and are valuable in side-stepping the risk of overlooking unexpected or even treatable genetic disease diagnoses.

Background: Mitochondrial disorders are a group of rare diseases, caused by nuclear or mitochondrial DNA mutations. Their marked clinical and genetic heterogeneity as well as referral and ascertainment biases render phenotype-based prevalence estimations difficult. Here we calculated the lifetime risk of all known autosomal recessive mitochondrial disorders on basis of genetic data. Methods: We queried the publicly available Genome Aggregation Database (gnomAD) and our in-house exome database to assess the allele frequency of disease-causing variants in genes associated with autosomal recessive mitochondrial disorders. Based on this, we estimated the lifetime risk of 249 autosomal recessive mitochondrial disorders. Three of these disorders and phenylketonuria (PKU) served as a proof of concept since calculations could be aligned with known birth prevalence data from newborn screening reports. Findings: The estimated lifetime risks are very close to newborn screening data (where available), supporting the validity of the approach. For example, calculated lifetime risk of PKU (16·0/100,000) correlates well with known birth prevalence data (18·7/100,000). The combined estimated lifetime risk of 249 investigated mitochondrial disorders is 31·8 (20·9–50·6)/100,000 in our in-house database, 48·4 (40·3–58·5)/100,000 in the European gnomAD dataset, and 31·1 (26·7–36·3)/100,000 in the global gnomAD dataset. The disorders with the highest lifetime risk (> 3 per 100,000) were, in all datasets, those caused by mutations in the SPG7, ACADM, POLG and SLC22A5 genes. Interpretation: We provide a population-genetic estimation on the lifetime risk of an entire class of monogenic disorders. Our findings reveal the substantial cumulative prevalence of autosomal recessive mitochondrial disorders, far above previous estimates. These data will be very important for assigning diagnostic a priori probabilities, and for resource allocation in therapy development, public health management and biomedical research. Funding: German Federal Ministry of Education and Research. Keywords: Autosomal recessive mitochondrial disorders, Population genetics, Prevalence, Lifetime risk, POLG, SPG7

The regulation of mitochondrial proteome is unique in that its components have origins in both mitochondria and nucleus. With the development of OMICS technologies, emerging evidence indicates an interaction between mitochondria and nucleus based not only on the proteins but also on the non-coding RNAs (ncRNAs). It is now accepted that large parts of the non‐coding genome are transcribed into various ncRNA species. Although their characterization has been a hot topic in recent years, the function of the majority remains unknown. Recently, ncRNA species microRNA (miRNA) and long-non coding RNAs (lncRNA) have been gaining attention as direct or indirect modulators of the mitochondrial proteome homeostasis. These ncRNA can impact mitochondria indirectly by affecting transcripts encoding for mitochondrial proteins in the cytoplasm. Furthermore, reports of mitochondria-localized miRNAs, termed mitomiRs, and lncRNAs directly regulating mitochondrial gene expression suggest the import of RNA to mitochondria, but also transcription from the mitochondrial genome. Interestingly, ncRNAs have been also shown to hide small open reading frames (sORFs) encoding for small functional peptides termed micropeptides, with several examples reported with a role in mitochondria. In this review, we provide a literature overview on ncRNAs and micropeptides found to be associated with mitochondrial biology in the context of both health and disease. Although reported, small study overlap and rare replications by other groups make the presence, transport, and role of ncRNA in mitochondria an attractive, but still challenging subject. Finally, we touch the topic of their potential as prognosis markers and therapeutic targets.

Abstract Background Breast cancer is the most prevalent tumor entity in Li-Fraumeni syndrome. Up to 80% of individuals with a Li-Fraumeni-like phenotype do not harbor detectable causative germline TP53 variants. Yet, no systematic panel analyses for a wide range of cancer predisposition genes have been conducted on cohorts of women with breast cancer fulfilling Li-Fraumeni(-like) clinical diagnostic criteria. Methods To specifically help explain the diagnostic gap of TP53 wild-type Li-Fraumeni(-like) breast cancer cases, we performed array-based CGH (comparative genomic hybridization) and panel-based sequencing of 94 cancer predisposition genes on 83 breast cancer patients suggestive of Li-Fraumeni syndrome who had previously had negative test results for causative BRCA1, BRCA2, and TP53 germline variants. Results We identified 13 pathogenic or likely pathogenic germline variants in ten patients and in nine genes, including four copy number aberrations and nine single-nucleotide variants or small indels. Three patients presented as double-mutation carriers involving two different genes each. In five patients (5 of 83; 6% of cohort), we detected causative pathogenic variants in established hereditary breast cancer susceptibility genes (i.e., PALB2, CHEK2, ATM). Five further patients (5 of 83; 6% of cohort) were found to harbor pathogenic variants in genes lacking a firm association with breast cancer susceptibility to date (i.e., Fanconi pathway genes, RECQ family genes, CDKN2A/p14ARF, and RUNX1). Conclusions Our study details the mutational spectrum in breast cancer patients suggestive of Li-Fraumeni syndrome and indicates the need for intensified research on monoallelic variants in Fanconi pathway and RECQ family genes. Notably, this study further reveals a large portion of still unexplained Li-Fraumeni(-like) cases, warranting comprehensive investigation of recently described candidate genes as well as noncoding regions of the TP53 gene in patients with Li-Fraumeni(-like) syndrome lacking TP53 variants in coding regions.

Abstract Background Mitochondrial acyl-CoA dehydrogenase family member 9 (ACAD9) is essential for the assembly of mitochondrial respiratory chain complex I. Disease causing biallelic variants in ACAD9 have been reported in individuals presenting with lactic acidosis and cardiomyopathy. Results We describe the genetic, clinical and biochemical findings in a cohort of 70 patients, of whom 29 previously unpublished. We found 34 known and 18 previously unreported variants in ACAD9. No patients harbored biallelic loss of function mutations, indicating that this combination is unlikely to be compatible with life. Causal pathogenic variants were distributed throughout the entire gene, and there was no obvious genotype-phenotype correlation. Most of the patients presented in the first year of life. For this subgroup the survival was poor (50% not surviving the first 2 years) comparing to patients with a later presentation (more than 90% surviving 10 years). The most common clinical findings were cardiomyopathy (85%), muscular weakness (75%) and exercise intolerance (72%). Interestingly, severe intellectual deficits were only reported in one patient and severe developmental delays in four patients. More than 70% of the patients were able to perform the same activities of daily living when compared to peers. Conclusions Our data show that riboflavin treatment improves complex I activity in the majority of patient-derived fibroblasts tested. This effect was also reported for most of the treated patients and is mirrored in the survival data. In the patient group with disease-onset below 1 year of age, we observed a statistically-significant better survival for patients treated with riboflavin.

Biallelic pathogenic variants in CLPP, encoding mitochondrial matrix peptidase ClpP, cause a rare autosomal recessive condition, Perrault syndrome type 3 (PRLTS3). It is characterized by primary ovarian insufficiency and early sensorineural hearing loss, often associated with progressive neurological deficits. Mouse models showed that accumulations of (i) its main protein interactor, the substrate-selecting AAA+ ATPase ClpX, (ii) mitoribosomes, and (iii) mtDNA nucleoids are the main cellular consequences of ClpP absence. However, the sequence of these events and their validity in human remain unclear. Here, we studied global proteome profiles to define ClpP substrates among mitochondrial ClpX interactors, which accumulated consistently in ClpP-null mouse embryonal fibroblasts and brains. Validation work included novel ClpP-mutant patient fibroblast proteomics. ClpX co-accumulated in mitochondria with the nucleoid component POLDIP2, the mitochondrial poly(A) mRNA granule element LRPPRC, and tRNA processing factor GFM1 (in mouse, also GRSF1). Only in mouse did accumulated ClpX, GFM1, and GRSF1 appear in nuclear fractions. Mitoribosomal accumulation was minor. Consistent accumulations in murine and human fibroblasts also affected multimerizing factors not known as ClpX interactors, namely, OAT, ASS1, ACADVL, STOM, PRDX3, PC, MUT, ALDH2, PMPCB, UQCRC2, and ACADSB, but the impact on downstream metabolites was marginal. Our data demonstrate the primary impact of ClpXP on the assembly of proteins with nucleic acids and show nucleoid enlargement in human as a key consequence.

Genome sequencing alone fails to provide a genetic diagnosis for many Mendelian disorder patients. Here, the authors utilize RNA sequencing to complement genotyping of patients with a rare mitochondrial disease by detecting aberrant RNA expression, splicing and allele-specific expression.

Abstract Coenzyme Q (CoQ) is a key component of the mitochondrial respiratory chain, but it also has several other functions in the cellular metabolism. One of them is to function as an electron carrier in the reaction catalyzed by sulfide:quinone oxidoreductase (SQR), which catalyzes the first reaction in the hydrogen sulfide oxidation pathway. Therefore, SQR may be affected by CoQ deficiency. Using human skin fibroblasts and two mouse models with primary CoQ deficiency, we demonstrate that severe CoQ deficiency causes a reduction in SQR levels and activity, which leads to an alteration of mitochondrial sulfide metabolism. In cerebrum of Coq9R239X mice, the deficit in SQR induces an increase in thiosulfate sulfurtransferase and sulfite oxidase, as well as modifications in the levels of thiols. As a result, biosynthetic pathways of glutamate, serotonin, and catecholamines were altered in the cerebrum, and the blood pressure was reduced. Therefore, this study reveals the reduction in SQR activity as one of the pathomechanisms associated with CoQ deficiency syndrome.

The post-transcriptional 5-methylcytosine (m5C) modification occurs in a wide range of nuclear-encoded RNAs. Here the authors identify the mitochondrial tRNA-Met as a target for the m5C methyltransferase NSun3—found mutated in a mitochondrial disease patient—and link mitochondrial tRNA modifications with energy metabolism.

Karin Tuschl, Philippa Mills and colleagues report mutations in the manganese (Mn) transporter gene SLC39A14in childhood-onset parkinsonism-dystonia. Using functional recapitulation, the authors also show that slc39A14 loss-of-function in zebrafish can lead to Mn dysregulation and locomotor impairment.

The accurate quantification of cellular and mitochondrial bioenergetic activity is of great interest in medicine and biology. Mitochondrial stress tests performed with Seahorse Bioscience XF Analyzers allow the estimation of different bioenergetic measures by monitoring the oxygen consumption rates (OCR) of living cells in multi-well plates. However, studies of the statistical best practices for determining aggregated OCR measurements and comparisons have been lacking. Therefore, to understand how OCR behaves across different biological samples, wells, and plates, we performed mitochondrial stress tests in 126 96-well plates involving 203 fibroblast cell lines. We show that the noise of OCR is multiplicative, that outlier data points can concern individual measurements or all measurements of a well, and that the inter-plate variation is greater than the intra-plate variation. Based on these insights, we developed a novel statistical method, OCR-Stats, that: i) robustly estimates OCR levels modeling multiplicative noise and automatically identifying outlier data points and outlier wells; and ii) performs statistical testing between samples, taking into account the different magnitudes of the between- and within-plate variations. This led to a significant reduction of the coefficient of variation across plates of basal respiration by 45% and of maximal respiration by 29%. Moreover, using positive and negative controls, we show that our statistical test outperforms the existing methods, which suffer from an excess of either false positives (within-plate methods), or false negatives (between-plate methods). Altogether, this study provides statistical good practices to support experimentalists in designing, analyzing, testing, and reporting the results of mitochondrial stress tests using this high throughput platform.

Editorial summary The advent of mitochondrial replacement techniques poses many scientific, regulatory, and ethical questions. Previous studies suggest good safety and efficacy profiles of these techniques, but challenges remain for clinical implementation and international consensus is needed on the regulation of these approaches.

Epigenetic regulation has been postulated to affect glucose metabolism, insulin sensitivity and the risk of type 2 diabetes. Therefore, we performed an epigenome-wide association study for measures of glucose metabolism in whole blood samples of the population-based Cooperative Health Research in the Region of Augsburg F4 study using the Illumina HumanMethylation 450 BeadChip. We identified a total of 31 CpG sites where methylation level was associated with measures of glucose metabolism after adjustment for age, sex, smoking, and estimated white blood cell proportions and correction for multiple testing using the Benjamini-Hochberg (B-H) method (four for fasting glucose, seven for fasting insulin, 25 for homeostasis model assessment-insulin resistance [HOMA-IR]; B-H-adjusted p-values between 9.2x10(-5) and 0.047). In addition, DNA methylation at cg06500161 (annotated to ABCG1) was associated with all the aforementioned phenotypes and 2-hour glucose (B-H-adjusted p-values between 9.2x10(-5) and 3.0x10(-3)). Methylation status of additional three CpG sites showed an association with fasting insulin only after additional adjustment for body mass index (BMI) (B-H-adjusted p-values = 0.047). Overall, effect strengths were reduced by around 30% after additional adjustment for BMI, suggesting that this variable has an influence on the investigated phenotypes. Furthermore, we found significant associations between methylation status of 21 of the aforementioned CpG sites and 2-hour insulin in a subset of samples with seven significant associations persisting after additional adjustment for BMI. In a subset of 533 participants, methylation of the CpG site cg06500161 (ABCG1) was inversely associated with ABCG1 gene expression (B-H-adjusted p-value = 1.5x10(-9)). Additionally, we observed an enrichment of the top 1,000 CpG sites for diabetes-related canonical pathways using Ingenuity Pathway Analysis. In conclusion, our study indicates that DNA methylation and diabetes-related traits are associated and that these associations are partially BMI-dependent. Furthermore, the interaction of ABCG1 with glucose metabolism is modulated by epigenetic processes.

Elevated levels of C-reactive protein (CRP, determined by a high-sensitivity assay) indicate low-grade inflammation which is implicated in many age-related disorders. Epigenetic studies on CRP might discover molecular mechanisms underlying CRP regulation. We aimed to identify DNA methylation sites related to CRP concentrations in cells and tissues regulating low-grade inflammation.Genome-wide DNA methylation was measured in peripheral blood in 1,741 participants of the KORA F4 study using Illumina HumanMethylation450 BeadChip arrays. Four CpG sites (located at BCL3, AQP3, SOCS3, and cg19821297 intergenic at chromosome 19p13.2, P ≤ 1.01E-07) were significantly hypomethylated at high CRP concentrations independent of various confounders including age, sex, BMI, smoking, and white blood cell composition. Findings were not sex-specific. CRP-related top genes were enriched in JAK/STAT pathways (Benjamini-Hochberg corrected P < 0.05). Results were followed-up in three studies using DNA from peripheral blood (EPICOR, n = 503) and adipose tissue (TwinsUK, n = 368) measured as described above and from liver tissue (LMU liver cohort, n = 286) measured by MALDI-TOF mass spectrometry using EpiTYPER. CpG sites at the AQP3 locus (significant p-values in peripheral blood = 1.72E-03 and liver tissue = 1.51E-03) and the SOCS3 locus (p-values in liver < 2.82E-05) were associated with CRP in the validation panels.Epigenetic modifications seem to engage in low-grade inflammation, possibly via JAK/STAT mediated pathways. Results suggest a shared relevance across different tissues at the AQP3 locus and highlight a role of DNA methylation for CRP regulation at the SOCS3 locus.

Abstract Type 2 diabetes (T2D) can be prevented in pre‐diabetic individuals with impaired glucose tolerance (IGT). Here, we have used a metabolomics approach to identify candidate biomarkers of pre‐diabetes. We quantified 140 metabolites for 4297 fasting serum samples in the population‐based Cooperative Health Research in the Region of Augsburg (KORA) cohort. Our study revealed significant metabolic variation in pre‐diabetic individuals that are distinct from known diabetes risk indicators, such as glycosylated hemoglobin levels, fasting glucose and insulin. We identified three metabolites (glycine, lysophosphatidylcholine (LPC) (18:2) and acetylcarnitine) that had significantly altered levels in IGT individuals as compared to those with normal glucose tolerance, with P‐values ranging from 2.4 × 10−4 to 2.1 × 10−13. Lower levels of glycine and LPC were found to be predictors not only for IGT but also for T2D, and were independently confirmed in the European Prospective Investigation into Cancer and Nutrition (EPIC)‐Potsdam cohort. Using metabolite–protein network analysis, we identified seven T2D‐related genes that are associated with these three IGT‐specific metabolites by multiple interactions with four enzymes. The expression levels of these enzymes correlate with changes in the metabolite concentrations linked to diabetes. Our results may help developing novel strategies to prevent T2D.

Biological systems consist of multiple organizational levels all densely interacting with each other to ensure function and flexibility of the system. Simultaneous analysis of cross-sectional multi-omics data from large population studies is a powerful tool to comprehensively characterize the underlying molecular mechanisms on a physiological scale. In this study, we systematically analyzed the relationship between fasting serum metabolomics and whole blood transcriptomics data from 712 individuals of the German KORA F4 cohort. Correlation-based analysis identified 1,109 significant associations between 522 transcripts and 114 metabolites summarized in an integrated network, the 'human blood metabolome-transcriptome interface' (BMTI). Bidirectional causality analysis using Mendelian randomization did not yield any statistically significant causal associations between transcripts and metabolites. A knowledge-based interpretation and integration with a genome-scale human metabolic reconstruction revealed systematic signatures of signaling, transport and metabolic processes, i.e. metabolic reactions mainly belonging to lipid, energy and amino acid metabolism. Moreover, the construction of a network based on functional categories illustrated the cross-talk between the biological layers at a pathway level. Using a transcription factor binding site enrichment analysis, this pathway cross-talk was further confirmed at a regulatory level. Finally, we demonstrated how the constructed networks can be used to gain novel insights into molecular mechanisms associated to intermediate clinical traits. Overall, our results demonstrate the utility of a multi-omics integrative approach to understand the molecular mechanisms underlying both normal physiology and disease.

The functional consequences of trait associated SNPs are often investigated using expression quantitative trait locus (eQTL) mapping. While trait-associated variants may operate in a cell-type specific manner, eQTL datasets for such cell-types may not always be available. We performed a genome-environment interaction (GxE) meta-analysis on data from 5,683 samples to infer the cell type specificity of whole blood cis-eQTLs. We demonstrate that this method is able to predict neutrophil and lymphocyte specific cis-eQTLs and replicate these predictions in independent cell-type specific datasets. Finally, we show that SNPs associated with Crohn's disease preferentially affect gene expression within neutrophils, including the archetypal NOD2 locus.

Genome-wide association studies (GWAS) have uncovered numerous genetic variants (SNPs) that are associated with blood pressure (BP). Genetic variants may lead to BP changes by acting on intermediate molecular phenotypes such as coded protein sequence or gene expression, which in turn affect BP variability. Therefore, characterizing genes whose expression is associated with BP may reveal cellular processes involved in BP regulation and uncover how transcripts mediate genetic and environmental effects on BP variability. A meta-analysis of results from six studies of global gene expression profiles of BP and hypertension in whole blood was performed in 7017 individuals who were not receiving antihypertensive drug treatment. We identified 34 genes that were differentially expressed in relation to BP (Bonferroni-corrected p

It has been suggested that mitochondrial dysfunction has an influence on lipid metabolism. The fact that mitochondrial defects can be accumulated over time as a normal part of aging may explain why cholesterol levels often are altered with age. To test the hypothesis whether mitochondrial variants are associated with lipid profile (total cholesterol, LDL, HDL, and triglycerides) we analyzed a total number of 978 mitochondrial single nucleotide polymorphisms (mtSNPs) in a sample of 2,815 individuals participating in the population-based KORA F4 study. To assess mtSNP association while taking the presence of heteroplasmy into account we used the raw signal intensity values measured on the microarray and applied linear regression. Ten mtSNPs (mt3285, mt3336, mt5285, mt6591, mt6671, mt9163, mt13855, mt13958, mt14000, and mt14580) were significantly associated with HDL cholesterol and one mtSNP (mt15074) with triglycerides levels. These results highlight the importance of the mitochondrial genome among the factors that contribute to the regulation of lipid levels. Focusing on mitochondrial variants may lead to further insights regarding the underlying physiological mechanisms, or even to the development of innovative treatments. Since this is the first mitochondrial genome-wide association analysis (mtGWAS) for lipid profile, further analyses are needed to follow up on the present findings.

Autoimmune thyroid diseases (AITD) are common, affecting 2-5% of the general population. Individuals with positive thyroid peroxidase antibodies (TPOAbs) have an increased risk of autoimmune hypothyroidism (Hashimoto's thyroiditis), as well as autoimmune hyperthyroidism (Graves' disease). As the possible causative genes of TPOAbs and AITD remain largely unknown, we performed GWAS meta-analyses in 18,297 individuals for TPOAb-positivity (1769 TPOAb-positives and 16,528 TPOAb-negatives) and in 12,353 individuals for TPOAb serum levels, with replication in 8,990 individuals. Significant associations (P

Recently, mutations in the mitochondrial translation optimization factor 1 gene (MTO1) were identified as causative in children with hypertrophic cardiomyopathy, lactic acidosis and respiratory chain defect. Here, we describe an MTO1-deficient mouse model generated by gene trap mutagenesis that mirrors the human phenotype remarkably well. As in patients, the most prominent signs and symptoms were cardiovascular and included bradycardia and cardiomyopathy. In addition, the mutant mice showed a marked worsening of arrhythmias during induction and reversal of anaesthesia. The detailed morphological and biochemical workup of murine hearts indicated that the myocardial damage was due to complex I deficiency and mitochondrial dysfunction. In contrast, neurological examination was largely normal in Mto1-deficient mice. A translational consequence of this mouse model may be to caution against anaesthesia-related cardiac arrhythmias which may be fatal in patients.

Restless legs syndrome (RLS) is a common neurologic disorder characterized by nightly dysesthesias affecting the legs primarily during periods of rest and relieved by movement. RLS is a complex genetic disease and susceptibility factors in six genomic regions have been identified by means of genome-wide association studies (GWAS). For some complex genetic traits, expression quantitative trait loci (eQTLs) are enriched among trait-associated single nucleotide polymorphisms (SNPs). With the aim of identifying new genetic susceptibility factors for RLS, we assessed the 332 best-associated SNPs from the genome-wide phase of the to date largest RLS GWAS for cis-eQTL effects in peripheral blood from individuals of European descent. In 740 individuals belonging to the KORA general population cohort, 52 cis-eQTLs with pnominal

BackgroundWe aimed to assess whether whole blood expression quantitative trait loci (eQTLs) with effects in cis and trans are robust and can be used to identify regulatory pathways affecting disease susceptibility.Materials and methodsWe performed whole-genome eQTL analyses in 890 participants of the KORA F4 study and in two independent replication samples (SHIP-TREND, N = 976 and EGCUT, N = 842) using linear regression models and Bonferroni correction.ResultsIn the KORA F4 study, 4,116 cis-eQTLs (defined as SNP-probe pairs where the SNP is located within a 500 kb window around the transcription unit) and 94 trans-eQTLs reached genome-wide significance and overall 91% (92% of cis-, 84% of trans-eQTLs) were confirmed in at least one of the two replication studies. Different study designs including distinct laboratory reagents (PAXgene™ vs. Tempus™ tubes) did not affect reproducibility (separate overall replication overlap: 78% and 82%). Immune response pathways were enriched in cis- and trans-eQTLs and significant cis-eQTLs were partly coexistent in other tissues (cross-tissue similarity 40-70%). Furthermore, four chromosomal regions displayed simultaneous impact on multiple gene expression levels in trans, and 746 eQTL-SNPs have been previously reported to have clinical relevance. We demonstrated cross-associations between eQTL-SNPs, gene expression levels in trans, and clinical phenotypes as well as a link between eQTLs and human metabolic traits via modification of gene regulation in cis.ConclusionsOur data suggest that whole blood is a robust tissue for eQTL analysis and may be used both for biomarker studies and to enhance our understanding of molecular mechanisms underlying gene-disease associations.