Anchored phylogenomics unravels the evolution of spider flies (Diptera, Acroceridae) and reveals discordance between nucleotides and amino acids

Resolving Marine-Freshwater Transitions by Diatoms Through a Fog of Gene Tree Discordance

Despite the obstacles facing marine colonists, most lineages of aquatic organisms have colonized and diversified in freshwaters repeatedly. These transitions can trigger rapid morphological or physiological change and, on longer timescales, lead to increased rates of speciation and extinction. Diatoms are a lineage of ancestrally marine microalgae that have diversified throughout freshwater habitats worldwide. We generated a phylogenomic data set of genomes and transcriptomes for 59 diatom taxa to resolve freshwater transitions in one lineage, the Thalassiosirales. Although most parts of the species tree were consistently resolved with strong support, we had difficulties resolving a Paleocene radiation, which affected the placement of one freshwater lineage. This and other parts of the tree were characterized by high levels of gene tree discordance caused by incomplete lineage sorting and low phylogenetic signal. Despite differences in species trees inferred from concatenation versus summary methods and codons versus amino acids, traditional methods of ancestral state reconstruction supported six transitions into freshwaters, two of which led to subsequent species diversification. Evidence from gene trees, protein alignments, and diatom life history together suggest that habitat transitions were largely the product of homoplasy rather than hemiplasy, a condition where transitions occur on branches in gene trees not shared with the species tree. Nevertheless, we identified a set of putatively hemiplasious genes, many of which have been associated with shifts to low salinity, indicating that hemiplasy played a small but potentially important role in freshwater adaptation. Accounting for differences in evolutionary outcomes, in which some taxa became locked into freshwaters while others were able to return to the ocean or become salinity generalists, might help further distinguish different sources of adaptive mutation in freshwater diatoms.

Evolutionary history of interactions among terrestrial arthropods

The study of terrestrial arthopod fossils preserved with microscopic fidelity in amber and as permineralized replicas has been revolutionized by CT scanning. Fine preservation facilitates phylogenetic interpretation of fossils, but molecular divergence-time models still commonly use insufficient fossil calibrations, skewing estimates away from the direct (i.e. fossil, morphological) evidence. Interactions among terrestrial arthropods (predation, parasitoidism; phoresy, social symbionts) are briefly reviewed from the fossil record. Predation is the oldest and most widespread, originating with arachnids since probably the Silurian. The first phoretic arthropods were probably mites (Acari). Parasitoidism extends to the early Jurassic ~200 mya, with four main episodes proposed by [1•]. 100-myo Burmese amber, the most diverse Cretaceous paleobiota, is unique for our understanding of insect eusociality and interrelationships among terrestrial arthropods. Eusocial insect colonies are ecological sinks for thousands of symbiont species; ages of the major eusocial groups and some of their nest symbionts are discussed. Fossilized arthropod interrelationships in Miocene Dominican amber are presented as visual exemplars.

Host-trailing satellite flight behaviour is associated with greater investment in peripheral visual sensory system in miltogrammine flies

Insect sensory systems are the subjects of different selective pressures that shape their morphology. In many species of the flesh fly subfamily Miltogramminae (Diptera: Sarcophagidae) that are kleptoparasitic on bees and wasps, females perch on objects close to the host nests and, once a returning host is detected, they follow it in flight at a fixed distance behind until reaching the nest. We hypothesized that such satellite (SAT) flight behaviour, which implies a finely coordinated trailing flight, is associated with an improved visual system, compared to species adopting other, non-satellite (NON-SAT) strategies. After looking at body size and common ancestry, we found that SAT species have a greater number of ommatidia and a greater eye surface area when compared to NON-SAT species. Ommatidium area is only affected by body size, suggesting that selection changes disproportionately (relative to body size variation) the number of ommatidia and as a consequence the eye area, instead of ommatidium size. SAT species also tend to have larger ocelli, but their role in host-finding was less clear. This suggests that SAT species may have a higher visual acuity by increasing ommatidia number, as well as better stability during flight and motion perception through larger ocelli. Interestingly, antennal length was significantly reduced in SAT species, and ommatidia number negatively correlated with antennal length. While this finding does not imply a selection pressure of improved antennal sensory system in species adopting NON-SAT strategies, it suggests an inverse resource (i.e. a single imaginal disc) allocation between eyes and antennae in this fly subfamily.

Review of Chilean Cyrtinae (Diptera: Acroceridae) with the Description of Three New Species and the First Record of Villalus inanis from Argentina

Acrocerid flies are endoparasitoids of spiders that are distributed in all ecosystems except deserts and the poles. This group was historically classified into three subfamilies. However, molecular analyses have indicated two new subfamilies. One of these is Cyrtinae Newman, which is made up of two clades: the north, with Palearctic species, and the south, with Chilean species (except for the fossil genus Villalites Hennig). In the southern clade, endemic genera Holops Philippi and Villalus Cole have been poorly described over time by different authors. The aim of this study is to re-describe the known species of Holops and Villalus and describe three new Holops species based on collected material: Holops grezi sp. nov., Holops anarayae sp. nov., and Holops pullomen sp. nov. In addition, an identification species key with updated distributions for all species and the first record of Villalus inanis from Argentina is also provided. New species are mainly characterized by the length of flagellum; colors of the postocellar lobe and thoracic segments, pilosity and legs; wing cells and veins, and shape of the abdomen. Holops and Villalus do not share morphological traits with any other known extant genus of Acroceridae in the southern hemisphere. The real biodiversity of this group must be evaluated with a molecular phylogeny in the future. In addition, the new species in the Pehuén and Valdivian forest provinces, areas considered to have been sheltered during the last glacial period, promoting speciation in several taxa. New and cryptic Holops species open the doors to disentangling the real biodiversity of spider flies in Chile.

Whole-Genome Phylogenetic Reconstruction as a Powerful Tool to Reveal Homoplasy and Ancient Rapid Radiation in Waterflea Evolution

Although phylogeny estimation is notoriously difficult in radiations that occurred several hundred million years ago, phylogenomic approaches offer new ways to examine relationships among ancient lineages and evaluate hypotheses that are key to evolutionary biology. Here, we reconstruct the deep-rooted relationships of one of the oldest living arthropod clades, the branchiopod crustaceans, using a kaleidoscopic approach. We use concatenation and coalescent tree-building methods to analyze a large multigene data set at the nucleotide and amino acid level and examine gene tree versus species tree discordance. We unequivocally resolve long-debated relationships among extant orders of the Cladocera, the waterfleas, an ecologically relevant zooplankton group in global aquatic and marine ecosystems that is famous for its model systems in ecology and evolution. To build the data set, we assembled eight de novo genomes of key taxa including representatives of all extant cladoceran orders and suborders. Our phylogenetic analysis focused on a BUSCO-based set of 823 conserved single-copy orthologs shared among 23 representative taxa spanning all living branchiopod orders, including 11 cladoceran families. Our analysis supports the monophyly of the Cladocera and reveals remarkable homoplasy in their body plans. We found large phylogenetic distances between lineages with similar ecological specializations, indicating independent evolution in major body plans, such as in the pelagic predatory orders Haplopoda and Onychopoda (the “Gymnomera”). In addition, we assessed rapid cladogenesis by estimating relative timings of divergence in major lineages using reliable fossil-calibrated priors on eight nodes in the branchiopod tree, suggesting a Paleozoic origin around 325 Ma for the cladoceran ancestor and an ancient rapid radiation around 252 Ma at the Perm/Triassic boundary. These findings raise new questions about the roles of homoplasy and rapid radiation in the diversification of the cladocerans and help examine trait evolution from a genomic perspective in a functionally well understood, ancient arthropod group. [Cladocera; Daphnia; evolution; homoplasy; molecular clock; phylogenomics; systematics; waterfleas.]

New insights into the phylogeny and evolution of Podocarpaceae inferred from transcriptomic data

Phylogenies of an increasing number of taxa have been resolved with the development of phylogenomics. However, the intergeneric relationships of Podocarpaceae, the second largest family of conifers comprising 19 genera and approximately 187 species mainly distributed in the Southern Hemisphere, have not been well disentangled in previous studies, even when genome-scale data sets were used. Here we used 993 nuclear orthologous groups (OGs) and 54 chloroplast OGs (genes), which were generated from 47 transcriptomes of Podocarpaceae and its sister group Araucariaceae, to reconstruct the phylogeny of Podocarpaceae. Our study completely resolved the intergeneric relationships of Podocarpaceae represented by all extant genera and revealed that topological conflicts among phylogenetic trees could be attributed to synonymous substitutions. Moreover, we found that two morphological traits, fleshy seed cones and flattened leaves, might be important for Podocarpaceae to adapt to angiosperm-dominated forests and thus could have promoted its species diversification. In addition, our results indicate that Podocarpaceae originated in Gondwana in the late Triassic and both vicariance and dispersal have contributed to its current biogeographic patterns. Our study provides the first robust transcriptome-based phylogeny of Podocarpaceae, an evolutionary framework important for future studies of this family.

Phylogenomics reveals accelerated late Cretaceous diversification of bee flies (Diptera: Bombyliidae)

Bombyliidae is a very species-rich and widespread family of parasitoid flies with more than 250 genera classified into 17 extant subfamilies. However, little is known about their evolutionary history or how their present-day diversity was shaped. Transcriptomes of 15 species and anchored hybrid enrichment (AHE) sequence captures of 86 species, representing 94 bee fly species and 14 subfamilies, were used to reconstruct the phylogeny of Bombyliidae. We integrated data from transcriptomes across each of the main lineages in our AHE tree to build a data set with more genes (550 loci versus 216 loci) and higher support levels. Our overall results show strong congruence with the current classification of the family, with 11 out of 14 included subfamilies recovered as monophyletic. Heterotropinae and Mythicomyiinae are successive sister groups to the remainder of the family. We examined the evolution of key morphological characters through our phylogenetic hypotheses and show that neither the "sand chamber subfamilies" nor the "Tomophthalmae" are monophyletic in our phylogenomic analyses. Based on our results, we reinstate two tribes at the subfamily level (Phthiriinae stat. rev. and Ecliminae stat. rev.) and we include the genus Sericosoma Macquart (previously incertae sedis) in the subfamily Oniromyiinae, bringing the total number of bee fly subfamilies to 19. Our dating analyses indicate a Jurassic origin of the family (165-194 Ma), with the sand chamber evolving early in bee fly evolution, in the late Jurassic or mid-Cretaceous (100-165 Ma). We hypothesize that the angiosperm radiation and the hothouse climate established during the late Cretaceous accelerated the diversification of bee flies, by providing an expanded range of resources for the parasitoid larvae and nectarivorous adults.

From e-voucher to genomic data: Preserving archive specimens as demonstrated with medically important mosquitoes (Diptera: Culicidae) and kissing bugs (Hemiptera: Reduviidae)

Scientific collections such as the U.S. National Museum (USNM) are critical to filling knowledge gaps in molecular systematics studies. The global taxonomic impediment has resulted in a reduction of expert taxonomists generating new collections of rare or understudied taxa and these large historic collections may be the only reliable source of material for some taxa. Integrated systematics studies using both morphological examinations and DNA sequencing are often required for resolving many taxonomic issues but as DNA methods often require partial or complete destruction of a sample, there are many factors to consider before implementing destructive sampling of specimens within scientific collections. We present a methodology for the use of archive specimens that includes two crucial phases: 1) thoroughly documenting specimens destined for destructive sampling—a process called electronic vouchering, and 2) the pipeline used for whole genome sequencing of archived specimens, from extraction of genomic DNA to assembly of putative genomes with basic annotation. The process is presented for eleven specimens from two different insect subfamilies of medical importance to humans: Anophelinae (Diptera: Culicidae)—mosquitoes and Triatominae (Hemiptera: Reduviidae)—kissing bugs. Assembly of whole mitochondrial genome sequences of all 11 specimens along with the results of an ortholog search and BLAST against the NCBI nucleotide database are also presented.

New insights into the phylogeny and evolution of lady beetles (Coleoptera: Coccinellidae) by extensive sampling of genes and species

Ladybirds (family Coccinellidae) are one of the most diverse groups of beetles and globally comprise over 6000 species. Despite their scientific and economic significance, the taxonomy of Coccinellidae remains unstable, and we still know little about their evolutionary history. By using a small number of genes, previous phylogenetic analyses have not reliably resolved the relationships among major ladybird lineages. In this study, we sequenced 94 nuclear protein-coding genes for 214 species of Coccinellidae and 14 outgroups, covering 90 genera and 35 tribes. We found that nucleotide compositional heterogeneity is present among ladybird tribes so that phylogenetic inference at the amino acid level is more reliable than at the DNA level. Based on the maximum likelihood analyses of the amino acid dataset, we recognize three subfamilies in Coccinellidae: Microweiseinae, Monocoryninae stat. nov., and Coccinellinae. The subfamily relationships are strongly supported as (Microweiseinae, (Monocoryninae stat. nov., Coccinellinae)). The tribes of ladybirds are mostly monophyletic, except Ortaliini, Sticholotidini, Scymnini, and Coccidulini. The phylogenetic relationships among tribes of Coccinellinae are still not well resolved, with many nodes weakly supported. Our divergence time analysis suggests that the crown group of extant lady beetles arose in the Early Cretaceous ~ 143 million years ago (Mya) and experienced a rapid diversification during the Late Cretaceous (120-70 Mya). We hypothesize that the boom of angiosperms in the Late Cretaceous promoted the diversification of herbivorous Sternorrhyncha insects, especially aphids, which in turn drove the rapid radiation of predatory lady beetles. In summary, our work provides a comprehensive time-calibrated phylogeny of Coccinellidae that provides a sound framework for revising their classification and understanding the origin of their biodiversity.

Data Types and the Phylogeny of Neoaves

The phylogeny of Neoaves, the largest clade of extant birds, has remained unclear despite intense study. The difficulty associated with resolving the early branches in Neoaves is likely driven by the rapid radiation of this group. However, conflicts among studies may be exacerbated by the data type analyzed. For example, analyses of coding exons typically yield trees that place Strisores (nightjars and allies) sister to the remaining Neoaves, while analyses of non-coding data typically yield trees where Mirandornites (flamingos and grebes) is the sister of the remaining Neoaves. Our understanding of data type effects is hampered by the fact that previous analyses have used different taxa, loci, and types of non-coding data. Herein, we provide strong corroboration of the data type effects hypothesis for Neoaves by comparing trees based on coding and non-coding data derived from the same taxa and gene regions. A simple analytical method known to minimize biases due to base composition (coding nucleotides as purines and pyrimidines) resulted in coding exon data with increased congruence to the non-coding topology using concatenated analyses. These results improve our understanding of the resolution of neoavian phylogeny and point to a challenge—data type effects—that is likely to be an important factor in phylogenetic analyses of birds (and many other taxonomic groups). Using our results, we provide a summary phylogeny that identifies well-corroborated relationships and highlights specific nodes where future efforts should focus.

Phylogenies from mitochondrial genomes of 120 species of ticks: Insights into the evolution of the families of ticks and of the genus Amblyomma

The evolution and phylogenetic relationships of the ticks at both the family and genus levels are contested. The genus Amblyomma and its subgenera are in a state of flux; moreover, the relationships among the three tick families are controversial due to conflicting phylogenetic support for different arrangements of the three families of living ticks. With 18 newly sequenced mitochondrial (mt) genomes of ticks included, we executed the largest mt genome phylogenetic study of ticks so far. Phylogenetic trees were inferred from one sea spider mt genome, one horseshoe crab, five mite mt genomes and 146 tick mt genomes from 120 species: 153 mt genomes in total. Sixteen phylogenetic trees were inferred from 10 datasets using both maximum likelihood and Bayesian inference methods. We describe the first novel mt gene-arrangement for the metastriate Ixodidae in Amblyomma (Africaniella) transversale. Also, three unusual partial 16S rRNA gene inserts were found in the mt genome of Haemaphysalis (Alloceraea) kitaokai: we consider the possible role of past genome translocation events in the formation of these inserts. Our phylogenies revealed evidence that: (i) the genus Amblyomma is polyphyletic with respect to Amblyomma (Africaniella) transversale; (ii) the subgenus Aponomma is apparently embedded in the genus Amblyomma; (iii) Haemaphysalis (Segalia) parva and Haemaphysalis (Alloceraea) kitaokai form a clade to the exclusion of other Haemaphysalis species; and (iv) the phylogenetic position of the family Nuttalliellidae is unstable among phylogenies from different datasets.

An integrative phylogenomic approach to elucidate the evolutionary history and divergence times of Neuropterida (Insecta: Holometabola)

BACKGROUND

The latest advancements in DNA sequencing technologies have facilitated the resolution of the phylogeny of insects, yet parts of the tree of Holometabola remain unresolved. The phylogeny of Neuropterida has been extensively studied, but no strong consensus exists concerning the phylogenetic relationships within the order Neuroptera. Here, we assembled a novel transcriptomic dataset to address previously unresolved issues in the phylogeny of Neuropterida and to infer divergence times within the group. We tested the robustness of our phylogenetic estimates by comparing summary coalescent and concatenation-based phylogenetic approaches and by employing different quartet-based measures of phylogenomic incongruence, combined with data permutations.

RESULTS

Our results suggest that the order Raphidioptera is sister to Neuroptera + Megaloptera. Coniopterygidae is inferred as sister to all remaining neuropteran families suggesting that larval cryptonephry could be a ground plan feature of Neuroptera. A clade that includes Nevrorthidae, Osmylidae, and Sisyridae (i.e. Osmyloidea) is inferred as sister to all other Neuroptera except Coniopterygidae, and Dilaridae is placed as sister to all remaining neuropteran families. Ithonidae is inferred as the sister group of monophyletic Myrmeleontiformia. The phylogenetic affinities of Chrysopidae and Hemerobiidae were dependent on the data type analyzed, and quartet-based analyses showed only weak support for the placement of Hemerobiidae as sister to Ithonidae + Myrmeleontiformia. Our molecular dating analyses suggest that most families of Neuropterida started to diversify in the Jurassic and our ancestral character state reconstructions suggest a primarily terrestrial environment of the larvae of Neuropterida and Neuroptera.

CONCLUSION

Our extensive phylogenomic analyses consolidate several key aspects in the backbone phylogeny of Neuropterida, such as the basal placement of Coniopterygidae within Neuroptera and the monophyly of Osmyloidea. Furthermore, they provide new insights into the timing of diversification of Neuropterida. Despite the vast amount of analyzed molecular data, we found that certain nodes in the tree of Neuroptera are not robustly resolved. Therefore, we emphasize the importance of integrating the results of morphological analyses with those of sequence-based phylogenomics. We also suggest that comparative analyses of genomic meta-characters should be incorporated into future phylogenomic studies of Neuropterida.

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.

Putting the genome in insect phylogenomics

Next-generation sequencing technologies provide a substantial increase in the size of molecular phylogenetic datasets that can be obtained for studies of insect systematics. Several new genome reduction approaches are leveraging these technologies to generate large phylogenomic datasets: targeted amplicon sequencing, target capture, and transcriptome sequencing. Although cost effective, these approaches provide limited data for questions outside of phylogenetics. For many groups of insects, sequencing the entire genome at modest coverage is feasible. Using these genomic reads, an automated Target Restricted Assembly Method (aTRAM) can use the results of blast searches to assemble thousands of single copy ortholog genes across a group of interest. These locally assembled genes can then be compiled into very large phylogenomic datasets. These genomic libraries have the advantage in that they also contain reads from the mitochondrial genome and symbiont genomes, as well the entire insect genome, and can be leveraged for additional studies beyond phylogenetics.

Whole Genome Shotgun Phylogenomics Resolves the Pattern and Timing of Swallowtail Butterfly Evolution

Evolutionary relationships have remained unresolved in many well-studied groups, even though advances in next-generation sequencing and analysis, using approaches such as transcriptomics, anchored hybrid enrichment, or ultraconserved elements, have brought systematics to the brink of whole genome phylogenomics. Recently, it has become possible to sequence the entire genomes of numerous nonbiological models in parallel at reasonable cost, particularly with shotgun sequencing. Here, we identify orthologous coding sequences from whole-genome shotgun sequences, which we then use to investigate the relevance and power of phylogenomic relationship inference and time-calibrated tree estimation. We study an iconic group of butterflies—swallowtails of the family Papilionidae—that has remained phylogenetically unresolved, with continued debate about the timing of their diversification. Low-coverage whole genomes were obtained using Illumina shotgun sequencing for all genera. Genome assembly coupled to BLAST-based orthology searches allowed extraction of 6621 orthologous protein-coding genes for 45 Papilionidae species and 16 outgroup species (with 32% missing data after cleaning phases). Supermatrix phylogenomic analyses were performed with both maximum-likelihood (IQ-TREE) and Bayesian mixture models (PhyloBayes) for amino acid sequences, which produced a fully resolved phylogeny providing new insights into controversial relationships. Species tree reconstruction from gene trees was performed with ASTRAL and SuperTriplets and recovered the same phylogeny. We estimated gene site concordant factors to complement traditional node-support measures, which strengthens the robustness of inferred phylogenies. Bayesian estimates of divergence times based on a reduced data set (760 orthologs and 12% missing data) indicate a mid-Cretaceous origin of Papilionoidea around 99.2 Ma (95% credibility interval: 68.6–142.7 Ma) and Papilionidae around 71.4 Ma (49.8–103.6 Ma), with subsequent diversification of modern lineages well after the Cretaceous-Paleogene event. These results show that shotgun sequencing of whole genomes, even when highly fragmented, represents a powerful approach to phylogenomics and molecular dating in a group that has previously been refractory to resolution.

Gene Selection and Evolutionary Modeling Affect Phylogenomic Inference of Neuropterida Based on Transcriptome Data

Neuropterida is a super order of Holometabola that consists of the orders Megaloptera (dobsonflies, fishflies, and alderflies), Neuroptera (lacewings) and Raphidioptera (snakeflies). Several proposed higher-level relationships within Neuropterida, such as the relationships between the orders or between the families, have been extensively debated. To further understand the evolutionary history of Neuropterida, we conducted phylogenomic analyses of all 13 published transcriptomes of the neuropterid species, as well as of a new transcriptome of the fishfly species Ctenochauliodes similis of Liu and Yang, 2006 (Megaloptera: Corydalidae: Chauliodinae) that we sequenced. Our phylogenomic data matrix contained 1392 ortholog genes from 22 holometabolan species representing six families from Neuroptera, two families from Raphidioptera, and two families from Megaloptera as the ingroup taxa, and nine orders of Holometabola as outgroups. Phylogenetic reconstruction was performed using both concatenation and coalescent-based approaches under a site-homogeneous model as well as under a site-heterogeneous model. Surprisingly, analyses using the site-homogeneous model strongly supported a paraphyletic Neuroptera, with Coniopterygidae assigned as the sister group of all other Neuropterida. In contrast, analyses using the site-heterogeneous model recovered Neuroptera as monophyletic. The monophyly of Neuroptera was also recovered in concatenation and coalescent-based analyses using genes with stronger phylogenetic signals [i.e., higher average bootstrap support (ABS) values and higher relative tree certainty including all conflicting bipartitions (RTCA) values] under the site-homogeneous model. The present study illustrated how both data selection and model selection influence phylogenomic analyses of large-scale data matrices comprehensively.