A novel assembly pipeline and functional annotations for targeted sequencing: A case study on the globally threatened Margaritiferidae (Bivalvia: Unionida)

The proliferation of genomic sequencing approaches has significantly impacted the field of phylogenetics. Target capture approaches provide a cost-effective, fast and easily applied strategy for phylogenetic inference of non-model organisms. However, several existing target capture processing pipelines are incapable of incorporating whole genome sequencing (WGS). Here, we develop a new pipeline for capture and de novo assembly of the targeted regions using whole genome re-sequencing reads. This new pipeline captured targeted loci accurately, and given its unbiased nature, can be used with any target capture probe set. Moreover, due to its low computational demand, this new pipeline may be ideal for users with limited resources and when high-coverage sequencing outputs are required. We demonstrate the utility of our approach by incorporating WGS data into the first comprehensive phylogenomic reconstruction of the freshwater mussel family Margaritiferidae. We also provide a catalogue of well-curated functional annotations of these previously uncharacterized freshwater mussel-specific target regions, representing a complementary tool for scrutinizing phylogenetic inferences while expanding future applications of the probe set.

Corrigendum: Phylogenomics of the Andean tetraploid clade of the American Amaryllidaceae (subfamily Amaryllidoideae): Unlocking a polyploid generic radiation abetted by continental geodynamics

[This corrects the article DOI: 10.3389/fpls.2020.582422.].

Sage Insights Into the Phylogeny of Salvia: Dealing With Sources of Discordance Within and Across Genomes

Next-generation sequencing technologies have facilitated new phylogenomic approaches to help clarify previously intractable relationships while simultaneously highlighting the pervasive nature of incongruence within and among genomes that can complicate definitive taxonomic conclusions. Salvia L., with ∼1,000 species, makes up nearly 15% of the species diversity in the mint family and has attracted great interest from biologists across subdisciplines. Despite the great progress that has been achieved in discerning the placement of Salvia within Lamiaceae and in clarifying its infrageneric relationships through plastid, nuclear ribosomal, and nuclear single-copy genes, the incomplete resolution has left open major questions regarding the phylogenetic relationships among and within the subgenera, as well as to what extent the infrageneric relationships differ across genomes. We expanded a previously published anchored hybrid enrichment dataset of 35 exemplars of Salvia to 179 terminals. We also reconstructed nearly complete plastomes for these samples from off-target reads. We used these data to examine the concordance and discordance among the nuclear loci and between the nuclear and plastid genomes in detail, elucidating both broad-scale and species-level relationships within Salvia. We found that despite the widespread gene tree discordance, nuclear phylogenies reconstructed using concatenated, coalescent, and network-based approaches recover a common backbone topology. Moreover, all subgenera, except for Audibertia, are strongly supported as monophyletic in all analyses. The plastome genealogy is largely resolved and is congruent with the nuclear backbone. However, multiple analyses suggest that incomplete lineage sorting does not fully explain the gene tree discordance. Instead, horizontal gene flow has been important in both the deep and more recent history of Salvia. Our results provide a robust species tree of Salvia across phylogenetic scales and genomes. Future comparative analyses in the genus will need to account for the impacts of hybridization/introgression and incomplete lineage sorting in topology and divergence time estimation.

A rodent anchored hybrid enrichment probe set for a range of phylogenetic utility: From order to species

Rodents are the largest order of mammals and contain several model organisms important to scientific research in a variety of fields, yet no large set of genomic markers have been designed for this group to date, hindering evolutionary studies into relationships of the group as a whole. Here we present a genomic probe set designed and optimized for rodents with a protocol that is easy to replicate with little laboratory investment. This design utilizes an anchored hybrid enrichment approach specifically targeting rodents to generate longer loci with a higher substitution rate than existing vertebrate probes to provide utility at various taxonomic levels. Using a test set of rodents from all five suborders, we successfully obtained alignments for 416 of the 418 target loci with an average of 1379 bp per locus and a total alignment of more than half a million base pairs. This genomic data set performed well in all phylogenetic analyses, especially in recent phylogenetic splits, with ample parsimony-informative sites within genera and even within species, showing more than four times as many single nucleotide polymorphisms per locus than a recent vertebrate ultraconserved elements study. Additional support is provided in resolving deeper clades in Rodentia. By providing this probe design, we hope that more laboratories can easily generate data for answering questions in rodents from species delimitation to understanding relationships among families in rapid radiations.

Identification of Potential Host Plants of Sap-Sucking Insects (Hemiptera: Cicadellidae) Using Anchored Hybrid By-Catch Data

Reliable host plant records are available for only a small fraction of herbivorous insect species, despite their potential agricultural importance. Most available data on insect-plant associations have been obtained through field observations of occurrences of insects on particular plants. Molecular methods have more recently been used to identify potential host plants using DNA extracted from insects, but most prior studies using these methods have focused on chewing insects that ingest tissues expected to contain large quantities of plant DNA. Screening of Illumina data obtained from sap feeders of the hemipteran family Cicadellidae (leafhoppers) using anchored hybrid enrichment indicates that, despite feeding on plant fluids, these insects often contain detectable quantities of plant DNA. Although inclusion of probes for bacterial 16S in the original anchored hybrid probe kit yielded relatively high detection rates for chloroplast 16S, the Illumina short reads also, in some cases, included DNA for various plant barcode genes as "by-catch". Detection rates were generally only slightly higher for Typhlocybinae, which feed preferentially on parenchyma cell contents, compared to other groups of leafhoppers that feed preferentially on phloem or xylem. These results indicate that next-generation sequencing provides a powerful tool to investigate the specific association between individual insect and plant species.

Food Plant Shifts Drive the Diversification of Sack-Bearer Moths

AbstractLepidoptera are a highly diverse group of herbivorous insects; however, some superfamilies have relatively few species. Two alternative hypotheses for drivers of Lepidoptera diversity are shifts in food plant use or shifts from concealed to external feeding as larvae. Many studies address the former hypothesis but with bias toward externally feeding taxa. One of the most striking examples of species disparity between sister lineages in Lepidoptera is between the concealed-feeding sack-bearer moths (Mimallonoidea), which contain about 300 species, and externally feeding Macroheterocera, which have over 74,000 species. We provide the first dated tree of Mimallonidae to understand the diversification dynamics of these moths in order to fill a knowledge gap pertaining to drivers of diversity within an important concealed-feeding clade. We find that Mimallonidae is an ancient Lepidoptera lineage that originated in the Cretaceous ∼105 million years ago and has had a close association with the plant order Myrtales for the past 40 million years. Diversification dynamics are tightly linked with food plant usage in this group. Reliance on Myrtales may have influenced diversification of Mimallonidae because clades that shifted away from the ancestral condition of feeding on Myrtales have the highest speciation rates in the family.

Anchored Phylogenomics, Evolution and Systematics of Elateridae: Are All Bioluminescent Elateroidea Derived Click Beetles?

Simple Summary

In the era of phylogenomics, new molecular sequencing and computational techniques can aid in resolving phylogenetic relationships that were previously intractable by morphological or limited molecular data. In this study, we used anchored hybrid enrichment—designed to recover DNA sequences from hundreds of single-copy orthologous genes—to resolve the phylogeny of the Elateridae (click-beetles) and establish their placement within superfamily Elateroidea. The resulting data were compatible with published transcriptomes, allowing for integrating our dataset with previously published data. Using a wide range of analyses on these molecular data, we tested hypotheses long-debated in the morphological literature and also the robustness of our phylogenetic inferences. Our results placed the bioluminescent lampyroids (fireflies and relatives) within the click-beetles, challenging the current classification of Elateridae, Lampyridae, Phengodidae, and Rhagophthalmidae. However, despite the large amount of molecular data analyzed, a few nodes with conflicting phylogenetic signals could not be unambiguously resolved. Overall, we recovered well-resolved tree topologies that will serve as a framework for further systematic and evolutionary studies of click-beetles. This work further demonstrates that the click-beetle lineage contains not only pest wireworms, but also many species that benefit agriculture.

Abstract

Click-beetles (Coleoptera: Elateridae) are an abundant, diverse, and economically important beetle family that includes bioluminescent species. To date, molecular phylogenies have sampled relatively few taxa and genes, incompletely resolving subfamily level relationships. We present a novel probe set for anchored hybrid enrichment of 2260 single-copy orthologous genes in Elateroidea. Using these probes, we undertook the largest phylogenomic study of Elateroidea to date (99 Elateroidea, including 86 Elateridae, plus 5 non-elateroid outgroups). We sequenced specimens from 88 taxa to test the monophyly of families, subfamilies and tribes. Maximum likelihood and coalescent phylogenetic analyses produced well-resolved topologies. Notably, the included non-elaterid bioluminescent families (Lampyridae + Phengodidae + Rhagophthalmidae) form a clade within the otherwise monophyletic Elateridae, and Sinopyrophoridae may not warrant recognition as a family. All analyses recovered the elaterid subfamilies Elaterinae, Agrypninae, Cardiophorinae, Negastriinae, Pityobiinae, and Tetralobinae as monophyletic. Our results were conflicting on whether the hypnoidines are sister to Dendrometrinae or Cardiophorinae + Negastriinae. Moreover, we show that fossils with the eucnemid-type frons and elongate cylindrical shape may belong to Eucnemidae, Elateridae: Thylacosterninae, ancestral hard-bodied cantharoids or related extinct groups. Proposed taxonomic changes include recognition of Plastocerini as a tribe in Dendrometrinae and Hypnoidinae stat. nov. as a subfamily within Elateridae.

Phylogenomics of the Andean Tetraploid Clade of the American Amaryllidaceae (Subfamily Amaryllidoideae): Unlocking a Polyploid Generic Radiation Abetted by Continental Geodynamics

One of the two major clades of the endemic American Amaryllidaceae subfam. Amaryllidoideae constitutes the tetraploid-derived (n = 23) Andean-centered tribes, most of which have 46 chromosomes. Despite progress in resolving phylogenetic relationships of the group with plastid and nrDNA, certain subclades were poorly resolved or weakly supported in those previous studies. Sequence capture using anchored hybrid enrichment was employed across 95 species of the clade along with five outgroups and generated sequences of 524 nuclear genes and a partial plastome. Maximum likelihood phylogenetic analyses were conducted on concatenated supermatrices, and coalescent-based species tree analyses were run on the gene trees, followed by hybridization network, age diversification and biogeographic analyses. The four tribes Clinantheae, Eucharideae, Eustephieae, and Hymenocallideae (sister to Clinantheae) are resolved in all analyses with > 90 and mostly 100% support, as are almost all genera within them. Nuclear gene supermatrix and species tree results were largely in concordance; however, some instances of cytonuclear discordance were evident. Hybridization network analysis identified significant reticulation in Clinanthus, Hymenocallis, Stenomesson and the subclade of Eucharideae comprising Eucharis, Caliphruria, and Urceolina. Our data support a previous treatment of the latter as a single genus, Urceolina, with the addition of Eucrosia dodsonii. Biogeographic analysis and penalized likelihood age estimation suggests an origin in the Cauca, Desert and Puna Neotropical bioprovinces for the complex in the mid-Oligocene, with more dispersals than vicariances in its history, but no extinctions. Hymenocallis represents the only instance of long-distance vicariance from the tropical Andean origin of its tribe Hymenocallideae. The absence of extinctions correlates with the lack of diversification rate shifts within the clade. The Eucharideae experienced a sudden lineage radiation ca. 10 Mya. We tie much of the divergences in the Andean-centered lineages to the rise of the Andes, and suggest that the Amotape-Huancabamba Zone functioned as both a corridor (dispersal) and a barrier to migration (vicariance). Several taxonomic changes are made. This is the largest DNA sequence data set to be applied within Amaryllidaceae to date.

Uncovering the genomic signature of ancient introgression between white oak lineages (Quercus)

Botanists have long recognised interspecific gene flow as a common occurrence within white oaks (Quercus section Quercus). Historical allele exchange, however, has not been fully characterised and the complex genomic signals resulting from the combination of vertical and horizontal gene transmission may confound phylogenetic inference and obscure our ability to accurately infer the deep evolutionary history of oaks. Using anchored enrichment, we obtained a phylogenomic dataset consisting of hundreds of single-copy nuclear loci. Concatenation, species-tree and network analyses were carried out in an attempt to uncover the genomic signal of ancient introgression and infer the divergent phylogenetic topology for the white oak clade. Locus and site-level likelihood comparisons were then conducted to further explore the introgressed signal within our dataset. Historical, intersectional gene flow is suggested to have occurred between an ancestor of the Eurasian Roburoid lineage and Quercus pontica and North American Dumosae and Prinoideae lineages. Despite extensive time past, our approach proved successful in detecting the genomic signature of ancient introgression. Our results, however, highlight the importance of sampling and the use of a plurality of analytical tools and methods to sufficiently explore genomic datasets, uncover this signal, and accurately infer evolutionary history.

A phylogenomic framework, evolutionary timeline and genomic resources for comparative studies of decapod crustaceans

Comprising over 15 000 living species, decapods (crabs, shrimp and lobsters) are the most instantly recognizable crustaceans, representing a considerable global food source. Although decapod systematics have received much study, limitations of morphological and Sanger sequence data have yet to produce a consensus for higher-level relationships. Here, we introduce a new anchored hybrid enrichment kit for decapod phylogenetics designed from genomic and transcriptomic sequences that we used to capture new high-throughput sequence data from 94 species, including 58 of 179 extant decapod families, and 11 of 12 major lineages. The enrichment kit yields 410 loci (greater than 86 000 bp) conserved across all lineages of Decapoda, more clade-specific molecular data than any prior study. Phylogenomic analyses recover a robust decapod tree of life strongly supporting the monophyly of all infraorders, and monophyly of each of the reptant, 'lobster' and 'crab' groups, with some results supporting pleocyemate monophyly. We show that crown decapods diverged in the Late Ordovician and most crown lineages diverged in the Triassic-Jurassic, highlighting a cryptic Palaeozoic history, and post-extinction diversification. New insights into decapod relationships provide a phylogenomic window into morphology and behaviour, and a basis to rapidly and cheaply expand sampling in this economically and ecologically significant invertebrate clade.

Optimizing Phylogenomics with Rapidly Evolving Long Exons: Comparison with Anchored Hybrid Enrichment and Ultraconserved Elements

Marker selection has emerged as an important component of phylogenomic study design due to rising concerns of the effects of gene tree estimation error, model misspecification, and data-type differences. Researchers must balance various trade-offs associated with locus length and evolutionary rate among other factors. The most commonly used reduced representation data sets for phylogenomics are ultraconserved elements (UCEs) and Anchored Hybrid Enrichment (AHE). Here, we introduce Rapidly Evolving Long Exon Capture (RELEC), a new set of loci that targets single exons that are both rapidly evolving (evolutionary rate faster than RAG1) and relatively long in length (>1,500 bp), while at the same time avoiding paralogy issues across amniotes. We compare the RELEC data set to UCEs and AHE in squamate reptiles by aligning and analyzing orthologous sequences from 17 squamate genomes, composed of 10 snakes and 7 lizards. The RELEC data set (179 loci) outperforms AHE and UCEs by maximizing per-locus genetic variation while maintaining presence and orthology across a range of evolutionary scales. RELEC markers show higher phylogenetic informativeness than UCE and AHE loci, and RELEC gene trees show greater similarity to the species tree than AHE or UCE gene trees. Furthermore, with fewer loci, RELEC remains computationally tractable for full Bayesian coalescent species tree analyses. We contrast RELEC to and discuss important aspects of comparable methods, and demonstrate how RELEC may be the most effective set of loci for resolving difficult nodes and rapid radiations. We provide several resources for capturing or extracting RELEC loci from other amniote groups.

Target Nuclear and Off-Target Plastid Hybrid Enrichment Data Inform a Range of Evolutionary Depths in the Orchid Genus Epidendrum

Universal angiosperm enrichment probe sets designed to enrich hundreds of putatively orthologous nuclear single-copy loci are increasingly being applied to infer phylogenetic relationships of different lineages of angiosperms at a range of evolutionary depths. Studies applying such probe sets have focused on testing the universality and performance of the target nuclear loci, but they have not taken advantage of off-target data from other genome compartments generated alongside the nuclear loci. Here we do so to infer phylogenetic relationships in the orchid genus Epidendrum and closely related genera of subtribe Laeliinae. Our aims are to: 1) test the technical viability of applying the plant anchored hybrid enrichment (AHE) method (Angiosperm v.1 probe kit) to our focal group, 2) mine plastid protein coding genes from off-target reads; and 3) evaluate the performance of the target nuclear and off-target plastid loci in resolving and supporting phylogenetic relationships along a range of taxonomical depths. Phylogenetic relationships were inferred from the nuclear data set through coalescent summary and site-based methods, whereas plastid loci were analyzed in a concatenated partitioned matrix under maximum likelihood. The usefulness of target and flanking non-target nuclear regions and plastid loci was assessed through the estimation of their phylogenetic informativeness. Our study successfully applied the plant AHE probe kit to Epidendrum, supporting the universality of this kit in angiosperms. Moreover, it demonstrated the feasibility of mining plastome loci from off-target reads generated with the Angiosperm v.1 probe kit to obtain additional, uniparentally inherited sequence data at no extra sequencing cost. Our analyses detected some strongly supported incongruences between nuclear and plastid data sets at shallow divergences, an indication of potential lineage sorting, hybridization, or introgression events in the group. Lastly, we found that the per site phylogenetic informativeness of the ycf1 plastid gene surpasses that of all other plastid genes and several nuclear loci, making it an excellent candidate for assessing phylogenetic relationships at medium to low taxonomic levels in orchids.

Phylogenomic Analysis of a Putative Missing Link Sparks Reinterpretation of Leech Evolution

Leeches (Hirudinida) comprise a charismatic, yet often maligned group of worms. Despite their ecological, economic, and medical importance, a general consensus on the phylogenetic relationships of major hirudinidan lineages is lacking. This absence of a consistent, robust phylogeny of early-diverging lineages has hindered our understanding of the underlying processes that enabled evolutionary diversification of this clade. Here, we used an anchored hybrid enrichment-based phylogenomic approach, capturing hundreds of loci to investigate phylogenetic relationships among major hirudinidan lineages and their closest living relatives. We recovered Branchiobdellida as sister to a clade that includes all major lineages of hirudinidans and Acanthobdella, casting doubt on the utility of Acanthobdella as a “missing link” between hirudinidans and the clitellate group formerly known as Oligochaeta. Further, our results corroborate the reciprocal monophyly of jawed and proboscis-bearing leeches. Our phylogenomic resolution of early-diverging leeches provides a useful framework for illuminating the evolution of key adaptations and host–symbiont associations that have allowed leeches to colonize a wide diversity of habitats worldwide.

Phylogenomic Evidence Overturns Current Conceptions of Social Evolution in Wasps (Vespidae)

The hypothesis that eusociality originated once in Vespidae has shaped interpretation of social evolution for decades and has driven the supposition that preimaginal morphophysiological differences between castes were absent at the outset of eusociality. Many researchers also consider casteless nest-sharing an antecedent to eusociality. Together, these ideas endorse a stepwise progression of social evolution in wasps (solitary → casteless nest-sharing → eusociality with rudimentary behavioral castes → eusociality with preimaginal caste-biasing (PCB) → morphologically differentiated castes). Here, we infer the phylogeny of Vespidae using sequence data generated via anchored hybrid enrichment from 378 loci across 136 vespid species and perform ancestral state reconstructions to test whether rudimentary and monomorphic castes characterized the initial stages of eusocial evolution. Our results reject the single origin of eusociality hypothesis, contest the supposition that eusociality emerged from a casteless nest-sharing ancestor, and suggest that eusociality in Polistinae + Vespinae began with castes having morphological differences. An abrupt appearance of castes with ontogenetically established morphophysiological differences conflicts with the current conception of stepwise social evolution and suggests that the climb up the ladder of sociality does not occur through sequential mutation. Phenotypic plasticity and standing genetic variation could explain how cooperative brood care evolved in concert with nest-sharing and how morphologically dissimilar castes arose without a rudimentary intermediate. Furthermore, PCB at the outset of eusociality implicates a subsocial route to eusociality in Polistinae + Vespinae, emphasizing the role of mother-daughter interactions and subfertility (i.e. the cost component of kin selection) in the origin of workers.

Are 100 enough? Inferring acanthomorph teleost phylogeny using Anchored Hybrid Enrichment

BACKGROUND

The past decade has witnessed remarkable progress towards resolution of the Tree of Life. However, despite the increased use of genomic scale datasets, some phylogenetic relationships remain difficult to resolve. Here we employ anchored phylogenomics to capture 107 nuclear loci in 29 species of acanthomorph teleost fishes, with 25 of these species sampled from the recently delimited clade Ovalentaria. Previous studies employing multilocus nuclear exon datasets have not been able to resolve the nodes at the base of the Ovalentaria tree with confidence. Here we test whether a phylogenomic approach will provide better support for these nodes, and if not, why this may be.

RESULTS

After using a novel method to account for paralogous loci, we estimated phylogenies with maximum likelihood and species tree methods using DNA sequence alignments of over 80,000 base pairs. Several key relationships within Ovalentaria are well resolved, including 1) the sister taxon relationship between Cichlidae and Pholidichthys, 2) a clade containing blennies, grammas, clingfishes, and jawfishes, and 3) monophyly of Atherinomorpha (topminnows, flyingfishes, and silversides). However, many nodes in the phylogeny associated with the early diversification of Ovalentaria are poorly resolved in several analyses. Through the use of rarefaction curves we show that limited phylogenetic resolution among the earliest nodes in the Ovalentaria phylogeny does not appear to be due to a deficiency of data, as average global node support ceases to increase when only 1/3rd of the sampled loci are used in analyses. Instead this lack of resolution may be driven by model misspecification as a Bayesian mixed model analysis of the amino acid dataset provided good support for parts of the base of the Ovalentaria tree.

CONCLUSIONS

Although it does not appear that the limited phylogenetic resolution among the earliest nodes in the Ovalentaria phylogeny is due to a deficiency of data, it may be that both stochastic and systematic error resulting from model misspecification play a role in the poor resolution at the base of the Ovalentaria tree as a Bayesian approach was able to resolve some of the deeper nodes, where the other methods failed.