Mitochondrial DNA intra-individual variation in a bumblebee species: A challenge for evolutionary studies and molecular identification

Environmental DNA: The next chapter

Molecular tools are an indispensable part of ecology and biodiversity sciences and implemented across all biomes. About a decade ago, the use and implementation of environmental DNA (eDNA) to detect biodiversity signals extracted from environmental samples opened new avenues of research. Initial eDNA research focused on understanding population dynamics of target species. Its scope thereafter broadened, uncovering previously unrecorded biodiversity via metabarcoding in both well-studied and understudied ecosystems across all taxonomic groups. The application of eDNA rapidly became an established part of biodiversity research, and a research field by its own. Here, we revisit key expectations made in a land-mark special issue on eDNA in Molecular Ecology in 2012 to frame the development in six key areas: (1) sample collection, (2) primer development, (3) biomonitoring, (4) quantification, (5) behaviour of DNA in the environment and (6) reference database development. We pinpoint the success of eDNA, yet also discuss shortfalls and expectations not met, highlighting areas of research priority and identify the unexpected developments. In parallel, our retrospective couples a screening of the peer-reviewed literature with a survey of eDNA users including academics, end-users and commercial providers, in which we address the priority areas to focus research efforts to advance the field of eDNA. With the rapid and ever-increasing pace of new technical advances, the future of eDNA looks bright, yet successful applications and best practices must become more interdisciplinary to reach its full potential. Our retrospect gives the tools and expectations towards concretely moving the field forward.

The complete mitochondrial genome of Trigonisca nataliae (Hymenoptera, Apidae) assemblage reveals heteroplasmy in the control region

The evolution of mitochondrial genomes in the stingless bees is surprisingly dynamic, making them a model system to understand mitogenome structure, function, and evolution. Out of the seven mitogenomes available in this group, five exhibit atypical characteristics, including extreme rearrangements, rapid evolution and complete mitogenome duplication. To further explore the mitogenome diversity in these bees, we utilized isolated mtDNA and Illumina sequencing to assemble the complete mitogenome of Trigonisca nataliae, a species found in Northern Brazil. The mitogenome of T. nataliae was highly conserved in gene content and structure when compared to Melipona species but diverged in the control region (CR). Using PCR amplification, cloning and Sanger sequencing, six different CR haplotypes, varying in size and content, were recovery. These findings indicate that heteroplasmy, where different mitochondrial haplotypes coexist within individuals, occurs in T. nataliae. Consequently, we argue that heteroplasmy might indeed be a common phenomenon in bees that could be associated with variations in mitogenome size and challenges encountered during the assembly process.

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.

Genetic Differentiation and Widespread Mitochondrial Heteroplasmy among Geographic Populations of the Gourmet Mushroom Thelephora ganbajun from Yunnan, China

The mitochondrial genomes are generally considered non-recombining and homoplasmic in nature. However, our previous study provided the first evidence of extensive and stable mitochondrial heteroplasmy in natural populations of the basidiomycete fungus Thelephora ganbajun from Yunnan province, China. The heteroplasmy was characterized by the presence of two types of introns residing at adjacent but different sites in the cytochrome oxidase subunits I (cox1) gene within an individual strain. However, the frequencies of these two introns among isolates from different geographical populations and the implications for the genetic structure in natural populations have not been investigated. In this study, we analyzed DNA sequence variation at the internal transcribed spacer (ITS) regions of the nuclear ribosomal RNA gene cluster among 489 specimens from 30 geographic locations from Yunnan and compared that variation with distribution patterns of the two signature introns in the cox1 gene that are indicative of heteroplasmy in this species. In our samples, evidence for gene flow, abundant genetic diversity, and genotypic uniqueness among geographic samples in Yunnan were revealed by ITS sequence variation. While there was insignificant positive correlation between geographic distance and genetic differentiation among the geographic samples based on ITS sequences, a moderate significant correlation was found between ITS sequence variation, geographical distance of sampling sites, and distribution patterns of the two heteroplasmic introns in the cox1 gene. Interestingly, there was a significantly negative correlation between the copy numbers of the two co-existing introns. We discussed the implications of our results for a better understanding of the spread of stable mitochondrial heteroplasmy, mito-nuclear interactions, and conservation of this important gourmet mushroom.

Unparalleled mitochondrial heteroplasmy and Wolbachia co-infection in the non-model bee, Amphylaeus morosus

Mitochondrial heteroplasmy is the occurrence of more than one type of mitochondrial DNA within a single individual. Although generally reported to occur in a small subset of individuals within a species, there are some instances of widespread heteroplasmy across entire populations. Amphylaeus morosus is an Australian native bee species in the diverse and cosmopolitan bee family Colletidae. This species has an extensive geographical range along the eastern Australian coast, from southern Queensland to western Victoria, covering approximately 2,000 km. Seventy individuals were collected from five localities across this geographical range and sequenced using Sanger sequencing for the mitochondrial cytochrome c oxidase subunit I (COI) gene. These data indicate that every individual had the same consistent heteroplasmic sites but no other nucleotide variation, suggesting two conserved and widespread heteroplasmic mitogenomes. Ion Torrent shotgun sequencing revealed that heteroplasmy occurred across multiple mitochondrial protein-coding genes and is unlikely explained by transposition of mitochondrial genes into the nuclear genome (NUMTs). DNA sequence data also demonstrated a consistent co-infection of Wolbachia across the A. morosus distribution with every individual infected with both bacterial strains. Our data are consistent with the presence of two mitogenomes within all individuals examined in this species and suggest a major divergence from standard patterns of mitochondrial inheritance. Because the host's mitogenome and the Wolbachia genome are genetically linked through maternal inheritance, we propose three possible hypotheses that could explain maintenance of the widespread and conserved co-occurring bacterial and mitochondrial genomes in this species.

Next-generation sequencing reveals the mitogenomic heteroplasmy in the topmouth culter (Culter alburnus Basilewsky, 1855)

BACKGROUND

The mitogenomic heteroplasmy is the presence of multiple haplotypes in the mitochondria, which could cause genetic diseases and is also associated with many critical biological functions. The topmouth culter (Culter alburnus Basilewsky, 1855) is one of the most important freshwater fish in the family of Cyprinidae in China. At present, there are no reports on the topmouth culter's mtDNA heteroplasmy and the existence of which is not known.

METHODS AND RESULTS

This study aimed to analyze the mitogenomic heteroplasmy in the topmouth culter by the next-generation sequencing of the fins' total DNA. The results confirmed the existence of the heteroplasmy and indicated the presence of the extensive heteroplasmy in the topmouth culter's mitogenome. There were 38 heteroplasmic variations in the protein-coding genes from the three specimens, with 33 non-synonymous substitutions accounting for 86.84% and five synonymous substitutions accounting for 13.16%. Among them, the ND6 had the most heteroplasmic variations but only one synonymous substitution. After removing the putative nuclear mitochondrial DNA fragments, the ratio of primary haplotype in the three specimens was 43.89%, 74.72%, and 32.76%, respectively. The three specimens contained 21, 7, and 21 haplotypes of the mitogenomes, respectively. Due to the extensive heteroplasmy, we reconstructed the phylogenetic tree of the topmouth culter using the RY-coding method, which improved the performance of the phylogenetic tree to some extent.

CONCLUSIONS

This study reported the mitogenomic heteroplasmy in the topmouth culter and enhanced the knowledge regarding the mitogenomic heteroplasmy in phylogenetic studies. As the topmouth culter is a commercial species, the mitogenomic heteroplasmy is crucial for the fisheries management of the topmouth culter.

Genus-Wide Characterization of Nuclear Mitochondrial DNAs in Bumblebee (Hymenoptera: Apidae) Genomes

Simple Summary

The DNA of mitochondria can be transferred into the nucleus and form nuclear mitochondrial DNAs (NUMTs). In this study, we identified and characterized NUMTs in genus-wide bumblebee species. The number of NUMTs in bumblebee is much lower than those in its closely related taxon, honeybee. The insertion sites of NUMTs in bumblebee are not random, with AT-rich regions harboring more NUMTs. In addition, NUMTs derived from the mitochondrial COX1 gene are most abundant in the nuclear genome. While the majority of NUMTs seem unfunctional in the bumblebee, some NUMTs show functional clues, which could fuse with their flanking sequences to form novel proteins. Our results shed light on the molecular features of NUMTs and uncover their contribution to genome innovation in the bumblebee.

Abstract

In eukaryotes, DNA of mitochondria is transferred into the nucleus and forms nuclear mitochondrial DNAs (NUMTs). Taking advantage of the abundant genomic resources for bumblebees, in this study, we de novo generated mitochondrial genomes (mitogenomes) for 11 bumblebee species. Then, we identified and characterized NUMTs in genus-wide bumblebee species. The number of identified NUMTs varies across those species, with numbers ranging from 32 to 72, and nuclear genome size is not positively related to NUMT number. The insertion sites of NUMTs in the nuclear genome are not random, with AT-rich regions harboring more NUMTs. In addition, our results suggest that NUMTs derived from the mitochondrial COX1 gene are most abundant in the bumblebee nuclear genome. Although the majority of NUMTs are found within intergenic regions, some NUMTs do reside within genic regions. Transcripts that contain both the NUMT sequence and its flanking non-NUMT sequences could be found in the bumblebee transcriptome, suggesting a potential domestication of NUMTs in the bumblebee. Taken together, our results shed light on the molecular features of NUMTs in the bumblebee and uncover their contribution to genome innovation.