Using ultraconserved elements to reconstruct the termite tree of life

Phylogenomic analyses of Blattodea combining traditional methods, incremental tree-building, and quality-aware support

Despite the many advances of the genomic era, there is a persistent problem in assessing the uncertainty of phylogenomic hypotheses. We see this in the recent history of phylogenetics for cockroaches and termites (Blattodea), where huge advances have been made, but there are still major inconsistencies between studies. To address this, we present a phylogenetic analysis of Blattodea that emphasizes identification and quantification of uncertainty. We analyze 1183 gene domains using three methods (multi-species coalescent inference, concatenation, and a supermatrix-supertree hybrid approach) and assess support for controversial relationships while considering data quality. The hybrid approach-here dubbed "tiered phylogenetic inference"-incorporates information about data quality into an incremental tree building framework. Leveraging this method, we are able to identify cases of low or misleading support that would not be possible otherwise, and explore them more thoroughly with follow-up tests. In particular, quality annotations pointed towards nodes with high bootstrap support that later turned out to have large ambiguities, sometimes resulting from low-quality data. We also clarify issues related to some recalcitrant nodes: Anaplectidae's placement lacks unbiased signal, Ectobiidae s.s. and Anaplectoideini need greater taxon sampling, the deepest relationships among most Blaberidae lack signal. As a result, several previous phylogenetic uncertainties are now closer to being resolved (e.g., African and Malagasy "Rhabdoblatta" spp. are the sister to all other Blaberidae, and Oxyhaloinae is sister to the remaining Blaberidae). Overall, we argue for more approaches to quantifying support that take data quality into account to uncover the nature of recalcitrant nodes.

Dating in the Dark: Elevated Substitution Rates in Cave Cockroaches (Blattodea: Nocticolidae) Have Negative Impacts on Molecular Date Estimates

Rates of nucleotide substitution vary substantially across the Tree of Life, with potentially confounding effects on phylogenetic and evolutionary analyses. A large acceleration in mitochondrial substitution rate occurs in the cockroach family Nocticolidae, which predominantly inhabit subterranean environments. To evaluate the impacts of this among-lineage rate heterogeneity on estimates of phylogenetic relationships and evolutionary timescales, we analyzed nuclear ultraconserved elements (UCEs) and mitochondrial genomes from nocticolids and other cockroaches. Substitution rates were substantially elevated in nocticolid lineages compared with other cockroaches, especially in mitochondrial protein-coding genes. This disparity in evolutionary rates is likely to have led to different evolutionary relationships being supported by phylogenetic analyses of mitochondrial genomes and UCE loci. Furthermore, Bayesian dating analyses using relaxed-clock models inferred much deeper divergence times compared with a flexible local clock. Our phylogenetic analysis of UCEs, which is the first genome-scale study to include all 13 major cockroach families, unites Corydiidae and Nocticolidae and places Anaplectidae as the sister lineage to the rest of Blattoidea. We uncover an extraordinary level of genetic divergence in Nocticolidae, including two highly distinct clades that separated ~115 million years ago despite both containing representatives of the genus Nocticola. The results of our study highlight the potential impacts of high among-lineage rate variation on estimates of phylogenetic relationships and evolutionary timescales.

Genomic data provide insights into the classification of extant termites

The higher classification of termites requires substantial revision as the Neoisoptera, the most diverse termite lineage, comprise many paraphyletic and polyphyletic higher taxa. Here, we produce an updated termite classification using genomic-scale analyses. We reconstruct phylogenies under diverse substitution models with ultraconserved elements analyzed as concatenated matrices or within the multi-species coalescence framework. Our classification is further supported by analyses controlling for rogue loci and taxa, and topological tests. We show that the Neoisoptera are composed of seven family-level monophyletic lineages, including the Heterotermitidae Froggatt, Psammotermitidae Holmgren, and Termitogetonidae Holmgren, raised from subfamilial rank. The species-rich Termitidae are composed of 18 subfamily-level monophyletic lineages, including the new subfamilies Crepititermitinae, Cylindrotermitinae, Forficulitermitinae, Neocapritermitinae, Protohamitermitinae, and Promirotermitinae; and the revived Amitermitinae Kemner, Microcerotermitinae Holmgren, and Mirocapritermitinae Kemner. Building an updated taxonomic classification on the foundation of unambiguously supported monophyletic lineages makes it highly resilient to potential destabilization caused by the future availability of novel phylogenetic markers and methods. The taxonomic stability is further guaranteed by the modularity of the new termite classification, designed to accommodate as-yet undescribed species with uncertain affinities to the herein delimited monophyletic lineages in the form of new families or subfamilies.

Ultraconserved elements from transcriptome and genome data provide insight into the phylogenomics of Sternorrhyncha (Insecta: Hemiptera)

Sternorrhyncha, one of the four major suborders of Hemiptera, is a phytophagous taxon inclusive of nearly 18 000 described species. The phylogenetic relationships within the taxon and the earliest-branching lineage of its infraorders remain incompletely understood. This study attempted to illuminate the phylogenetic relationships within Sternorrhyncha through the use of maximum likelihood, Bayesian inference and maximum parsimony analyses, employing ultraconserved element (UCE) data from 39 genomic and 62 transcriptomic datasets and thereby representing most families within the taxon. The probe set Hemiptera 2.7Kv1 was used to recover a total of 2731 UCE loci: from 547 to 1699 (with an average of 1084) across all genomic datasets and from 108 to 849 (with an average of 329) across all transcriptomic datasets. All three types of phylogenetic analyses employed in this study produced robust statistical support for Sternorrhyncha being a monophyletic group. The different methods of phylogenetic analysis produced inconsistent descriptions of topological structure at the infraorder level: while maximum likelihood and Bayesian inference analyses produced strong statistical evidence (100%) indicating the clade Psylloidea + Aleyrodoidea to be a sister of the clade Aphidoidea (Aphidomorpha) + Coccoidea (Coccomorpha), the maximum parsimony analysis failed to recover a similar result. Our results also provide detail on the phylogenetic relationships within each infraorder. This study presents the first use of UCE data to investigate the phylogeny of Sternorrhyncha. It also shows the viability of amalgamating genomic and transcriptomic data in studies of phylogenetic relationships, potentially highlighting a resource-efficient approach for future inquiries into diverse taxa through the integration of varied data sources.

Comparative analysis of adipokinetic hormones and their receptors in Blattodea reveals novel patterns of gene evolution

Adipokinetic hormone (AKH) is a neuropeptide produced in the insect corpora cardiaca that plays an essential role in mobilising carbohydrates and lipids from the fat body to the haemolymph. AKH acts by binding to a rhodopsin-like G protein-coupled receptor (GPCR), the adipokinetic hormone receptor (AKHR). In this study, we tackle AKH ligand and receptor gene evolution as well as the evolutionary origins of AKH gene paralogues from the order Blattodea (termites and cockroaches). Phylogenetic analyses of AKH precursor sequences point to an ancient AKH gene duplication event in the common ancestor of Blaberoidea, yielding a new group of putative decapeptides. In total, 16 different AKH peptides from 90 species were obtained. Two octapeptides and seven putatively novel decapeptides are predicted for the first time. AKH receptor sequences from 18 species, spanning solitary cockroaches and subsocial wood roaches as well as lower and higher termites, were subsequently acquired using classical molecular methods and in silico approaches employing transcriptomic data. Aligned AKHR open reading frames revealed 7 highly conserved transmembrane regions, a typical arrangement for GPCRs. Phylogenetic analyses based on AKHR sequences support accepted relationships among termite, subsocial (Cryptocercus spp.) and solitary cockroach lineages to a large extent, while putative post-translational modification sites do not greatly differ between solitary and subsocial roaches and social termites. Our study provides important information not only for AKH and AKHR functional research but also for further analyses interested in their development as potential candidates for biorational pest control agents against invasive termites and cockroaches.

Evidence of cospeciation between termites and their gut bacteria on a geological time scale

Termites host diverse communities of gut microbes, including many bacterial lineages only found in this habitat. The bacteria endemic to termite guts are transmitted via two routes: a vertical route from parent colonies to daughter colonies and a horizontal route between colonies sometimes belonging to different termite species. The relative importance of both transmission routes in shaping the gut microbiota of termites remains unknown. Using bacterial marker genes derived from the gut metagenomes of 197 termites and one Cryptocercus cockroach, we show that bacteria endemic to termite guts are mostly transferred vertically. We identified 18 lineages of gut bacteria showing cophylogenetic patterns with termites over tens of millions of years. Horizontal transfer rates estimated for 16 bacterial lineages were within the range of those estimated for 15 mitochondrial genes, suggesting that horizontal transfers are uncommon and vertical transfers are the dominant transmission route in these lineages. Some of these associations probably date back more than 150 million years and are an order of magnitude older than the cophylogenetic patterns between mammalian hosts and their gut bacteria. Our results suggest that termites have cospeciated with their gut bacteria since first appearing in the geological record.

Cryptic subterranean diversity: regional phylogeography of the sand termite Psammotermes allocerus Silvestri, 1908 in the wider Namib region

Psammotermes allocerus Silvestri, 1908 is the only described species representing the genus Psammotermes Desneux, 1902 in Southern Africa. The large geographical range of this subterranean termite covers both summer and winter rainfall regimes. Deadwood is the preferred food when available, but in more arid habitats, both live and dead grasses form the major dietary component. Along the Namib Desert margins, the species’ localised herbivory creates circular bare patches known as fairy circles. For a regional phylogeographic study of this species, we sampled 65 sand termite populations within drier parts of Namibia, South Africa, and Angola. Based on combined molecular and ecological data, we found considerable genetic diversification within P. allocerus. Analyses of two mitochondrial markers (COI, COII), including a Bayesian inference tree, haplotype analysis and genetic distances suggest a delineation into seven highly differentiated genetic groups. The ‘Succulent Karoo’ group is additionally characterised by unique features of the royal chamber, nest and tunnel system. In conclusion, our data suggest that P. allocerus should be not regarded as one species but as a species complex. Termites of each analysed group ‘Northern Namib’, ‘Western Kalahari Basin’, ‘Nama’, ‘Southwestern Kalahari’, ‘East Gariep’, ‘Southern Namib’ and ‘Succulent Karoo’ should be considered as distinct species. The species name P. allocerus should be used for termites of the ‘Succulent Karoo’.