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

Abstract

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 3kb. 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.

Comparative mitogenome analysis of two ectomycorrhizal fungi (Paxillus) reveals gene rearrangement, intron dynamics, and phylogeny of basidiomycetes

“…During the evolution of eukaryotes, most mitochondrial genes have been transferred to the nucleus (Adams and Palmer 2003 ). However, some nuclear genes have also been found transferred to the mitogenome, and this has affected the evolution of eukaryotes (Li et al 2019f ; Srirattana and St John 2018 ). In our study, 10.64 and 36.50% of gene fragments were found potentially transferred between nuclear and mitochondrial genomes of the two Paxillus , respectively, which is a higher proportion than reported from Ganoderma lucidum (Li et al 2013 ) and Pleurotus eryngii (Li et al 2018a ).…”

Potential Predictive Factors for Breast Cancer Subtypes from a North Cyprus Cohort Analysis

“…The t-distributed stochastic neighbor embedding (t-SNE) [13] is a popular method to represent high-dimensional data in 2 or 3 dimensions. We previously use the technique in other applications in biology/genomics [14, 15].…”

COVID-19 outpatients and surviving inpatients exhibit comparable blood test results that are distinct from non-surviving inpatients

“…To check this, we first consider each person as a point in the 9-dimensional space (for nine blood test variables, and note the age factor is not used), then project the points to low-dimensional space by a dimension-reduction technique. We use both t-SNE and UMAP, which are popular in single-cell expression data analysis as well as other fields (14,15,(19)(20)(21)(22)(23).…”