Uneven Missing Data Skew Phylogenomic Relationships within the Lories and Lorikeets

The effect of missing data on evolutionary analysis of sequence capture bycatch, with application to an agricultural pest

Sequence capture is a genomic technique that selectively enriches target sequences before high throughput next-generation sequencing, to generate specific sequences of interest. Off-target or 'bycatch' data are often discarded from capture experiments, but can be leveraged to address evolutionary questions under some circumstances. Here, we investigated the effects of missing data on a variety of evolutionary analyses using bycatch from an exon capture experiment on the global pest moth, Helicoverpa armigera. We added > 200 new samples from across Australia in the form of mitogenomes obtained as bycatch from targeted sequence capture, and combined these into an additional larger dataset to total > 1000 mitochondrial cytochrome c oxidase subunit I (COI) sequences across the species' global distribution. Using discriminant analysis of principal components and Bayesian coalescent analyses, we showed that mitogenomes assembled from bycatch with up to 75% missing data were able to return evolutionary inferences consistent with higher coverage datasets and the broader literature surrounding H. armigera. For example, low-coverage sequences broadly supported the delineation of two H. armigera subspecies and also provided new insights into the potential for geographic turnover among these subspecies. However, we also identified key effects of dataset coverage and composition on our results. Thus, low-coverage bycatch data can offer valuable information for population genetic and phylodynamic analyses, but caution is required to ensure the reduced information does not introduce confounding factors, such as sampling biases, that drive inference. We encourage more researchers to consider maximizing the potential of the targeted sequence approach by examining evolutionary questions with their off-target bycatch where possible-especially in cases where no previous mitochondrial data exists-but recommend stratifying data at different genome coverage thresholds to separate sampling effects from genuine genomic signals, and to understand their implications for evolutionary research.

Relicts in the mist: Two new frog families, genera and species highlight the role of Pantepui as a biodiversity museum throughout the Cenozoic

The iconic mountains of the Pantepui biogeographical region host many early-diverging endemic animal and plant lineages, concurring with Conan Doyle's novel about an ancient "Lost World". While this is the case of several frog lineages, others appear to have more recent origins, adding to the controversy around the diversification processes in this region. Due to its remoteness, Pantepui is challenging for biological surveys, and only a glimpse of its biodiversity has been described, which hampers comprehensive evolutionary studies in many groups. During a recent expedition to the Neblina massif on the Brazil-Venezuela border, we sampled two new frog species that could not be assigned to any known genus. Here, we perform phylogenetic analyses of mitogenomic and nuclear loci to infer the evolutionary relationships of the new taxa and support their description. We find that both species represent single lineages deeply nested within Brachycephaloidea, a major Neotropical clade of direct-developing frogs. Both species diverged >45 Ma from their closest relatives: the first is sister to all other Brachycephaloidea except for Ceuthomantis, another Pantepui endemic, and the second is sister to Brachycephalidae, endemic to the Brazilian Atlantic Forest. In addition to these considerable phylogenetic and biogeographic divergences, external morphology and osteological features support the proposition of two new family and genus-level taxa to accommodate these new branches of the amphibian tree of life. These findings add to other recently described ancient vertebrate lineages from the Neblina massif, providing a bewildering reminder that our perception of the Pantepui's biodiversity remains vastly incomplete. It also provides insights into how these mountains acted as "museums" during the diversification of Brachycephaloidea and of Neotropical biotas more broadly, in line with the influential "Plateau theory". Finally, these discoveries point at the yet unknown branches of the tree of life that may go extinct, due to global climate change and zoonotic diseases, before we even learn about their existence, amphibians living at higher elevations being particularly at risk.

Incongruence in the phylogenomics era

Genome-scale data and the development of novel statistical phylogenetic approaches have greatly aided the reconstruction of a broad sketch of the tree of life and resolved many of its branches. However, incongruence - the inference of conflicting evolutionary histories - remains pervasive in phylogenomic data, hampering our ability to reconstruct and interpret the tree of life. Biological factors, such as incomplete lineage sorting, horizontal gene transfer, hybridization, introgression, recombination and convergent molecular evolution, can lead to gene phylogenies that differ from the species tree. In addition, analytical factors, including stochastic, systematic and treatment errors, can drive incongruence. Here, we review these factors, discuss methodological advances to identify and handle incongruence, and highlight avenues for future research.

An evaluation of DNA extraction methods on historical and roadkill mammalian specimen

Guidelines identifying appropriate DNA extraction methods for both museum and modern biological samples are scarce or non-existent for mammalian species. Yet, obtaining large-scale genetic material collections are vital for conservation and management purposes. In this study, we evaluated five protocols making use of either spin-column, organic solvents, or magnetic bead-based methods for DNA extraction on skin samples from both modern, traffic-killed (n = 10) and museum (n = 10) samples of European hedgehogs, Ericaneus europaeus. We showed that phenol-chloroform or silica column (NucleoSpin Tissue) protocols yielded the highest amount of DNA with satisfactory purity compared with magnetic bead-based protocols, especially for museum samples. Furthermore, extractions using the silica column protocol appeared to produce longer DNA fragments on average than the other methods tested. Our investigation demonstrates that both commercial extraction kits and phenol-chloroform protocol retrieve acceptable DNA concentrations for downstream processes, from degraded remnants of traffic-killed and museum samples of mammalian specimens. Although all the tested methods could be applied depending on the research questions and laboratory conditions, commercial extraction kits may be preferred due to their effectiveness, safety and the higher quality of the DNA extractions.

Appendage-Bearing Sordariomycetes from Dipterocarpus alatus Leaf Litter in Thailand

Leaf litter is an essential functional aspect of forest ecosystems, acting as a source of organic matter, a protective layer in forest soils, and a nurturing habitat for micro- and macro-organisms. Through their successional occurrence, litter-inhabiting microfungi play a key role in litter decomposition and nutrient recycling. Despite their importance in terrestrial ecosystems and their abundance and diversity, information on the taxonomy, diversity, and host preference of these decomposer taxa is scarce. This study aims to clarify the taxonomy and phylogeny of four saprobic fungal taxa inhabiting Dipterocarpus alatus leaf litter. Leaf litter samples were collected from Doi Inthanon National Park in Chiang Mai, northern Thailand. Fungal isolates were characterized based on morphology and molecular phylogeny of the nuclear ribosomal DNA (ITS, LSU) and protein-coding genes (tub2, tef1-α, rpb2). One novel saprobic species, Ciliochorella dipterocarpi, and two new host records, Pestalotiopsis dracontomelon and Robillarda australiana, are introduced. The newly described taxa are compared with similar species, and comprehensive descriptions, micrographs, and phylogenetic trees are provided.

Museomics and phylogenomics of lovebirds (Psittaciformes, Psittaculidae, Agapornis) using low-coverage whole-genome sequencing

Natural history collections contain specimens that provide important insights into studies of ecology and evolution. With the advancement of high-throughput sequencing, historical DNA (hDNA) from museum specimens has become a valuable source of genomic data to study the evolutionary history of organisms. Low-coverage whole genome sequencing (WGS) has been increasingly applied to museum specimens for analyzing organelle genomes, but is still uncommon for genotyping the nuclear DNA fraction. In this study, we applied low-coverage WGS to phylogenomic analyses of parrots in the genus Agapornis by including both modern samples and historical specimens of ∼100-year-old. Agapornis are small-sized African and Malagasy parrots with diverse characters. Earlier phylogenetic studies failed to resolve the positions of some key lineages, prohibiting a robust interpretation of the biogeography and evolution of these African parrots. Here, we demonstrated the use of low-coverage WGS for generating both mitochondrial and nuclear genomic data, and evaluated data quality differences between modern and historical samples. Our resolved Agapornis phylogeny indicates the ancestor of Agapornis likely colonized Madagascar from Australasia by trans-oceanic dispersal events before dispersing to the African continent. Genome-wide SNPs also allowed us to identify the parental origins of hybrid Agapornis individuals. This study demonstrates the potential of applying low-coverage WGS to phylogenomics and population genomics analyses and illustrates how including historical museum specimens can address outstanding questions regarding the evolutionary history of contemporary lineages.

Phylogeny of sea spiders (Arthropoda: Pycnogonida) inferred from mitochondrial genome and 18S ribosomal RNA gene sequences

The phylogeny of sea spiders has been debated for more than a century. Despite several molecular studies in the last twenty years, interfamilial relationships remain uncertain. In the present study, relationships within Pycnogonida are examined in the light of a new dataset composed of 160 mitochondrial genomes (including 152 new sequences) and 130 18S rRNA gene sequences (including 120 new sequences), from 141 sea spider morphospecies representing 26 genera and 9 families. Node congruence between mitochondrial and nuclear markers was analysed to identify the most reliable relationships. We also reanalysed a multilocus dataset previously published and showed that the high percentages of missing data make phylogenetic conclusions difficult and uncertain. Our results support the monophyly of most families currently accepted, except Callipallenidae and Nymphonidae, the monophyly of the superfamilies Ammotheoidea (Ammotheidae + Pallenopsidae), Nymphonoidea (Nymphonidae + Callipallenidae), Phoxichilidioidea (Phoxichilidiidae + Endeidae) and Colossendeoidea (Colossendeidae + Pycnogonidae + Rhynchothoracidae), and the sister-group relationship between Ammotheoidea and Phoxichilidioidea. We discuss the morphological evolution of sea spiders, identifying homoplastic characters and possible synapomorphies. We also discuss the palaeontological and phylogenetic arguments supporting either a radiation of sea spiders prior to Jurassic or a progressive diversification from Ordovician or Cambrian.

Phylogenomic analysis of the parrots of the world distinguishes artifactual from biological sources of gene tree discordance

Gene tree discordance is expected in phylogenomic trees and biological processes are often invoked to explain it. However, heterogeneous levels of phylogenetic signal among individuals within datasets may cause artifactual sources of topological discordance. We examined how the information content in tips and subclades impacts topological discordance in the parrots (Order: Psittaciformes), a diverse and highly threatened clade of nearly 400 species. Using ultraconserved elements from 96% of the clade's species-level diversity, we estimated concatenated and species trees for 382 ingroup taxa. We found that discordance among tree topologies was most common at nodes dating between the late Miocene and Pliocene, and often at the taxonomic level of genus. Accordingly, we used two metrics to characterize information content in tips and assess the degree to which conflict between trees was being driven by lower quality samples. Most instances of topological conflict and non-monophyletic genera in the species tree could be objectively identified using these metrics. For subclades still discordant after tip-based filtering, we used a machine learning approach to determine whether phylogenetic signal or noise was the more important predictor of metrics supporting the alternative topologies. We found that when signal favored one of the topologies, noise was the most important variable in poorly performing models that favored the alternative topology. In sum, we show that artifactual sources of gene tree discordance, which are likely a common phenomenon in many datasets, can be distinguished from biological sources by quantifying the information content in each tip and modeling which factors support each topology.

Historical specimens and the limits of subspecies phylogenomics in the New World quails (Odontophoridae)

As phylogenomics focuses on comprehensive taxon sampling at the species and population/subspecies levels, incorporating genomic data from historical specimens has become increasingly common. While historical samples can fill critical gaps in our understanding of the evolutionary history of diverse groups, they also introduce additional sources of phylogenomic uncertainty, making it difficult to discern novel evolutionary relationships from artifacts caused by sample quality issues. These problems highlight the need for improved strategies to disentangle artifactual patterns from true biological signal as historical specimens become more prevalent in phylogenomic datasets. Here, we tested the limits of historical specimen-driven phylogenomics to resolve subspecies-level relationships within a highly polytypic family, the New World quails (Odontophoridae), using thousands of ultraconserved elements (UCEs). We found that relationships at and above the species-level were well-resolved and highly supported across all analyses, with the exception of discordant relationships within the two most polytypic genera which included many historical specimens. We examined the causes of discordance and found that inferring phylogenies from subsets of taxa resolved the disagreements, suggesting that analyzing subclades can help remove artifactual causes of discordance in datasets that include historical samples. At the subspecies-level, we found well-resolved geographic structure within the two most polytypic genera, including the most polytypic species in this family, Northern Bobwhites (Colinus virginianus), demonstrating that variable sites within UCEs are capable of resolving phylogenetic structure below the species level. Our results highlight the importance of complete taxonomic sampling for resolving relationships among polytypic species, often through the inclusion of historical specimens, and we propose an integrative strategy for understanding and addressing the uncertainty that historical samples sometimes introduce to phylogenetic analyses.

Evaluation of Arabian Vascular Plant Barcodes (rbcL and matK): Precision of Unsupervised and Supervised Learning Methods towards Accurate Identification

Arabia is the largest peninsula in the world, with >3000 species of vascular plants. Not much effort has been made to generate a multi-locus marker barcode library to identify and discriminate the recorded plant species. This study aimed to determine the reliability of the available Arabian plant barcodes (>1500; rbcL and matK) at the public repository (NCBI GenBank) using the unsupervised and supervised methods. Comparative analysis was carried out with the standard dataset (FINBOL) to assess the methods and markers' reliability. Our analysis suggests that from the unsupervised method, TaxonDNA's All Species Barcode criterion (ASB) exhibits the highest accuracy for rbcL barcodes, followed by the matK barcodes using the aligned dataset (FINBOL). However, for the Arabian plant barcode dataset (GBMA), the supervised method performed better than the unsupervised method, where the Random Forest and K-Nearest Neighbor (gappy kernel) classifiers were robust enough. These classifiers successfully recognized true species from both barcode markers belonging to the aligned and alignment-free datasets, respectively. The multi-class classifier showed high species resolution following the two classifiers, though its performance declined when employed to recognize true species. Similar results were observed for the FINBOL dataset through the supervised learning approach; overall, matK marker showed higher accuracy than rbcL. However, the lower rate of species identification in matK in GBMA data could be due to the higher evolutionary rate or gaps and missing data, as observed for the ASB criterion in the FINBOL dataset. Further, a lower number of sequences and singletons could also affect the rate of species resolution, as observed in the GBMA dataset. The GBMA dataset lacks sufficient species membership. We would encourage the taxonomists from the Arabian Peninsula to join our campaign on the Arabian Barcode of Life at the Barcode of Life Data (BOLD) systems. Our efforts together could help improve the rate of species identification for the Arabian Vascular plants.

A phylogenomic perspective on the evolutionary history of the stonefly genus Suwallia (Plecoptera: Chloroperlidae) revealed by ultraconserved genomic elements

Evolutionary biologists have long sought to disentangle phylogenetic relationships among taxa spanning the tree of life, an increasingly important task as anthropogenic influences accelerate population declines and species extinctions, particularly in insects. Phylogenetic analyses are commonly used to identify unique evolutionary lineages, to clarify taxonomic designations of the focal taxa, and to inform conservation decisions. Advances in DNA sequencing techniques have increasingly facilitated the ability of researchers to apply genomic methods to phylogenetic analyses, even for non-model organisms. Stoneflies are non-model insects that are important bioindicators of the quality of freshwater habitats and landscape disturbance as they spend the immature stages of their life cycles in fresh water, and the adult stages in terrestrial environments. Phylogenetic relationships within the stonefly genus Suwallia (Insecta: Plecoptera: Chloroperlidae) are poorly understood, and have never been assessed using molecular data. We used DNA sequence data from genome-wide ultraconserved element loci to generate the first molecular phylogeny for the group and assess its monophyly. We found that Palearctic and Nearctic Suwallia do not form reciprocally monophyletic clades, and that a biogeographic history including dispersal, vicariance, and founder event speciation via jump dispersal best explains the geographic distribution of this group. Our results also strongly suggest that Neaviperla forcipata (Neave, 1929) is nested within Suwallia, and the concept of the genus Suwallia should be revised to include it. Thus, we formally propose a new taxonomic combination wherein Neaviperla forcipata (Neave, 1929) is reclassified as Suwallia forcipata (Neave, 1929). Moreover, some Suwallia species (e.g., S. amoenacolens, S. kerzhneri, S. marginata, S. pallidula, and S. starki) exhibit pronounced cryptic diversity that is worthy of further investigation. These findings provide a first glimpse into the evolutionary history of Suwallia, improve our understanding of stonefly diversity in the tribe Suwallini, and highlight areas where additional research is needed.

Mining museums for historical DNA: advances and challenges in museomics

Historical DNA (hDNA), obtained from museum and herbarium specimens, has yielded spectacular new insights into the history of organisms. This includes documenting historical genetic erosion and extinction, discovering species new to science, resolving evolutionary relationships, investigating epigenetic effects, and determining origins of infectious diseases. However, the development of best-practices in isolating, processing, and analyzing hDNA remain under-explored, due to the substantial diversity of specimen preparation types, tissue sources, archival ages, and collecting histories. Thus, for hDNA to reach its full potential, and justify the destructive sampling of the rarest specimens, more experimental work using time-series collections, and the development of improved methods to correct for data asymmetries and biases due to DNA degradation are required.

Genomic library preparation and hybridization capture of formalin-fixed tissues and allozyme supernatant for population genomics and considerations for combining capture- and RADseq-based single nucleotide polymorphism data sets

Until recently many historical museum specimens were largely inaccessible to genomic inquiry, but high-throughput sequencing (HTS) approaches have allowed researchers to successfully sequence genomic DNA from dried and fluid-preserved museum specimens. In addition to preserved specimens, many museums contain large series of allozyme supernatant samples, but the amenability of these samples to HTS has not yet been assessed. Here, we compared the performance of a target-capture approach using alternative sources of genomic DNA from 10 specimens of spring salamanders (Plethodontidae: Gyrinophilus porphyriticus) collected between 1985 and 1990: allozyme supernatants, allozyme homogenate pellets and formalin-fixed tissues. We designed capture probes based on double-digest restriction-site associated sequencing (RADseq) derived loci from frozen blood samples available for seven of the specimens and assessed the success and consistency of capture and RADseq approaches. This study design enabled direct comparisons of data quality and potential biases among the different data sets for phylogenomic and population genomic analyses. We found that in phylogenetic analyses, all enrichment types for a given specimen clustered together. In principal component space all capture-based samples clustered together, but RADseq samples did not cluster with corresponding capture-based samples. Single nucleotide polymorphism calls were on average 18.3% different between enrichment types for a given individual, but these discrepancies were primarily due to differences in heterozygous/homozygous single nucleotide polymorphism calls. We demonstrate that both allozyme supernatant and formalin-fixed samples can be successfully used for population genomic analyses and we discuss ways to identify and reduce biases associated with combining capture and RADseq data.