Children born after assisted reproduction more commonly carry a mitochondrial genotype associating with low birthweight

Children conceived through assisted reproductive technologies (ART) have an elevated risk of lower birthweight, yet the underlying cause remains unclear. Our study explores mitochondrial DNA (mtDNA) variants as contributors to birthweight differences by impacting mitochondrial function during prenatal development. We deep-sequenced the mtDNA of 451 ART and spontaneously conceived (SC) individuals, 157 mother-child pairs and 113 individual oocytes from either natural menstrual cycles or after ovarian stimulation (OS) and find that ART individuals carried a different mtDNA genotype than SC individuals, with more de novo non-synonymous variants. These variants, along with rRNA variants, correlate with lower birthweight percentiles, independent of conception mode. Their higher occurrence in ART individuals stems from de novo mutagenesis associated with maternal aging and OS-induced oocyte cohort size. Future research will establish the long-term health consequences of these changes and how these findings will impact the clinical practice and patient counselling in the future.

Low mitochondrial DNA copy number in buffy coat DNA of primary open-angle glaucoma patients

Primary open-angle glaucoma (POAG) is characterized by optic nerve degeneration and irreversible loss of retinal ganglion cells (RGCs). The pathophysiology is not fully understood. Since RGCs have a high energy demand, suboptimal mitochondrial function may put the survival of these neurons at risk. In the present study, we explored whether mtDNA copy number or mtDNA deletions could reveal a mitochondrial component in POAG pathophysiology. Buffy coat DNA was isolated from EDTA blood of age- and sex-matched study groups, namely POAG patients with high intraocular pressure (IOP) at diagnosis (high tension glaucoma: HTG; n = 97), normal tension glaucoma patients (NTG, n = 37), ocular hypertensive controls (n = 9), and cataract controls (without glaucoma; n = 32), all without remarkable comorbidities. The number of mtDNA copies was assessed through qPCR quantification of the mitochondrial D-loop and nuclear B2M gene. Presence of the common 4977 base pair mtDNA deletion was assessed by a highly sensitive breakpoint PCR. Analysis showed that HTG patients had a lower number of mtDNA copies per nuclear DNA than NTG patients (p-value <0.01, Dunn test) and controls (p-value <0.001, Dunn test). The common 4977 base pair mtDNA deletion was not detected in any of the participants. A lower mtDNA copy number in blood of HTG patients suggests a role for a genetically defined, deficient mtDNA replication in the pathology of HTG. This may cause a low number of mtDNA copies in RGCs, which together with aging and high IOP, may lead to mitochondrial dysfunction, and contribute to glaucoma pathology.

Genetic Profiling of Sodium Channels in Diabetic Painful and Painless and Idiopathic Painful and Painless Neuropathies

Neuropathic pain is a frequent feature of diabetic peripheral neuropathy (DPN) and small fiber neuropathy (SFN). Resolving the genetic architecture of these painful neuropathies will lead to better disease management strategies, counselling and intervention. Our aims were to profile ten sodium channel genes (SCG) expressed in a nociceptive pathway in painful and painless DPN and painful and painless SFN patients, and to provide a perspective for clinicians who assess patients with painful peripheral neuropathy. Between June 2014 and September 2016, 1125 patients with painful-DPN (n = 237), painless-DPN (n = 309), painful-SFN (n = 547) and painless-SFN (n = 32), recruited in four different centers, were analyzed for SCN3A, SCN7A-SCN11A and SCN1B-SCN4B variants by single molecule Molecular inversion probes-Next Generation Sequence. Patients were grouped based on phenotype and the presence of SCG variants. Screening of SCN3A, SCN7A-SCN11A, and SCN1B-SCN4B revealed 125 different (potential) pathogenic variants in 194 patients (17.2%, n = 194/1125). A potential pathogenic variant was present in 18.1% (n = 142/784) of painful neuropathy patients vs. 15.2% (n = 52/341) of painless neuropathy patients (17.3% (n = 41/237) for painful-DPN patients, 14.9% (n = 46/309) for painless-DPN patients, 18.5% (n = 101/547) for painful-SFN patients, and 18.8% (n = 6/32) for painless-SFN patients). Of the variants detected, 70% were in SCN7A, SCN9A, SCN10A and SCN11A. The frequency of SCN9A and SCN11A variants was the highest in painful-SFN patients, SCN7A variants in painful-DPN patients, and SCN10A variants in painless-DPN patients. Our findings suggest that rare SCG genetic variants may contribute to the development of painful neuropathy. Genetic profiling and SCG variant identification should aid in a better understanding of the genetic variability in patients with painful and painless neuropathy, and may lead to better risk stratification and the development of more targeted and personalized pain treatments.

Peripheral Ion Channel Genes Screening in Painful Small Fiber Neuropathy

Neuropathic pain is a characteristic feature of small fiber neuropathy (SFN), which in 18% of the cases is caused by genetic variants in voltage-gated sodium ion channels. In this study, we assessed the role of fifteen other ion channels in neuropathic pain. Patients with SFN (n = 414) were analyzed for ANO1, ANO3, HCN1, KCNA2, KCNA4, KCNK18, KCNN1, KCNQ3, KCNQ5, KCNS1, TRPA1, TRPM8, TRPV1, TRPV3 and TRPV4 variants by single-molecule molecular inversion probes-next-generation sequencing. These patients did not have genetic variants in SCN3A, SCN7A-SCN11A and SCN1B-SCN4B. In twenty patients (20/414, 4.8%), a potentially pathogenic heterozygous variant was identified in an ion-channel gene (ICG). Variants were present in seven genes, for two patients (0.5%) in ANO3, one (0.2%) in KCNK18, two (0.5%) in KCNQ3, seven (1.7%) in TRPA1, three (0.7%) in TRPM8, three (0.7%) in TRPV1 and two (0.5%) in TRPV3. Variants in the TRP genes were the most frequent (n = 15, 3.6%), partly in patients with high mean maximal pain scores VAS = 9.65 ± 0.7 (n = 4). Patients with ICG variants reported more severe pain compared to patients without such variants (VAS = 9.36 ± 0.72 vs. VAS = 7.47 ± 2.37). This cohort study identified ICG variants in neuropathic pain in SFN, complementing previous findings of ICG variants in diabetic neuropathy. These data show that ICG variants are central in neuropathic pain of different etiologies and provides promising gene candidates for future research.

Lead-exposure associated miRNAs in humans and Alzheimer’s disease as potential biomarkers of the disease and disease processes

Alzheimer’s disease (AD) is a neurodegenerative disease that eventually affects memory and behavior. The identification of biomarkers based on risk factors for AD provides insight into the disease since the exact cause of AD remains unknown. Several studies have proposed microRNAs (miRNAs) in blood as potential biomarkers for AD. Exposure to heavy metals is a potential risk factor for onset and development of AD. Blood cells of subjects that are exposed to lead detected in the circulatory system, potentially reflect molecular responses to this exposure that are similar to the response of neurons. In this study we analyzed blood cell-derived miRNAs derived from a general population as proxies of potentially AD-related mechanisms triggered by lead exposure. Subsequently, we analyzed these mechanisms in the brain tissue of AD subjects and controls. A total of four miRNAs were identified as lead exposure-associated with hsa-miR-3651, hsa-miR-150-5p and hsa-miR-664b-3p being negatively and hsa-miR-627 positively associated. In human brain derived from AD and AD control subjects all four miRNAs were detected. Moreover, two miRNAs (miR-3651, miR-664b-3p) showed significant differential expression in AD brains versus controls, in accordance with the change direction of lead exposure. The miRNAs’ gene targets were validated for expression in the human brain and were found enriched in AD-relevant pathways such as axon guidance. Moreover, we identified several AD relevant transcription factors such as CREB1 associated with the identified miRNAs. These findings suggest that the identified miRNAs are involved in the development of AD and might be useful in the development of new, less invasive biomarkers for monitoring of novel therapies or of processes involved in AD development.

Peripheral Ion Channel Gene Screening in Painful- and Painless-Diabetic Neuropathy

Neuropathic pain is common in diabetic peripheral neuropathy (DN), probably caused by pathogenic ion channel gene variants. Therefore, we performed molecular inversion probes-next generation sequencing of 5 transient receptor potential cation channels, 8 potassium channels and 2 calcium-activated chloride channel genes in 222 painful- and 304 painless-DN patients. Twelve painful-DN (5.4%) patients showed potentially pathogenic variants (five nonsense/frameshift, seven missense, one out-of-frame deletion) in ANO3 (n = 3), HCN1 (n = 1), KCNK18 (n = 2), TRPA1 (n = 3), TRPM8 (n = 3) and TRPV4 (n = 1) and fourteen painless-DN patients (4.6%-three nonsense/frameshift, nine missense, one out-of-frame deletion) in ANO1 (n = 1), KCNK18 (n = 3), KCNQ3 (n = 1), TRPA1 (n = 2), TRPM8 (n = 1), TRPV1 (n = 3) and TRPV4 (n = 3). Missense variants were present in both conditions, presumably with loss- or gain-of-functions. KCNK18 nonsense/frameshift variants were found in painless/painful-DN, making a causal role in pain less likely. Surprisingly, premature stop-codons with likely nonsense-mediated RNA-decay were more frequent in painful-DN. Although limited in number, painful-DN patients with ion channel gene variants reported higher maximal pain during the night and day. Moreover, painful-DN patients with TRP variants had abnormal thermal thresholds and more severe pain during the night and day. Our results suggest a role of ion channel gene variants in neuropathic pain, but functional validation is required.

RRM1 variants cause a mitochondrial DNA maintenance disorder via impaired de novo nucleotide synthesis

Mitochondrial DNA (mtDNA) depletion/deletions syndromes (MDDS) encompass a clinically and etiologically heterogenous group of mitochondrial disorders due to impaired mtDNA maintenance. Among the most frequent causes of MDDS are defects in nucleoside/nucleotide metabolism, which is critical for synthesis and homeostasis of the deoxynucleoside triphosphate (dNTP) substrates of mtDNA replication. A central enzyme for generating dNTPs is ribonucleotide reductase, a critical mediator of de novo nucleotide synthesis composed of catalytic RRM1 subunits in complex with RRM2 or p53R2. Here, we report five probands from four families who presented with ptosis and ophthalmoplegia, plus other manifestations and multiple mtDNA deletions in muscle. We identified three RRM1 loss-of-function variants, including a dominant catalytic site variant (NP_001024.1: p.N427K) and two homozygous recessive variants at p.R381, which has evolutionarily conserved interactions with the specificity site. Atomistic molecular dynamics simulations indicate mechanisms by which RRM1 variants affect protein structure. Cultured primary skin fibroblasts of probands manifested mtDNA depletion under cycling conditions, indicating impaired de novo nucleotide synthesis. Fibroblasts also exhibited aberrant nucleoside diphosphate and dNTP pools and mtDNA ribonucleotide incorporation. Our data reveal primary RRM1 deficiency and, by extension, impaired de novo nucleotide synthesis are causes of MDDS.

Lacosamide Inhibition of NaV1.7 Channels Depends on its Interaction With the Voltage Sensor Domain and the Channel Pore

Lacosamide, developed as an anti-epileptic drug, has been used for the treatment of pain. Unlike typical anticonvulsants and local anesthetics which enhance fast-inactivation and bind within the pore of sodium channels, lacosamide enhances slow-inactivation of these channels, suggesting different binding mechanisms and mode of action. It has been reported that lacosamide's effect on Na_V1.5 is sensitive to a mutation in the local anesthetic binding site, and that it binds with slow kinetics to the fast-inactivated state of Na_V1.7. We recently showed that the Na_V1.7-W1538R mutation in the voltage-sensing domain 4 completely abolishes Na_V1.7 inhibition by clinically-achievable concentration of lacosamide. Our molecular docking analysis suggests a role for W1538 and pore residues as high affinity binding sites for lacosamide. Aryl sulfonamide sodium channel blockers are also sensitive to substitutions of the W1538 residue but not of pore residues. To elucidate the mechanism by which lacosamide exerts its effects, we used voltage-clamp recordings and show that lacosamide requires an intact local anesthetic binding site to inhibit Na_V1.7 channels. Additionally, the W1538R mutation does not abrogate local anesthetic lidocaine-induced blockade. We also show that the naturally occurring arginine in Na_V1.3 (Na_V1.3-R1560), which corresponds to Na_V1.7-W1538R, is not sufficient to explain the resistance of Na_V1.3 to clinically-relevant concentrations of lacosamide. However, the Na_V1.7-W1538R mutation conferred sensitivity to the Na_V1.3-selective aryl-sulfonamide blocker ICA-121431. Together, the W1538 residue and an intact local anesthetic site are required for lacosamide's block of Na_V1.7 at a clinically-achievable concentration. Moreover, the contribution of W1538 to lacosamide inhibitory effects appears to be isoform-specific.

OUP accepted manuscript

Mutations of the mitochondrial DNA are an important cause of inherited diseases that can severely affect the tissue’s homeostasis and integrity. The m.3243A > G mutation is the most commonly observed across mitochondrial disorders and is linked to multisystemic complications, including cognitive deficits. In line with in vitro experiments demonstrating the m.3243A > G’s negative impact on neuronal energy production and integrity, m.3243A > G patients show cerebral grey matter tissue changes. However, its impact on the most neuron dense, and therefore energy-consuming brain structure—the cerebellum—remains elusive. In this work, we used high-resolution structural and functional data acquired using 7 T MRI to characterize the neurodegenerative and functional signatures of the cerebellar cortex in m.3243A > G patients. Our results reveal altered tissue integrity within distinct clusters across the cerebellar cortex, apparent by their significantly reduced volume and longitudinal relaxation rate compared with healthy controls, indicating macroscopic atrophy and microstructural pathology. Spatial characterization reveals that these changes occur especially in regions related to the frontoparietal brain network that is involved in information processing and selective attention. In addition, based on resting-state functional MRI data, these clusters exhibit reduced functional connectivity to frontal and parietal cortical regions, especially in patients characterized by (i) a severe disease phenotype and (ii) reduced information-processing speed and attention control. Combined with our previous work, these results provide insight into the neuropathological changes and a solid base to guide longitudinal studies aimed to track disease progression.

Pathogenic SLIRP variants as a novel cause of autosomal recessive mitochondrial encephalomyopathy with complex I and IV deficiency

In a Dutch non-consanguineous patient having mitochondrial encephalomyopathy with complex I and complex IV deficiency, whole exome sequencing revealed two compound heterozygous variants in SLIRP. SLIRP gene encodes a stem-loop RNA-binding protein that regulates mitochondrial RNA expression and oxidative phosphorylation (OXPHOS). A frameshift and a deep-intronic splicing variant reduced the amount of functional wild-type SLIRP RNA to 5%. Consequently, in patient fibroblasts, MT-ND1, MT-ND6, and MT-CO1 expression was reduced. Lentiviral transduction of wild-type SLIRP cDNA in patient fibroblasts increased MT-ND1, MT-ND6, and MT-CO1 expression (2.5-7.2-fold), whereas mutant cDNAs did not. A fourfold decrease of citrate synthase versus total protein ratio in patient fibroblasts indicated that the resulting reduced mitochondrial mass caused the OXPHOS deficiency. Transduction with wild-type SLIRP cDNA led to a 2.4-fold increase of this ratio and partly restored OXPHOS activity. This confirmed causality of the SLIRP variants. In conclusion, we report SLIRP variants as a novel cause of mitochondrial encephalomyopathy with OXPHOS deficiency.

Whole exome sequencing reveals a homozygous C1QBP deletion as the cause of progressive external ophthalmoplegia and multiple mtDNA deletions

Whole exome sequencing (WES), analyzed with GENESIS and WeGET, revealed a homozygous deletion in the C1QBP gene in a patient with progressive external ophthalmoplegia (PEO) and multiple mtDNA deletions. The gene encodes the mitochondria-located complementary 1 Q subcomponent-binding protein, involved in mitochondrial homeostasis. Biallelic mutations in C1QBP cause mitochondrial cardiomyopathy and/or PEO with variable age of onset. Our patient showed only late-onset PEO-plus syndrome without overt cardiac involvement. Available data suggest that early-onset cardiomyopathy variants localize in important structural domains and PEO-plus variants in the coiled-coil region. Our patient demonstrates that C1QBP mutations should be considered in individuals with PEO with or without cardiomyopathy.

Natural history, outcome measures and trial readiness in LAMA2-related muscular dystrophy and SELENON-related myopathy in children and adults: protocol of the LAST STRONG study

BACKGROUND

SELENON (SEPN1)-related myopathy (SELENON-RM) is a rare congenital myopathy characterized by slowly progressive proximal muscle weakness, early onset spine rigidity and respiratory insufficiency. A muscular dystrophy caused by mutations in the LAMA2 gene (LAMA2-related muscular dystrophy, LAMA2-MD) has a similar clinical phenotype, with either a severe, early-onset due to complete Laminin subunit α2 deficiency (merosin-deficient congenital muscular dystrophy type 1A (MDC1A)), or a mild, childhood- or adult-onset due to partial Laminin subunit α2 deficiency. For both muscle diseases, no curative treatment options exist, yet promising preclinical studies are ongoing. Currently, there is a paucity on natural history data and appropriate clinical and functional outcome measures are needed to reach trial readiness.

METHODS

LAST STRONG is a natural history study in Dutch-speaking patients of all ages diagnosed with SELENON-RM or LAMA2-MD, starting August 2020. Patients have four visits at our hospital over a period of 1.5 year. At all visits, they undergo standardized neurological examination, hand-held dynamometry (age ≥ 5 years), functional measurements, questionnaires (patient report and/or parent proxy; age ≥ 2 years), muscle ultrasound including diaphragm, pulmonary function tests (spirometry, maximal inspiratory and expiratory pressure, sniff nasal inspiratory pressure; age ≥ 5 years), and accelerometry for 8 days (age ≥ 2 years); at visit one and three, they undergo cardiac evaluation (electrocardiogram, echocardiography; age ≥ 2 years), spine X-ray (age ≥ 2 years), dual-energy X-ray absorptiometry (DEXA-)scan (age ≥ 2 years) and full body magnetic resonance imaging (MRI) (age ≥ 10 years). All examinations are adapted to the patient's age and functional abilities. Correlation between key parameters within and between subsequent visits will be assessed.

DISCUSSION

Our study will describe the natural history of patients diagnosed with SELENON-RM or LAMA2-MD, enabling us to select relevant clinical and functional outcome measures for reaching clinical trial-readiness. Moreover, our detailed description (deep phenotyping) of the clinical features will optimize clinical management and will establish a well-characterized baseline cohort for prospective follow-up.

CONCLUSION

Our natural history study is an essential step for reaching trial readiness in SELENON-RM and LAMA2-MD.

TRIAL REGISTRATION

This study has been approved by medical ethical reviewing committee Region Arnhem-Nijmegen (NL64269.091.17, 2017-3911) and is registered at ClinicalTrial.gov ( NCT04478981 ).

Plasma GDF-15 concentration is not elevated in open-angle glaucoma

Aim

Recently, the level of growth differentiation factor 15 (GDF-15) in blood, was proposed as biomarker to detect mitochondrial dysfunction. In the current study, we evaluate this biomarker in open-angle glaucoma (OAG), as there is increasing evidence that mitochondrial dysfunction plays a role in the pathophysiology of this disease.

Methods

Plasma GDF-15 concentrations were measured with ELISA in 200 OAG patients and 61 age-matched controls (cataract without glaucoma). The OAG patient group consisted of high tension glaucoma (HTG; n = 162) and normal tension glaucoma (NTG; n = 38). Groups were compared using the Kruskal-Wallis nonparametric test with Dunn’s multiple comparison post-hoc correction. GDF-15 concentration was corrected for confounders identified with forward linear regression models.

Results

Before correcting for confounders, median plasma GDF-15 levels was significantly lower in the combined OAG group (p = 0.04), but not when analysing HTG and NTG patients separately. Forward linear regression analysis showed that age, gender, smoking and systemic hypertension were significant confounders affecting GDF-15 levels. After correction for these confounders, GDF-15 levels in OAG patients were no longer significantly different from controls. Subgroup analysis of the glaucoma patients did not show a correlation between disease severity and plasma GDF-15, but did reveal that for NTG patients, intake of dietary supplements, which potentially improve mitochondrial function, correlated with lower plasma GDF-15.

Conclusion

The present study suggests that plasma GDF-15 is not suited as biomarker of mitochondrial dysfunction in OAG patients.

Distal muscle weakness and optic atrophy without central nervous system involvement in a patient with a homozygous missense mutation in the C19ORF12-gene

Variants of the C19ORF12-gene have been described in patients with spastic paraplegia type 43 and in patients with mitochondrial membrane protein-associated neurodegeneration (MPAN), a subtype of neurodegeneration associated with brain iron accumulation (NBIA). In both subtypes optic atrophy and neuropathy have been frequently described. This case report describes a patient with bilateral optic atrophy and severe distal muscle weakness based on motor neuropathy without involvement of the central nervous system. Exome sequencing revealed a homozygous pathogenic missense variant (c.187G>C;p.Ala63Pro) of the C19ORF12-gene while iron deposits were absent on repeat MR-imaging of the brain, thus showing that peripheral neuropathy and optic neuropathy can be the sole manifestations of the C19ORF12-related disease spectrum whereby iron accumulation in the brain may be absent.

Evaluation of molecular inversion probe versus TruSeq® custom methods for targeted next-generation sequencing

Resolving the genetic architecture of painful neuropathy will lead to better disease management strategies. We aimed to develop a reliable method to re-sequence multiple genes in a large cohort of painful neuropathy patients at low cost. In this study, we compared sensitivity, specificity, targeting efficiency, performance and cost effectiveness of Molecular Inversion Probes-Next generation sequencing (MIPs-NGS) and TruSeq® Custom Amplicon-Next generation sequencing (TSCA-NGS). Capture probes were designed to target nine sodium channel genes (SCN3A, SCN8A-SCN11A, and SCN1B-SCN4B). One hundred sixty-six patients with diabetic and idiopathic neuropathy were tested by both methods, 70 patients were validated by Sanger sequencing. Sensitivity, specificity and performance of both techniques were comparable, and in agreement with Sanger sequencing. The average targeted regions coverage for MIPs-NGS was 97.3% versus 93.9% for TSCA-NGS. MIPs-NGS has a more versatile assay design and is more flexible than TSCA-NGS. The cost of MIPs-NGS is >5 times cheaper than TSCA-NGS when 500 or more samples are tested. In conclusion, MIPs-NGS is a reliable, flexible, and relatively inexpensive method to detect genetic variations in a large cohort of patients. In our centers, MIPs-NGS is currently implemented as a routine diagnostic tool for screening of sodium channel genes in painful neuropathy patients.

Differential effect of lacosamide on Nav1.7 variants from responsive and non-responsive patients with small fibre neuropathy

Small fibre neuropathy is a common pain disorder, which in many cases fails to respond to treatment with existing medications. Gain-of-function mutations of voltage-gated sodium channel Na_v1.7 underlie dorsal root ganglion neuronal hyperexcitability and pain in a subset of patients with small fibre neuropathy. Recent clinical studies have demonstrated that lacosamide, which blocks sodium channels in a use-dependent manner, attenuates pain in some patients with Na_v1.7 mutations; however, only a subgroup of these patients responded to the drug. Here, we used voltage-clamp recordings to evaluate the effects of lacosamide on five Na_v1.7 variants from patients who were responsive or non-responsive to treatment. We show that, at the clinically achievable concentration of 30 μM, lacosamide acts as a potent sodium channel inhibitor of Na_v1.7 variants carried by responsive patients, via a hyperpolarizing shift of voltage-dependence of both fast and slow inactivation and enhancement of use-dependent inhibition. By contrast, the effects of lacosamide on slow inactivation and use-dependence in Na_v1.7 variants from non-responsive patients were less robust. Importantly, we found that lacosamide selectively enhances fast inactivation only in variants from responders. Taken together, these findings begin to unravel biophysical underpinnings that contribute to responsiveness to lacosamide in patients with small fibre neuropathy carrying select Na_v1.7 variants.

Tfam Knockdown Results in Reduction of mtDNA Copy Number, OXPHOS Deficiency and Abnormalities in Zebrafish Embryos

High mitochondrial DNA (mtDNA) copy numbers are essential for oogenesis and embryogenesis and correlate with fertility of oocytes and viability of embryos. To understand the pathology and mechanisms associated with low mtDNA copy numbers, we knocked down mitochondrial transcription factor A (tfam), a regulator of mtDNA replication, during early zebrafish development. Reduction of tfam using a splice-modifying morpholino (MO) resulted in a 42 ± 17% decrease in mtDNA copy number in embryos at 4 days post fertilization. Morphant embryos displayed abnormal development of the eye, brain, heart, and muscle, as well as a 50 ± 22% decrease in ATP production. Transcriptome analysis revealed a decrease in protein-encoding transcripts from the heavy strand of the mtDNA, and down-regulation of genes involved in haem production and the metabolism of metabolites, which appear to trigger increased rRNA and tRNA synthesis in the nucleoli. However, this stress or compensatory response appears to fall short as pathology emerges and expression of genes related to eye development are severely down-regulated. Taken together, this study highlights the importance of sufficient mtDNA copies for early zebrafish development. Zebrafish is an excellent model to manipulate the mtDNA bottleneck and study its effect on embryogenesis rapidly and in large numbers of offspring.

Parental repeat length instability in myotonic dystrophy type 1 pre- and protomutations

Myotonic dystrophy type 1 (DM1) is caused by a CTG trinucleotide repeat expansion on chromosome 19q13.3. While DM1 premutation (36-50 repeats) and protomutation (51-80 repeats) allele carriers are mostly asymptomatic, offspring is at risk of inheriting expanded, symptom-associated, (CTG)n repeats of n > 80. In this study we aimed to evaluate the intergenerational instability of DM1 pre- and protomutation alleles, focussing on the influence of parental gender. One hundred and forty-six parent-child pairs (34 parental premutations, 112 protomutations) were retrospectively selected from the DM1 patient cohort of the Maastricht University Medical Center+. CTG repeat size of parents and children was determined by (triplet-primed) PCR followed by fragment length analysis and Southern blot analysis. Fifty-eight out of eighty-one (71.6%) paternal transmissions led to a (CTG)n repeat of n > 80 in offspring, compared with 15 out of 65 (23.1%) maternal transmissions (p < 0.001). Repeat length instability occurred for paternal (CTG)n repeats of n ≥ 45, while maternal instability did not occur until (CTG)n repeats reached a length of n ≥ 71. Transmission of premutations caused (CTG)n repeats of n > 80 in offspring only when paternally transmitted (two cases), while protomutations caused (CTG)n repeats of n > 80 in offspring in 71 cases, of which 56 (78.9%) were paternally transmitted. In conclusion, our data show that paternally transmitted pre- and protomutations were more unstable than maternally transmitted pre- and protomutations. For genetic counseling, this implies that males with a small DMPK mutation have a higher risk of symptomatic offspring compared with females. Consequently, we suggest addressing sex-dependent factors in genetic counseling of small-sized CTG repeat carriers.

Mutation-specific effects in germline transmission of pathogenic mtDNA variants

STUDY QUESTION

Does germline selection (besides random genetic drift) play a role during the transmission of heteroplasmic pathogenic mitochondrial DNA (mtDNA) mutations in humans?

SUMMARY ANSWER

We conclude that inheritance of mtDNA is mutation-specific and governed by a combination of random genetic drift and negative and/or positive selection.

WHAT IS KNOWN ALREADY

mtDNA inherits maternally through a genetic bottleneck, but the underlying mechanisms are largely unknown. Although random genetic drift is recognized as an important mechanism, selection mechanisms are thought to play a role as well.

STUDY DESIGN, SIZE, DURATION

We determined the mtDNA mutation loads in 160 available oocytes, zygotes, and blastomeres of five carriers of the m.3243A>G mutation, one carrier of the m.8993T>G mutation, and one carrier of the m.14487T>C mutation.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Mutation loads were determined in PGD samples using PCR assays and analysed mathematically to test for random sampling effects. In addition, a meta-analysis has been performed on mutation load transmission data in the literature to confirm the results of the PGD samples.

MAIN RESULTS AND THE ROLE OF CHANCE

By applying the Kimura distribution, which assumes random mechanisms, we found that mtDNA segregations patterns could be explained by variable bottleneck sizes among all our carriers (moment estimates ranging from 10 to 145). Marked differences in the bottleneck size would determine the probability that a carrier produces offspring with mutations markedly different than her own. We investigated whether bottleneck sizes might also be influenced by non-random mechanisms. We noted a consistent absence of high mutation loads in all our m.3243A>G carriers, indicating non-random events. To test this, we fitted a standard and a truncated Kimura distribution to the m.3243A>G segregation data. A Kimura distribution truncated at 76.5% heteroplasmy has a significantly better fit (P-value = 0.005) than the standard Kimura distribution. For the m.8993T>G mutation, we suspect a skewed mutation load distribution in the offspring. To test this hypothesis, we performed a meta-analysis on published blood mutation levels of offspring-mother (O-M) transmission for the m.3243A>G and m.8993T>G mutations. This analysis revealed some evidence that the O-M ratios for the m.8993T>G mutation are different from zero (P-value <0.001), while for the m.3243A>G mutation there was little evidence that the O-M ratios are non-zero. Lastly, for the m.14487T>G mutation, where the whole range of mutation loads was represented, we found no indications for selective events during its transmission.

LARGE SCALE DATA

All data are included in the Results section of this article.

LIMITATIONS, REASON FOR CAUTION

The availability of human material for the mutations is scarce, requiring additional samples to confirm our findings.

WIDER IMPLICATIONS OF THE FINDINGS

Our data show that non-random mechanisms are involved during mtDNA segregation. We aimed to provide the mechanisms underlying these selection events. One explanation for selection against high m.3243A>G mutation loads could be, as previously reported, a pronounced oxidative phosphorylation (OXPHOS) deficiency at high mutation loads, which prohibits oogenesis (e.g. progression through meiosis). No maximum mutation loads of the m.8993T>G mutation seem to exist, as the OXPHOS deficiency is less severe, even at levels close to 100%. In contrast, high mutation loads seem to be favoured, probably because they lead to an increased mitochondrial membrane potential (MMP), a hallmark on which healthy mitochondria are being selected. This hypothesis could provide a possible explanation for the skewed segregation pattern observed. Our findings are corroborated by the segregation pattern of the m.14487T>C mutation, which does not affect OXPHOS and MMP significantly, and its transmission is therefore predominantly determined by random genetic drift. Our conclusion is that mutation-specific selection mechanisms occur during mtDNA inheritance, which has implications for PGD and mitochondrial replacement therapy.

STUDY FUNDING/COMPETING INTEREST(S)

This work has been funded by GROW-School of Oncology and Developmental Biology. The authors declare no competing interests.

Leigh syndrome caused by mutations in MTFMT is associated with a better prognosis

Objectives

Mitochondrial methionyl‐tRNA formyltransferase (MTFMT) is required for the initiation of translation and elongation of mitochondrial protein synthesis. Pathogenic variants in MTFMT have been associated with Leigh syndrome (LS) and mitochondrial multiple respiratory chain deficiencies. We sought to elucidate the spectrum of clinical, neuroradiological and molecular genetic findings of patients with bi‐allelic pathogenic variants in MTFMT.

Methods

Retrospective cohort study combining new cases and previously published cases.

Results

Thirty‐eight patients with pathogenic variants in MTFMT were identified, including eight new cases. The median age of presentation was 14 months (range: birth to 17 years, interquartile range [IQR] 4.5 years), with developmental delay and motor symptoms being the most frequent initial manifestation. Twenty‐nine percent of the patients survived into adulthood. MRI headings in MTFMT pathogenic variants included symmetrical basal ganglia changes (62%), periventricular and subcortical white matter abnormalities (55%), and brainstem lesions (48%). Isolated complex I and combined respiratory chain deficiencies were identified in 31% and 59% of the cases, respectively. Reduction of the mitochondrial complex I and complex IV subunits was identified in the fibroblasts (13/13). Sixteen pathogenic variants were identified, of which c.626C>T was the most common. Seventy‐four percent of the patients were alive at their last clinical review (median 6.8 years, range: 14 months to 31 years, IQR 14.5 years).

Interpretation

Patients that harbour pathogenic variants in MTFMT have a milder clinical phenotype and disease progression compared to LS caused by other nuclear defects. Fibroblasts may preclude the need for muscle biopsy, to prove causality of any novel variant.

Whole Exome Sequencing Is the Preferred Strategy to Identify the Genetic Defect in Patients With a Probable or Possible Mitochondrial Cause

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.

Clinical, biochemical and genetic spectrum of 70 patients with ACAD9 deficiency: is riboflavin supplementation effective?

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.

Genetic defects in mtDNA-encoded protein translation cause pediatric, mitochondrial cardiomyopathy with early-onset brain disease

This study aims to identify gene defects in pediatric cardiomyopathy and early-onset brain disease with oxidative phosphorylation (OXPHOS) deficiencies. We applied whole-exome sequencing in three patients with pediatric cardiomyopathy and early-onset brain disease with OXPHOS deficiencies. The brain pathology was studied by MRI analysis. In consanguineous patient 1, we identified a homozygous intronic variant (c.850-3A > G) in the QRSL1 gene, which was predicted to cause abnormal splicing. The variant segregated with the disease and affected the protein function, which was confirmed by complementation studies, restoring OXPHOS function only with wild-type QRSL1. Patient 2 was compound heterozygous for two novel affected and disease-causing variants (c.[253G > A];[938G > A]) in the MTO1 gene. In patient 3, we detected one unknown affected and disease-causing variants (c.2872C > T) and one known disease-causing variant (c.1774C > T) in the AARS2 gene. The c.1774C > T variant was present in the paternal copy of the AARS2 gene, the c.2872C > T in the maternal copy. All genes were involved in translation of mtDNA-encoded proteins. Defects in mtDNA-encoded protein translation lead to severe pediatric cardiomyopathy and brain disease with OXPHOS abnormalities. This suggests that the heart and brain are particularly sensitive to defects in mitochondrial protein synthesis during late embryonic or early postnatal development, probably due to the massive mitochondrial biogenesis occurring at that stage. If both the heart and brain are involved, the prognosis is poor with a likely fatal outcome at young age.

PGD for the m.14487 T>C mitochondrial DNA mutation resulted in the birth of a healthy boy

We report on the first PGD performed for the m.14487 T>C mitochondrial DNA (mtDNA) mutation in the MT-ND6 gene, associated with Leigh syndrome. The female carrier gave birth to a healthy baby boy at age 42. This case adds to the successes of PGD for mtDNA mutations.

Anatomic & metabolic brain markers of the m.3243A>G mutation: A multi-parametric 7T MRI study

One of the most common mitochondrial DNA (mtDNA) mutations, the A to G transition at base pair 3243, has been linked to changes in the brain, in addition to commonly observed hearing problems, diabetes and myopathy. However, a detailed quantitative description of m.3243A>G patients' brains has not been provided so far. In this study, ultra-high field MRI at 7T and volume- and surface-based data analyses approaches were used to highlight morphology (i.e. atrophy)-, microstructure (i.e. myelin and iron concentration)- and metabolism (i.e. cerebral blood flow)-related differences between patients (N = 22) and healthy controls (N = 15). The use of quantitative MRI at 7T allowed us to detect subtle changes of biophysical processes in the brain with high accuracy and sensitivity, in addition to typically assessed lesions and atrophy. Furthermore, the effect of m.3243A>G mutation load in blood and urine epithelial cells on these MRI measures was assessed within the patient population and revealed that blood levels were most indicative of the brain's state and disease severity, based on MRI as well as on neuropsychological data. Morphometry MRI data showed a wide-spread reduction of cortical, subcortical and cerebellar gray matter volume, in addition to significantly enlarged ventricles. Moreover, surface-based analyses revealed brain area-specific changes in cortical thickness (e.g. of the auditory cortex), and in T₁, T₂* and cerebral blood flow as a function of mutation load, which can be linked to typically m.3243A>G-related clinical symptoms (e.g. hearing impairment). In addition, several regions linked to attentional control (e.g. middle frontal gyrus), the sensorimotor network (e.g. banks of central sulcus) and the default mode network (e.g. precuneus) were characterized by alterations in cortical thickness, T₁, T₂* and/or cerebral blood flow, which has not been described in previous MRI studies. Finally, several hypotheses, based either on vascular, metabolic or astroglial implications of the m.3243A>G mutation, are discussed that potentially explain the underlying pathobiology. To conclude, this is the first 7T and also the largest MRI study on this patient population that provides macroscopic brain correlates of the m.3243A>G mutation indicating potential MRI biomarkers of mitochondrial diseases and might guide future (longitudinal) studies to extensively track neuropathological and clinical changes.

Selection and Characterization of Palmitic Acid Responsive Patients with an OXPHOS Complex I Defect

Mitochondrial disorders are genetically and clinically heterogeneous, mainly affecting high energy-demanding organs due to impaired oxidative phosphorylation (OXPHOS). Currently, effective treatments for OXPHOS defects, with complex I deficiency being the most prevalent, are not available. Yet, clinical practice has shown that some complex I deficient patients benefit from a high-fat or ketogenic diet, but it is unclear how these therapeutic diets influence mitochondrial function and more importantly, which complex I patients could benefit from such treatment. Dietary studies in a complex I deficient patient with exercise intolerance showed increased muscle endurance on a high-fat diet compared to a high-carbohydrate diet. We performed whole-exome sequencing to characterize the genetic defect. A pathogenic homozygous p.G212V missense mutation was identified in the TMEM126B gene, encoding an early assembly factor of complex I. A complementation study in fibroblasts confirmed that the p.G212V mutation caused the complex I deficiency. The mechanism turned out to be an incomplete assembly of the peripheral arm of complex I, leading to a decrease in the amount of mature complex I. The patient clinically improved on a high-fat diet, which was supported by the 25% increase in maximal OXPHOS capacity in TMEM126B defective fibroblast by the saturated fatty acid palmitic acid, whereas oleic acid did not have any effect in those fibroblasts. Fibroblasts of other patients with a characterized complex I gene defect were tested in the same way. Patient fibroblasts with complex I defects in NDUFS7 and NDUFAF5 responded to palmitic acid, whereas ACAD9, NDUFA12, and NDUFV2 defects were non-responding. Although the data are too limited to draw a definite conclusion on the mechanism, there is a tendency that protein defects involved in early assembly complexes, improve with palmitic acid, whereas proteins defects involved in late assembly, do not. Our data show at a clinical and biochemical level that a high fat diet can be beneficial for complex I patients and that our cell line assay will be an easy tool for the selection of patients, who might potentially benefit from this therapeutic diet.

Preimplantation genetic diagnosis for mitochondrial DNA mutations: analysis of one blastomere suffices

BACKGROUND

Preimplantation genetic diagnosis (PGD) is a reproductive strategy for mitochondrial DNA (mtDNA) mutation carriers, strongly reducing their risk of affected offspring. Embryos either without the mutation or with mutation load below the phenotypic threshold are transferred to the uterus. Because of incidental heteroplasmy deviations in single blastomere and the relatively limited data available, we so far preferred relying on two blastomeres rather than one. Considering the negative effect of a two-blastomere biopsy protocol compared with a single-blastomere biopsy protocol on live birth delivery rate, we re-evaluated the error rate in our current dataset.

METHODS

For the m.3243A>G mutation, sufficient embryos/blastomeres were available for a powerful analysis. The diagnostic error rate, defined as a potential false-negative result, based on a threshold of 15%, was determined in 294 single blastomeres analysed in 73 embryos of 9 female m.3243A>G mutation carriers.

RESULTS

Only one out of 294 single blastomeres (0.34%) would have resulted in a false-negative diagnosis. False-positive diagnoses were not detected.

CONCLUSION

Our findings support a single-blastomere biopsy PGD protocol for the m.3243A>G mutation as the diagnostic error rate is very low. As in the early preimplantation embryo no mtDNA replication seems to occur and the mtDNA is divided randomly among the daughter cells, we conclude this result to be independent of the specific mutation and therefore applicable to all mtDNA mutations.

Rapid Resolution of Blended or Composite Multigenic Disease in Infants by Whole-Exome Sequencing

Whole-exome sequencing identified multiple genetic causes in 2 infants with heterogeneous disease. Three gene defects in the first patient explained all symptoms, but manifestations were overlapping (blended phenotype). Two gene defects in the second patient explained nonoverlapping symptoms (composite phenotype). Whole-exome sequencing rapidly and comprehensively resolves heterogeneous genetic disease.

Rare novel variants in the ZIC3 gene cause X-linked heterotaxy

Variants in the ZIC3 gene are rare, but have demonstrated their profound clinical significance in X-linked heterotaxy, affecting in particular male patients with abnormal arrangement of thoracic and visceral organs. Several reports have shown relevance of ZIC3 gene variants in both familial and sporadic cases and with a predominance of mutations detected in zinc-finger domains. No studies so far have assessed the functional consequences of ZIC3 variants in an in vivo model organism. A study population of 348 patients collected over more than 10 years with a large variety of congenital heart disease including heterotaxy was screened for variants in the ZIC3 gene. Functional effects of three variants were assessed both in vitro and in vivo in the zebrafish. We identified six novel pathogenic variants (1,7%), all in either male patients with heterotaxy (n=5) or a female patient with multiple male deaths due to heterotaxy in the family (n=1). All variants were located within the zinc-finger domains or leading to a truncation before these domains. Truncating variants showed abnormal trafficking of mutated ZIC3 proteins, whereas the missense variant showed normal trafficking. Overexpression of wild-type and mutated ZIC protein in zebrafish showed full non-functionality of the two frame-shift variants and partial activity of the missense variant compared with wild-type, further underscoring the pathogenic character of these variants. Concluding, we greatly expanded the number of causative variants in ZIC3 and delineated the functional effects of three variants using in vitro and in vivo model systems.

Specific MRI Abnormalities Reveal Severe Perrault Syndrome due to CLPP Defects

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.

De novo mtDNA point mutations are common and have a low recurrence risk

Background

Severe, disease-causing germline mitochondrial (mt)DNA mutations are maternally inherited or arise de novo. Strategies to prevent transmission are generally available, but depend on recurrence risks, ranging from high/unpredictable for many familial mtDNA point mutations to very low for sporadic, large-scale single mtDNA deletions. Comprehensive data are lacking for de novo mtDNA point mutations, often leading to misconceptions and incorrect counselling regarding recurrence risk and reproductive options. We aim to study the relevance and recurrence risk of apparently de novo mtDNA point mutations.

Methods

Systematic study of prenatal diagnosis (PND) and recurrence of mtDNA point mutations in families with de novo cases, including new and published data. ‘De novo’ based on the absence of the mutation in multiple (postmitotic) maternal tissues is preferred, but mutations absent in maternal blood only were also included.

Results

In our series of 105 index patients (33 children and 72 adults) with (likely) pathogenic mtDNA point mutations, the de novo frequency was 24.6%, the majority being paediatric. PND was performed in subsequent pregnancies of mothers of four de novo cases. A fifth mother opted for preimplantation genetic diagnosis because of a coexisting Mendelian genetic disorder. The mtDNA mutation was absent in all four prenatal samples and all 11 oocytes/embryos tested. A literature survey revealed 137 de novo cases, but PND was only performed for 9 (including 1 unpublished) mothers. In one, recurrence occurred in two subsequent pregnancies, presumably due to germline mosaicism.

Conclusions

De novo mtDNA point mutations are a common cause of mtDNA disease. Recurrence risk is low. This is relevant for genetic counselling, particularly for reproductive options. PND can be offered for reassurance.

Differences in Strength and Timing of the mtDNA Bottleneck between Zebrafish Germline and Non-germline Cells

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.

Polyhydramnios and cerebellar atrophy: a prenatal presentation of mitochondrial encephalomyopathy caused by mutations in the FBXL4 gene

Severe recessive mitochondrial myopathy caused by FBXL4 gene mutations may present prenatally with polyhydramnios and cerebellar hypoplasia. Characteristic dysmorphic features are: high and arched eyebrows, triangular face, a slight upslant of palpebral fissures, and a prominent pointed chin. Metabolic investigations invariably show increased serum lactate and pyruvate levels.

Leigh syndrome: Resolving the clinical and genetic heterogeneity paves the way for treatment options

Leigh syndrome is a progressive neurodegenerative disorder, affecting 1 in 40,000 live births. Most patients present with symptoms between the ages of three and twelve months, but adult onset Leigh syndrome has also been described. The disease course is characterized by a rapid deterioration of cognitive and motor functions, in most cases resulting in death due to respiratory failure. Despite the high genetic heterogeneity of Leigh syndrome, patients present with identical, symmetrical lesions in the basal ganglia or brainstem on MRI, while additional clinical manifestations and age of onset varies from case to case. To date, mutations in over 60 genes, both nuclear and mitochondrial DNA encoded, have been shown to cause Leigh syndrome, still explaining only half of all cases. In most patients, these mutations directly or indirectly affect the activity of the mitochondrial respiratory chain or pyruvate dehydrogenase complex. Exome sequencing has accelerated the discovery of new genes and pathways involved in Leigh syndrome, providing novel insights into the pathophysiological mechanisms. This is particularly important as no general curative treatment is available for this devastating disorder, although several recent studies imply that early treatment might be beneficial for some patients depending on the gene or process affected. Timely, gene-based personalized treatment may become an important strategy in rare, genetically heterogeneous disorders like Leigh syndrome, stressing the importance of early genetic diagnosis and identification of new genes/pathways. In this review, we provide a comprehensive overview of the most important clinical manifestations and genes/pathways involved in Leigh syndrome, and discuss the current state of therapeutic interventions in patients.

Mitochondrial DNA sequence characteristics modulate the size of the genetic bottleneck

With a combined carrier frequency of 1:200, heteroplasmic mitochondrial DNA (mtDNA) mutations cause human disease in ∼1:5000 of the population. Rapid shifts in the level of heteroplasmy seen within a single generation contribute to the wide range in the severity of clinical phenotypes seen in families transmitting mtDNA disease, consistent with a genetic bottleneck during transmission. Although preliminary evidence from human pedigrees points towards a random drift process underlying the shifting heteroplasmy, some reports describe differences in segregation pattern between different mtDNA mutations. However, based on limited observations and with no direct comparisons, it is not clear whether these observations simply reflect pedigree ascertainment and publication bias. To address this issue, we studied 577 mother–child pairs transmitting the m.11778G>A, m.3460G>A, m.8344A>G, m.8993T>G/C and m.3243A>G mtDNA mutations. Our analysis controlled for inter-assay differences, inter-laboratory variation and ascertainment bias. We found no evidence of selection during transmission but show that different mtDNA mutations segregate at different rates in human pedigrees. m.8993T>G/C segregated significantly faster than m.11778G>A, m.8344A>G and m.3243A>G, consistent with a tighter mtDNA genetic bottleneck in m.8993T>G/C pedigrees. Our observations support the existence of different genetic bottlenecks primarily determined by the underlying mtDNA mutation, explaining the different inheritance patterns observed in human pedigrees transmitting pathogenic mtDNA mutations.

Hypertrophic remodelling in cardiac regulatory myosin light chain (MYL2) founder mutation carriers

AIMS

Phenotypic heterogeneity and incomplete penetrance are common in patients with hypertrophic cardiomyopathy (HCM). We aim to improve the understanding in genotype-phenotype correlations in HCM, particularly the contribution of an MYL2 founder mutation and risk factors to left ventricular hypertrophic remodelling.

METHODS AND RESULTS

We analysed 14 HCM families of whom 38 family members share the MYL2 c.64G > A [p.(Glu22Lys)] mutation and a common founder haplotype. In this unique cohort, we investigated factors influencing phenotypic outcome in addition to the primary mutation. The mutation alone showed benign disease manifestation with low penetrance. The co-presence of additional risk factors for hypertrophy such as hypertension, obesity, or other sarcomeric gene mutation increased disease penetrance substantially and caused HCM in 89% of MYL2 mutation carriers (P = 0.0005). The most prominent risk factor was hypertension, observed in 71% of mutation carriers with HCM and an additional risk factor.

CONCLUSION

The MYL2 mutation c.64G > A on its own is incapable of triggering clinical HCM in most carriers. However, the presence of an additional risk factor for hypertrophy, particularly hypertension, adds to the development of HCM. Early diagnosis of risk factors is important for early treatment of MYL2 mutation carriers and close monitoring should be guaranteed in this case. Our findings also suggest that the presence of hypertension or another risk factor for hypertrophy should not be an exclusion criterion for genetic studies.

Mitochondrial DNA copy number in colorectal cancer: between tissue comparisons, clinicopathological characteristics and survival

Low mitochondrial DNA (mtDNA) copy number in tumors has been associated with worse prognosis in colorectal cancer (CRC). This study further deciphers the role of mtDNA copy number in CRC by comparing mtDNA copy number between healthy, adenoma and carcinoma tissue, by investigating its association according to several clinicopathological characteristics in CRC, and by relating it to CRC-specific survival in CRC patients. A hospital-based series of samples including cancer, adenoma and adjacent histologically normal tissue from primary CRC patients (n = 56) and recurrent CRC (n = 16) was studied as well as colon mucosa samples from healthy subjects (n = 76). Furthermore, mtDNA copy number was assessed in carcinomas of 693 CRC cases identified from the population-based Netherlands Cohort Study (NLCS). MtDNA copy number was significantly lower in carcinoma tissue (P = 0.011) and adjacent tissue (P < 0.001) compared to earlier resected adenoma tissue and in primary CRC tissue compared to recurrent CRC tissue (P = 0.011). Within both study populations, mtDNA copy number was significantly lower in mutated BRAF (P = 0.027 and P = 0.006) and in microsatellite unstable (MSI) tumors (P = 0.033 and P < 0.001) and higher in KRAS mutated tumors (P = 0.004). Furthermore, the association between mtDNA and survival seemed to follow an inverse U-shape with the highest HR observed in the second quintile of mtDNA copy number (HR = 1.70, 95% CI = 1.18, 2.44) compared to the first quintile. These results might reflect an association of mtDNA copy number with various malignant processes in cancer cells and warrants further research on tumor energy metabolism in CRC prognosis.

Transcriptome analysis of complex I-deficient patients reveals distinct expression programs for subunits and assembly factors of the oxidative phosphorylation system

Background

Transcriptional control of mitochondrial metabolism is essential for cellular function. A better understanding of this process will aid the elucidation of mitochondrial disorders, in particular of the many genetically unsolved cases of oxidative phosphorylation (OXPHOS) deficiency. Yet, to date only few studies have investigated nuclear gene regulation in the context of OXPHOS deficiency. In this study we performed RNA sequencing of two control and two complex I-deficient patient cell lines cultured in the presence of compounds that perturb mitochondrial metabolism: chloramphenicol, AICAR, or resveratrol. We combined this with a comprehensive analysis of mitochondrial and nuclear gene expression patterns, co-expression calculations and transcription factor binding sites.

Results

Our analyses show that subsets of mitochondrial OXPHOS genes respond opposingly to chloramphenicol and AICAR, whereas the response of nuclear OXPHOS genes is less consistent between cell lines and treatments. Across all samples nuclear OXPHOS genes have a significantly higher co-expression with each other than with other genes, including those encoding mitochondrial proteins. We found no evidence for complex-specific mRNA expression regulation: subunits of different OXPHOS complexes are similarly (co-)expressed and regulated by a common set of transcription factors. However, we did observe significant differences between the expression of nuclear genes for OXPHOS subunits versus assembly factors, suggesting divergent transcription programs. Furthermore, complex I co-expression calculations identified 684 genes with a likely role in OXPHOS biogenesis and function. Analysis of evolutionarily conserved transcription factor binding sites in the promoters of these genes revealed almost all known OXPHOS regulators (including GABP, NRF1/2, SP1, YY1, E-box factors) and a set of novel candidates (ELK1, KLF7, SP4, EHF, ZNF143, and TEL2).

Conclusions

OXPHOS genes share an expression program distinct from other genes encoding mitochondrial proteins, indicative of targeted nuclear regulation of a mitochondrial sub-process. Within the subset of OXPHOS genes we established a difference in expression between mitochondrial and nuclear genes, and between nuclear genes encoding subunits and assembly factors. Most transcription regulators of genes that co-express with complex I are well-established factors for OXPHOS biogenesis. For the remaining six factors we here suggest for the first time a link with transcription regulation in OXPHOS deficiency.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1883-8) contains supplementary material, which is available to authorized users.

Pathogenic CWF19L1 variants as a novel cause of autosomal recessive cerebellar ataxia and atrophy

Autosomal recessive cerebellar ataxia (ARCA) is a group of neurological disorders characterized by degeneration or abnormal development of the cerebellum and spinal cord. ARCA is clinically and genetically highly heterogeneous, with over 20 genes involved. Exome sequencing of a girl with ARCA from non-consanguineous Dutch parents revealed two pathogenic variants c.37G>C; p.D13H and c.946A>T; p.K316* in CWF19L1, a gene with an unknown function, recently reported to cause ARCA in a Turkish family. Sanger sequencing showed that the c.37G>C variant was inherited from the father and the c.946A>T variant from the mother. Pathogenicity was based on the damaging effect on protein function as the c.37G>C variant changed the highly conserved, negatively charged aspartic acid to the positively charged histidine and the c.946A>T variant introduced a premature stop codon. In addition, 27 patients with ARCA were tested for pathogenic variants in CWF19L1, however, no pathogenic variants were identified. Our data confirm CWF19L1 as a novel but rare gene causing ARCA.

Evolutionary defined role of the mitochondrial DNA in fertility, disease and ageing

BACKGROUND

The endosymbiosis of an alpha-proteobacterium and a eubacterium a billion years ago paved the way for multicellularity and enabled eukaryotes to flourish. The selective advantage for the host was the acquired ability to generate large amounts of intracellular hydrogen-dependent adenosine triphosphate. The price was increased reactive oxygen species (ROS) inside the eukaryotic cell, causing high mutation rates of the mitochondrial DNA (mtDNA). According to the Muller's ratchet theory, this accumulation of mutations in asexually transmitted mtDNA would ultimately lead to reduced reproductive fitness and eventually extinction. However, mitochondria have persisted over the course of evolution, initially due to a rapid, extreme evolutionary reduction of the mtDNA content. After the phylogenetic divergence of eukaryotes into animals, fungi and plants, differences in evolution of the mtDNA occurred with different adaptations for coping with the mutation burden within these clades. As a result, mitochondrial evolutionary mechanisms have had a profound effect on human adaptation, fertility, healthy reproduction, mtDNA disease manifestation and transmission and ageing. An understanding of these mechanisms might elucidate novel approaches for treatment and prevention of mtDNA disease.

METHODS

The scientific literature was investigated to determine how mtDNA evolved in animals, plants and fungi. Furthermore, the different mechanisms of mtDNA inheritance and of balancing Muller's ratchet in these species were summarized together with the consequences of these mechanisms for human health and reproduction.

RESULTS

Animal, plant and fungal mtDNA have evolved differently. Animals have compact genomes, little recombination, a stable number of genes and a high mtDNA copy number, whereas plants have larger genomes with variable gene counts, a low mtDNA copy number and many recombination events. Fungal mtDNA is somewhere in between. In plants, the mtDNA mutation rate is kept low by effective ROS defence and efficient recombination-mediated mtDNA repair. In animal mtDNA, these mechanisms are not or less well-developed and the detrimental mutagenesis events are controlled by a high mtDNA copy number in combination with a genetic bottleneck and purifying selection during transmission. The mtDNA mutation rates in animals are higher than in plants, which allow mobile animals to adapt more rapidly to various environmental conditions in terms of energy production, whereas static plants do not have this need. Although at the level of the species, these mechanisms have been extremely successful, they can have adverse effects for the individual, resulting, in humans, in severe or unpredictably segregating mtDNA diseases, as well as fertility problems and unhealthy ageing.

CONCLUSIONS

Understanding the forces and processes that underlie mtDNA evolution among different species increases our knowledge on the detrimental consequences that individuals can have from these evolutionary end-points. Alternative outcomes in animals, fungi and plants will lead to a better understanding of the inheritance of mtDNA disorders and mtDNA-related fertility problems. These will allow the development of options to ameliorate, cure and/or prevent mtDNA diseases and mtDNA-related fertility problems.

Preventing the transmission of mitochondrial DNA disorders: selecting the good guys or kicking out the bad guys

Mitochondrial disorders represent the most common group of inborn errors of metabolism. Clinical manifestations can be extremely variable, ranging from single affected tissues to multisystemic syndromes. Maternally inherited mitochondrial DNA (mtDNA) mutations are a frequent cause, affecting about one in 5000 individuals. The expression of mtDNA mutations differs from nuclear gene defects. Mutations are either homoplasmic or heteroplasmic, and in the latter case disease becomes manifest when the mutation load exceeds a tissue-specific threshold. Mutation load can vary between tissues and in time, and often an exact correlation between mutation load and clinical manifestations is lacking. Because of the possible clinical severity, the lack of treatment and the high recurrence risk of affected offspring for female carriers, couples request prevention of transmission of mtDNA mutations. Previously, choices have been limited due to a segregational bottleneck, which makes the mtDNA mutation load in embryos highly variable and the consequences largely unpredictable. However, recently it was shown that preimplantation genetic diagnosis offers a fair chance of unaffected offspring to carriers of heteroplasmic mtDNA mutations. Technically and ethically challenging possibilities, such maternal spindle transfer and pronuclear transfer, are emerging and providing carriers additional prospects of giving birth to a healthy child.

Myelin-derived lipids modulate macrophage activity by liver X receptor activation

Multiple sclerosis is a chronic, inflammatory, demyelinating disease of the central nervous system in which macrophages and microglia play a central role. Foamy macrophages and microglia, containing degenerated myelin, are abundantly found in active multiple sclerosis lesions. Recent studies have described an altered macrophage phenotype after myelin internalization. However, it is unclear by which mechanisms myelin affects the phenotype of macrophages and how this phenotype can influence lesion progression. Here we demonstrate, by using genome wide gene expression analysis, that myelin-phagocytosing macrophages have an enhanced expression of genes involved in migration, phagocytosis and inflammation. Interestingly, myelin internalization also induced the expression of genes involved in liver-X-receptor signaling and cholesterol efflux. In vitro validation shows that myelin-phagocytosing macrophages indeed have an increased capacity to dispose intracellular cholesterol. In addition, myelin suppresses the secretion of the pro-inflammatory mediator IL-6 by macrophages, which was mediated by activation of liver-X-receptor β. Our data show that myelin modulates the phenotype of macrophages by nuclear receptor activation, which may subsequently affect lesion progression in demyelinating diseases such as multiple sclerosis.

Physical activity is the key determinant of skeletal muscle mitochondrial function in type 2 diabetes

CONTEXT

Conflicting data exist on mitochondrial function and physical activity in type 2 diabetes mellitus (T2DM) development.

OBJECTIVE

The aim was to assess mitochondrial function at different stages during T2DM development in combination with physical exercise in longstanding T2DM patients.

DESIGN AND METHODS

We performed cross-sectional analysis of skeletal muscle from 12 prediabetic 11 longstanding T2DM male subjects and 12 male controls matched by age and body mass index.

INTERVENTION

One-year intrasubject controlled supervised exercise training intervention was done in longstanding T2DM patients.

MAIN OUTCOME MEASUREMENTS

Extensive ex vivo analyses of mitochondrial quality, quantity, and function were collected and combined with global gene expression analysis and in vivo ATP production capacity after 1 yr of training.

RESULTS

Mitochondrial density, complex I activity, and the expression of Krebs cycle and oxidative phosphorylation system-related genes were lower in longstanding T2DM subjects but not in prediabetic subjects compared with controls. This indicated a reduced capacity to generate ATP in longstanding T2DM patients only. Gene expression analysis in prediabetic subjects suggested a switch from carbohydrate toward lipid as an energy source. One year of exercise training raised in vivo skeletal muscle ATP production capacity by 21 ± 2% with an increased trend in mitochondrial density and complex I activity. In addition, expression levels of β-oxidation, Krebs cycle, and oxidative phosphorylation system-related genes were higher after exercise training.

CONCLUSIONS

Mitochondrial dysfunction is apparent only in inactive longstanding T2DM patients, which suggests that mitochondrial function and insulin resistance do not depend on each other. Prolonged exercise training can, at least partly, reverse the mitochondrial impairments associated with the longstanding diabetic state.

Patient-derived fibroblasts indicate oxidative stress status and may justify antioxidant therapy in OXPHOS disorders

Oxidative phosphorylation disorders are often associated with increased oxidative stress and antioxidant therapy is frequently given as treatment. However, the role of oxidative stress in oxidative phosphorylation disorders or patients is far from clear and consequently the preventive or therapeutic effect of antioxidants is highly anecdotic. Therefore, we performed a systematic study of a panel of oxidative stress parameters (reactive oxygen species levels, damage and defense) in fibroblasts of twelve well-characterized oxidative phosphorylation patients with a defect in the POLG1 gene, in the mitochondrial DNA-encoded tRNA-Leu gene (m.3243A>G or m.3302A>G) and in one of the mitochondrial DNA-encoded NADH dehydrogenase complex I (CI) subunits. All except two cell lines (one POLG1 and one tRNA-Leu) showed increased reactive oxygen species levels compared with controls, but only four (two CI and two tRNA-Leu) cell lines provided evidence for increased oxidative protein damage. The absence of a correlation between reactive oxygen species levels and oxidative protein damage implies differences in damage prevention or correction. This was investigated by gene expression studies, which showed adaptive and compensating changes involving antioxidants and the unfolded protein response, especially in the POLG1 group. This study indicated that patients display individual responses and that detailed analysis of fibroblasts enables the identification of patients that potentially benefit from antioxidant therapy. Furthermore, the fibroblast model can also be used to search for and test novel, more specific antioxidants or explore ways to stimulate compensatory mechanisms.

Large scale mtDNA sequencing reveals sequence and functional conservation as major determinants of homoplasmic mtDNA variant distribution

The mitochondrial DNA (mtDNA) is highly variable, containing large numbers of pathogenic mutations and neutral polymorphisms. The spectrum of homoplasmic mtDNA variation was characterized in 730 subjects and compared with known pathogenic sites. The frequency and distribution of variants in protein coding genes were inversely correlated with conservation at the amino acid level. Analysis of tRNA secondary structures indicated a preference of variants for the loops and some acceptor stem positions. This comprehensive overview of mtDNA variants distinguishes between regions and positions which are likely not critical, mainly conserved regions with pathogenic mutations and essential regions containing no mutations at all.

Biochemical abnormalities in a patient with thymidine phosphorylase deficiency with fatal outcome

Deficiency of the cytosolic enzyme thymidine phosphorylase (TP) causes a multisystem disorder called mitochondrial neurogastrointestinal encephalomyopathy (MNGIE) syndrome. Clinical symptoms are gastrointestinal dysfunction, muscle involvement and neurological deterioration. TP deficiency is biochemically characterised by accumulation of thymidine and deoxyuridine in body fluids and compromised mitochondrial deoxyribose nucleic acid (mtDNA) integrity (depletion and multiple deletions). In this report we describe a patient with the clinical and biochemical features related to the end stage of the disease. Home parenteral nutrition had started to improve the clinical condition and preparations were initiated for stem cell transplantation (SCT) as a last resort treatment. Unfortunately, the patient died during the induction phase of SCT. This report shows that TP deficiency is a severe clinical condition with a broad spectrum of affected tissues. TP deficiency can be easily determined by the measurement of pyrimidine metabolites in body fluids and TP activity in peripheral blood leucocytes. Early detection and treatment may prevent the progress of the clinical symptoms and, therefore, should be considered for inclusion in newborn screening programmes.