Mitochondrial data are not suitable for resolving placental mammal phylogeny

Museum specimens shedding light on the evolutionary history and cryptic diversity of the hedgehog family Erinaceidae

The family Erinaceidae encompasses 27 extant species in two subfamilies: Erinaceinae, which includes spiny hedgehogs, and Galericinae, which comprises silky-furred gymnures and moonrats. Although they are commonly recognized by the general public, their phylogenetic history remains incompletely understood, and several species have never been included in any molecular analyses. Additionally, previous research suggested that the species diversity of Erinaceidae might be underestimated. In this study, we sequenced the mitochondrial genomes of 29 individuals representing 18 erinaceid species using 18 freshly collected tissue and 11 historical museum specimens. We also integrated previously published data for a concatenated analysis. We aimed to elucidate the evolutionary relationships within Erinaceidae, estimate divergence times, and uncover potential underestimated species diversity. Our data finely resolved intergeneric and interspecific relationships and presented the first molecular evidence for the phylogenetic position of Mesechinus wangi, Paraechinus micropus, and P. nudiventris. Our results revealed a sister relationship between Neotetracus and Neohylomys gymnures, as well as a sister relationship between Hemiechinus and Mesechinus, supporting previous hypotheses. Additionally, our findings provided a novel phylogenetic position for Paraechinus aethiopicus, placing it in a basal position within the genus. Furthermore, our study uncovered cryptic species diversity within Hylomys suillus as well as in Neotetracus sinensis, Atelerix albiventris, P. aethiopicus, and Hemiechinus auratus, most of which have been previously overlooked.

Comparative Genome Microsynteny Illuminates the Fast Evolution of Nuclear Mitochondrial Segments (NUMTs) in Mammals

The escape of DNA from mitochondria into the nuclear genome (nuclear mitochondrial DNA, NUMT) is an ongoing process. Although pervasively observed in eukaryotic genomes, their evolutionary trajectories in a mammal-wide context are poorly understood. The main challenge lies in the orthology assignment of NUMTs across species due to their fast evolution and chromosomal rearrangements over the past 200 million years. To address this issue, we systematically investigated the characteristics of NUMT insertions in 45 mammalian genomes and established a novel, synteny-based method to accurately predict orthologous NUMTs and ascertain their evolution across mammals. With a series of comparative analyses across taxa, we revealed that NUMTs may originate from nonrandom regions in mtDNA, are likely found in transposon-rich and intergenic regions, and unlikely code for functional proteins. Using our synteny-based approach, we leveraged 630 pairwise comparisons of genome-wide microsynteny and predicted the NUMT orthology relationships across 36 mammals. With the phylogenetic patterns of NUMT presence-and-absence across taxa, we constructed the ancestral state of NUMTs given the mammal tree using a coalescent method. We found support on the ancestral node of Fereuungulata within Laurasiatheria, whose subordinal relationships are still controversial. This study broadens our knowledge on NUMT insertion and evolution in mammalian genomes and highlights the merit of NUMTs as alternative genetic markers in phylogenetic inference.

Phylogenomics reveals an almost perfect polytomy among the almost ungulates (Paenungulata)

Phylogenetic studies have resolved most relationships among Eutherian Orders. However, the branching order of elephants (Proboscidea), hyraxes (Hyracoidea), and sea cows (Sirenia) (i.e., the Paenungulata) has remained uncertain since at least 1758, when Linnaeus grouped elephants and manatees into a single Order (Bruta) to the exclusion of hyraxes. Subsequent morphological, molecular, and large-scale phylogenomic datasets have reached conflicting conclusions on the branching order within Paenungulates. We use a phylogenomic dataset of alignments from 13,388 protein-coding genes across 261 Eutherian mammals to infer phylogenetic relationships within Paenungulates. We find that gene trees almost equally support the three alternative resolutions of Paenungulate relationships and that despite strong support for a Proboscidea+Hyracoidea split in the multispecies coalescent (MSC) tree, there is significant evidence for gene tree uncertainty, incomplete lineage sorting, and introgression among Proboscidea, Hyracoidea, and Sirenia. Indeed, only 8-10% of genes have statistically significant phylogenetic signal to reject the hypothesis of a Paenungulate polytomy. These data indicate little support for any resolution for the branching order Proboscidea, Hyracoidea, and Sirenia within Paenungulata and suggest that Paenungulata may be as close to a real, or at least unresolvable, polytomy as possible.

Revised phylogeny from complete mitochondrial genomes of phyllostomid bats resolves subfamilial classification

Classically, molecular phylogenetic trees of Phyllostomidae have been inferred using a combination of a few mitochondrial and nuclear markers. However, there is still uncertainty in the relationships, especially among deep clades within the family. In this study, we provide newly sequenced complete mitochondrial genomes from 26 bat species, including genomes of 23 species reported here for the first time. By carefully analysing these genomes using maximum likelihood and Bayesian methods and different ingroup and outgroup samples, partition schemes and data types, we investigated the robustness and sensitivity of our phylogenetic results. The optimal topologies were those inferred from the complete data matrix of nucleotides, with complex and highly parameterized substitution models and partition schemes. Our results show a statistically robust picture of the evolutionary relationships between phyllostomid subfamilies and clarify hitherto uncertain relationships of Lonchorhininae and Macrotinae.

The value of updating GenBank accessions for supermatrix phylogeny: The case of the New Guinean marsupial carnivore genus Myoictis

Erroneous taxonomic attributions in GenBank accessions can mislead phylogenetic inference and appear to be widespread within genera. We investigate the influence of taxonomic misattributions for reconstructing the phylogeny of three-striped dasyures, which include four recognized Myoictis species (Marsupialia: Dasyuridae) that are distributed across New Guinea and nearby islands. Molecular phylogenetic studies that have focused on dasyurids consistently resolve the interrelationships of these small carnivores, grouping M. leucura with M. wavicus, and placing M. wallacei and M. melas as successively deeper divergences from these. Two recent marsupial and mammalian supermatrix phylogenies instead favour an alternative Myoictis topology that is discordant with each of these relationships. We add new nuclear and mitochondrial sequences and employ randomized accession resampling that shows the supermatrix topologies are an artefact of several outdated taxonomic attributions in GenBank. Updating these accessions brings agreement across Myoictis phylogenies with randomly resampled accessions. We encourage authors to update GenBank taxonomic attributions and we argue that an option is needed for flagging accessions that are not demonstrably incorrect, but that provide anomalous results. This would serve both as a caution for future supermatrix construction and to highlight accessions of potentially significant biological interest for further study.

Unlocking inaccessible historical genomes preserved in formalin

Museum specimens represent an unparalleled record of historical genomic data. However, the widespread practice of formalin preservation has thus far impeded genomic analysis of a large proportion of specimens. Limited DNA sequencing from formalin-preserved specimens has yielded low genomic coverage with unpredictable success. We set out to refine sample processing methods and to identify specimen characteristics predictive of sequencing success. With a set of taxonomically diverse specimens collected between 1962 and 2006 and ranging in preservation quality, we compared the efficacy of several end-to-end whole genome sequencing workflows alongside a k-mer-based trimming-free read alignment approach to maximize mapping of endogenous sequence. We recovered complete mitochondrial genomes and up to 3× nuclear genome coverage from formalin-preserved tissues. Hot alkaline lysis coupled with phenol-chloroform extraction out-performed proteinase K digestion in recovering DNA, while library preparation method had little impact on sequencing success. The strongest predictor of DNA yield was overall specimen condition, which additively interacts with preservation conditions to accelerate DNA degradation. Here, we demonstrate a significant advance in capability beyond limited recovery of a small number of loci via PCR or target-capture sequencing. To facilitate strategic selection of suitable specimens for genomic sequencing, we present a decision-making framework that utilizes independent and nondestructive assessment criteria. Sequencing of formalin-preserved specimens will contribute to a greater understanding of temporal trends in genetic adaptation, including those associated with a changing climate. Our work enhances the value of museum collections worldwide by unlocking genomes of specimens that have been disregarded as a valid molecular resource.

Coverage and quality of DNA barcode references for Central and Northern European Odonata

Background

Dragonflies and damselflies (Odonata) are important components in biomonitoring due to their amphibiotic lifecycle and specific habitat requirements. They are charismatic and popular insects, but can be challenging to identify despite large size and often distinct coloration, especially the immature stages. DNA-based assessment tools rely on validated DNA barcode reference libraries evaluated in a supraregional context to minimize taxonomic incongruence and identification mismatches.

Methods

This study reports on findings from the analysis of the most comprehensive DNA barcode dataset for Central European Odonata to date, with 103 out of 145 recorded European species included and publicly deposited in the Barcode of Life Data System (BOLD). The complete dataset includes 697 specimens (548 adults, 108 larvae) from 274 localities in 16 countries with a geographic emphasis on Central Europe. We used BOLD to generate sequence divergence metrics and to examine the taxonomic composition of the DNA barcode clusters within the dataset and in comparison with all data on BOLD.

Results

Over 88% of the species included can be readily identified using their DNA barcodes and the reference dataset provided. Considering the complete European dataset, unambiguous identification is hampered in 12 species due to weak mitochondrial differentiation and partial haplotype sharing. However, considering the known species distributions only two groups of five species possibly co-occur, leading to an unambiguous identification of more than 95% of the analysed Odonata via DNA barcoding in real applications. The cases of small interspecific genetic distances and the observed deep intraspecific variation in Cordulia aenea (Linnaeus, 1758) are discussed in detail and the corresponding taxa in the public reference database are highlighted. They should be considered in future applications of DNA barcoding and metabarcoding and represent interesting evolutionary biological questions, which call for in depth analyses of the involved taxa throughout their distribution ranges.

Enhancing mitogenomic phylogeny and resolving the relationships of extinct megafaunal placental mammals

Mitochondrial genomes provided the first widely used sequences that were sufficiently informative to resolve relationships among animals across a wide taxonomic domain, from within species to between phyla. However, mitogenome studies supported several anomalous relationships and fell partly out of favour as sequencing multiple, independent nuclear loci proved to be highly effective. A tendency to blame mitochondrial DNA (mtDNA) has overshadowed efforts to understand and ameliorate underlying model misspecification. Here we find that influential assessments of the infidelity of mitogenome phylogenies have often been overstated, but nevertheless, substitution saturation and compositional non-stationarity substantially mislead reconstruction. We show that RY coding the mtDNA, excluding protein-coding 3rd codon sites, partitioning models based on amino acid hydrophobicity and enhanced taxon sampling improve the accuracy of mitogenomic phylogeny reconstruction for placental mammals, almost to the level of multi-gene nuclear datasets. Indeed, combined analysis of mtDNA with 3-fold longer nuclear sequence data either maintained or improved upon the nuclear support for all generally accepted clades, even those that mtDNA alone did not favour, thus indicating "hidden support". Confident mtDNA phylogeny reconstruction is especially important for understanding the evolutionary dynamics of mitochondria themselves, and for merging extinct taxa into the tree of life, with ancient DNA often only accessible as mtDNA. Our ancient mtDNA analyses lend confidence to the relationships of three extinct megafaunal taxa: glyptodonts are nested within armadillos, the South American ungulate, Macrauchenia is sister to horses and rhinoceroses, and sabre-toothed and scimitar cats are the monophyletic sister-group of modern cats.

Novel Molecular Synapomorphies Demarcate Different Main Groups/Subgroups of Plasmodium and Piroplasmida Species Clarifying Their Evolutionary Relationships

The class Hematozoa encompasses several clinically important genera, including Plasmodium, whose members cause the major life-threating disease malaria. Hence, a good understanding of the interrelationships of organisms from this class and reliable means for distinguishing them are of much importance. This study reports comprehensive phylogenetic and comparative analyses on protein sequences on the genomes of 28 hematozoa species to understand their interrelationships. In addition to phylogenetic trees based on two large datasets of protein sequences, detailed comparative analyses were carried out on the genomes of hematozoa species to identify novel molecular synapomorphies consisting of conserved signature indels (CSIs) in protein sequences. These studies have identified 79 CSIs that are exclusively present in specific groups of Hematozoa/Plasmodium species, also supported by phylogenetic analysis, providing reliable means for the identification of these species groups and understanding their interrelationships. Of these CSIs, six CSIs are specifically shared by all hematozoa species, two CSIs serve to distinguish members of the order Piroplasmida, five CSIs are uniquely found in all Piroplasmida species except B. microti and two CSIs are specific for the genus Theileria. Additionally, we also describe 23 CSIs that are exclusively present in all genome-sequenced Plasmodium species and two, nine, ten and eight CSIs which are specific for members of the Plasmodium subgenera Haemamoeba, Laverania, Vinckeia and Plasmodium (excluding P. ovale and P. malariae), respectively. Additionally, our work has identified several CSIs that support species relationships which are not evident from phylogenetic analysis. Of these CSIs, one CSI supports the ancestral nature of the avian-Plasmodium species in comparison to the mammalian-infecting groups of Plasmodium species, four CSIs strongly support a specific relationship of species between the subgenera Plasmodium and Vinckeia and three CSIs each that reliably group P. malariae with members of the subgenus Plasmodium and P. ovale within the subgenus Vinckeia, respectively. These results provide a reliable framework for understanding the evolutionary relationships among the Plasmodium/Piroplasmida species. Further, in view of the exclusivity of the described molecular markers for the indicated groups of hematozoa species, particularly large numbers of unique characteristics that are specific for all Plasmodium species, they provide important molecular tools for biochemical/genetic studies and for developing novel diagnostics and therapeutics for these organisms.

A Novel Approach to Clustering Genome Sequences Using Inter-nucleotide Covariance

Classification of DNA sequences is an important issue in the bioinformatics study, yet most existing methods for phylogenetic analysis including Multiple Sequence Alignment (MSA) are time-consuming and computationally expensive. The alignment-free methods are popular nowadays, whereas the manual intervention in those methods usually decreases the accuracy. Also, the interactions among nucleotides are neglected in most methods. Here we propose a new Accumulated Natural Vector (ANV) method which represents each DNA sequence by a point in ℝ¹⁸. By calculating the Accumulated Indicator Functions of nucleotides, we can further find an Accumulated Natural Vector for each sequence. This new Accumulated Natural Vector not only can capture the distribution of each nucleotide, but also provide the covariance among nucleotides. Thus global comparison of DNA sequences or genomes can be done easily in ℝ¹⁸. The tests of ANV of datasets of different sizes and types have proved the accuracy and time-efficiency of the new proposed ANV method.

Highly Resolved Phylogenetic Relationships within Order Acipenseriformes According to Novel Nuclear Markers

Order Acipenseriformes contains 27 extant species distributed across the northern hemisphere, including so-called "living fossil" species of garfish and sturgeons. Previous studies have focused on their mitochondrial genetics and have rarely used nuclear genetic data, leaving questions as to their phylogenetic relationships. This study aimed to utilize a bioinformatics approach to screen for candidate single-copy nuclear genes, using transcriptomic data from sturgeon species and genomic data from the spotted gar, Lepisosteus oculatus. We utilized nested polymerase chain reaction (PCR) and degenerate primers to identify nuclear protein-coding (NPC) gene markers to determine phylogenetic relationships among the Acipenseriformes. We identified 193 nuclear single-copy genes, selected from 1850 candidate genes with at least one exon larger than 700 bp. Forty-three of these genes were used for primer design and development of 30 NPC markers, which were sequenced for at least 14 Acipenseriformes species. Twenty-seven NPC markers were found completely in 16 species. Gene trees according to Bayesian inference (BI) and maximum likelihood (ML) were calculated based on the 30 NPC markers (20,946 bp total). Both gene and species trees produced very similar topologies. A molecular clock model estimated the divergence time between sturgeon and paddlefish at 204.1 Mya, approximately 10% later than previous estimates based on cytochrome b data (184.4 Mya). The successful development and application of NPC markers provides a new perspective and insight for the phylogenetic relationships of Acipenseriformes. Furthermore, the newly developed nuclear markers may be useful in further studies on the conservation, evolution, and genomic biology of this group.