Functional recurrent mutations in the human mitochondrial phylogeny: dual roles in evolution and disease

Natural selection shaped the protective effect of the mtDNA lineage against obesity in Han Chinese populations

Mitochondria play a key role in lipid metabolism, and mitochondrial DNA (mtDNA) mutations are thus considered to affect obesity susceptibility by altering oxidative phosphorylation and mitochondrial function. In this study, we investigated mtDNA variants that may affect obesity risk in 2877 Han Chinese individuals from three independent populations. The association analysis of 16 basal mtDNA haplogroups with body mass index, waist circumference, and waist-to-hip ratio revealed that only haplogroup M7 was significantly negatively correlated with all three adiposity-related anthropometric traits in the overall cohort, verified by the analysis of a single population, i.e., the Zhengzhou population. Furthermore, subhaplogroup analysis suggested that M7b1a1 was the most likely haplogroup associated with a decreased obesity risk, and the variation T12811C (causing Y159H in ND5) harbored in M7b1a1 may be the most likely candidate for altering the mitochondrial function. Specifically, we found that proportionally more nonsynonymous mutations accumulated in M7b1a1 carriers, indicating that M7b1a1 was either under positive selection or subject to a relaxation of selective constraints. We also found that nuclear variants, especially in DACT2 and PIEZO1, may functionally interact with M7b1a1.

The Association Between Mitochondrial tRNAGlu Variants and Hearing Loss: A Case-Control Study

Purpose

This study aimed to examine the frequencies of mt-tRNAGlu variants in 180 pediatric patients with non-syndromic hearing loss (NSHL) and 100 controls.

Methods

Sanger sequencing was performed to screen for mt-tRNAGlu variants. These mitochondrial DNA (mtDNA) pathogenic mutations were further assessed using phylogenetic conservation and haplogroup analyses. We also traced the origins of the family history of probands carrying potential pathogenic mtDNA mutations. Mitochondrial functions including mtDNA content, ATP and reactive oxygen species (ROS) were examined in cells derived from patients carrying the mt-tRNAGlu A14692G and CO1/tRNASer(UCN) G7444A variants and controls.

Results

We identified four possible pathogenic variants: m.T14709C, m.A14683G, m.A14692G and m.A14693G, which were found in NSHL patients but not in controls. Genetic counseling suggested that one child with the m.A14692G variant had a family history of NSHL. Sequence analysis of mtDNA suggested the presence of the CO1/tRNASer(UCN) G7444A and mt-tRNAGlu A14692G variants. Molecular analysis suggested that, compared with the controls, patients with these variants exhibited much lower mtDNA copy numbers, ATP production, whereas ROS levels increased (p<0.05 for all), suggesting that the m.A14692G and m.G7444A variants led to mitochondrial dysfunction.

Conclusion

mt-tRNAGlu variants are important risk factors for NSHL.

Maternally Inherited Essential Hypertension May Be Associated with the Mutations in Mitochondrial tRNAGlu Gene

Background

Mitochondrial DNA (mtDNA) mutations are associated with essential hypertension (EH), but the molecular mechanism remains largely unknown.

Objective

The aim of this study is to explore the association between mtDNA mutations and EH.

Methods

Two maternally inherited families with EH are underwent clinical, genetic and biochemical assessments. mtDNA mutations are screened by PCR-Sanger sequencing and phylogenetic, and bioinformatics analyses are performed to evaluate the pathogenicity of mtDNA mutations. We also generate cytoplasmic hybrid (cybrid) cell lines to analysis mitochondrial functions.

Results

Matrilineal relatives exhibit variable degree of clinical phenotypes. Molecular analysis reveals the presence of m.A14693G and m.A14696G mutations in two pedigrees. Notably, the m.A14693G mutation occurs at position 54 in the TψC loop of tRNAGlu, a position which is critical for post-transcriptionally modification of tRNAGlu. While the m.A14696G mutation creates a novel base-pairing (51C-64G). Bioinformatic analysis shows that these mutations alter tRNAGlu secondary structure. Additionally, patients with tRNAGlu mutations exhibit markedly decreased in mtDNA copy number, mitochondrial membrane potential (MMP) and ATP, whereas the levels of reactive oxygen species (ROS) increase significantly.

Conclusion

The m.A14696G and m.A14693G mutations lead to failure in tRNAGlu metabolism and cause mitochondrial dysfunction that is responsible for EH.

The metazoan landscape of mitochondrial DNA gene order and content is shaped by selection and affects mitochondrial transcription

Mitochondrial DNA (mtDNA) harbors essential genes in most metazoans, yet the regulatory impact of the multiple evolutionary mtDNA rearrangements has been overlooked. Here, by analyzing mtDNAs from ~8000 metazoans we found high gene content conservation (especially of protein and rRNA genes), and codon preferences for mtDNA-encoded tRNAs across most metazoans. In contrast, mtDNA gene order (MGO) was selectively constrained within but not between phyla, yet certain gene stretches (ATP8-ATP6, ND4-ND4L) were highly conserved across metazoans. Since certain metazoans with different MGOs diverge in mtDNA transcription, we hypothesized that evolutionary mtDNA rearrangements affected mtDNA transcriptional patterns. As a first step to test this hypothesis, we analyzed available RNA-seq data from 53 metazoans. Since polycistron mtDNA transcripts constitute a small fraction of the steady-state RNA, we enriched for polycistronic boundaries by calculating RNA-seq read densities across junctions between gene couples encoded either by the same strand (SSJ) or by different strands (DSJ). We found that organisms whose mtDNA is organized in alternating reverse-strand/forward-strand gene blocks (mostly arthropods), displayed significantly reduced DSJ read counts, in contrast to organisms whose mtDNA genes are preferentially encoded by one strand (all chordates). Our findings suggest that mtDNA rearrangements are selectively constrained and likely impact mtDNA regulation.

Mitochondrial DNA haplogroup analysis in Saudi Arab patients with multiple sclerosis

Previous studies have suggested that mitochondrial DNA (mtDNA) variants are associated with multiple sclerosis (MS), a complex neurodegenerative immune-mediated disease of the central nervous system. Since mtDNA is maternally inherited without recombination, specific mtDNA variants defining genetic background are associated with the susceptibility to human diseases. To assess the contribution of mtDNA haplogroups to the predisposition of MS in an Arab population, we analysed sequencing data of mitochondrial genomes from 47 native Saudi Arab individuals including 23 patients with relapsing-remitting MS (RRMS) and 24 healthy controls. All patients and controls could be classified into ten haplogroups. The European-specific haplogroup U was more prevalent in patients than in the controls (26.1% vs. 4.2%), whereas haplogroup T was only present in patients and haplogroups HV and N were only found in controls. Haplogroup U was significantly association with increased risk of MS (odds ratio = 6.26, p<0.05), although the association did not maintain significance after adjustment for multiple comparisons. Haplotype U was more prevalent in patients with younger age of onset (p = 0.006), but there was no relationship between haplotype U and disease severity, disease duration or EDSS and age-matched carriers and non-carriers of haplogroup U (p>0.05). Definition site of haplogroup U include the variant m.12308A>G in MT-TL2 gene which was found to affect highly conserved position within the variable arm of tRNALeu(CUN) and thus may impact mitochondrial protein synthesis, and two other variants namely m.11467A>G in MT-ND4 gene and m.12372G>A in MT-ND5 gene which were previously linked with mitochondrial function. Despite the small number of subjects, which may limit the statistical power of the study, our results showed for the first time a possible contribution of haplogroup U to the predisposition to MS in an Arab population. These findings warrant further validation in a large cohort to distinguish a genuine effect specific to MS from a chance finding due to small sampling.

Human mitochondrial RNA modifications associate with tissue-specific changes in gene expression, and are affected by sunlight and UV exposure

RNA-DNA differences (RDD) have previously been identified in the human mitochondrial RNA (mt-RNA) transcripts, yet their functional impact is poorly understood. By analyzing 4928 RNA-seq samples from 23 body sites, we found that mtDNA gene expression negatively correlated with the levels of both m1A 947 16 S rRNA modification (mtDNA position 2617) and the m1A 1812 ND5 mRNA modification (mtDNA position 13,710) in 15 and 14 body sites, respectively. Such correlation was not evident in all tested brain tissues, thus suggesting a tissue-specific impact of these modifications on mtDNA gene expression. To assess the response of the tested modifications to environmental cues, we analyzed pairs of skin samples that were either exposed to the sun or not. We found that the correlations of mtDNA gene expression with both mt-RNA modifications were compromised upon sun exposure. As a first step to explore the underlying mechanism, we analyzed RNA-seq data from keratinocytes that were exposed to increasing doses of UV irradiation. Similar to sun exposure, we found a significant decrease in mtDNA gene expression upon increase in UV dosage. In contrast, there was a significant increase in the m1A 947 16 S rRNA modification levels upon elevation in UV dose. Finally, we identified candidate modulators of such responses. Taken together, our results indicate that mt-RNA modifications functionally correlate with mtDNA gene expression, and responds to environmental cues, hence supporting their physiological importance.

Mitochondrial tRNAGln 4394C>T Mutation May Contribute to the Clinical Expression of 1555A>G-Induced Deafness

The mitochondrial 1555A>G mutation plays a critical role in aminoglycoside-induced and non-syndromic hearing loss (AINSHL). Previous studies have suggested that mitochondrial secondary variants may modulate the clinical expression of m.1555A>G-induced deafness, but the molecular mechanism has remained largely undetermined. In this study, we investigated the contribution of a deafness-associated tRNAGln 4394C>T mutation to the clinical expression of the m.1555A>G mutation. Interestingly, a three-generation family with both the m.1555A>G and m.4394C>T mutations exhibited a higher penetrance of hearing loss than another family harboring only the m.1555A>G mutation. At the molecular level, the m.4394C>T mutation resides within a very conserved nucleotide of tRNAGln, which forms a new base-pairing (7T-66A) and may affect tRNA structure and function. Using trans-mitochondrial cybrid cells derived from three subjects with both the m.1555A>G and m.4394C>T mutations, three patients with only the m.1555A>G mutation and three control subjects without these primary mutations, we observed that cells with both the m.1555A>G and m.4394C>T mutations exhibited more severely impaired mitochondrial functions than those with only the m.1555A>G mutation. Furthermore, a marked decrease in mitochondrial RNA transcripts and respiratory chain enzymes was observed in cells harboring both the m.1555A>G and m.4394C>T mutations. Thus, our data suggest that the m.4394C>T mutation may play a synergistic role in the m.1555A>G mutation, enhancing mitochondrial dysfunctions and contributing to a high penetrance of hearing loss in families with both mtDNA pathogenic mutations.

Helena's Many Daughters: More Mitogenome Diversity behind the Most Common West Eurasian mtDNA Control Region Haplotype in an Extended Italian Population Sample

The high number of matching haplotypes of the most common mitochondrial (mt)DNA lineages are considered to be the greatest limitation for forensic applications. This study investigates the potential to solve this constraint by massively parallel sequencing a large number of mitogenomes that share the most common West Eurasian mtDNA control region (CR) haplotype motif (263G 315.1C 16519C). We augmented a pilot study on 29 to a total of 216 Italian mitogenomes that represents the largest set of the most common CR haplotype compiled from a single country. The extended population sample confirmed and extended the huge coding region diversity behind the most common CR motif. Complete mitogenome sequencing allowed for the detection of 163 distinct haplotypes, raising the power of discrimination from 0 (CR) to 99.6% (mitogenome). The mtDNAs were clustered into 61 named clades of haplogroup H and did not reveal phylogeographic trends within Italy. Rapid individualization approaches for investigative purposes are limited to the most frequent H clades of the dataset, viz. H1, H3, and H7.

Next-generation sequencing of the whole mitochondrial genome identifies functionally deleterious mutations in patients with multiple sclerosis

Multiple sclerosis (MS) is an immune-mediated disease of the central nervous system with genetics and environmental determinants. Studies focused on the neurogenetics of MS showed that mitochondrial DNA (mtDNA) mutations that can ultimately lead to mitochondrial dysfunction, alter brain energy metabolism and cause neurodegeneration. We analyzed the whole mitochondrial genome using next-generation sequencing (NGS) from 47 Saudi individuals, 23 patients with relapsing-remitting MS and 24 healthy controls to identify mtDNA disease-related mutations/variants. A large number of variants were detected in the D-loop and coding genes of mtDNA. While distinct unique variants were only present in patients or only occur in controls, a number of common variants were shared among the two groups. The prevalence of some common variants differed significantly between patients and controls, thus could be implicated in susceptibility to MS. Of the unique variants only present in the patients, 34 were missense mutations, located in different mtDNA-encoded genes. Seven of these mutations were not previously reported in MS, and predicted to be deleterious with considerable impacts on the functions and structures of encoded-proteins and may play a role in the pathogenesis of MS. These include two heteroplasmic mutations namely 10237T>C in MT-ND3 gene and 15884G>C in MT-CYB gene; and three homoplasmic mutations namely 9288A>G in MT-CO3 gene, 14484T>C in MT-ND6 gene, 15431G>A in MT-CYB gene, 8490T>C in MT-ATP8 gene and 5437C>T in MT-ND2 gene. Notably some patients harboured multiple mutations while other patients carried the same mutations. This study is the first to sequence the entire mitochondrial genome in MS patients in an Arab population. Our results expanded the mutational spectrum of mtDNA variants in MS and highlighted the efficiency of NGS in population-specific mtDNA variant discovery. Further investigations in a larger cohort are warranted to confirm the role of mtDNA MS.

Maternally transmitted nonsyndromic hearing impairment may be associated with mitochondrial tRNAAla 5601C>T and tRNALeu(CUN) 12311T>C mutations

BACKGROUND

Sequence alternations in mitochondrial genomes, especially in genes encoding mitochondrial tRNA (mt-tRNA), were the important contributors to nonsyndromic hearing loss (NSHL); however, the molecular mechanisms remained largely undetermined.

METHODS

A maternally transmitted Chinese pedigree with NSHL underwent clinical, genetic, and biochemical assessment. PCR and direct sequence analyses were performed to detect mitochondrial DNA (mtDNA), GJB2, and SLC26A4 gene mutations from matrilineal relatives of this family. Mitochondrial functions including mitochondrial membrane potential (MMP), ATP, and ROS were evaluated in polymononuclear leukocytes (PMNs) derived from three deaf patients and three controls from this pedigree.

RESULTS

Four of nine matrilineal relatives developed hearing loss at the variable age of onset. Two putative pathogenic mutations, m.5601C>T in tRNA^Ala and m.12311T>C in tRNA^Leu(CUN) , were identified via PCR-Sanger sequencing, as well as 34 variants that belonged to mtDNA haplogroup G2b2. Intriguingly, m.5601C>T mutation resided at very conserved nucleotide in the TψC loop of tRNA^Ala (position 59), while the T-to-C substitution at position 12311 located at position 48 in the variable stem of tRNA^Leu(CUN) and was believed to alter the aminoacylation and the steady-state level of tRNA. Biochemical analysis revealed the impairment of mitochondrial functions including the significant reductions of ATP and MMP, whereas markedly increased ROS levels were found in PMNs derived from NSHL patients with m.5601C>T and m.12311T>C mutations. However, we did not detect any mutations in GJB2 and SLC26A4 genes.

CONCLUSION

Our data indicated that mt-tRNA^Ala m.5601C>T and tRNA^Leu(CUN) 12311T>C mutations were associated with NSHL.

Maternally transmitted diabetes mellitus may be associated with mitochondrial ND5 T12338C and tRNAAla T5587C variants

INTRODUCTION

Mutations/variants in mitochondrial genomes are found to be associated with type 2 diabetes mellitus (T2DM), but the pathophysiology of this disease remains largely unknown.

AIM

The aim of this study is to investigate the relationship between mitochondrial DNA (mtDNA) variants and T2DM.

METHODOLOGY

A maternally inherited T2DM pedigree is underwent clinical, genetic, and molecular assessment. Moreover, the complete mitochondrial genomes of the matrilineal relatives of this family are PCR amplified and sequenced. We also utilize the phylogenetic conservation analysis, haplogroup classification, and the pathogenicity scoring system to determine the T2DM-associated potential pathogenic mtDNA variants.

RESULT

Four of seven matrilineal relatives of this pedigree suffered from T2DM with variable ages of onset. Screening for the entire mtDNA genes of matrilineal members reveals co-existence of ND5 T12338C and tRNA^Ala T5587C variants, as well as 21 genetic polymorphisms which belong to East Asian haplogroup F2. Interestingly, the T12338C variant causes the alternation of first amino acid Met to Thr, shortened two amino acids of ND5 protein. Furthermore, T5587C variant is located at position 73 in the 3'end of mt-tRNA^Ala and may have structural and functional consequences.

CONCLUSIONS

The co-occurrence of ND5 T12338C and tRNA^Ala T5587C variants may impair the mitochondrial function, which are associated with the development of T2DM in this family.

Contribution of mitochondrial gene variants in diabetes and diabetic kidney disease

OBJECTIVES

Mitochondrial DNA (mtDNA) plays an important role in the pathogenesis of diabetes. Variants in mtDNA have been reported in diabetes, but studies on the whole mtDNA variants were limited. Our study aims to explore the association of whole mtDNA variants with diabetes and diabetic kidney disease (DKD).

METHODS

The whole mitochondrial genome was screened by next-generation sequencing in cohort 1 consisting of 50 early-onset diabetes (EOD) patients with a maternally inherited diabetes (MID) family history. A total of 42 variants possibly associated with mitochondrial diseases were identified according to the filtering strategy. These variants were sequenced in cohort 2 consisting of 90 EOD patients with MID. The association between the clinical phenotype and these variants was analyzed. Then, these variants were genotyped in cohort 3 consisting of 1,571 type 2 diabetes mellitus patients and 496 subjects with normal glucose tolerance (NGT) to analyze the association between variants with diabetes and DKD.

RESULTS

Patients with variants in the non-coding region had a higher percentage of obesity and levels of fasting insulin (62.1% vs. 24.6%, P = 0.001; 80.0% vs. 26.5% P < 0.001). The patients with the variants in rRNA had a higher prevalence of obesity (71.4% vs. 30.3%, P = 0.007), and the patients with the variants in mitochondrial complex I had a higher percentage of the upper tertile of FINS (64.3% vs. 34.3%, P = 0.049). Among 20 homogeneous variants successfully captured, two known variants (m.A3943G, m.A10005G) associated with other mitochondrial diseases were only in the diabetic group, but not in the NGT group, which perhaps indicated its possible association with diabetes. The prevalence of DKD was significantly higher in the group with the 20 variants than those without these variants (18.7% vs. 14.6%, P = 0.049) in the participants with diabetes of cohort 3.

CONCLUSION

MtDNA variants are associated with MID and DKD, and our findings advance our understanding of mtDNA in diabetes and DKD. It will have important implications for the individual therapy of mitochondrial diabetes.

Mitochondrial Diabetes is Associated with tRNALeu(UUR) A3243G and ND6 T14502C Mutations

BACKGROUND

Mutations in mitochondrial DNA (mtDNA) are associated with type 2 diabetes mellitus (T2DM). In particular, m.A3243G is the most common T2DM-related mtDNA mutation in many families worldwide. However, the clinical features and pathophysiology of m.A3243G-induced T2DM are largely undefined.

METHODS

Two pedigrees with maternally inherited T2DM were underwent clinical, molecular and biochemical assessments. The mtDNA genes were PCR amplified and sequenced. Mitochondrial adenosine triphosphate (ATP) and reactive oxygen species (ROS) were measured in polymononuclear leukocytes derived from three patients with both the m.A3243G and m.T14502C mutations, three patients with only the m.A3243G mutation and three controls without these mutations. Moreover, GJB2, GJB3 and GJB6 mutations were screened by PCR-Sanger sequencing.

RESULTS

Members of the two pedigrees manifestated variable clinical phenotypes including diabetes and hearing and vision impairments. The age at onset of T2DM varied from 31 to 66 years, with an average of 41 years. Mutational analysis of mitochondrial genomes indicated the presence of the m.A3243G mutation in both pedigrees. Matrilineal relatives in one of the pedigrees harbored the coexisting of m.A3243G and m.T14502C mutations. Remarkably, the m.T14502C mutation, which causes the substitution of a conserved isoleucine for valine at position 58 in ND6 mRNA, may affect the mitochondrial respiratory chain functions. Biochemical analysis revealed that cell lines bearing both the m.A3243G and m.T14502C mutations exhibited greater reductions in ATP levels and increased ROS production compared with those carrying only the m.A3243G mutation. However, we did not find any mutations in the GJB2, GJB3 and GJB6 genes.

CONCLUSION

Our study indicated that mitochondrial diabetes is associated with the tRNA^Leu(UUR) A3243G and ND6 T14502C mutations.

Late onset of type 2 diabetes is associated with mitochondrial tRNATrp A5514G and tRNASer(AGY) C12237T mutations

Background

Mitochondrial dysfunctions caused by mitochondrial DNA (mtDNA) pathogenic mutations play putative roles in type 2 diabetes mellitus (T2DM) progression. But the underlying mechanism remains poorly understood.

Methods

A large Chinese family with maternally inherited diabetes and deafness (MIDD) underwent clinical, genetic, and molecular assessment. PCR and sequence analysis are carried out to detect mtDNA variants in affected family members, in addition, phylogenetic conservation analysis, haplogroup classification, and pathogenicity scoring system are performed. Moreover, the GJB2, GJB3, GJB6, and TRMU genes mutations are screened by PCR‐Sanger sequencing.

Results

Six of 18 matrilineal subjects manifested different clinical phenotypes of diabetes. The average age at onset of diabetic patients is 52 years. Screening for the entire mitochondrial genomes suggests the co‐existence of two possibly pathogenic mutations: tRNA^Trp A5514G and tRNA^Ser(AGY) C12237T, which belongs to East Asia haplogroup G2a. By molecular level, m.A5514G mutation resides at acceptor stem of tRNA^Trp (position 3), which is critical for steady‐state level of tRNA^Trp. Conversely, m.C12237T mutation occurs in the variable region of tRNA^Ser(AGY) (position 31), which creates a novel base‐pairing (11A‐31T). Thus, the mitochondrial dysfunctions caused by tRNA^Trp A5514G and tRNA^Ser(AGY) C12237T mutations, may be associated with T2DM in this pedigree. But we do not find any functional mutations in those nuclear genes.

Conclusion

Our findings suggest that m.A5514G and m.C12337T mutations are associated with T2DM, screening for mt‐tRNA mutations is useful for molecular diagnosis and prevention of mitochondrial diabetes.

Mutational Analysis of Mitochondrial tRNA Genes in 200 Patients with Type 2 Diabetes Mellitus

Objective

Previous studies showed that variants in mitochondrial DNA (mtDNA) are associated with type 2 diabetes mellitus (T2DM). However, the relationships between mitochondrial tRNA (mt-tRNA) variants and T2DM remain poorly understood.

Methods

In this study, we performed a mutational screening of 22 mt-tRNA genes in a cohort of 200 Han Chinese subjects with T2DM and 200 control subjects through PCR–Sanger sequencing. The identified mt-tRNA variants were assessed for their pathogenicity via the phylogenetic approach, structural and functional analysis. Furthermore, two Han Chinese pedigrees with maternally inherited diabetes and deafness (MIDD) were reported by clinical and genetic assessments.

Results

A total of 49 genetic variants in mt-tRNA genes were identified; among them, 31 variants (17 pathogenic/likely pathogenic) were absent in controls, located at extremely conserved nucleotides, may have potential structural and functional significance, thereby considered to be T2DM-associated variants. In addition, sequence analysis of entire mitochondrial genomes of the matrilineal relatives from two MIDD pedigrees revealed the occurrence of tRNA^Leu(UUR) A3243G and T3290C mutations, as well as sets of polymorphisms belonging to mitochondrial haplogroups F2 and D4. However, the lack of any functional variants in connexin 26 gene (GJB2) and tRNA 5-methylaminomethyl-2-thiouridylate (TRMU) suggested that nuclear genes may not play active roles in clinical expression of MIDD in these pedigrees.

Conclusion

Our data indicated that mt-tRNA variants were associated with T2DM, screening for mt-tRNA pathogenic mutations was recommended for early detection and prevention of mitochondrial diabetes.

Mitochondrial DNA variants modulate N-formylmethionine, proteostasis and risk of late-onset human diseases

Mitochondrial DNA (mtDNA) variants influence the risk of late-onset human diseases, but the reasons for this are poorly understood. Undertaking a hypothesis-free analysis of 5,689 blood-derived biomarkers with mtDNA variants in 16,220 healthy donors, here we show that variants defining mtDNA haplogroups Uk and H4 modulate the level of circulating N-formylmethionine (fMet), which initiates mitochondrial protein translation. In human cytoplasmic hybrid (cybrid) lines, fMet modulated both mitochondrial and cytosolic proteins on multiple levels, through transcription, post-translational modification and proteolysis by an N-degron pathway, abolishing known differences between mtDNA haplogroups. In a further 11,966 individuals, fMet levels contributed to all-cause mortality and the disease risk of several common cardiovascular disorders. Together, these findings indicate that fMet plays a key role in common age-related disease through pleiotropic effects on cell proteostasis.

Clinical and molecular features of two diabetes families carrying mitochondrial ND1 T3394C mutation

BACKGROUND

Mutations in mitochondrial DNA (mtDNA) are found to be associated with type 2 diabetes mellitus (T2DM). However, the molecular pathogenesis of these mutations in T2DM is still poorly understood.

METHODS

In this study, we report here the molecular features of two Han Chinese families with maternally transmitted T2DM. The matrilineal relatives are undergoing clinical, biochemical, genetic evaluations, and molecular analysis. Furthermore, the entire mitochondrial genomes of these matrilineal relatives are screened by PCR-Sanger sequencing.

RESULTS

The age at onset of T2DM of these participants varies from 28 to 71 years, with an average of 43 years. Molecular analysis of mitochondrial genomes identifies the existence of ND1 T3394C mutation in both families, together with sets of variants belonging to mitochondrial haplogroup Y2 and M9a. The m.T3394C mutation is localized at very conserved tyrosine at position 30 of ND1, may result the failure in ND1 mRNA metabolism, and lead to mitochondrial dysfunction. Moreover, sequence analysis of matrilineal relatives in Family 1 identifies the m.A14693G mutation which occurs in the TΨC-loop of tRNA^Glu (position 54), and is critical to the structural formation and stabilization of this tRNA. Thus, m.A14693G mutation may cause the impairment in tRNA metabolism, thereby worsens the mitochondrial dysfunction caused by ND1 T3394C mutation. However, no functional mtDNA variants are identified in Family 2 which suggest that mitochondrial haplogroup may not play an important role in diabetes expression.

CONCLUSIONS

Our study indicates that mitochondrial ND1 T3394C mutation is involved in the pathogenesis of maternally inherited T2DM in these families.

Molecular characterization of two Chinese pedigrees with maternally inherited hypertension

BACKGROUND

Mutations in mitochondrial tRNA (mt-tRNA) genes are associated with hypertension, although their pathogenic mechanisms remain poorly understood.

METHODS

In the present study, two Han Chinese families with maternally transmitted hypertension were interviewed. The mtDNA mutations of matrilineal relatives were screened by polymerase chain reaction-Sanger sequencing. Mitochondrial ATP, membrane potential and reactive oxygen species (ROS) were also analyzed in polymononuclear leukocytes carrying these mt-tRNA mutations. Additionally, the levels of oxidative stress-related biomarkers [malondialdehyde (MDA), superoxide dismutase (SOD) and glutathione peroxidase (GSH-Px) and 8-hydroxy-2-deoxyguanosine (8-OHdG)] were analyzed.

RESULTS

Nine of 13 adult matrilineal relatives of these pedigrees exhibited a wide range of severity of hypertension. The age at onset of hypertension was 30-62 years (average 46 years). Mutational screening of mitochondrial genomes revealed tRNA^Arg T10410C and T10454C mutations. Indeed, the m.T10454C and m.T10410C mutations occurred at conserved bases of TΨC-loop and acceptor arm of tRNA^Arg (positions 55 and 6), which are critical for tRNA^Arg post-transcriptional modification. Thus, the defects in tRNA modification may cause failure in tRNA metabolism, impairing mitochondrial translation. Biochemical analysis revealed that m.T10454C or m.T10410C mutation significantly reduced mitochondrial ATP and membrane potential and also increased ROS production in mutant cell lines (all p < 0.05). In addition, the levels of MDA and 8-OHdG in hypertensive patients markedly increased, whereas those of SOD and GSH-Px decreased (all p < 0.05).

CONCLUSIONS

These findings demonstrate that m.T10410C and m.T10454C mutations affect the structure and function of tRNA^Arg and consequently alter mitochondrial function and lead to oxidative stress, which are involved in the pathogenesis of maternally inherited hypertension.

Targeted Sequencing of Mitochondrial Genes Reveals Signatures of Molecular Adaptation in a Nearly Panmictic Small Pelagic Fish Species

Ongoing climatic changes, with predictable impacts on marine environmental conditions, are expected to trigger organismal responses. Recent evidence shows that, in some marine species, variation in mitochondrial genes involved in the aerobic conversion of oxygen into ATP at the cellular level correlate with gradients of sea surface temperature and gradients of dissolved oxygen. Here, we investigated the adaptive potential of the European sardine Sardina pilchardus populations offshore the Iberian Peninsula. We performed a seascape genetics approach that consisted of the high throughput sequencing of mitochondria’s ATP6, COI, CYTB and ND5 and five microsatellite loci on 96 individuals coupled with environmental information on sea surface temperature and dissolved oxygen across five sampling locations. Results show that, despite sardines forming a nearly panmictic population around Iberian Peninsula, haplotype frequency distribution can be explained by gradients of minimum sea surface temperature and dissolved oxygen. We further identified that the frequencies of the most common CYTB and ATP6 haplotypes negatively correlate with minimum sea surface temperature across the sampled area, suggestive of a signature of selection. With signatures of selection superimposed on highly connected populations, sardines may be able to follow environmental optima and shift their distribution northwards as a response to the increasing sea surface temperatures.

The roles of mitochondrial tRNA mutations in non-dystrophic myotonias

According a recent report by Heidari et al., a mutational screening for candidate pathogenic mitochondrial tRNA (mt-tRNA) mutations were performed in 45 Iranian patients with non-dystrophic myotonia (NDM) and 70 control subjects. Through PCR amplification and direct sequence analysis, nine mt-tRNA mutations were identified: tRNA^Met T4454C, tRNA^Trp A5568G, tRNA^Cys T5794C, tRNA^Arg A10438T and T10462C, tRNA^Leu(CUN) A12308G, tRNA^Thr A15907G, A15924G and G15928A. However, through the database searches and phylogenetic conservation analysis, we noticed that the tRNA^Thr A15924G, G15928A and tRNA^Leu(CUN) A12308G mutations should be classified 'pathogenic'. Thus, the roles of mt-tRNA mutations in clinical expression of NDM needed to be further experimentally addressed.

Ischemic Stroke Risk Associated with Mitochondrial Haplogroup F in the Asian Population

Mitochondrial dysfunction is involved in the pathogenesis of atherosclerosis, the primary risk factor for ischemic stroke. This study aims to explore the role of mitochondrial genomic variations in ischemic stroke, and to uncover the nuclear genes involved in this relationship. Eight hundred and thirty Taiwanese patients with a history of ischemic stroke and 966 normal controls were genotyped for their mitochondrial haplogroup (Mthapg). Cytoplasmic hybrid cells (cybrids) harboring different Mthapgs were used to observe functional differences under hypoxia-ischemia. RNA sequencing (RNASeq) was conducted to identify the particularly elevated mRNA. The patient study identified an association between Mthapg F1 and risk of ischemic stroke (OR 1.72:1.27-2.34, p = 0.001). The cellular study further demonstrated an impeded induction of hypoxic inducible factor 1α in the Mthapg F1 cybrid after hypoxia-ischemia. Additionally, the study demonstrated that Mthapg F cybrids were associated with an altered mitochondrial function, including decreased oxygen consumption, higher mitochondrial ROS production, and lower mitochondrial membrane potential. Mthapg F cybrids were also noted to be prone to inflammation, with increased expression of several inflammatory cytokines and elevated matrix metalloproteinase 9. The RNASeq identified significantly elevated expressions of angiopoietin-like 4 in Mthapg F1 cybrids after hypoxia-ischemia. Our study demonstrates an association between Mthapg F and susceptibility to ischemic stroke.

The Mitochondrial tRNAHis G12192A Mutation May Modulate the Clinical Expression of Deafness-Associated tRNAThr G15927A Mutation in a Chinese Pedigree

BACKGROUND

Mutations in mitochondrial tRNA (mt-tRNA) genes have been found to be associated with both syndromic and non-syndromic hearing impairment. However, the pathophysiology underlying mt-tRNA mutations in clinical expression of hearing loss remains poorly understood.

OBJECTIVE

The aim of this study was to explore the potential association between mttRNA mutations and hearing loss.

METHODS AND RESULTS

We reported here the molecular features of a pedigree with maternally transmitted non-syndromic hearing loss. Among 12 matrilineal relatives, five of them suffered variable degree of hearing impairment, but none of them had any medical history of using aminoglycosides antibiotics (AmAn). Genetic screening of the complete mitochondrial genomes from the matrilineal relatives identified the coexistence of mt-tRNAHis G12192A and mt-tRNAThr G15927A mutations, together with a set of polymorphisms belonging to human mitochondrial haplogroup B5b1b. Interestingly, the G12192A mutation occurred 2-bp from the 3' end of the TψC loop of mt-tRNAHis, which was evolutionarily conserved from various species. In addition, the well-known G15927A mutation, which disrupted the highly conserved C-G base-pairing at the anticodon stem of mt-tRNAThr, may lead to the failure in mt-tRNA metabolism. Furthermore, a significant decreased in ATP production and an increased ROS generation were observed in polymononuclear leukocytes (PMNs) which were isolated from the deaf patients carrying these mt-tRNA mutations, suggested that the G12192A and G15927A mutations may cause mitochondrial dysfunction that was responsible for deafness. However, the absence of any functional mutations/variants in GJB2, GJB3, GJB6 and TRMU genes suggested that the nuclear genes may not play important roles in the clinical expression of non-syndromic hearing loss in this family.

CONCLUSION

Our data indicated that mt-tRNAHis G12192A mutation may increase the penetrance and expressivity of deafness-associated m-tRNAThr G15927A mutation in this family.

NUMT Confounding Biases Mitochondrial Heteroplasmy Calls in Favor of the Reference Allele

Homology between mitochondrial DNA (mtDNA) and nuclear DNA of mitochondrial origin (nuMTs) causes confounding when aligning short sequence reads to the reference human genome, as the true sequence origin cannot be determined. Using a systematic in silico approach, we here report the impact of all potential mitochondrial variants on alignment accuracy and variant calling. A total of 49,707 possible mutations were introduced across the 16,569 bp reference mitochondrial genome (16,569 × 3 alternative alleles), one variant at-at-time. The resulting in silico fragmentation and alignment to the entire reference genome (GRCh38) revealed preferential mapping of mutated mitochondrial fragments to nuclear loci, as variants increased loci similarity to nuMTs, for a total of 807, 362, and 41 variants at 333, 144, and 27 positions when using 100, 150, and 300 bp single-end fragments. We subsequently modeled these affected variants at 50% heteroplasmy and carried out variant calling, observing bias in the reported allele frequencies in favor of the reference allele. Four variants (chrM:6023A, chrM:4456T, chrM:5147A, and chrM:7521A) including a possible hypertension factor, chrM:4456T, caused 100% loss of coverage at the mutated position (with all 100 bp single-end fragments aligning to homologous, nuclear positions instead of chrM), rendering these variants undetectable when aligning to the entire reference genome. Furthermore, four mitochondrial variants reported to be pathogenic were found to cause significant loss of coverage and select haplogroup-defining SNPs were shown to exacerbate the loss of coverage caused by surrounding variants. Increased fragment length and use of paired-end reads both improved alignment accuracy.

Sharing of heteroplasmies between human liver lobes varies across the mtDNA genome

Mitochondrial DNA (mtDNA) heteroplasmy (intra-individual variation) varies among different human tissues and increases with age, suggesting that the majority of mtDNA heteroplasmies are acquired, rather than inherited. However, the extent to which heteroplasmic sites are shared across a tissue remains an open question. We therefore investigated heteroplasmy in two liver samples (one from each primary lobe) from 83 Europeans, sampled at autopsy. Minor allele frequencies (MAF) at heteroplasmic sites were significantly correlated between the two liver samples from an individual, with significantly more sharing of heteroplasmic sites in the control region than in the non-control region. We show that this increased sharing for the control region cannot be explained by recent mutations at just a few specific heteroplasmic sites or by the possible presence of 7S DNA. Moreover, we carried out simulations to show that there is significantly more sharing than would be predicted from random genetic drift from a common progenitor cell. We also observe a significant excess of non-synonymous vs. synonymous heteroplasmies in the protein-coding region, but significantly more sharing of synonymous heteroplasmies. These contrasting patterns for the control vs. the non-control region, and for non-synonymous vs. synonymous heteroplasmies, suggest that selection plays a role in heteroplasmy sharing.

The mitochondrial tRNAAla 5587T>C and tRNALeu(CUN) 12280A>G mutations may be associated with hypertension in a Chinese family

Hypertension is a very common cardiovascular disorder, however, the molecular mechanism underlying this disease remains poorly understood. Recently, an increasing number of studies have demonstrated that mitochondrial (mt)DNA mutations serve important roles in the pathogenesis of hypertension. The current study reported the clinical and molecular characterization of a Chinese family with maternally inherited hypertension (the penetrance of hypertension was 71.4%). In addition, the entire mitochondrial transfer (mt-t)RNA genomes was amplified using a polymerase chain reaction (PCR) and identified through direct Sanger sequencing. Additionally, the mtDNA copy number in matrilineal relatives in this family was evaluated using quantitative PCR. The sequence analysis of the 22 mt-tRNA genes led to the identification of tRNA^Ala 5587T>C (thymine to cytosine) and tRNA^Leu(CUN) 12280A>G (adenine to guanine) mutations. Notably, the heteroplasmic 5587T>C mutation was located at the 3' end of tRNA^Ala (position 73), which is highly conserved from bacteria to human mitochondria. In addition, the 12280A>G mutation was revealed to occurs at the dihydrouridine loop of tRNA^Leu(CUN) (position 15) and may decrease the steady-state level of mt-tRNA. As a result, 5587T>C and 12280A>G mutations may lead to the failure of tRNAs metabolism and subsequently cause mitochondrial protein synthesis defects. Molecular analysis revealed that patients carrying the 5587T>C and 12280A>G mutations had a lower copy number of mtDNA compared with a control with hypertension, but without the mutations, suggesting that these mutations may cause mitochondrial dysfunctions that are responsible for hypertension. Therefore, mt-tRNA^Ala 5587T>C and tRNA^Leu(CUN) 12280A>G mutations may be involved in the pathogenesis of hypertension in this family.

A common pattern of DNase I footprinting throughout the human mtDNA unveils clues for a chromatin-like organization

Human mitochondrial DNA (mtDNA) is believed to lack chromatin and histones. Instead, it is coated solely by the transcription factor TFAM. We asked whether mtDNA packaging is more regulated than once thought. To address this, we analyzed DNase-seq experiments in 324 human cell types and found, for the first time, a pattern of 29 mtDNA Genomic footprinting (mt-DGF) sites shared by ∼90% of the samples. Their syntenic conservation in mouse DNase-seq experiments reflect selective constraints. Colocalization with known mtDNA regulatory elements, with G-quadruplex structures, in TFAM-poor sites (in HeLa cells) and with transcription pausing sites, suggest a functional regulatory role for such mt-DGFs. Altered mt-DGF pattern in interleukin 3-treated CD34⁺ cells, certain tissue differences, and significant prevalence change in fetal versus nonfetal samples, offer first clues to their physiological importance. Taken together, human mtDNA has a conserved protein-DNA organization, which is likely involved in mtDNA regulation.

The Mitonuclear Dimension of Neanderthal and Denisovan Ancestry in Modern Human Genomes

Some human populations interbred with Neanderthals and Denisovans, resulting in substantial contributions to modern-human genomes. Therefore, it is now possible to use genomic data to investigate mechanisms that shaped historical gene flow between humans and our closest hominin relatives. More generally, in eukaryotes, mitonuclear interactions have been argued to play a disproportionate role in generating reproductive isolation. There is no evidence of mtDNA introgression into modern human populations, which means that all introgressed nuclear alleles from archaic hominins must function on a modern-human mitochondrial background. Therefore, mitonuclear interactions are also potentially relevant to hominin evolution. We performed a detailed accounting of mtDNA divergence among hominin lineages and used population-genomic data to test the hypothesis that mitonuclear incompatibilities have preferentially restricted the introgression of nuclear genes with mitochondrial functions. We found a small but significant underrepresentation of introgressed Neanderthal alleles at such nuclear loci. Structural analyses of mitochondrial enzyme complexes revealed that these effects are unlikely to be mediated by physically interacting sites in mitochondrial and nuclear gene products. We did not detect any underrepresentation of introgressed Denisovan alleles at mitochondrial-targeted loci, but this may reflect reduced power because locus-specific estimates of Denisovan introgression are more conservative. Overall, we conclude that genes involved in mitochondrial function may have been subject to distinct selection pressures during the history of introgression from archaic hominins but that mitonuclear incompatibilities have had, at most, a small role in shaping genome-wide introgression patterns, perhaps because of limited functional divergence in mtDNA and interacting nuclear genes.

Cancer-specific SNPs originate from low-level heteroplasmic variants in human mitochondrial genomes of a matched cell line pair

Low-level mitochondrial heteroplasmy is a common phenomenon in both normal and cancer cells. Here, we investigate the link between low-level heteroplasmy and mitogenome mutations in a human breast cancer matched cell line by high-throughput sequencing. We identified 23 heteroplasmic sites, of which 15 were common between normal cells (Hs578Bst) and cancer cells (Hs578T). Most sites were clustered within the highly conserved Complex IV and ribosomal RNA genes. Two heteroplasmic variants in normal cells were found as fixed mutations in cancer cells. This indicates a positive selection of these variants in cancer cells. RNA-Seq analysis identified upregulated L-strand specific transcripts in cancer cells, which include three mitochondrial long non-coding RNA molecules. We hypothesize that this is due to two cancer cell-specific mutations in the control region.

Codon based co-occurrence network motifs in human mitochondria

The nucleotide polymorphism in the human mitochondrial genome (mtDNA) tolled by codon position bias plays an indispensable role in human population dispersion and expansion. Herein, genome-wide nucleotide co-occurrence networks were constructed using data comprised of five different geographical regions and around 3000 samples for each region. We developed a powerful network model to describe complex mitochondrial evolutionary patterns among codon and non-codon positions. We found evidence that the evolution of human mitochondria DNA is dominated by adaptive forces, particularly mutation and selection, which was supported by many previous studies. The diversity observed in the mtDNA was compared with mutations, co-occurring mutations, network motifs considering codon positions as causing agent. This comparison showed that long-range nucleotide co-occurrences have a large effect on genomic diversity. Most notably, codon motifs apparently underpinned the preferences among codon positions for co-evolution which is probably highly biased during the origin of the genetic code. Our analysis also showed that variable nucleotide positions of different human sub-populations implemented the independent mtDNA evolution to its geographical dispensation. Ergo, this study has provided both a network framework and a codon glance to investigate co-occurring genomic variations that are critical in underlying complex mitochondrial evolution.

Performance of the Early Access AmpliSeq™ Mitochondrial Panel with degraded DNA samples using the Ion Torrent™ platform

The Early Access AmpliSeq™ Mitochondrial Panel amplifies whole mitochondrial genomes for phylogenetic and kinship identifications, using Ion Torrent™ technology. There is currently limited information on its performance with degraded DNA, a common occurrence in forensic samples. This study evaluated the performance of the Panel with DNA samples degraded in vitro, to mimic conditions commonly found in forensic investigations. Purified DNA from five individuals was heat-treated at five time points each (125°C for 0, 30, 60, 120, and 240 min; total n = 25). The quality of DNA was assessed via a real-time DNA assay of genomic DNA and prepared for massively parallel sequencing on the Ion Torrent™ platform. Mitochondrial sequences were obtained for all samples and had an amplicon coverage averaging between 66X to 2803X. Most amplicons (157/162) displayed high coverages (452 ± 333X), while reads with less than 100X coverage were recorded in five amplicons only (90 ± 5X). Amplicon coverage was decreased with prolonged heating. At 72% strand balance, reads were well balanced between forward and reverse strands. Using a coverage threshold of ten reads per SNP, complete sequences were recovered in all samples and resolved kinship and, haplogroup relations. Additionally, the HV1 and HV2 regions of the reference and 240-min heat-treated samples (n = 10) were Sanger-sequenced for concordance. Overall, this study demonstrates the efficacy of a novel forensic Panel that recovers high quality mitochondrial sequences from degraded DNA samples.

Heteroplasmic shifts in tumor mitochondrial genomes reveal tissue-specific signals of relaxed and positive selection

Although mitochondrial genomes (mtDNA) accumulate elevated levels of mutations in cancer cells, the origin and functional impact of these mutations remain controversial. Here, we queried whole-genome sequence data from 1,916 patients across 24 cancer types to characterize patterns of mtDNA mutations and elucidate the selective constraints driving their fate. Given that mitochondrial genomes are polyploid, cells with advantageous levels of mtDNA mutations can be selected for depending on their cellular environment. Therefore, we tracked changes in per-cell abundances of mtDNA mutations from normal to tumor cells in the same patient. Tumor mitochondrial genomes show distinct mutational patterns and are disproportionately enriched for protein-altering changes. Moreover, protein-altering mtDNA variants that are initially present at low frequencies in normal cells preferentially expand in the altered tumor environment, suggesting selective advantage. We also perform these analyses with attention to the cancer's tissue of origin, which revealed tissue-specific differences in selective signals. The mitochondrial genomes in renal chromophobe and thyroid cancers show particularly strong signals of positive selection, indicated by higher proportions and per-cell abundances of truncating variants. Dramatic tumor- and tissue-specific variations in selective pressures suggest that cancer cells with advantageous levels of damaged mitochondrial genomes will selectively proliferate to facilitate the tumorigenic process.

Mitochondrial tRNALeu(UUR) C3275T, tRNAGln T4363C and tRNALys A8343G mutations may be associated with PCOS and metabolic syndrome

Polycystic ovary syndrome (PCOS) is a very prevalent endocrine disease affecting reproductive women. Clinically, patients with this disorder are more vulnerable to develop type 2 diabetes mellitus (T2DM), cardiovascular events, as well as metabolic syndrome (MetS). To date, the molecular mechanism underlying PCOS remains largely unknown. Previously, we showed that mitochondrial dysfunction caused by mitochondrial DNA (mtDNA) mutation was an important cause for PCOS. In the current study, we described the clinical and biochemical features of a three-generation pedigree with maternally transmitted MetS, combined with PCOS. A total of three matrilineal relatives exhibited MetS including obesity, high triglyceride (TG) and Hemoglobin A1c (HbA1c) levels, and hypertension. Whereas one patient from the third generation manifestated PCOS. Mutational analysis of the whole mitochondrial genes from the affected individuals identified a set of genetic variations belonging to East Asia haplogroup B4b1c. Among these variants, the homoplasmic C3275T mutation disrupted a highly evolutionary conserved base-pairing (28A-46C) on the variable region of tRNA^Leu(UUR), whereas the T4363C mutation created a new base-pairing (31T-37A) in the anticodon stem of tRNA^Gln, furthermore, the A8343G mutation occurred at the very conserved position of tRNA^Lys and may result the failure in mitochondrial tRNAs (mt-tRNAs) metabolism. Biochemical analysis revealed the deficiency in mitochondrial functions including lower levels of mitochondrial membrane potential (MMP), ATP production and mtDNA copy number, while a significantly increased reactive oxygen species (ROS) generation was observed in polymononuclear leukocytes (PMNs) from the individuals carrying these mt-tRNA mutations, suggesting that these mutations may cause mitochondrial dysfunction that was responsible for the clinical phenotypes. Taken together, our data indicated that mt-tRNA mutations were associated with MetS and PCOS in this family, which shaded additional light into the pathophysiology of PCOS that were manifestated by mitochondrial dysfunction.

Disease-associated mitochondrial mutations and the evolution of primate mitogenomes

Several human diseases have been associated with mutations in mitochondrial genes comprising a set of confirmed and reported mutations according to the MITOMAP database. An analysis of complete mitogenomes across 139 primate species showed that most confirmed disease-associated mutations occurred in aligned codon positions and gene regions under strong purifying selection resulting in a strong evolutionary conservation. Only two confirmed variants (7.1%), coding for the same amino acids accounting for severe human diseases, were identified without apparent pathogenicity in non-human primates, like the closely related Bornean orangutan. Conversely, reported disease-associated mutations were not especially concentrated in conserved codon positions, and a large fraction of them occurred in highly variable ones. Additionally, 88 (45.8%) of reported mutations showed similar variants in several non-human primates and some of them have been present in extinct species of the genus Homo. Considering that recurrent mutations leading to persistent variants throughout the evolutionary diversification of primates are less likely to be severely damaging to fitness, we suggest that these 88 mutations are less likely to be pathogenic. Conversely, 69 (35.9%) of reported disease-associated mutations occurred in extremely conserved aligned codon positions which makes them more likely to damage the primate mitochondrial physiology.

MtDNA meta-analysis reveals both phenotype specificity and allele heterogeneity: a model for differential association

Human mtDNA genetic variants have traditionally been considered markers for ancient population migrations. However, during the past three decades, these variants have been associated with altered susceptibility to various phenotypes, thus supporting their importance for human health. Nevertheless, mtDNA disease association has frequently been supported only in certain populations, due either to population stratification or differential epistatic compensations among populations. To partially overcome these obstacles, we performed meta-analysis of the multiple mtDNA association studies conducted until 2016, encompassing 53,975 patients and 63,323 controls. Our findings support the association of mtDNA haplogroups and recurrent variants with specific phenotypes such as Parkinson’s disease, type 2 diabetes, longevity, and breast cancer. Strikingly, our assessment of mtDNA variants’ involvement with multiple phenotypes revealed significant impact for Caucasian haplogroups H, J, and K. Therefore, ancient mtDNA variants could be divided into those that affect specific phenotypes, versus others with a general impact on phenotype combinations. We suggest that the mtDNA could serve as a model for phenotype specificity versus allele heterogeneity.

The genomic landscape of evolutionary convergence in mammals, birds and reptiles

Many lineage-defining (nodal) mutations possess high functionality. However, differentiating adaptive nodal mutations from those that are functionally compensated remains challenging. To address this challenge, we identified functional nodal mutations (fNMs) in ~3,400 nuclear DNA (nDNA) and 4 mitochondrial DNA (mtDNA) protein structures from 91 and 1,003 species, respectively, representing the entire mammalian, bird and reptile phylogeny. A screen for candidate compensatory mutations among co-occurring amino acid changes in close structural proximity revealed that such compensated fNMs encompass 37% and 27% of the mtDNA and nDNA datasets, respectively. Analysis of the remaining (non-compensated) mutations, which are enriched for adaptive mutations, showed that birds and mammals share most such recurrent fNMs (N = 51). Among the latter, we discovered mutations in thermoregulation-related genes. These represent the best candidates to explain the molecular basis of convergent body thermoregulation in birds and mammals. Our analysis reveals the landscape of possible mutational compensation and convergence in amniote phylogeny.

Initiation of mtDNA transcription is followed by pausing, and diverges across human cell types and during evolution

Mitochondrial DNA (mtDNA) genes are long known to be cotranscribed in polycistrones, yet it remains impossible to study nascent mtDNA transcripts quantitatively in vivo using existing tools. To this end, we used deep sequencing (GRO-seq and PRO-seq) and analyzed nascent mtDNA-encoded RNA transcripts in diverse human cell lines and metazoan organisms. Surprisingly, accurate detection of human mtDNA transcription initiation sites (TISs) in the heavy and light strands revealed a novel conserved transcription pausing site near the light-strand TIS. This pausing site correlated with the presence of a bacterial pausing sequence motif, with reduced SNP density, and with a DNase footprinting signal in all tested cells. Its location within conserved sequence block 3 (CSBIII), just upstream of the known transcription–replication transition point, suggests involvement in such transition. Analysis of nonhuman organisms enabled de novo mtDNA sequence assembly, as well as detection of previously unknown mtDNA TIS, pausing, and transcription termination sites with unprecedented accuracy. Whereas mammals (Pan troglodytes, Macaca mulatta, Rattus norvegicus, and Mus musculus) showed a human-like mtDNA transcription pattern, the invertebrate pattern (Drosophila melanogaster and Caenorhabditis elegans) profoundly diverged. Our approach paves the path toward in vivo, quantitative, reference sequence-free analysis of mtDNA transcription in all eukaryotes.

Late-onset Mitochondrial Cardiomyopathy Triggered by Anticancer Treatment

We report the case of a 62-year-old woman with a history of bilateral hearing impairment, who developed mitochondrial cardiomyopathy after chemotherapy. The patient underwent postoperative cisplatin chemotherapy after the surgical treatment of cervical cancer. The systolic function of her left ventricle decreased significantly. A tissue examination of the left ventricle revealed mitochondrial cardiomyopathy. Genetic testing revealed mutations in mitochondrial 3,243 A→G. Nine hundred fifty-five individual mutations were identified by next-generation sequencing. Since cardiovascular complications are the second leading cause of morbidity and mortality in patients undergoing cancer treatment, mitochondrial cardiomyopathy should be considered a potential cause of heart failure.

Mitochondrial variations in the MT-ND4 and MT-TL1 genes are associated with male infertility

Mitochondrial gene mutations have been reported to be associated with sperm motility and the quality of semen. The aim of this study was to investigate whether the two mitochondrial genes (MT-ND4 and MT-TL1) are involved in Chinese male infertility. A total of 97 asthenospermia patients and 80 fertile controls were recruited in this case-control study. Genomic DNA were extracted from the sperm of all participants. Two mitochondrial DNA genes (MT-ND4 and MT-TL1) were amplified by using polymerase chain reaction (PCR) with the gene-specific primers and sequenced on an ABI 3730XL DNA sequencer. For the MT-ND4 gene, we found a total of 64 and 54 nucleotide substitutions in patients and controls, respectively, with no discrepancy in the mutation rates (66.0% vs. 67.5%, p>0.05). However, one mutation (g.11084A>G, p.T109A) leading to an amino acid substitution in a highly conserved residue and predicted to be deleterious was detected only in the cases. For another gene MT-TL1, a novel mutation (g.3263C>T) near the anticodon TAA was identified in an asthenospermia patient and was absent from normal controls. However, the mutation positions in the cases varied from the controls and one highly conserved mutation (g.11084A>G, p.T109A) which was not found in the controls and probably caused damage to the protein structure might contribute to asthenospermia. For another gene MT-TL1, a highly conservative novel mutation which is located closely next to the anticodon also might contribute to asthenospermia. Our result suggests that the MT-ND4 and MT-TL1 genes might be associated with Chinese male infertility.

ABBREVIATIONS

MT-ND4: mitochondrially encoded NADH dehydrogenase 4; MT-TL1: mitochondrially encoded tRNA leucine 1 (UUA/G); PCR: polymerase chain reaction; OXPHOS: mitochondrial oxidative phosphorylation; ATP: adenosine triphosphate; mtDNA: mitochondrial DNA; SNPs: single nucleotide substitutions; AD: alzheimer's disease; PD: parkinson's disease; MELAS: mitochondrial encephalomyopathy with lactic acidosis and stroke-like episodes; ROS: reactive oxygen species.

Ancient Out-of-Africa Mitochondrial DNA Variants Associate with Distinct Mitochondrial Gene Expression Patterns

Mitochondrial DNA (mtDNA) variants have been traditionally used as markers to trace ancient population migrations. Although experiments relying on model organisms and cytoplasmic hybrids, as well as disease association studies, have served to underline the functionality of certain mtDNA SNPs, only little is known of the regulatory impact of ancient mtDNA variants, especially in terms of gene expression. By analyzing RNA-seq data of 454 lymphoblast cell lines from the 1000 Genomes Project, we found that mtDNA variants defining the most common African genetic background, the L haplogroup, exhibit a distinct overall mtDNA gene expression pattern, which was independent of mtDNA copy numbers. Secondly, intra-population analysis revealed subtle, yet significant, expression differences in four tRNA genes. Strikingly, the more prominent African mtDNA gene expression pattern best correlated with the expression of nuclear DNA-encoded RNA-binding proteins, and with SNPs within the mitochondrial RNA-binding proteins PTCD1 and MRPS7. Our results thus support the concept of an ancient regulatory transition of mtDNA-encoded genes as humans left Africa to populate the rest of the world.

The mitochondrion is an organelle found in all cells of our body and plays a significant role in the energy and heat production. This is the only organelle in animal cells harboring its own genome outside of the nucleus. Mitochondrial DNA (mtDNA) variants have been traditionally used as neutral markers to trace ancient population migrations. As a result, the functional impact of human mtDNA population variants on gene regulation is poorly understood. To address this question, we analyzed available data of mtDNA gene expression pattern in a large group of individuals (454) from diverse human populations. Here, we show for the first time that the ancient migration of humans out of Africa correlated with differences in mitochondrial gene expression patterns, and could be explained by the activity of certain RNA-binding proteins. These findings suggest a major mitochondrial regulatory transition, as humans left Africa to populate the rest of the world.

Analysis of mitochondrial alterations in Brazilian patients with sensorineural hearing loss using MALDI-TOF mass spectrometry

BACKGROUND

Mutations in the mitochondrial DNA (mtDNA) have been associated with aminoglycoside-induced and nonsyndromic deafness in different populations. In the present study, we investigated the contribution of mutations in mitochondrial genes to the etiology of hearing loss in a Brazilian sample.

METHODS

Using mass spectrometry genotyping technology, combined with direct sequencing, 50 alterations previously described in 14 mitochondrial genes were screened in 152 patients with sensorineural hearing loss and in104 normal hearing controls.

RESULTS

Fifteen known mitochondrial alterations were detected (G709A, A735G, A827G, G988A, A1555G, T4363C, T5628C, T5655C, G5821A, C7462T, G8363A, T10454C, G12236A, T1291C, G15927A). Pathogenic mutations in MT-RNR1 and MT-TK genes were detected in 3 % (5/152) of the patients with hearing loss.

CONCLUSIONS

This study contributed to show the spectrum of mitochondrial variants in Brazilian patients with hearing loss. Frequency of A1555G was relatively high (2.6 %), indicating that this mutation is an important cause of hearing loss in our population. This work reports for the first time the investigation and the detection of the tRNA(Lys) G8363A mutation in Brazilian patients with maternally inherited sensorineural hearing loss.

Mitochondrial cytochrome c oxidase deficiency

As with other mitochondrial respiratory chain components, marked clinical and genetic heterogeneity is observed in patients with a cytochrome c oxidase deficiency. This constitutes a considerable diagnostic challenge and raises a number of puzzling questions. So far, pathological mutations have been reported in more than 30 genes, in both mitochondrial and nuclear DNA, affecting either structural subunits of the enzyme or proteins involved in its biogenesis. In this review, we discuss the possible causes of the discrepancy between the spectacular advances made in the identification of the molecular bases of cytochrome oxidase deficiency and the lack of any efficient treatment in diseases resulting from such deficiencies. This brings back many unsolved questions related to the frequent delay of clinical manifestation, variable course and severity, and tissue-involvement often associated with these diseases. In this context, we stress the importance of studying different models of these diseases, but also discuss the limitations encountered in most available disease models. In the future, with the possible exception of replacement therapy using genes, cells or organs, a better understanding of underlying mechanism(s) of these mitochondrial diseases is presumably required to develop efficient therapy.

Mitochondrial Involvement in Vertebrate Speciation? The Case of Mito-nuclear Genetic Divergence in Chameleons

Compatibility between the nuclear (nDNA) and mitochondrial (mtDNA) genomes is important for organismal health. However, its significance for major evolutionary processes such as speciation is unclear, especially in vertebrates. We previously identified a sharp mtDNA-specific sequence divergence between morphologically indistinguishable chameleon populations (Chamaeleo chamaeleon recticrista) across an ancient Israeli marine barrier (Jezreel Valley). Because mtDNA introgression and gender-based dispersal were ruled out, we hypothesized that mtDNA spatial division was maintained by mito-nuclear functional compensation. Here, we studied RNA-seq generated from each of ten chameleons representing the north and south populations and identified candidate nonsynonymous substitutions (NSSs) matching the mtDNA spatial distribution. The most prominent NSS occurred in 14 nDNA-encoded mitochondrial proteins. Increased chameleon sample size (N = 70) confirmed the geographic differentiation in POLRMT, NDUFA5, ACO1, LYRM4, MARS2, and ACAD9. Structural and functionality evaluation of these NSSs revealed high functionality. Mathematical modeling suggested that this mito-nuclear spatial divergence is consistent with hybrid breakdown. We conclude that our presented evidence and mathematical model underline mito-nuclear interactions as a likely role player in incipient speciation in vertebrates.

Are mutagenic non D-loop direct repeat motifs in mitochondrial DNA under a negative selection pressure?

Non D-loop direct repeats (DRs) in mitochondrial DNA (mtDNA) have been commonly implicated in the mutagenesis of mtDNA deletions associated with neuromuscular disease and ageing. Further, these DRs have been hypothesized to put a constraint on the lifespan of mammals and are under a negative selection pressure. Using a compendium of 294 mammalian mtDNA, we re-examined the relationship between species lifespan and the mutagenicity of such DRs. Contradicting the prevailing hypotheses, we found no significant evidence that long-lived mammals possess fewer mutagenic DRs than short-lived mammals. By comparing DR counts in human mtDNA with those in selectively randomized sequences, we also showed that the number of DRs in human mtDNA is primarily determined by global mtDNA properties, such as the bias in synonymous codon usage (SCU) and nucleotide composition. We found that SCU bias in mtDNA positively correlates with DR counts, where repeated usage of a subset of codons leads to more frequent DR occurrences. While bias in SCU and nucleotide composition has been attributed to nucleotide mutational bias, mammalian mtDNA still exhibit higher SCU bias and DR counts than expected from such mutational bias, suggesting a lack of negative selection against non D-loop DRs.

Underlying Data for Sequencing the Mitochondrial Genome with the Massively Parallel Sequencing Platform Ion Torrent™ PGM™

BACKGROUND

Massively parallel sequencing (MPS) technologies have the capacity to sequence targeted regions or whole genomes of multiple nucleic acid samples with high coverage by sequencing millions of DNA fragments simultaneously. Compared with Sanger sequencing, MPS also can reduce labor and cost on a per nucleotide basis and indeed on a per sample basis. In this study, whole genomes of human mitochondria (mtGenome) were sequenced on the Personal Genome Machine (PGMTM) (Life Technologies, San Francisco, CA), the out data were assessed, and the results were compared with data previously generated on the MiSeqTM (Illumina, San Diego, CA). The objectives of this paper were to determine the feasibility, accuracy, and reliability of sequence data obtained from the PGM.

RESULTS

24 samples were multiplexed (in groups of six) and sequenced on the at least 10 megabase throughput 314 chip. The depth of coverage pattern was similar among all 24 samples; however the coverage across the genome varied. For strand bias, the average ratio of coverage between the forward and reverse strands at each nucleotide position indicated that two-thirds of the positions of the genome had ratios that were greater than 0.5. A few sites had more extreme strand bias. Another observation was that 156 positions had a false deletion rate greater than 0.15 in one or more individuals. There were 31-98 (SNP) mtGenome variants observed per sample for the 24 samples analyzed. The total 1237 (SNP) variants were concordant between the results from the PGM and MiSeq. The quality scores for haplogroup assignment for all 24 samples ranged between 88.8%-100%.

CONCLUSIONS

In this study, mtDNA sequence data generated from the PGM were analyzed and the output evaluated. Depth of coverage variation and strand bias were identified but generally were infrequent and did not impact reliability of variant calls. Multiplexing of samples was demonstrated which can improve throughput and reduce cost per sample analyzed. Overall, the results of this study, based on orthogonal concordance testing and phylogenetic scrutiny, supported that whole mtGenome sequence data with high accuracy can be obtained using the PGM platform.

Mito-nuclear co-evolution: the positive and negative sides of functional ancient mutations

Most cell functions are carried out by interacting factors, thus underlying the functional importance of genetic interactions between genes, termed epistasis. Epistasis could be under strong selective pressures especially in conditions where the mutation rate of one of the interacting partners notably differs from the other. Accordingly, the order of magnitude higher mitochondrial DNA (mtDNA) mutation rate as compared to the nuclear DNA (nDNA) of all tested animals, should influence systems involving mitochondrial-nuclear (mito-nuclear) interactions. Such is the case of the energy producing oxidative phosphorylation (OXPHOS) and mitochondrial translational machineries which are comprised of factors encoded by both the mtDNA and the nDNA. Additionally, the mitochondrial RNA transcription and mtDNA replication systems are operated by nDNA-encoded proteins that bind mtDNA regulatory elements. As these systems are central to cell life there is strong selection toward mito-nuclear co-evolution to maintain their function. However, it is unclear whether (A) mito-nuclear co-evolution befalls only to retain mitochondrial functions during evolution or, also, (B) serves as an adaptive tool to adjust for the evolving energetic demands as species' complexity increases. As the first step to answer these questions we discuss evidence of both negative and adaptive (positive) selection acting on the mtDNA and nDNA-encoded genes and the effect of both types of selection on mito-nuclear interacting factors. Emphasis is given to the crucial role of recurrent ancient (nodal) mutations in such selective events. We apply this point-of-view to the three available types of mito-nuclear co-evolution: protein-protein (within the OXPHOS system), protein-RNA (mainly within the mitochondrial ribosome), and protein-DNA (at the mitochondrial replication and transcription machineries).

A genetic basis for coronary artery disease

CAD and cancer account for over one-half of all deaths in the world. It is claimed that the 21st century is the last century for CAD. This is, in part, because CAD is preventable based on randomized, placebo-controlled clinical trials, which show modifying known risk factors such as cholesterol is associated consistently with 40-60% reduction in morbidity and mortality from CAD. Comprehensive prevention will require modifying genetic risk factors that are claimed to account for 40-60% of predisposition to CAD. The 21st century is meeting this challenge with over 50 genetic risk variants discovered and replicated in large genome-wide association studies involving over 200,000 cases and controls. Similarly, 157 genetic variants have been discovered that regulate plasma lipids including, LDL-C, HDL-C, triglycerides, and total cholesterol. A major finding from these studies is that only 15 of the 50 genetic variants for CAD act through known risk factors. Hence, the pathogenesis of CAD in addition to cholesterol and other known risk factors is due to various other factors, many of which remain unknown. Secondly, genes regulating the plasma triglyceride levels are strongly associated with the pathogenesis of CAD. Thirdly, Mendelian randomization studies show no protection from genes that increase plasma HDL cholesterol. This is contrary to current opinion. These genetic risk variants have provided new targets for the development of novel therapies to prevent CAD. Already a new and potent drug has been developed targeting PCSK9, which is in phase 3 clinical trials and shows great efficacy and safety for prevention of CAD. The 21st century is looking very bright for the prevention of CAD.