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

Ancestry affects the transcription of small mitochondrial RNAs in human lymphocytes

Mitochondrial mutations have been linked to changes in phenotypes such as fertility or longevity, however, these changes have been often inconsistent across populations for unknown reasons. A hypothesis that could explain this inconsistency is that some still uncharacterized mitochondrial products are mediating the phenotypic changes across populations. It has been hypothesized that one such product could be the small RNAs encoded in the mitochondrial genome, thus this work will provide new evidence for their existence and function. By using data from the 1000 genome project and knowledge from previously characterized nuclear small RNAs, this study found that 10 small RNAs encoded in tRNA fragments are consistently expressed in 450 individuals from five different populations. Furthermore, this study demonstrated that the expression of some small mitochondrial RNAs is different in individuals of African ancestry, similar to what was observed before in nuclear and mitochondria mRNAs. Lastly, we investigate the causes behind these differences in expression, showing that at least one of the mt-tRFs might be regulated by TRMT10B. The analyses presented in this work further support the small mitochondrial RNAs as functional molecules, and their population-specific expression supports the hypothesis that they act as a mediator between the nucleus and mitochondria differently across populations.

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.

Exploring the Effects of Mitonuclear Interactions on Mitochondrial DNA Gene Expression in Humans

Most mitochondrial protein complexes include both nuclear and mitochondrial gene products, which coevolved to work together. This coevolution can be disrupted due to disparity in genetic ancestry between the nuclear and mitochondrial genomes in recently admixed populations. Such mitonuclear DNA discordance might result in phenotypic effects. Several nuclear-encoded proteins regulate expression of mitochondrial DNA (mtDNA) genes. We hypothesized that mitonuclear DNA discordance affects expression of genes encoded by mtDNA. To test this, we utilized the data from the GTEx project, which contains expression levels for ∼100 African Americans and >600 European Americans. The varying proportion of African and European ancestry in recently admixed African Americans provides a range of mitonuclear discordance values, which can be correlated with mtDNA gene expression levels (adjusted for age and ischemic time). In contrast, European Americans did not undergo recent admixture. We demonstrated that, for most mtDNA protein-coding genes, expression levels in energetically-demanding tissues were lower in African Americans than in European Americans. Furthermore, gene expression levels were lower in individuals with higher mitonuclear discordance, independent of population. Moreover, we found a negative correlation between mtDNA gene expression and mitonuclear discordance. In African Americans, the average value of African ancestry was higher for nuclear-encoded mitochondrial than non-mitochondrial genes, facilitating a match in ancestry with the mtDNA and more optimal interactions. These results represent an example of a phenotypic effect of mitonuclear discordance on human admixed populations, and have potential biomedical applications.

Genome-wide analysis of mitochondrial DNA copy number reveals loci implicated in nucleotide metabolism, platelet activation, and megakaryocyte proliferation

Mitochondrial DNA copy number (mtDNA-CN) measured from blood specimens is a minimally invasive marker of mitochondrial function that exhibits both inter-individual and intercellular variation. To identify genes involved in regulating mitochondrial function, we performed a genome-wide association study (GWAS) in 465,809 White individuals from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium and the UK Biobank (UKB). We identified 133 SNPs with statistically significant, independent effects associated with mtDNA-CN across 100 loci. A combination of fine-mapping, variant annotation, and co-localization analyses was used to prioritize genes within each of the 133 independent sites. Putative causal genes were enriched for known mitochondrial DNA depletion syndromes (p = 3.09 × 10-15) and the gene ontology (GO) terms for mtDNA metabolism (p = 1.43 × 10-8) and mtDNA replication (p = 1.2 × 10-7). A clustering approach leveraged pleiotropy between mtDNA-CN associated SNPs and 41 mtDNA-CN associated phenotypes to identify functional domains, revealing three distinct groups, including platelet activation, megakaryocyte proliferation, and mtDNA metabolism. Finally, using mitochondrial SNPs, we establish causal relationships between mitochondrial function and a variety of blood cell-related traits, kidney function, liver function and overall (p = 0.044) and non-cancer mortality (p = 6.56 × 10-4).

Immune system cells from COVID-19 patients display compromised mitochondrial-nuclear expression co-regulation and rewiring toward glycolysis

Mitochondria are pivotal for bioenergetics, as well as in cellular response to viral infections. Nevertheless, their role in COVID-19 was largely overlooked. Here, we analyzed available bulk RNA-seq datasets from COVID-19 patients and corresponding healthy controls (three blood datasets, N = 48 healthy, 119 patients; two respiratory tract datasets, N = 157 healthy, 524 patients). We found significantly reduced mtDNA gene expression in blood, but not in respiratory tract samples from patients. Next, analysis of eight single-cells RNA-seq datasets from peripheral blood mononuclear cells, nasopharyngeal samples, and Bronchoalveolar lavage fluid (N = 1,192,243 cells), revealed significantly reduced mtDNA gene expression especially in immune system cells from patients. This is associated with elevated expression of nuclear DNA-encoded OXPHOS subunits, suggesting compromised mitochondrial-nuclear co-regulation. This, together with elevated expression of ROS-response genes and glycolysis enzymes in patients, suggest rewiring toward glycolysis, thus generating beneficial conditions for SARS-CoV-2 replication. Our findings underline the centrality of mitochondrial dysfunction in COVID-19.

•

mtDNA gene expression is downregulated in COVID-19 blood, but not in respiratory tract

•

Decreased mtDNA gene expression disrupts mito-nuclear coordination

•

mtDNA is downregulated and rewired toward glycolysis especially in immune system cells

•

Mitochondrial dysfunction is central to the etiology of COVID19

Emerging methods for and novel insights gained by absolute quantification of mitochondrial DNA copy number and its clinical applications

The past thirty years have seen a surge in interest in pathophysiological roles of mitochondria, and the accurate quantification of mitochondrial DNA copy number (mCN) in cells and tissue samples is a fundamental aspect of assessing changes in mitochondrial health and biogenesis. Quantification of mCN between studies is surprisingly variable due to a combination of physiological variability and diverse protocols being used to measure this endpoint. The advent of novel methods to quantify nucleic acids like digital polymerase chain reaction (dPCR) and high throughput sequencing offer the ability to measure absolute values of mCN. We conducted an in-depth survey of articles published between 1969 -- 2020 to create an overview of mCN values, to assess consensus values of tissue-specific mCN, and to evaluate consistency between methods of assessing mCN. We identify best practices for methods used to assess mCN, and we address the impact of using specific loci on the mitochondrial genome to determine mCN. Current data suggest that clinical measurement of mCN can provide diagnostic and prognostic value in a range of diseases and health conditions, with emphasis on cancer and cardiovascular disease, and the advent of means to measure absolute mCN should improve future clinical applications of mCN measurements.

Possible role of the mitochondrial genome in the pathogenesis of autosomal dominant polycystic kidney disease

Autosomal dominant polycystic kidney disease (ADPKD) is the most common monogenic renal disease in adults and is due to heterozygous germ line variants in either PKD1, PKD2 or rarely other genes. It is characterized by marked intra-familial disease variability suggesting that other genetic and/or environmental factors are involved in determining the lifetime course ADPKD. Recently, research indicates that polycystin-mediated mitochondrial dysfunction and metabolic re-programming contributes to the progression of ADPKD. Although biochemical abnormalities have gained the most interest, variants in the mitochondrial genome could be one of the mechanisms underlying the phenotypic variability in ADPKD. This narrative review aims to evaluate the role of the mitochondrial genome in the pathogenesis of APDKD.

Genetic basis of mitochondrial diseases

Mitochondrial disorders are monogenic disorders characterized by a defect in oxidative phosphorylation and caused by pathogenic variants in one of over 340 different genes. The implementation of whole-exome sequencing has led to a revolution in their diagnosis, duplicated the number of associated disease genes, and significantly increased the diagnosed fraction. However, the genetic etiology of a substantial fraction of patients exhibiting mitochondrial disorders remains unknown, highlighting limitations in variant detection and interpretation, which calls for improved computational and DNA sequencing methods, as well as the addition of OMICS tools. More intriguingly, this also suggests that some pathogenic variants lie outside of the protein-coding genes and that the mechanisms beyond the Mendelian inheritance and the mtDNA are of relevance. This review covers the current status of the genetic basis of mitochondrial diseases, discusses current challenges and perspectives, and explores the contribution of factors beyond the protein-coding regions and monogenic inheritance in the expansion of the genetic spectrum of disease.

Nuclear genome-wide associations with mitochondrial heteroplasmy

The role of the nuclear genome in maintaining the stability of the mitochondrial genome (mtDNA) is incompletely known. mtDNA sequence variants can exist in a state of heteroplasmy, which denotes the coexistence of organellar genomes with different sequences. Heteroplasmic variants that impair mitochondrial capacity cause disease, and the state of heteroplasmy itself is deleterious. However, mitochondrial heteroplasmy may provide an intermediate state in the emergence of novel mitochondrial haplogroups. We used genome-wide genotyping data from 982,072 European ancestry individuals to evaluate variation in mitochondrial heteroplasmy and to identify the regions of the nuclear genome that affect it. Age, sex, and mitochondrial haplogroup were associated with the extent of heteroplasmy. GWAS identified 20 loci for heteroplasmy that exceeded genome-wide significance. This included a region overlapping mitochondrial transcription factor A (TFAM), which has multiple roles in mtDNA packaging, replication, and transcription. These results show that mitochondrial heteroplasmy has a heritable nuclear component.

Mitochondrial gene expression in single cells shape pancreatic beta cells' sub-populations and explain variation in insulin pathway

Mitochondrial gene expression is pivotal to cell metabolism. Nevertheless, it is unknown whether it diverges within a given cell type. Here, we analysed single-cell RNA-seq experiments from human pancreatic alpha (N = 3471) and beta cells (N = 1989), as well as mouse beta cells (N = 1094). Cluster analysis revealed two distinct human beta cells populations, which diverged by mitochondrial (mtDNA) and nuclear DNA (nDNA)-encoded oxidative phosphorylation (OXPHOS) gene expression in healthy and diabetic individuals, and in newborn but not in adult mice. Insulin gene expression was elevated in beta cells with higher mtDNA gene expression in humans and in young mice. Such human beta cell populations also diverged in mitochondrial RNA mutational repertoire, and in their selective signature, thus implying the existence of two previously overlooked distinct and conserved beta cell populations. While applying our approach to human alpha cells, two sub-populations of cells were identified which diverged in mtDNA gene expression, yet these cellular populations did not consistently diverge in nDNA OXPHOS genes expression, nor did they correlate with the expression of glucagon, the hallmark of alpha cells. Thus, pancreatic beta cells within an individual are divided into distinct groups with unique metabolic-mitochondrial signature.

Mitochondrial DNA haplogroups affect physical performances in Han older adults: an 8-year follow-up prospective cohort study

AIM

The objective of this study was to explore the association between mitochondrial DNA (mtDNA) haplogroups and physical performances in Han older adults.

METHODS

This study was an 8-year follow-up prospective cohort study. A total of 104 Han older adults completed the measurements of the 6-min walk test, grip strength and mitochondrial DNA sequencing. The mtDNA haplogroups were classified by using HaploGrep2 software. We used the mixed model to analyze the longitudinal effect of mtDNA haplogroups on physical performance tests among three waves of data.

RESULTS

The mean age at wave 3 among men and women were 78.3 and 77.2 years, respectively. The overall proportions of mtDNA haplogroups in this study population was 26.9% F, 21.2% M, 15.4% R, 14.4% D, 8.7% B and others. After adjusted for age, sex, ethnicity, body mass index and exercise, the interaction of mtDNA haplogroup M and waves significantly affected the 6-min walk distance. Notably, the adjusted mean of the 6-min walk distance among the group of mtDNA haplogroup M was significantly lower than other haplogroups at wave 3. The adjusted mean of grip strength among the group of mtDNA haplogroup R was significantly higher than other haplogroups at wave 3.

CONCLUSIONS

This finding suggests that mtDNA haplogroups might have effects on the 6-min walk test and grip strength in Han older adults, although studies of the physical performance of older adults with larger sample sizes are necessary to further substantiate these findings. Geriatr Gerontol Int 2021; 21: 166-171.

Association between mitochondrial genetic variation and breast cancer risk: The Multiethnic Cohort

Background

The mitochondrial genome encodes for thirty-seven proteins, among them thirteen are essential for the oxidative phosphorylation (OXPHOS) system. Inherited variation in mitochondrial genes may influence cancer development through changes in mitochondrial proteins, altering the OXPHOS process and promoting the production of reactive oxidative species.

Methods

To investigate the association between mitochondrial genetic variation and breast cancer risk, we tested 314 mitochondrial SNPs (mtSNPs), capturing four complexes of the mitochondrial OXPHOS pathway and mtSNP groupings for rRNA and tRNA, in 2,723 breast cancer cases and 3,260 controls from the Multiethnic Cohort Study.

Results

We examined the collective set of 314 mtSNPs as well as subsets of mtSNPs grouped by mitochondrial OXPHOS pathway, complexes, and genes, using the sequence kernel association test and adjusting for age, sex, and principal components of global ancestry. We also tested haplogroup associations using unconditional logistic regression and adjusting for the same covariates. Stratified analyses were conducted by self-reported maternal race/ethnicity. No significant mitochondrial OXPHOS pathway, gene, and haplogroup associations were observed in African Americans, Asian Americans, Latinos, and Native Hawaiians. In European Americans, a global test of all genetic variants of the mitochondrial genome identified an association with breast cancer risk (P = 0.017, q = 0.102). In mtSNP-subset analysis, the gene MT-CO2 (P = 0.001, q = 0.09) in Complex IV (cytochrome c oxidase) and MT-ND2 (P = 0.004, q = 0.19) in Complex I (NADH dehydrogenase (ubiquinone)) were significantly associated with breast cancer risk.

Conclusions

In summary, our findings suggest that collective mitochondrial genetic variation and particularly in the MT-CO2 and MT-ND2 may play a role in breast cancer risk among European Americans. Further replication is warranted in larger populations and future studies should evaluate the contribution of mitochondrial proteins encoded by both the nuclear and mitochondrial genomes to breast cancer risk.

Quantitative mitochondrial DNA copy number determination using droplet digital PCR with single-cell resolution

Mitochondria are involved in a number of diverse cellular functions, including energy production, metabolic regulation, apoptosis, calcium homeostasis, cell proliferation, and motility, as well as free radical generation. Mitochondrial DNA (mtDNA) is present at hundreds to thousands of copies per cell in a tissue-specific manner. mtDNA copy number also varies during aging and disease progression and therefore might be considered as a biomarker that mirrors alterations within the human body. Here, we present a new quantitative, highly sensitive droplet digital PCR (ddPCR) method, droplet digital mitochondrial DNA measurement (ddMDM), to measure mtDNA copy number not only from cell populations but also from single cells. Our developed assay can generate data in as little as 3 h, is optimized for 96-well plates, and also allows the direct use of cell lysates without the need for DNA purification or nuclear reference genes. We show that ddMDM is able to detect differences between samples whose mtDNA copy number was close enough as to be indistinguishable by other commonly used mtDNA quantitation methods. By utilizing ddMDM, we show quantitative changes in mtDNA content per cell across a wide variety of physiological contexts including cancer progression, cell cycle progression, human T cell activation, and human aging.

Mitochondrial DNA, nuclear context, and the risk for carcinogenesis

The inheritance of mitochondrial DNA (mtDNA) from mother to child is complicated by differences in the stability of the mitochondrial genome. Although the germ line mtDNA is protected through the minimization of replication between generations, sequence variation can occur either through mutation or due to changes in the ratio between distinct genomes that are present in the mother (known as heteroplasmy). Thus, the unpredictability in transgenerational inheritance of mtDNA may cause the emergence of pathogenic mitochondrial and cellular phenotypes in offspring. Studies of the role of mitochondrial metabolism in cancer have a long and rich history, but recent evidence strongly suggests that changes in mitochondrial genotype and phenotype play a significant role in the initiation, progression and treatment of cancer. At the intersection of these two fields lies the potential for emerging mtDNA mutations to drive carcinogenesis in the offspring. In this review, we suggest that this facet of transgenerational carcinogenesis remains underexplored and is a potentially important contributor to cancer. Environ. Mol. Mutagen. 60:455-462, 2019. © 2018 Wiley Periodicals, Inc.

Nuclear genetic regulation of the human mitochondrial transcriptome

Mitochondria play important roles in cellular processes and disease, yet little is known about how the transcriptional regime of the mitochondrial genome varies across individuals and tissues. By analyzing >11,000 RNA-sequencing libraries across 36 tissue/cell types, we find considerable variation in mitochondrial-encoded gene expression along the mitochondrial transcriptome, across tissues and between individuals, highlighting the importance of cell-type specific and post-transcriptional processes in shaping mitochondrial-encoded RNA levels. Using whole-genome genetic data we identify 64 nuclear loci associated with expression levels of 14 genes encoded in the mitochondrial genome, including missense variants within genes involved in mitochondrial function (TBRG4, MTPAP and LONP1), implicating genetic mechanisms that act in trans across the two genomes. We replicate ~21% of associations with independent tissue-matched datasets and find genetic variants linked to these nuclear loci that are associated with cardio-metabolic phenotypes and Vitiligo, supporting a potential role for variable mitochondrial-encoded gene expression in complex disease.

Whole exome sequencing for the identification of CYP3A7 variants associated with tacrolimus concentrations in kidney transplant patients

The purpose of this study was to identify genotypes associated with dose-adjusted tacrolimus trough concentrations (C₀/D) in kidney transplant recipients using whole-exome sequencing (WES). This study included 147 patients administered tacrolimus, including seventy-five patients in the discovery set and seventy-two patients in the replication set. The patient genomes in the discovery set were sequenced using WES. Also, known tacrolimus pharmacokinetics-related intron variants were genotyped. Tacrolimus C₀/D was log-transformed. Sixteen variants were identified including novel CYP3A7 rs12360 and rs10211 by ANOVA. CYP3A7 rs2257401 was found to be the most significant variant among the periods by ANOVA. Seven variants including CYP3A7 rs2257401, rs12360, and rs10211 were analyzed by SNaPshot in the replication set and the effects on tacrolimus C₀/D were verified. A linear mixed model (LMM) was further performed to account for the effects of the variants and clinical factors. The combined set LMM showed that only CYP3A7 rs2257401 was associated with tacrolimus C₀/D after adjusting for patient age, albumin, and creatinine. The CYP3A7 rs2257401 genotype variant showed a significant difference on the tacrolimus C₀/D in those expressing CYP3A5, showing its own effect. The results suggest that CYP3A7 rs2257401 may serve as a significant genetic marker for tacrolimus pharmacokinetics in kidney transplantation.

Mitochondrial gene sequence variants in children with severe malaria anaemia with or without lactic acidosis: a case control study

BACKGROUND

Evolutionary pressure by Plasmodium falciparum malaria is known to have favoured a large number of human gene adaptations, but there is surprisingly little investigation of the effect of malaria on human mitochondrial sequence variation. Plasmodium falciparum infection can cause severe malaria anaemia (SMA) with insufficient tissue oxygenation, lactic acidosis and death. Despite equal degrees of severe anaemia, some individuals develop lactic acidosis while others do not. A case-control study design was used to investigate whether differences in host mitochondrial gene sequences were associated with lactic acidosis in SMA. Full mitochondrial sequences were obtained from 36 subjects with SMA complicated by lactic acidosis and 37 subjects with SMA without lactic acidosis. The two groups were matched for age, sex, and degree of anaemia.

RESULTS

Compared with the reference sequence, a median of 60 nucleotide variants per individual (interquartile range 4-91) was found, with an average frequency of 3.97 variants per 1000 nucleotides. The frequency and distribution of non-synonymous DNA variants in genes associated with oxidative phosphorylation were not statistically different between the two groups. Non-synonymous variants predicted to have the most disruptive effect on proteins responsible for oxidative phosphorylation were present at a similar frequency in both groups.

CONCLUSIONS

Lactic acidosis in SMA does not appear to be consistently associated with the high prevalence of any mitochondrial gene variant.

Whole sequence of the mitochondrial DNA genome of Kearns Sayre Syndrome patients: Identification of deletions and variants

Mitochondria both produce the energy of the cell as ATP via respiration and regulate cellular metabolism. Accordingly, any deletion or mutation in the mitochondrial DNA (mtDNA) may result in a disease. One of these diseases is Kearns Sayre syndrome (KSS), described for the first time in 1958, where different large-scale deletions of different sizes and at different positions have been reported in the mitochondrial genome of patients with similar clinical symptoms. In this study, sequences of the mitochondrial genome of three patients with clinic features of KSS were analyzed. Our results revealed the position, heteroplasmy percentage, size of deletions, and their haplogroups. Two patients contained deletions reported previously and one patient showed a new deletion not reported previously. These results display for the first time a systematic analysis of mtDNA variants in the whole mtDNA genome of patients with KSS to help to understand their association with the disease.

Meta-analysis identifies mitochondrial DNA sequence variants associated with walking speed

Declines in walking speed are associated with a variety of poor health outcomes including disability, comorbidity, and mortality. While genetic factors are putative contributors to variability in walking, few genetic loci have been identified for this trait. We examined the role of mitochondrial genomic variation on walking speed by sequencing the entire mitochondrial DNA (mtDNA). Data were meta-analyzed from 1758 Lifestyle Interventions and Independence for Elders (LIFE) Study and replication data from 730 Health, Aging, and Body Composition (HABC) Study participants with baseline walking speed information. Participants were 69+ years old of diverse racial backgrounds (African, European, and other race/ethnic groups) and had a wide range of mean walking speeds [4-6 m (0.78-1.09 m/s) and 400 m (0.83-1.24 m/s)]. Meta-analysis across studies and racial groups showed that m.12705C>T, ND5 variant was significantly associated (p < 0.0001) with walking speed at both short and long distances. Replication and meta-analysis also identified statistically significant walking speed associations (p < 0.0001) between the m.5460.G>A, ND2 and m.309C>CT, HV2 variants at short and long distances, respectively. All results remained statistically significant after multiple comparisons adjustment for 499 mtDNA variants. The m.12705C>T variant can be traced to the beginnings of human global migration and that cells carrying this variant display altered tRNA expression. Significant pooled effects related to stopping during the long-distance walk test were observed across OXPHOS complexes I (p = 0.0017) and III (p = 0.0048). These results suggest that mtDNA-encoded variants are associated with differences in walking speed among older adults, potentially identifying those at risk of developing mobility impairments.

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.

Divergent Patterns of Mitochondrial and Nuclear Ancestry Are Associated with the Risk for Preterm Birth

OBJECTIVE

To examine linkages between mitochondrial genetics and preterm birth by assessing the risk for preterm birth associated with the inheritance of nuclear haplotypes that are ancestrally distinct from mitochondrial haplogroup.

STUDY DESIGN

Genome-wide genotyping studies of cohorts of preterm and term individuals were evaluated. We determined the mitochondrial haplogroup and nuclear ancestry for individuals and developed a scoring for the degree to which mitochondrial ancestry is divergent from nuclear ancestry.

RESULTS

Infants with higher degrees of divergent mitochondrial ancestry were at increased risk for preterm birth (0.124 for preterm vs 0.105 for term infants; P< .05). This finding was validated in 1 of 2 replication cohorts. We also observed that greater degrees of divergent ancestry correlated with earlier delivery within the primary study population, but this finding was not replicated in secondary cohorts born preterm.

CONCLUSIONS

Individuals with divergent patterns of mitochondrial and nuclear ancestry are at increased risk for preterm birth. These findings may in part explain the higher rates of preterm birth in African Americans and in individuals with a matrilineal family history of preterm birth.

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.