Orthograph: a versatile tool for mapping coding nucleotide sequences to clusters of orthologous genes

Evolution across the adaptive landscape in a hyperdiverse beetle radiation

The extraordinary diversification of beetles on Earth is a textbook example of adaptive evolution. Yet, the tempo and drivers of this super-radiation remain largely unclear. Here, we address this problem by investigating macroevolutionary dynamics in darkling beetles (Coleoptera: Tenebrionidae), one of the most ecomorphologically diverse beetle families (with over 30,000 species). Using multiple genomic datasets and analytical approaches, we resolve the long-standing inconsistency over deep relationships in the family. In conjunction with a landmark-based dataset of body shape morphology, we show that the evolutionary history of darkling beetles is marked by ancient rapid radiations, frequent ecological transitions, and rapid bursts of morphological diversification. On a global scale, our analyses uncovered a significant pulse of phenotypic diversification proximal to the Cretaceous-Palaeogene (K/Pg) mass extinction and convergence of body shape associated with recurrent ecological specializations. On a regional scale, two major Australasian radiations, the Adeliini and the Heleine clade, exhibited contrasting patterns of ecomorphological diversification, representing phylogenetic niche conservatism versus adaptive radiation. Our findings align with the Simpsonian model of adaptive evolution across the macroevolutionary landscape and highlight a significant role of ecological opportunity in driving the immense ecomorphological diversity in a hyperdiverse beetle group.

PhyloAln: A Convenient Reference-Based Tool to Align Sequences and High-Throughput Reads for Phylogeny and Evolution in the Omic Era

The current trend in phylogenetic and evolutionary analyses predominantly relies on omic data. However, prior to core analyses, traditional methods typically involve intricate and time-consuming procedures, including assembly from high-throughput reads, decontamination, gene prediction, homology search, orthology assignment, multiple sequence alignment, and matrix trimming. Such processes significantly impede the efficiency of research when dealing with extensive data sets. In this study, we develop PhyloAln, a convenient reference-based tool capable of directly aligning high-throughput reads or complete sequences with existing alignments as a reference for phylogenetic and evolutionary analyses. Through testing with simulated data sets of species spanning the tree of life, PhyloAln demonstrates consistently robust performance compared with other reference-based tools across different data types, sequencing technologies, coverages, and species, with percent completeness and identity at least 50 percentage points higher in the alignments. Additionally, we validate the efficacy of PhyloAln in removing a minimum of 90% foreign and 70% cross-contamination issues, which are prevalent in sequencing data but often overlooked by other tools. Moreover, we showcase the broad applicability of PhyloAln by generating alignments (completeness mostly larger than 80%, identity larger than 90%) and reconstructing robust phylogenies using real data sets of transcriptomes of ladybird beetles, plastid genes of peppers, or ultraconserved elements of turtles. With these advantages, PhyloAln is expected to facilitate phylogenetic and evolutionary analyses in the omic era. The tool is accessible at https://github.com/huangyh45/PhyloAln.

Phylogenomics resolves long-standing questions about the affinities of an endangered Corsican endemic fly

Recent studies on oestroidean Diptera (Brachycera) are providing a comprehensive and nuanced understanding of the evolutionary history of this remarkably diverse clade of holometabolous insects. The Oestroidea, which includes formidable pests such as various blowflies, botflies, and flesh flies that infest livestock, pets and humans, are mostly composed of beneficial species that act as scavengers or parasitoids on various pest insects. In our research, we used genomic methods to elucidate the phylogenetic position of Nesodexia corsicana Villeneuve, 1911 (Diptera: Calliphoridae), a mysterious oestroid species endemic to Corsica and characterized by distinctive morphological features that have puzzled taxonomists for years. Contrary to initial hypotheses, our results place Nesodexia Villeneuve, 1911 within the Calliphoridae subfamily Rhinophorinae, a small lineage of terrestrial isopod parasitoids. Through detailed morphological analysis of adults of both sexes and eggs, we uncovered significant insights consistent with our phylogenomic reconstruction. The unique morphological features of the species, coupled with its restricted and fragmented habitat, highlight its potential conservation importance. We delineated the area of occupancy for N. corsicana and assessed its "threatened" category using specific IUCN Red List criteria. In addition, we mapped the available habitat within its range and determined potential key biodiversity areas (KBA) triggered by N. corsicana. New potential KBAs are only partially covered by the Corsican Regional Park. Finally, we mapped the distribution of habitats on the island to assess the potential distribution of the species beyond its currently known geographic range.

Phylogenomics recovers multiple origins of portable case making in caddisflies (Insecta: Trichoptera), nature's underwater architects

Caddisflies (Trichoptera) are among the most diverse groups of freshwater animals with more than 16 000 described species. They play a fundamental role in freshwater ecology and environmental engineering in streams, rivers and lakes. Because of this, they are frequently used as indicator organisms in biomonitoring programmes. Despite their importance, key questions concerning the evolutionary history of caddisflies, such as the timing and origin of larval case making, remain unanswered owing to the lack of a well-resolved phylogeny. Here, we estimated a phylogenetic tree using a combination of transcriptomes and targeted enrichment data for 207 species, representing 48 of 52 extant families and 174 genera. We calibrated and dated the tree with 33 carefully selected fossils. The first caddisflies originated approximately 295 million years ago in the Permian, and major suborders began to diversify in the Triassic. Furthermore, we show that portable case making evolved in three separate lineages, and shifts in diversification occurred in concert with key evolutionary innovations beyond case making.

Metazoa-level USCOs as markers in species delimitation and classification

Metazoa-level universal single-copy orthologs (mzl-USCOs) are universally applicable markers for DNA taxonomy in animals that can replace or supplement single-gene barcodes. Previously, mzl-USCOs from target enrichment data were shown to reliably distinguish species. Here, we tested whether USCOs are an evenly distributed, representative sample of a given metazoan genome and therefore able to cope with past hybridization events and incomplete lineage sorting. This is relevant for coalescent-based species delimitation approaches, which critically depend on the assumption that the investigated loci do not exhibit autocorrelation due to physical linkage. Based on 239 chromosome-level assembled genomes, we confirmed that mzl-USCOs are genetically unlinked for practical purposes and a representative sample of a genome in terms of reciprocal distances between USCOs on a chromosome and of distribution across chromosomes. We tested the suitability of mzl-USCOs extracted from genomes for species delimitation and phylogeny in four case studies: Anopheles mosquitos, Drosophila fruit flies, Heliconius butterflies and Darwin's finches. In almost all instances, USCOs allowed delineating species and yielded phylogenies that corresponded to those generated from whole genome data. Our phylogenetic analyses demonstrate that USCOs may complement single-gene DNA barcodes and provide more accurate taxonomic inferences. Combining USCOs from sources that used different versions of ortholog reference libraries to infer marker orthology may be challenging and, at times, impact taxonomic conclusions. However, we expect this problem to become less severe as the rapidly growing number of reference genomes provides a better representation of the number and diversity of organismal lineages.

Large-scale genomic data reveal the phylogeny and evolution of owlet moths (Noctuoidea)

The owlet moths (Noctuoidea; ~43-45K described species) are one of the most ecologically diverse and speciose superfamilies of animals. Moreover, they comprise some of the world's most notorious pests of agriculture and forestry. Despite their contributions to terrestrial biodiversity and impacts on ecosystems and economies, the evolutionary history of Noctuoidea remains unclear because the superfamily lacks a statistically robust phylogenetic and temporal framework. We reconstructed the phylogeny of Noctuoidea using data from 1234 genes (946.4 kb nucleotides) obtained from the genome and transcriptome sequences of 76 species. The relationships among the six families of Noctuoidea were well resolved and consistently recovered based on both concatenation and gene coalescence approaches, supporting the following relationships: Oenosandridae + (Notodontidae + (Erebidae + (Nolidae + (Euteliidae + Noctuidae)))). A Yule tree prior with three unlinked molecular clocks was identified as the preferred BEAST analysis using marginal-likelihood estimations. The crown age of Noctuoidea was estimated at 74.5 Ma, with most families originating before the end of the Paleogene (23 Ma). Our study provides the first statistically robust phylogenetic and temporal framework for Noctuoidea, including all families of owlet moths, based on large-scale genomic data.

Phylogenomics reveals the history of host use in mosquitoes

Mosquitoes have profoundly affected human history and continue to threaten human health through the transmission of a diverse array of pathogens. The phylogeny of mosquitoes has remained poorly characterized due to difficulty in taxonomic sampling and limited availability of genomic data beyond the most important vector species. Here, we used phylogenomic analysis of 709 single copy ortholog groups from 256 mosquito species to produce a strongly supported phylogeny that resolves the position of the major disease vector species and the major mosquito lineages. Our analyses support an origin of mosquitoes in the early Triassic (217 MYA [highest posterior density region: 188-250 MYA]), considerably older than previous estimates. Moreover, we utilize an extensive database of host associations for mosquitoes to show that mosquitoes have shifted to feeding upon the blood of mammals numerous times, and that mosquito diversification and host-use patterns within major lineages appear to coincide in earth history both with major continental drift events and with the diversification of vertebrate classes.

Potential Contribution of Ancient Introgression to the Evolution of a Derived Reproductive Strategy in Ricefishes

Transitions from no parental care to extensive care are costly and involve major changes in life history, behavior, and morphology. Nevertheless, in Sulawesi ricefishes, pelvic brooding evolved from transfer brooding in two distantly related lineages within the genera Adrianichthys and Oryzias, respectively. Females of pelvic brooding species carry their eggs attached to their belly until the fry hatches. Despite their phylogenetic distance, both pelvic brooding lineages share a set of external morphological traits. A recent study found no direct gene flow between pelvic brooding lineages, suggesting independent evolution of the derived reproductive strategy. Convergent evolution can, however, also rely on repeated sorting of preexisting variation of an admixed ancestral population, especially when subjected to similar external selection pressures. We thus used a multispecies coalescent model and D-statistics to identify gene-tree-species-tree incongruencies, to evaluate the evolution of pelvic brooding with respect to interspecific gene flow not only between pelvic brooding lineages but also between pelvic brooding lineages and other Sulawesi ricefish lineages. We found a general network-like evolution in Sulawesi ricefishes, and as previously reported, we detected no gene flow between the pelvic brooding lineages. Instead, we found hybridization between the ancestor of pelvic brooding Oryzias and the common ancestor of the Oryzias species from the Lake Poso area. We further detected signs of introgression within the confidence interval of a quantitative trait locus associated with pelvic brooding in O. eversi. Our results hint toward a contribution of ancient standing genetic variation to the evolution of pelvic brooding in Oryzias.

The genome of the glasshouse plant noble rhubarb (Rheum nobile) provides a window into alpine adaptation

Glasshouse plants are species that trap warmth via specialized morphology and physiology, mimicking a human glasshouse. In the Himalayan alpine region, the highly specialized glasshouse morphology has independently evolved in distinct lineages to adapt to intensive UV radiation and low temperature. Here we demonstrate that the glasshouse structure - specialized cauline leaves - is highly effective in absorbing UV light but transmitting visible and infrared light, creating an optimal microclimate for the development of reproductive organs. We reveal that this glasshouse syndrome has evolved at least three times independently in the rhubarb genus Rheum. We report the genome sequence of the flagship glasshouse plant Rheum nobile and identify key genetic network modules in association with the morphological transition to specialized glasshouse leaves, including active secondary cell wall biogenesis, upregulated cuticular cutin biosynthesis, and suppression of photosynthesis and terpenoid biosynthesis. The distinct cell wall organization and cuticle development might be important for the specialized optical property of glasshouse leaves. We also find that the expansion of LTRs has likely played an important role in noble rhubarb adaptation to high elevation environments. Our study will enable additional comparative analyses to identify the genetic basis underlying the convergent occurrence of glasshouse syndrome.

300 Million years of coral treaders (Insecta: Heteroptera: Hermatobatidae) back to the ocean in the phylogenetic context of Arthropoda

Among hundreds of insect families, Hermatobatidae (commonly known as coral treaders) is one of the most unique. They are small, wingless predaceous bugs in the suborder Heteroptera. Adults are almost black in colour, measuring about 5 mm in body length and 3 mm in width. Thirteen species are known from tropical coral reefs or rocky shores, but their origin and evolutionary adaptation to their unusual marine habitat were unexplored. We report here the genome and metagenome of Hermatobates lingyangjiaoensis, hitherto known only from its type locality in the South China Sea. We further reconstructed the evolutionary history and origin of these marine bugs in the broader context of Arthropoda. The dated phylogeny indicates that Hexapoda diverged from their marine sister groups approximately 498 Ma and that Hermatobatidae originated 192 Ma, indicating that they returned to an oceanic life some 300 Myr after their ancestors became terrestrial. Their origin is consistent with the recovery of tropical reef ecosystems after the end-Triassic mass extinction, which might have provided new and open niches for them to occupy and thrive. Our analyses also revealed that both the genome changes and the symbiotic bacteria might have contributed to adaptations necessary for life in the sea.

Advances in the understanding of Blattodea evolution: Insights from phylotranscriptomics and spermathecae

Cockroaches, an ancient and diverse group of insects on earth that originated in the Carboniferous, displays a wide array of morphology or biology diversity. The spermatheca is an organ of the insect reproductive system; the diversity of spermathecae might be the adaption to different mating and sperm storage strategies. Yet a consensus about the phylogenetic relationships among the main lineages of Blattodea and the evolution of spermatheca has not been reached until now. Here we added the transcriptome data of Anaplectidae for the first time and supplemented other family level groups (such as Blaberidae, Corydiidae) to address the pending issues. Our results showed that Blattoidea was recovered as sister to Corydioidea, which was strongly supported by molecular evidence. In Blattoidea, (Lamproblattidae + Anaplectidae) + (Cryptocercidae + Termitoidae) was strongly supported by our molecular data. In Blaberoidea, Pseudophyllodromiidae and Blaberidae were recovered to be monophyletic, while Blattellidae was found to be paraphyletic with respect to Malaccina. Ectobius sylvestris + Malaccina discoidalis formed the sister group to other Blaberoidea; Blattellidae (except Malaccina discoidalis) + Nyctiboridae was found as the sister of Blaberidae. Corydiidae was recovered to be non-monophyletic due to the embedding of Nocticola sp. Our ASR analysis of spermatheca suggested that primary spermathecae were present in the common ancestor, and it transformed at least six times during the evolutionary history of Blattodea. The evolution of spermatheca could be described as a unidirectional trend: the increased size to accommodate more sperm. Furthermore, major splits within the existing genera of cockroaches occurred in the Upper Paleogene to Neogene. Our study provides strong support for the relationship among three superfamilies and offers some new insights into the phylogeny of cockroaches. Meanwhile, this study also provides basic knowledge on the evolution of spermathecae and reproductive patterns.

The phylogeny and divergence times of leaf-mining flies (Diptera: Agromyzidae) from anchored phylogenomics

Leaf-mining flies (Diptera: Agromyzidae) are a diverse clade of phytophagous Diptera known largely for their economic impact as leaf- or stem-miners on vegetable and ornamental plants. Higher-level phylogenetic relationships of Agromyzidae have remained uncertain because of challenges in sampling of both taxa and characters for morphology and PCR-based Sanger-era molecular systematics. Here, we used hundreds of orthologous single-copy nuclear loci obtained from anchored hybrid enrichment (AHE) to reconstruct phylogenetic relationships among the major lineages of leaf-mining flies. The resulting phylogenetic trees are highly congruent and well-supported, except for a few deep nodes, when using different molecular data types and phylogenetic methods. Based on divergence time dating using a relaxed clock and model-based historical biogeography analysis, leaf-mining flies are shown to have diversified in multiple lineages since the early Paleocene, approximately 65 million years ago. Our study not only reveals a revised classification system of leaf-mining flies, but also provides a new phylogenetic framework to understand their macroevolution.

OrthoPhy: A Program to Construct Ortholog Data Sets Using Taxonomic Information

Species phylogenetic trees represent the evolutionary processes of organisms, and they are fundamental in evolutionary research. Therefore, new methods have been developed to obtain more reliable species phylogenetic trees. A highly reliable method is the construction of an ortholog data set based on sequence information of genes, which is then used to infer the species phylogenetic tree. However, although methods for constructing an ortholog data set for species phylogenetic analysis have been developed, they cannot remove some paralogs, which is necessary for reliable species phylogenetic inference. To address the limitations of current methods, we developed OrthoPhy, a program that excludes paralogs and constructs highly accurate ortholog data sets using taxonomic information dividing analyzed species into monophyletic groups. OrthoPhy can remove paralogs, detecting inconsistencies between taxonomic information and phylogenetic trees of candidate ortholog groups clustered by sequence similarity. Performance tests using evolutionary simulated sequences and real sequences of 40 bacteria revealed that the precision of ortholog inference by OrthoPhy is higher than that of existing programs. Additionally, the phylogenetic analysis of species was more accurate when performed using ortholog data sets constructed by OrthoPhy than that performed using data sets constructed by existing programs. Furthermore, we performed a benchmark test of the Quest for Orthologs using real sequence data and found that the concordance rate between the phylogenetic trees of orthologs inferred by OrthoPhy and those of species was higher than the rates obtained by other ortholog inference programs. Therefore, ortholog data sets constructed using OrthoPhy enabled a more accurate phylogenetic analysis of species than those constructed using the existing programs, and OrthoPhy can be used for the phylogenetic analysis of species even for distantly related species that have experienced many evolutionary events.

Phylogenomics of darkling beetles (Coleoptera: Tenebrionidae) from the Atacama Desert

Background

Tenebrionidae (Insecta: Coleoptera) are a conspicuous component of desert fauna worldwide. In these ecosystems, they are significantly responsible for nutrient cycling and show remarkable morphological and physiological adaptations. Nevertheless, Tenebrionidae colonizing individual deserts have repeatedly emerged from different lineages. The goal of our study was to gain insights into the phylogenetic relationships of the tenebrionid genera from the Atacama Desert and how these taxa are related to the globally distributed Tenebrionidae.

Methods

We used newly generated transcriptome data (47 tribes, 7 of 11 subfamilies) that allowed for a comprehensive phylogenomic analysis of the tenebrionid fauna of this hyperarid desert and fills a gap in our knowledge of the highly diversified Tenebrionidae. We examined two independent data sets known to be suitable for phylogenomic reconstructions. One is based on 35 neuropeptide precursors, the other on 1,742 orthologous genes shared among Coleoptera.

Results

The majority of Atacama genera are placed into three groups, two of which belong to typical South American lineages within the Pimeliinae. While the data support the monophyly of the Physogasterini, Nycteliini and Scotobiini, this does not hold for the Atacama genera of Edrotini, Epitragini, Evaniosomini, Praociini, Stenosini, Thinobatini, and Trilobocarini. A suggested very close relationship of Psammetichus with the Mediterranean Leptoderis also could not be confirmed. We also provide hints regarding the phylogenetic relationships of the Caenocrypticini, which occur both in South America and southern Africa. Apart from the focus on the Tenebrionidae from the Atacama Desert, we found a striking synapomorphy grouping Alleculinae, Blaptinae, Diaperinae, Stenochinae, and several taxa of Tenebrioninae, but not Tenebrio and Tribolium. This character, an insertion in the myosuppressin gene, defines a higher-level monophyletic group within the Tenebrionidae.

Conclusion

Transcriptome data allow a comprehensive phylogenomic analysis of the tenebrionid fauna of the Atacama Desert, which represents one of the seven major endemic tribal areas in the world for Tenebrionidae. Most Atacama genera could be placed in three lineages typical of South America; monophyly is not supported for several tribes based on molecular data, suggesting that a detailed systematic revision of several groups is necessary.

Phylogenetic resolution of the fly superfamily Ephydroidea-Molecular systematics of the enigmatic and diverse relatives of Drosophilidae

The schizophoran superfamily Ephydroidea (Diptera: Cyclorrhapha) includes eight families, ranging from the well-known vinegar flies (Drosophilidae) and shore flies (Ephydridae), to several small, relatively unusual groups, the phylogenetic placement of which has been particularly challenging for systematists. An extraordinary diversity in life histories, feeding habits and morphology are a hallmark of fly biology, and the Ephydroidea are no exception. Extreme specialization can lead to "orphaned" taxa with no clear evidence for their phylogenetic position. To resolve relationships among a diverse sample of Ephydroidea, including the highly modified flies in the families Braulidae and Mormotomyiidae, we conducted phylogenomic sampling. Using exon capture from Anchored Hybrid Enrichment and transcriptomics to obtain 320 orthologous nuclear genes sampled for 32 species of Ephydroidea and 11 outgroups, we evaluate a new phylogenetic hypothesis for representatives of the superfamily. These data strongly support monophyly of Ephydroidea with Ephydridae as an early branching radiation and the placement of Mormotomyiidae as a family-level lineage sister to all remaining families. We confirm placement of Cryptochetidae as sister taxon to a large clade containing both Drosophilidae and Braulidae-the latter a family of honeybee ectoparasites. Our results reaffirm that sampling of both taxa and characters is critical in hyperdiverse clades and that these factors have a major influence on phylogenomic reconstruction of the history of the schizophoran fly radiation.

A Practical Guide to Design and Assess a Phylogenomic Study

Over the last decade, molecular systematics has undergone a change of paradigm as high-throughput sequencing now makes it possible to reconstruct evolutionary relationships using genome-scale datasets. The advent of "big data" molecular phylogenetics provided a battery of new tools for biologists but simultaneously brought new methodological challenges. The increase in analytical complexity comes at the price of highly specific training in computational biology and molecular phylogenetics, resulting very often in a polarized accumulation of knowledge (technical on one side and biological on the other). Interpreting the robustness of genome-scale phylogenetic studies is not straightforward, particularly as new methodological developments have consistently shown that the general belief of "more genes, more robustness" often does not apply, and because there is a range of systematic errors that plague phylogenomic investigations. This is particularly problematic because phylogenomic studies are highly heterogeneous in their methodology, and best practices are often not clearly defined. The main aim of this article is to present what I consider as the ten most important points to take into consideration when planning a well-thought-out phylogenomic study and while evaluating the quality of published papers. The goal is to provide a practical step-by-step guide that can be easily followed by nonexperts and phylogenomic novices in order to assess the technical robustness of phylogenomic studies or improve the experimental design of a project.

Reconstruction of the rRNA Sequences of LUCA, with Bioinformatic Implication of the Local Similarities Shared by Them

The theory of the RNA world, especially with the catalytic capability of RNA, provides a reasonable framework explaining the evolution of molecular genetics system before the scenario of the central dogma. However, it remains a challenge to deduce the origin mechanism of rRNAs. Here we reconstructed the phylogenetic relationships of archaea and bacteria with bootstrap values of most nodes, especially the deep ones, higher than 90%. Based on the well-resolved tree, the full lengths of 16S, 5S, and 23S rRNA sequences of the last universal common ancestor (LUCA) were reconstructed for the first time. The potential similarities shared by the three ancestral rRNA sequences were further explored by searching for repeat short fragments in the level of purine–pyrimidine (RY) with certain lengths and arrangements. With the lengths ranging from 2 to 14, functional short fragments could be found in the three RNAs. As a representative, a set with a total of 75 short fragments of 11 nucleotides in length can recover all types of the known functional sites of ribosomes in a most concise manner. The 75 short fragments cluster around the functional center of the ribosome, among which 18 of them are highly conserved across five or six kingdoms and still contain all types of known functional sites except one. Alternatively, according to the strategy using the level of AUGC instead of RY, a similar pattern can be recovered. Such results indicate the local similarities shared by 16S, 5S, and 23S rRNAs and thus suggest a possible general mechanism in the formation of the LUCA rRNAs.

Phylogenomics resolves timing and patterns in the evolution of Australasian Cerambycinae (Coleoptera: Cerambycidae), and reveals new insights into the subfamily-level classification and historical biogeography of longhorn beetles

Cerambycinae is the second-largest subfamily of longhorn beetles in the Southern Hemisphere. The phylogeny of Cerambycinae is poorly known, resulting in a highly artificial tribal-level classification and a largely speculative evolutionary history. We reconstructed the phylogenetic relationships of Cerambycinae at the generic level using anchored hybrid enrichment data from hundreds of nuclear genes, with a primary focus on the extraordinarily diverse faunas of Australia and New Zealand. We also estimated divergence times by incorporating fossil calibrations in our analyses. We identified two main clades within Cerambycinae, which can also be separated morphologically by a distinct type of antennal foramen. We recovered a Late Jurassic origin of crown Cerambycinae. Dorcasominae, which was newly found to have representatives in Australia, was notably derived from within Cerambycinae. We recovered two independent origins of Australian Cerambycinae: one clade originated in the Early Cretaceous and is likely endemic to the Southern Hemisphere, while the other clade appears to have immigrated to Australia, perhaps from the Northern Hemisphere. Within the Australian lineages were multiple independent origins of New Zealand taxa, all of which are relative host-plant generalists. Tribal relationships and assignments are discussed, and based on our results, the following major nomenclatural acts were made: Dorcasominae Lacordaire, 1868 is downgraded to a tribe Dorcasomini of Cerambycinae Latreille, 1804; Neostenini Lacordaire, 1868syn. nov. is treated as a junior synonym of Uracanthini Blanchard, 1851.

Diversification of the phytophagous lineages of true bugs (Insecta: Hemiptera: Heteroptera) shortly after that of the flowering plants

More than 95% of phytophagous true bug (Hemiptera: Heteroptera) species belong to four superfamilies: Miroidea (Cimicomorpha), Pentatomoidea, Coreoidea, and Lygaeoidea (all Pentatomomorpha). These iconic groups of highly diverse, overwhelmingly phytophagous insects include several economically prominent agricultural and silvicultural pest species, though their evolutionary history has not yet been well resolved. In particular, superfamily- and family-level phylogenetic relationships of these four lineages have remained controversial, and the divergence times of some crucial nodes for phytophagous true bugs have hitherto been little known, which hampers a better understanding of the evolutionary processes and patterns of phytophagous insects. In the present study, we used 150 species and concatenated nuclear and mitochondrial protein-coding genes and rRNA genes to infer the phylogenetic relationships within the Terheteroptera (Cimicomorpha + Pentatomomorpha) and estimated their divergence times. Our results support the monophyly of Cimicomorpha, Pentatomomorpha, Miroidea, Pentatomoidea, Pyrrhocoroidea, Coreoidea, and Lygaeoidea. The phylogenetic relationships across phytophagous lineages are largely congruent at deep nodes across the analyses based on different datasets and tree-reconstructing methods with just a few exceptions. Estimated divergence times and ancestral state reconstructions for feeding habit indicate that phytophagous true bugs explosively radiated in the Early Cretaceous-shortly after the angiosperm radiation-with the subsequent diversification of the most speciose clades (Mirinae, Pentatomidae, Coreinae, and Rhyparochromidae) in the Late Cretaceous.

Phylotranscriptomic insights into a Mesoproterozoic-Neoproterozoic origin and early radiation of green seaweeds (Ulvophyceae)

The Ulvophyceae, a major group of green algae, is of particular evolutionary interest because of its remarkable morphological and ecological diversity. Its phylogenetic relationships and diversification timeline, however, are still not fully resolved. In this study, using an extensive nuclear gene dataset, we apply coalescent- and concatenation-based approaches to reconstruct the phylogeny of the Ulvophyceae and to explore the sources of conflict in previous phylogenomic studies. The Ulvophyceae is recovered as a paraphyletic group, with the Bryopsidales being a sister group to the Chlorophyceae, and the remaining taxa forming a clade (Ulvophyceae sensu stricto). Molecular clock analyses with different calibration strategies emphasize the large impact of fossil calibrations, and indicate a Meso-Neoproterozoic origin of the Ulvophyceae (sensu stricto), earlier than previous estimates. The results imply that ulvophyceans may have had a profound influence on oceanic redox structures and global biogeochemical cycles at the Mesoproterozoic-Neoproterozoic transition.

“Ulvophyceae is a remarkably morphologically and ecologically diverse clade of green algae. Here, the authors reconstruct the Ulvophyceae phylogeny, showing that these algae originated earlier than expected and may have influenced biogeochemical cycles at the Mesoproterozoic-Neoproterozoic transition.”

Transcriptomics provides a robust framework for the relationships of the major clades of cladobranch sea slugs (Mollusca, Gastropoda, Heterobranchia), but fails to resolve the position of the enigmatic genus Embletonia

Background

The soft-bodied cladobranch sea slugs represent roughly half of the biodiversity of marine nudibranch molluscs on the planet. Despite their global distribution from shallow waters to the deep sea, from tropical into polar seas, and their important role in marine ecosystems and for humans (as targets for drug discovery), the evolutionary history of cladobranch sea slugs is not yet fully understood.

Results

To enlarge the current knowledge on the phylogenetic relationships, we generated new transcriptome data for 19 species of cladobranch sea slugs and two additional outgroup taxa (Berthella plumula and Polycera quadrilineata). We complemented our taxon sampling with previously published transcriptome data, resulting in a final data set covering 56 species from all but one accepted cladobranch superfamilies. We assembled all transcriptomes using six different assemblers, selecting those assemblies that provided the largest amount of potentially phylogenetically informative sites. Quality-driven compilation of data sets resulted in four different supermatrices: two with full coverage of genes per species (446 and 335 single-copy protein-coding genes, respectively) and two with a less stringent coverage (667 genes with 98.9% partition coverage and 1767 genes with 86% partition coverage, respectively). We used these supermatrices to infer statistically robust maximum-likelihood trees. All analyses, irrespective of the data set, indicate maximal statistical support for all major splits and phylogenetic relationships at the family level. Besides the questionable position of Noumeaella rubrofasciata, rendering the Facelinidae as polyphyletic, the only notable discordance between the inferred trees is the position of Embletonia pulchra. Extensive testing using Four-cluster Likelihood Mapping, Approximately Unbiased tests, and Quartet Scores revealed that its position is not due to any informative phylogenetic signal, but caused by confounding signal.

Conclusions

Our data matrices and the inferred trees can serve as a solid foundation for future work on the taxonomy and evolutionary history of Cladobranchia. The placement of E. pulchra, however, proves challenging, even with large data sets and various optimization strategies. Moreover, quartet mapping results show that confounding signal present in the data is sufficient to explain the inferred position of E. pulchra, again leaving its phylogenetic position as an enigma.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01944-0.

Phylogenomic and mitogenomic data can accelerate inventorying of tropical beetles during the current biodiversity crisis

Conservation efforts must be evidence-based, so rapid and economically feasible methods should be used to quantify diversity and distribution patterns. We have attempted to overcome current impediments to the gathering of biodiversity data by using integrative phylogenomic and three mtDNA fragment analyses. As a model, we sequenced the Metriorrhynchini beetle fauna, sampled from ~700 localities in three continents. The species-rich dataset included ~6,500 terminals, ~1,850 putative species delimited at 5% uncorrected pairwise threshold, possibly ~1,000 of them unknown to science. Neither type of data could alone answer our questions on biodiversity and phylogeny. The phylogenomic backbone enabled the integrative delimitation of robustly defined natural genus-group units that will inform future research. Using constrained mtDNA analysis, we identified the spatial structure of species diversity, very high species-level endemism, and a biodiversity hotspot in New Guinea. We suggest that focused field research and subsequent laboratory and bioinformatic workflow steps would substantially accelerate the inventorying of any hyperdiverse tropical group with several thousand species. The outcome would be a scaffold for the incorporation of further data from environmental sequencing and ecological studies. The database of sequences could set a benchmark for the spatiotemporal evaluation of biodiversity, would support evidence-based conservation planning, and would provide a robust framework for systematic, biogeographic, and evolutionary studies.

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.

Lactitol Supplementation Modulates Intestinal Microbiome in Liver Cirrhotic Patients

Background: Cirrhosis is a common chronic liver disease characterized by irreversible diffuse liver damage. Intestinal microbiome dysbiosis and metabolite dysfunction contribute to the development of cirrhosis. Lactitol (4-β-D-galactopyranosyl-D-glucitol) was previously reported to promote the growth of intestinal Bifidobacteria. However, the effect of lactitol on the intestinal microbiome and fecal short-chain fatty acids (SCFAs) and bile acids (BAs) and the interactions among these factors in cirrhotic patients pre- and post-lactitol treatment remain poorly understood.

Methods: Here, using shotgun metagenomics and targeted metabolomics methods.

Results: we found that health-promoting lactic acid bacteria, including Bifidobacterium longum, B.pseudocatenulatum, and Lactobacillus salivarius, were increased after lactitol intervention, and significant decrease of pathogen Klebsiella pneumonia and associated antibiotic resistant genes /virulence factors. Functionally, pathways including Pseudomonas aeruginosa biofilm formation, endotoxin biosynthesis, and horizontal transfer of pathogenic genes were decreased in cirrhotic patients after 4-week lactitol intervention compared with before treatment.

Conclusion: We identified lactitol-associated metagenomic changes, and provide insight into the understanding of the roles of lactitol in modulating gut microbiome in cirrhotic patients.

Aspergillus sp. A31 and Curvularia geniculata P1 mitigate mercury toxicity to Oryza sativa L

Aspergillus sp. A31 and Curvularia geniculata P1 are endophytes that colonize the roots of Aeschynomene fluminensis Vell. and Polygonum acuminatum Kunth. in humid environments contaminated with mercury. The two strains mitigated mercury toxicity and promoted Oryza sativa L growth. C. geniculata P1 stood out for increasing the host biomass by fourfold and reducing the negative effects of the metal on photosynthesis. Assembling and annotation of Aspergillus sp. A31 and C. geniculata P1 genomes resulted in 28.60 Mb (CG% 53.1; 10,312 coding DNA sequences) and 32.92 Mb (CG% 50.72; 8,692 coding DNA sequences), respectively. Twelve and 27 genomes of Curvularia/Bipolaris and Aspergillus were selected for phylogenomic analyzes, respectively. Phylogenetic analysis inferred the separation of species from the genus Curvularia and Bipolaris into different clades, and the separation of species from the genus Aspergillus into three clades; the species were distinguished by occupied niche. The genomes had essential gene clusters for the adaptation of microorganisms to high metal concentrations, such as proteins of the phytoquelatin-metal complex (GO: 0090423), metal ion binders (GO: 0046872), ABC transporters (GO: 0042626), ATPase transporters (GO: 0016887), and genes related to response to reactive oxygen species (GO: 0000302) and oxidative stress (GO: 0006979). The results reported here help to understand the unique regulatory mechanisms of mercury tolerance and plant development.

Monophyletic blowflies revealed by phylogenomics

BACKGROUND

Blowflies are ubiquitous insects, often shiny and metallic, and the larvae of many species provide important ecosystem services (e.g., recycling carrion) and are used in forensics and debridement therapy. Yet, the taxon has repeatedly been recovered to be para- or polyphyletic, and the lack of a well-corroborated phylogeny has prevented a robust classification.

RESULTS

We here resolve the relationships between the different blowfly subclades by including all recognized subfamilies in a phylogenomic analysis using 2221 single-copy nuclear protein-coding genes of Diptera. Maximum likelihood (ML), maximum parsimony (MP), and coalescent-based phylogeny reconstructions all support the same relationships for the full data set. Based on this backbone phylogeny, blowflies are redefined as the most inclusive monophylum within the superfamily Oestroidea not containing Mesembrinellidae, Mystacinobiidae, Oestridae, Polleniidae, Sarcophagidae, Tachinidae, and Ulurumyiidae. The constituent subfamilies are re-classified as Ameniinae (including the Helicoboscinae, syn. nov.), Bengaliinae, Calliphorinae (including Aphyssurinae, syn. nov., Melanomyinae, syn. nov., and Toxotarsinae, syn. nov.), Chrysomyinae, Luciliinae, Phumosiinae, Rhiniinae stat. rev., and Rhinophorinae stat. rev. Metallic coloration in the adult is shown to be widespread but does not emerge as the most likely ground plan feature.

CONCLUSIONS

Our study provides the first phylogeny of oestroid calyptrates including all blowfly subfamilies. This allows settling a long-lasting controversy in Diptera by redefining blowflies as a well-supported monophylum, and blowfly classification is adjusted accordingly. The archetypical blowfly trait of carrion-feeding maggots most likely evolved twice, and the metallic color may not belong to the blowfly ground plan.

A phylotranscriptomic framework for flesh fly evolution (Diptera, Calyptratae, Sarcophagidae)

The Sarcophagidae (flesh flies) comprise a large and widely distributed radiation within the Calyptratae (Diptera). Larval feeding habits are ecologically diverse and include sarcosaprophagy, coprophagy, herbivory, invertebrate and vertebrate predation, and kleptoparasitism. To elucidate the geographic origin and evolution of flesh fly life-history, we inferred a backbone phylogeny based on transcriptomic data from 26 sarcophagid species covering all three subfamilies plus 15 outgroups. The phylogeny was inferred using maximum parsimony and maximum likelihood methods based on a series of supermatrices, one set with overall information content improved by MARE (2290 loci), one set with 100% gene coverage for all included species (587 loci), and the last set including mitochondrial and nuclear genes (589 loci) and additional taxa. In order to obtain a more detailed hypothesis, we utilized the supertree approach to combine results from the present study with previously published hypotheses. This resulted supertree covers 84 of the one hundred currently recognized sarcophagid genera and formed the basis for the ancestral state reconstructions. The monophyletic Sarcophagidae is well-supported as sister to {Mystacinobiidae + Oestridae}, and relationships at the subfamily level are inferred as {Sarcophaginae, (Paramacronychiinae + Miltogramminae)}. The Sarcophagidae and each subfamily originated in the Americas, with Sarcophaginae diversifying mainly in the Neotropics, whereas the major radiation of both Miltogramminae and Paramacronychiinae occurred in the Palaearctic. Sarcosaprophagy is reconstructed as the ancestral larval feeding habit of the family Sarcophagidae and each subfamily. The ancestral sarcophagid larva probably utilized dead invertebrates as food, and the food spectrum expanded together with the diversification of breeding strategies. Particularly, kleptoparasitism in Miltogramminae is derived from sarcosaprophagy and may be seen as having derived from the breeding biology of 'lower' miltogrammines, the larvae of which feed on buried vertebrate carrion.

Development and evaluation of a custom bait design based on 469 single-copy protein-coding genes for exon capture of isopods (Philosciidae: Haloniscus)

Transcriptome-based exon capture approaches, along with next-generation sequencing, are allowing for the rapid and cost-effective production of extensive and informative phylogenomic datasets from non-model organisms for phylogenetics and population genetics research. These approaches generally employ a reference genome to infer the intron-exon structure of targeted loci and preferentially select longer exons. However, in the absence of an existing and well-annotated genome, we applied this exon capture method directly, without initially identifying intron-exon boundaries for bait design, to a group of highly diverse Haloniscus (Philosciidae), paraplatyarthrid and armadillid isopods, and examined the performance of our methods and bait design for phylogenetic inference. Here, we identified an isopod-specific set of single-copy protein-coding loci, and a custom bait design to capture targeted regions from 469 genes, and analysed the resulting sequence data with a mapping approach and newly-created post-processing scripts. We effectively recovered a large and informative dataset comprising both short (<100 bp) and longer (>300 bp) exons, with high uniformity in sequencing depth. We were also able to successfully capture exon data from up to 16-year-old museum specimens along with more distantly related outgroup taxa, and efficiently pool multiple samples prior to capture. Our well-resolved phylogenies highlight the overall utility of this methodological approach and custom bait design, which offer enormous potential for application to future isopod, as well as broader crustacean, molecular studies.

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.

Fifty million years of beetle evolution along the Antarctic Polar Front

The Antarctic environment is famously inhospitable to most terrestrial biodiversity, traditionally viewed as a driver of species extinction. Combining population- and species-level molecular data, we show that beetles on islands along the Antarctic Polar Front diversified in response to major climatic events over the last 50 Ma in surprising synchrony with the region’s marine organisms. Unique algae- and moss-feeding habits enabled beetles to capitalize on cooling conditions, which resulted in a decline in flowering plants—the typical hosts for beetles elsewhere. Antarctica’s cooling paleoclimate thus fostered the diversification of both terrestrial and marine life. Climatically driven evolutionary processes since the Miocene may underpin much of the region’s diversity, are still ongoing, and should be further investigated among Antarctic biota.

Global cooling and glacial–interglacial cycles since Antarctica’s isolation have been responsible for the diversification of the region’s marine fauna. By contrast, these same Earth system processes are thought to have played little role terrestrially, other than driving widespread extinctions. Here, we show that on islands along the Antarctic Polar Front, paleoclimatic processes have been key to diversification of one of the world’s most geographically isolated and unique groups of herbivorous beetles—Ectemnorhinini weevils. Combining phylogenomic, phylogenetic, and phylogeographic approaches, we demonstrate that these weevils colonized the sub-Antarctic islands from Africa at least 50 Ma ago and repeatedly dispersed among them. As the climate cooled from the mid-Miocene, diversification of the beetles accelerated, resulting in two species-rich clades. One of these clades specialized to feed on cryptogams, typical of the polar habitats that came to prevail under Miocene conditions yet remarkable as a food source for any beetle. This clade’s most unusual representative is a marine weevil currently undergoing further speciation. The other clade retained the more common weevil habit of feeding on angiosperms, which likely survived glaciation in isolated refugia. Diversification of Ectemnorhinini weevils occurred in synchrony with many other Antarctic radiations, including penguins and notothenioid fishes, and coincided with major environmental changes. Our results thus indicate that geo-climatically driven diversification has progressed similarly for Antarctic marine and terrestrial organisms since the Miocene, potentially constituting a general biodiversity paradigm that should be sought broadly for the region’s taxa.

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.

Is Phylotranscriptomics as Reliable as Phylogenomics?

Phylogenomics, the study of phylogenetic relationships among taxa based on their genome sequences, has emerged as the preferred phylogenetic method because of the wealth of phylogenetic information contained in genome sequences. Genome sequencing, however, can be prohibitively expensive, especially for taxa with huge genomes and when many taxa need sequencing. Consequently, the less costly phylotranscriptomics has seen an increased use in recent years. Phylotranscriptomics reconstructs phylogenies using DNA sequences derived from transcriptomes, which are often orders of magnitude smaller than genomes. However, in the absence of corresponding genome sequences, comparative analyses of transcriptomes can be challenging and it is unclear whether phylotranscriptomics is as reliable as phylogenomics. Here, we respectively compare the phylogenomic and phylotranscriptomic trees of 22 mammals and 15 plants that have both sequenced nuclear genomes and publicly available RNA sequencing data from multiple tissues. We found that phylotranscriptomic analysis can be sensitive to orthologous gene identification. When a rigorous method for identifying orthologs is employed, phylogenomic and phylotranscriptomic trees are virtually identical to each other, regardless of the tissue of origin of the transcriptomes and whether the same tissue is used across species. These findings validate phylotranscriptomics, brighten its prospect, and illustrate the criticality of reliable ortholog detection in such practices.

Large Phylogenomic Data sets Reveal Deep Relationships and Trait Evolution in Chlorophyte Green Algae

The chlorophyte green algae (Chlorophyta) are species-rich ancient groups ubiquitous in various habitats with high cytological diversity, ranging from microscopic to macroscopic organisms. However, the deep phylogeny within core Chlorophyta remains unresolved, in part due to the relatively sparse taxon and gene sampling in previous studies. Here we contribute new transcriptomic data and reconstruct phylogenetic relationships of core Chlorophyta based on four large data sets up to 2,698 genes of 70 species, representing 80% of extant orders. The impacts of outgroup choice, missing data, bootstrap-support cutoffs, and model misspecification in phylogenetic inference of core Chlorophyta are examined. The species tree topologies of core Chlorophyta from different analyses are highly congruent, with strong supports at many relationships (e.g., the Bryopsidales and the Scotinosphaerales-Dasycladales clade). The monophyly of Chlorophyceae and of Trebouxiophyceae as well as the uncertain placement of Chlorodendrophyceae and Pedinophyceae corroborate results from previous studies. The reconstruction of ancestral scenarios illustrates the evolution of the freshwater-sea and microscopic–macroscopic transition in the Ulvophyceae, and the transformation of unicellular→colonial→multicellular in the chlorophyte green algae. In addition, we provided new evidence that serine is encoded by both canonical codons and noncanonical TAG code in Scotinosphaerales, and stop-to-sense codon reassignment in the Ulvophyceae has originated independently at least three times. Our robust phylogenetic framework of core Chlorophyta unveils the evolutionary history of phycoplast, cyto-morphology, and noncanonical genetic codes in chlorophyte green algae.

Large-Scale Phylogenomic Analyses Reveal the Monophyly of Bryophytes and Neoproterozoic Origin of Land Plants

The relationships among the four major embryophyte lineages (mosses, liverworts, hornworts, vascular plants) and the timing of the origin of land plants are enigmatic problems in plant evolution. Here, we resolve the monophyly of bryophytes by improving taxon sampling of hornworts and eliminating the effect of synonymous substitutions. We then estimate the divergence time of crown embryophytes based on three fossil calibration strategies, and reveal that maximum calibration constraints have a major effect on estimating the time of origin of land plants. Moreover, comparison of priors and posteriors provides a guide for evaluating the optimal calibration strategy. By considering the reliability of fossil calibrations and the influences of molecular data, we estimate that land plants originated in the Precambrian (980–682 Ma), much older than widely recognized. Our study highlights the important contribution of molecular data when faced with contentious fossil evidence, and that fossil calibrations used in estimating the timescale of plant evolution require critical scrutiny.

New alignment-based sequence extraction software (ALiBaSeq) and its utility for deep level phylogenetics

Despite many bioinformatic solutions for analyzing sequencing data, few options exist for targeted sequence retrieval from whole genomic sequencing (WGS) data with the ultimate goal of generating a phylogeny. Available tools especially struggle at deep phylogenetic levels and necessitate amino-acid space searches, which may increase rates of false positive results. Many tools are also difficult to install and may lack adequate user resources. Here, we describe a program that uses freely available similarity search tools to find homologs in assembled WGS data with unparalleled freedom to modify parameters. We evaluate its performance compared to other commonly used bioinformatics tools on two divergent insect species (>200 My) for which annotated genomes exist, and on one large set each of highly conserved and more variable loci. Our software is capable of retrieving orthologs from well-curated or unannotated, low or high depth shotgun, and target capture assemblies as well or better than other software as assessed by recovering the most genes with maximal coverage and with a low rate of false positives throughout all datasets. When assessing this combination of criteria, ALiBaSeq is frequently the best evaluated tool for gathering the most comprehensive and accurate phylogenetic alignments on all types of data tested. The software (implemented in Python), tutorials, and manual are freely available at https://github.com/AlexKnyshov/alibaseq.

DeepNOG: fast and accurate protein orthologous group assignment

MOTIVATION

Protein orthologous group databases are powerful tools for evolutionary analysis, functional annotation or metabolic pathway modeling across lineages. Sequences are typically assigned to orthologous groups with alignment-based methods, such as profile hidden Markov models, which have become a computational bottleneck.

RESULTS

We present DeepNOG, an extremely fast and accurate, alignment-free orthology assignment method based on deep convolutional networks. We compare DeepNOG against state-of-the-art alignment-based (HMMER, DIAMOND) and alignment-free methods (DeepFam) on two orthology databases (COG, eggNOG 5). DeepNOG can be scaled to large orthology databases like eggNOG, for which it outperforms DeepFam in terms of precision and recall by large margins. While alignment-based methods still provide the most accurate assignments among the investigated methods, computing time of DeepNOG is an order of magnitude lower on CPUs. Optional GPU usage further increases throughput massively. A command-line tool enables rapid adoption by users.

AVAILABILITYAND IMPLEMENTATION

Source code and packages are freely available at https://github.com/univieCUBE/deepnog. Install the platform-independent Python program with $pip install deepnog.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Ancestral sequence reconstruction - An underused approach to understand the evolution of gene function in plants?

Whilst substantial research effort has been placed on understanding the interactions of plant proteins with their molecular partners, relatively few studies in plants - by contrast to work in other organisms - address how these interactions evolve. It is thought that ancestral proteins were more promiscuous than modern proteins and that specificity often evolved following gene duplication and subsequent functional refining. However, ancestral protein resurrection studies have found that some modern proteins have evolved de novo from ancestors lacking those functions. Intriguingly, the new interactions evolved as a consequence of just a few mutations and, as such, acquisition of new functions appears to be neither difficult nor rare, however, only a few of them are incorporated into biological processes before they are lost to subsequent mutations. Here, we detail the approach of ancestral sequence reconstruction (ASR), providing a primer to reconstruct the sequence of an ancestral gene. We will present case studies from a range of different eukaryotes before discussing the few instances where ancestral reconstructions have been used in plants. As ASR is used to dig into the remote evolutionary past, we will also present some alternative genetic approaches to investigate molecular evolution on shorter timescales. We argue that the study of plant secondary metabolism is particularly well suited for ancestral reconstruction studies. Indeed, its ancient evolutionary roots and highly diverse landscape provide an ideal context in which to address the focal issue around the emergence of evolutionary novelties and how this affects the chemical diversification of plant metabolism.

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.

Beyond Drosophila: resolving the rapid radiation of schizophoran flies with phylotranscriptomics

BACKGROUND

The most species-rich radiation of animal life in the 66 million years following the Cretaceous extinction event is that of schizophoran flies: a third of fly diversity including Drosophila fruit fly model organisms, house flies, forensic blow flies, agricultural pest flies, and many other well and poorly known true flies. Rapid diversification has hindered previous attempts to elucidate the phylogenetic relationships among major schizophoran clades. A robust phylogenetic hypothesis for the major lineages containing these 55,000 described species would be critical to understand the processes that contributed to the diversity of these flies. We use protein encoding sequence data from transcriptomes, including 3145 genes from 70 species, representing all superfamilies, to improve the resolution of this previously intractable phylogenetic challenge.

RESULTS

Our results support a paraphyletic acalyptrate grade including a monophyletic Calyptratae and the monophyly of half of the acalyptrate superfamilies. The primary branching framework of Schizophora is well supported for the first time, revealing the primarily parasitic Pipunculidae and Sciomyzoidea stat. rev. as successive sister groups to the remaining Schizophora. Ephydroidea, Drosophila's superfamily, is the sister group of Calyptratae. Sphaeroceroidea has modest support as the sister to all non-sciomyzoid Schizophora. We define two novel lineages corroborated by morphological traits, the 'Modified Oviscapt Clade' containing Tephritoidea, Nerioidea, and other families, and the 'Cleft Pedicel Clade' containing Calyptratae, Ephydroidea, and other families. Support values remain low among a challenging subset of lineages, including Diopsidae. The placement of these families remained uncertain in both concatenated maximum likelihood and multispecies coalescent approaches. Rogue taxon removal was effective in increasing support values compared with strategies that maximise gene coverage or minimise missing data.

CONCLUSIONS

Dividing most acalyptrate fly groups into four major lineages is supported consistently across analyses. Understanding the fundamental branching patterns of schizophoran flies provides a foundation for future comparative research on the genetics, ecology, and biocontrol.

Comparative transcriptomics of ice-crawlers demonstrates cold specialization constrains niche evolution in a relict lineage

Key changes in ecological niche space are often critical to understanding how lineages diversify during adaptive radiations. However, the converse, or understanding why some lineages are depauperate and relictual, is more challenging, as many factors may constrain niche evolution. In the case of the insect order Grylloblattodea, highly conserved thermal breadth is assumed to be closely tied to their relictual status, but has not been formerly tested. Here, we investigate whether evolutionary constraints in the physiological tolerance of temperature can help explain relictualism in this lineage. Using a comparative transcriptomics approach, we investigate gene expression following acute heat and cold stress across members of Grylloblattodea and their sister group, Mantophasmatodea. We additionally examine patterns of protein evolution, to identify candidate genes of positive selection. We demonstrate that cold specialization in Grylloblattodea has been accompanied by the loss of the inducible heat shock response under both acute heat and cold stress. Additionally, there is widespread evidence of selection on protein-coding genes consistent with evolutionary constraints due to cold specialization. This includes positive selection on genes involved in trehalose transport, metabolic function, mitochondrial function, oxygen reduction, oxidative stress, and protein synthesis. These patterns of molecular adaptation suggest that Grylloblattodea have undergone evolutionary trade-offs to survive in cold habitats and should be considered highly vulnerable to climate change. Finally, our transcriptomic data provide a robust backbone phylogeny for generic relationships within Grylloblattodea and Mantophasmatodea. Major phylogenetic splits in each group relate to arid conditions driving biogeographical patterns, with support for a sister-group relationship between North American Grylloblatta and Altai-Sayan Grylloblattella, and a range disjunction in Namibia splitting major clades within Mantophasmatodea.

Horizontally acquired antibacterial genes associated with adaptive radiation of ladybird beetles

BACKGROUND

Horizontal gene transfer (HGT) has been documented in many herbivorous insects, conferring the ability to digest plant material and promoting their remarkable ecological diversification. Previous reports suggest HGT of antibacterial enzymes may have contributed to the insect immune response and limit bacterial growth. Carnivorous insects also display many evolutionary successful lineages, but in contrast to the plant feeders, the potential role of HGTs has been less well-studied.

RESULTS

Using genomic and transcriptomic data from 38 species of ladybird beetles, we identified a set of bacterial cell wall hydrolase (cwh) genes acquired by this group of beetles. Infection with Bacillus subtilis led to upregulated expression of these ladybird cwh genes, and their recombinantly produced proteins limited bacterial proliferation. Moreover, RNAi-mediated cwh knockdown led to downregulation of other antibacterial genes, indicating a role in antibacterial immune defense. cwh genes are rare in eukaryotes, but have been maintained in all tested Coccinellinae species, suggesting that this putative immune-related HGT event played a role in the evolution of this speciose subfamily of predominant predatory ladybirds.

CONCLUSION

Our work demonstrates that, in a manner analogous to HGT-facilitated plant feeding, enhanced immunity through HGT might have played a key role in the prey adaptation and niche expansion that promoted the diversification of carnivorous beetle lineages. We believe that this represents the first example of immune-related HGT in carnivorous insects with an association with a subsequent successful species radiation.

Four myriapod relatives - but who are sisters? No end to debates on relationships among the four major myriapod subgroups

Background

Phylogenetic relationships among the myriapod subgroups Chilopoda, Diplopoda, Symphyla and Pauropoda are still not robustly resolved. The first phylogenomic study covering all subgroups resolved phylogenetic relationships congruently to morphological evidence but is in conflict with most previously published phylogenetic trees based on diverse molecular data. Outgroup choice and long-branch attraction effects were stated as possible explanations for these incongruencies. In this study, we addressed these issues by extending the myriapod and outgroup taxon sampling using transcriptome data.

Results

We generated new transcriptome data of 42 panarthropod species, including all four myriapod subgroups and additional outgroup taxa. Our taxon sampling was complemented by published transcriptome and genome data resulting in a supermatrix covering 59 species. We compiled two data sets, the first with a full coverage of genes per species (292 single-copy protein-coding genes), the second with a less stringent coverage (988 genes). We inferred phylogenetic relationships among myriapods using different data types, tree inference, and quartet computation approaches. Our results unambiguously support monophyletic Mandibulata and Myriapoda. Our analyses clearly showed that there is strong signal for a single unrooted topology, but a sensitivity of the position of the internal root on the choice of outgroups. However, we observe strong evidence for a clade Pauropoda+Symphyla, as well as for a clade Chilopoda+Diplopoda.

Conclusions

Our best quartet topology is incongruent with current morphological phylogenies which were supported in another phylogenomic study. AU tests and quartet mapping reject the quartet topology congruent to trees inferred with morphological characters. Moreover, quartet mapping shows that confounding signal present in the data set is sufficient to explain the weak signal for the quartet topology derived from morphological characters. Although outgroup choice affects results, our study could narrow possible trees to derivatives of a single quartet topology. For highly disputed relationships, we propose to apply a series of tests (AU and quartet mapping), since results of such tests allow to narrow down possible relationships and to rule out confounding signal.

Phylogenomic analysis sheds light on the evolutionary pathways towards acoustic communication in Orthoptera

Acoustic communication is enabled by the evolution of specialised hearing and sound producing organs. In this study, we performed a large-scale macroevolutionary study to understand how both hearing and sound production evolved and affected diversification in the insect order Orthoptera, which includes many familiar singing insects, such as crickets, katydids, and grasshoppers. Using phylogenomic data, we firmly establish phylogenetic relationships among the major lineages and divergence time estimates within Orthoptera, as well as the lineage-specific and dynamic patterns of evolution for hearing and sound producing organs. In the suborder Ensifera, we infer that forewing-based stridulation and tibial tympanal ears co-evolved, but in the suborder Caelifera, abdominal tympanal ears first evolved in a non-sexual context, and later co-opted for sexual signalling when sound producing organs evolved. However, we find little evidence that the evolution of hearing and sound producing organs increased diversification rates in those lineages with known acoustic communication.

Measuring Genome Sizes Using Read-Depth, k-mers, and Flow Cytometry: Methodological Comparisons in Beetles (Coleoptera)

Measuring genome size across different species can yield important insights into evolution of the genome and allow for more informed decisions when designing next-generation genomic sequencing projects. New techniques for estimating genome size using shallow genomic sequence data have emerged which have the potential to augment our knowledge of genome sizes, yet these methods have only been used in a limited number of empirical studies. In this project, we compare estimation methods using next-generation sequencing (k-mer methods and average read depth of single-copy genes) to measurements from flow cytometry, a standard method for genome size measures, using ground beetles (Carabidae) and other members of the beetle suborder Adephaga as our test system. We also present a new protocol for using read-depth of single-copy genes to estimate genome size. Additionally, we report flow cytometry measurements for five previously unmeasured carabid species, as well as 21 new draft genomes and six new draft transcriptomes across eight species of adephagan beetles. No single sequence-based method performed well on all species, and all tended to underestimate the genome sizes, although only slightly in most samples. For one species, Bembidion sp. nr. transversale, most sequence-based methods yielded estimates half the size suggested by flow cytometry.

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.

Phylogenetic tree building in the genomic age

Knowing phylogenetic relationships among species is fundamental for many studies in biology. An accurate phylogenetic tree underpins our understanding of the major transitions in evolution, such as the emergence of new body plans or metabolism, and is key to inferring the origin of new genes, detecting molecular adaptation, understanding morphological character evolution and reconstructing demographic changes in recently diverged species. Although data are ever more plentiful and powerful analysis methods are available, there remain many challenges to reliable tree building. Here, we discuss the major steps of phylogenetic analysis, including identification of orthologous genes or proteins, multiple sequence alignment, and choice of substitution models and inference methodologies. Understanding the different sources of errors and the strategies to mitigate them is essential for assembling an accurate tree of life.

The Parallel Molecular Adaptations to the Antarctic Cold Environment in Two Psychrophilic Green Algae

Psychrophilic green algae from independent phylogenetic lines thrive in the polar extreme environments, but the hypothesis that their psychrophilic characteristics appeared through parallel routes of molecular evolution remains untested. The recent surge of transcriptome data enables large-scale evolutionary analyses to investigate the genetic basis for the adaptations to the Antarctic extreme environment, and the identification of the selective forces that drive molecular evolution is the foundation to understand the strategies of cold adaptation. Here, we conducted transcriptome sequencing of two Antarctic psychrophilic green algae (Chlamydomonas sp. ICE-L and Tetrabaena socialis) and performed positive selection and convergent substitution analyses to investigate their molecular convergence and adaptive strategies against extreme cold conditions. Our results revealed considerable shared positively selected genes and significant evidence of molecular convergence in two Antarctic psychrophilic algae. Significant evidence of positive selection and convergent substitution were detected in genes associated with photosynthetic machinery, multiple antioxidant systems, and several crucial translation elements in Antarctic psychrophilic algae. Our study reveals that the psychrophilic algae possess more stable photosynthetic apparatus and multiple protective mechanisms and provides new clues of parallel adaptive evolution in Antarctic psychrophilic green algae.

Phylogenomics reveals the evolutionary timing and pattern of butterflies and moths

Lepidoptera play key roles in many biological systems. Butterflies are hypothesized to have evolved contemporaneously with flowering plants, and moths are thought to have gained anti-bat defenses in response to echolocating predatory bats, but these hypotheses have largely gone untested. Using a transcriptomic, dated evolutionary tree of Lepidoptera, we demonstrate that the most recent common ancestor of Lepidoptera is considerably older than previously hypothesized. The oldest moths in crown Lepidoptera were present in the Carboniferous, some 300 million years ago, and began to diversify largely in synchrony with angiosperms. We show that multiple lineages of moths independently evolved hearing organs well before the origin of bats, rejecting the hypothesis that lepidopteran hearing organs arose in response to these predators.

Butterflies and moths (Lepidoptera) are one of the major superradiations of insects, comprising nearly 160,000 described extant species. As herbivores, pollinators, and prey, Lepidoptera play a fundamental role in almost every terrestrial ecosystem. Lepidoptera are also indicators of environmental change and serve as models for research on mimicry and genetics. They have been central to the development of coevolutionary hypotheses, such as butterflies with flowering plants and moths’ evolutionary arms race with echolocating bats. However, these hypotheses have not been rigorously tested, because a robust lepidopteran phylogeny and timing of evolutionary novelties are lacking. To address these issues, we inferred a comprehensive phylogeny of Lepidoptera, using the largest dataset assembled for the order (2,098 orthologous protein-coding genes from transcriptomes of 186 species, representing nearly all superfamilies), and dated it with carefully evaluated synapomorphy-based fossils. The oldest members of the Lepidoptera crown group appeared in the Late Carboniferous (∼300 Ma) and fed on nonvascular land plants. Lepidoptera evolved the tube-like proboscis in the Middle Triassic (∼241 Ma), which allowed them to acquire nectar from flowering plants. This morphological innovation, along with other traits, likely promoted the extraordinary diversification of superfamily-level lepidopteran crown groups. The ancestor of butterflies was likely nocturnal, and our results indicate that butterflies became day-flying in the Late Cretaceous (∼98 Ma). Moth hearing organs arose multiple times before the evolutionary arms race between moths and bats, perhaps initially detecting a wide range of sound frequencies before being co-opted to specifically detect bat sonar. Our study provides an essential framework for future comparative studies on butterfly and moth evolution.