Nuclear-embedded mitochondrial DNA sequences in 66,083 human genomes

Evolution and maintenance of mtDNA gene content across eukaryotes

Across eukaryotes, most genes required for mitochondrial function have been transferred to, or otherwise acquired by, the nucleus. Encoding genes in the nucleus has many advantages. So why do mitochondria retain any genes at all? Why does the set of mtDNA genes vary so much across different species? And how do species maintain functionality in the mtDNA genes they do retain? In this review, we will discuss some possible answers to these questions, attempting a broad perspective across eukaryotes. We hope to cover some interesting features which may be less familiar from the perspective of particular species, including the ubiquity of recombination outside bilaterian animals, encrypted chainmail-like mtDNA, single genes split over multiple mtDNA chromosomes, triparental inheritance, gene transfer by grafting, gain of mtDNA recombination factors, social networks of mitochondria, and the role of mtDNA dysfunction in feeding the world. We will discuss a unifying picture where organismal ecology and gene-specific features together influence whether organism X retains mtDNA gene Y, and where ecology and development together determine which strategies, importantly including recombination, are used to maintain the mtDNA genes that are retained.

Nanophotonic waveguide-based sensing of circulating cell-free mitochondrial DNA: implications for personalized medicine

Circulating cell-free mitochondrial DNA (ccf-mtDNA) has emerged as a promising biomarker, with potential implications for disease diagnosis. Changes in mtDNA, such as deletions, mutations or variations in the number of copies, have been associated with mitochondrial disorders, heart diseases, cancer and age-related non-communicable diseases. Previous methods, such as polymerase chain reaction-based approaches, next-generation sequencing and imaging-based techniques, have shown improved accuracy in identifying rare mtDNA variants or mutations, but they have limitations. This article explains the basic principles and benefits of using planar optical waveguide-based detection devices, which represent an advanced approach in the field of sensing.

Single-mitochondrion sequencing uncovers distinct mutational patterns and heteroplasmy landscape in mouse astrocytes and neurons

BACKGROUND

Mitochondrial (mt) heteroplasmy can cause adverse biological consequences when deleterious mtDNA mutations accumulate disrupting "normal" mt-driven processes and cellular functions. To investigate the heteroplasmy of such mtDNA changes, we developed a moderate throughput mt isolation procedure to quantify the mt single-nucleotide variant (SNV) landscape in individual mouse neurons and astrocytes. In this study, we amplified mt-genomes from 1645 single mitochondria isolated from mouse single astrocytes and neurons to (1) determine the distribution and proportion of mt-SNVs as well as mutation pattern in specific target regions across the mt-genome, (2) assess differences in mtDNA SNVs between neurons and astrocytes, and (3) study co-segregation of variants in the mouse mtDNA.

RESULTS

(1) The data show that specific sites of the mt-genome are permissive to SNV presentation while others appear to be under stringent purifying selection. Nested hierarchical analysis at the levels of mitochondrion, cell, and mouse reveals distinct patterns of inter- and intra-cellular variation for mt-SNVs at different sites. (2) Further, differences in the SNV incidence were observed between mouse neurons and astrocytes for two mt-SNV 9027:G > A and 9419:C > T showing variation in the mutational propensity between these cell types. Purifying selection was observed in neurons as shown by the Ka/Ks statistic, suggesting that neurons are under stronger evolutionary constraint as compared to astrocytes. (3) Intriguingly, these data show strong linkage between the SNV sites at nucleotide positions 9027 and 9461.

CONCLUSIONS

This study suggests that segregation as well as clonal expansion of mt-SNVs is specific to individual genomic loci, which is important foundational data in understanding of heteroplasmy and disease thresholds for mutation of pathogenic variants.

The rate and nature of mitochondrial DNA mutations in human pedigrees

We examined the rate and nature of mitochondrial DNA (mtDNA) mutations in humans using sequence data from 64,806 contemporary Icelanders from 2,548 matrilines. Based on 116,663 mother-child transmissions, 8,199 mutations were detected, providing robust rate estimates by nucleotide type, functional impact, position, and different alleles at the same position. We thoroughly document the true extent of hypermutability in mtDNA, mainly affecting the control region but also some coding-region variants. The results reveal the impact of negative selection on viable deleterious mutations, including rapidly mutating disease-associated 3243A>G and 1555A>G and pre-natal selection that most likely occurs during the development of oocytes. Finally, we show that the fate of new mutations is determined by a drastic germline bottleneck, amounting to an average of 3 mtDNA units effectively transmitted from mother to child.

The extent and characteristics of DNA transfer between plasmids and chromosomes

Plasmids are extrachromosomal genetic elements that reside in prokaryotes. The acquisition of plasmids encoding beneficial traits can facilitate short-term survival in harsh environmental conditions or long-term adaptation of new ecological niches. Due to their ability to transfer between cells, plasmids are considered agents of gene transfer. Nonetheless, the frequency of DNA transfer between plasmids and chromosomes remains understudied. Using a novel approach for detection of homologous loci between genome pairs, we uncover gene sharing with the chromosome in 1,974 (66%) plasmids residing in 1,016 (78%) taxonomically diverse isolates. The majority of homologous loci correspond to mobile elements, which may be duplicated in the host chromosomes in tens of copies. Neighboring shared genes often encode similar functional categories, indicating the transfer of multigene functional units. Rare transfer events of antibiotics resistance genes are observed mainly with mobile elements. The frequent erosion of sequence similarity in homologous regions indicates that the transferred DNA is often devoid of function. DNA transfer between plasmids and chromosomes thus generates genetic variation that is akin to workings of endosymbiotic gene transfer in eukaryotic evolution. Our findings imply that plasmid contribution to gene transfer most often corresponds to transfer of the plasmid entity rather than transfer of protein-coding genes between plasmids and chromosomes.

A novel protein CYTB-187AA encoded by the mitochondrial gene CYTB modulates mammalian early development

The mitochondrial genome transcribes 13 mRNAs coding for well-known proteins essential for oxidative phosphorylation. We demonstrate here that cytochrome b (CYTB), the only mitochondrial-DNA-encoded transcript among complex III, also encodes an unrecognized 187-amino-acid-long protein, CYTB-187AA, using the standard genetic code of cytosolic ribosomes rather than the mitochondrial genetic code. After validating the existence of this mtDNA-encoded protein arising from cytosolic translation (mPACT) using mass spectrometry and antibodies, we show that CYTB-187AA is mainly localized in the mitochondrial matrix and promotes the pluripotent state in primed-to-naive transition by interacting with solute carrier family 25 member 3 (SLC25A3) to modulate ATP production. We further generated a transgenic knockin mouse model of CYTB-187AA silencing and found that reduction of CYTB-187AA impairs females' fertility by decreasing the number of ovarian follicles. For the first time, we uncovered the novel mPACT pattern of a mitochondrial mRNA and demonstrated the physiological function of this 14th protein encoded by mtDNA.

Single-Mitochondrion Sequencing Uncovers Distinct Mutational Patterns and Heteroplasmy Landscape in Mouse Astrocytes and Neurons

Background

Mitochondrial (mt) heteroplasmy can cause adverse biological consequences when deleterious mtDNA mutations accumulate disrupting 'normal' mt-driven processes and cellular functions. To investigate the heteroplasmy of such mtDNA changes we developed a moderate throughput mt isolation procedure to quantify the mt single-nucleotide variant (SNV) landscape in individual mouse neurons and astrocytes In this study we amplified mt-genomes from 1,645 single mitochondria (mts) isolated from mouse single astrocytes and neurons to 1. determine the distribution and proportion of mt-SNVs as well as mutation pattern in specific target regions across the mt-genome, 2. assess differences in mtDNA SNVs between neurons and astrocytes, and 3. Study cosegregation of variants in the mouse mtDNA.

Results

1. The data show that specific sites of the mt-genome are permissive to SNV presentation while others appear to be under stringent purifying selection. Nested hierarchical analysis at the levels of mitochondrion, cell, and mouse reveals distinct patterns of inter- and intra-cellular variation for mt-SNVs at different sites. 2. Further, differences in the SNV incidence were observed between mouse neurons and astrocytes for two mt-SNV 9027:G>A and 9419:C>T showing variation in the mutational propensity between these cell types. Purifying selection was observed in neurons as shown by the Ka/Ks statistic, suggesting that neurons are under stronger evolutionary constraint as compared to astrocytes. 3. Intriguingly, these data show strong linkage between the SNV sites at nucleotide positions 9027 and 9461.

Conclusion

This study suggests that segregation as well as clonal expansion of mt-SNVs is specific to individual genomic loci, which is important foundational data in understanding of heteroplasmy and disease thresholds for mutation of pathogenic variants.

Nanobiopsy investigation of the subcellular mtDNA heteroplasmy in human tissues

Mitochondrial function is critical to continued cellular vitality and is an important contributor to a growing number of human diseases. Mitochondrial dysfunction is typically heterogeneous, mediated through the clonal expansion of mitochondrial DNA (mtDNA) variants in a subset of cells in a given tissue. To date, our understanding of the dynamics of clonal expansion of mtDNA variants has been technically limited to the single cell-level. Here, we report the use of nanobiopsy for subcellular sampling from human tissues, combined with next-generation sequencing to assess subcellular mtDNA mutation load in human tissue from mitochondrial disease patients. The ability to map mitochondrial mutation loads within individual cells of diseased tissue samples will further our understanding of mitochondrial genetic diseases.

THAP3 recruits SMYD3 to OXPHOS genes and epigenetically promotes mitochondrial respiration in hepatocellular carcinoma

Mitochondria harbor the oxidative phosphorylation (OXPHOS) system to sustain cellular respiration. However, the transcriptional regulation of OXPHOS remains largely unexplored. Through the cancer genome atlas (TCGA) transcriptome analysis, transcription factor THAP domain-containing 3 (THAP3) was found to be strongly associated with OXPHOS gene expression. Mechanistically, THAP3 recruited the histone methyltransferase SET and MYND domain-containing protein 3 (SMYD3) to upregulate H3K4me3 and promote OXPHOS gene expression. The levels of THAP3 and SMYD3 were altered by metabolic cues. They collaboratively supported liver cancer cell proliferation and colony formation. In clinical human liver cancer, both of them were overexpressed. THAP3 positively correlated with OXPHOS gene expression. Together, THAP3 cooperates with SMYD3 to epigenetically upregulate cellular respiration and liver cancer cell proliferation.

Connecting Species-Specific Extents of Genome Reduction in Mitochondria and Plastids

Mitochondria and plastids have both dramatically reduced their genomes since the endosymbiotic events that created them. The similarities and differences in the evolution of the two organelle genome types have been the target of discussion and investigation for decades. Ongoing work has suggested that similar mechanisms may modulate the reductive evolution of the two organelles in a given species, but quantitative data and statistical analyses exploring this picture remain limited outside of some specific cases like parasitism. Here, we use cross-eukaryote organelle genome data to explore evidence for coevolution of mitochondrial and plastid genome reduction. Controlling for differences between clades and pseudoreplication due to relatedness, we find that extents of mtDNA and ptDNA gene retention are related to each other across taxa, in a generally positive correlation that appears to differ quantitatively across eukaryotes, for example, between algal and nonalgal species. We find limited evidence for coevolution of specific mtDNA and ptDNA gene pairs, suggesting that the similarities between the two organelle types may be due mainly to independent responses to consistent evolutionary drivers.

HiFi long-read amplicon sequencing for full-spectrum variants of human mtDNA

BACKGROUND

Mitochondrial diseases (MDs) can be caused by single nucleotide variants (SNVs) and structural variants (SVs) in the mitochondrial genome (mtDNA). Presently, identifying deletions in small to medium-sized fragments and accurately detecting low-percentage variants remains challenging due to the limitations of next-generation sequencing (NGS).

METHODS

In this study, we integrated targeted long-range polymerase chain reaction (LR-PCR) and PacBio HiFi sequencing to analyze 34 participants, including 28 patients and 6 controls. Of these, 17 samples were subjected to both targeted LR-PCR and to compare the mtDNA variant detection efficacy.

RESULTS

Among the 28 patients tested by long-read sequencing (LRS), 2 patients were found positive for the m.3243 A > G hotspot variant, and 20 patients exhibited single or multiple deletion variants with a proportion exceeding 4%. Comparison between the results of LRS and NGS revealed that both methods exhibited similar efficacy in detecting SNVs exceeding 5%. However, LRS outperformed NGS in detecting SNVs with a ratio below 5%. As for SVs, LRS identified single or multiple deletions in 13 out of 17 cases, whereas NGS only detected single deletions in 8 cases. Furthermore, deletions identified by LRS were validated by Sanger sequencing and quantified in single muscle fibers using real-time PCR. Notably, LRS also effectively and accurately identified secondary mtDNA deletions in idiopathic inflammatory myopathies (IIMs).

CONCLUSIONS

LRS outperforms NGS in detecting various types of SNVs and SVs in mtDNA, including those with low frequencies. Our research is a significant advancement in medical comprehension and will provide profound insights into genetics.

Hypoxic adaptation of mitochondria and its impact on tumor cell function

Mitochondria are the major sink for oxygen in the cell, consuming it during ATP production. Therefore, when environmental oxygen levels drop in the tumor, significant adaptation is required. Mitochondrial activity is also a major producer of biosynthetic precursors and a regulator of cellular oxidative and reductive balance. Because of the complex biochemistry, mitochondrial adaptation to hypoxia occurs through multiple mechanisms and has significant impact on other cellular processes such as macromolecule synthesis and gene regulation. In tumor hypoxia, mitochondria shift their location in the cell and accelerate the fission and quality control pathways. Hypoxic mitochondria also undergo significant changes to fundamental metabolic pathways of carbon metabolism and electron transport. These metabolic changes further impact the nuclear epigenome because mitochondrial metabolites are used as enzymatic substrates for modifying chromatin. This coordinated response delivers physiological flexibility and increased tumor cell robustness during the environmental stress of low oxygen.

A PCR-independent approach for mtDNA enrichment and next-generation sequencing: comprehensive evaluation and clinical application

BACKGROUND

Sequencing the mitochondrial genome has been increasingly important for the investigation of primary mitochondrial diseases (PMD) and mitochondrial genetics. To overcome the limitations originating from PCR-based mtDNA enrichment, we set out to develop and evaluate a PCR-independent approach in this study, named Pime-Seq (PCR-independent mtDNA enrichment and next generation Sequencing).

RESULTS

By using the optimized mtDNA enrichment procedure, the mtDNA reads ratio reached 88.0 ± 7.9% in the sequencing library when applied on human PBMC samples. We found the variants called by Pime-Seq were highly consistent among technical repeats. To evaluate the accuracy and reliability of this method, we compared Pime-Seq with lrPCR based NGS by performing both methods simultaneously on 45 samples, yielding 1677 concordant variants, as well as 146 discordant variants with low-level heteroplasmic fraction, in which Pime-Seq showed higher reliability. Furthermore, we applied Pime-Seq on 4 samples of PMD patients retrospectively, and successfully detected all the pathogenic mtDNA variants. In addition, we performed a prospective study on 192 apparently healthy pregnant women during prenatal screening, in which Pime-Seq identified pathogenic mtDNA variants in 4 samples, providing extra information for better health monitoring in these cases.

CONCLUSIONS

Pime-Seq can obtain highly enriched mtDNA in a PCR-independent manner for high quality and reliable mtDNA deep-sequencing, which provides us an effective and promising tool for detecting mtDNA variants for both clinical and research purposes.

Long-read sequencing and structural variant characterization in 1,019 samples from the 1000 Genomes Project

Structural variants (SVs) contribute significantly to human genetic diversity and disease 1-4 . Previously, SVs have remained incompletely resolved by population genomics, with short-read sequencing facing limitations in capturing the whole spectrum of SVs at nucleotide resolution 5-7 . Here we leveraged nanopore sequencing 8 to construct an intermediate coverage resource of 1,019 long-read genomes sampled within 26 human populations from the 1000 Genomes Project. By integrating linear and graph-based approaches for SV analysis via pangenome graph-augmentation, we uncover 167,291 sequence-resolved SVs in these samples, considerably advancing SV characterization compared to population-wide short-read sequencing studies 3,4 . Our analysis details diverse SV classes-deletions, duplications, insertions, and inversions-at population-scale. LINE-1 and SVA retrotransposition activities frequently mediate transductions 9,10 of unique sequences, with both mobile element classes transducing sequences at either the 3'- or 5'-end, depending on the source element locus. Furthermore, analyses of SV breakpoint junctions suggest a continuum of homology-mediated rearrangement processes are integral to SV formation, and highlight evidence for SV recurrence involving repeat sequences. Our open-access dataset underscores the transformative impact of long-read sequencing in advancing the characterisation of polymorphic genomic architectures, and provides a resource for guiding variant prioritisation in future long-read sequencing-based disease studies.

Leveraging new methods for comprehensive characterization of mitochondrial DNA in esophageal squamous cell carcinoma

BACKGROUND

Mitochondria play essential roles in tumorigenesis; however, little is known about the contribution of mitochondrial DNA (mtDNA) to esophageal squamous cell carcinoma (ESCC). Whole-genome sequencing (WGS) is by far the most efficient technology to fully characterize the molecular features of mtDNA; however, due to the high redundancy and heterogeneity of mtDNA in regular WGS data, methods for mtDNA analysis are far from satisfactory.

METHODS

Here, we developed a likelihood-based method dMTLV to identify low-heteroplasmic mtDNA variants. In addition, we described fNUMT, which can simultaneously detect non-reference nuclear sequences of mitochondrial origin (non-ref NUMTs) and their derived artifacts. Using these new methods, we explored the contribution of mtDNA to ESCC utilizing the multi-omics data of 663 paired tumor-normal samples.

RESULTS

dMTLV outperformed the existing methods in sensitivity without sacrificing specificity. The verification using Nanopore long-read sequencing data showed that fNUMT has superior specificity and more accurate breakpoint identification than the current methods. Leveraging the new method, we identified a significant association between the ESCC overall survival and the ratio of mtDNA copy number of paired tumor-normal samples, which could be potentially explained by the differential expression of genes enriched in pathways related to metabolism, DNA damage repair, and cell cycle checkpoint. Additionally, we observed that the expression of CBWD1 was downregulated by the non-ref NUMTs inserted into its intron region, which might provide precursor conditions for the tumor cells to adapt to a hypoxic environment. Moreover, we identified a strong positive relationship between the number of mtDNA truncating mutations and the contribution of signatures linked to tumorigenesis and treatment response.

CONCLUSIONS

Our new frameworks promote the characterization of mtDNA features, which enables the elucidation of the landscapes and roles of mtDNA in ESCC essential for extending the current understanding of ESCC etiology. dMTLV and fNUMT are freely available from https://github.com/sunnyzxh/dMTLV and https://github.com/sunnyzxh/fNUMT , respectively.

Vitamin B12 status and folic acid supplementation influence mitochondrial heteroplasmy levels in mice

One-carbon metabolism is a complex network of metabolic reactions that are essential for cellular function including DNA synthesis. Vitamin B12 and folate are micronutrients that are utilized in this pathway and their deficiency can result in the perturbation of one-carbon metabolism and subsequent perturbations in DNA replication and repair. This effect has been well characterized in nuclear DNA but to date, mitochondrial DNA (mtDNA) has not been investigated extensively. Mitochondrial variants have been associated with several inherited and age-related disease states; therefore, the study of factors that impact heteroplasmy are important for advancing our understanding of the mitochondrial genome's impact on human health. Heteroplasmy studies require robust and efficient mitochondrial DNA enrichment to carry out in-depth mtDNA sequencing. Many of the current methods for mtDNA enrichment can introduce biases and false-positive results. Here, we use a method that overcomes these limitations and have applied it to assess mitochondrial heteroplasmy in mouse models of altered one-carbon metabolism. Vitamin B12 deficiency was found to cause increased levels of mitochondrial DNA heteroplasmy across all tissues that were investigated. Folic acid supplementation also contributed to elevated mitochondrial DNA heteroplasmy across all mouse tissues investigated. Heteroplasmy analysis of human data from the Framingham Heart Study suggested a potential sex-specific effect of folate and vitamin B12 status on mitochondrial heteroplasmy. This is a novel relationship that may have broader consequences for our understanding of one-carbon metabolism, mitochondrial-related disease and the influence of nutrients on DNA mutation rates.

Mitochondrial genome plasticity of mammalian species

There is an ongoing process in which mitochondrial sequences are being integrated into the nuclear genome. The importance of these sequences has already been revealed in cancer biology, forensic, phylogenetic studies and in the evolution of the eukaryotic genetic information. Human and numerous model organisms' genomes were described from those sequences point of view. Furthermore, recent studies were published on the patterns of these nuclear localised mitochondrial sequences in different taxa.However, the results of the previously released studies are difficult to compare due to the lack of standardised methods and/or using few numbers of genomes. Therefore, in this paper our primary goal is to establish a uniform mining pipeline to explore these nuclear localised mitochondrial sequences.Our results show that the frequency of several repetitive elements is higher in the flanking regions of these sequences than expected. A machine learning model reveals that the flanking regions' repetitive elements and different structural characteristics are highly influential during the integration process.In this paper, we introduce a general mining pipeline for all mammalian genomes. The workflow is publicly available and is believed to serve as a validated baseline for future research in this field. We confirm the widespread opinion, on - as to our current knowledge - the largest dataset, that structural circumstances and events corresponding to repetitive elements are highly significant. An accurate model has also been trained to predict these sequences and their corresponding flanking regions.

Mitochondrial DNA integrity and metabolome profile are preserved in the human induced pluripotent stem cell reference line KOLF2.1J

Quality control of human induced pluripotent stem cells (iPSCs) is critical to ensure reproducibility of research. Recently, KOLF2.1J was characterized and published as a male iPSC reference line to study neurological disorders. Emerging evidence suggests potential negative effects of mtDNA mutations, but its integrity was not analyzed in the original publication. To assess mtDNA integrity, we conducted a targeted mtDNA analysis followed by untargeted metabolomics analysis. We found that KOLF2.1J mtDNA integrity was intact at the time of publication and is still preserved in the commercially distributed cell line. In addition, the basal KOLF2.1J metabolome profile was similar to that of the two commercially available iPSC lines IMR90 and iPSC12, but clearly distinct from an in-house-generated ERCC6R683X/R683X iPSC line modeling Cockayne syndrome. Conclusively, we validate KOLF2.1J as a reference iPSC line, and encourage scientists to conduct mtDNA analysis and unbiased metabolomics whenever feasible.

RPA and Rad27 limit templated and inverted insertions at DNA breaks

Formation of templated insertions at DNA double-strand breaks (DSBs) is very common in cancer cells. The mechanisms and enzymes regulating these events are largely unknown. Here, we investigated templated insertions in yeast at DSBs using amplicon sequencing across a repaired locus. We document very short (most ∼5-34 bp), templated inverted duplications at DSBs. They are generated through a foldback mechanism that utilizes microhomologies adjacent to the DSB. Enzymatic requirements suggest a hybrid mechanism wherein one end requires Polδ-mediated synthesis while the other end is captured by nonhomologous end joining (NHEJ). This process is exacerbated in mutants with low levels or mutated RPA ( rtt105 Δ; rfa1 -t33) or extensive resection mutant ( sgs1 Δ exo1 Δ). Templated insertions from various distant genomic locations also increase in these mutants as well as in rad27 Δ and originate from fragile regions of the genome. Among complex insertions, common events are insertions of two sequences, originating from the same locus and with inverted orientation. We propose that these inversions are also formed by microhomology-mediated template switching. Taken together, we propose that a shortage of RPA typical in cancer cells is one possible factor stimulating the formation of templated insertions.

A complete workflow for single cell mtDNAseq in CHO cells, from cell culture to bioinformatic analysis

Chinese hamster ovary (CHO) cells have a long history in the biopharmaceutical industry and currently produce the vast majority of recombinant therapeutic proteins. A key step in controlling the process and product consistency is the development of a producer cell line derived from a single cell clone. However, it is recognized that genetic and phenotypic heterogeneity between individual cells in a clonal CHO population tends to arise over time. Previous bulk analysis of CHO cell populations revealed considerable variation within the mtDNA sequence (heteroplasmy), which could have implications for the performance of the cell line. By analyzing the heteroplasmy of single cells within the same population, this heterogeneity can be characterized with greater resolution. Such analysis may identify heterogeneity in the mitochondrial genome, which impacts the overall phenotypic performance of a producer cell population, and potentially reveal routes for genetic engineering. A critical first step is the development of robust experimental and computational methods to enable single cell mtDNA sequencing (termed scmtDNAseq). Here, we present a protocol from cell culture to bioinformatic analysis and provide preliminary evidence of significant mtDNA heteroplasmy across a small panel of single CHO cells.

Screening for Rare Mitochondrial Genome Variants Reveals a Potentially Novel Association between MT-CO1 and MT-TL2 Genes and Diabetes Phenotype

Variations in several nuclear genes predisposing humans to the development of MODY diabetes have been very well characterized by modern genetic diagnostics. However, recent reports indicate that variants in the mtDNA genome may also be associated with the diabetic phenotype. As relatively little research has addressed the entire mitochondrial genome in this regard, the aim of the present study is to evaluate the genetic variations present in mtDNA among individuals susceptible to MODY diabetes. In total, 193 patients with a MODY phenotype were tested with a custom panel with mtDNA enrichment. Heteroplasmic variants were selected for further analysis via further sequencing based on long-range PCR to evaluate the potential contribution of frequent NUMTs (acronym for nuclear mitochondrial DNA) insertions. Twelve extremely rare variants with a potential damaging character were selected, three of which were likely to be the result of NUMTs from the nuclear genome. The variant m.3243A>G in MT-TL1 was responsible for 3.5% of MODY cases in our study group. In addition, a novel, rare, and possibly pathogenic leucine variant m.12278T>C was found in MT-TL2. Our findings also found the MT-CO1 gene to be over-represented in the study group, with a clear phenotype-genotype correlation observed in one family. Our data suggest that heteroplasmic variants in MT-COI and MT-TL2 genes may play a role in the pathophysiology of glucose metabolism in humans.

The known unknowns of mitochondrial carcinogenesis: de novo NUMTs and intercellular mitochondrial transfer

The translocation of mitochondrial DNA (mtDNA) sequences into the nuclear genome, resulted in the occurrence of nuclear sequences of mitochondrial origin (NUMTs) which can be detected in nearly all sequenced eukaryotes. However, de novo mtDNA insertions can contribute to the development of pathological conditions including cancer. Recent data indicate that de novo mtDNA translocation into chromosomes can occur due to genotoxic influence of DNA double-strand break-inducing environmental mutagens. This confirms the hypothesis of the involvement of genome instability in the occurrence of mtDNA fragments in chromosomes. Mounting evidence indicates that mitochondria can be transferred from normal cells to cancer cells and recover cellular respiration. These exchanged mitochondria can facilitate cancer progression and metastasis. This review article provides a comprehensive overview of the potential carcinogenicity of mtDNA insertions, and the relevance of mtDNA escape in cancer progression, metastasis, and treatment resistance in humans. Potential molecular targets involved in mtDNA escape and exchange of mitochondria that can be of possible clinical benefits are presented and discussed. Understanding these processes could lead to improved diagnostic approaches, novel therapeutic strategies, and a deeper understanding of the intricate relationship between mitochondria, nuclear DNA, and cancer biology.

Mitochondrial genomes revisited: why do different lineages retain different genes?

The mitochondria contain their own genome derived from an alphaproteobacterial endosymbiont. From thousands of protein-coding genes originally encoded by their ancestor, only between 1 and about 70 are encoded on extant mitochondrial genomes (mitogenomes). Thanks to a dramatically increasing number of sequenced and annotated mitogenomes a coherent picture of why some genes were lost, or relocated to the nucleus, is emerging. In this review, we describe the characteristics of mitochondria-to-nucleus gene transfer and the resulting varied content of mitogenomes across eukaryotes. We introduce a 'burst-upon-drift' model to best explain nuclear-mitochondrial population genetics with flares of transfer due to genetic drift.

Yeast EndoG prevents genome instability by degrading cytoplasmic DNA

In metazoans release of mitochondrial DNA or retrotransposon cDNA to cytoplasm can cause sterile inflammation and disease 1. Cytoplasmic nucleases degrade these DNA species to limit inflammation 2,3. It remains unknown whether degradation these DNA also prevents nuclear genome instability. To address this question, we decided to identify the nuclease regulating transfer of these cytoplasmic DNA species to the nucleus. We used an amplicon sequencing-based method in yeast enabling analysis of millions of DSB repair products. Nuclear mtDNA (NUMTs) and retrotransposon cDNA insertions increase dramatically in nondividing stationary phase cells. Yeast EndoG (Nuc1) nuclease limits insertions of cDNA and transfer of very long mtDNA (>10 kb) that forms unstable circles or rarely insert in the genome, but it promotes formation of short NUMTs (~45-200 bp). Nuc1 also regulates transfer of cytoplasmic DNA to nucleus in aging or during meiosis. We propose that Nuc1 preserves genome stability by degrading retrotransposon cDNA and long mtDNA, while short NUMTs can originate from incompletely degraded mtDNA. This work suggests that nucleases eliminating cytoplasmic DNA play a role in preserving genome stability.

Comparative Genome Microsynteny Illuminates the Fast Evolution of Nuclear Mitochondrial Segments (NUMTs) in Mammals

The escape of DNA from mitochondria into the nuclear genome (nuclear mitochondrial DNA, NUMT) is an ongoing process. Although pervasively observed in eukaryotic genomes, their evolutionary trajectories in a mammal-wide context are poorly understood. The main challenge lies in the orthology assignment of NUMTs across species due to their fast evolution and chromosomal rearrangements over the past 200 million years. To address this issue, we systematically investigated the characteristics of NUMT insertions in 45 mammalian genomes and established a novel, synteny-based method to accurately predict orthologous NUMTs and ascertain their evolution across mammals. With a series of comparative analyses across taxa, we revealed that NUMTs may originate from nonrandom regions in mtDNA, are likely found in transposon-rich and intergenic regions, and unlikely code for functional proteins. Using our synteny-based approach, we leveraged 630 pairwise comparisons of genome-wide microsynteny and predicted the NUMT orthology relationships across 36 mammals. With the phylogenetic patterns of NUMT presence-and-absence across taxa, we constructed the ancestral state of NUMTs given the mammal tree using a coalescent method. We found support on the ancestral node of Fereuungulata within Laurasiatheria, whose subordinal relationships are still controversial. This study broadens our knowledge on NUMT insertion and evolution in mammalian genomes and highlights the merit of NUMTs as alternative genetic markers in phylogenetic inference.

Numerous insertions of mitochondrial DNA in the genome of the northern mole vole, Ellobius talpinus

BACKGROUND

Ellobius talpinus is a subterranean rodent representing an attractive model in population ecology studies due to its highly special lifestyle and sociality. In such studies, mitochondrial DNA (mtDNA) is widely used. However, if nuclear copies of mtDNA, aka NUMTs, are present, they may co-amplify with the target mtDNA fragment, generating misleading results. The aim of this study was to determine whether NUMTs are present in E. talpinus.

METHODS AND RESULTS

PCR amplification of the putative mtDNA CytB-D-loop fragment using 'universal' primers from 56 E. talpinus samples produced multiple double peaks in 90% of the sequencing chromatograms. To reveal NUMTs, molecular cloning and sequencing of PCR products of three specimens was conducted, followed by phylogenetic analysis. The pseudogene nature of three out of the seven detected haplotypes was confirmed by their basal positions in relation to other Ellobius haplotypes in the phylogenetic tree. Additionally, 'haplotype B' was basal in relation to other E. talpinus haplotypes and found present in very distant sampling sites. BLASTN search revealed 195 NUMTs in the E. talpinus nuclear genome, including fragments of all four PCR amplified pseudogenes. Although the majority of the NUMTs studied were short, the entire mtDNA had copies in the nuclear genome. The most numerous NUMTs were found for rrnL, COXI, and D-loop.

CONCLUSIONS

Numerous NUMTs are present in E. talpinus and can be difficult to discriminate against mtDNA sequences. Thus, in future population or phylogenetic studies in E. talpinus, the possibility of cryptic NUMTs amplification should always be taken into account.

Yeast EndoG prevents genome instability by degrading cytoplasmic DNA

In metazoans release of mitochondrial DNA or retrotransposon cDNA to cytoplasm can cause sterile inflammation and disease. Cytoplasmic nucleases degrade these DNA species to limit inflammation. It remains unknown whether degradation these DNA also prevents nuclear genome instability. To address this question, we decided to identify the nuclease regulating transfer of these cytoplasmic DNA species to the nucleus. We used an amplicon sequencing-based method in yeast enabling analysis of millions of DSB repair products. Nu clear mt DNA (NUMTs) and retrotransposon cDNA insertions increase dramatically in nondividing stationary phase cells. Yeast EndoG (Nuc1) nuclease limits insertions of cDNA and transfer of very long mtDNA (>10 kb) that forms unstable circles or rarely insert in the genome, but it promotes formation of short NUMTs (∼45-200 bp). Nuc1 also regulates transfer of cytoplasmic DNA to nucleus in aging or during meiosis. We propose that Nuc1 preserves genome stability by degrading retrotransposon cDNA and long mtDNA, while short NUMTs can originate from incompletely degraded mtDNA. This work suggests that nucleases eliminating cytoplasmic DNA play a role in preserving genome stability.

North and East African mitochondrial genetic variation needs further characterization towards precision medicine

INTRODUCTION

Mitochondria are maternally inherited cell organelles with their own genome, and perform various functions in eukaryotic cells such as energy production and cellular homeostasis. Due to their inheritance and manifold biological roles in health and disease, mitochondrial genetics serves a dual purpose of tracing the history as well as disease susceptibility of human populations across the globe. This work requires a comprehensive catalogue of commonly observed genetic variations in the mitochondrial DNAs for all regions throughout the world. So far, however, certain regions, such as North and East Africa have been understudied.

OBJECTIVES

To address this shortcoming, we have created the most comprehensive quality-controlled North and East African mitochondrial data set to date and use it for characterizing mitochondrial genetic variation in this region.

METHODS

We compiled 11 published cohorts with novel data for mitochondrial genomes from 159 Sudanese individuals. We combined these 641 mitochondrial sequences with sequences from the 1000 Genomes (n = 2504) and the Human Genome Diversity Project (n = 828) and used the tool haplocheck for extensive quality control and detection of in-sample contamination, as well as Nanopore long read sequencing for haplogroup validation of 18 samples.

RESULTS

Using a subset of high-coverage mitochondrial sequences, we predict 15 potentially novel haplogroups in North and East African subjects and observe likely phylogenetic deviations from the established PhyloTree reference for haplogroups L0a1 and L2a1.

CONCLUSION

Our findings demonstrate common hitherto unexplored variants in mitochondrial genomes of North and East Africa that lead to novel phylogenetic relationships between haplogroups present in these regions. These observations call for further in-depth population genetic studies in that region to enable the prospective use of mitochondrial genetic variation for precision medicine.

Comprehensive Identification of Mitochondrial Pseudogenes (NUMTs) in the Human Telomere-to-Telomere Reference Genome

Practices related to mitochondrial research have long been hindered by the presence of mitochondrial pseudogenes within the nuclear genome (NUMTs). Even though partially assembled human reference genomes like hg38 have included NUMTs compilation, the exhaustive NUMTs within the only complete reference genome (T2T-CHR13) remain unknown. Here, we comprehensively identified the fixed NUMTs within the reference genome using human pan-mitogenome (HPMT) from GeneBank. The inclusion of HPMT serves the purpose of establishing an authentic mitochondrial DNA (mtDNA) mutational spectrum for the identification of NUMTs, distinguishing it from the polymorphic variations found in NUMTs. Using HPMT, we identified approximately 10% of additional NUMTs in three human reference genomes under stricter thresholds. And we also observed an approximate 6% increase in NUMTs in T2T-CHR13 compared to hg38, including NUMTs on the short arms of chromosomes 13, 14, and 15 that were not assembled previously. Furthermore, alignments based on 20-mer from mtDNA suggested the presence of more mtDNA-like short segments within the nuclear genome, which should be avoided for short amplicon or cell free mtDNA detection. Finally, through the assay of transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) on cell lines before and after mtDNA elimination, we concluded that NUMTs have a minimal impact on bulk ATAC-seq, even though 16% of sequencing data originated from mtDNA.

Massive invasion of organellar DNA drives nuclear genome evolution in Toxoplasma

Toxoplasma gondii is a zoonotic protist pathogen that infects up to one third of the human population. This apicomplexan parasite contains three genome sequences: nuclear (65 Mb); plastid organellar, ptDNA (35 kb); and mitochondrial organellar, mtDNA (5.9 kb of non-repetitive sequence). We find that the nuclear genome contains a significant amount of NUMTs (nuclear integrants of mitochondrial DNA) and NUPTs (nuclear integrants of plastid DNA) that are continuously acquired and represent a significant source of intraspecific genetic variation. NUOT (nuclear DNA of organellar origin) accretion has generated 1.6% of the extant T. gondii ME49 nuclear genome-the highest fraction ever reported in any organism. NUOTs are primarily found in organisms that retain the non-homologous end-joining repair pathway. Significant movement of organellar DNA was experimentally captured via amplicon sequencing of a CRISPR-induced double-strand break in non-homologous end-joining repair competent, but not ku80 mutant, Toxoplasma parasites. Comparisons with Neospora caninum, a species that diverged from Toxoplasma ~28 mya, revealed that the movement and fixation of five NUMTs predates the split of the two genera. This unexpected level of NUMT conservation suggests evolutionary constraint for cellular function. Most NUMT insertions reside within (60%) or nearby genes (23% within 1.5 kb), and reporter assays indicate that some NUMTs have the ability to function as cis-regulatory elements modulating gene expression. Together, these findings portray a role for organellar sequence insertion in dynamically shaping the genomic architecture and likely contributing to adaptation and phenotypic changes in this important human pathogen.

The constraints of allotopic expression

Allotopic expression is the functional transfer of an organellar gene to the nucleus, followed by synthesis of the gene product in the cytosol and import into the appropriate organellar sub compartment. Here, we focus on mitochondrial genes encoding OXPHOS subunits that were naturally transferred to the nucleus, and critically review experimental evidence that claim their allotopic expression. We emphasize aspects that may have been overlooked before, i.e., when modifying a mitochondrial gene for allotopic expression━besides adapting the codon usage and including sequences encoding mitochondrial targeting signals━three additional constraints should be considered: (i) the average apparent free energy of membrane insertion (μΔGapp) of the transmembrane stretches (TMS) in proteins earmarked for the inner mitochondrial membrane, (ii) the final, functional topology attained by each membrane-bound OXPHOS subunit; and (iii) the defined mechanism by which the protein translocator TIM23 sorts cytosol-synthesized precursors. The mechanistic constraints imposed by TIM23 dictate the operation of two pathways through which alpha-helices in TMS are sorted, that eventually determine the final topology of membrane proteins. We used the biological hydrophobicity scale to assign an average apparent free energy of membrane insertion (μΔGapp) and a "traffic light" color code to all TMS of OXPHOS membrane proteins, thereby predicting which are more likely to be internalized into mitochondria if allotopically produced. We propose that the design of proteins for allotopic expression must make allowance for μΔGapp maximization of highly hydrophobic TMS in polypeptides whose corresponding genes have not been transferred to the nucleus in some organisms.

Complete mitochondrial DNA profile in stroke: A geographical matched case-control study in Spanish population

INTRODUCTION

Stroke, the second leading cause of death worldwide, is a complex disease influenced by many risk factors among which we can find reactive oxygen species (ROS). Since mitochondria are the main producers of cellular ROS, nowadays studies are trying to elucidate the role of these organelles and its DNA (mtDNA) variation in stroke risk. The aim of the present study was to perform a comprehensive evaluation of the association between mtDNA mutations and mtDNA content and stroke risk.

MATERIAL AND METHODS

Homoplasmic and heteroplasmic mutations of the mtDNA were analysed in a case-controls study using 110 S cases and their corresponding control individuals. Mitochondrial DNA copy number (mtDNA-CN) was analysed in 73 of those case-control pairs.

RESULTS

Our results suggest that haplogroup V, specifically variants m.72C > T, m.4580G > A, m.15904C > T and m.16298 T > C have a protective role in relation to stroke risk. On the contrary, variants m.73A > G, m.11719G > A and m.14766C > T appear to be genetic risk factors for stroke. In this study, we found no statistically significant association between stroke risk and mitochondrial DNA copy number.

CONCLUSIONS

These results demonstrate the possible role of mtDNA genetics on the pathogenesis of stroke, probably through alterations in mitochondrial ROS production.

A human identification system for hair shaft using RNA polymorphism

Hair is one of the common pieces of evidence at crime scenes, with abundant mitochondrial DNA but limited nuclear DNA in its shaft. It also helps to narrow the investigation scope to maternal lineage but fails to provide unique individual information. We assumed that RNA in hair shafts would be an alternative resource used to perform human identification based on the facts that (1) RNA retains the polymorphic information; (2) the multi-copy of RNA in a cell resists degradation as compared to the one-copy of nuclear DNA. In this study, we explored the potential of RNA polymorphism in hair shafts for forensic individual identification. A SNaPshot typing system was constructed using 18 SNPs located on 11 genes (ABCA13, AHNAK, EXPH5, KMT2D, KRT35, PPP1R15A, RBM33, S100A5, TBC1D4, TMC5, TRPV2). The RNA typing system was evaluated for sensitivity, species specificity, and feasibility for aged hair samples. Hair samples from a Shanxi population in China were used for the population study of the system. The detection limit of the assay was 2 ng RNA. The CDP of these 11 genes was 0.999969 in the Shanxi population. We also identified the concordance of the RNA and DNA typing results. In summary, we developed an RNA typing method to perform human identification from hair shafts, which performed as accurately as nuclear DNA typing. Our method provides a potential basis for solving the human identification problem from hair shafts, as well as other biological materials that lack nuclear DNA.

Target-Induced Stepwise Disintegration of Starlike Branched and Multiplex Embedded Systems for Amplified Detection of Serum MicroRNA

DNA nanotechnology has shown great promise for biosensing and molecular recognition. However, the practical application of conventional DNA biosensors is constrained by inadequate target stimuli, intricate design schemes, multicomponent systems, and susceptibility to nuclease degradation. To overcome these limitations, we present a class of starlike branched and multiplex embedded system (SBES) with an integrated functional design and cascade exponential amplification for serum microRNA (miRNA) detection. The DNA arms can be integrated into an all-in-one system by surrounding a branch point, with each arm endowed with specific functionalities by embedding different DNA fragments. These fragments include a segment complementary to the target miRNA for the recognition element, palindromic tails for self-primed polymerization, and a region with the same sequences as the target serving as the target analogue. Upon exposure to a target miRNA, the DNA arms unwind in a stepwise manner through palindrome-mediated dimerization and polymerization. This enables target recycling for subsequent reactions while releasing the target analogue to generate a secondary response in a feedback manner. A comparative analysis illustrates that the signal-to-noise ratio (SNR) of a full SBES with a feedback strategy is approximately 250% higher than the system without a feedback design. We demonstrate that the four-arm 4pSBES has the benefits of multifunctional integration, enhanced sensitivity, and low false-positive signals, which makes this approach ideally suited for clinical diagnosis. Moreover, an upgraded SBES with additional DNA arms (e.g., 6pSBES) can be constructed to allow multifunctional extension, offering unprecedented opportunities to build versatile DNA nanostructures for biosensing.

Dysregulation of circular RNAs in inflammation and cancers

Emerging lines of evidence have shown that the production of the covalently closed single-stranded circular RNAs is not splicing errors, but rather a regulated process with distinct biogenesis and turnover. Circular RNAs are expressed in a cell type- and tissue-specific manner and often localize to specific subcellular regions or organelles for functions. The dysregulation of circular RNAs from birth to death is linked to the pathogenesis and progression of diverse diseases. This review outlines how aberrant circular RNA biogenesis, subcellular location, and degradation are linked to disease progression, focusing on metaflammation and cancers. We also discuss potential therapeutic strategies and obstacles in targeting such disease-related circular RNAs.

Interplay between endogenous and exogenous human retroviruses

In humans, retroviruses thrive more as symbionts than as parasites. Apart from the only two modern exogenous human retroviruses (human T-cell lymphotropic and immunodeficiency viruses; HTLV and HIV, respectively), ~8% of the human genome is occupied by ancient retroviral DNA [human endogenous retroviruses (HERVs)]. Here, we review the recent discoveries about the interactions between the two groups, the impact of infection by exogenous retroviruses on the expression of HERVs, the effect of HERVs on the pathogenicity of HIV and HTLV and on the severity of the diseases caused by them, and the antiviral protection that HERVs can allegedly provide to the host. Tracing the crosstalk between contemporary retroviruses and their endogenized ancestors will provide better understanding of the retroviral world.

Mitochondrial polymorphism m.3017C>T of SHLP6 relates to heterothermy

Heterothermic thermoregulation requires intricate regulation of metabolic rate and activation of pro-survival factors. Eliciting these responses and coordinating the necessary energy shifts likely involves retrograde signalling by mitochondrial-derived peptides (MDPs). Members of the group were suggested before to play a role in heterothermic physiology, a key component of hibernation and daily torpor. Here we studied the mitochondrial single-nucleotide polymorphism (SNP) m.3017C>T that resides in the evolutionarily conserved gene MT-SHLP6. The substitution occurring in several mammalian orders causes truncation of SHLP6 peptide size from twenty to nine amino acids. Public mass spectrometric (MS) data of human SHLP6 indicated a canonical size of 20 amino acids, but not the use of alternative translation initiation codons that would expand the peptide. The shorter isoform of SHLP6 was found in heterothermic rodents at higher frequency compared to homeothermic rodents (p < 0.001). In heterothermic mammals it was associated with lower minimal body temperature (T b, p < 0.001). In the thirteen-lined ground squirrel, brown adipose tissue-a key organ required for hibernation, showed dynamic changes of the steady-state transcript level of mt-Shlp6. The level was significantly higher before hibernation and during interbout arousal and lower during torpor and after hibernation. Our finding argues to further explore the mode of action of SHLP6 size isoforms with respect to mammalian thermoregulation and possibly mitochondrial retrograde signalling.

Optimized bisulfite sequencing analysis reveals the lack of 5-methylcytosine in mammalian mitochondrial DNA

BACKGROUND

DNA methylation is one of the best characterized epigenetic modifications in the mammalian nuclear genome and is known to play a significant role in various biological processes. Nonetheless, the presence of 5-methylcytosine (5mC) in mitochondrial DNA remains controversial, as data ranging from the lack of 5mC to very extensive 5mC have been reported.

RESULTS

By conducting comprehensive bioinformatic analyses of both published and our own data, we reveal that previous observations of extensive and strand-biased mtDNA-5mC are likely artifacts due to a combination of factors including inefficient bisulfite conversion, extremely low sequencing reads in the L strand, and interference from nuclear mitochondrial DNA sequences (NUMTs). To reduce false positive mtDNA-5mC signals, we establish an optimized procedure for library preparation and data analysis of bisulfite sequencing. Leveraging our modified workflow, we demonstrate an even distribution of 5mC signals across the mtDNA and an average methylation level ranging from 0.19% to 0.67% in both cell lines and primary cells, which is indistinguishable from the background noise.

CONCLUSIONS

We have developed a framework for analyzing mtDNA-5mC through bisulfite sequencing, which enables us to present multiple lines of evidence for the lack of extensive 5mC in mammalian mtDNA. We assert that the data available to date do not support the reported presence of mtDNA-5mC.

Mitochondrial DNA-targeted therapy: A novel approach to combat cancer

Mitochondrial DNA (mtDNA) encodes proteins and RNAs that are essential for mitochondrial function and cellular homeostasis, and participates in important processes of cellular bioenergetics and metabolism. Alterations in mtDNA are associated with various diseases, especially cancers, and are considered as biomarkers for some types of tumors. Moreover, mtDNA alterations have been found to affect the proliferation, progression and metastasis of cancer cells, as well as their interactions with the immune system and the tumor microenvironment (TME). The important role of mtDNA in cancer development makes it a significant target for cancer treatment. In recent years, many novel therapeutic methods targeting mtDNA have emerged. In this study, we first discussed how cancerogenesis is triggered by mtDNA mutations, including alterations in gene copy number, aberrant gene expression and epigenetic modifications. Then, we described in detail the mechanisms underlying the interactions between mtDNA and the extramitochondrial environment, which are crucial for understanding the efficacy and safety of mtDNA-targeted therapy. Next, we provided a comprehensive overview of the recent progress in cancer therapy strategies that target mtDNA. We classified them into two categories based on their mechanisms of action: indirect and direct targeting strategies. Indirect targeting strategies aimed to induce mtDNA damage and dysfunction by modulating pathways that are involved in mtDNA stability and integrity, while direct targeting strategies utilized molecules that can selectively bind to or cleave mtDNA to achieve the therapeutic efficacy. This study highlights the importance of mtDNA-targeted therapy in cancer treatment, and will provide insights for future research and development of targeted drugs and therapeutic strategies.

Nuclear genetic control of mtDNA copy number and heteroplasmy in humans

Mitochondrial DNA (mtDNA) is a maternally inherited, high-copy-number genome required for oxidative phosphorylation1. Heteroplasmy refers to the presence of a mixture of mtDNA alleles in an individual and has been associated with disease and ageing. Mechanisms underlying common variation in human heteroplasmy, and the influence of the nuclear genome on this variation, remain insufficiently explored. Here we quantify mtDNA copy number (mtCN) and heteroplasmy using blood-derived whole-genome sequences from 274,832 individuals and perform genome-wide association studies to identify associated nuclear loci. Following blood cell composition correction, we find that mtCN declines linearly with age and is associated with variants at 92 nuclear loci. We observe that nearly everyone harbours heteroplasmic mtDNA variants obeying two principles: (1) heteroplasmic single nucleotide variants tend to arise somatically and accumulate sharply after the age of 70 years, whereas (2) heteroplasmic indels are maternally inherited as mixtures with relative levels associated with 42 nuclear loci involved in mtDNA replication, maintenance and novel pathways. These loci may act by conferring a replicative advantage to certain mtDNA alleles. As an illustrative example, we identify a length variant carried by more than 50% of humans at position chrM:302 within a G-quadruplex previously proposed to mediate mtDNA transcription/replication switching2,3. We find that this variant exerts cis-acting genetic control over mtDNA abundance and is itself associated in-trans with nuclear loci encoding machinery for this regulatory switch. Our study suggests that common variation in the nuclear genome can shape variation in mtCN and heteroplasmy dynamics across the human population.

De novo assembled mitogenome analysis of Trichuris trichiura from Korean individuals using nanopore-based long-read sequencing technology

Knowledge about mitogenomes has been proven to be essential in human parasite diagnostics and understanding of their diversity. However, the lack of substantial data for comparative analysis is still a challenge in Trichuris trichiura research. To provide high quality mitogenomes, we utilized long-read sequencing technology of Oxford Nanopore Technologies (ONT) to better resolve repetitive regions and to construct de novo mitogenome assembly minimizing reference biases. In this study, we got three de novo assembled mitogenomes of T. trichiura isolated from Korean individuals. These circular complete mitogenomes of T. trichiura are 14,508 bp, 14,441 bp, and 14,440 bp in length. A total of 37 predicted genes were identified consisting of 13 protein-coding genes (PCGs), 22 transfer RNA (tRNAs) genes, two ribosomal RNA (rRNA) genes (rrnS and rrnL), and two non-coding regions. Interestingly, the assembled mitogenome has up to six times longer AT-rich regions than previous reference sequences, thus proving the advantage of long-read sequencing in resolving unreported non-coding regions. Furthermore, variant detection and phylogenetic analysis using concatenated protein coding genes, cox1, rrnL, and nd1 genes confirmed the distinct molecular identity of this newly assembled mitogenome while at the same time showing high genetic relationship with sequences from China or Tanzania. Our study provided a new set of reference mitogenome with better contiguity and resolved repetitive regions that could be used for meaningful phylogenetic analysis to further understand disease transmission and parasite biology.

Mitochondrial DNA in Human Diversity and Health: From the Golden Age to the Omics Era

Mitochondrial DNA (mtDNA) is a small fraction of our hereditary material. However, this molecule has had an overwhelming presence in scientific research for decades until the arrival of high-throughput studies. Several appealing properties justify the application of mtDNA to understand how human populations are-from a genetic perspective-and how individuals exhibit phenotypes of biomedical importance. Here, I review the basics of mitochondrial studies with a focus on the dawn of the field, analysis methods and the connection between two sides of mitochondrial genetics: anthropological and biomedical. The particularities of mtDNA, with respect to inheritance pattern, evolutionary rate and dependence on the nuclear genome, explain the challenges of associating mtDNA composition and diseases. Finally, I consider the relevance of this single locus in the context of omics research. The present work may serve as a tribute to a tool that has provided important insights into the past and present of humankind.

The Vital Role Played by Deferiprone in the Transition of Thalassaemia from a Fatal to a Chronic Disease and Challenges in Its Repurposing for Use in Non-Iron-Loaded Diseases

The iron chelating orphan drug deferiprone (L1), discovered over 40 years ago, has been used daily by patients across the world at high doses (75-100 mg/kg) for more than 30 years with no serious toxicity. The level of safety and the simple, inexpensive synthesis are some of the many unique properties of L1, which played a major role in the contribution of the drug in the transition of thalassaemia from a fatal to a chronic disease. Other unique and valuable clinical properties of L1 in relation to pharmacology and metabolism include: oral effectiveness, which improved compliance compared to the prototype therapy with subcutaneous deferoxamine; highly effective iron removal from all iron-loaded organs, particularly the heart, which is the major target organ of iron toxicity and the cause of mortality in thalassaemic patients; an ability to achieve negative iron balance, completely remove all excess iron, and maintain normal iron stores in thalassaemic patients; rapid absorption from the stomach and rapid clearance from the body, allowing a greater frequency of repeated administration and overall increased efficacy of iron excretion, which is dependent on the dose used and also the concentration achieved at the site of drug action; and its ability to cross the blood-brain barrier and treat malignant, neurological, and microbial diseases affecting the brain. Some differential pharmacological activity by L1 among patients has been generally shown in relation to the absorption, distribution, metabolism, elimination, and toxicity (ADMET) of the drug. Unique properties exhibited by L1 in comparison to other drugs include specific protein interactions and antioxidant effects, such as iron removal from transferrin and lactoferrin; inhibition of iron and copper catalytic production of free radicals, ferroptosis, and cuproptosis; and inhibition of iron-containing proteins associated with different pathological conditions. The unique properties of L1 have attracted the interest of many investigators for drug repurposing and use in many pathological conditions, including cancer, neurodegenerative conditions, microbial conditions, renal conditions, free radical pathology, metal intoxication in relation to Fe, Cu, Al, Zn, Ga, In, U, and Pu, and other diseases. Similarly, the properties of L1 increase the prospects of its wider use in optimizing therapeutic efforts in many other fields of medicine, including synergies with other drugs.

Characterization of Mitochondrial DNA Methylation of Alzheimer's Disease in Plasma Cell-Free DNA

Noninvasive diagnosis of Alzheimer's disease (AD) is important for patients. Significant differences in the methylation of mitochondrial DNA (mtDNA) were found in AD brain tissue. Cell-free DNA (cfDNA) is a noninvasive and economical diagnostic tool. We aimed to characterize mtDNA methylation alterations in the plasma cfDNA of 31 AD patients and 26 age- and sex-matched cognitively normal control subjects. We found that the mtDNA methylation patterns differed between AD patients and control subjects. The mtDNA was predominantly hypomethylated in the plasma cfDNA of AD patients. The hypomethylation sites or regions were mainly located in mt-rRNA, mt-tRNA, and D-Loop regions. The hypomethylation of the D-Loop region in plasma cfDNA of AD patients was consistent with that in previous studies. This study presents evidence that hypomethylation in the non-protein coding region of mtDNA may contribute to the pathogenesis of AD and potential application for the diagnosis of AD.

TREC mediated oncogenesis in human immature T lymphoid malignancies preferentially involves ZFP36L2

The reintegration of excised signal joints resulting from human V(D)J recombination was described as a potent source of genomic instability in human lymphoid cancers. However, such molecular events have not been recurrently reported in clinical patient lymphoma/leukemia samples. Using a specifically designed NGS-capture pipeline, we here demonstrated the reintegration of T-cell receptor excision circles (TRECs) in 20/1533 (1.3%) patients with T-cell acute lymphoblastic leukemia (T-ALL) and T-cell lymphoblastic lymphoma (T-LBL). Remarkably, the reintegration of TREC recurrently targeted the tumor suppressor gene, ZFP36L2, in 17/20 samples. Thus, our data identified a new and hardly detectable mechanism of gene deregulation in lymphoid cancers providing new insights in human oncogenesis.

Enhanced mitochondrial genome analysis: bioinformatic and long-read sequencing advances and their diagnostic implications

INTRODUCTION

Primary mitochondrial diseases (PMDs) comprise a large and heterogeneous group of genetic diseases that result from pathogenic variants in either nuclear DNA (nDNA) or mitochondrial DNA (mtDNA). Widespread adoption of next-generation sequencing (NGS) has improved the efficiency and accuracy of mtDNA diagnoses; however, several challenges remain.

AREAS COVERED

In this review, we briefly summarize the current state of the art in molecular diagnostics for mtDNA and consider the implications of improved whole genome sequencing (WGS), bioinformatic techniques, and the adoption of long-read sequencing, for PMD diagnostics.

EXPERT OPINION

We anticipate that the application of PCR-free WGS from blood DNA will increase in diagnostic laboratories, while for adults with myopathic presentations, WGS from muscle DNA may become more widespread. Improved bioinformatic strategies will enhance WGS data interrogation, with more accurate delineation of mtDNA and NUMTs (nuclear mitochondrial DNA segments) in WGS data, superior coverage uniformity, indirect measurement of mtDNA copy number, and more accurate interpretation of heteroplasmic large-scale rearrangements (LSRs). Separately, the adoption of diagnostic long-read sequencing could offer greater resolution of complex LSRs and the opportunity to phase heteroplasmic variants.

Nanopore long-read next-generation sequencing for detection of mitochondrial DNA large-scale deletions

Primary mitochondrial diseases are progressive genetic disorders affecting multiple organs and characterized by mitochondrial dysfunction. These disorders can be caused by mutations in nuclear genes coding proteins with mitochondrial localization or by genetic defects in the mitochondrial genome (mtDNA). The latter include point pathogenic variants and large-scale deletions/rearrangements. MtDNA molecules with the wild type or a variant sequence can exist together in a single cell, a condition known as mtDNA heteroplasmy. MtDNA single point mutations are typically detected by means of Next-Generation Sequencing (NGS) based on short reads which, however, are limited for the identification of structural mtDNA alterations. Recently, new NGS technologies based on long reads have been released, allowing to obtain sequences of several kilobases in length; this approach is suitable for detection of structural alterations affecting the mitochondrial genome. In the present work we illustrate the optimization of two sequencing protocols based on long-read Oxford Nanopore Technology to detect mtDNA structural alterations. This approach presents strong advantages in the analysis of mtDNA compared to both short-read NGS and traditional techniques, potentially becoming the method of choice for genetic studies on mtDNA.

Independent and sensory human mitochondrial functions reflecting symbiotic evolution

The bacterial origin of mitochondria has been a widely accepted as an event that occurred about 1.45 billion years ago and endowed cells with internal energy producing organelle. Thus, mitochondria have traditionally been viewed as subcellular organelle as any other - fully functionally dependent on the cell it is a part of. However, recent studies have given us evidence that mitochondria are more functionally independent than other organelles, as they can function outside the cells, engage in complex "social" interactions, and communicate with each other as well as other cellular components, bacteria and viruses. Furthermore, mitochondria move, assemble and organize upon sensing different environmental cues, using a process akin to bacterial quorum sensing. Therefore, taking all these lines of evidence into account we hypothesize that mitochondria need to be viewed and studied from a perspective of a more functionally independent entity. This view of mitochondria may lead to new insights into their biological function, and inform new strategies for treatment of disease associated with mitochondrial dysfunction.

Mitochondrial Deoxyribonucleic Acid (mtDNA), Maternal Inheritance, and Their Role in the Development of Cancers: A Scoping Review

Mitochondrial DNA (mtDNA) is a small, circular, double-stranded DNA inherited from the mother during fertilization. Evolutionary evidence supported by the endosymbiotic theory identifies mitochondria as an organelle that could have descended from prokaryotes. This may be the reason for the independent function and inheritance pattern shown by mtDNA. The unstable nature of mtDNA due to the lack of protective histones, and effective repair systems make it more vulnerable to mutations. The mtDNA and its mutations could be maternally inherited thereby predisposing the offspring to various cancers like breast and ovarian cancers among others. Although mitochondria are considered heteroplasmic wherein variations among the multiple mtDNA genomes are noticed, mothers can have mitochondrial populations that are homoplasmic for a given mitochondrial mutation. Homoplasmic mitochondrial mutations may be transmitted to all maternal offspring. However, due to the complex interplay between the mitochondrial and nuclear genomes, it is often difficult to predict disease outcomes, even with homoplasmic mitochondrial populations. Heteroplasmic mtDNA mutations can be maternally inherited, but the proportion of mutated alleles differs markedly between offspring within one generation. This led to the genetic bottleneck hypothesis, explaining the rapid changes in allele frequency witnessed during the transmission of mtDNA from one generation to the next. Although a physical reduction in mtDNA has been demonstrated in several species, a comprehensive understanding of the molecular mechanisms is yet to be demonstrated. Despite initially thought to be limited to the germline, there is evidence that blockages exist in different cell types during development, perhaps explaining why different tissues in the same organism contain different levels of mutated mtDNA. In this review, we comprehensively discuss the potential mechanisms through which mtDNA undergoes mutations and the maternal mode of transmission that contributes to the development of tumors, especially breast and ovarian cancers.

Massive invasion of organellar DNA drives nuclear genome evolution in Toxoplasma

Toxoplasma gondii is a zoonotic protist pathogen that infects up to 1/3 of the human population. This apicomplexan parasite contains three genome sequences: nuclear (63 Mb); plastid organellar, ptDNA (35 kb); and mitochondrial organellar, mtDNA (5.9 kb of non-repetitive sequence). We find that the nuclear genome contains a significant amount of NUMTs (nuclear DNA of mitochondrial origin) and NUPTs (nuclear DNA of plastid origin) that are continuously acquired and represent a significant source of intraspecific genetic variation. NUOT (nuclear DNA of organellar origin) accretion has generated 1.6% of the extant T. gondii ME49 nuclear genome; the highest fraction ever reported in any organism. NUOTs are primarily found in organisms that retain the non-homologous end-joining repair pathway. Significant movement of organellar DNA was experimentally captured via amplicon sequencing of a CRISPR-induced double-strand break in non-homologous end-joining repair competent, but not ku80 mutant, Toxoplasma parasites. Comparisons with Neospora caninum, a species that diverged from Toxoplasma ~28 MY ago, revealed that the movement and fixation of 5 NUMTs predates the split of the two genera. This unexpected level of NUMT conservation suggests evolutionary constraint for cellular function. Most NUMT insertions reside within (60%) or nearby genes (23% within 1.5 kb) and reporter assays indicate that some NUMTs have the ability to function as cis-regulatory elements modulating gene expression. Together these findings portray a role for organellar sequence insertion in dynamically shaping the genomic architecture and likely contributing to adaptation and phenotypic changes in this important human pathogen.