
==== Front
BMC Genomics
BMC Genomics
BMC Genomics
1471-2164
BioMed Central London

38486136
10201
10.1186/s12864-024-10201-9
Research
Mitochondrial genome plasticity of mammalian species
Biró Bálint biro.balint@uni-mate.hu

12
Gál Zoltán 1
Fekete Zsófia 3
Klecska Eszter 4
Hoffmann Orsolya Ivett hoffmann.orsolya.ivett@uni-mate.hu

1
1 https://ror.org/01394d192 grid.129553.9 0000 0001 1015 7851 Agribiotechnology and Precision Breeding for Food Security National Laboratory, Department of Animal Biotechnology, Institute of Genetics and Biotechnology, Hungarian University of Agriculture and Life Sciences, Szent-Györgyi Albert str. 4, 2100 Gödöllő, Hungary
2 Group BM, Data Insights Team, _VOIS, Kerepesi str. 35, 1087 Budapest, Hungary
3 https://ror.org/01394d192 grid.129553.9 0000 0001 1015 7851 Department of Genetics and Genomics, Institute of Genetics and Biotechnology, Hungarian University of Agriculture and Life Sciences, Szent-Györgyi Albert str. 4, 2100 Gödöllő, Hungary
4 FamiCord Group, Krio Institute, Kelemen László str, 1026 Budapest, Hungary
14 3 2024
14 3 2024
2024
25 27827 6 2023
8 3 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-024-10201-9.

Keywords

NUMT
Mammals
Genome
Bioinformatics
Machine learning
Hungarian University of Agriculture and Life SciencesOpen access funding provided by Hungarian University of Agriculture and Life Sciences.

issue-copyright-statement© BioMed Central Ltd., part of Springer Nature 2024
==== Body
pmcIntroduction

At the beginning of the evolution of multicellular organisms, an intracellular cooperation occurred between alpha-proteobacteria and Archaea [1, 2]. One of the organelles that have formed during this cooperation is the mitochondria, which is primarily responsible for oxidative phosphorylation, however it also participates in several other intracellular processes [2]. One of the unique phenotypic characteristics of the mitochondria is that it has its own genome (mtDNA), which is considered as the most important evidence for the endosymbiotic theory [1]. During the evolution of eukaryotes, endosymbiotic gene transfer (EGT) occurred, resulting in a large amount of genetic material being transferred from the organellar genomes to the host cell’s nuclear genome [3, 4]. The process in which mitochondrial sequences transfer to the nuclear genome is referred to as NUMTogenesis [5].

The integration of certain parts of organelle genomes into other genomes was first described in the case of the maize mitochondria and chloroplasts [6]. The presence of NUMTs in animals was first proven in the genome of the domestic cat. The NUMTs found in the cat genome are unique in several ways. On one hand, the cat nuclear genome contains almost half of the corresponding mtDNA, a 7.9 kb sequence. On the other hand, this large NUMT is present in a copy number that is several tens of times [7]. Larger NUMTs that cover almost the entire mtDNA (megaNUMTs) have also been identified in human samples [8]. Approximately 140 NUMTs have been described in the human genome [9]. NUMTs are also present in the genome of the honeybee, and their number is almost ten times higher than that of NUMTs found in the human genome [10]. The genome of organelles of endosymbiotic origin is usually less than 0.05% of their independent ancestors. Therefore, a significant part of the products of the host cell’s genome must be redirected to the mitochondria [11]. The molecular driving force behind EGT is explained by Müller’s theory [12]. According to this theory, deletions will occur in a genome that is sexually isolated, i.e. without recombination (in the case of NUMTogenesis, the mtDNA itself), which contributes to the loss of genetic material of the given genome. This results in the erosion of the genome in the short term and leading to the disappearance of the genome in the long term [13–15], mtDNA is sexually isolated because it is uniparentally inherited in most eukaryotes [16]. Therefore, the mtDNA is influenced by the effect defined by Müller [17], and the operation of EGT, by directing mitochondrial-derived genetic material into the cell nucleus, rescues the mtDNA from the degrading effect of Müller’s theory [18].

Another assumption is that maintaining multiple organelles per cell and multiple genomes per organelle requires too much energy investment from the cell. Therefore, the transfer of organelle genomes through EGT to the cell nucleus, and then redirecting the gene products back to the organelles is a less energy-intensive process than if each organelle had to do it on its own [11].

Environmental factors can lead to the transfer of genetic material from mitochondria to the nucleus, referred to as NUMTogens [5]. NUMTogens are physical [19], chemical [20] and biological stresses that increase the level of mitochondrial toxicity or stress, causing damage to the mitochondrial membrane which is the first step in the NUMTogenesis process. The disruption of the integrity of the mitochondrial membrane, resulting from events such as the excessive production of reactive oxygen species (ROS), the release of cytochrome c, and mitophagy, can allow mtDNA to escape from mitochondria. While the production of reactive oxygen species is generally considered as a random occurrence, the latter two take place during gametogenesis and are therefore largely controlled [5]. On top of that, the disruption of the integrity of the mitochondrial membrane can be caused by both external factors that induce mitostress (ionizing radiation, aging, endotoxins, ROS, endonucleases) and mutations [21]. These mitostressors damage the mtDNA to a greater extent than the nuclear genomic DNA (nDNAs) when exposed to radiation, probably due to the higher efficiency of the repair mechanisms in the cell nucleus [22]. The induction of NUMTogenesis by ionizing radiation has been proven in chicken embryos and yeast cells [19, 23]. In another study, the presence of NUMTs in the brain and liver tissue of rats was investigated over time [24]. Results showed that the number of NUMTs increased in older tissues. The increased frequency of NUMTs in aging cells has also been proven in yeast [25].

When the mitochondrial membrane is damaged, the organelle receives a degradation signal which induces the phenomenon of mitophagy. Mitophagy is a special case of autophagy that takes place in the mitophagosome [26]. This organelle is responsible for the degradation of damaged mitochondria and the recycling of its components [27]. In case of inadequate mitophagy, mitochondrial-originated sequences enter the cytoplasm. According to other theories, the entry of these sequences into the cytoplasm mainly occurs due to inadequate division and membrane fusion events [20, 27]. NUMT precursors are protected in the cytoplasm from the digestion of nucleases by a vesicle-mediated route or through the formation of a complex with histone-like proteins that bind to DNA. Mitochondrial-originated sequences located in the cytoplasm enter the nucleus through membrane fusion and/or pores [20]. After entering the nucleus, NUMT precursors can be incorporated into the nuclear genome during Double Stranded Breaks (DSB) through the non-homologous end joining (NHEJ) DNA repair mechanism. From the moment of incorporation, these sequences can be referred to as NUMTs. Mitostressors also increase the frequency of DSBs in nDNA [22]. In the case of NHEJ, a nuclease-mediated deletion always occurs in the absence of template DNA. NHEJ often results in a long single-stranded DNA segment, which increases the risk of longer deletions and translocations. According to some explanations, the cell uses NUMT precursors as template DNA to prevent larger damage caused by deletions and translocations [27].

Here we performed a systematic investigation of NUMTs across all the NCBI mammalian genomes.

We have performed this investigation in order to test the hypothesis whether all mammalian species indeed harbour NUMTs integrated within their nuclear genomes.

Based on this hypothesis our goal was also to explore the characteristics of the described NUMTs.

Furthermore, we aimed to identify the factors influencing NUMTogenesis. The motivation behind that particular hypothesis testing was to generalize our knowledge about this biological phenomenon. Previous studies have described NUMTogenesis in an isolated form with small sample sizes and different approaches, our investigation is intended to contribute to the field by providing comparable insights.

Our other hypothesis was that it is possible to predict NUMTs by utilizing modern machine learning based methods. This part is solely motivated by its practicality given that having NUMTs in an experimental setup (mainly during genome assembly and genome modification) often leads to biased, unreliable results [28–30]. Furthermore, by utilizing predictive techniques, one can eliminate the necessity of a complete reference genome when it comes to NUMT detection, hence there is no need to perform sequence alignment which is considered a computationally intensive process.

Results

Repetitive elements and NUMTs

Several repetitive elements were found that have different frequency in the 200 bp flanking regions of NUMTs than it should be expected. Three distinct behaviours could be observed when it comes to repetitive elements.

Simple repeat, Long Terminal Repeat Endogenous RetroVirus-related eLement (LTR ERVL) and low complexity elements are getting more frequent in close proximity of the upstream flanking region of NUMTs. However, the frequency of those elements is going down in distant downstream flanking regions. In the above cases, the repetitive elements have the highest density in NUMTs (Fig. 1).

Fig. 1 The frequency patterns of repetitive elements in the 200 bp flanking regions of NUMTs. Blue shaded areas are +-20 bps from NUMTs. Asterisk indicates p value < 0.05

While as the second particular behaviour, the frequencies of Short Interspersed Nuclear Element Mammalian-wide Interspersed Repeats (SINE MIR), Long Interspersed Nuclear Element L2 (LINE L2), DNA hAT-Charlie, LTR ERV1, tRNA, LINE L1, LTR ERVL-MaLR, SINE Alu and SINE B4 elements drop drastically in close upstream proximity of NUMTs, the frequency of these elements are getting sparser in distant downstream flanking of NUMTs. These elements also have their corresponding peak frequency in NUMTs. The distributions of the repetitive elements mentioned above significantly differ from random distribution (p value < 0.05).

In the third behaviour, there is no pattern in the case of sparse repetitive elements (DNA hAT-Tip100, LINE CR1, LTR ERVK, SINE B2), hence the distribution of these elements does not differ from random distributions (p value < 0.05) (Fig. 1). However, several repetitive elements did not show a particular pattern in the close proximity of NUMTs hence only 4 of them are displayed in the last row of Fig. 1.

Descriptive analysis of NUMTs

Most of the analysed nuclear genomes were between 2000 and 3000 Mb in size and had 0.4–0.42 GC ratio. We found a moderated positive Spearman correlation between nuclear genome sizes and number/cumulative length of NUMTs (0.38 and 0.33 correlation coefficients respectively with p values < 0.01) (Fig. 2/a). This means that the larger the nuclear genome the more NUMTs will be inserted.

By contrast, we found that the GC content of the nuclear genome has an opposite effect. Namely, the higher the GC content, the smaller the number and cumulative length of NUMTs (-0.29 and − 0.43 Spearman correlation coefficient respectively with p values < 0.01) (Fig. 2/b).

NUMTs tend to have lower GC content than their host nuclear genome (lower mean than 1.0). Meanwhile their corresponding mitochondrial sequences have a 1.0 centered distribution (mean ∼ 1.0) compared to their parent mtDNA genome when it comes to GC content. However, a subpopulation of NUMTs with around 1.25 relative GC content can be observed (Fig. 2/c). Since the GC content of the mtDNA is usually lower than nDNA’s GC content, this result can be considered as a feature that strengthens the theory of the mitochondrial origin of NUMTs.

The distribution of relative NUMT sizes is highly skewed toward smaller NUMTs (Fig. 2/d). The absolute NUMT sizes can be described with a 121, 632 interquartile range and a 248 median. We found nine species (Delphinapterus leucas, Tursiops truncatus, Castor canadensis, Cavia porcellus, Globicephala melas, Lagenorhynchus obliquidens, Monodelphis domestica, Orcinus orca, Theropithecus gelada and Tursiops truncatus) with their total mtDNA genomes integrated into their nuclear genomes as NUMTs. In the Castor canadensis nuclear genome, two concatenated mtDNA genomes are present as a megaNUMT.

Fig. 2 Descriptive statistics of the NCBI mammalian NUMTs. (a) Correlation between NUMT length/count and nuclear genome size in Mb. (b) Correlation between NUMT length/count and total nuclear genome GC content (%). (c) The relative GC contents of NUMTs and their corresponding mitochondrial sequences compared to the GC content of total nuclear genome and total mtDNA genome GC contents respectively. (d) Distributions of relative and absolute NUMT sizes. Relativization was performed based on the corresponding mtDNA genome size

Order specific patterns of NUMTogenesis

The NUMTs of orders that have several genomes involved, are clustered together based on NAC, kmers, NMBroto, Z-curve and mismatch profile iFeatureOmegaCLI features. For example, Artiodactyla related NUMTs which are from 25 genomes are tightly clustered together. The same applies to Rodentia and Primates. By contrast, Didelphimorphia related NUMTs that are derived from only two genomes, are scatteredand do not form such distinct clusters (Fig. 3/a, b).

Primates is the only order that is clustered tightly together when it comes to the number of nucleotide involvement in the process of NUMTogenesis (Fig. 3/c).

In general, most of the mtDNA genome is involved in NUMTogenesis throughout mammalian species. However, position specific discrepancies can be observed. Namely that the ends of the mtDNA genomes show huge differences in terms of the NUMTogenesis involvement of particular nucleotides (Fig. 3/c).

Fig. 3 NUMT features integrated with taxonomical data. (a) UMAP dimension reduction of NUMT features coloured by taxonomical order. (b) Consensus phylogenetic tree with the General Time Reversible (GTR) model and 100 bootstrap iterations of the corresponding mammalian mtDNA genomes. (c) Heatmap shows how many times a given nucleotide contributed to NUMTogenesis. Columns represent species while rows represent mitochondrial nucleotides

Machine learning approach for NUMT classification

As expected, NUMTs and random sequences form distinctly separated UMAP clusters based on the previously mentioned features despite the very few numbers of data points closer to the opposite cluster. However, NUMTs are divided into two subclusters (Fig. 4/a).

Nevertheless, no distinct pattern could be detected between these subclusters in terms of taxonomical order i.e. order specific UMAP points did not form distinct clusters.

Random search hyperparameter optimisation of a random forest model resulted in a maximum AUROC of 0.94. During the classification, 7 decision trees with a maximum depth of 3 and 50 features gave the best result (Fig. 4/b).

k-fold (k = 10) cross validation was run with the best hyperparameters (Fig. 4/c).

In every split, AUROC was calculated which resulted in a mean AUROC of 0.95 with a +- 0.01 standard deviation on the test dataset. The accuracy of the random forest classifier was 0.878 in the test set. The model correctly predicted 83% and 89% of NUMTs (17 804 correctly predicted instances) and random sequences (20 133 correctly predicted instances) of the test set respectively.

NMBroto proved to be the most important feature, while NAC and Z-curve features were the least important features in the classification of NUMTs and random sequences of the test set. Besides NMBroto, Kmer type1, RCKmer type1 and mismatch profile features were also important (Fig. 4/d).

Fig. 4 NUMT and random sequence classification with random forest algorithm. (a) UMAP dimension reduction of NUMT and random sequence features. (b) Parallel coordinates diagram of the hyperparameter optimisation. (c) k-fold cross validation ROC analysis. (d) Feature importances of iFeatureOmegaCLI’s features

Machine learning approach for flanking sequence classification

Surprisingly, the flanking regions of NUMTs and random sequences also form separated UMAP clusters. On top of that, using UMAP, the flanking positions (whether it is an upstream or downstream flanking of a random sequence etc.) can also be separated. However, some small upstream flanking clusters of random sequences are closer to the clusters of downstream flanking of random sequences (Fig. 5/a).

In the case of flanking sequence classification, the random forest algorithm’s performance was worse than the observation made during the classification of NUMTs and random sequences. However, the performance is still good, since the random search hyperparameter optimisation resulted in a maximum AUROC of 0.85 (Fig. 5/b).

Throughout the classification, 7 decision trees with a maximum depth of 3 and 50 features gave the best result (Fig. 5/b).

During k-fold (k = 10) cross validation, the mean AUROC was 0.86 with +- 0.01 standard deviation on the test dataset (Fig. 5/c).

The accuracy of the random forest classifier was about 0.774 in the test set. The model correctly predicted 68% and 73% of NUMTs’ flankings (29 374 correctly predicted instances) and random sequences’ flankings (30 373 correctly predicted instances) of the test set respectively.

It turned out that NAC and Z-curve features do not contribute to the classification, hence these two features had 0.0 feature importance. However, there is no such positive outlier importance of a particular feature as we saw in the case of NUMT and random sequence classification (Fig. 5/d).

Fig. 5 NUMT and random sequence flanking classification with random forest algorithm. (a) UMAP dimension reduction of NUMT and random flanking sequence features. (b) k-fold cross validation ROC analysis. (c) Parallel coordinates diagram of the hyperparameter optimisation. (d) Feature importance of iFeatureOmegaCLI’s features

Discussion

We observed repetitive elements that have altered presence next to NUMTs (Fig. 1).

The uneven distribution of certain repetitive elements in the close proximity of NUMTs have already been explored in several mammalian species. For example, it has been reported that in the case of several bat nuclear DNAs, the frequency of repetitive elements around NUMTs differ (higher density in the proximity of NUMT, lower density in the proximity of NUMT or no patterns at all) from other ‘non-NUMT’ regions of the genome [31]. In that study the authors came up with a theory for the co-occurrence of NUMTs and repetitive elements. Their explanation is that the repetitive elements that are enriched next to NUMTs are mainly mobile DNA related elements and so those elements can be responsible for NUMT propagation due to a copy-paste mechanism [31]. In the bovine genome, SINE and LINE elements proved to be the most frequent repetitive elements in NUMT regions and in NUMTs themselves too. These repetitive element integration events can contribute to NUMT evolution [32]. SINEs and simple repeats are also enriched in the NUMTs of the human genome [33]. However, based on this study, repetitive element integrations into NUMTs are extremely rare events (0.1% of the 66 000 genomes investigated) which mainly occur in tumour specific NUMTs. And so, the germline transmission of repetitive elements containing NUMTs is highly unlikely hence the repetitive elements cannot contribute to the evolution of NUMTs. Our result show that it is quite common that repetitive elements have their corresponding highest frequencies inside NUMTs. However, we did not investigate the tumour-specificity of these NUMTs. From these repetitive elements SINE MIR, SINE Alu and SINE B4 are non-autonomous, Class I transposable elements (TEs). One of their structural characteristics is a recognition site for a LINE mediated retrotransposition [34, 35]. As to our current knowledge, there is no evidence of SINE mediated NUMTogenesis. However, SINE elements (especially Alu) were reported as contributors to mitochondrial ROS generation and transition pore opening during ageing which are prerequisites of NUMT integration [36, 37]. Furthermore, Alu elements are frequently enriched within genes that are associated with mitochondrial transport processes [38] which supports our findings (Fig. 1).

From the LINE repetitive element class, LINE L1 repetitive element has its peak frequency in NUMTs. LINE L1 is the only autonomously active retrotransposon in the human genome. The components of LINE L1 code for open reading frame 1 and 2 proteins (ORF1p and ORF2p respectively). ORF1p is associated with the L1 RNA to form a chaperon, while ORF2p has retrotransposon related endonuclease and reverse transcriptase activities [39]. Mitochondrial inner and outer membrane translocases TIMM13 and TOMM40 do interact with ORF2p, while the mitochondrial GTPase interacts with ORF1p [40]. Additionally, it has also been described that at DSB sites, LINE L1 retrotransposition occurs more frequently [39] just as it has also been previously proved in case of NUMTs [41]. LTRs are retrotransposons that contain protein coding regions in between two long terminal repeat ends. One of their superfamilies is the superfamily of endogenous retroviruses (ERV) [42]. ERV derived transcripts are strongly connected to highly complex processes which can facilitate mitochondrial membrane permeabilization, hence can facilitate the escape of mitochondrial content into the cytosol [43]. Furthermore, cancerous human cells tend to accumulate ERV transcripts [44].

Based on our results, there is a moderate positive correlation between the size of the nuclear genome and number/total length of NUMTs (Fig. 2/a).

This means that to some extent, the bigger the nuclear genome, the more NUMTs are integrated. The same pattern was found in the case of eukaryotic genomes [27, 45]. This phenomenon is possibly due to the elevated number of DSBs in bigger genomes [27]. However, this motif is non general since in some scenarios these two features (genome size and NUMTs) are just weakly correlated [29], or not correlated at all [46]. We found a modest negative correlation between total nuclear genome GC content and number/total length of NUMTs (Fig. 2/b).

This indicates that the smaller the GC content of a genome, the more NUMTs are integrated into it. However, on the other hand, it is worth noting that the previously described results regarding the correlations were conducted on different taxonomical groups with different methods. It has been previously reported, that during certain types of cancerous transformation, NUMTs tend to integrate into gene rich regions with elevated GC contents [21]. This means that there is going to be a negative selection against NUMT integrations into genomic parts with high GC contents since those NUMTs can disrupt gene functions and cause cancer [33]. The sizes of NUMTs proved to be skewed towards shorter sequences with very few longer ones (Fig. 2/d).

Our results are consistent with the ones published before regarding the human genome when it comes to size distribution of NUMTs [33]. However, megaNUMTs, which we have described in the cases of Delphinapterus leucas, Tursiops truncatus, Castor canadensis, Cavia porcellus, Globicephala melas, Lagenorhynchus obliquidens, Monodelphis domestica, Orcinus orca, Theropithecus gelada and Tursiops truncatus have been already reported beforehand [7, 8, 47].

NUMTs form order specific clusters based on the extracted features. The taxonomically characteristic patterns of NUMTs make these sequences applicable in phylogenetic studies [48–50] (Fig. 3/a).

Despite the taxonomically well-defined inner structure of NUMT sequences, the process of NUMTogenesis itself seems to be largely universal across the species that were investigated in this study. Meaning that NUMTs are from the whole mtDNA genome, however local ‘hotspots’ (subsequences that contributed to NUMTogenesis more frequently) do exist (Fig. 3/c). Interestingly, we found out that from the dataset investigated in this study, Primates was the only order that displayed a well-defined cluster when it came to nucleotide involvement in NUMTogenesis. A possible explanation to this phenomenon is the conserved regions of the mtDNA genomes across Primates [51–53]. However, mtDNA is considered to be conserved, environmental factors seem to be highly influential to the structure and plasticity of it. For example, energy need and behaviour have high impact on the selective exposure of different parts of mtDNA [54]. Hence, mammals with bigger brain sizes tend to have different mtDNA compared to others. This difference is due to altered nucleotide substitution rates that makes catalytic activity modifiable through mitochondrial coded enzymes [55]. In our dataset, diverse spectra of mammals were included. In this case, the word ‘diverse’ can be interpreted on many levels such as body composition, energy need or even habitat [56]. As a speculative explanation, this diversity could be considered as a causing factor for the heterogeneities through the organisms investigated here.

NUMTs and random sequences form separated UMAP clusters (Fig. 4/a).

The observation is not surprising considering the fact that NUMTs are from the mtDNA genome which operates with a different genetic code compared to the nuclear genome [57]. The random forest classifier effectively distinguishes between NUMTs and random sequences, resulting in a 0.95 mean AUROC value during k-fold cross validation on the test set (Fig. 4/a, b, c).

NMBroto features had the highest feature importances, i.e. NMBroto contributed to the classification to the largest extent (Fig. 4/d).

This feature is a kind of a normalized autocorrelation measure, which symbolizes the statistical pattern of a biological sequence [58]. In other words, this feature tells us whether a property of a given nucleotide is independent from the same property of the neighbouring nucleotides [59]. NAC and Z-curve features do not contribute to the classification of NUMTs and random sequences. Both features reflect the composition of a sequence, which makes this result more interesting considering the previously mentioned fact that nuclear and mitochondrial DNA (hence NUMTs also) use different genetic codes.

The flanking regions of NUMTs and random sequences also formed well distinguished clusters. Even more remarkable is the result of the upstream and downstream flanking regions of NUMTs and random sequences clustered separately (Fig. 5/a).

It is also possible to classify the flanking regions (whether the given flanking region belongs to a NUMT or a random sequence) using a random forest-based machine learning approach (Fig. 4/c).

There is no such positive outlier in the feature importances as we have previously seen in the case of NUMT and random sequence classification. Moreover, NAC and Z-curve features did not contribute to the classification in this task either (Fig. 4/d).

It is worth mentioning that the NUMT classification algorithm that we have used in this experiment has its drawbacks. This can be seen when it comes to flanking region classification (considering the lower AUROC values). However, the random forest algorithm is still considered a reliable and widely used method in biological studies [60–63], mainly due it’s results interpretability and “non-black-box” nature.

Conclusions

In this article we characterised the NUMTs of all the available mammalian genomes at the NCBI database for the first time as to our current knowledge. We described several repetitive elements that show altered presence in the proximity of NUMTs. With the use of different features of nucleic acid sequences, we were able to classify NUMTs and also their corresponding flanking regions. Our results on the large dataset of mammalian species contribute to the theory that NUMT insertion is non-random.

To achieve our goal, we used machine learning methods that have not yet appeared in the NUMT literature hence can be considered as cutting-edge technology. Due to their ability to analyse large amounts of biological data and make accurate decisions and predictions without being explicitly programmed, machine learning methods become increasingly popular in bioinformatics. With machine learning methods, one can identify patterns and relationships in complex and high-dimensional data that may be missed with traditional statistical methods. Furthermore, NUMT prediction also has a practical benefit since NUMT integration can lead to bias during molecular works or even non accurate genome assembly. With the help of our predictive model, NUMTs can be accurately identified to eliminate their impact on the downstream analyses. Another important result of our experiments is that we have built and published here a platform that can be used by anyone, with the help of which cross-species NUMT characterization projects will finally become comparable and more manageable. This workflow can be used easily and quickly on large data, which greatly facilitates the mining of NUMTs in all genomes.

Integration of NUMTs into the nuclear genome de novo may play important roles in the development of various diseases and ageing. The increased NUMT integration with elevated temperature or increased ROS exposure, raises the question of how the amount of NUMTs in genomes will change due to climate change and the rising toxin levels found in food. A thorough understanding of the process of NUMTogenesis would give us the potential to provide greater insight into the biological relevance of the role of NUMTs in ageing, cancer, and genome integrity. By elucidating the exact mechanisms underlying NUMTogenesis, we can better understand and interpret the role of NUMTs in the genomes. This study provides a valuable contribution in two aspects. Firstly, it presents an analysis conducted on mammalian species, shedding light on the NUMTogenesis processes occurring in their genomes. Secondly, it introduces a novel workflow that facilitates the comparative analysis of these processes across diverse genomes.

Methods

Genome curation

A total number of 153 nuclear and mtDNA genomes were analysed. All the reference nuclear and mtDNA genome sequences were acquired from the FTP site of NCBI. Taxonomical related data is also from NCBI data source. The exact URLs are in the Availability of data and materials section.

Sequence alignments

The sequence alignment of the corresponding genomes was performed using LASTAL (v.: v1219) [64] with the settings and e value threshold as it has been described beforehand [65]. We chose LASTAL since it is believed to provide more accurate results for e value calculations than other widely used methods [65]. The sequence alignments of the 153 genomes resulted in 79 645 NUMTs.

While the multiple sequence alignment of the mitochondrial sequences was performed with ClustalO (v.: 1.2.4) [66] with default parameters.

Repetitive elements were investigated with RepeatMasker (v.: 4.1.2.p1) [67] using species specific settings.

Feature extraction

Nucleic acid composition (NAC), kmers (Kmer type 1 and 2, RCKmer type 1 and 2), normalized Moreau-Broto autocorrelation (NMBroto), Z-curve (geometrical features of a nucleic acid sequence [68]) and mismatch profile features were extracted with iFeatureOmegaCLI (v.: 1.0.2) [58] for further classification. IFeatureOmegaCLI provides a high-throughput, robust, easy to automate workflow for extracting meaningful features from biological sequences.

Machine learning

Uniform Manifold Approximation and Projection (UMAP), classification, grid search, random search and k-fold (k = 10) cross validation algorithms were implemented in Scikit-learn (v.: 1.0.2) [69]. During grid and random search, the corresponding hyper parameters were optimised. The best hyper parameter combination was selected based on the area under receiver operating characteristic curve (AUROC score). We decided to use the Scikit-learn implementation of the above-mentioned algorithms since these models are scalable and the resulting models are easy to share cross platforms.

Negative labelled sequence identification

For negative labelled samples acquisition, we generated as many random positions for a given genomic id as many NUMTs were located on the given genomic part. In this way we downsampled the majority class (“negative” sample sequences at random positions). This resulted in the exact same number of random sequences as many NUMTs we explored for the sake of a balanced dataset for the classification. Then, a sequence was extracted using Samtools (v.: 1.6) [70] starting from the previously determined random position in a length of the corresponding NUMTs that were integrated into the same genomic part. For instance, if the X chromosome of the mouse genome contains 2 NUMTs in a length of 123 and 304 bps, then this chromosome is going to be sampled two times at random positions in the same lengths i.e., 123 and 304 bps, respectively. As we have pointed out previously, this sampling technique provides a totally balanced dataset which is crucial when it comes to training ML models without overfitting or other biased behaviour.

Phylogenetic tree construction

Ape (v.: 5.6-2) [71] and phangorn (v.: 2.10.0) [72] libraries were used during phylogenetic analysis. The consensus phylogenetic tree was constructed from a bootstrap dataset with 100 bootstrap iterations using the GTR nucleic acid substitution model for optimizing the JC69 model.

General statistics

The statistical analysis (Spearman correlation, Kolmogorov-Smirnov test) was carried out in Scipy (v.: 1.7.3) [73]. Statistics with p value < 0.05 were considered as significant results.

Relative GC contents were calculated using whole genome GC contents derived from NCBI. The exact URL is in the Availability of data and materials section.

The NUMT mining application and the trained model are available at: https://github.com/balintbiro/NUMT_finder.

The workflow of our experiments can be seen on Fig. 6.

Fig. 6 Graphical summary of the methods

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Acknowledgements

The authors’ sincere thanks go to Endre Barta for providing the computing cluster.

Author contributions

OIH and BB conceptualized the study. All authors conceived and designed the final study. OIH coordinated the research. BB wrote the bioinformatics pipeline, ZG, EK, ZSF contributed to the analysis of the results. OIH, ZG and BB drafted the manuscript. All authors read and approved the final version of the manuscript.

Funding

BB received the NTP-NFTÖ-21-B-0127 grant. OIH received NKFIH OTKA 11964 and NKFIH OTKA 124708 grants and was supported by the National Research, Development, and Innovation Office through the Agribiotechnology and Precision Breeding for Food Security National Laboratory grant under the reference number RRF-2.3.1-21-2022-00007.

Open access funding provided by Hungarian University of Agriculture and Life Sciences.

Data availability

All the nuclear and mtDNA genome sequences are available on the FTP site of NCBI (https://ftp.ncbi.nlm.nih.gov/genomes/refseq/vertebrate_mammalian/; https://ftp.ncbi.nlm.nih.gov/genomes/refseq/mitochondrion/). Taxonomical related data is also from NCBI data source (https://ftp.ncbi.nlm.nih.gov/genomes/GENOME_REPORTS/eukaryotes.txt). The NUMT mining application and custom BASH, R and Python codes are available at GitHub repositories (https://github.com/balintbiro/NUMT_finder; https://github.com/balintbiro/numt_mining). Supplementary Information contains a file with the species names involved in the analysis.

Declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Abbreviations

AUROC Area Under the Receiver Operating characteristic Curve

DNA hAT-Charlie hAT-Charlie DNA transposon

DSB Double Stranded DNA Break

EGT Endosymbiotic Gene Transfer

ERVL Endogenous RetroVirus-related eLement

GC Guanine-Cytosine ratio

GTR General Time Reversal nucleic acid substitution model

nDNA Nuclear genomic DNA

LINE Long Interspersed Nuclear Element

LTR Long Terminal Repeat

Mb Mega base

NCBI National Centre for Biotechnology Information

NAC Nucleic Acid Composition

NHEJ Non-Homologous End Joining DNA repair mechanism

NUMT Nuclear Mitochondrial Sequence

NUMTogenesis The process in which NUMTs are created

mtDNA Mitochondrial DNA

ROS Reactive Oxygen Species

SINE Short Interspersed Nuclear Element

TRNA Transfer RNA

UMAP Uniform Manifold Approximation and Projection for dimension reduction

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Martin WF Garg S Zimorski V Endosymbiotic theories for eukaryote origin Philos Trans R Soc B Biol Sci 2015 370 1678 20140330 10.1098/rstb.2014.0330
Martin WF, Garg S, Zimorski V. Endosymbiotic theories for eukaryote origin. Philos Trans R Soc B Biol Sci. 2015;370(1678):20140330.10.1098/rstb.2014.0330
2. Roger AJ Muñoz-Gómez SA Kamikawa R The origin and diversification of mitochondria Curr Biol 2017 27 21 R1177 92 10.1016/j.cub.2017.09.015 29112874
Roger AJ, Muñoz-Gómez SA, Kamikawa R. The origin and diversification of mitochondria. Curr Biol. 2017;27(21):R1177–92.29112874 10.1016/j.cub.2017.09.015
3. Kelly S. The economics of endosymbiotic gene transfer and the evolution of organellar genomes. bioRxiv. 2020.
4. Lopez JV, Yuhki N, Masuda R, Modi W, O’Brien SJ. Numt, a recent transfer and tandem amplification of mitochondrial DNA to the nuclear genome of the domestic cat. J Mol Evol. 1994 [cited 2023 Aug 28];39(2):174–90. Available from: https://link.springer.com/article/10.1007/BF00163806.
5. Singh KK, Choudhury AR, Tiwari HK. Numtogenesis as a mechanism for development of cancer. In: Seminars in cancer biology. 2017. p. 101–9.
6. Stern DB Lonsdale DM Mitochondrial and chloroplast genomes of maize have a 12-kilobase DNA sequence in common Nature 1982 299 5885 698 702 10.1038/299698a0 6889685
Stern DB, Lonsdale DM. Mitochondrial and chloroplast genomes of maize have a 12-kilobase DNA sequence in common. Nature. 1982;299(5885):698–702.6889685 10.1038/299698a0
7. Lopez JV Cevario S O’Brien SJ Complete nucleotide sequences of the domestic cat (Felis catus) mitochondrial genome and a transposed mtDNA tandem repeat (numt) in the nuclear genome Genomics 1996 33 2 229 46 10.1006/geno.1996.0188 8660972
Lopez JV, Cevario S, O’Brien SJ. Complete nucleotide sequences of the domestic cat (Felis catus) mitochondrial genome and a transposed mtDNA tandem repeat (numt) in the nuclear genome. Genomics. 1996;33(2):229–46.8660972 10.1006/geno.1996.0188
8. Lutz-Bonengel S Niederstätter H Naue J Koziel R Yang F Sänger T Evidence for multi-copy Mega-NUMT s in the human genome Nucleic Acids Res 2021 49 3 1517 31 10.1093/nar/gkaa1271 33450006
Lutz-Bonengel S, Niederstätter H, Naue J, Koziel R, Yang F, Sänger T, et al. Evidence for multi-copy Mega-NUMT s in the human genome. Nucleic Acids Res. 2021;49(3):1517–31.33450006 10.1093/nar/gkaa1271
9. Dayama G Emery SB Kidd JM Mills RE The genomic landscape of polymorphic human nuclear mitochondrial insertions Nucleic Acids Res 2014 42 20 12640 9 10.1093/nar/gku1038 25348406
Dayama G, Emery SB, Kidd JM, Mills RE. The genomic landscape of polymorphic human nuclear mitochondrial insertions. Nucleic Acids Res. 2014;42(20):12640–9.25348406 10.1093/nar/gku1038
10. Pamilo P Viljakainen L Vihavainen A Exceptionally high density of NUMTs in the honeybee genome Mol Biol Evol 2007 24 6 1340 6 10.1093/molbev/msm055 17383971
Pamilo P, Viljakainen L, Vihavainen A. Exceptionally high density of NUMTs in the honeybee genome. Mol Biol Evol. 2007;24(6):1340–6.17383971 10.1093/molbev/msm055
11. Kelly S The economics of organellar gene loss and endosymbiotic gene transfer Genome Biol 2021 22 1 1 22 10.1186/s13059-021-02567-w 33397451
Kelly S. The economics of organellar gene loss and endosymbiotic gene transfer. Genome Biol. 2021;22(1):1–22.33397451 10.1186/s13059-021-02567-w
12. Muller HJ The relation of recombination to mutational advance Mutat Res Mol Mech Mutagen 1964 1 1 2 9 10.1016/0027-5107(64)90047-8
Muller HJ. The relation of recombination to mutational advance. Mutat Res Mol Mech Mutagen. 1964;1(1):2–9.10.1016/0027-5107(64)90047-8
13. Metzger JJ Eule S Distribution of the fittest individuals and the rate of Muller’s ratchet in a model with overlapping generations PLoS Comput Biol 2013 9 11 e1003303 10.1371/journal.pcbi.1003303 24244123
Metzger JJ, Eule S. Distribution of the fittest individuals and the rate of Muller’s ratchet in a model with overlapping generations. PLoS Comput Biol. 2013;9(11):e1003303.24244123 10.1371/journal.pcbi.1003303
14. Naito M Pawlowska TE Defying Muller’s Ratchet: ancient heritable endobacteria escape extinction through retention of recombination and genome plasticity MBio 2016 7 3 e02057 15 10.1128/mBio.02057-15 27329757
Naito M, Pawlowska TE. Defying Muller’s Ratchet: ancient heritable endobacteria escape extinction through retention of recombination and genome plasticity. MBio. 2016;7(3):e02057–15.27329757 10.1128/mBio.02057-15
15. Singh LN, Kao SH, Wallace DC. Unlocking the complexity of mitochondrial dna: a key to understanding neurodegenerative disease caused by injury. Cells. 2021;10(12).
16. Breton S Stewart DT Atypical mitochondrial inheritance patterns in eukaryotes Genome 2015 58 10 423 31 10.1139/gen-2015-0090 26501689
Breton S, Stewart DT. Atypical mitochondrial inheritance patterns in eukaryotes. Genome. 2015;58(10):423–31.26501689 10.1139/gen-2015-0090
17. Howe DK Denver DR Muller’s Ratchet and compensatory mutation in Caenorhabditis briggsae mitochondrial genome evolution BMC Evol Biol 2008 8 1 1 13 10.1186/1471-2148-8-62 18179683
Howe DK, Denver DR. Muller’s Ratchet and compensatory mutation in Caenorhabditis briggsae mitochondrial genome evolution. BMC Evol Biol. 2008;8(1):1–13.18179683 10.1186/1471-2148-8-62
18. Martin W Herrmann RG Gene transfer from organelles to the nucleus: how much, what happens, and why? Plant Physiol 1998 118 1 9 17 10.1104/pp.118.1.9 9733521
Martin W, Herrmann RG. Gene transfer from organelles to the nucleus: how much, what happens, and why? Plant Physiol. 1998;118(1):9–17.9733521 10.1104/pp.118.1.9
19. Abdullaev SA Fomenko LA Kuznetsova EA Gaziev AI Experimental detection of integration of mTDNA in the nuclear genome induced by ionizing radiation Radiatsionnaia Biol Radioecol 2013 53 4 380 8
Abdullaev SA, Fomenko LA, Kuznetsova EA, Gaziev AI. Experimental detection of integration of mTDNA in the nuclear genome induced by ionizing radiation. Radiatsionnaia Biol Radioecol. 2013;53(4):380–8.
20. Puertas MJ, González-Sánchez M. Insertions of mitochondrial DNA into the nucleus-effects and role in cell evolution. Genome. 2020 [cited 2023 Aug 28];63(8):365–74. Available from: https://pubmed.ncbi.nlm.nih.gov/32396758/.
21. Srinivasainagendra V Sandel MW Singh B Sundaresan A Mooga VP Bajpai P Migration of mitochondrial DNA in the nuclear genome of colorectal adenocarcinoma Genome Med 2017 9 1 1 15 10.1186/s13073-017-0420-6 28081715
Srinivasainagendra V, Sandel MW, Singh B, Sundaresan A, Mooga VP, Bajpai P, et al. Migration of mitochondrial DNA in the nuclear genome of colorectal adenocarcinoma. Genome Med. 2017;9(1):1–15.28081715 10.1186/s13073-017-0420-6
22. Gaziev AI Shaikhaev GO Ionizing radiation can activate the insertion of mitochondrial DNA fragments in the nuclear genome Radiatsionnaia Biol Radioecol 2007 47 6 673 83
Gaziev AI, Shaikhaev GO. Ionizing radiation can activate the insertion of mitochondrial DNA fragments in the nuclear genome. Radiatsionnaia Biol Radioecol. 2007;47(6):673–83.
23. Chan CY Kiechle M Manivasakam P Schiestl RH Ionizing radiation and restriction enzymes induce microhomology-mediated illegitimate recombination in Saccharomyces cerevisiae Nucleic Acids Res 2007 35 15 5051 9 10.1093/nar/gkm442 17652322
Chan CY, Kiechle M, Manivasakam P, Schiestl RH. Ionizing radiation and restriction enzymes induce microhomology-mediated illegitimate recombination in Saccharomyces cerevisiae. Nucleic Acids Res. 2007;35(15):5051–9.17652322 10.1093/nar/gkm442
24. Caro P Gómez J Arduini A González-Sánchez M González-Garc\’\ia M Borrás C Mitochondrial DNA sequences are present inside nuclear DNA in rat tissues and increase with age Mitochondrion 2010 10 5 479 86 10.1016/j.mito.2010.05.004 20546951
Caro P, Gómez J, Arduini A, González-Sánchez M, González-Garc\’\ia M, Borrás C, et al. Mitochondrial DNA sequences are present inside nuclear DNA in rat tissues and increase with age. Mitochondrion. 2010;10(5):479–86.20546951 10.1016/j.mito.2010.05.004
25. Cheng X Ivessa AS The migration of mitochondrial DNA fragments to the nucleus affects the chronological aging process of Saccharomyces cerevisiae Aging Cell 2010 9 5 919 23 10.1111/j.1474-9726.2010.00607.x 20626726
Cheng X, Ivessa AS. The migration of mitochondrial DNA fragments to the nucleus affects the chronological aging process of Saccharomyces cerevisiae. Aging Cell. 2010;9(5):919–23.20626726 10.1111/j.1474-9726.2010.00607.x
26. Goldman SJ Taylor R Zhang Y Jin S Autophagy and the degradation of mitochondria Mitochondrion 2010 10 4 309 15 10.1016/j.mito.2010.01.005 20083234
Goldman SJ, Taylor R, Zhang Y, Jin S. Autophagy and the degradation of mitochondria. Mitochondrion. 2010;10(4):309–15.20083234 10.1016/j.mito.2010.01.005
27. Hazkani-Covo E Zeller RM Martin W Molecular poltergeists: mitochondrial DNA copies (numts) in sequenced nuclear genomes PLoS Genet 2010 6 2 e1000834 10.1371/journal.pgen.1000834 20168995
Hazkani-Covo E, Zeller RM, Martin W. Molecular poltergeists: mitochondrial DNA copies (numts) in sequenced nuclear genomes. PLoS Genet. 2010;6(2):e1000834.20168995 10.1371/journal.pgen.1000834
28. Maude H, Davidson M, Charitakis N, Diaz L, Bowers T, Gradovich WH. E, NUMT confounding biases mitochondrial heteroplasmy calls in Favor of the reference allele. Front Cell Dev Biol. 2019 [cited 2023 Nov 20];7:201. Available from: www.frontiersin.org.
29. Triant DA, Pearson WR. Comparison of detection methods and genome quality when quantifying nuclear mitochondrial insertions in vertebrate genomes. Front Genet. 2022;13.
30. Martínez M, Harms L, Abele D, Held C. Mitochondrial Heteroplasmy and PCR Amplification Bias Lead to Wrong Species Delimitation with High Confidence in the South American and Antarctic Marine Bivalve Aequiyoldia eightsii Species Complex. Genes (Basel). 2023 Apr 1 [cited 2023 Nov 20];14(4):935. Available from: https://www.mdpi.com/2073-4425/14/4/935/htm.
31. Zhang G Geng D Guo Q Liu W Li S Gao W Genomic landscape of mitochondrial DNA insertions in 23 bat genomes: characteristics, loci, phylogeny, and polymorphism Integr Zool 2022 17 5 890 903 10.1111/1749-4877.12582 34496458
Zhang G, Geng D, Guo Q, Liu W, Li S, Gao W, et al. Genomic landscape of mitochondrial DNA insertions in 23 bat genomes: characteristics, loci, phylogeny, and polymorphism. Integr Zool. 2022;17(5):890–903.34496458 10.1111/1749-4877.12582
32. Grau ET Charles M Féménia M Rebours E Vaiman A Rocha D Survey of mitochondrial sequences integrated into the bovine nuclear genome Sci Rep 2020 10 1 1 11 10.1038/s41598-020-59155-4 31913322
Grau ET, Charles M, Féménia M, Rebours E, Vaiman A, Rocha D. Survey of mitochondrial sequences integrated into the bovine nuclear genome. Sci Rep. 2020;10(1):1–11.31913322 10.1038/s41598-020-59155-4
33. Wei W Schon KR Elgar G Orioli A Tanguy M Giess A Nuclear-embedded mitochondrial DNA sequences in 66,083 human genomes Nature 2022 611 7934 105 14 10.1038/s41586-022-05288-7 36198798
Wei W, Schon KR, Elgar G, Orioli A, Tanguy M, Giess A, et al. Nuclear-embedded mitochondrial DNA sequences in 66,083 human genomes. Nature. 2022;611(7934):105–14.36198798 10.1038/s41586-022-05288-7
34. Deininger P Alu elements: know the SINEs Genome Biol 2011 12 12 1 12 10.1186/gb-2011-12-12-236
Deininger P. Alu elements: know the SINEs. Genome Biol. 2011;12(12):1–12.10.1186/gb-2011-12-12-236
35. Han G Zhang N Jiang H Meng X Qian K Zheng Y Diversity of short interspersed nuclear elements (SINEs) in lepidopteran insects and evidence of horizontal SINE transfer between baculovirus and lepidopteran hosts BMC Genomics 2021 22 1 1 16 10.1186/s12864-021-07543-z 33388042
Han G, Zhang N, Jiang H, Meng X, Qian K, Zheng Y, et al. Diversity of short interspersed nuclear elements (SINEs) in lepidopteran insects and evidence of horizontal SINE transfer between baculovirus and lepidopteran hosts. BMC Genomics. 2021;22(1):1–16.33388042 10.1186/s12864-021-07543-z
36. Tarallo V Hirano Y Gelfand BD Dridi S Kerur N Kim Y DICER1 loss and Alu RNA induce age-related macular degeneration via the NLRP3 inflammasome and MyD88 Cell 2012 149 4 847 59 10.1016/j.cell.2012.03.036 22541070
Tarallo V, Hirano Y, Gelfand BD, Dridi S, Kerur N, Kim Y, et al. DICER1 loss and Alu RNA induce age-related macular degeneration via the NLRP3 inflammasome and MyD88. Cell. 2012;149(4):847–59.22541070 10.1016/j.cell.2012.03.036
37. Kerur N Fukuda S Banerjee D Kim Y Fu D Apicella I cGAS drives noncanonical-inflammasome activation in age-related macular degeneration Nat Med 2018 24 1 50 61 10.1038/nm.4450 29176737
Kerur N, Fukuda S, Banerjee D, Kim Y, Fu D, Apicella I, et al. cGAS drives noncanonical-inflammasome activation in age-related macular degeneration. Nat Med. 2018;24(1):50–61.29176737 10.1038/nm.4450
38. Larsen PA Hunnicutt KE Larsen RJ Yoder AD Saunders AM Warning SINEs: Alu elements, evolution of the human brain, and the spectrum of neurological disease Chromosom Res 2018 26 1 93 111 10.1007/s10577-018-9573-4
Larsen PA, Hunnicutt KE, Larsen RJ, Yoder AD, Saunders AM. Warning SINEs: Alu elements, evolution of the human brain, and the spectrum of neurological disease. Chromosom Res. 2018;26(1):93–111.10.1007/s10577-018-9573-4
39. Tao J Wang Q Mendez-Dorantes C Burns KH Chiarle R Frequency and mechanisms of LINE-1 retrotransposon insertions at CRISPR/Cas9 sites Nat Commun 2022 13 1 1 17 10.1038/s41467-022-31322-3 34983933
Tao J, Wang Q, Mendez-Dorantes C, Burns KH, Chiarle R. Frequency and mechanisms of LINE-1 retrotransposon insertions at CRISPR/Cas9 sites. Nat Commun. 2022;13(1):1–17.34983933 10.1038/s41467-022-31322-3
40. Taylor MS, Altukhov I, Molloy KR, Mita P, Jiang H, Adney EM et al. Dissection of affinity captured LINE-1 macromolecular complexes. Elife. 2018;7.
41. Hazkani-Covo E Covo S Numt-mediated double-strand break repair mitigates deletions during primate genome evolution PLoS Genet 2008 4 10 e1000237 10.1371/journal.pgen.1000237 18949041
Hazkani-Covo E, Covo S. Numt-mediated double-strand break repair mitigates deletions during primate genome evolution. PLoS Genet. 2008;4(10):e1000237.18949041 10.1371/journal.pgen.1000237
42. Kojima KK Structural and sequence diversity of eukaryotic transposable elements Genes \& Genet Syst 2019 94 6 233 52 10.1266/ggs.18-00024
Kojima KK. Structural and sequence diversity of eukaryotic transposable elements. Genes \& Genet Syst. 2019;94(6):233–52.10.1266/ggs.18-00024
43. Di Giorgio E, Xodo LE. Endogenous retroviruses (ERVs): does RLR (RIG-I-Like Receptors)-MAVS Pathway directly Control Senescence and Aging as a consequence of ERV De-repression? Front Immunol. 2022;13.
44. Díaz-Carballo D Klein J Acikelli AH Wilk C Saka S Jastrow H Cytotoxic stress induces transfer of mitochondria-associated human endogenous retroviral RNA and proteins between cancer cells Oncotarget 2017 8 56 95945 10.18632/oncotarget.21606 29221178
Díaz-Carballo D, Klein J, Acikelli AH, Wilk C, Saka S, Jastrow H, et al. Cytotoxic stress induces transfer of mitochondria-associated human endogenous retroviral RNA and proteins between cancer cells. Oncotarget. 2017;8(56):95945.29221178 10.18632/oncotarget.21606
45. Ko YJ Kim S Analysis of nuclear mitochondrial DNA segments of nine plant species: size, distribution, and insertion loci Genomics \& Inf 2016 14 3 90 10.5808/GI.2016.14.3.90
Ko YJ, Kim S. Analysis of nuclear mitochondrial DNA segments of nine plant species: size, distribution, and insertion loci. Genomics \& Inf. 2016;14(3):90.10.5808/GI.2016.14.3.90
46. Song S Jiang F Yuan J Guo W Miao Y Exceptionally high cumulative percentage of NUMTs originating from linear mitochondrial DNA molecules in the Hydra magnipapillata genome BMC Genomics 2013 14 1 1 13 10.1186/1471-2164-14-447 23323973
Song S, Jiang F, Yuan J, Guo W, Miao Y. Exceptionally high cumulative percentage of NUMTs originating from linear mitochondrial DNA molecules in the Hydra magnipapillata genome. BMC Genomics. 2013;14(1):1–13.23323973 10.1186/1471-2164-14-447
47. Balciuniene J Balciunas D A nuclear mtDNA concatemer (mega-NUMT) could mimic paternal inheritance of mitochondrial genome Front Genet 2019 10 518 10.3389/fgene.2019.00518 31244882
Balciuniene J, Balciunas D. A nuclear mtDNA concatemer (mega-NUMT) could mimic paternal inheritance of mitochondrial genome. Front Genet. 2019;10:518.31244882 10.3389/fgene.2019.00518
48. Nacer DF do Amaral FR Striking pseudogenization in avian phylogenetics: numts are large and common in falcons Mol Phylogenet Evol 2017 115 1 6 10.1016/j.ympev.2017.07.002 28690127
Nacer DF, do Amaral FR. Striking pseudogenization in avian phylogenetics: numts are large and common in falcons. Mol Phylogenet Evol. 2017;115:1–6.28690127 10.1016/j.ympev.2017.07.002
49. Parakatselaki ME Zhu CT Rand D Ladoukakis ED NUMTs can imitate biparental transmission of mtDNA—A case in Drosophila melanogaster Genes (Basel) 2022 13 6 1023 10.3390/genes13061023 35741785
Parakatselaki ME, Zhu CT, Rand D, Ladoukakis ED. NUMTs can imitate biparental transmission of mtDNA—A case in Drosophila melanogaster. Genes (Basel). 2022;13(6):1023.35741785 10.3390/genes13061023
50. Lucas T Vincent B Eric P Translocation of mitochondrial DNA into the nuclear genome blurs phylogeographic and conservation genetic studies in seabirds R Soc Open Sci 2022 9 6 211888 10.1098/rsos.211888 35719890
Lucas T, Vincent B, Eric P. Translocation of mitochondrial DNA into the nuclear genome blurs phylogeographic and conservation genetic studies in seabirds. R Soc Open Sci. 2022;9(6):211888.35719890 10.1098/rsos.211888
51. Baral K, Rotwein P. ZMAT2 in Humans and Other Primates: A Highly Conserved and Understudied Gene. Evol Bioinform Online. 2020 [cited 2023 Aug 28];16. Available from: https://pubmed.ncbi.nlm.nih.gov/32952394/.
52. Cartault F, Munier P, Benko E, Desguerre I, Hanein S, Boddaert N, S A. Mutation in a primate-conserved retrotransposon reveals a noncoding RNA as a mediator of infantile encephalopathy. Proc Natl Acad Sci U. 2012 Mar 27 [cited 2023 Aug 28];109(13):4980–5. Available from: https://www.pnas.org/doi/abs/10.1073/pnas.1111596109.
53. Raina SZ, Faith JJ, Disotell TR, Seligmann H, Stewart CB, Pollock DD. Evolution of base-substitution gradients in primate mitochondrial genomes. Genome Res. 2005 May 1 [cited 2023 Aug 28];15(5):665–73. Available from: https://genome.cshlp.org/content/15/5/665.full.
54. Lavrov DV, Pett W. Animal Mitochondrial DNA as We Do Not Know It: mt-Genome Organization and Evolution in Nonbilaterian Lineages. Genome Biol Evol. 2016 Sep 1 [cited 2023 Aug 28];8(9):2896–913. Available from: https://pubmed.ncbi.nlm.nih.gov/27557826/.
55. Piganeau G, Eyre-Walker A. Evidence for variation in the effective population size of animal mitochondrial DNA. PLoS One. 2009 Feb 9 [cited 2023 Aug 28];4(2). Available from: https://pubmed.ncbi.nlm.nih.gov/19198657/.
56. Mori S, Matsunami M. Signature of positive selection in mitochondrial DNA in Cetartiodactyla. Genes Genet Syst. 2018 Apr 1 [cited 2023 Aug 28];93(2):65–73. 10.1266/ggs.17-00015.
57. Barrell BG Bankier AT Drouin J A different genetic code in human mitochondria Nature 1979 282 5735 189 94 10.1038/282189a0 226894
Barrell BG, Bankier AT, Drouin J. A different genetic code in human mitochondria. Nature. 1979;282(5735):189–94.226894 10.1038/282189a0
58. Chen Z, Liu X, Zhao P, Li C, Wang Y, Li F et al. iFeatureOmega: an integrative platform for engineering, visualization and analysis of features from molecular sequences, structural and ligand data sets. Nucleic Acids Res. 2022.
59. Caballero J Fernandez L Abreu JI Fernández M Amino acid sequence autocorrelation vectors and ensembles of bayesian-regularized genetic neural networks for prediction of conformational stability of human lysozyme mutants J Chem Inf Model 2006 46 3 1255 68 10.1021/ci050507z 16711745
Caballero J, Fernandez L, Abreu JI, Fernández M. Amino acid sequence autocorrelation vectors and ensembles of bayesian-regularized genetic neural networks for prediction of conformational stability of human lysozyme mutants. J Chem Inf Model. 2006;46(3):1255–68.16711745 10.1021/ci050507z
60. Palimkar P, Shaw R, Ghosh A. In. Machine learning technique to Prognosis Diabetes Disease: Random Forest Classifier Approach. 2021. p. 219–44.
61. Xin Y, Ren X. Predicting depression among rural and urban disabled elderly in China using a random forest classifier. BMC Psychiatry. 2022;22.
62. Gao Y, Zhu Z, Sun F. Increasing prediction performance of colorectal cancer disease status using random forests classification based on metagenomic shotgun sequencing data. Synth Syst Biotechnol. 2022 [cited 2023 Nov 20];7:574–85. 10.1016/j.synbio.2022.01.005.
63. Fekete JT, Győrffy B. New Transcriptomic Biomarkers of 5-Fluorouracil Resistance. Int J Mol Sci. 2023 Jan 1 [cited 2023 Nov 20];24(2):1508. Available from: https://www.mdpi.com/1422-0067/24/2/1508/htm.
64. Kiełbasa SM Wan R Sato K Horton P Frith MC Adaptive seeds tame genomic sequence comparison Genome Res 2011 21 3 487 93 10.1101/gr.113985.110 21209072
Kiełbasa SM, Wan R, Sato K, Horton P, Frith MC. Adaptive seeds tame genomic sequence comparison. Genome Res. 2011;21(3):487–93.21209072 10.1101/gr.113985.110
65. Tsuji J Frith MC Tomii K Horton P Mammalian NUMT insertion is non-random Nucleic Acids Res 2012 40 18 9073 88 10.1093/nar/gks424 22761406
Tsuji J, Frith MC, Tomii K, Horton P. Mammalian NUMT insertion is non-random. Nucleic Acids Res. 2012;40(18):9073–88.22761406 10.1093/nar/gks424
66. Sievers F, Higgins DG. Clustal Omega, accurate alignment of very large numbers of sequences. Multiple sequence alignment methods. Springer; 2014. pp. 105–16.
67. Chen N Using repeat Masker to identify repetitive elements in genomic sequences Curr Protoc Bioinforma 2004 5 1 4 10 10.1002/0471250953.bi0410s05
Chen N. Using repeat Masker to identify repetitive elements in genomic sequences. Curr Protoc Bioinforma. 2004;5(1):4–10.10.1002/0471250953.bi0410s05
68. Zhang R Zhang CT A brief review: the Z-curve theory and its application in genome analysis Curr Genomics 2014 15 78 94 10.2174/1389202915999140328162433 24822026
Zhang R, Zhang CT. A brief review: the Z-curve theory and its application in genome analysis. Curr Genomics. 2014;15:78–94.24822026 10.2174/1389202915999140328162433
69. Kramer O. Scikit-learn. Machine learning for evolution strategies. Springer; 2016. pp. 45–53.
70. Li H Handsaker B Wysoker A Fennell T Ruan J Homer N The sequence alignment/map format and SAMtools Bioinformatics 2009 25 16 2078 9 10.1093/bioinformatics/btp352 19505943
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.19505943 10.1093/bioinformatics/btp352
71. Paradis E Claude J Strimmer K APE: analyses of phylogenetics and evolution in R language Bioinformatics 2004 20 2 289 90 10.1093/bioinformatics/btg412 14734327
Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20(2):289–90.14734327 10.1093/bioinformatics/btg412
72. Schliep KP Phangorn: phylogenetic analysis in R Bioinformatics 2011 27 4 592 3 10.1093/bioinformatics/btq706 21169378
Schliep KP. Phangorn: phylogenetic analysis in R. Bioinformatics. 2011;27(4):592–3.21169378 10.1093/bioinformatics/btq706
73. Virtanen P Gommers R Oliphant TE Haberland M Reddy T Cournapeau D SciPy 1.0: fundamental algorithms for scientific computing in Python Nat Methods 2020 17 3 261 72 10.1038/s41592-019-0686-2 32015543
Virtanen P, Gommers R, Oliphant TE, Haberland M, Reddy T, Cournapeau D, et al. SciPy 1.0: fundamental algorithms for scientific computing in Python. Nat Methods. 2020;17(3):261–72.32015543 10.1038/s41592-019-0686-2
