Mitochondrial disorders, characterized by clinical symptoms and/or OXPHOS deficiencies, are caused by pathogenic variants in mitochondrial genes. However, pathogenic variants in some of these genes can lead to clinical manifestations which overlap with other neuromuscular diseases, which can be caused by pathogenic variants in non-mitochondrial genes as well. Mitochondrial pathogenic variants can be found in the mitochondrial DNA (mtDNA) or in any of the 1,500 nuclear genes with a mitochondrial function. We have performed a two-step next-generation sequencing approach in a cohort of 117 patients, mostly children, in whom a mitochondrial disease-cause could likely or possibly explain the phenotype. A total of 86 patients had a mitochondrial disorder, according to established clinical and biochemical criteria. The other 31 patients had neuromuscular symptoms, where in a minority a mitochondrial genetic cause is present, but a non-mitochondrial genetic cause is more likely. All patients were screened for pathogenic variants in the mtDNA and, if excluded, analyzed by whole exome sequencing (WES). Variants were filtered for being pathogenic and compatible with an autosomal or X-linked recessive mode of inheritance in families with multiple affected siblings and/or consanguineous parents. Non-consanguineous families with a single patient were additionally screened for autosomal and X-linked dominant mutations in a predefined gene-set. We identified causative pathogenic variants in the mtDNA in 20% of the patient-cohort, and in nuclear genes in 49%, implying an overall yield of 68%. We identified pathogenic variants in mitochondrial and non-mitochondrial genes in both groups with, obviously, a higher number of mitochondrial genes affected in mitochondrial disease patients. Furthermore, we show that 31% of the disease-causing genes in the mitochondrial patient group were not included in the MitoCarta database, and therefore would have been missed with MitoCarta based gene-panels. We conclude that WES is preferable to panel-based approaches for both groups of patients, as the mitochondrial gene-list is not complete and mitochondrial symptoms can be secondary. Also, clinically and genetically heterogeneous disorders would require sequential use of multiple different gene panels. We conclude that WES is a comprehensive and unbiased approach to establish a genetic diagnosis in these patients, able to resolve multi-genic disease-causes.

We studied the mtDNA bottleneck in zebrafish to elucidate size, timing, and variation in germline and non-germline cells. Mature zebrafish oocytes contain, on average, 19.0 × 10(6) mtDNA molecules with high variation between oocytes. During embryogenesis, the mtDNA copy number decreases to ∼170 mtDNA molecules per primordial germ cell (PGC), a number similar to that in mammals, and to ∼50 per non-PGC. These occur at the same developmental stage, implying considerable variation in mtDNA copy number in (non-)PGCs of the same female, dictated by variation in the mature oocyte. The presence of oocytes with low mtDNA numbers, if similar in humans, could explain how (de novo) mutations can reach high mutation loads within a single generation. High mtDNA copy numbers in mature oocytes are established by mtDNA replication during oocyte development. Bottleneck differences between germline and non-germline cells, due to early differentiation of PGCs, may account for different distribution patterns of familial mutations.

In a 51-year-old patient of consanguineous parents with a severe neuromuscular phenotype of early-onset ataxia, myoclonia, dysarthria, muscle weakness and exercise intolerance, exome sequencing revealed a novel homozygous variant (c.-264_31delinsCTCACAAATGCTCA) in the mitochondrial FAD-transporter gene SLC25A32. Flavin adenine dinucleotide (FAD) is an essential co-factor for many mitochondrial enzymes and impaired mitochondrial FAD-transport was supported by a reduced oxidative phosphorylation complex II activity in the patient's muscle, decreased ATP production in fibroblasts, and a deficiency of mitochondrial FAD-dependent enzymes. Clinically, the patient showed improvement upon riboflavin treatment, which is a precursor of FAD. Our results confirm the recently reported case of SLC25A32 as a cause of riboflavin-responsive disease. Our patient showed a more severe clinical phenotype compared with the reported patient, corresponding with the (most likely) complete absence of the SLC25A32-encoding MFT (Mitochondrial Folate Transporter) protein.

Hirschsprung disease (HSCR) is the most frequent developmental anomaly of the enteric nervous system with an incidence of 1/5000 live births. Chronic intestinal pseudo-obstruction (CIPO) is less frequent and classified as neurogenic or myogenic. Isolated HSCR has an oligogenic inheritance with RET as the major disease-causing gene, while CIPO is genetically heterogeneous, caused by mutations in smooth muscle-specific genes. Here, we describe a series of patients with developmental disorders including gastrointestinal dysmotility, and investigate the underlying molecular bases. Trio-exome sequencing led to the identification of biallelic variants in ERBB3 and ERBB2 in eight individuals variably associating HSCR, CIPO, peripheral neuropathy and arthrogryposis. Thorough gut histology revealed aganglionosis, hypoganglionosis and intestinal smooth muscle abnormalities. The cell-type-specific ErbB3 and ErbB2 function was further analysed in mouse single-cell RNA sequencing data and in a conditional ErbB3-deficient mouse model, revealing a primary role for ERBB3 in enteric progenitors. The consequences of the identified variants were evaluated using RT-qPCR on patient-derived fibroblasts or immunoblot assays on Neuro-2a cells overexpressing either wild-type or mutant proteins, revealing either decreased expression or altered phosphorylation of the mutant receptors. Our results demonstrate that dysregulation of ERBB3 or ERBB2 leads to a broad spectrum of developmental anomalies including intestinal dysmotility.

In establishing a genetic diagnosis in heterogeneous neurological disease, clinical characterization and whole exome sequencing (WES) go hand-in-hand. Clinical data are essential, not only to guide WES variant selection and define the clinical severity of a genetic defect but also to identify other patients with defects in the same gene. In an infant patient with sensorineural hearing loss, psychomotor retardation, and epilepsy, WES resulted in identification of a novel homozygous CLPP frameshift mutation (c.21delA). Based on the gene defect and clinical symptoms, the diagnosis Perrault syndrome type 3 (PRLTS3) was established. The patient’s brain-MRI revealed specific abnormalities of the subcortical and deep cerebral white matter and the middle blade of the corpus callosum, which was used to identify similar patients in the Amsterdam brain-MRI database, containing over 3000 unclassified leukoencephalopathy cases. In three unrelated patients with similar MRI abnormalities the CLPP gene was sequenced, and in two of them novel missense mutations were identified together with a large deletion that covered part of the CLPP gene on the other allele. The severe neurological and MRI abnormalities in these young patients were due to the drastic impact of the CLPP mutations, correlating with the variation in clinical manifestations among previously reported patients. Our data show that similarity in brain-MRI patterns can be used to identify novel PRLTS3 patients, especially during early disease stages, when only part of the disease manifestations are present. This seems especially applicable to the severely affected cases in which CLPP function is drastically affected and MRI abnormalities are pronounced.