Pervasive heteroplasmy in an invasive ambrosia beetle (Scolytinae) in southern California

Heteroplasmy, the presence of multiple mitochondrial genotypes (mitotypes) within an individual, has long been thought to be a rare aberrance that is quickly removed by selection or drift. However, heteroplasmy is being reported in natural populations of eukaryotes with increasing frequency, in part due to improved diagnostic methods. Here, we report a seemingly stable heteroplasmic state in California populations of the polyphagous shothole borer (PSHB), Euwallacea fornicatus; an invasive ambrosia beetle that is causing significant tree dieback. We develop and validate a qPCR assay utilizing locked nucleic acid probes to detect different mitotypes, and qualitatively assess heteroplasmy in individual PSHB. We prove the utility of this assay by: (1) mitotyping field-collected PSHB, documenting the prevalence of heteroplasmy across its range in California; and, (2) measuring relative titers of each mitotype across multiple generations of heteroplasmic laboratory colonies to assess the stability of transmission through the maternal germline. We show that our findings are unlikely to be explained by the existence of NUMTs by next generation sequencing of contiguous sections of mitochondrial DNA, where each of the observed heteroplasmic sites are found within fully functional coding regions of mtDNA. Subsequently, we find heteroplasmic individuals are common in Californian field populations, and that heteroplasmy persists for at least 10 generations in experimental colonies. We also looked for evidence of the common occurrence of paternal leakage, but found none. In light of our results, we discuss competing hypotheses as to how heteroplasmy may have arisen, and continues to perpetuate, in Californian PSHB populations.

Assembly and comparative analysis of the mitochondrial genome in diploid potatoes

KEY MESSAGE

We report the mitochondrial genome of 39 diploid potatoes and identify a candidate ORF potentially linked to cytoplasmic male sterility in potatoes. Potato (Solanum tuberosum L.) holds a critical position as the foremost non-grain food crop, playing a pivotal role in ensuring global food security. Diploid potatoes constitute a vital genetic resource pool, harboring the potential to revolutionize modern potato breeding. Nevertheless, diploid potatoes are relatively understudied, and mitochondrial DNA can provide valuable insights into key potato breeding traits such as CMS. In this study, we examine and assemble the mitochondrial genome evolution and diversity of 39 accessions of diploid potatoes using high-fidelity (HiFi) sequencing. We annotated 54 genes for all the investigated accessions, comprising 34 protein-coding genes, 3 rRNA genes, and 17 tRNA genes. Our analyses revealed differences in repeats sequences between wild and cultivated landraces. To understand the evolution of diploid maternal lineage inheritance, we conducted phylogenetic analysis, which clearly distinguished mitochondrial from nuclear gene trees, further supporting the evidence-based of clustering between wild and cultivated landraces accessions. Our study discovers new candidate ORFs associated with CMS in potatoes, including ORF137, which is homologous to other CMS in Solanaceae. Ultimately, this work bridges the gap in mitochondrial genome research for diploid potatoes, providing a steppingstone into evolutionary studies and potato breeding.

Liberties of the genome: insertions of mitochondrial DNA fragments into nuclear genome

The transition of detached fragments of mitochondrial DNA into the nucleus and their integration into chromosomal DNA is a special kind of genetic variability that highlights the relation between the two genomes and their interaction in a eukaryotic cell. The human genome contains several hundreds of insertions of mtDNA fragments (NUMTS). This paper presents an overview of the current state of research in this area. To date, evidence has been obtained that the occurrence of new mtDNA insertions in the nuclear genome is a seldom but not exceptionally rare event. The integration of new mtDNA fragments into the nuclear genome occurs during double-strand DNA break repair through the non-homologous end joining mechanism. Along with evolutionarily stable "genetic fossils" that were integrated into the nuclear genome millions of years ago and are shared by many species, there are NUMTS that could be species-specific, polymorphic in a species, or "private". Partial copies of mitochondrial DNA in the human nuclear genome can interfere with mtDNA during experimental studies of the mitochondrial genome, such as genotyping, heteroplasmy assessment, mtDNA methylation analysis, and mtDNA copy number estimation. In some cases, the insertion of multiple copies of the complete mitochondrial genome sequence may mimic paternal inheritance of mtDNA. The functional significance of NUMTS is poorly understood. For instance, they may be a source of variability for expression and splicing modulation. The role of NUMTS as a cause of hereditary diseases is negligible, since only a few cases of diseases caused by NUMTS have been described so far. In addition, NUMTS can serve as markers for evolutionary genetic studies. Of particular interest is the meaning of NUMTS in eukaryotic genome evolution. The constant flow of functionally inactive DNA sequences from mitochondria into the nucleus and its significance could be studied in view of the modern concepts of evolutionary theory suggesting non-adaptive complexity and the key role of stochastic processes in the formation of genomic structure.

Strong mitonuclear discordance in the phylogeny of Neodermata and evolutionary rates of Polyopisthocotylea

The genomic evolution of Polyopisthocotylea remains poorly understood in comparison to the remaining three classes of Neodermata: Monopisthocotylea, Cestoda, and Trematoda. Moreover, the evolutionary sequence of major events in the phylogeny of Neodermata remains unresolved. Herein we sequenced the mitogenome and transcriptome of the polyopisthocotylean Diplorchis sp., and conducted comparative evolutionary analyses using nuclear (nDNA) and mitochondrial (mtDNA) genomic datasets of Neodermata. We found strong mitonuclear discordance in the phylogeny of Neodermata. Polyopisthocotylea exhibited striking mitonuclear discordance in relative evolutionary rates: the fastest-evolving mtDNA in Neodermata and a comparatively slowly-evolving nDNA genome. This was largely attributable to its very long stem branch in mtDNA topologies, not exhibited by the nDNA data. We found indications that the fast evolution of mitochondrial genomes of Polyopisthocotylea may be driven both by relaxed purifying selection pressures and elevated levels of directional selection. We identified mitochondria-associated genes encoded in the nuclear genome: they exhibited unique evolutionary rates, but not correlated with the evolutionary rate of mtDNA, and there is no evidence for compensatory evolution (they evolved slower than the rest of the genome). Finally, there appears to exist an exceptionally large (≈6.3 kb) nuclear mitochondrial DNA segment (numt) in the nuclear genome of newly sequenced Diplorchis sp. A 3'-end segment of the 16S rRNA gene encoded by the numt was expressed, suggesting that this gene acquired novel, regulatory functions after the transposition to the nuclear genome. In conclusion, Polyopisthocotylea appears to be the lineage with the fastest-evolving mtDNA sequences among all of Bilateria, but most of the substitutions were accumulated deep in the evolutionary history of this lineage. As the nuclear genome does not exhibit a similar pattern, the circumstances underpinning this evolutionary phenomenon remain a mystery.

Epigenetic signatures of social status in wild female spotted hyenas (Crocuta crocuta)

In mammalian societies, dominance hierarchies translate into inequalities in health, reproductive performance and survival. DNA methylation is thought to mediate the effects of social status on gene expression and phenotypic outcomes, yet a study of social status-specific DNA methylation profiles in different age classes in a wild social mammal is missing. We tested for social status signatures in DNA methylation profiles in wild female spotted hyenas (Crocuta crocuta), cubs and adults, using non-invasively collected gut epithelium samples. In spotted hyena clans, female social status influences access to resources, foraging behavior, health, reproductive performance and survival. We identified 149 differentially methylated regions between 42 high- and low-ranking female spotted hyenas (cubs and adults). Differentially methylated genes were associated with energy conversion, immune function, glutamate receptor signalling and ion transport. Our results provide evidence that socio-environmental inequalities are reflected at the molecular level in cubs and adults in a wild social mammal.

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

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

HaploCart: Human mtDNA haplogroup classification using a pangenomic reference graph human mtDNA haplogroup inference

Current mitochondrial DNA (mtDNA) haplogroup classification tools map reads to a single reference genome and perform inference based on the detected mutations to this reference. This approach biases haplogroup assignments towards the reference and prohibits accurate calculations of the uncertainty in assignment. We present HaploCart, a probabilistic mtDNA haplogroup classifier which uses a pangenomic reference graph framework together with principles of Bayesian inference. We demonstrate that our approach significantly outperforms available tools by being more robust to lower coverage or incomplete consensus sequences and producing phylogenetically-aware confidence scores that are unbiased towards any haplogroup. HaploCart is available both as a command-line tool and through a user-friendly web interface. The C++ program accepts as input consensus FASTA, FASTQ, or GAM files, and outputs a text file with the haplogroup assignments of the samples along with the level of confidence in the assignments. Our work considerably reduces the amount of data required to obtain a confident mitochondrial haplogroup assignment.

Bringing to light nuclear-mitochondrial insertions in the genomes of nocturnal predatory birds

Mito-nuclear insertions, or NUMTs, relate to genetic material of mitochondrial origin that have been transferred to the nuclear DNA molecule. The increasing amounts of genomic data currently being produced presents an opportunity to investigate this type of patterns in genome evolution of non-model organisms. Identifying NUMTs across a range of closely related taxa allows one to generalize patterns of insertion and maintenance in autosomes, which is ultimately relevant to the understanding of genome biology and evolution. Here we collected existing pairwise genome-mitogenome data of the order Strigiformes, a group that includes all the nocturnal bird predators. We identified NUMTs by applying percent similarity thresholds after blasting mitochondrial genomes against nuclear genome assemblies. We identified NUMTsin all genomes with numbers ranging from 4 in Bubo bubo to 24 in Ciccaba nigrolineata. Statistical analyses revealed NUMT size to negatively correlate with NUMT's sequence similarity to with original mtDNA region. Lastly, characterizing these nuclear insertions of mitochondrial origin in a comparative genomics framework produced variable phylogenetic patterns, suggesting in some cases that insertions might pre-date speciation events within Strigiformes.

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.

Contribution bias of parental genomes to the hybrid lineages of black Amur bream and topmouth culter revealed by low-coverage whole-genome sequencing

Interspecific hybridization has the potential to increase animal genetic diversity, enable the introgression of advantageous genetic variation, and even create hybrid lineages. Using the low-coverage whole-genome sequencing data of black Amur bream (Megalobrama terminalis [SJ]), topmouth culter (Culter alburnus [QZ]), and their hybrid progeny, we focused on the sequence variation and contribution bias to the genomic composition of the hybrid lineages of SJ♀ × QZ♂ [ZJ] and QZ♀ × SJ♂ [FJ]. Both mash-based genetic distance analysis and variant analysis based on various reference genomes revealed that the hybrid lineages' genomes are more prone to QZ composition. Furthermore, we assess the degree of subgenomic intermixing between parents QZ and SJ in the FJ and ZJ hybrid lineages. In the Jaccard index-zero genomic area, 89.41% of the genomic region in ZJ and 90.23% in FJ are entirely from one of the parental genomes. This genomic area is primarily from QZ, regardless of FJ or ZJ. In FJ, the proportion is 55.70%; in ZJ, the proportion is 54.60%. Using the reference-free variant discovery tool, we identified two InDel markers (Both were SSR-type InDel markers) and one (T/A)-type SNP marker that displayed differences in agarose electrophoresis. Our current findings suggest that parental genomic contribution bias leads to an imbalance in the genomic composition of hybrid lineages. The genomic composition bias is connected to species but not maternal effects, giving insights into the development of vertebrate genomes shortly after hybridization.

Mitochondrial signal transduction

The analogy of mitochondria as powerhouses has expired. Mitochondria are living, dynamic, maternally inherited, energy-transforming, biosynthetic, and signaling organelles that actively transduce biological information. We argue that mitochondria are the processor of the cell, and together with the nucleus and other organelles they constitute the mitochondrial information processing system (MIPS). In a three-step process, mitochondria (1) sense and respond to both endogenous and environmental inputs through morphological and functional remodeling; (2) integrate information through dynamic, network-based physical interactions and diffusion mechanisms; and (3) produce output signals that tune the functions of other organelles and systemically regulate physiology. This input-to-output transformation allows mitochondria to transduce metabolic, biochemical, neuroendocrine, and other local or systemic signals that enhance organismal adaptation. An explicit focus on mitochondrial signal transduction emphasizes the role of communication in mitochondrial biology. This framework also opens new avenues to understand how mitochondria mediate inter-organ processes underlying human health.

Frequent Paternal Mitochondrial Inheritance and Rapid Haplotype Frequency Shifts in Copepod Hybrids

Mitochondria are assumed to be maternally inherited in most animal species, and this foundational concept has fostered advances in phylogenetics, conservation, and population genetics. Like other animals, mitochondria were thought to be solely maternally inherited in the marine copepod Tigriopus californicus, which has served as a useful model for studying mitonuclear interactions, hybrid breakdown, and environmental tolerance. However, we present PCR, Sanger sequencing, and Illumina Nextera sequencing evidence that extensive paternal mitochondrial DNA (mtDNA) transmission is occurring in inter-population hybrids of T. californicus. PCR on four types of crosses between three populations (total sample size of 376 F1 individuals) with 20% genome-wide mitochondrial divergence showed 2% to 59% of F1 hybrids with both paternal and maternal mtDNA, where low and high paternal leakage values were found in different cross directions of the same population pairs. Sequencing methods further verified nucleotide similarities between F1 mtDNA and paternal mtDNA sequences. Interestingly, the paternal mtDNA in F1s from some crosses inherited haplotypes that were uncommon in the paternal population. Compared to some previous research on paternal leakage, we employed more rigorous methods to rule out contamination and false detection of paternal mtDNA due to non-functional nuclear mitochondrial DNA fragments. Our results raise the potential that other animal systems thought to only inherit maternal mitochondria may also have paternal leakage, which would then affect the interpretation of past and future population genetics or phylogenetic studies that rely on mitochondria as uniparental markers.

Comparative hologenomics of two Ixodes scapularis tick populations in New Jersey

Tick-borne diseases, such as those transmitted by the blacklegged tick Ixodes scapularis, are a significant and growing public health problem in the US. There is mounting evidence that co-occurring non-pathogenic microbes can also impact tick-borne disease transmission. Shotgun metagenome sequencing enables sampling of the complete tick hologenome—the collective genomes of the tick and all of the microbial species contained therein, whether pathogenic, commensal or symbiotic. This approach simultaneously uncovers taxonomic composition and allows the detection of intraspecific genetic variation, making it a useful tool to compare spatial differences across tick populations. We evaluated this approach by comparing hologenome data from two tick samples (N = 6 ticks per location) collected at a relatively fine spatial scale, approximately 23 km apart, within a single US county. Several intriguing variants in the data between the two sites were detected, including polymorphisms in both in the tick’s own mitochondrial DNA and that of a rickettsial endosymbiont. The two samples were broadly similar in terms of the microbial species present, including multiple known tick-borne pathogens (Borrelia burgdorferi, Babesia microti, and Anaplasma phagocytophilum), filarial nematodes, and Wolbachia and Babesia species. We assembled the complete genome of the rickettsial endosymbiont (most likely Rickettsia buchneri) from both populations. Our results provide further evidence for the use of shotgun metagenome sequencing as a tool to compare tick hologenomes and differentiate tick populations across localized spatial scales.

Detection of heteroplasmy and nuclear mitochondrial pseudogenes in the Japanese spiny lobster Panulirus japonicus

Partial mtDNA cytochrome oxidase subunit I (COI) fragments and near entire stretch of 12S rDNA (12S) and control region (Dloop) of the Japanese spiny lobster (Panulirus japonicus) (n = 3) were amplified by PCR and used for direct nucleotide sequencing and for clone library-based nucleotide sequence analysis. Nucleotide sequences of a total of 75 clones in COI, 77 in 12S and 92 in Dloop were determined. Haplotypes of the clones matched with those obtained by direct sequencing were determined to be genuine mtDNA sequence of the individual. Phylogenetic analysis revealed several distinct groups of haplotypes in all three regions. Genuine mtDNA sequences were observed to form a group with their closely related variables, and most of these variables may be due to amplification error but a few to be heteroplasmy. Haplotypes determined as nuclear mitochondrial pseudogenes (NUMTs) formed distinct groups. Nucleotide sequence divergence (K2P distance) between genuine haplotypes and NUMTs were substantial (7.169-23.880% for COI, 1.336-23.434% for 12S, and 7.897-71.862% for Dloop). These values were comparable to or smaller than those between species of the genus Panulirus, indicating that integration of mtDNA into the nuclear genome is a continuous and dynamic process throughout pre- and post-speciation events. Double peaks in electropherograms obtained by direct nucleotide sequencing were attributed to common nucleotides shared by multiple NUMTs. Information on the heteroplasmy and NUMTs would be very important for addressing their impact on direct nucleotide sequencing and for quality control of nucleotide sequences obtained.

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.

Mechanisms for sperm mitochondrial removal in embryos

BACKGROUND

Different animal species have different characteristics regarding the transmission of mitochondrial DNA. While some species have biparental mitochondrial inheritance, others have developed pathways to remove paternal mtDNA. These pathways guarantee the uniparental mitochondrial inheritance, so far well known in mammals, avoiding heteroplasmy, which may have the potential to cause certain mitochondrial diseases in the offspring.

SCOPE OF REVIEW

This review aims to address the main mechanisms that involve mitochondrial degradation in different animal species, as well as to describe what is present in the literature on the mechanisms involved in mitochondrial inheritance.

MAJOR CONCLUSIONS

Two theories are proposed to explain the uniparental inheritance of mtDNA: (i) active degradation, where mechanisms for paternal mitochondrial DNA elimination involve mitochondrial degradation pathway by autophagy and, in some species, may also involve the endocytic degradation pathway; and (ii) passive dilution, where the paternal mitochondria are diluted in the cells of the embryo according to cell division, until becoming undetectable.

GENERAL SIGNIFICANCE

This work brings a wide review of the already published evidence on mitochondrial inheritance in the animal kingdom and the possible mechanisms to mtDNA transmission already described in literature.

Identification of Mitochondrial DNA (NUMTs) in the Nuclear Genome of Daphnia magna

This is the first study in which the Daphnia magna (D. magna) nuclear genome (nDNA) obtained from the GenBank database was analyzed for pseudogene sequences of mitochondrial origin. To date, there is no information about pseudogenes localized in D. magna genome. This study aimed to identify NUMTs, their length, homology, and location for potential use in evolutionary studies and to check whether their occurrence causes co-amplification during mitochondrial genome (mtDNA) analyses. Bioinformatic analysis showed 1909 fragments of the mtDNA of D. magna, of which 1630 were located in ten linkage groups (LG) of the nDNA. The best-matched NUMTs covering >90% of the gene sequence have been identified for two mt-tRNA genes, and they may be functional nuclear RNA molecules. Isolating the total DNA in mtDNA studies, co-amplification of nDNA fragments is unlikely in the case of amplification of the whole tRNA genes as well as fragments of other genes. It was observed that TRNA-MET fragments had the highest level of sequence homology, thus they could be evolutionarily the youngest. The lowest homology was found in the D-loop-derived pseudogene. It may probably be the oldest NUMT incorporated into the nDNA; however, further analysis is necessary.

Fish-Based Baby Food Concern-From Species Authentication to Exposure Risk Assessment

In this work, two different but complementary approaches were used to evaluate the reliability of fish-based baby foods as a source of safe nourishment for babies. More specifically, barcoding analysis based on the Cytochrome Oxidase I sequences was used for fish species authentication and an analysis of metal/metalloid levels was performed to estimate the exposure risk assessment derived from consumption of selected fish-based baby food in infants and toddlers. COI DNA barcoding revealed that in three samples the species detected did not match the common name of the species shown on the label. In particular, G. chalcogrammus and M. australis were found in place of M. merluccius and O. mykiss was found in place of S. salar. The analysis of exposure risk assessment indicated a low risk for developing chronic systemic and carcinogenic effects in infants and toddler, under an exposure scenario based on daily consumption of a single box of fish-based baby food. However, it is important to highlight that in order to provide a comprehensive risk assessment it would be important to supplement the levels of exposure resulting from the total diet. Overall, our results suggest that more attention should be paid by authorities to ensure the safety of food for infants and toddlers.

Mitochondrial DNA enrichment reduced NUMT contamination in porcine NGS analyses

Genetic associations between mitochondrial DNA (mtDNA) and economic traits have been widely reported for pigs, which indicate the importance of mtDNA. However, studies on mtDNA heteroplasmy in pigs are rare. Next generation sequencing (NGS) methodologies have emerged as a promising genomic approach for detection of mitochondrial heteroplasmy. Due to the short reads, flexible bioinformatic analyses and the contamination of nuclear mitochondrial sequences (NUMTs), NGS was expected to increase false-positive detection of heteroplasmy. In this study, Sanger sequencing was performed as a gold standard to detect heteroplasmy with a detection sensitivity of 5% in pigs and then one whole-genome sequencing method (WGS) and two mtDNA enrichment sequencing methods (Capture and LongPCR) were carried out. The aim of this study was to determine whether mitochondrial heteroplasmy identification from NGS data was affected by NUMTs. We find that WGS generated more false intra-individual polymorphisms and less mapping specificity than the two enrichment sequencing methods, suggesting NUMTs indeed led to false-positive mitochondrial heteroplasmies from NGS data. In addition, to accurately detect mitochondrial diversity, three commonly used tools-SAMtools, VarScan and GATK-with different parameter values were compared. VarScan achieved the best specificity and sensitivity when considering the base alignment quality re-computation and the minimum variant frequency of 0.25. It also suggested bioinformatic workflow interfere in the identification of mtDNA SNPs. In conclusion, intra-individual polymorphism in pig mitochondria from NGS data was confused with NUMTs, and mtDNA-specific enrichment is essential before high-throughput sequencing in the detection of mitochondrial genome sequences.

The Genomic Origins of Small Mitochondrial RNAs: Are They Transcribed by the Mitochondrial DNA or by Mitochondrial Pseudogenes within the Nucleus (NUMTs)?

Several studies have linked mitochondrial genetic variation to phenotypic modifications; albeit the identity of the mitochondrial polymorphisms involved remains elusive. The search for these polymorphisms led to the discovery of small noncoding RNAs, which appear to be transcribed by the mitochondrial DNA ("small mitochondrial RNAs"). This contention is, however, controversial because the nuclear genome of most animals harbors mitochondrial pseudogenes (NUMTs) of identical sequence to regions of mtDNA, which could alternatively represent the source of these RNAs. To discern the likely contributions of the mitochondrial and nuclear genome to transcribing these small mitochondrial RNAs, we leverage data from six vertebrate species exhibiting markedly different levels of NUMT sequence. We explore whether abundances of small mitochondrial RNAs are associated with levels of NUMT sequence across species, or differences in tissue-specific mtDNA content within species. Evidence for the former would support the hypothesis these RNAs are primarily transcribed by NUMT sequence, whereas evidence for the latter would provide strong evidence for the counter hypothesis that these RNAs are transcribed directly by the mtDNA. No association exists between the abundance of small mitochondrial RNAs and NUMT levels across species. Moreover, a sizable proportion of transcripts map exclusively to the mtDNA sequence, even in species with highest NUMT levels. Conversely, tissue-specific abundances of small mitochondrial RNAs are strongly associated with the mtDNA content. These results support the hypothesis that small mitochondrial RNAs are primarily transcribed by the mitochondrial genome and that this capacity is conserved across Amniota and, most likely, across most metazoan lineages.

Beyond Biodiversity: Can Environmental DNA (eDNA) Cut It as a Population Genetics Tool?

Population genetic data underpin many studies of behavioral, ecological, and evolutionary processes in wild populations and contribute to effective conservation management. However, collecting genetic samples can be challenging when working with endangered, invasive, or cryptic species. Environmental DNA (eDNA) offers a way to sample genetic material non-invasively without requiring visual observation. While eDNA has been trialed extensively as a biodiversity and biosecurity monitoring tool with a strong taxonomic focus, it has yet to be fully explored as a means for obtaining population genetic information. Here, we review current research that employs eDNA approaches for the study of populations. We outline challenges facing eDNA-based population genetic methodologies, and suggest avenues of research for future developments. We advocate that with further optimizations, this emergent field holds great potential as part of the population genetics toolkit.

Pseudogenes as Functionally Significant Elements of the Genome

Pseudogene is a gene copy that has lost its original function. For a long time, pseudogenes have been considered as "junk DNA" that inevitably arises as a result of ongoing evolutionary process. However, experimental data obtained during recent years indicate this understanding of the nature of pseudogenes is not entirely correct, and many pseudogenes perform important genetic functions. In the review, we have addressed classification of pseudogenes, methods of their detection in the genome, and the problem of their evolutionary conservatism and prevalence among species belonging to different taxonomic groups in the light of modern data. The mechanisms of gene expression regulation by pseudogenes and the role of pseudogenes in pathogenesis of various human diseases are discussed.

Between the Balkans and the Baltic: Phylogeography of a Common Vole Mitochondrial DNA Lineage Limited to Central Europe

The common vole (Microtus arvalis) has been a model species of small mammal for studying end-glacial colonization history. In the present study we expanded the sampling from central and eastern Europe, analyzing contemporary genetic structure to identify the role of a potential 'northern glacial refugium', i.e. a refugium at a higher latitude than the traditional Mediterranean refugia. Altogether we analyzed 786 cytochrome b (cytb) sequences (representing mitochondrial DNA; mtDNA) from the whole of Europe, adding 177 new sequences from central and eastern Europe, and we conducted analyses on eight microsatellite loci for 499 individuals (representing nuclear DNA) from central and eastern Europe, adding data on 311 new specimens. Our new data fill gaps in the vicinity of the Carpathian Mountains, the potential northern refugium, such that there is now dense sampling from the Balkans to the Baltic Sea. Here we present evidence that the Eastern mtDNA lineage of the common vole was present in the vicinity of this Carpathian refugium during the Last Glacial Maximum and the Younger Dryas. The Eastern lineage expanded from this refugium to the Baltic and shows low cytb nucleotide diversity in those most northerly parts of the distribution. Analyses of microsatellites revealed a similar pattern but also showed little differentiation between all of the populations sampled in central and eastern Europe.

Functions of noncoding sequences in mammalian genomes

Most of the mammalian genome consists of nucleotide sequences not coding for proteins. Exons of genes make up only 3% of the human genome, while the significance of most other sequences remains unknown. Recent genome studies with high-throughput methods demonstrate that the so-called noncoding part of the genome may perform important functions. This hypothesis is supported by three groups of experimental data: 1) approximately 10% of the sequences, most of which are located in noncoding parts of the genome, is evolutionarily conserved and thus can be of functional importance; 2) up to 99% of the mammalian genome is being transcribed forming short and long noncoding RNAs in addition to common mRNA; and 3) mutations in noncoding parts of the genome can be accompanied by progression of pathological states of the organism. In the light of these data, in the review we consider the functional role of numerous known sequences of noncoding parts of the genome including introns, DNA methylation regions, enhancers and locus control regions, insulators, S/MAR sequences, pseudogenes, and genes of noncoding RNAs, as well as transposons and simple repeats of centromeric and telomeric regions of chromosomes. The assumption is made that the intergenic noncoding sequences without definite/clear functions can be involved in spatial organization of genetic loci in interphase nuclei.

Could pseudogenes be widespread in ants? Evidence of numts in the leafcutter ant Acromyrmex striatus (Roger, 1863) (Formicidae: Attini)

The incorporation of fragments of mitochondrial DNA (mtDNA) in the nuclear genome, known as numts (nuclear mitochondrial pseudogenes), undermines general assumptions concerning the use of mtDNA in phylogenetic and phylogeographic studies. Accidental amplifications of these nuclear copies instead of the mitochondrial target can lead to crucial misinterpretations, thus the correct identification of numts and their differentiation from true mitochondrial sequences are important in preventing this kind of error. Our goal was to describe the existence of cytochrome b (cytb) numts in the leafcutter ant Acromyrmex striatus (Roger, 1863). PCR products were directly sequenced using a pair of universal primers designed to amplify the cytb gene of these insects. Other species of leafcutter ants were also sequenced. The sequences were analyzed and the numts were identified by the presence of double peaks, indels and premature stop codons. Only A. striatus clearly showed the presence of numts, while the other species displayed the expected amplification of the mtDNA cytb gene target using the same primer pair. We hope that our report will highlight the benefits and challenges of using mtDNA in the molecular phylogenetic reconstruction and phylogeographic studies of ants, while establishing the importance of numts reports for future studies.

Primates and mouse NumtS in the UCSC Genome Browser

BACKGROUND

NumtS (Nuclear MiTochondrial Sequences) are mitochondrial DNA sequences that, after stress events involving the mitochondrion, colonized the nuclear genome. Accurate mapping of NumtS avoids contamination during mtDNA PCR amplification, thus supplying reliable bases for detecting false heteroplasmies. In addition, since they commonly populate mammalian genomes (especially primates) and are polymorphic, in terms of presence/absence and content of SNPs, they may be used as evolutionary markers in intra- and inter-species population analyses.

RESULTS

The need for an exhaustive NumtS annotation led us to produce the Reference Human NumtS compilation, followed, as reported in this paper, by those for chimpanzee, rhesus macaque and mouse ones. Identification of NumtS inside the UCSC Genome Browser and their inter-species comparison required the design and the implementation of NumtS tracks, starting from the compilation data. NumtS retrieval through the UCSC Genome Browser, in the species examined, is now feasible at a glance.

CONCLUSIONS

Analyses involving NumtS tracks, together with other genome element tracks publicly available at the UCSC Genome Browser, can provide deep insight into genome evolution and comparative genomics, thus improving studies dealing with the mechanisms that drove the generation of NumtS. In addition, the NumtS tracks constitute a useful tool in the design of mitochondrial DNA primers.

Highly conserved d-loop sequences in woolly mouse opossums Marmosa (Micoureus)

This study reports the occurrence of highly conserved d-loop sequences in the mitochondrial genome of the woolly mouse opossum genus Marmosa subgenus Micoureus (Mammalia, Didelphimorphia, Didelphidae). Sixty-six sequences of Marmosa (Micoureus) demerarae, Marmosa (Micoureus) constantiae, and Marmosa (Micoureus) paraguayanus were amplified using universal d-loop primers and virtually no genetic differences were detected within and among species. These sequences matched the control region of the mitochondrial marsupial genome. Analyses of qualitative aspects of these sequences revealed that their structural composition is very similar to the d-loop region of other didelphid species. However, the total lack of variability has not been reported from other closely related species. The data analyzed here support the occurrence of highly conserved d-loop sequences, and we found no support for the hypothesis that these sequences are d-loop-like nuclear pseudogenes. Furthermore, the control and flanking regions obtained with different primers corroborate the lack of variability of the d-loop sequences in the mitochondrial genome of Marmosa (Micoureus).

Complete sequence and analysis of the mitochondrial genome of the methylotrophic yeast Hansenula polymorpha DL-1

We determined the complete nucleotide sequence of the 41 719 bp mitochondrial genome of the methylotrophic yeast Hansenula polymorpha strain DL-1. It contains genes for three subunits of cytochrome oxidase (cox1, cox2 and cox3), three subunits of ATP synthase (atp6, atp8 and atp9), seven subunits of NADH dehydrogenase (nad1-6 and nad4L), apocytochrome b (cob), four endonuclease/maturase homologs, a ribosomal protein (rps3), large and small rRNAs and a complete set of tRNAs. The structural genes are organized in two major transcriptional units. Phylogenetic, gene content and gene order analyses revealed the close phylogenetic relationship between H. polymorpha and Brettanomyces custersianus, and support the assignment of strain DL-1 to a separate genus rather than including it in the polyphyletic genus Pichia.