Deep ancestral introgression shapes evolutionary history of dragonflies and damselflies

Ultrafast learning of four-node hybridization cycles in phylogenetic networks using algebraic invariants

Motivation

The abundance of gene flow in the Tree of Life challenges the notion that evolution can be represented with a fully bifurcating process which cannot capture important biological realities like hybridization, introgression, or horizontal gene transfer. Coalescent-based network methods are increasingly popular, yet not scalable for big data, because they need to perform a heuristic search in the space of networks as well as numerical optimization that can be NP-hard. Here, we introduce a novel method to reconstruct phylogenetic networks based on algebraic invariants. While there is a long tradition of using algebraic invariants in phylogenetics, our work is the first to define phylogenetic invariants on concordance factors (frequencies of four-taxon splits in the input gene trees) to identify level-1 phylogenetic networks under the multispecies coalescent model.

Results

Our novel hybrid detection methodology is optimization-free as it only requires the evaluation of polynomial equations, and as such, it bypasses the traversal of network space, yielding a computational speed at least 10 times faster than the fastest-to-date network methods. We illustrate our method's performance on simulated and real data from the genus Canis.

Availability and implementation

We present an open-source publicly available Julia package PhyloDiamond.jl available at https://github.com/solislemuslab/PhyloDiamond.jl with broad applicability within the evolutionary community.

Integrated phylogenomic approaches in insect systematics

The increased accessibility of genomic and imaging methods, and the improved access to ecological, spatial, and other natural history-related data is allowing for insect systematics to grow and find answers to central evolutionary and taxonomic questions. Today, integrated studies in insect phylogenomics and systematics are combining natural history, behavior, developmental biology, morphology, fossils, geographic range data, and ecological interactions. This integration is contributing to the clarification of evolutionary relationships, and the recognition of the role played by these factors on the evolution of insects. Future work should continue to build on these advances, seeking to further increase open-access databasing and support for natural history research, as well as expand its analytical palettes.

Exploring chromosome evolution in 250 million year old groups of dragonflies and damselflies (Insecta:Odonata)

Using recently published chromosome-length genome assemblies of two damselfly species, Ischnura elegans and Platycnemis pennipes, and two dragonfly species, Pantala flavescens and Tanypteryx hageni, we demonstrate that the autosomes of Odonata have undergone few fission, fusion, or inversion events, despite 250 million years of separation. In the four genomes discussed here, our results show that all autosomes have a clear ortholog in the ancestral karyotype. Despite this clear chromosomal orthology, we demonstrate that different factors, including concentration of repeat dynamics, GC content, relative position on the chromosome, and the relative proportion of coding sequence all influence the density of syntenic blocks across chromosomes. However, these factors do not interact to influence synteny the same way in any two pairs of species, nor is any one factor retained in all four species. Furthermore, it was previously unknown whether the micro-chromosomes in Odonata are descended from one ancestral chromosome. Despite structural rearrangements, our evidence suggests that the micro-chromosomes in the sampled Odonata do indeed descend from an ancestral chromosome, and that the micro-chromosome in P. flavescens was lost through fusion with autosomes.

PhyloMatcher: a tool for resolving conflicts in taxonomic nomenclature

Summary

Large-scale comparative studies rely on the application of both phylogenetic trees and phenotypic data, both of which come from a variety of sources, but due to the changing nature of phylogenetic classification over time, many taxon names in comparative datasets do not match the nomenclature in phylogenetic trees. Manual curation of taxonomic synonyms in large comparative datasets can be daunting. To address this issue, we introduce PhyloMatcher, a tool which allows for programmatic querying of the National Center for Biotechnology Information Taxonomy and Global Biodiversity Information Facility databases to find associated synonyms with given target species names.

Availability and implementation

PhyloMatcher is easily installed as a Python package with pip, or as a standalone GUI application. PhyloMatcher source code and documentation are freely available at https://github.com/Lswhiteh/PhyloMatcher, the GUI application can be downloaded from the Releases page.

Deciphering complex reticulate evolution of Asian Buddleja (Scrophulariaceae): insights into the taxonomy and speciation of polyploid taxa in the Sino-Himalayan region

BACKGROUND AND AIMS

Species of the genus Buddleja in Asia are mainly distributed in the Sino-Himalayan region and form a challenging taxonomic group, with extensive hybridization and polyploidization. A phylogenetic approach to unravelling the history of reticulation in this lineage will deepen our understanding of the speciation in biodiversity hotspots.

METHODS

For this study, we obtained 80 accessions representing all the species in the Asian Buddleja clade, and the ploidy level of each taxon was determined by flow cytometry analyses. Whole plastid genomes, nuclear ribosomal DNA, single nucleotide polymorphisms and a large number of low-copy nuclear genes assembled from genome skimming data were used to investigate the reticulate evolutionary history of Asian Buddleja. Complex cytonuclear conflicts were detected through a comparison of plastid and species trees. Gene tree incongruence was also analysed to detect any reticulate events in the history of this lineage.

KEY RESULTS

Six hybridization events were detected, which are able to explain the cytonuclear conflict in Asian Buddleja. Furthermore, PhyloNet analysis combining species ploidy data indicated several allopolyploid speciation events. A strongly supported species tree inferred from a large number of low-copy nuclear genes not only corrected some earlier misinterpretations, but also indicated that there are many Asian Buddleja species that have been lumped mistakenly. Divergent time estimation shows two periods of rapid diversification (8-10 and 0-3 Mya) in the Asian Buddleja clade, which might coincide with the final uplift of the Hengduan Mountains and Quaternary climate fluctuations, respectively.

CONCLUSIONS

This study presents a well-supported phylogenetic backbone for the Asian Buddleja species, elucidates their complex and reticulate evolutionary history and suggests that tectonic activity, climate fluctuations, polyploidization and hybridization together promoted the diversification of this lineage.

PhyloMatcher: a tool for resolving conflicts in taxonomic nomenclature

Summary

Large-scale comparative studies rely on the application of both phylogenetic trees and phenotypic data, both of which come from a variety of sources, but due to the changing nature of phylogenetic classification over time, many taxon names in comparative datasets do not match the nomenclature in phylogenetic trees. Manual curation of taxonomic synonyms in large comparative datasets can be daunting. To address this issue, we introduce PhyloMatcher, a tool which allows for programmatic querying of two commonly used taxonomic databases to find associated synonyms with given target species names.

Availability and implementation

PhyloMatcher is easily installed as a Python package with pip, or as a standalone GUI application. PhyloMatcher source code and documentation are freely available at https://github.com/Lswhiteh/PhyloMatcher, the GUI application can be downloaded from the Releases page.

Contact

Lswhiteh@unc.edu.

Supplemental Information

We provide documentation for PhyloMatcher, including walkthrough instructions for the GUI application on the Releases page of https://github.com/Lswhiteh/PhyloMatcher.

The damselfly genus Furagrion Petrulevičius et al. (Odonata, Zygoptera) from the early Eocene Fur Formation of Denmark and the dysagrionoid grade

The earliest Eocene odonate genus Furagrion Petrulevičius et al. from the Danish Fur Formation is revised based on eighteen specimens, two of which apparently have been lost since their publication. The holotype of Phenacolestes jutlandicus Henriksen, type species of Furagrion, is incomplete and lacks the characters currently used to differentiate species, genera and higher taxa in Odonata. We, therefore, propose that the holotype is set aside and a recently discovered nearly complete Fur Formation fossil is designated as neotype. Furagrion possesses all of the nine wing character states currently used along with head shape for diagnosing the Dysagrionidae; however, Furagrion has a characteristically zygopteran head, not the distinctive head shape of the suborder Cephalozygoptera. We, therefore, treat it as a zygopteran unassigned to family. These nine wing character states appear in different combinations not only in various Zygoptera and Cephalozygoptera, but also in the Frenguelliidae, an Eocene family of Argentina that may represent an unnamed suborder. We recognise these taxa as constituting a dysagrionoid grade, in which these character states appear either convergently or as symplesiomorphies. Furagrion morsi Zessin is synonymized with Phenacolestes jutlandicus Henriksen, syn. nov. and Morsagrion Zessin with Furagrion Petrulevičius, Wappler, Wedmann, Rust, and Nel, syn. nov.

Phylogeny and historical biogeography of the water boatmen (Insecta: Hemiptera: Heteroptera: Nepomorpha: Corixoidea)

The water boatmen of Corixoidea, a group of aquatic bugs with more than 600 extant species, is one of the largest superfamilies of Nepomorpha. Contrary to the other nepomorphan lineages, the Corixoidea are most diverse in the Laurasian remnant Holarctic region. To explicitly test whether the present-day Holarctic distribution of diverse corixids is associated with the arising of the Laurasian landmass that was separated from Gondwana, we investigated the phylogeny, divergence times and historical biogeography of Corixoidea based on morphological and molecular characters sampled from 122 taxa representing all families, subfamilies, tribes and approximately 54 % of the genera. Our results were largely congruent with the phylogenetic relationships within the established nepomorphan phylogenetic context. The fossil calibrated chronogram, diversification analysis and ancestral ranges reconstruction indicated that Corixoidea began to diversify in Gondwana in the late Triassic approximately at 224 Ma and the arising of the most diverse subfamily Corixinae in Corixidae in the Holarctic region was largely congruent with the time of separation of Laurasia from Gondwana. The large-scale expansion of the temperate and cold zones on the northward-moving Laurasian landmass after the breakup of the Pangea provided new aquatic niches and ecological opportunities for promoting rapid diversification for the Holarctic corixid lineage.

Species diversity, hotspot congruence, and conservation of North American damselflies (Odonata: Zygoptera)

The rapid extinction of species is of considerable concern for biodiversity conservation. Identifying the drivers of species diversity and hotspots is beneficial for developing conservation strategies. Studies on insects have mainly focused on terrestrial species and rarely on semiaquatic species. Using 135,208 georeferenced occurrence records of 296 damselflies across North America, their species richness and endemism (represented by weighted endemism) patterns were mapped in a 100 × 100-km grid size, and the effects of environmental variables on species richness and endemism were investigated using generalized linear models and hierarchical partitioning. Subsequently, the top 5% grids with species richness and weighted endemism were separately selected as hotspots and their congruence was evaluated. Finally, species diversity hotspots were identified by integrating two types of hotspot grids, and gap analysis was performed to evaluate their conservation status. Temperature conditions and water availability had the strongest influence on species richness and endemism, respectively. Low congruence among species richness and endemism hotspots was observed. Moreover, four species diversity hotspots were identified, namely, region of the eastern United States and southeastern Canada, southwestern United States, central Mexico, and southernmost North America. Approximately 69.31% of the hotspot grids are not a part of the existing protected areas, presenting a significant conservation gap. The habitats of taxonomic groups should be considered while identifying the most common driving mechanisms of endemism. Strengthening the establishment of protected areas in regions with conservation gaps is urgently needed to promote the conservation of damselflies in North America.

Evolutionary history and divergence times of Odonata (dragonflies and damselflies) revealed through transcriptomics

Dragonflies and damselflies are among the earliest flying insects with extant representatives. However, unraveling details of their long evolutionary history, such as egg laying (oviposition) strategies, is impeded by unresolved phylogenetic relationships, particularly in damselflies. Here we present a transcriptome-based phylogenetic reconstruction of Odonata, analyzing 2,980 protein-coding genes in 105 species representing nearly all the order's families. All damselfly and most dragonfly families are recovered as monophyletic. Our data suggest a sister relationship between dragonfly families of Gomphidae and Petaluridae. According to our divergence time estimates, both crown-Zygoptera and -Anisoptera arose during the late Triassic. Egg-laying with a reduced ovipositor apparently evolved in dragonflies during the late Jurassic/early Cretaceous. Lastly, we also test the impact of fossil choice and placement, particularly, of the extinct fossil species, †Triassolestodes asiaticus, and †Proterogomphus renateae on divergence time estimates. We find placement of †Proterogomphus renateae to be much more impactful than †Triassolestodes asiaticus.