Arguably, the most bizarre mitochondrial DNA (mtDNA) is that of the euglenozoan eukaryote Diplonema papillatum. The genome consists of numerous small circular chromosomes none of which appears to encode a complete gene. For instance, the cox1 coding sequence is spread out over nine different chromosomes in non-overlapping pieces (modules), which are transcribed separately and joined to a contiguous mRNA by trans-splicing. Here, we examine how many genes are encoded by Diplonema mtDNA and whether all are fragmented and their transcripts trans-spliced. Module identification is challenging due to the sequence divergence of Diplonema mitochondrial genes. By employing most sensitive protein profile search algorithms and comparing genomic with cDNA sequence, we recognize a total of 11 typical mitochondrial genes. The 10 protein-coding genes are systematically chopped up into three to 12 modules of 60-350 bp length. The corresponding mRNAs are all trans-spliced. Identification of ribosomal RNAs is most difficult. So far, we only detect the 3'-module of the large subunit ribosomal RNA (rRNA); it does not trans-splice with other pieces. The small subunit rRNA gene remains elusive. Our results open new intriguing questions about the biochemistry and evolution of mitochondrial trans-splicing in Diplonema.

The mitochondrial genome exists in numerous structural conformations, complicating the study of mitochondrial DNA (mtDNA) metabolism. Here, we describe the development of 2D intact mtDNA agarose gel electrophoresis (2D-IMAGE) for the separation and detection of approximately two-dozen distinct topoisomers. Although the major topoisomers were well conserved across many cell and tissue types, unique differences in certain cells and tissues were also observed. RNase treatment revealed that partially hybridized RNAs associated primarily with covalently closed circular DNA, consistent with this structure being the template for transcription. Circular structures composed of RNA:DNA hybrids contained only heavy-strand DNA sequences, implicating them as lagging-strand replication intermediates. During recovery from replicative arrest, 2D-IMAGE showed changes in both template selection and replication products. These studies suggest that discrete topoisomers are associated with specific mtDNA-directed processes. Because of the increased resolution, 2D-IMAGE has the potential to identify novel mtDNA intermediates involved in replication or transcription, or pathology including oxidative linearization, deletions or depletion.

Yeast Qri7 and human OSGEPL are members of the orthologous Kae1(OSGEP)/YgjD protein family, the last class of universally conserved proteins without assigned function. Phylogenetic analyses indicate that the eukaryotic Qri7(OSGEPL) proteins originated from bacterial YgjD proteins. We have recently shown that the archaeal Kae1 protein is a DNA-binding protein that exhibits apurinic endonuclease activity in vitro. We show here that the Qri7/OSGEPL proteins localize in mitochondria and are involved in mitochondrial genome maintenance in two model eukaryotic organisms, Saccharomyces cerevisiae and Caenorhabditis elegans. Furthermore, S. cerevisiae Qri7 complements the loss of the bacterial YgjD protein in Escherichia coli, suggesting that Qri7/OSGEPL and YgjD proteins have retained similar functions in modern organisms. We suggest to name members of the Kae1(OSGEP)/YgjD family UGMP, for Universal Genome Maintenance Proteins.

The mitochondrial genome (mtDNA) is assembled into nucleo-protein structures termed nucleoids and maintained differently compared to nuclear DNA, the involved molecular basis remaining poorly understood. In yeast (Saccharomyces cerevisiae), mtDNA is a ∼80 kbp linear molecule and Abf2p, a double HMG-box protein, packages and maintains it. The protein binds DNA in a non-sequence-specific manner, but displays a distinct 'phased-binding' at specific DNA sequences containing poly-adenine tracts (A-tracts). We present here two crystal structures of Abf2p in complex with mtDNA-derived fragments bearing A-tracts. Each HMG-box of Abf2p induces a 90° bend in the contacted DNA, causing an overall U-turn. Together with previous data, this suggests that U-turn formation is the universal mechanism underlying mtDNA compaction induced by HMG-box proteins. Combining this structural information with mutational, biophysical and computational analyses, we reveal a unique DNA binding mechanism for Abf2p where a characteristic N-terminal flag and helix are crucial for mtDNA maintenance. Additionally, we provide the molecular basis for A-tract mediated exclusion of Abf2p binding. Due to high prevalence of A-tracts in yeast mtDNA, this has critical relevance for nucleoid architecture. Therefore, an unprecedented A-tract mediated protein positioning mechanism regulates DNA packaging proteins in the mitochondria, and in combination with DNA-bending and U-turn formation, governs mtDNA compaction.

Mitochondrial genome (mitogenome) plays important roles in evolutionary and ecological studies. It becomes routine to utilize multiple genes on mitogenome or the entire mitogenomes to investigate phylogeny and biodiversity of focal groups with the onset of High Throughput Sequencing (HTS) technologies. We developed a mitogenome toolkit MitoZ, consisting of independent modules of de novo assembly, findMitoScaf (find Mitochondrial Scaffolds), annotation and visualization, that can generate mitogenome assembly together with annotation and visualization results from HTS raw reads. We evaluated its performance using a total of 50 samples of which mitogenomes are publicly available. The results showed that MitoZ can recover more full-length mitogenomes with higher accuracy compared to the other available mitogenome assemblers. Overall, MitoZ provides a one-click solution to construct the annotated mitogenome from HTS raw data and will facilitate large scale ecological and evolutionary studies. MitoZ is free open source software distributed under GPLv3 license and available at https://github.com/linzhi2013/MitoZ.

Genomic instability is related to a wide-range of human diseases. Here, we show that mitochondrial iron-sulfur cluster biosynthesis is important for the maintenance of nuclear genome stability in Saccharomyces cerevisiae. Cells lacking the mitochondrial chaperone Zim17 (Tim15/Hep1), a component of the iron-sulfur biosynthesis machinery, have limited respiration activity, mimic the metabolic response to iron starvation and suffer a dramatic increase in nuclear genome recombination. Increased oxidative damage or deficient DNA repair do not account for the observed genomic hyperrecombination. Impaired cell-cycle progression and genetic interactions of ZIM17 with components of the RFC-like complex involved in mitotic checkpoints indicate that replicative stress causes hyperrecombination in zim17Δ mutants. Furthermore, nuclear accumulation of pre-ribosomal particles in zim17Δ mutants reinforces the importance of iron-sulfur clusters in normal ribosome biosynthesis. We propose that compromised ribosome biosynthesis and cell-cycle progression are interconnected, together contributing to replicative stress and nuclear genome instability in zim17Δ mutants.

In quiescent fibroblasts, the expression levels of cytosolic enzymes for thymidine triphosphate (dTTP) synthesis are down-regulated, causing a marked reduction in the dTTP pool. In this study, we provide evidence that mitochondrial thymidylate synthesis via thymidine kinase 2 (TK2) is a limiting factor for the repair of ultraviolet (UV) damage in the nuclear compartment in quiescent fibroblasts. We found that TK2 deficiency causes secondary DNA double-strand breaks formation in the nuclear genome of quiescent cells at the late stage of recovery from UV damage. Despite slower repair of quiescent fibroblast deficient in TK2, DNA damage signals eventually disappeared, and these cells were capable of re-entering the S phase after serum stimulation. However, these cells displayed severe genome stress as revealed by the dramatic increase in 53BP1 nuclear body in the G1 phase of the successive cell cycle. Here, we conclude that mitochondrial thymidylate synthesis via TK2 plays a role in facilitating the quality repair of UV damage for the maintenance of genome integrity in the cells that are temporarily arrested in the quiescent state.

Nucleic acid-dependent ATPases are involved in nearly all aspects of DNA and RNA metabolism. Previous studies have described a number of mitochondrial helicases. However, double-stranded DNA-dependent ATPases, including translocases or enzymes remodeling DNA-protein complexes, have not been identified in mitochondria of the yeast Saccharomyces cerevisae. Here, we demonstrate that Irc3p is a mitochondrial double-stranded DNA-dependent ATPase of the Superfamily II. In contrast to the other mitochondrial Superfamily II enzymes Mss116p, Suv3p and Mrh4p, which are RNA helicases, Irc3p has a direct role in mitochondrial DNA (mtDNA) maintenance. Specific Irc3p-dependent mtDNA metabolic intermediates can be detected, including high levels of double-stranded DNA breaks that accumulate in irc3Δ mutants. irc3Δ-related topology changes in rho- mtDNA can be reversed by the deletion of mitochondrial RNA polymerase RPO41, suggesting that Irc3p counterbalances adverse effects of transcription on mitochondrial genome stability.

Xeroderma pigmentosum group D (XPD/ERCC2) encodes an ATP-dependent helicase that plays essential roles in both transcription and nucleotide excision repair of nuclear DNA, however, whether or not XPD exerts similar functions in mitochondria remains elusive. In this study, we provide the first evidence that XPD is localized in the inner membrane of mitochondria, and cells under oxidative stress showed an enhanced recruitment of XPD into mitochondrial compartment. Furthermore, mitochondrial reactive oxygen species production and levels of oxidative stress-induced mitochondrial DNA (mtDNA) common deletion were significantly elevated, whereas capacity for oxidative damage repair of mtDNA was markedly reduced in both XPD-suppressed human osteosarcoma (U2OS) cells and XPD-deficient human fibroblasts. Immunoprecipitation-mass spectrometry analysis was used to identify interacting factor(s) with XPD and TUFM, a mitochondrial Tu translation elongation factor was detected to be physically interacted with XPD. Similar to the findings in XPD-deficient cells, mitochondrial common deletion and oxidative damage repair capacity in U2OS cells were found to be significantly altered after TUFM knock-down. Our findings clearly demonstrate that XPD plays crucial role(s) in protecting mitochondrial genome stability by facilitating an efficient repair of oxidative DNA damage in mitochondria.

Transfer RNAs (tRNAs) are present in all types of cells as well as in organelles. tRNAs of animal mitochondria show a low level of primary sequence conservation and exhibit 'bizarre' secondary structures, lacking complete domains of the common cloverleaf. Such sequences are hard to detect and hence frequently missed in computational analyses and mitochondrial genome annotation. Here, we introduce an automatic annotation procedure for mitochondrial tRNA genes in Metazoa based on sequence and structural information in manually curated covariance models. The method, applied to re-annotate 1876 available metazoan mitochondrial RefSeq genomes, allows to distinguish between remaining functional genes and degrading 'pseudogenes', even at early stages of divergence. The subsequent analysis of a comprehensive set of mitochondrial tRNA genes gives new insights into the evolution of structures of mitochondrial tRNA sequences as well as into the mechanisms of genome rearrangements. We find frequent losses of tRNA genes concentrated in basal Metazoa, frequent independent losses of individual parts of tRNA genes, particularly in Arthropoda, and wide-spread conserved overlaps of tRNAs in opposite reading direction. Direct evidence for several recent Tandem Duplication-Random Loss events is gained, demonstrating that this mechanism has an impact on the appearance of new mitochondrial gene orders.

To study the mechanisms involved in the maintenance of a linear mitochondrial genome we investigated the biochemical properties of the recombination protein Mgm101 from Candida parapsilosis. We show that CpMgm101 complements defects associated with the Saccharomyces cerevisiae mgm101-1(ts) mutation and that it is present in both the nucleus and mitochondrial nucleoids of C. parapsilosis. Unlike its S. cerevisiae counterpart, CpMgm101 is associated with the entire nucleoid population and is able to bind to a broad range of DNA substrates in a non-sequence specific manner. CpMgm101 is also able to catalyze strand annealing and D-loop formation. CpMgm101 forms a roughly C-shaped trimer in solution according to SAXS. Electron microscopy of a complex of CpMgm101 with a model mitochondrial telomere revealed homogeneous, ring-shaped structures at the telomeric single-stranded overhangs. The DNA-binding properties of CpMgm101, together with its DNA recombination properties, suggest that it can play a number of possible roles in the replication of the mitochondrial genome and the maintenance of its telomeres.

Uridine insertion/deletion (U-indel) editing of mitochondrial mRNA, unique to the protistan class Kinetoplastea, generates canonical as well as potentially non-productive editing events. While the molecular machinery and the role of the guide (g) RNAs that provide required information for U-indel editing are well understood, little is known about the forces underlying its apparently error-prone nature. Analysis of a gRNA:mRNA pair allows the dissection of editing events in a given position of a given mitochondrial transcript. A complete gRNA dataset, paired with a fully characterized mRNA population that includes non-canonically edited transcripts, would allow such an analysis to be performed globally across the mitochondrial transcriptome. To achieve this, we have assembled 67 minicircles of the insect parasite Leptomonas pyrrhocoris, with each minicircle typically encoding one gRNA located in one of two similar-sized units of different origin. From this relatively narrow set of annotated gRNAs, we have dissected all identified mitochondrial editing events in L. pyrrhocoris, the strains of which dramatically differ in the abundance of individual minicircle classes. Our results support a model in which a multitude of editing events are driven by a limited set of gRNAs, with individual gRNAs possessing an inherent ability to guide canonical and non-canonical editing.

Kinetoplastids are protists defined by one of the most complex mitochondrial genomes in nature, the kinetoplast. In the sleeping sickness parasite Trypanosoma brucei, the kinetoplast is a chain mail-like network of two types of interlocked DNA molecules: a few dozen ∼23-kb maxicircles (homologs of the mitochondrial genome of other eukaryotes) and thousands of ∼1-kb minicircles. Maxicircles encode components of respiratory chain complexes and the mitoribosome. Several maxicircle-encoded mRNAs undergo extensive post-transcriptional RNA editing via addition and deletion of uridines. The process is mediated by hundreds of species of minicircle-encoded guide RNAs (gRNAs), but the precise number of minicircle classes and gRNA genes was unknown. Here we present the first essentially complete assembly and annotation of the kinetoplast genome of T. brucei. We have identified 391 minicircles, encoding not only ∼930 predicted 'canonical' gRNA genes that cover nearly all known editing events (accessible via the web at http://hank.bio.ed.ac.uk), but also ∼370 'non-canonical' gRNA genes of unknown function. Small RNA transcriptome data confirmed expression of the majority of both categories of gRNAs. Finally, we have used our data set to refine definitions for minicircle structure and to explore dynamics of minicircle copy numbers.

Pathological conditions impairing functions of mitochondria often lead to compensatory upregulation of the mitochondrial DNA (mtDNA) replisome machinery, and the replicative DNA helicase appears to be a key factor in regulating mtDNA copy number. Moreover, mtDNA helicase mutations have been associated with structural rearrangements of the mitochondrial genome. To evaluate the effects of elevated levels of the mtDNA helicase on the integrity and replication of the mitochondrial genome, we overexpressed the helicase in Drosophila melanogaster Schneider cells and analyzed the mtDNA by two-dimensional neutral agarose gel electrophoresis and electron microscopy. We found that elevation of mtDNA helicase levels increases the quantity of replication intermediates and alleviates pausing at the replication slow zones. Though we did not observe a concomitant alteration in mtDNA copy number, we observed deletions specific to the segment of repeated elements in the immediate vicinity of the origin of replication, and an accumulation of species characteristic of replication fork stalling. We also found elevated levels of RNA that are retained in the replication intermediates. Together, our results suggest that upregulation of mtDNA helicase promotes the process of mtDNA replication but also results in genome destabilization.

N6-methyldeoxyadenosine (6mA) is a well-characterized DNA modification in prokaryotes but reports on its presence and function in mammals have been controversial. To address this issue, we established the capacity of 6mA-Crosslinking-Exonuclease-sequencing (6mACE-seq) to detect genome-wide 6mA at single-nucleotide-resolution, demonstrating this by accurately mapping 6mA in synthesized DNA and bacterial genomes. Using 6mACE-seq, we generated a human-genome-wide 6mA map that accurately reproduced known 6mA enrichment at active retrotransposons and revealed mitochondrial chromosome-wide 6mA clusters asymmetrically enriched on the heavy-strand. We identified a novel putative 6mA-binding protein in single-stranded DNA-binding protein 1 (SSBP1), a mitochondrial DNA (mtDNA) replication factor known to coat the heavy-strand, linking 6mA with the regulation of mtDNA replication. Finally, we characterized AlkB homologue 1 (ALKBH1) as a mitochondrial protein with 6mA demethylase activity and showed that its loss decreases mitochondrial oxidative phosphorylation. Our results show that 6mA clusters play a previously unappreciated role in regulating human mitochondrial function, despite 6mA being an uncommon DNA modification in the human genome.

The maintenance of the mitochondrial genomic integrity is a prerequisite for proper mitochondrial function. Due to the high concentration of reactive oxygen species (ROS) generated by the oxidative phosphorylation pathway, the mitochondrial genome is highly exposed to oxidative stress leading to mitochondrial DNA injury. Accordingly, mitochondrial DNA damage was shown to be associated with ageing as well as with numerous human diseases including neurodegenerative disorders and cancer. To date, several methods have been described to detect damaged mitochondrial DNA, but those techniques are semi-quantitative and often require high amounts of genomic input DNA. We developed a rapid and quantitative method to evaluate the relative levels of damage in mitochondrial DNA by using the real time-PCR amplification of mitochondrial DNA fragments of different lengths. We investigated mitochondrial DNA damage in SH-SY5Y human neuroblastoma cells exposed to hydrogen peroxide or stressed by over-expression of the tyrosinase gene. In the past, there has been speculation about a variable vulnerability to oxidative stress along the mitochondrial genome. Our results indicate the existence of at least one mitochondrial DNA hot spot, namely the D-Loop, being more prone to ROS-derived damage.

Plant mitochondrial genomes show much more evolutionary plasticity than those of animals. We analysed the first mitochondrial DNA (mtDNA) of a lycophyte, the quillwort Isoetes engelmannii, which is separated from seed plants by more than 350 million years of evolution. The Isoetes mtDNA is particularly rich in recombination events, and chloroplast as well as nuclear DNA inserts document the incorporation of foreign sequences already in this most ancestral vascular plant lineage. On the other hand, particularly small group II introns and short intergenic regions reveal a tendency of evolution towards a compact mitochondrial genome. RNA editing reaches extreme levels exceeding 100 pyrimidine exchanges in individual mRNAs and, hitherto unobserved in such frequency, also in tRNAs with 18 C-to-U conversions in the tRNA for proline. In total, some 1500 sites of RNA editing can be expected for the Isoetes mitochondrial transcriptome. As a unique molecular novelty, the Isoetes cox1 gene requires trans-splicing via a discontinuous group I intron demonstrating disrupted, but functional, RNAs for yet another class of natural ribozymes.

DNA damage is tightly associated with various biological and pathological processes, such as aging and tumorigenesis. Although detection of DNA damage is attracting increasing attention, only a limited number of methods are available to quantify DNA lesions, and these techniques are tedious or only detect global DNA damage. In this study, we present a high-sensitivity long-run real-time PCR technique for DNA-damage quantification (LORD-Q) in both the mitochondrial and nuclear genome. While most conventional methods are of low-sensitivity or restricted to abundant mitochondrial DNA samples, we established a protocol that enables the accurate sequence-specific quantification of DNA damage in >3-kb probes for any mitochondrial or nuclear DNA sequence. In order to validate the sensitivity of this method, we compared LORD-Q with a previously published qPCR-based method and the standard single-cell gel electrophoresis assay, demonstrating a superior performance of LORD-Q. Exemplarily, we monitored induction of DNA damage and repair processes in human induced pluripotent stem cells and isogenic fibroblasts. Our results suggest that LORD-Q provides a sequence-specific and precise method to quantify DNA damage, thereby allowing the high-throughput assessment of DNA repair, genotoxicity screening and various other processes for a wide range of life science applications.

Mutations in mitochondrial DNA (mtDNA) cause maternally inherited diseases, while somatic mutations are linked to common diseases of aging. Although mtDNA mutations impact health, the processes that give rise to them are under considerable debate. To investigate the mechanism by which de novo mutations arise, we analyzed the distribution of naturally occurring somatic mutations across the mouse and human mtDNA obtained by Duplex Sequencing. We observe distinct mutational gradients in G→A and T→C transitions delimited by the light-strand origin and the mitochondrial Control Region (mCR). The gradient increases unequally across the mtDNA with age and is lost in the absence of DNA polymerase γ proofreading activity. In addition, high-resolution analysis of the mCR shows that important regulatory elements exhibit considerable variability in mutation frequency, consistent with them being mutational 'hot-spots' or 'cold-spots'. Collectively, these patterns support genome replication via a deamination prone asymmetric strand-displacement mechanism as the fundamental driver of mutagenesis in mammalian DNA. Moreover, the distribution of mtDNA single nucleotide polymorphisms in humans and the distribution of bases in the mtDNA across vertebrate species mirror this gradient, indicating that replication-linked mutations are likely the primary source of inherited polymorphisms that, over evolutionary timescales, influences genome composition during speciation.

Abundant ribonucleoside-triphosphate (rNTP) incorporation into DNA by DNA polymerases in the form of ribonucleoside monophosphates (rNMPs) is a widespread phenomenon in nature, resulting in DNA-structural change and genome instability. The rNMP distribution, characteristics, hotspots and association with DNA metabolic processes in human mitochondrial DNA (hmtDNA) remain mostly unknown. Here, we utilize the ribose-seq technique to capture embedded rNMPs in hmtDNA of six different cell types. In most cell types, the rNMPs are preferentially embedded on the light strand of hmtDNA with a strong bias towards rCMPs; while in the liver-tissue cells, the rNMPs are predominately found on the heavy strand. We uncover common rNMP hotspots and conserved rNMP-enriched zones across the entire hmtDNA, including in the control region, which links the rNMP presence to the frequent hmtDNA replication-failure events. We show a strong correlation between coding-sequence size and rNMP-embedment frequency per nucleotide on the non-template, light strand in all cell types, supporting the presence of transient RNA-DNA hybrids preceding light-strand replication. Moreover, we detect rNMP-embedment patterns that are only partly conserved across the different cell types and are distinct from those found in yeast mtDNA. The study opens new research directions to understand the biology of hmtDNA and genomic rNMPs.