Mammalian NUMT insertion is non-random

Nuclear sequences of mitochondrial origin in domestic yak

Mitochondrial DNA sequences are frequently transferred into the nuclear genome, generating nuclear mitochondrial DNA sequences (NUMTs). Here, we analysed, for the first time, NUMTs in the domestic yak genome. We obtained 499 alignment matches covering 340.2 kbp of the yak nuclear genome. After a merging step, we identified 167 NUMT regions with a total length of ~ 503 kbp, representing 0.02% of the nuclear genome. We discovered copies of all mitochondrial regions and found that most NUMT regions are intergenic or intronic and mostly untranscribed. 98 different NUMT regions from domestic yak showed high homology with cow and/or wild yak genomes, suggesting selection or hybridization between domestic/wild yak and cow. To rule out the possibility that the identified NUMTs could be artifacts of the domestic yak genome assembly, we validated experimentally five NUMT regions by PCR amplification. As NUMT regions show high similarity to the mitochondrial genome can potentially pose a risk to domestic yak DNA mitochondrial studies, special care is therefore needed to select primers for PCR amplification of mitochondrial DNA sequences.

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.

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.

The mitogenomes of Leptographium aureum, Leptographium sp., and Grosmannia fruticeta: expansion by introns

Introduction

Many members of the Ophiostomatales are of economic importance as they are bark-beetle associates and causative agents for blue stain on timber and in some instances contribute towards tree mortality. The taxonomy of these fungi has been challenging due to the convergent evolution of many traits associated with insect dispersal and a limited number of morphological characters that happen to be highly pleomorphic. This study examines the mitochondrial genomes for three members of Leptographium sensu lato [Leptographium aureum (also known as Grosmannia aurea), Grosmannia fruticeta (also known as Leptographium fruticetum), and Leptographium sp. WIN(M)1376)].

Methods

Illumina sequencing combined with gene and intron annotations and phylogenetic analysis were performed.

Results

Sequence analysis showed that gene content and gene synteny are conserved but mitochondrial genome sizes were variable: G. fruticeta at 63,821 bp, Leptographium sp. WIN(M)1376 at 81,823 bp and L. aureum at 104,547 bp. The variation in size is due to the number of introns and intron-associated open reading frames. Phylogenetic analysis of currently available mitochondrial genomes for members of the Ophiostomatales supports currently accepted generic arrangements within this order and specifically supports the separation of members with Leptographium-like conidiophores into two genera, with L. aureum grouping with Leptographium and G. fruticeta aligning with Grosmannia.

Discussion

Mitochondrial genomes are promising sequences for resolving evolutionary relationships within the Ophiostomatales.

The Mighty NUMT: Mitochondrial DNA Flexing Its Code in the Nuclear Genome

Nuclear-mitochondrial DNA segments (NUMTs) are mitochondrial DNA (mtDNA) fragments that have been inserted into the nuclear genome. Some NUMTs are common within the human population but most NUMTs are rare and specific to individuals. NUMTs range in size from 24 base pairs to encompassing nearly the entire mtDNA and are found throughout the nuclear genome. Emerging evidence suggests that the formation of NUMTs is an ongoing process in humans. NUMTs contaminate sequencing results of the mtDNA by introducing false positive variants, particularly heteroplasmic variants present at a low variant allele frequency (VAF). In our review, we discuss the prevalence of NUMTs in the human population, the potential mechanisms of de novo NUMT insertion via DNA repair mechanisms, and provide an overview of the existing approaches for minimizing NUMT contamination. Apart from filtering known NUMTs, both wet lab-based and computational methods can be used to minimize the contamination of NUMTs in analyses of human mtDNA. Current approaches include: (1) isolating mitochondria to enrich for mtDNA; (2) applying basic local alignment to identify NUMTs for subsequent filtering; (3) bioinformatic pipelines for NUMT detection; (4) k-mer-based NUMT detection; and (5) filtering candidate false positive variants by mtDNA copy number, VAF, or sequence quality score. Multiple approaches must be applied in order to effectively identify NUMTs in samples. Although next-generation sequencing is revolutionizing our understanding of heteroplasmic mtDNA, it also raises new challenges with the high prevalence and individual-specific NUMTs that need to be handled with care in studies of mitochondrial genetics.

Extensive mitogenomic heteroplasmy and its implications in the phylogeny of the fish genus Megalobrama

Megalobrama is China's most economically valuable fish genera. Four species make up this genus: M. amblycephala (MA), M. terminalis (MT), M. pellegrini (MP), and M. hoffmanni (MH). Many researchers have investigated the genetic relationship of Megalobrama based on mitochondrial DNA (mtDNA) and discovered that the branches of the phylogenetic tree for MT and MP are intertwined. We hypothesized that this occurs because mitogenomic heteroplasmy is overlooked when working with mtDNA, which causes MP and MT positions to intersect in phylogenetic trees. To eliminate the influence of nuclear mitochondrial DNA fragments (NUMTs) before analyzing mitogenomic heteroplasmy, we used PLastZ to identify NUMTs, which were then removed from the samples for the subsequent heteroplasmy analysis. Using the heteroplasmy caller icHET, we discovered 126, 339, 135, and 203 heteroplasmic variants in six MA, MT, MP, and MH samples. We reconstructed the Megalobrama fish genus's phylogenetic tree using the RY coding method and rejecting the third position on codons, which improved the performance of the phylogenetic tree by increasing the ratio of treeness to relative component variability from 100.02 ± 1.76 to 688.59 ± 190.56. Despite this, the RY coding method cannot alter the intersection of MP and MT positions in phylogenetic trees. We hypothesize that gene flow between MT and MP leads to intertwining mtDNA-based phylogenetic trees. In conclusion, our findings on the mitogenomic heteroplasmy of Megalobrama provide new insights into mtDNA-based phylogenetic studies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03523-0.

Comparison of detection methods and genome quality when quantifying nuclear mitochondrial insertions in vertebrate genomes

The integration of mitochondrial genome fragments into the nuclear genome is well documented, and the transfer of these mitochondrial nuclear pseudogenes (numts) is thought to be an ongoing evolutionary process. With the increasing number of eukaryotic genomes available, genome-wide distributions of numts are often surveyed. However, inconsistencies in genome quality can reduce the accuracy of numt estimates, and methods used for identification can be complicated by the diverse sizes and ages of numts. Numts have been previously characterized in rodent genomes and it was postulated that they might be more prevalent in a group of voles with rapidly evolving karyotypes. Here, we examine 37 rodent genomes, and an additional 26 vertebrate genomes, while also considering numt detection methods. We identify numts using DNA:DNA and protein:translated-DNA similarity searches and compare numt distributions among rodent and vertebrate taxa to assess whether some groups are more susceptible to transfer. A combination of protein sequence comparisons (protein:translated-DNA) and BLASTN genomic DNA searches detect 50% more numts than genomic DNA:DNA searches alone. In addition, higher-quality RefSeq genomes produce lower estimates of numts than GenBank genomes, suggesting that lower quality genome assemblies can overestimate numts abundance. Phylogenetic analysis shows that mitochondrial transfers are not associated with karyotypic diversity among rodents. Surprisingly, we did not find a strong correlation between numt counts and genome size. Estimates using DNA: DNA analyses can underestimate the amount of mitochondrial DNA that is transferred to the nucleus.

Trichoderma - genomes and genomics as treasure troves for research towards biology, biotechnology and agriculture

The genus Trichoderma is among the best studied groups of filamentous fungi, largely because of its high relevance in applications from agriculture to enzyme biosynthesis to biofuel production. However, the physiological competences of these fungi, that led to these beneficial applications are intriguing also from a scientific and ecological point of view. This review therefore summarizes recent developments in studies of fungal genomes, updates on previously started genome annotation efforts and novel discoveries as well as efforts towards bioprospecting for enzymes and bioactive compounds such as cellulases, enzymes degrading xenobiotics and metabolites with potential pharmaceutical value. Thereby insights are provided into genomes, mitochondrial genomes and genomes of mycoviruses of Trichoderma strains relevant for enzyme production, biocontrol and mycoremediation. In several cases, production of bioactive compounds could be associated with responsible genes or clusters and bioremediation capabilities could be supported or predicted using genome information. Insights into evolution of the genus Trichoderma revealed large scale horizontal gene transfer, predominantly of CAZyme genes, but also secondary metabolite clusters. Investigation of sexual development showed that Trichoderma species are competent of repeat induced point mutation (RIP) and in some cases, segmental aneuploidy was observed. Some random mutants finally gave away their crucial mutations like T. reesei QM9978 and QM9136 and the fertility defect of QM6a was traced back to its gene defect. The Trichoderma core genome was narrowed down to 7000 genes and gene clustering was investigated in the genomes of multiple species. Finally, recent developments in application of CRISPR/Cas9 in Trichoderma, cloning and expression strategies for the workhorse T. reesei as well as the use genome mining tools for bioprospecting Trichoderma are highlighted. The intriguing new findings on evolution, genomics and physiology highlight emerging trends and illustrate worthwhile perspectives in diverse fields of research with Trichoderma.

Nuclear mitochondrial DNA sequences in the rabbit genome

Numtogenesis is observable in the mammalian genomes resulting in the integration of mitochondrial segments into the nuclear genomes (numts). To identify numts in rabbit, we aligned mitochondrial and nuclear genomes. Alignment significance threshold was calculated and individual characteristics of numts were analysed. We found 153 numts in the nuclear genome. The GC content of numts were significantly lower than the GC content of their genomic flanking regions or the genome itself. The frequency of three mammalian-wide interspersed repeats were increased in the proximity of numts. The decreased GC content around numts strengthen the theory which supposes a link between DNA structural instability and numt integration.

Signatures of selection on mitonuclear integrated genes uncover hidden mitogenomic variation in fur seals

Nuclear copies of mitochondrial genes (numts) are commonplace in vertebrate genomes and have been characterized in many species. However, relatively little attention has been paid to understanding their evolutionary origins and to disentangling alternative sources of insertions. Numts containing genes with intact mitochondrial reading frames represent good candidates for this purpose. The sequences of the genes they contain can be compared to their mitochondrial homologs to characterize synonymous to non-synonymous substitution rates, which can shed light on the selection pressures these genes have been subjected to. Here, we characterise 25 numts in the Antarctic fur seal (Arctocephalus gazella) genome. Among those containing genes with intact mitochondrial reading frames, three carry multiple substitutions in comparison to their mitochondrial homologs. Our analyses reveal that one represents a historic insertion subjected to strong purifying selection since it colonized the Otarioidea in a genomic region enriched in retrotransposons. By contrast, the other two numts appear to be more recent and their large number of substitutions can be attributed to non-canonical insertions, either the integration of heteroplasmic mtDNA or hybridization. Our study sheds new light on the evolutionary history of pinniped numts and uncovers the presence of hidden sources of mitonuclear variation.

Orphan gene in Littorina: An unexpected role of symbionts in the host evolution

Mechanisms of reproductive isolation between closely related sympatric species are of high evolutionary significance as they may function as initial drivers of speciation and protect species integrity afterwards. Proteins involved in the establishment of reproductive barriers often evolve fast and may be key players in cessation of gene flow between the incipient species. The five Atlantic Littorina (Neritrema) species represent a notable example of recent radiation. The geographic ranges of these young species largely overlap and the mechanisms of reproductive isolation are poorly understood. In this study, we performed a detailed analysis of the reproductive protein LOSP, previously identified in Littorina. We showed that this protein is evolutionary young and taxonomically restricted to the genus Littorina. It has high sequence variation both within and between Littorina species, which is compatible with its presumable role in the reproductive isolation. The strongest differences in the LOSP structure were detected between Littorina subgenera with distinctive repetitive motifs present exclusively in the Neritrema species, but not in L. littorea. Moreover, the sequence of these repetitive structural elements demonstrates a high homology with genetic elements of bacteria, identified as components of Littorina associated microbiomes. We suggest that these elements were acquired from a symbiotic bacterial donor via horizontal genetic transfer (HGT), which is indirectly confirmed by the presence of multiple transposable elements in the LOSP flanking and intronic regions. Furthermore, we hypothesize that this HGT-driven evolutionary innovation promoted LOSP function in reproductive isolation, which might be one of the factors determining the intensive cladogenesis in the Littorina (Neritrema) lineage in contrast to the anagenesis in the L. littorea clade.

Diet, Trimethylamine Metabolism, and Mitochondrial DNA: An Observational Study

SCOPE

Mitochondrial DNA copy number (mtDNAcn) and its methylation level in the D-loop area have been correlated with metabolic health and are suggested to vary in response to environmental stimuli, including diet. Circulating levels of trimethylamine-n-oxide (TMAO), which is an oxidative derivative of the trimethylamine (TMA) produced by the gut microbiome from dietary precursors, have been associated with chronic diseases and are suggested to have an impact on mitochondrial dynamics. This study is aimed to investigate the relationship between diet, TMA, TMAO, and mtDNAcn, as well as DNA methylation.

METHODS AND RESULTS

Two hundred subjects with extreme (healthy and unhealthy) dietary patterns are recruited. Dietary records are collected to assess their nutrient intake and diets' quality (Healthy Eating Index). Blood levels of TMA and TMAO, circulating levels of TMA precursors and their dietary intakes are measured. MtDNAcn, nuclear DNA methylation long interspersed nuclear element 1 (LINE-1), and strand-specific D-loop methylation levels are assessed. There is no association between dietary patterns and mtDNAcn. The TMAO/TMA ratio is negatively correlated with d-loop methylation levels but positively with mtDNAcn.

CONCLUSIONS

These findings suggest a potential association between TMA metabolism and mitochondrial dynamics (and mtDNA), indicating a new avenue for further research.

A Burst of Numt Insertion in the Dasyuridae Family During Marsupial Evolution

Nuclear pseudogenes of mitochondrial origin (numts) are common in all eukaryotes. Our previous scan of numts in sequenced nuclear genomes suggested that the highest numt content currently known in animals is that in the gray short-tailed opossum. The present work sought to determine numt content in marsupials and to compare it to those in placental and monothematic mammals as well as in non-mammalian vertebrates. To achieve this, 70 vertebrate species with available nuclear and mitochondrial genomes were scanned for numt content. An extreme numt content was found in the Dasyuridae, with 3,450 in Sarcophilus harrisii (1,955 kb) and 2,813 in Antechinus flavipes (847 kb). The evolutionarily closest species analyzed, the extinct Thylacinus cynocephalus belonging to the Thylacindae family, had only 435 numts (238 kb). These two Dasyuridae genomes featured the highest numt content identified in animals to date. A phylogenetic analysis of numts longer than 300 bp, using a Diprotodonita mitochondrial tree, indicated a burst of numt insertion that began before the divergence of the Dasyurini and Phascogalini, reaching a peak in the early evolution of the two tribes. No comparable increase was found in the early divergent species T. cynocephalus. Divergence of the Dasyuridae tribes has been previously dated to shortly after the Miocene climate transition, characterized by a rapid temperature decline. Interestingly, deviation from optimal growth temperature is one of the environmental factors reported to increase numt insertions in a laboratory setting.

Mitochondrial DNA and Epigenetics: Investigating Interactions with the One-Carbon Metabolism in Obesity

Mitochondrial DNA copy number (mtDNAcn) has been proposed for use as a surrogate biomarker of mitochondrial health, and evidence suggests that mtDNA might be methylated. Intermediates of the one-carbon cycle (1CC), which is duplicated in the cytoplasm and mitochondria, have a major role in modulating the impact of diet on the epigenome. Moreover, epigenetic pathways and the redox system are linked by the metabolism of glutathione (GSH). In a cohort of 101 normal-weight and 97 overweight/obese subjects, we evaluated mtDNAcn and methylation levels in both mitochondrial and nuclear areas to test the association of these marks with body weight, metabolic profile, and availability of 1CC intermediates associated with diet. Body composition was associated with 1CC intermediate availability. Reduced levels of GSH were measured in the overweight/obese group (p = 1.3∗10⁻⁵). A high BMI was associated with lower LINE-1 (p = 0.004) and nominally lower methylenetetrahydrofolate reductase (MTHFR) gene methylation (p = 0.047). mtDNAcn was lower in overweight/obese subjects (p = 0.004) and independently correlated with MTHFR methylation levels (p = 0.005) but not to LINE-1 methylation levels (p = 0.086). DNA methylation has been detected in the light strand but not in the heavy strand of the mtDNA. Although mtDNA methylation in the light strand did not differ between overweight/obese and normal-weight subjects, it was nominally correlated with homocysteine levels (p = 0.035) and MTHFR methylation (p = 0.033). This evidence suggests that increased body weight might perturb mitochondrial-nuclear homeostasis affecting the availability of nutrients acting as intermediates of the one-carbon cycle.

Complete Genome Sequences and Genome-Wide Characterization of Trichoderma Biocontrol Agents Provide New Insights into their Evolution and Variation in Genome Organization, Sexual Development, and Fungal-Plant Interactions

Trichoderma spp. represent one of the most important fungal genera to mankind and in natural environments. The genus harbors prolific producers of wood-decaying enzymes, biocontrol agents against plant pathogens, plant-growth-promoting biofertilizers, as well as model organisms for studying fungal-plant-plant pathogen interactions. Pursuing highly accurate, contiguous, and chromosome-level reference genomes has become a primary goal of fungal research communities. Here, we report the chromosome-level genomic sequences and whole-genome annotation data sets of four strains used as biocontrol agents or biofertilizers (Trichoderma virens Gv29-8, Trichoderma virens FT-333, Trichoderma asperellum FT-101, and Trichoderma atroviride P1). Our results provide comprehensive categorization, correct positioning, and evolutionary detail of both nuclear and mitochondrial genomes, including telomeres, AT-rich blocks, centromeres, transposons, mating-type loci, nuclear-encoded mitochondrial sequences, as well as many new secondary metabolic and carbohydrate-active enzyme gene clusters. We have also identified evolutionarily conserved core genes contributing to plant-fungal interactions, as well as variations potentially linked to key behavioral traits such as sex, genome defense, secondary metabolism, and mycoparasitism. The genomic resources we provide herein significantly extend our knowledge not only of this economically important fungal genus, but also fungal evolution and basic biology in general. IMPORTANCE Telomere-to-telomere and gapless reference genome assemblies are necessary to ensure that all genomic variants are studied and discovered, including centromeres, telomeres, AT-rich blocks, mating type loci, biosynthetic, and metabolic gene clusters. Here, we applied long-range sequencing technologies to determine the near-completed genome sequences of four widely used biocontrol agents or biofertilizers: Trichoderma virens Gv29-8 and FT-333, Trichoderma asperellum FT-101, and Trichoderma atroviride P1. Like those of three Trichoderma reesei wild isolates [QM6a, CBS999.97(MAT1-1) and CBS999.97(MAT1-2)] we reported previously, these four biocontrol agent genomes each contain seven nuclear chromosomes and a circular mitochondrial genome. Substantial intraspecies and intragenus diversities are also discovered, including single nucleotide polymorphisms, chromosome shuffling, as well as genomic relics derived from historical transposition events and repeat-induced point (RIP) mutations.

Global Characterization of Fungal Mitogenomes: New Insights on Genomic Diversity and Dynamism of Coding Genes and Accessory Elements

Fungi comprise a great diversity of species with distinct ecological functions and lifestyles. Similar to other eukaryotes, fungi rely on interactions with prokaryotes and one of the most important symbiotic events was the acquisition of mitochondria. Mitochondria are organelles found in eukaryotic cells whose main function is to generate energy through aerobic respiration. Mitogenomes (mtDNAs) are double-stranded circular or linear DNA from mitochondria that may contain core genes and accessory elements that can be replicated, transcribed, and independently translated from the nuclear genome. Despite their importance, investigative studies on the diversity of fungal mitogenomes are scarce. Herein, we have evaluated 788 curated fungal mitogenomes available at NCBI database to assess discrepancies and similarities among them and to better understand the mechanisms involved in fungal mtDNAs variability. From a total of 12 fungal phyla, four do not have any representative with available mitogenomes, which highlights the underrepresentation of some groups in the current available data. We selected representative and non-redundant mitogenomes based on the threshold of 90% similarity, eliminating 81 mtDNAs. Comparative analyses revealed considerable size variability of mtDNAs with a difference of up to 260 kb in length. Furthermore, variation in mitogenome length and genomic composition are generally related to the number and length of accessory elements (introns, HEGs, and uORFs). We identified an overall average of 8.0 (0–39) introns, 8.0 (0–100) HEGs, and 8.2 (0–102) uORFs per genome, with high variation among phyla. Even though the length of the core protein-coding genes is considerably conserved, approximately 36.3% of the mitogenomes evaluated have at least one of the 14 core coding genes absent. Also, our results revealed that there is not even a single gene shared among all mitogenomes. Other unusual genes in mitogenomes were also detected in many mitogenomes, such as dpo and rpo, and displayed diverse evolutionary histories. Altogether, the results presented in this study suggest that fungal mitogenomes are diverse, contain accessory elements and are absent of a conserved gene that can be used for the taxonomic classification of the Kingdom Fungi.

Unlocking the Complexity of Mitochondrial DNA: A Key to Understanding Neurodegenerative Disease Caused by Injury

Neurodegenerative disorders that are triggered by injury typically have variable and unpredictable outcomes due to the complex and multifactorial cascade of events following the injury and during recovery. Hence, several factors beyond the initial injury likely contribute to the disease progression and pathology, and among these are genetic factors. Genetics is a recognized factor in determining the outcome of common neurodegenerative diseases. The role of mitochondrial genetics and function in traditional neurodegenerative diseases, such as Alzheimer’s and Parkinson’s diseases, is well-established. Much less is known about mitochondrial genetics, however, regarding neurodegenerative diseases that result from injuries such as traumatic brain injury and ischaemic stroke. We discuss the potential role of mitochondrial DNA genetics in the progression and outcome of injury-related neurodegenerative diseases. We present a guide for understanding mitochondrial genetic variation, along with the nuances of quantifying mitochondrial DNA variation. Evidence supporting a role for mitochondrial DNA as a risk factor for neurodegenerative disease is also reviewed and examined. Further research into the impact of mitochondrial DNA on neurodegenerative disease resulting from injury will likely offer key insights into the genetic factors that determine the outcome of these diseases together with potential targets for treatment.

MitoScape: A big-data, machine-learning platform for obtaining mitochondrial DNA from next-generation sequencing data

The growing number of next-generation sequencing (NGS) data presents a unique opportunity to study the combined impact of mitochondrial and nuclear-encoded genetic variation in complex disease. Mitochondrial DNA variants and in particular, heteroplasmic variants, are critical for determining human disease severity. While there are approaches for obtaining mitochondrial DNA variants from NGS data, these software do not account for the unique characteristics of mitochondrial genetics and can be inaccurate even for homoplasmic variants. We introduce MitoScape, a novel, big-data, software for extracting mitochondrial DNA sequences from NGS. MitoScape adopts a novel departure from other algorithms by using machine learning to model the unique characteristics of mitochondrial genetics. We also employ a novel approach of using rho-zero (mitochondrial DNA-depleted) data to model nuclear-encoded mitochondrial sequences. We showed that MitoScape produces accurate heteroplasmy estimates using gold-standard mitochondrial DNA data. We provide a comprehensive comparison of the most common tools for obtaining mtDNA variants from NGS and showed that MitoScape had superior performance to compared tools in every statistically category we compared, including false positives and false negatives. By applying MitoScape to common disease examples, we illustrate how MitoScape facilitates important heteroplasmy-disease association discoveries by expanding upon a reported association between hypertrophic cardiomyopathy and mitochondrial haplogroup T in men (adjusted p-value = 0.003). The improved accuracy of mitochondrial DNA variants produced by MitoScape will be instrumental in diagnosing disease in the context of personalized medicine and clinical diagnostics.

Recent studies have highlighted the importance of mitochondrial DNA variation in both primary mitochondrial disease and complex, human pathology including COVID-19, and space-flight stress. The vast amount of existing, next-generation sequencing (NGS) data can be leveraged to interrogate both nuclear and mitochondrial DNA (mtDNA) sequence simultaneously, allowing for analysis of the interplay between mitochondrial and nuclear encoded genes in mitochondrial function. Identifying mtDNA sequence accurately is complicated by the presence of nuclear encoded mitochondrial sequences (NUMTs), which are homologous to mtDNA. Current software for analyzing mtDNA from NGS do not accurately model the unique characteristics of mitochondrial genetics. We introduce MitoScape, a novel, big-data, software which models mitochondrial genetics through machine learning to accurately identify mtDNA sequence from NGS data. MitoScape takes advantage of rho-zero cell data to model the characteristics of NUMTs. We show that MitoScape produces more accurate heteroplasmy estimates compared to published software. We provide an example of applying MitoScape in replicating an association between hypertrophic cardiomyopathy and mitochondrial haplogroup T in men. MitoScape is an important contribution to mitochondrial genomics allowing for accurate mtDNA variants, and the ability to tailor mtDNA analysis in different population and disease contexts, which is not available in other software.

Pt-ttpy, a G-quadruplex binding platinum complex, induces telomere dysfunction and G-rich regions DNA damage

Pt-ttpy (tolyl terpyridin-Pt complex) covalently binds to G-quadruplex (G4) structures in vitro and to telomeres in cellulo via its Pt moiety. Here, we identified its targets in the human genome, in comparison to Pt-tpy, its derivative without G4 affinity, and cisplatin. Pt-ttpy, but not Pt-tpy, induces the release of the shelterin protein TRF2 from telomeres concomitantly to the formation of DNA damage foci at telomeres but also at other chromosomal locations. γ-H2AX chromatin immunoprecipitation (ChIP-seq) after treatment with Pt-ttpy or cisplatin revealed accumulation in G- and A-rich tandemly repeated sequences, but not particularly in potential G4 forming sequences. Collectively, Pt-ttpy presents dual targeting efficiency on DNA, by inducing telomere dysfunction and genomic DNA damage at specific loci.

Chromosome-length genome assembly and structural variations of the primal Basenji dog (Canis lupus familiaris) genome

BACKGROUND

Basenjis are considered an ancient dog breed of central African origins that still live and hunt with tribesmen in the African Congo. Nicknamed the barkless dog, Basenjis possess unique phylogeny, geographical origins and traits, making their genome structure of great interest. The increasing number of available canid reference genomes allows us to examine the impact the choice of reference genome makes with regard to reference genome quality and breed relatedness.

RESULTS

Here, we report two high quality de novo Basenji genome assemblies: a female, China (CanFam_Bas), and a male, Wags. We conduct pairwise comparisons and report structural variations between assembled genomes of three dog breeds: Basenji (CanFam_Bas), Boxer (CanFam3.1) and German Shepherd Dog (GSD) (CanFam_GSD). CanFam_Bas is superior to CanFam3.1 in terms of genome contiguity and comparable overall to the high quality CanFam_GSD assembly. By aligning short read data from 58 representative dog breeds to three reference genomes, we demonstrate how the choice of reference genome significantly impacts both read mapping and variant detection.

CONCLUSIONS

The growing number of high-quality canid reference genomes means the choice of reference genome is an increasingly critical decision in subsequent canid variant analyses. The basal position of the Basenji makes it suitable for variant analysis for targeted applications of specific dog breeds. However, we believe more comprehensive analyses across the entire family of canids is more suited to a pangenome approach. Collectively this work highlights the importance the choice of reference genome makes in all variation studies.

Intronic tRNAs of mitochondrial origin regulate constitutive and alternative splicing

BACKGROUND

The presence of nuclear mitochondrial DNA (numtDNA) has been reported within several nuclear genomes. Next to mitochondrial protein-coding genes, numtDNA sequences also encode for mitochondrial tRNA genes. However, the biological roles of numtDNA remain elusive.

RESULTS

Employing in silico analysis, we identify 281 mitochondrial tRNA homologs in the human genome, which we term nimtRNAs (nuclear intronic mitochondrial-derived tRNAs), being contained within introns of 76 nuclear host genes. Despite base changes in nimtRNAs when compared to their mtRNA homologs, a canonical tRNA cloverleaf structure is maintained. To address potential functions of intronic nimtRNAs, we insert them into introns of constitutive and alternative splicing reporters and demonstrate that nimtRNAs promote pre-mRNA splicing, dependent on the number and positioning of nimtRNA genes and splice site recognition efficiency. A mutational analysis reveals that the nimtRNA cloverleaf structure is required for the observed splicing increase. Utilizing a CRISPR/Cas9 approach, we show that a partial deletion of a single endogenous nimtRNALys within intron 28 of the PPFIBP1 gene decreases inclusion of the downstream-located exon 29 of the PPFIBP1 mRNA. By employing a pull-down approach followed by mass spectrometry, a 3'-splice site-associated protein network is identified, including KHDRBS1, which we show directly interacts with nimtRNATyr by an electrophoretic mobility shift assay.

CONCLUSIONS

We propose that nimtRNAs, along with associated protein factors, can act as a novel class of intronic splicing regulatory elements in the human genome by participating in the regulation of splicing.

Characterization of nuclear mitochondrial insertions in the whole genomes of primates

The transfer and integration of whole and partial mitochondrial genomes into the nuclear genomes of eukaryotes is an ongoing process that has facilitated the transfer of genes and contributed to the evolution of various cellular pathways. Many previous studies have explored the impact of these insertions, referred to as NumtS, but have focused primarily on older events that have become fixed and are therefore present in all individual genomes for a given species. We previously developed an approach to identify novel Numt polymorphisms from next-generation sequence data and applied it to thousands of human genomes. Here, we extend this analysis to 79 individuals of other great ape species including chimpanzee, bonobo, gorilla, orang-utan and also an old world monkey, macaque. We show that recent Numt insertions are prevalent in each species though at different apparent rates, with chimpanzees exhibiting a significant increase in both polymorphic and fixed Numt sequences as compared to other great apes. We further assessed positional effects in each species in terms of evolutionary time and rate of insertion and identified putative hotspots on chromosome 5 for Numt integration, providing insight into both recent polymorphic and older fixed reference NumtS in great apes in comparison to human events.

Whole-genome de novo assemblies reveal extensive structural variations and dynamic organelle-to-nucleus DNA transfers in African and Asian rice

Asian cultivated rice (Oryza sativa) and African cultivated rice (Oryza glaberrima) originated from the wild rice species Oryza rufipogon and Oryza barthii, respectively. The genomes of both cultivated species have undergone profound changes during domestication. Whole-genome de novo assemblies of O. barthii, O. glaberrima, O. rufipogon and Oryza nivara, produced using PacBio single-molecule real-time (SMRT) and next-generation sequencing (NGS) technologies, showed that Gypsy-like retrotransposons are the major contributors to genome size variation in African and Asian rice. Through the detection of genome-wide structural variations (SVs), we observed that besides 28 shared SV hot spots, another 67 hot spots existed in either the Asian or African rice genomes. Based on gene annotation information of the SVs, we established that organelle-to-nucleus DNA transfers resulted in numerous SVs that participated in the nuclear genome divergence of rice species and subspecies. We detected 52 giant nuclear integrants of organelle DNA (NORGs, defined as >10 kb) in six Oryza AA genomes. In addition, we developed an effective method to genotype giant NORGs, based on genome assembly, and first showed the dynamic change in the distribution of giant NORGs in rice natural population. Interestingly, 16 highly differentiated giant NORGs tended to accumulate in natural populations of Asian rice from higher latitude regions, grown at lower temperatures and light intensities. Our study provides new insight into the genome divergence of African and Asian rice, and establishes that organelle-to-nucleus DNA transfers, as potentially powerful contributors to environmental adaptation during rice evolution, play a major role in producing SVs in rice genomes.

Tracking the Distribution and Burst of Nuclear Mitochondrial DNA Sequences (NUMTs) in Fig Wasp Genomes

Simple Summary

Nuclear mitochondrial DNA sequences (NUMTs), which result from the insertion of exogenous mtDNA into the nuclear genome, are widely distributed in eukaryotes. However, how NUMTs are inserted into the nuclear genome and their post-insertion fates remain a mystery. Previous studies have suggested that Hymenoptera may be a group rich in NUMTs, which will be helpful to study the biological issues of NUMTs. We here select 11 species of fig wasps (Chalcidoidea, Hymenoptera) to analyze the distribution and evolution of NUMTs at the genomic level. The results show that the distributions of NUMTs are species- or lineage-specific. Furthermore, genomic environmental factors such as genome size, the damage-prone regions, and the mode of TE dynamics can determine the insertion and post-insertion fate of NUMTs. Especially because of TEs, the fragmentation and duplication of NUMTs, and thus their burst, are common. This is a relatively comprehensive investigation of the specific distribution of NUMTs and its influencing factors. Our study will help people to understand the evolution of exogenous fragments in the nuclear genome.

Abstract

Mitochondrial DNA sequences can be transferred into the nuclear genome, giving rise to nuclear mitochondrial DNA sequences (NUMTs). NUMTs have been described in numerous eukaryotes. However, the studies on the distribution of NUMTs and its influencing factors are still inadequate and even controversial. Previous studies have suggested that Hymenoptera may be a group rich in NUMTs, in which we selected 11 species of fig wasps (Chalcidoidea, Hymenoptera) to analyze the distribution and evolution of NUMTs at the genomic level. The results showed that the contents of NUMTs varied greatly in these species, and bursts of NUMTs existed in some species or lineages. Further detailed analyses showed that the large number of NUMTs might be related to the large genomes; NUMTs tended to be inserted into unstable regions of the genomes; and the inserted NUMTs might also be affected by transposable elements (TEs) in the neighbors, leading to fragmentations and duplications, followed by bursts of NUMTs. In summary, our results suggest that a variety of genomic environmental factors can determine the insertion and post-insertion fate of NUMTs, resulting in their species- or lineage-specific distribution patterns, and that studying the evolution of NUMTs can provide good evidence and theoretical basis for exploring the dynamics of exogenous DNA entering into the nuclear genome.

A pipeline for complete characterization of complex germline rearrangements from long DNA reads

BACKGROUND

Many genetic/genomic disorders are caused by genomic rearrangements. Standard methods can often characterize these variations only partly, e.g., copy number changes or breakpoints. It is important to fully understand the order and orientation of rearranged fragments, with precise breakpoints, to know the pathogenicity of the rearrangements.

METHODS

We performed whole-genome-coverage nanopore sequencing of long DNA reads from four patients with chromosomal translocations. We identified rearrangements relative to a reference human genome, subtracted rearrangements shared by any of 33 control individuals, and determined the order and orientation of rearranged fragments, with our newly developed analysis pipeline.

RESULTS

We describe the full characterization of complex chromosomal rearrangements, by filtering out genomic rearrangements seen in controls without the same disease, reducing the number of loci per patient from a few thousand to a few dozen. Breakpoint detection was very accurate; we usually see ~ 0 ± 1 base difference from Sanger sequencing-confirmed breakpoints. For one patient with two reciprocal chromosomal translocations, we find that the translocation points have complex rearrangements of multiple DNA fragments involving 5 chromosomes, which we could order and orient by an automatic algorithm, thereby fully reconstructing the rearrangement. A rearrangement is more than the sum of its parts: some properties, such as sequence loss, can be inferred only after reconstructing the whole rearrangement. In this patient, the rearrangements were evidently caused by shattering of the chromosomes into multiple fragments, which rejoined in a different order and orientation with loss of some fragments.

CONCLUSIONS

We developed an effective analytic pipeline to find chromosomal aberration in congenital diseases by filtering benign changes, only from long read sequencing. Our algorithm for reconstruction of complex rearrangements is useful to interpret rearrangements with many breakpoints, e.g., chromothripsis. Our approach promises to fully characterize many congenital germline rearrangements, provided they do not involve poorly understood loci such as centromeric repeats.

Exploring the Relationship Among Divergence Time and Coding and Non-coding Elements in the Shaping of Fungal Mitochondrial Genomes

The order Hypocreales (Ascomycota) is composed of ubiquitous and ecologically diverse fungi such as saprobes, biotrophs, and pathogens. Despite their phylogenetic relationship, these species exhibit high variability in biomolecules production, lifestyle, and fitness. The mitochondria play an important role in the fungal biology, providing energy to the cells and regulating diverse processes, such as immune response. In spite of its importance, the mechanisms that shape fungal mitogenomes are still poorly understood. Herein, we investigated the variability and evolution of mitogenomes and its relationship with the divergence time using the order Hypocreales as a study model. We sequenced and annotated for the first time Trichoderma harzianum mitochondrial genome (mtDNA), which was compared to other 34 mtDNAs species that were publicly available. Comparative analysis revealed a substantial structural and size variation on non-coding mtDNA regions, despite the conservation of copy number, length, and structure of protein-coding elements. Interestingly, we observed a highly significant correlation between mitogenome length, and the number and size of non-coding sequences in mitochondrial genome. Among the non-coding elements, group I and II introns and homing endonucleases genes (HEGs) were the main contributors to discrepancies in mitogenomes structure and length. Several intronic sequences displayed sequence similarity among species, and some of them are conserved even at gene position, and were present in the majority of mitogenomes, indicating its origin in a common ancestor. On the other hand, we also identified species-specific introns that advocate for the origin by different mechanisms. Investigation of mitochondrial gene transfer to the nuclear genome revealed that nuclear copies of the nad5 are the most frequent while atp8, atp9, and cox3 could not be identified in any of the nuclear genomes analyzed. Moreover, we also estimated the divergence time of each species and investigated its relationship with coding and non-coding elements as well as with the length of mitogenomes. Altogether, our results demonstrated that introns and HEGs are key elements on mitogenome shaping and its presence on fast-evolving mtDNAs could be mostly explained by its divergence time, although the intron sharing profile suggests the involvement of other mechanisms on the mitochondrial genome evolution, such as horizontal transference.

Survey of mitochondrial sequences integrated into the bovine nuclear genome

Nuclear copies of the mitochondrial DNA (NUMTs) have already been described in several species. In this context, we identified and analysed 166 bovine NUMT regions with a total length of 430 kbp, representing about 0.02% of the cattle nuclear genome. Copies of all mitochondrial regions were detected in the nuclear genome, with distinct degrees of sequence similarity to the mitogenome. Some NUMT regions include large mitogenome segments and show high similarity to the organelle genome sequence. NUMT regions are frequently modified by insertion of repetitive sequences and by sequence rearrangements. We confirmed the existence of 29 NUMT regions by PCR amplification using DNA from the cow (Dominette) which was used to generate the bovine genome reference sequence, ruling out the possibility that these NUMTs could be artifacts of the genome assembly. As there are NUMT regions with high similarity to the mitogenome, special care is needed when designing primers for mitochondrial DNA amplification. Our results can therefore be used to avoid co-amplification of bovine nuclear sequences similar to mitochondrial DNA.

Transposable Elements Cross Kingdom Boundaries and Contribute to Inflammation and Ageing: Somatic Acquisition of Foreign Transposable Elements as a Catalyst of Genome Instability, Epigenetic Dysregulation, Inflammation, Senescence, and Ageing

The de-repression of transposable elements (TEs) in mammalian genomes is thought to contribute to genome instability, inflammation, and ageing, yet is viewed as a cell-autonomous event. In contrast to mammalian cells, prokaryotes constantly exchange genetic material through TEs, crossing both cell and species barriers, contributing to rapid microbial evolution and diversity in complex communities such as the mammalian gut. Here, it is proposed that TEs released from prokaryotes in the microbiome or from pathogenic infections regularly cross the kingdom barrier to the somatic cells of their eukaryotic hosts. It is proposed this horizontal transfer of TEs from microbe to host is a stochastic, ongoing catalyst of genome destabilization, resulting in structural and epigenetic variations, and activation of well-evolved host defense mechanisms contributing to inflammation, senescence, and biological ageing. It is proposed that innate immunity pathways defend against the horizontal acquisition of microbial TEs, and that activation of this pathway during horizontal transposon transfer promotes chronic inflammation during ageing. Finally, it is suggested that horizontal acquisition of prokaryotic TEs into mammalian genomes has been masked and subsequently under-reported due to flaws in current sequencing pipelines, and new strategies to uncover these events are proposed.

Archaic mitochondrial DNA inserts in modern day nuclear genomes

Background

Traces of interbreeding of Neanderthals and Denisovans with modern humans in the form of archaic DNA have been detected in the genomes of present-day human populations outside sub-Saharan Africa. Up to now, only nuclear archaic DNA has been detected in modern humans; we therefore attempted to identify archaic mitochondrial DNA (mtDNA) residing in modern human nuclear genomes as nuclear inserts of mitochondrial DNA (NUMTs).

Results

We analysed 221 high-coverage genomes from Oceania and Indonesia using an approach which identifies reads that map both to the nuclear and mitochondrial DNA. We then classified reads according to the source of the mtDNA, and found one NUMT of Denisovan mtDNA origin, present in 15 analysed genomes; analysis of the flanking region suggests that this insertion is more likely to have happened in a Denisovan individual and introgressed into modern humans with the Denisovan nuclear DNA, rather than in a descendant of a Denisovan female and a modern human male.

Conclusions

Here we present our pipeline for detecting introgressed NUMTs in next generation sequencing data that can be used on genomes sequenced in the future. Further discovery of such archaic NUMTs in modern humans can be used to detect interbreeding between archaic and modern humans and can reveal new insights into the nature of such interbreeding events.

Mitochondrial DNA methylation and copy number predict body composition in a young female population

Background

Since both genomic and environmental factors are involved in obesity etiology, several studies about the influence of adiposity on both nuclear DNA and mitochondrial DNA methylation patterns have been carried out. Nevertheless, few evidences exploring the usage of buccal swab samples to study mitochondrial DNA epigenetics can be found in literature.

Methods

In this study, mitochondrial DNA from buccal swabs collected from a young Caucasian population (n = 69) have been used to examine potential correlation between mitochondrial DNA copy number and methylation with body composition (BMI, WHtR and bioimpedance measurements).

Results

A negative correlation between mitochondrial DNA copy number and BMI was measured in females (p = 0.028), but not in males. The mean percentage of D-loop methylation is significantly higher in overweight than in lean female subjects (p = 0.003), and a specific CpG located in the D-loop shows per se an association with impaired body composition (p = 0.004). Body composition impairment is predicted by a combined variable including mtDNA copy number and the D-loop methylation (AUC = 0.785; p = 0.009).

Conclusions

This study corroborates the hypothesis that mitochondrial DNA carries relevant information about body composition. However, wider investigations able to validate the usage of mtDNA methylation from buccal swabs as a biomarker are warranted.

Alignment-free approaches for predicting novel Nuclear Mitochondrial Segments (NUMTs) in the human genome

The nuclear human genome harbors sequences of mitochondrial origin, indicating an ancestral transfer of DNA from the mitogenome. Several Nuclear Mitochondrial Segments (NUMTs) have been detected by alignment-based sequence similarity search, as implemented in the Basic Local Alignment Search Tool (BLAST). Identifying NUMTs is important for the comprehensive annotation and understanding of the human genome. Here we explore the possibility of detecting NUMTs in the human genome by alignment-free sequence similarity search, such as k-mers (k-tuples, k-grams, oligos of length k) distributions. We find that when k=6 or larger, the k-mer approach and BLAST search produce almost identical results, e.g., detect the same set of NUMTs longer than 3 kb. However, when k=5 or k=4, certain signals are only detected by the alignment-free approach, and these may indicate yet unrecognized, and potentially more ancestral NUMTs. We introduce a "Manhattan plot" style representation of NUMT predictions across the genome, which are calculated based on the reciprocal of the Jensen-Shannon divergence between the nuclear and mitochondrial k-mer frequencies. The further inspection of the k-mer-based NUMT predictions however shows that most of them contain long-terminal-repeat (LTR) annotations, whereas BLAST-based NUMT predictions do not. Thus, similarity of the mitogenome to LTR sequences is recognized, which we validate by finding the mitochondrial k-mer distribution closer to those for transposable sequences and specifically, close to some types of LTR.

Genetic variation of the nuclear sequences of mitochondrial origin associated with retrotransposon Tv1 insertions in Drosophila species of the virilis group

Mitochondrial DNA sequences integrated into chromosomes are a promising object for designing genetic markers for studies of phylogenesis and genomic instability. Mitochondrial genomes of D. virilis and other Drosophila species of the virilis group contain (AT)nmicrosatellites in the spacer region between the atp6 and cox3 genes, and this microsatellite sequence is one of the hallmarks of the virilis group. The nuclear genome of D. virilis contains many extended fragments of mitochondrial DNA, which in total are several times longer than the mitochondrial genome. These nuclear sequences of mitochondrial origin contain all types of mitochondrial sequences, including mitochondrial genes and the aforementioned microsatellite sequence. The presence of the (AT)nmicrosatellite allows insertion of retrotransposon Tv1, which can transpose into the (AT)n microsatellite in a site-specific manner. The Tv1 insertion into (AT)n, close to the atp6 or cox3 pseudogenes produces a unique sequence. This sequence is formed by retrotransposon Tv1 and pseudogenes atp6 or cox3. This unique sequence can be detected in the genome by a PCR-based method. We applied this method to the detection and analysis of the nucleotide variability of the pseudogenes atp6 and cox3 associated with Tv1 insertions in a D. virilis cell culture and in the genomes of four Drosophila species of the virilis group: D. virilis, D. montana, D. borealis, and D. lacicola. We discovered new events of mitochondrial sequence transfer to the nucleus in the transplanted cell culture of D. virilis, and new Tv1 insertions, having emerged during the passage of this cell line were detected in the genome of the D. virilis transplanted cell culture. We found atp6 and cox3 pseudogenes associated with insertions of retrotransposon Tv1 in the nuclear genomes of four Drosophila species from the virilis group. These chimeric sequences proved to be species-specific. The age of the Tv1 insertion into the atp6 and cox3 pseudogenes is estimated at 1.50 Ma for D. virilis, 1.31 Ma for D. lacicola, and 1.56 Ma for D. borealis. A specific situation was revealed for D. montana, in which Tv1 insertions with nearly identical 5' and 3' long terminal repeats (LTRs) were present in accessions of flies from Europe and Asia. The age of this insertion was about 300 thousand years, and the insertion was absent from the D. montana fly line from North America.

Mitochondrial DNA and their nuclear copies in the parasitic wasp Pteromalus puparum: A comparative analysis in Chalcidoidea

Chalcidoidea (chalcidoid wasps) are an abundant and megadiverse insect group with both ecological and economical importance. Here we report a complete mitochondrial genome in Chalcidoidea from Pteromalus puparum (Pteromalidae). Eight tandem repeats followed by 6 reversed repeats were detected in its 3308 bp control region. This long and complex control region may explain failures of amplifying and sequencing of complete mitochondrial genomes in some chalcidoids. In addition to 37 typical mitochondrial genes, an extra identical isoleucine tRNA (trnI) was detected at the opposite end of the control region. This recent mitochondrial gene duplication indicates that gene arrangements in chalcidoids are ongoing. A comparison among available chalcidoid mitochondrial genomes reveals rapid gene order rearrangements overall and high protein substitution rates in most chalcidoid taxa. In addition, we identified 24 nuclear sequences of mitochondrial origin (NUMTs) in P. puparum, summing up to 9989 bp, with 3617 bp of these NUMTs originating from mitochondrial coding regions. NUMTs abundance in P. puparum is only one-twelfth of that in its relative, Nasonia vitripennis. Based on phylogenetic analysis, we provide evidence that a faster nuclear degradation rate contributes to the reduced NUMT numbers in P. puparum. Overall, our study shows unusually high rates of mitochondrial evolution and considerable variation in NUMT accumulation in Chalcidoidea.

The landscape of mitochondrial small non-coding RNAs in the PGCs of male mice, spermatogonia, gametes and in zygotes

Background

Mitochondria are organelles that fulfill a fundamental role in cell bioenergetics, as well as in other processes like cell signaling and death. Small non-coding RNAs (sncRNA) are now being considered as pivotal post-transcriptional regulators, widening the landscape of their diversity and functions. In mammalian cells, small RNAs encoded by the mitochondrial genome, mitosRNAs were discovered recently, although their biological role remains uncertain.

Results

Here, using specific bioinformatics analyses, we have defined the diversity of mitosRNAs present in early differentiated germ cells of male mice (PGCs and spermatogonia), and in the gametes of both sexes and in zygotes. We found strong transcription of mitosRNAs relative to the size of the mtDNA, and classifying these mitosRNAs into different functional sncRNA groups highlighted the predominance of Piwi-interacting RNAs (piRNAs) relative to the other types of mitosRNAs. Mito-piRNAs were more abundant in oocytes and zygotes, where mitochondria fulfill key roles in fecundation process. Functional analysis of some particular mito-piRNAs (mito-piR-7,456,245), also expressed in 3T3-L1 cells, was assessed after exposure to RNA antagonists.

Conclusions

As far as we are aware, this is the first integrated analysis of sncRNAs encoded by mtDNA in germ cells and zygotes. The data obtained suggesting that mitosRNAs fulfill key roles in gamete differentiation and fertilization.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5020-3) contains supplementary material, which is available to authorized users.

Few mitochondrial DNA sequences are inserted into the turkey (Meleagris gallopavo) nuclear genome: evolutionary analyses and informativity in the domestic lineage

Mitochondrial DNA (mtDNA) insertions have been detected in the nuclear genome of many eukaryotes. These sequences are pseudogenes originated by horizontal transfer of mtDNA fragments into the nuclear genome, producing nuclear DNA sequences of mitochondrial origin (numt). In this study we determined the frequency and distribution of mtDNA-originated pseudogenes in the turkey (Meleagris gallopavo) nuclear genome. The turkey reference genome (Turkey_2.01) was aligned with the reference linearized mtDNA sequence using last. A total of 32 numt sequences (corresponding to 18 numt regions derived by unique insertional events) were identified in the turkey nuclear genome (size ranging from 66 to 1415 bp; identity against the modern turkey mtDNA corresponding region ranging from 62% to 100%). Numts were distributed in nine chromosomes and in one scaffold. They derived from parts of 10 mtDNA protein-coding genes, ribosomal genes, the control region and 10 tRNA genes. Seven numt regions reported in the turkey genome were identified in orthologues positions in the Gallus gallus genome and therefore were present in the ancestral genome that in the Cretaceous originated the lineages of the modern crown Galliformes. Five recently integrated turkey numts were validated by PCR in 168 turkeys of six different domestic populations. None of the analysed numts were polymorphic (i.e. absence of the inserted sequence, as reported in numts of recent integration in other species), suggesting that the reticulate speciation model is not useful for explaining the origin of the domesticated turkey lineage.

A genomic landscape of mitochondrial DNA insertions in the pig nuclear genome provides evolutionary signatures of interspecies admixture

Nuclear DNA sequences of mitochondrial origin (numts) are derived by insertion of mitochondrial DNA (mtDNA), into the nuclear genome. In this study, we provide, for the first time, a genome picture of numts inserted in the pig nuclear genome. The Sus scrofa reference nuclear genome (Sscrofa10.2) was aligned with circularized and consensus mtDNA sequences using LAST software. A total of 430 numt sequences that may represent 246 different numt integration events (57 numt regions determined by at least two numt sequences and 189 singletons) were identified, covering about 0.0078% of the nuclear genome. Numt integration events were correlated (0.99) to the chromosome length. The longest numt sequence (about 11 kbp) was located on SSC2. Six numts were sequenced and PCR amplified in pigs of European commercial and local pig breeds, of the Chinese Meishan breed and in European wild boars. Three of them were polymorphic for the presence or absence of the insertion. Surprisingly, the estimated age of insertion of two of the three polymorphic numts was more ancient than that of the speciation time of the Sus scrofa, supporting that these polymorphic sites were originated from interspecies admixture that contributed to shape the pig genome.

Quantifying the Number of Independent Organelle DNA Insertions in Genome Evolution and Human Health

Fragments of organelle genomes are often found as insertions in nuclear DNA. These fragments of mitochondrial DNA (numts) and plastid DNA (nupts) are ubiquitous components of eukaryotic genomes. They are, however, often edited out during the genome assembly process, leading to systematic underestimation of their frequency. Numts and nupts, once inserted, can become further fragmented through subsequent insertion of mobile elements or other recombinational events that disrupt the continuity of the inserted sequence relative to the genuine organelle DNA copy. Because numts and nupts are typically identified through sequence comparison tools such as BLAST, disruption of insertions into smaller fragments can lead to systematic overestimation of numt and nupt frequencies. Accurate identification of numts and nupts is important, however, both for better understanding of their role during evolution, and for monitoring their increasingly evident role in human disease. Human populations are polymorphic for 141 numt loci, five numts are causal to genetic disease, and cancer genomic studies are revealing an abundance of numts associated with tumor progression. Here, we report investigation of salient parameters involved in obtaining accurate estimates of numt and nupt numbers in genome sequence data. Numts and nupts from 44 sequenced eukaryotic genomes reveal lineage-specific differences in the number, relative age and frequency of insertional events as well as lineage-specific dynamics of their postinsertional fragmentation. Our findings outline the main technical parameters influencing accurate identification and frequency estimation of numts in genomic studies pertinent to both evolution and human health.

Cytoplasmic transfer of heritable elements other than mtDNA from SAMP1 mice into mouse tumor cells suppresses their ability to form tumors in C57BL6 mice

In a previous study, we generated transmitochondrial P29mtSAMP1 cybrids, which had nuclear DNA from the C57BL6 (referred to as B6) mouse strain-derived P29 tumor cells and mitochondrial DNA (mtDNA) exogenously-transferred from the allogeneic strain SAMP1. Because P29mtSAMP1 cybrids did not form tumors in syngeneic B6 mice, we proposed that allogeneic SAMP1 mtDNA suppressed tumor formation of P29mtSAMP1 cybrids. To test this hypothesis, current study generated P29mt(sp)B6 cybrids carrying all genomes (nuclear DNA and mtDNA) from syngeneic B6 mice by eliminating SAMP1 mtDNA from P29mtSAMP1 cybrids and reintroducing B6 mtDNA. However, the P29mt(sp)B6 cybrids did not form tumors in B6 mice, even though they had no SAMP1 mtDNA, suggesting that SAMP1 mtDNA is not involved in tumor suppression. Then, we examined another possibility of whether SAMP1 mtDNA fragments potentially integrated into the nuclear DNA of P29mtSAMP1 cybrids are responsible for tumor suppression. We generated P29^H(sp)B6 cybrids by eliminating nuclear DNA from P29mt(sp)B6 cybrids and reintroducing nuclear DNA with no integrated SAMP1 mtDNA fragment from mtDNA-less P29 cells resistant to hygromycin in selection medium containing hygromycin. However, the P29^H(sp)B6 cybrids did not form tumors in B6 mice, even though they carried neither SAMP1 mtDNA nor nuclear DNA with integrated SAMP1 mtDNA fragments. Moreover, overproduction of reactive oxygen species (ROS) and bacterial infection were not involved in tumor suppression. These observations suggest that tumor suppression was caused not by mtDNA with polymorphic mutations or infection of cytozoic bacteria but by hypothetical heritable cytoplasmic elements other than mtDNA from SAMP1 mice.

Data on the time of integration of the human mitochondrial pseudogenes (NUMTs) into the nuclear genome

The data and methods presented in this article are supplementing the research article "Integration of mtDNA pseudogenes into the nuclear genome coincides with speciation of the human genus. A hypothesis", DOI: 10.1016/j.mito.2016.12.001 (Gunbin et al., 2017) [1]. Mitochondrial DNA is known to get inserted into nuclear DNA to form NUMTs, i.e. nuclear DNA pseudogenes of the mtDNA. We present here the sequences of selected NUMTs, in which time of integration can be determined with sufficient precision. We report their chromosomal positions , their position within the great ape mtDNA phylogeny, and their times of integration into the nuclear genome. The methods used to generate the data and to control their quality are also presented. The dataset is made publicly available to enable critical or extended analyzes.

Numtogenesis as a mechanism for development of cancer

Transfer of genetic material from cytoplasmic organelles to the nucleus, an ongoing process, has implications in evolution, aging, and human pathologies such as cancer. The transferred mitochondrial DNA (mtDNA) fragments in the nuclear genome are called nuclear mtDNA or NUMTs. We have named the process numtogenesis, defining the term as the transfer of mtDNA into the nuclear genome, or, less specifically, the transfer of mitochondria or mitochondrial components into the nucleus. There is increasing evidence of the involvement of NUMTs in human biology and pathology. Although information pertaining to NUMTs and numtogenesis is sparse, the role of this aspect of mitochondrial biology to human cancers is apparent. In this review, we present available knowledge about the origin and mechanisms of numtogenesis, with special emphasis on the role of NUMTs in human malignancies. We describe studies undertaken in our laboratory and in others and discuss the influence of NUMTs in tumor initiation and progression and in survival of cancer patients. We describe suppressors of numtogenesis and evolutionary conserved mechanisms underlying numtogenesis in cancer. An understanding the emerging field of numtogenesis should allow comprehension of this process in various malignancies and other diseases and, more generally, in human health.

Migration of mitochondrial DNA in the nuclear genome of colorectal adenocarcinoma

Background

Colorectal adenocarcinomas are characterized by abnormal mitochondrial DNA (mtDNA) copy number and genomic instability, but a molecular interaction between mitochondrial and nuclear genome remains unknown. Here we report the discovery of increased copies of nuclear mtDNA (NUMT) in colorectal adenocarcinomas, which supports link between mtDNA and genomic instability in the nucleus. We name this phenomenon of nuclear occurrence of mitochondrial component as numtogenesis. We provide a description of NUMT abundance and distribution in tumor versus matched blood-derived normal genomes.

Methods

Whole-genome sequence data were obtained for colon adenocarcinoma and rectum adenocarcinoma patients participating in The Cancer Genome Atlas, via the Cancer Genomics Hub, using the GeneTorrent file acquisition tool. Data were analyzed to determine NUMT proportion and distribution on a genome-wide scale. A NUMT suppressor gene was identified by comparing numtogenesis in other organisms.

Results

Our study reveals that colorectal adenocarcinoma genomes, on average, contains up to 4.2-fold more somatic NUMTs than matched normal genomes. Women colorectal tumors contained more NUMT than men. NUMT abundance in tumor predicted parallel abundance in blood. NUMT abundance positively correlated with GC content and gene density. Increased numtogenesis was observed with higher mortality. We identified YME1L1, a human homolog of yeast YME1 (yeast mitochondrial DNA escape 1) to be frequently mutated in colorectal tumors. YME1L1 was also mutated in tumors derived from other tissues. We show that inactivation of YME1L1 results in increased transfer of mtDNA in the nuclear genome.

Conclusions

Our study demonstrates increased somatic transfer of mtDNA in colorectal tumors. Our study also reveals sex-based differences in frequency of NUMT occurrence and that NUMT in blood reflects NUMT in tumors, suggesting NUMT may be used as a biomarker for tumorigenesis. We identify YME1L1 as the first NUMT suppressor gene in human and demonstrate that inactivation of YME1L1 induces migration of mtDNA to the nuclear genome. Our study reveals that numtogenesis plays an important role in the development of cancer.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-017-0420-6) contains supplementary material, which is available to authorized users.