Abstract
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.
Similar content being viewed by others
Introduction
Termites, or Isoptera, are a clade of eusocial cockroaches in a sister-group relationship with the subsocial Cryptocercus wood roaches1,2,3,4,5,6. With over 3000 described species7, the diversity of termites is relatively low compared to hyper-diverse insect clades (e.g., Coleoptera or Lepidoptera). Termites achieve critical ecosystem functions, such as decomposing plant material, providing bioturbation of soils, and influencing water infiltration in soils of tropical and subtropical terrestrial ecosystems, where they are amongst the most abundant animals8,9,10,11,12,13,14,15. Termites are also well-known as major pests of wood in service, such as in buildings and constructions, utility poles, fencing, etc.7,16. In regions where they are found, they account for over half of total urban pest control costs17,18,19, excluding costs from insect-vectored diseases. Consequently, public opinion on termites has been shaped by, and the bulk of termite biology research has focused on, a handful of pestiferous termite genera. The combination of moderate species diversity and high abundance, leading to economic and ecological importance, has made termites of considerable interest to insect taxonomists. The discovery of new extant taxa is rare at the supra-generic level (but see Romero Arias et al.20), although it is relatively common at genus and species levels.
The higher classification of termites has undergone substantial changes since its inception more than a century ago21,22,23,24,25,26,27,28,29,30,31. In the Treatise on the Isoptera of the World, the last comprehensive catalogue published on the global termite fauna, extant termites were divided into nine families: Archotermopsidae, Hodotermitidae, Kalotermitidae, Mastotermitidae, Rhinotermitidae, Serritermitidae, Stolotermitidae, Stylotermitidae, and Termitidae7. An additional family was subsequently erected for Hodotermopsis, the Hodotermopsidae32,33, resolving the paraphyly of Archotermopsidae with respect to Hodotermitidae identified by preceding molecular phylogenies3,34,35,36,37. Molecular phylogenetic analyses have identified additional non-monophyletic termite families and subfamilies, indicating that additional augmentations of termite classification are required.
All ten currently recognised extant termite families (sensu Krishna et al.7, modified by Wang et al.33) are monophyletic, except for the Rhinotermitidae, which is paraphyletic with respect to Serritermitidae and Termitidae. The non-monophyletic status of Rhinotermitidae has been evidenced in numerous phylogenetic studies based on morphological characters38,39, a few molecular markers34,40,41, or a combination of both3,35,42. Phylogenetic reconstructions performed with complete mitogenomes36,37 and, more recently, with thousands of nuclear markers obtained from transcriptomes43 or Ultraconserved Elements (UCEs)44, confirmed the paraphyly of Rhinotermitidae. In addition, molecular phylogenies indicate that the rhinotermitid subfamily Heterotermitinae is paraphyletic with respect to Coptotermitinae43,44. A taxonomic revision of these family and subfamily names is required, and phylogenies inferred from many independent nuclear markers, which yield robust reconstructions, provide an ideal framework.
One additional polyphyletic subfamily requires attention within the Termitidae. Nine subfamilies of Termitidae are currently recognised: Apicotermitinae, Cubitermitinae, Engelitermitinae, Foraminitermitinae, Macrotermitinae, Nasutitermitinae, Sphaerotermitinae, Syntermitinae, and Termitinae7,20. All subfamilies were consistently retrieved as monophyletic, with the notable exception of Termitinae, which has consistently been found to be polyphyletic34,37,41,43,44,45. Indeed, the Cubitermitinae, Engelitermitinae, Nasutitermitinae, and Syntermitinae are all nested within the traditional concept of Termitinae. The early diversification of crown Termitidae presumably occurred at a fast pace, leaving termitid lineages separated by short, unresolved internodes, along with significant discordance among individual gene trees43,44, suggestive of incomplete lineage sorting and hybridisation. While the relationships among the main termitid lineages are unclear in many cases, UCE phylogenies may allow for the identification of monophyletic lineages and their elevation to a subfamilial rank.
The purpose of this study is to revise termites at the family and subfamily levels and to produce a robust classification only composed of monophyletic groups. To achieve this, we thoroughly inspect the termite tree of life through the computation of 30 phylogenetic trees reconstructed from concatenated sequence supermatrices and 21 trees reconstructed within the multi-species coalescence framework. Our maximum-likelihood reconstructions use UCEs and diverse nucleotide, codon, and amino acid substitution models. UCEs are ideal markers to build robust phylogenetic trees, as they are highly conserved nuclear markers spread across all chromosomes. In arthropods, most UCEs are found in exons44,46,47,48,49. Our sampling includes 135 termite species, covering the breadth of the termite tree of life, and is especially exhaustive for the Rhinotermitidae and Termitinae. We delineate new circumscriptions for all families and subfamilies to abide by the monophyly criterion and corroborate our taxonomic amendments with analyses of rogue loci and taxa and topological tests. These comprehensive phylogenetic analyses provide a robust framework to create a lasting global taxonomic classification of termites.
Results
Phylogenetic analyses and their implication on the higher classification of termites
Our phylogenetic reconstructions support the monophyly of all currently recognised termite families (sensu Krishna et al.7, modified by Wang et al.33), except Rhinotermitidae, which is paraphyletic with respect to Serritermitidae and Termitidae, as previously acknowledged3,34,35,36,37,38,39,40,41,42,43,44.
The Rhinotermitidae (sensu Krishna et al.7) can be divided into four monophyletic lineages diverging early alongside Serritermitidae and Termitidae (Figs. 1, 2): Heterotermitinae Froggattâ+âCoptotermitinae Holmgren, Psammotermitinae Holmgrenâ+âProrhinotermitinae Quennedey & Deligne, Rhinotermitinae Froggatt, and Termitogetoninae Holmgren. All our analyses supported the monophyly of these four clades, and our analyses with RogueNaRok found no evidence of rogue taxa within these clades. Assigning each of these four clades a familial rank resolves the challenges embodied in Rhinotermitidae as traditionally conceived. Accordingly, we raise these four lineages to familial rank. We restrict Rhinotermitidae to Acorhinotermes, Dolichorhinotermes, Parrhinotermes, Rhinotermes, and Schedorhinotermes, the five genera composing the Rhinotermitinae (i.e., Rhinotermitidae Froggatt, sensu novo). We synonymise the two largely unused names Arrhinotermitinae and Leucotermitinae with Coptotermitinae and Heterotermitinae, respectively. We further synonymise Coptotermitinae with Heterotermitinae. The synonymy of the monogeneric Coptotermitinae (Coptotermes) with Heterotermitinae is justified by the paraphyly of Heterotermitinae, within which Coptotermitinae are nested. Heterotermitinae, having precedence over the other names, is here elevated to familial rank (i.e., elevated as Heterotermitidae Froggatt, stat. nov.). The monogeneric subfamilies Prorhinotermitinae and Psammotermitinae are together raised to familial rank, the latter having precedence (i.e., Psammotermitidae Holmgren, stat. nov.). Lastly, we elevate Termitogetoninae to Termitogetonidae Holmgren, stat. nov. (see Systematics, infra).
For readability, we only provided support values (UFBs, ultrafast bootstraps) for nodes representing families and subfamilies. Node labels are site concordance factors (sCFs) calculated on the majority rule consensus tree and the unfiltered UCE supermatrix. Species whose names are in bold were identified as rogue taxa by RogueNaRok (number of times under parentheses; for details, see Supplementary Data 5). The source data used for this figure can be found in the Dryad repository at https://doi.org/10.5061/dryad.02v6wwqbm (trees 1-30 in âFile 4: trees.tarâ).
The tree presented herein (tree 75 in File 4 on Dryad) was constructed from gene trees for which outliers were pruned with TreeShrink (trees 52â72). The summary topology without pruning outliers is identical for all but one intrasubfamilial node (see tree 74 built from trees 31â51 on Dryad: File 4). For readability, we only provided support values (LPPs, local posterior probabilities) for nodes representing families and subfamilies. Node labels are site concordance factors (sCFs) calculated on the majority rule consensus tree and the unfiltered UCE supermatrix. The source data used for this figure can be found in the Dryad repository at https://doi.org/10.5061/dryad.02v6wwqbm (âFile 4: trees.tarâ).
The phylogenetic relationships of Psammotermitidae and Termitogetonidae, which diverged shortly after the origin of the clade Serritermitidaeâ+âRhinotermitidaeâ+âTermitogetonidaeâ+âPsammotermitidaeâ+âHeterotermitidaeâ+âTermitidae, are unresolved and variable among analyses (Fig. 1). Our phylogenetic analyses placed Termitogetonidae in four potential topological conformations, being sister to: Psammotermitidae (analyses on concatenated data: UFBsâ=â49â93, nâ=â7/30 analyses; coalescence analyses: LPPâ=â0.49, nâ=â1/21); Serritermitidaeâ+âRhinotermitidae (UFBsâ=â60â99, nâ=â13/30; LPPsâ=â0.35â0.54, 6/21); Heterotermitidaeâ+âTermitidae (LPPsâ=â0.49â0.90, nâ=â6/21); and Psammotermitidaeâ+âHeterotermitidaeâ+âTermitidae (UFBsâ=â52â99, nâ=â10/30; LPPsâ=â0.38â0.88, nâ=â8/21). The latter branching was recovered in both majority rule consensus trees, but with sCFs close to 33.3% (32.1â32.5%; Figs. 1, 2), the values obtained for alignments generated at random, indicating the absence of signal for this node in the alignments50. To assess the anomaly zone surrounding the origin of Termitogetonidae and Psammotermitidae, we performed a series of topology tests displacing Termitogetonidae. Topology tests indicated that Termitogetonidae could be sister to non-Stylotermitidae Neoisoptera and Psammotermitidaeâ+âHeterotermitidaeâ+âTermitidae, as reported in some mitogenome-based Bayesian phylogenies51, but it was never in a direct sister relationship with any of these families (AU tests; Supplementary Data 6).
Our sampling of Termitidae included a total of 77 representative genera. Notably, the Syntermitinae were found to be paraphyletic and include the enigmatic monotypic genus Genuotermes Emerson (see Fig. 1). Although sharing key features with the Syntermitinae, Genuotermes was first tentatively assigned to the Termitinae52, based on similar soldier morphologies (Fig. 3 and Supplementary Figs. 17, 18). Genuotermes was identified as a rogue taxon by RogueNaRok in one analysis performed on concatenated sequence data (Fig. 1 and Supplementary Data 5). Genuotermes was also retrieved sister to Onkotermes in one coalescence analysis (Fig. 2; tree 36 in âFile 4: trees.tarâ, available on Dryad at https://doi.org/10.5061/dryad.02v6wwqbm; LPPâ=â0.66), both genera forming a clade outside the Syntermitinae and the Microcerotermes clade. The phylogenetic position of Genuotermes is, therefore, not entirely resolved; however, our analyses generally placed Genuotermes within the Syntermitinae53, and we hereby formally transfer this genus to the subfamily Syntermitinae.
This cladogram integrates topologies summarised from analyses performed on concatenated data and coalescence analyses. The preferred topology within anomaly zones estimated from approximately unbiased topological tests is indicated in dashed grey. The biogeographic matrix indicates the realms occupied by extant members of all isopteran families and termitid subfamilies (adapted from Krishna et al. 7). Here, we recognise nine realms (sensu Holt et al.54): Afrotropical (including Madagascan), Neotropical (incl. Panamanian), Oriental, Australian, Oceanian, Saharo-Arabian, Sino-Japanese, Nearctic, and Palaearctic. For termitid subfamilies, we report the main diagnostic features based on workersâ digestive tubes and soldiersâ heads. These features are summarised from Supplementary Notes 1 and 3, and Supplementary Figs. 1â18. (i) For workersâ midgut-hindgut junction, we indicate: the presence of a mixed segment (square; stippled section), the presence of Malpighian nodules (circle); the presence of a Malpighian knot (full diamond) or pseudo-knot (full diamond with empty circle); the shape of the junction between the mixed segment and the ileum (P1) either elbowed (empty oval) or arched (full oval). For the hindgut, we indicate the presence of a ventral loop (star) formed under the rectum (P5) by P1; and the position (triangle) of the enteric valve (P2, indicated by an arrow) at the insertion of P1 into the paunch (P3). (ii) For termitid soldiers, we indicate the presence of a frontal projection (hexagon) and whether the fontanelle opens at the tip of the projection (full hexagon with empty circle), or at its base (full hexagon). Greyed shapes indicate that both trait states occur within the considered subfamily. Diagnostic features for all subfamilies are extensively presented in Supplementary Figs. 1â18. Soldier pictures, from left to right: (1) Microcerotermes sp. (Microcerotermitinae); (2) Termes fatalis (Termitinae); (3) Silvestritermes heyeri (Syntermitinae), (4) Genuotermes spinifer (Syntermitinae); (5) Nasutitermes sp. (Nasutitermitinae). Picture credits: M. M. Rocha (4); R. H. Scheffrahn (2,3); J. Å obotnÃk (1,5).
Our phylogenetic analyses indicate that all other subfamilies of Termitidae are monophyletic, with the exception of the polyphyletic Termitinae, within which four other subfamilies are nested: Cubitermitinae, Engelitermitinae, Nasutitermitinae, and Syntermitinae. The relationships within the Termitinaeâ+âCubitermitinaeâ+âEngelitermitinaeâ+âNasutitermitinaeâ+âSyntermitinae clade are largely unresolved. This clade is composed of a dozen monophyletic lineages with uncertain relationships and separated by exceedingly short internodes (see Figs. 1, 2). Based on our phylogenetic reconstructions, we recognise 18 major monophyletic lineages within Termitidae, eight of which already have subfamilial status. Giving subfamilial ranks to the other ten lineages resolves the polyphyly of Termitinae while simultaneously retaining the subfamilial status of four subfamilies supported as monophyletic: Cubitermitinae (UFBsâ=â100; LPPsâ=â1; sCFsâ=â64â64.9%), Engelitermitinae, Nasutitermitinae (UFBsâ=â100; LPPsâ=â0.99â1; sCFsâ=â52.1â53%), and Syntermitinae (UFBsâ=â100; LPPsâ=â0.92â1; sCFsâ=â62.9â64%). Accordingly, we redefine Termitinae Latreille sensu novo and restrict it to Termes and related genera, which unambiguously formed a monophyletic group in all our analyses (UFBsâ=â100; LPPsâ=â1; sCFsâ=â56.5â57%). We advocate the resurrection of three monophyletic subfamilies: Microcerotermitinae Holmgren stat. rev. (UFBsâ=â100; LPPsâ=â0.74â1; sCFsâ=â48.2â49.7%), Mirocapritermitinae Kemner stat. rev. (UFBsâ=â99â100; LPPsâ=â1; sCFsâ=â55â56.1%), and Amitermitinae Kemner stat. rev. (UFBsâ=â99â100; LPPsâ=â0.98â1; sCFsâ=â45.8â46.9%). The aforementioned subfamilies were found monophyletic across all analyses, with the exception of Microcerotermitinae, as Onkotermes was placed sister to Microcerotermesâ+âSyntermitinae in two coalescence analyses (LPPsâ=â1; trees 34, 50 in âFile 4: trees.tarâ, available on Dryad at https://doi.org/10.5061/dryad.02v6wwqbm) and sister to Genuotermes in one coalescence analysis (see above). Although not entirely supported by our phylogenetic analyses, the placement of Onkotermes within Microcerotermitinae is also justified by their similar gut morphologies, with a characteristic elbowed junction between the mixed segment and the ileum (Fig. 3). Our topology tests did not reject Microcerotermitinaeâ+âSyntermitinae as sisters to other non-Sphaerotermitinae, non-Macrotermitinae, non-Foraminitermitinae and non-Apicotermitinae Termitidae (Supplementary Data 6), while coalescence analyses placed the clade on a polytomy in 8 of 21 coalescence analyses run with diverse sequence matrices and substitution models (Fig. 2).
Lastly, we propose to erect six monophyletic lineages to new subfamilies in Termitidae: Crepititermitinae subfam. nov., Cylindrotermitinae subfam. nov. (UFBsâ=â78â100; LPPsâ=â0.70â1; sCFsâ=â44.4â44.7%), Forficulitermitinae subfam. nov., Neocapritermitinae subfam. nov. (UFBsâ=â100; LPPsâ=â1; sCFsâ=â76.6â77%), Protohamitermitinae subfam. nov., and Promirotermitinae subfam. nov. (UFBsâ=â100; LPPsâ=â1; sCFsâ=â77.4â79.5%). These six subfamilies were placed on a polytomy, reflecting their uncertain phylogenetic positions, varying among analyses. The monotypic subfamilies Crepititermitinae and Forficulitermitinae were identified as rogue taxa in, respectively, two and ten analyses performed on concatenated sequence data from edge-trimmed protein-coding loci (Fig. 1 and Supplementary Data 5). Our topology tests did not reject Engelitermitinae and Forficulitermitinae as forming a clade sister to other non-Sphaerotermitinae, non-Macrotermitinae, non-Foraminitermitinae, non-Apicotermitinae, non-Microcerotermitinae and non-Syntermitinae Termitidae (AU tests; Supplementary Data 6). The Crepititermitinae were retrieved as sister to Cylindrotermitinaeâ+âProtohamitermitinae in 18 of 30 analyses performed on concatenated sequence data (Fig. 1; UFBsâ=â51â87), a scenario absent in our coalescence analyses (Fig. 2). The uncertain phylogenetic positions of Crepititermitinae and Forficulitermitinae, together with their unique combinations of morphological and biological characteristics (see the Systematics section below; Fig. 3), justified their placements in two separate subfamilies.
Systematics
We present the revised classification of extant Isoptera in Table 1. The diagnoses and distribution maps of Neoisoptera families and subfamilies presented herein are based on and modified from Krishna et al.7. The biogeographic boundaries follow Holt et al.54. The descriptions of digestive tubes follow Noirot55.
Family-level classification changes for Neoisoptera
In the following, we provide abridged diagnoses for families with changes. Please refer to Supplementary Note 1 for extended diagnoses of all neoisopteran families and Supplementary Note 2 for a revised identification key to all families of Isoptera.
Family Rhinotermitidae Froggatt, sensu novo
The subfamily Rhinotermitinae Froggatt21, 1897: 518, is herein elevated to the familial rank.
Type genus: Rhinotermes Hagen56, 1858.
Included genera (5): Acorhinotermes Emerson, Dolichorhinotermes Snyder & Emerson, Parrhinotermes Holmgren, Rhinotermes Hagen, and Schedorhinotermes Silvestri (Macrorhinotermes Holmgren, syn. nov.).
Diagnosis: Imagoes with inflated clypeus, slightly or not at all drawn out over the labrum âreminiscent of a ânose-likeâ projection; and a narrow groove running from the fontanelle to the clypeus. Imago-worker left mandible with an apical tooth, followed by three distinct (unequal) marginal teeth. Imagoes, workers, and soldiers with anterior margin of postclypeus convex. Soldiers without setae surrounding the fontanelle. Soldier labrum elongated, forming a brush apically and endowed with a groove in the centre along which the chemical secretion of the frontal gland flows from the fontanelle.
Remarks: Macrorhinotermes was initially described as a subgenus of Rhinotermes57. Snyder26 considered Macrorhinotermes as a potential synonym of Schedorhinotermes but did not unambiguously synonymize the two names, and Macrorhinotermes has remained a valid genus name. Here, we formally synonymize Macrorhinotermes and Schedorhinotermes, the latter having precedence over the former. There is an interesting diversity of soldier morphologies among Rhinotermitidae. Soldiers can be monomorphic as in Acorhinotermes and Parrhinotermes; dimorphic as in Dolichorhinotermes, Rhinotermes, and most Schedorhinotermes; or trimorphic as in some Schedorhinotermes. Rhinotermitidae occur in the Australian, Afrotropical, Neotropical, Oriental, Palaearctic (Sino-Japanese), and Oceanian realms.
Family Termitogetonidae Holmgren, stat. nov
The subfamily Termitogetoninae Holmgren23, 1910: 286, is herein elevated to the familial rank.
Type genus: Termitogeton Desneux58, 1904
Included genus: monogeneric.
Diagnosis: Imagoes densely hairy. Head subtriangular with posterior margin sinuate. Pronotum small, with median projections from anterior margins. Soldiers are densely hairy with a heart-shaped head. Soldier mandibles elongated without marginal teeth. Soldier legs with two protibial apical spurs. All castes dorso-ventrally flattened.
Remarks: Termitogeton is a wood-feeding termite known from the Oriental and Oceanian realms. It is composed of only two known extant species.
Family Psammotermitidae Holmgren, stat. nov
The family Psammotermitidae regroups the Psammotermitinae and Prorhinotermitinae subfamilies.
Subfamily Psammotermitinae Holmgren24, 1911: 64.
Type genus: Psammotermes Desneux59, 1902
Subfamily Prorhinotermitinae Quennedey & Deligne60, 1975: 265.
Type genus: Prorhinotermes Silvestri61, 1909.
Included extant genera (2): Psammotermes Desneux and Prorhinotermes Silvestri.
Diagnosis: Imagoes with pronotum flat; tarsi tetramerous. Wings with all veins arising independently from inside the wing scales. Nymphs with fused wing buds forming two shield-looking structures on the mesothorax and metathorax. Imago-worker left mandible with an apical tooth followed by three distinct marginal teeth. Soldiers with tibial spurs formula 3:2:2; presenting a more or less marked groove from the fontanelle to clypeus60.
Remarks: The extinct genus â Zophotermes Engel, 2011 is included in the Psammotermitidae. The distribution of Prorhinotermes is mostly insular. This distribution has fuelled speculation that colonies can be transported in driftwood by ocean currents, thus accounting for such a wide distribution on islands of the South Pacific. In sharp contrast, Psammotermes is a desert dweller occurring in the Afrotropical (incl. Madagascan), Saharo-Arabian and Oriental realms.
Clade Geoisoptera Engel, Hellemans, & Bourguignon, nov
This is an important group of termites as it includes the Termitidae, the most species-rich family of Isoptera, and the Heterotermitidae, which include some of the most widespread and critical pest species. In this clade, the number of chromosomes is fixed to 2nâ=â42 (reviewed in Jankásek et al.62), and their mitogenomes are characterised by a complex repeat structure in their control region36. Additionally, all species of this clade have a true worker caste (as in Mastotermitidae, Hodotermitidae, and Rhinotermitidae sensu novo) and typically live in large colonies.
Etymology: The clade name is a combination of the Ancient Greek prefix geoâ (γεÏâ, from gê / γá¿, meaning âsoil,â âland,â or âearthâ) and Isoptera. The name references the growing reliance on organically rich soil as a source of nitrogen63.
Family Heterotermitidae Froggatt, stat. nov
The subfamily Heterotermitinae Froggatt21,1897: 550, is herein elevated to the familial rank.
Heterotermitinae Froggatt21, 1897: 550. Type genus: Heterotermes Froggatt21.
Leucotermitinae Holmgren23, 1910: 285. Type genus: Leucotermes Silvestri64, 1901 [=âHeterotermes Froggatt21, 1897].
Coptotermitinae Holmgren23, 1910: 285. Type genus: Coptotermes Wasmann65, 1896. Syn. nov.
Arrhinotermitinae Sjöstedt66, 1926: 8. Type genus: Arrhinotermes Wasmann67, 1902 [=âCoptotermes Wasmann]. Syn. nov.
Type genus: Heterotermes Froggatt21, 1897
Included genera (3): Coptotermes Wasmann, Heterotermes Froggatt, and Reticulitermes Holmgren. â Lukotermes Perkovsky & Nel, syn. nov. is an extinct member of Heterotermitidae.
Diagnosis: Imagoes with postclypeus devoid of nose-like projection or groove extending from the fontanelle. Head round in Coptotermes and oval with even posterior margin in Heterotermes and Reticulitermes. Pronotum without median projection on the anterior margin; tarsi tetramerous. Workers and soldiers with anterior margin of postclypeus flat or concave. Soldiers with fontanelle surrounded by setae and katana-shaped mandibles devoid of marginal teeth (though with basal serrations); protibia with more than three protibial apical spurs. Soldiers of Coptotermes displaying an egg-shaped head with a large fontanelle directed forward and oozing latex-appearing secretions. Soldiers of Heterotermes and Reticulitermes exhibit a rectangular head with a small indistinct fontanelle.
Remarks: Some species exhibit slightly dimorphic soldiers68. The Heterotermitidae now encompass the three genera often referred to as âsubterranean termitesâ, a term historically associated with âRhinotermitidaeâ with an extensive pest status in the literature. The Heterotermitidae occur across most temperate and tropical regions.
Subfamily-level classification changes for Termitidae Latreille
The number of subfamilies within the Termitidae is expanded given that the current Termitinae are demonstrably polyphyletic. We describe the new subfamilies of Termitidae erected in this paper and redescribe the Termitinae. We do not address the Syntermitinae, as the inclusion of Genuotermes does not alter the diagnosis of the subfamily. We do not address the other seven subfamilies of Termitinae, as they remain unaltered (i.e., Sphaerotermitinae Engel & Krishna; Macrotermitinae Kemner; Foraminitermitinae Holmgren; Apicotermitinae Grassé & Noirot; Engelitermitinae Romero Arias, Roisin, & Scheffrahn; Nasutitermitinae Hare; and Cubitermitinae Weidner). Diagnoses for these subfamilies can be found in Krishna et al.7 and Romero Arias et al.20. An updated key for all subfamilies of Termitidae is given in Supplementary Note 3. An updated list of all genera of extant Isoptera currently recognised is given in Supplementary Note 4. A list of the top 20 UCEs diagnostic for each subfamily of Termitidae (excluding subfamilies represented by only one sample, i.e., Crepititermitinae, Forficulitermitinae, and Protohamitermitinae) is provided in Supplementary Data 7. An accompanying diagnostic tool is available on the Termite UCE Database (in âFile 5: termitidae_diagnosing_database_v1.fasta.gzâ available on Dryad at https://doi.org/10.5061/dryad.02v6wwqbm; also available at: https://github.com/sihellem/TER-UCE-DB/). Diagnostic features for each of the 18 subfamilies are given in Supplementary Figs. 1â18.
Subfamily Amitermitinae Kemner, stat. nov
Amitermitinae Kemner69, 1934: 110.
Type genus: Amitermes Silvestri64, 1901.
Included genera (14): Ahamitermes Mjöberg, Amitermes Silvestri, Dentispicotermes Emerson, Drepanotermes Silvestri, Eremotermes Silvestri, Globitermes Holmgren, Gnathamitermes Light, Hoplotermes Light, Incolitermes Gay, Invasitermes Miller, Orthognathotermes Holmgren, Prohamitermes Holmgren, Pseudhamitermes Holmgren, and Synhamitermes Holmgren.
Diagnosis: Imago-worker left mandible with apical tooth (LAt) shorter than marginal teeth (LMt), except in a few genera such as: Ahamitermes and Invasitermes70, or Orthognathotermes (see Supplementary Fig. 1i). Worker gut with ileum (P1) tubular, dilated, or globular, without loop ventrally; each pair of Malpighian tubules united in one base (or nodule), except in Prohamitermes in which each tubule ends in its own nodule. Soldiers generally have a pear-shaped head and are devoid of frontal projection (except in Eremotermes and Dentispicotermes). Diagnostic UCE loci are listed in Supplementary Data 7. Diagnostic features are given in Supplementary Fig. 1.
The subfamily is most similar to the Promirotermitinae, whose enteric valve (P2) position, i.e. the insertion of P1 into the paunch (P3), occurs on the right side of the abdomen (rather than the left side, as in Microcerotermitinae, Mirocapritermitinae, and Neocapritermitinae; Fig. 3). The Amitermitinae can be distinguished from the promirotermitines by: P1 tubular to inflated (always inflated in Promirotermitinae); soldier mandibles of various forms, ranging from crushing, piercing, slashing, or symmetrical snapping (always symmetrical snapping in Promirotermitinae); and soldier head pear-shaped without frontal projection (a notable frontal hump or projection is present in Promirotermitinae).
Remarks: Amitermitinae are widespread across every tropical and subtropical biogeographic realm. Species feed on various substrates, including mounds, soil, wood (incl. dry grasses), and the soil-wood interface. The genus Invasitermes is reported as soldierless.
Subfamily Crepititermitinae Hellemans, Engel, & Bourguignon, subfam. nov
Zoobank registration: https://zoobank.org/urn:lsid:zoobank.org:act:5F7D1841-7543-4914-B552-10318E2224C2.
Type genus: Crepititermes Emerson71, 1925.
Included genera: monogeneric.
Diagnosis: Imago-worker left mandible lacks a third marginal tooth (LMt3). Right mandible lacking second marginal tooth (RMt2). Worker gut with a blind diverticulum or caecum on P3; Malpighian tubules forming two pairs attached individually. Soldiers with subrectangular heads, labrums with concave apical margins and sharp apicolateral points, flattened frons, and symmetrical snapping mandibles. Diagnostic features are given in Supplementary Fig. 3.
The presence of a blind diverticulum on P3 and the absence of an indentation between LMt3 and LMt1+2 in Crepititermitinae are reminiscent of the Cubitermitinae, but the former can be distinguished by: the presence of a columnar belt (of the generalised type) in the gizzard (absent in Cubitermitinae); P1 tubular (P1 inflated in Cubitermitinae); the absence of intermediate cushions on the enteric valve armature (intermediate cushions present in Cubitermitinae; see Supplementary Fig. 4j, k); and the blind diverticulum on P3 occurring dorsally (ventrally in Cubitermitinae).
Remarks: Crepititermes was previously assigned to the Termes-group. A comprehensive re-description of Crepititermes was provided by Rocha & Cuezzo72. The genus is exclusively Neotropical, with species that feed on humus. Crepititermes workers use an autothysis mechanism that involves the rupture of the upper side of the abdomen, near the junction with the thorax, to expose their sticky gut for colony defence.
Subfamily Cylindrotermitinae Hellemans, Engel, & Bourguignon, subfam. nov
Zoobank registration: https://zoobank.org/urn:lsid:zoobank.org:act:DEEDCE51-6696-4C97-BB2E-654EF86C8F9F.
Type genus: Cylindrotermes Holmgren73, 1906.
Included genera (2): Cephalotermes Silvestri and Cylindrotermes Holmgren.
Diagnosis: Imago-worker left mandible with apical tooth (LAt) roughly equal or shorter than LMt1+2 and LMt3 well-developed. Worker gut with long, P1 tubular, forming a loop ventrally under the rectum (P5). Enteric valve bearing six unequal cushions. Soldier head elongated, subrectangular, devoid of frontal projection, covered by short bristles; mandibles short, robust, of biting/slashing type; weak humps on the frons located immediately between the antennae above the basolateral corners of clypeus. Diagnostic UCE loci are listed in Supplementary Data 7. Diagnostic features are given in Supplementary Fig. 5.
The Cylindrotermitinae can be differentiated from the Amitermitinae as in the latter, P1 forms a loop ventrally under P5 and P2 occurs on the left side of the abdomen (Fig. 3).
Remarks: The two genera of Cylindrotermitinae were previously assigned to the Amitermes-group sensu Krishna et al.7. Cylindrotermitinae are wood-feeding termite genera from separate biogeographic regions: Cephalotermes is Afrotropical, while Cylindrotermes is Neotropical. Cephalotermes builds mounds, and forms populous colonies with few soldiers, headed by one king and one extremely physogastric queen that loses the ability to crawl as the colony matures (J. Å obotnÃk, pers. obs.). Cylindrotermes are found in dead wood and form small colonies.
Subfamily Forficulitermitinae Hellemans, Engel, & Bourguignon, subfam. nov
Zoobank registration: https://zoobank.org/urn:lsid:zoobank.org:act:429D8F1D-1D78-44AD-BF91-4AAE9492251F.
Type genus: Forficulitermes Emerson74, 1960.
Included genera: monogeneric.
Diagnosis: Imago-worker right mandible with well-developed second marginal tooth (RMt2). Worker gut with four Malpighian tubules attached in pairs at the junction of the mixed segment. Soldier head subquadrangular with fontanelle located in a small depression, with two flat longitudinal ridges located in front of the opening on each side and forming a faint groove-like structure leading to the base of the clypeus; small, elevated humps on each side between the fontanelle and antennal base; mandibles elongate and slender, without marginal teeth, used to bite and unable to snap; apical margin of the labrum with medioapical convexity and sharply pointed apicolateral angles. Diagnostic features are given in Supplementary Fig. 8.
The Forficulitermitinae can be distinguished from termitines in their newly restricted circumscription by the absence of a Malpighian nodule (generally present in Termitinae sensu novo). In addition, the soldier mandibles are elongated and not of snapping type, while in Termitinae the mandibles are generally snapping (except in three genera with biting/slashing mandibles: Divinotermes, Spinitermes, and Tuberculitermes). Lastly, the soldier head has a faint groove running from the fontanelle, which sits in a shallow depression, to the clypeal base and with small humps on either side of the channel between the fontanelle and the antennal toruli, while termitines lack such a groove, depression, or lateral humps and instead have a medial frontal prominence of varying degrees of development.
Remarks: Forficulitermes was previously assigned to the Termes-group sensu Krishna et al.7. It is an exclusively Afrotropical genus, including species that feed on soil. A comprehensive re-description of Forficulitermes was provided by Scheffrahn & KÅeÄek75.
Subfamily Microcerotermitinae Holmgren, stat. nov
Microcerotermitinae Holmgren76, 1910: 145.
Type genus: Microcerotermes Silvestri64, 1901.
Included genera (2): Microcerotermes Silvestri and Onkotermes Constantino.
Diagnosis: Imago-worker left mandible with third marginal tooth (LMt3) well developed, right mandible with posterior margin of LMt2 long and straight. Worker gut with a conspicuous Malpighian knot; junction of the mixed segment with P1 constricted and elbowed, not fully visible in ventral view; P1 tubular, dilated anteriorly, constricted posteriorly at the connection with P2 on the left side of the abdomen, without loop ventrally. P2 is constricted and poorly ornamented. Soldier devoid of frontal projection, head generally long and narrow. Imago head subrectangular. Diagnostic UCE loci are listed in Supplementary Data 7. Diagnostic features are given in Supplementary Fig. 10.
Microcerotermitines are similar to Amitermitinae but have P1 inserted on P3 on the left side of the abdomen rather than the right in the latter clade (Fig. 3). In addition, the soldier head of Microcerotermitinae is subrectangular rather than the pear-shaped form of Amitermitinae.
Remarks: Species of this group were previously assigned to the Amitermes-group but were found to consistently form a sister relationship with Syntermitinae. The microcerotermitines are pantropical wood feeders (sensu lato: wood and grass).
Subfamily Mirocapritermitinae Kemner, stat. nov
Mirocapritermitinae Kemner69, 1934: 166.
Type genus: Mirocapritermes Holmgren76, 1914.
Included genera (14): Dicuspiditermes Krishna, Homallotermes John, Indocapritermes Chhotani, Kemneritermes Ahmad & Akhtar, Krishnacapritermes Chhotani, Labiocapritermes Krishna, Mirocapritermes Holmgren, Oriencapritermes Ahmad & Akhtar, Pericapritermes Silvestri, Procapritermes Holmgren, Pseudocapritermes Kemner, Rinacapritermes Amina & Rajmohana, Sinocapritermes Ping & Xu, and Syncapritermes Ahmad & Akhtar.
Diagnosis: Imago-worker right mandible with anterior margin of RMt2 convex or straight, posterior margin concave. Worker gut with Malpighian tubules attached in two separate pairs fused proximally, devoid of Malpighian nodules. Soldiers with rectangular heads and snapping mandibles are generally asymmetrical (symmetrical in some species of Procapritermes). Soldier head devoid of frontal projection (except in Mirocapritermes). Diagnostic UCE loci are listed in Supplementary Data 7. Diagnostic features are given in Supplementary Fig. 11.
Mirocapritermitinae can be distinguished from the termitines in their restricted sense by the absence of a Malpighian nodule (present in Termitinae sensu novo) and the absence of a frontal projection (generally present in Termitinae). In the imago-worker mandible, RMt2 is distinctly present and with a concave posterior margin, while it is absent or vestigial with a straight posterior margin in Termitinae sensu novo.
Remarks: Mirocapritermitinae include species previously assigned to the Pericapritermes-group. Their distribution is restricted to the Oriental region, with the exception of Pericapritermes, which is also found in the Australian, Afrotropical, Palaearctic (Sino-Japanese), and Oceanian regions. Mirocapritermitines feed on soil or at the soil-wood interface.
Subfamily Neocapritermitinae Hellemans, Engel, & Bourguignon, subfam. nov
Zoobank registration: https://zoobank.org/urn:lsid:zoobank.org:act:B9567537-33EF-4B4D-992F-34E0C16B8701.
Type genus: Neocapritermes Holmgren25, 1912.
Included genera (3): Neocapritermes Holmgren, Planicapritermes Emerson, and Schievitermes Roisin.
Diagnosis: Imago-worker left mandible with apical tooth (LAt) shorter than LMt1+2; right mandible with RMt2 slightly reduced, posterior margin straight. Worker gut with two pairs of Malpighian tubules fused anteriorly and inserted in a single point at the mixed segment junction, devoid of Malpighian nodule; mixed segment long, distally bilobed into mesenteric sacs in Neocapritermes and Schievitermes (see Roisin77); junction of the mixed segment with P1 constricted and elbowed; P3 inflated, with wall bearing numerous spines. Soldiers with asymmetrical snapping mandibles and subrectangular heads devoid of frontal projection. Diagnostic UCE loci are listed in Supplementary Data 7. Diagnostic features are given in Supplementary Fig. 13.
Like the Termitinae, the imago-worker mandible has RMt2 reduced and with a straight posterior margin but can be distinguished by the absence of a Malpighian nodule (present in Termitinae sensu novo) and the absence of a frontal prominence on the soldier head (generally present to varying degrees in Termitinae). The Neocapritermitinae differ from the Old World Mirocapritermitinae in three ways: the posterior margin of RMt2 is straight (rather than concave in mirocapritermitines); the molar area has distinct ridges (molar area smooth to faint and reduced striations in mirocapritermitines); the enteric valve armature is composed of three folds (six in mirocapritermitines). The Neocapritermitinae exhibit an elbowed midgut-P1 junction as the Microcerotermitinae but can be distinguished by their soldiers (snapping in Neocapritermitinae, piercing/slashing in Microcerotermitinae).
Remarks: The species of Neocapritermitinae were previously assigned to the Termes-group. The subfamily is exclusively Neotropical, with species that feed on decayed wood or at the wood-soil interface.
Subfamily Promirotermitinae Hellemans, Engel, & Bourguignon, subfam. nov
Zoobank registration: https://zoobank.org/urn:lsid:zoobank.org:act:1F08A9B3-6F94-498C-A8B4-78ABF508FDA6.
Type genus: Promirotermes Silvestri78, 1914.
Included genera (2): Angulitermes Sjöstedt and Promirotermes Silvestri.
Diagnosis: Imago-worker right mandible with well-developed RMt2. Worker gut with two pairs of Malpighian tubules fused in a common base. Soldier head is subrectangular to egg-shaped (Promirotermes), with symmetrical snapping mandibles and either a pointed frontal projection (Angulitermes) or a blunt frontal hump (Promirotermes). Diagnostic features are given in Supplementary Fig. 14.
The subfamily can be distinguished from Termitinae sensu novo by its imago-worker mandible with well-developed RMt2, which is otherwise absent or at least considerably reduced in termitines. The subfamily is most easily distinguished from Amitermitinae by its subrectangular soldier heads with a frontal process (generally absent in Amitermitinae) and bifurcate labrum (rounded, or with slightly pointed tips in Amitermitinae).
Remarks: The species of this group were previously assigned to the Termes-group (sensu Krishna et al.7) and are humus-feeders. Promirotermes is exclusively Afrotropical, while Angulitermes is also distributed across the Oriental region, the Middle East, and Africa.
Subfamily Protohamitermitinae Hellemans, Engel, & Bourguignon, subfam. nov
Zoobank registration: https://zoobank.org/urn:lsid:zoobank.org:act:2C4DDA82-D86C-40CA-BCB9-F5CD5AA60E7B.
Type genus: Protohamitermes Holmgren25, 1912.
Included genera (2): Orientotermes Ahmad and Protohamitermes Holmgren.
Diagnosis: Imago-worker right mandible with well-developed subsidiary tooth at the base of the first and second marginal teeth (RMt1 and RMt2), with anterior margin angular and posterior margin curved. Worker gut with two pairs of Malpighian tubules connected to the mixed segment by two nodules, one for each pair; with long, tubular ileum (P1), forming a loop ventrally under the rectum (P5) (as in Cylindrotermitinae and part of the Apicotermitinae). Diagnostic features are given in Supplementary Fig. 15.
The Protohamitermitinae can be distinguished from the Amitermitinae by the presence of a long P1 that loops ventrally under P5 (ventral loop absent in Amitermitinae) and by the presence of a subsidiary tooth on the right imago-worker mandible (absent in Amitermitinae). The subfamily differs from Cylindrotermitinae by the presence of the subsidiary tooth and the absence of soldiers.
Remarks: The Protohamitermitinae include two Oriental soldierless species feeding on wood. They were previously assigned to the Amitermes-group sensu Krishna et al.7. Although the type-genus Protohamitermes was not included in our analyses, its relationship with Orientotermes is known from previous phylogenetic analyses3.
Subfamily Termitinae Latreille, sensu novo
Termitina Latreille80, 1802: 293. Type genus: Termes Linnaeus81.
Mirotermitini Weidner82, 1956: 99. Type genus: Mirotermes Wasmann83, 1897 [=âTermes Linnaeus81, 1758].
Capritermitini Weidner82, 1956: 100. Type genus: Capritermes Wasmann83, 1897.
Type genus: Termes Linnaeus81, 1758
Included genera (23): Apsenterotermes Miller, Capritermes Wasmann, Cavitermes Emerson, Cornicapritermes Emerson, Cristatitermes Miller, Dihoplotermes Araujo, Divinotermes Carrijo & Cancello, Ekphysotermes Miller, Ephelotermes Miller, Hapsidotermes Miller, Hesperotermes Gay, Inquilinitermes Mathews, Lophotermes Miller, Macrognathotermes Silvestri, Palmitermes Hellemans & Roisin, Paracapritermes Hill, Protocapritermes Holmgren, Quasitermes Emerson, Saxatilitermes Miller, Spinitermes Wasmann, Termes Linnaeus, Tuberculitermes Holmgren, and Xylochomitermes Miller.
Diagnosis: Imago-worker right mandible with RMt2 reduced or absent and posterior margin straight. Left mandible with apical tooth (LAt) longer or slightly shorter than LMt1+2. Worker gut with two pairs of Malpighian tubules attached on a nodule connected to the mixed segment. P2 is on the left side of the abdomen (except in Cornicapritermes, whose P1 forms a ventral loop and P2 is on the right; M.M. Rocha and T.F. Carrijo, pers. obs.). Soldiers with head subrectangular, a frontal prominence varying in size among genera, and snapping mandibles, except in Divinotermes and Spinitermes, which exhibit biting/slashing mandibles; and Tuberculitermes, which has snapping/slashing mandibles (A. BuÄek, pers. obs.). Diagnostic UCE loci are listed in Supplementary Data 7. Diagnostic features are given in Supplementary Fig. 18. Additional diagnoses are provided in the systematic sections of Promirotermitinae, Neocapritermitinae, Mirocapritermitinae, and Forficulitermitinae (above).
Remarks: The Termitinae have a pantropical distribution, found in the Afrotropical, Australian, Neotropical, Oceanian and Oriental regions. Termitinae comprises species feeding on diverse substrates, including soil, mounds, wood, and the soil-wood interface.
Discussion
The Rhinotermitidae and Termitinae, as defined before this study, have been recognised as polyphyletic groups since the advent of the molecular phylogenetics era3,37,38,40. Our robust phylogenies inferred with UCE data confirm that, unlike all other families and subfamilies of termites, the Rhinotermitidae and Termitinae were non-monophyletic groups in need of a taxonomic revision. Herein, we carried out the needed taxonomic revision and fixed the termite classification. Our comprehensive phylogenetic analyses unambiguously identified four neoisopteran family-level clades diverging alongside the Serritermitidae and Termitidae: the Heterotermitidae (Coptotermes, Heterotermes, and Reticulitermes), Psammotermitidae (Psammotermes and Prorhinotermes), Termitogetonidae (Termitogeton), and Rhinotermitidae (now restricted to Acorhinotermes, Dolichorhinotermes, Parrhinotermes, Rhinotermes, and Schedorhinotermes). Our analyses also recognised ten subfamily-level clades forming monophyletic groups together with the Cubitermitinae, Engelitermitinae, Nasutitermitinae, and Syntermitinae: i.e., the Amitermitinae, Crepititermitinae, Cylindrotermitinae, Forficulitermitinae, Microcerotermitinae, Mirocapritermitinae, Neocapritermitinae, Promirotermitinae, Protohamitermitinae, and Termitinae (see Table 1). The erection of these new families and subfamilies enriches the termite nomenclature and allows for a predictive classification entirely comprised of monophyletic groups at the supra-generic level.
One potential concern of this newly proposed termite classification is the inflation of the number of family and subfamily names. However, the alternative strategy of synonymizing a number of families and subfamilies has important caveats. Synonymizing the Rhinotermitidae, Serritermitidae, and Termitidae would require downgrading the Rhinotermitidae and Termitidae to subfamily ranks. This would profoundly affect the nomenclature because the terms Rhinotermitidae and Termitidae have been widely used, and both Rhinotermitinae and Termitinae are already valid subfamilies. Such a downgrade would bring unnecessary confusion and undermine our goal to resolve paraphyly in Neoisoptera. Similarly, synonymizing the Cubitermitinae, Engelitermitinae, Nasutitermitinae, Syntermitinae, and Termitinae would strongly affect nomenclature, especially due to the loss of the Nasutitermitinae Hare84. The Nasutitermitinae are the largest subfamily of Termitidae with the obvious synapomorphy of soldiers bearing conical nasi, and this name has been widely accepted and used ever since its erection. Therefore, expanding family and subfamily names largely maintain the stability of the termite classification while also avoiding confusion in the vast termite literature. Our classification preserves the stability of termite nomenclature by maintaining some key families and subfamilies that are featured in the literature, and synonymizing them would have negative impacts on many disciplines.
Despite the above-listed advantages, favouring splitting the formerly polyphyletic Rhinotermitidae and Termitinae over lumping also has consequences. Many genera previously placed in the Rhinotermitidae and Termitinae now belong to different families and subfamilies, requiring updates to databases and specimen collections. The new proposed termite classification will also have impacts outside the scientific community. For example, âRhinotermitidaeâ is well known in the pest control industry for containing the two dominant genera with urban pest species, Coptotermes and Reticulitermes, which now belong to the Heterotermitidae. While preserving the stability of termite nomenclature, our new termite classification will require termitologists at large to adjust their vocabulary.
Increasing the number of families and subfamilies of Neoisoptera is justifiable considering the distribution of species amongst the main termite lineages. Termite diversity is disproportionately biased toward the Neoisoptera and, more specifically, toward the Termitidae, which include approximately 75% of all extant termite species7. Thus, the classification of Termitidae has received more attention than that of the smaller families. Similarly, it justifies that the systems of families within Neoisoptera and subfamilies within Termitidae are enriched compared to less speciose lineages. Several new subfamilies of Termitidae are species-poor lineages. While some of these subfamilies are widespread, such as the cylindrotermitine genera Cephalotermes and Cylindrotermes in the Afrotropical and Neotropical regions, respectively; the Afrotropical Engelitermitinae and Forficulitermitinae are rare monotypic subfamilies75. Notably, the recent description of Engelitermes zambo, the only representative of Engelitermitinae20, indicates that additional termitid subfamilies may be erected in the future to accommodate new species that cannot be placed in existing subfamilies.
Uncertainties in our phylogenetic reconstructions and the morphological diversity of the Neoisoptera were the principal challenges to delimiting new families and subfamilies. The exact phylogenetic relationships amongst the families of Neoisoptera and subfamilies of Termitidae are often unclear. Although we used thousands of UCE markers44, we failed to resolve the phylogenetic relationships amongst some families of Neoisoptera and many subfamilies of Termitidae because the internodes separating subfamilies are very short. This suggests insufficient time for alleles to coalesce, leading to incomplete lineage sorting85,86. Although future phylogenetic analyses may resolve some of these relationships, uncertainties will probably remain due to the apparently rapid radiations within the Termitidae. In addition, some of the new subfamilies of Termitidae include species that do not share clear morphological synapomorphies, complicating their delimitation based on morphological characters. Consequently, some new subfamilies are morphologically diverse. To ensure the robustness of the proposed classification of termites, each subfamily we described is supported as a monophyletic group by all our phylogenetic analyses. The placements of Genuotermes and Onkotermes in the Syntermitinae and Microcerotermitinae, respectively, represent the only two exceptions, as a few analyses placed them outside these subfamilies. Although these placements are favoured in our analyses, we cannot entirely rule out alternative scenarios. We acknowledge that the presently proposed classification is likely to be amended in the future; however, we expect it to serve as a robust backbone toward a definitive classification of termites.
Methods
Biological samples and sequencing
We gathered genome assemblies of 135 termites and three wood-feeding subsocial cockroaches from the genus Cryptocercus to reconstruct the phylogenies. We used previously published UCE datasets44,87, which we supplemented with sequences of 56 samples (see Supplementary Data 1). Our dataset includes representative species from all genera of Rhinotermitidae sensu lato and every subfamily of Termitidae, including all the new subfamilies described here. Within the Termitidae, we included representatives of all termitine lineages that arose shortly after the origin of the first Termitinae.
All samples reported in this study were housed in collections, and no new field collection was performed. Specimens from collections were obtained in accordance with local and worldwide regulations at the time of their collection. Specifically, this involved the terms of the following permits: Australia, no sampling permits required; Bolivia, No. MMAyA-VMABCCGDF-DGBAP No. 1052/2013; Brazil, No. IBAMA #144/2010 and SISBIO #40673; Cameroon, No. 000000010/MINRESI/B00/C00/C10/C12 and No. 00000075/MINRESI/B00/C00/C10/C12; Democratic Republic of Congo (DRC), collections permitted under Belgian Science Policy (BELSPO) institutional project COBIMFO (2010â2015), as well as during the Boyekoli-Ebale-Congo Expedition in 2010 (organised by: the University of Kisangani (DRC), the Royal Museum of Central Africa (Tervuren, Belgium), the Royal Belgian Institute of Natural Sciences, the National Botanical Garden of Belgium; supported by: the Belgian Development Cooperation, BELSPO); Ecuador, No. 06-2011-FAU-DPAP-MA; French Guiana, No. TREL1902817S/136; Israel, no permits required; Ivory Coast, permits granted for protected areas (Taï, Marahoué, Comoé, Banco) by Office Ivoirien des Parcs et Réserves; Kenya, No. NACOSTI/P/15/7983/7214; Madagascar, No. 035/19/MEDD/SG/DGF/DSAP/SCB.Re and No. 202-19/MEDD/SG/DIREDD.AATS.ANS.AND; Martinique, no permits required; Mexico, permits provided to the late Paul M. Ban in 1996; Panama, No. SEX/A-36-10; Paraguay, permits provided to Rudolf Scheffrahn through the Museo Nacional de Historia Natural del Paraguay in 2012; Peru, Carta No. 592-2014-MINAGRI-DGFFS/DGEFFS; Singapore, No. NP/RP12-063a; Thailand, sampling carried out through the Khao Chong Insect Laboratory; USA, no permits required; Venezuela, permits provided to Rudolf Scheffrahn in 2008.
Samples were either preserved in 80â100% ethanol and stored at room temperature or preserved in RNA-later®, initially stored for up to several weeks at room temperature, and subsequently at temperatures varying between ââ20â°C and ââ80â°C. Whole genomic DNA extraction was performed using the DNeasy Blood & Tissue extraction kit (Qiagen). Tissue samples were disrupted in 2âmL microtubes using two 3âmm steel beads and the TissueLyser II (Qiagen). Library preparation was achieved using the NEBNext® UltraTM II FS DNA Library Preparation Kit (New England Biolabs) and the Unique Dual Indexing Kit (New England Biolabs), with reagent volumes reduced to one-fifteenth of recommended volumes. Because the DNA of samples preserved in ~â80% ethanol at room temperature is typically highly fragmented, the incubation time of the enzymatic fragmentation step was set to a maximum of five minutes. This duration ensures a proper end repair, 5´ phosphorylation and dA-tailing, but limits any further significant fragmentation. Concentrations of final libraries were measured with a Qubit 3.0 Fluorometer (Invitrogen) using the QubitTM 1X dsDNA HS Assay Kit (Invitrogen). The fragment size distribution of libraries was assessed with a 4200 TapeStation System (Agilent) using the High Sensitivity D1000 Kit (Agilent). Libraries were pooled in equimolar concentration, and paired-end reads were generated using the NovaSeq 6000 and HiSeq X Illumina platforms at a read length of 150âbp.
Genome assemblies and extraction of UCEs
Raw reads were quality-trimmed using fastp v0.20.188 and assembled with metaSPAdes v3.1389. We extracted the core UCEs from these draft assemblies with flanking 200âbp at both 5â and 3â ends (~â600âbp loci) using PHYLUCE v1.6.690 and LASTZ91. We used the termite-specific bait set targeting 50,616 loci published by Hellemans et al.44. The species identification codes (âFile 1: ids_to_species.txtâ) and UCE dataset produced in this study (âFile 2: TER_UCE_DB_CONTRIB_4.fasta.gzâ) are available on the Dryad Digital Repository, at https://doi.org/10.5061/dryad.02v6wwqbm. This paper represents Contribution #4 to the Termite UCE Database, available at: https://github.com/oist/TER-UCE-DB/.
Alignments and supermatrices
Extracted UCEs were combined in a total of seven supermatrices (see below) used for subsequent phylogenetic analyses. All UCE loci were aligned with MAFFT92, as implemented in phyluce_align_seqcap_align. Alignments were either trimmed internally using phyluce_align_get_gblocks_trimmed_alignments_from_untrimmed that implements Gblocks93,94 with default parameters or only edge-trimmed with phyluce_align_seqcap_align. Loci absent in more than 30% of taxa were filtered out using phyluce_align_get_only_loci_with_min_taxa.
We assessed the compositional bias of each locus after both internal- and edge-trimming treatments (3234 and 3233 loci, respectively) to identify low-quality and rogue loci. We computed the relative composition frequency variability (RCFV) and its normalised version (nRCFV) using RCFV_Reader v195. Chi-square tests of nucleotide homogeneity were performed with BaCoCa v1.10596. None of the internally trimmed loci significantly deviated from homogeneity, while five edge-trimmed loci did (at the default αâ=â0.01; for details, see Supplementary Data 2, 3). However, the nRCFV values of these five loci were below 0.005, indicating the deviations from base composition homogeneity were marginal. Consequently, we kept all loci in our analyses.
From the internally trimmed alignments, we constructed three supermatrices using phyluce_align_format_nexus_files_for_raxml: the unfiltered, cogenic, and intergenic datasets. The cogenic and intergenic datasets were obtained using the functional annotation of UCEs from the GFF file (NCBI Annotation Release 100) accompanying the draft genome of the archotermopsid Zootermopsis nevadensis (GCF_000696155)97. The cogenic dataset consisted of merged non-overlapping UCEs (from the pool of the 40,966 singly annotated UCEs) occurring within the same gene.
From the edge-trimmed alignments, we constructed four supermatrices using the loci coding for proteins along their full length. The position of the UCE loci and their strand orientation was determined by mapping Z. nevadensis loci (TER-4-UCEDB) against their corresponding genome (accession GCF_000696155) and GFF annotation file (NCBI Annotation Release 100) using BLASTn v2.10.0+ and bedtools v2.29.298. The reading frame and corresponding amino acid sequences were identified and retrieved using the transeq function of EMBOSS v6.6.099. We constructed two datasets by filtering out sequences with stop codons: one retaining all loci and the other retaining only loci displaying fewer than three stop codons. For each dataset, we aligned the protein sequences of each locus with MAFFT v7.305 and back-translated the protein alignments into codon alignments with PAL2NAL v14100. Subsequent analyses were carried out on both the protein and the codon matrices for each dataset (i.e., four supermatrices). All concatenated sequence supermatrices are described in Supplementary Data 4 and provided in the Dryad Digital Repository (âFile 3: supermatrices.tarâ), at https://doi.org/10.5061/dryad.02v6wwqbm.
Phylogenetic reconstructions
We reconstructed 51 phylogenetic trees using diverse models and two main approaches: (1) maximum-likelihood phylogenetic analyses performed on concatenated sequence data and (2) multi-species coalescence analyses, which infer the species tree by summarising a collection of maximum-likelihood phylogenetic trees constructed for each UCE locus.
A total of 30 maximum-likelihood phylogenetic trees were reconstructed using IQ-TREE v1.6.12101 on concatenated sequence data with 1000 ultrafast bootstrap replicates (UFB) to assess branch supports102. A summary of the parameters used for all 30 analyses performed in this study is available in Supplementary Data 5. Phylogenetic trees were reconstructed using nucleotide, codon, or protein sequences. We either used the best-fit substitution model selected with the Bayesian Information Criterion using ModelFinder103 implemented in IQ-TREE or directly specified one of several predefined substitution models. Briefly, we used a GTRâ+âFâ+âIâ+âG4 nucleotide substitution model on five types of nucleotide alignments: unfiltered, cogenic, intergenic, and protein-coding with and without third codon positions. For alignments of protein-coding nucleotide sequences, we used the empirical ECMK07â+âFU model and the mechanistic GY2Kâ+âF and MG2Kâ+âF3X4 codon models104. For amino acid alignments, we used the improved general amino acid replacement matrix (LG) LGâ+âFâ+âG4 alone and combined with a C20 mixture of profiles (LGâ+âC20â+âFâ+âG)105,106, the revised Dayhoff matrix (DCMut)107, and the insect-based Q matrix (Q.insect)108,109. The models were either applied to the unpartitioned supermatrix or the loci-partitioned supermatrix by examining the top 10% of partitioning schemes using the fast relaxed clustering approach110,111. For each reconstruction, we used RogueNaRok v1.0.1-3112 to identify potential rogue taxa based on their bootstrap tree sets.
As maximum-likelihood phylogenetic trees reconstructed with concatenated data may fail to properly resolve short internodes113, we also ran a total of 21 multi-species coalescence reconstructions with ASTRAL-III v5.7.7114. ASTRAL was run on gene trees reconstructed independently for each locus with IQ-TREE using the same models and parameters as in the analyses on concatenated sequence data (for details, see Supplementary Data 5). For gene tree reconstructions, we used the option â-czbâ to allow polytomies and reduce gene tree biases. To control for poor-quality loci and rogue taxa, we ran all ASTRAL analyses once more after pruning outliers (at quantile 0.05) independently in each gene tree with TreeShrink v1.3.9115. Local posterior probabilities (LPPs) estimated from gene tree quartet frequencies116 were calculated as branch support values for each tree.
All phylogenetic trees reconstructed in this study are provided in the Dryad Digital Repository (âFile 4: trees.tarâ), at https://doi.org/10.5061/dryad.02v6wwqbm. We used phytools117 implemented in R v4.0.2118 to generate 50% majority rule consensus summary trees, independently for analyses performed on concatenated sequence data and for coalescence analyses (with or without pruning by TreeShrink). Branch support values are UFBs and LPPs for concatenation and coalescence analyses, respectively. We also calculated site concordance factors (sCFs) as an alternative branch support metric50. sCFs are defined as the percentage of decisive alignment sites supporting a branch. Branches with sCFs >â50% (i.e., supported by a majority of alignment sites) are considered robust. We calculated sCFs on each consensus tree with IQ-TREE using the option â-scfâ and the unfiltered supermatrix. Finally, we assessed alternative topologies for unresolved interfamilial nodes using the approximately unbiased (AU) test119. The tests were performed with IQ-TREE, using the unfiltered supermatrix as data input and 10,000 resampling estimated log-likelihood (RELL) replicates (option â-zb 10000â). Alternative topologies were tested against the majority rule consensus tree summarising the 30 analyses performed on concatenated sequence data.
Nomenclatural acts
This published work and the nomenclatural acts it contains have been registered in ZooBank, the proposed online registration system for the International Code of Zoological Nomenclature (ICZN). The ZooBank LSIDs (Life Science Identifiers) and the associated information can be accessed through a web browser by appending the LSID to the prefix âhttp://zoobank.org/â. The LSIDs for this publication are: urn:lsid:zoobank.org:act:5F7D1841-7543-4914-B552-10318E2224C2; urn:lsid:zoobank.org:act:DEEDCE51-6696-4C97-BB2E-654EF86C8F9F; urn:lsid:zoobank.org:act:429D8F1D-1D78-44AD-BF91-4AAE9492251F; urn:lsid:zoobank.org:act:B9567537-33EF-4B4D-992F-34E0C16B8701; urn:lsid:zoobank.org:act:1F08A9B3-6F94-498C-A8B4-78ABF508FDA6; urn:lsid:zoobank.org:act:2C4DDA82-D86C-40CA-BCB9-F5CD5AA60E7B.
Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
Data availability
Samples used in this study are described in Supplementary Data 1. The raw sequence reads generated in this study have been deposited in the NCBI Sequence Read Archive under the BioProject accession code PRJNA1122986. The raw data from prior publications used in this study are available under BioProject accession codes PRJNA560101, PRJNA1122307 and PRJNA1122966; and copys are available under accession codes GCF_000696155, GCF_002891405, GCA_013340265 and GigaDB100057. The raw sequence data can also be accessed using individual accession numbers given in Supplementary Data 1. The UCE data from prior publications used in this study are available in the Dryad Digital Repository at https://doi.org/10.5061/dryad.x0k6djhn0120 (file: âSupplementary_Data_6_uces_TER_UCE_DB_CONTRIB_1.fastaâ); https://doi.org/10.5061/dryad.tmpg4f53w121 (file: âTER_UCE_DB_CONTRIB_3.tar.gzâ). The source data for this study (UCE sequences, supermatrices, and trees) are available from the Dryad Digital Repository at https://doi.org/10.5061/dryad.02v6wwqbm122. This study is registered under the ZooBank LSID: https://zoobank.org/urn:lsid:zoobank.org:pub:B8A8464B-2FB2-4CA2-8330-A907111FD62E.
Code availability
No new custom code is published with this article. The Termite UCE database and associated codes are available at: https://github.com/sihellem/TER-UCE-DB/.
References
Lo, N. et al. Evidence from multiple gene sequences indicates that termites evolved from wood-feeding cockroaches. Curr. Biol. 10, 801â804 (2000).
Grimaldi, D. A. & Engel, M. S. Evolution of the Insects. (Cambridge University Press, 2005).
Inward, D. J. G., Vogler, A. P. & Eggleton, P. A comprehensive phylogenetic analysis of termites (Isoptera) illuminates key aspects of their evolutionary biology. Mol. Phylogenet. Evol. 44, 953â967 (2007).
Engel, M. S., Barden, P., Riccio, M. L. & Grimaldi, D. A. Morphologically specialized termite castes and advanced sociality in the early Cretaceous. Curr. Biol. 26, 522â530 (2016).
Bourguignon, T. et al. Transoceanic dispersal and plate tectonics shaped global cockroach distributions: evidence from mitochondrial phylogenomics. Mol. Biol. Evol. 35, 970â983 (2018).
Evangelista, D. A. et al. An integrative phylogenomic approach illuminates the evolutionary history of cockroaches and termites (Blattodea). Proc. R. Soc. B. 286, 20182076 (2019).
Krishna, K., Grimaldi, D. A., Krishna, V. & Engel, M. S. Treatise on the Isoptera of the World. 1. Introduction. Bull. Am. Mus. Nat. Hist. 377, 1â200 (2013).
Eggleton, P. et al. The diversity, abundance and biomass of termites under differing levels of disturbance in the Mbalmayo Forest Reserve, southern Cameroon. Philos. Trans. R. Soc. B. 351, 51â68 (1996).
Bignell, D. E. & Eggleton, P. Termites in ecosystems. in Termites: Evolution, Sociality, Symbioses, Ecology (eds. Abe, T., Bignell, D. E. & Higashi, M.) 363â387 (Kluwer Academic Publishers, 2000).
Holt, J. A. & Lepage, M. Termite and soil properties. in Termites: Evolution, Sociality, Symbioses, Ecology (eds. Abe, T., Bignell, D. E. & Higashi, M.) 389â407 (Kluwer Academic Publishers, 2000).
Evans, T. A., Dawes, T. Z., Ward, P. R. & Lo, N. Ants and termites increase crop yield in a dry climate. Nat. Commun. 2, 262 (2011).
Jouquet, P., Traoré, S., Choosai, C., Hartmann, C. & Bignell, D. Influence of termites on ecosystem functioning. Ecosystem services provided by termites. Eur. J. Soil Biol. 47, 215â222 (2011).
Bonachela, J. A. et al. Termite mounds can increase the robustness of dryland ecosystems to climatic change. Science 347, 651â655 (2015).
Ashton, L. A. et al. Termites mitigate the effects of drought in tropical rainforest. Science 363, 174â177 (2019).
Elizalde, L. et al. The ecosystem services provided by social insects: traits, management tools and knowledge gaps. Biol. Rev. 95, 1418â1441 (2020).
Evans, T. A., Forschler, B. T. & Kenneth Grace, J. Biology of invasive termites: a worldwide review. Annu. Rev. Entomol. 58, 455â474 (2013).
Gerozisis, J., Hadlington, P. & Staunton, I. Urban Pest Management in Australia. (University of New South Wales Press, 2008).
Dhang, P. Urban Pest Management: An Environmental Perspective. (CABI, 2011).
Rust, M. K. & Su, N. Y. Managing social insects of urban importance. Annu. Rev. Entomol. 57, 355â375 (2012).
Romero Arias, J. et al. Mitochondrial phylogenetics position a new Afrotropical termite species into its own subfamily, the Engelitermitinae (Blattodea: Termitidae). Syst. Entomol. 49, 72â83 (2024).
Froggatt, W. W. Australian Termitidae. Part II. Proc. Linn. Soc. N. South Wales 21, 510â552 (1897).
Desneux, J. à propos de la phylogénie des termitides. Ann. Soc. Entomol. Belgique 48, 278â286 (1904).
Holmgren, N. Das System der Termiten. Zool. Anz. 35, 284â286 (1910).
Holmgren, N. Termitenstudien. 2. Systematik der Termiten. Die Familien Mastotermitidae, Protermitidae und Mesotermitidae. K. Sven. Vetensk. Akad. Handl. 46, 1â86 (1911).
Holmgren, N. Termitenstudien. 3. Systematik der Termiten. Die Familie Metatermitidae. K. Sven. Vetensk. Akad. Handl. 48, 1â166 (1912).
Snyder, T. E. Catalog of the termites (Isoptera) of the world. Smithson. Misc. Collect. 112, 1â490 (1949).
Grassé, P.-P. Ordre des Isoptères ou Termites. in Traité de Zoologie IX (ed. Grassé, P.-P.) 408â544 (Masson, 1949).
Emerson, A. E. Geographical origins and dispersions of termite genera. Fieldiana Zool. 37, 465â521 (1955).
Krishna, K. Taxonomy, phylogeny, and distribution of termites. in Biology of Termites, Vol. 2 (eds. Krishna, K. & Weesner, F. M.) 127â152 (Academic Press, 1970).
Roonwal, M. L. & Chhotani, O. B. The Fauna of India and the Adjacent Countries. Isoptera (termites), Vol. 1. (Zoological Survey of India, 1989).
Engel, M. S., Grimaldi, D. A. & Krishna, K. Termites (Isoptera): their phylogeny, classification, and rise to ecological dominance. Am. Mus. Novit. 3650, 1â27 (2009).
Jiang, R.-X. et al. Further evidence of Cretaceous termitophily: description of new termite hosts of the trichopseniine Cretotrichopsenius (Coleoptera: Staphylinidae), with emendations to the classification of lower termites (Isoptera). Palaeoentomology 4, 374â389 (2021).
Wang, M. et al. Phylogeny, biogeography and classification of Teletisoptera (Blattaria: Isoptera). Syst. Entomol. 47, 581â590 (2022).
Legendre, F. et al. The phylogeny of termites (Dictyoptera: Isoptera) based on mitochondrial and nuclear markers: Implications for the evolution of the worker and pseudergate castes, and foraging behaviors. Mol. Phylogenet. Evol. 48, 615â627 (2008).
Ware, J. L., Grimaldi, D. A. & Engel, M. S. The effects of fossil placement and calibration on divergence times and rates: an example from the termites (Insecta: Isoptera). Arthropod Struct. Dev. 39, 204â219 (2010).
Cameron, S. L., Lo, N., Bourguignon, T., Svenson, G. J. & Evans, T. A. A mitochondrial genome phylogeny of termites (Blattodea: Termitoidae): robust support for interfamilial relationships and molecular synapomorphies define major clades. Mol. Phylogenet. Evol. 65, 163â173 (2012).
Bourguignon, T. et al. The evolutionary history of termites as inferred from 66 mitochondrial genomes. Mol. Biol. Evol. 32, 406â421 (2015).
Donovan, S. E., Jones, D. T., Sands, W. A. & Eggleton, P. Morphological phylogenetics of termites (Isoptera). Biol. J. Linn. Soc. 70, 467â513 (2000).
Jouault, C., Legendre, F., Grandcolas, P. & Nel, A. Revising dating estimates and the antiquity of eusociality in termites using the fossilized birthâdeath process. Syst. Entomol. 46, 592â610 (2021).
Lo, N., Kitade, O., Miura, T., Constantino, R. & Matsumoto, T. Molecular phylogeny of the Rhinotermitidae. Insectes Soc. 51, 365â371 (2004).
Legendre, F. et al. Phylogeny of Dictyoptera: dating the origin of cockroaches, praying mantises and termites with molecular data and controlled fossil evidence. PLoS ONE 10, e0130127 (2015).
Lo, N. & Eggleton, P. Termite phylogenetics and co-cladogenesis with symbionts. in Biology of Termites: A Modern Synthesis (eds. Bignell, D. E., Roisin, Y. & Lo, N.) 27â50 (Springer, 2011).
Bucek, A. et al. Evolution of termite symbiosis informed by transcriptome-based phylogenies. Curr. Biol. 29, 3728â3734 (2019).
Hellemans, S. et al. Using ultraconserved elements to reconstruct the termite tree of life. Mol. Phylogenet. Evol. 173, 107520 (2022).
Bourguignon, T. et al. Mitochondrial phylogenomics resolves the global spread of higher termites, ecosystem engineers of the tropics. Mol. Biol. Evol. 34, 589â597 (2017).
Bejerano, G. et al. Ultraconserved elements in the human genome. Science 304, 1321â1325 (2004).
Faircloth, B. C. et al. Ultraconserved elements anchor thousands of genetic markers spanning multiple evolutionary timescales. Syst. Biol. 61, 717â726 (2012).
Hedin, M., Derkarabetian, S., Alfaro, A., RamÃrez, M. J. & Bond, J. E. Phylogenomic analysis and revised classification of atypoid mygalomorph spiders (Araneae, Mygalomorphae), with notes on arachnid ultraconserved element loci. PeerJ 7, e6864 (2019).
Van Dam, M. H., Henderson, J. B., Esposito, L. & Trautwein, M. Genomic characterization and curation of UCEs improves species tree reconstruction. Syst. Biol. 70, 307â321 (2021).
Minh, B. Q., Hahn, M. W. & Lanfear, R. New methods to calculate concordance factors for phylogenomic datasets. Mol. Biol. Evol. 37, 2727â2733 (2020).
Wang, M. et al. Neoisoptera repeatedly colonised Madagascar after the Middle Miocene climatic optimum. Ecography 2023, e06463 (2023).
Rocha, M. M. Redescription of the enigmatic genus Genuotermes Emerson (Isoptera, Termitidae, Termitinae). Zookeys 340, 107â117 (2013).
Rocha, M. M., Morales-Corrêa e Castro, A. C., Cuezzo, C. & Cancello, E. M. Phylogenetic reconstruction of Syntermitinae (Isoptera, Termitidae) based on morphological and molecular data. PLoS ONE 12, e0174366 (2017).
Holt, B. G. et al. An update of Wallaceâs zoogeographic regions of the world. Science 339, 74â78 (2013).
Noirot, C. The gut of termites (Isoptera): comparative anatomy, systematics, phylogeny. II. - Higher termites (Termitidae). Ann. Soc. Entomol. France (N.S.)Â 37, 431â471 (2001).
Hagen, H. A. Monographie der Termiten. Linnaea Entomol. 12, 1â342 (1858).
Holmgren, N. Termitenstudien. 4. Versuch einer systematischen Monographie der Termiten der Orientalischen Region. K. Sven. Vetensk. Akad. Handl. 50, 1â276 (1913).
Desneux, J. Remarques critiques sur la phylogénie et la division systématique des Termitides (réponse à M. Wasmann). Ann. Soc. Entomol. Belgique 48, 372â378 (1904).
Desneux, J. Termites du Sahara algérien recueillis par M. le professeur Lameere. Ann. Soc. Entomol. Belgique 46, 436â440 (1902).
Quennedey, A. & Deligne, J. Lâarme frontale des soldats de termites. I. Rhinotermitidae. Insectes Soc. 22, 243â267 (1975).
Silvestri, F. Isoptera. in Die Fauna Südwest-Australiens. Ergebnisse der Hamburger Südwest-Australischen Forschungsreise 1905. Vol. 2 (eds. Michaelsen, W. & Hartmeyer, R.) 279â314 (Gustav Fischer, 1909).
Jankásek, M., Kotyková VaradÃnová, Z. & Stáhlavský, F. Blattodea karyotype database. Eur. J. Entomol. 118, 192â199 (2021).
Chouvenc, T., Å obotnÃk, J., Engel, M. S. & Bourguignon, T. Termite evolution: mutualistic associations, key innovations, and the rise of Termitidae. Cell. Mol. Life Sci. 78, 2749â2769 (2021).
Silvestri, F. Nota preliminare sui Termitidi sud-americani. Boll. Mus. Zool. Anat. Comp. R. Univ. Torino 16, 1â8 (1901).
Wasmann, E. Viaggio di Leonardo Fea in Birmania e regioni vicine LXXII. Neue Termitophilen und Termiten aus Indien. IâIII. Ann. del Mus. Civ. di Stor. Nat. di Genova (Ser. 2) 16, 613â630 (1896).
Sjöstedt, Y. Revision der Termiten Afrikas. 3. Monographie. K. Sven. Vetensk. Handl. (Ser. 3) 3, 1â419, 16pl. (1926).
Wasmann, E. Termiten, Termitophilen und Myrmekophilen, gesammelt auf Ceylon von Dr. W. Horn 1899, mit anderm ostindischen Material bearbeitet. Zool. Jahrb. Abt. Syst. Geog. Biol. Tiere 17, 99â164 (1902).
Constantino, R. Key to the soldiers of South American Heterotermes with a new species from Brazil (Isoptera: Rhinotermitidae). Insect Syst. Evol. 31, 463â472 (2000).
Kemner, N. A. Systematische und biologische Studien über die Termiten Javas und Celebesâ. K. Sven. Vetensk. Akad. Handl. 13, 1â241 (1934).
Miller, L. R. Invasitermes, a new genus of soldierless termites from Northern Australia (Isoptera: Termitidae). Aust. J. Entomol. 23, 33â37 (1984).
Emerson, A. E. The termites of Kartabo, Bartica District, British Guiana. Zoologica 6, 291â459 (1925).
Rocha, M. M. & Cuezzo, C. Redescription of the monotypic Neotropical genus Crepititermes Emerson (Termitidae: Termitinae). Neotrop. Entomol. 44, 457â465 (2015).
Holmgren, N. Studien über südamerikanische Termiten. Zool. Jahrb. Abt. Syst. Geog. Biol. Tiere 23, 521â676 (1906).
Emerson, A. E. Six new genera of Termitinae from the Belgian Congo (Isoptera, Termitidae). Am. Mus. Novit. 1988, 1â49 (1960).
Scheffrahn, R. H. & KÅeÄek, J. Redescription and reclassification of the African termite, Forficulitermes planifrons (Isoptera, Termitidae, Termitinae). Zootaxa 3946, 591â594 (2015).
Holmgren, N. The Percy Sladen Trust expedition to the Indian Ocean in 1905. Isoptera. Trans. Linn. Soc. Lond., Zool. 14, 135â148 (1910).
Holmgren, N. Wissenschaftliche Ergebnisse einer Forschungsreise nach Ostindien, ausgeführt im Auftrage der Kgl. Preuss. Akademie der Wissenschaften zu Berlin von H. v. Buttel-Reepen. III. Termiten aus Sumatra, Java, Malacca und Ceylon. Zool. Jahrb. Abt. Syst. Geog. Biol. Tiere 36, 229â290 (1914).
Roisin, Y. Schievitermes globicornis, a new genus and species of Termitinae (Blattodea, Termitidae) from French Guiana. ZooKeys 1125, 103â114 (2022).
Silvestri, F. Contribuzione alla conoscenza dei Termitidi e Termitofili dellâAfrica occidentale. I. Termitidi. Boll. Lab. Zool. Gen. Agrar. R. Sc. Super. Agricoltura Portici 9, 1â146 (1914).
Latreille, P. A. Histoire Naturelle, Générale et Particulière des Crustacés et des Insectes. Vol. 3. (F. Dufart, 1802).
Linnaeus, C. Systema Naturae per Regna Tria Natura, Secundum Classes, Ordines, Genera, Species, Cum Characteribus, Differentiis, Synonymis, Locis [10th ed. (revised), Vol. 1]. (Laurentius Salvius, 1758).
Weidner, H. Beiträge zur Kenntnis der Termiten Angolas, hauptsächlich auf Grund der Sammlungen und Beobachtungen von A. de Barros Machado (I. Beitrag). Publicações Cult. Cia. Diam. Angola 29, 55â106 (1956).
Wasmann, E. Termiten von Madagaskar und Ostafrika. Abh. Senckenbergi. Naturf. Ges. 21, 137â182 (1897).
Hare, L. Termite phylogeny as evidenced by soldier mandible development. Ann. Entomol. Soc. Am. 37, 459â486 (1937).
Whitfield, J. B. & Lockhart, P. J. Deciphering ancient rapid radiations. Trends Ecol. Evol. 22, 258â265 (2007).
Degnan, J. H. & Rosenberg, N. A. Gene tree discordance, phylogenetic inference and the multispecies coalescent. Trends Ecol. Evol. 24, 332â340 (2009).
Arora, J. et al. Evidence of cospeciation between termites and their gut bacteria on a geological time scale. Proc. R. Soc. B. 290, 20230619 (2023).
Chen, S., Zhou, Y., Chen, Y. & Gu, J. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884âi890 (2018).
Nurk, S., Meleshko, D., Korobeynikov, A. & Pevzner, P. A. metaSPAdes: a new versatile metagenomic assembler. Genome Res. 27, 824â834 (2017).
Faircloth, B. C. PHYLUCE is a software package for the analysis of conserved genomic loci. Bioinformatics 32, 786â788 (2016).
Harris, R. S. Improved Pairwise Alignment of Genomic DNA. (The Pennsylvania State University, 2007).
Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772â780 (2013).
Castresana, J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol. Biol. Evol. 17, 540â552 (2000).
Talavera, G. & Castresana, J. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst. Biol. 56, 564â577 (2007).
Fleming, J. F. & Struck, T. H. nRCFV: a new, dataset-size-independent metric to quantify compositional heterogeneity in nucleotide and amino acid datasets. BMC Bioinformatics 24, 145 (2023).
Kück, P. & Struck, T. H. BaCoCa - A heuristic software tool for the parallel assessment of sequence biases in hundreds of gene and taxon partitions. Mol. Phylogenet. Evol. 70, 94â98 (2014).
Terrapon, N. et al. Molecular traces of alternative social organization in a termite genome. Nat. Commun. 5, 3636 (2014).
Quinlan, A. R. & Hall, I. M. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26, 841â842 (2010).
Rice, P., Longden, L. & Bleasby, A. EMBOSS: the European molecular biology open software suite. Trends Genet. 16, 276â277 (2000).
Suyama, M., Torrents, D. & Bork, P. PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res. 34, W609âW612 (2006).
Nguyen, L. T., Schmidt, H. A., Von Haeseler, A. & Minh, B. Q. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 32, 268â274 (2015).
Hoang, D. T., Chernomor, O., von Haeseler, A., Minh, B. Q. & Vinh, L. S. UFBoot2: improving the ultrafast bootstrap approximation. Mol. Biol. Evol. 35, 518â522 (2018).
Kalyaanamoorthy, S., Minh, B. Q., Wong, T. K. F., Von Haeseler, A. & Jermiin, L. S. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat. Methods 14, 587â589 (2017).
Kosiol, C., Holmes, I. & Goldman, N. An empirical codon model for protein sequence evolution. Mol. Biol. Evol. 24, 1464â1479 (2007).
Le, S. Q. & Gascuel, O. An improved general amino acid replacement matrix. Mol. Biol. Evol. 25, 1307â1320 (2008).
Le, S. Q., Gascuel, O. & Lartillot, N. Empirical profile mixture models for phylogenetic reconstruction. Bioinformatics 24, 2317â2323 (2008).
Kosiol, C. & Goldman, N. Different versions of the Dayhoff rate matrix. Mol. Biol. Evol. 22, 193â199 (2005).
Minh, B. Q., Dang, C. C., Vinh, L. S. & Lanfear, R. QMaker: fast and accurate method to estimate empirical models of protein evolution. Syst. Biol. 70, 1046â1060 (2021).
Misof, B. et al. Phylogenomics resolves the timing and pattern of insect evolution. Science 346, 763â767 (2014).
Chernomor, O., Von Haeseler, A. & Minh, B. Q. Terrace aware data structure for phylogenomic inference from supermatrices. Syst. Biol. 65, 997â1008 (2016).
Lanfear, R., Frandsen, P. B., Wright, A. M., Senfeld, T. & Calcott, B. Partitionfinder 2: new methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses. Mol. Biol. Evol. 34, 772â773 (2017).
Aberer, A. J., Krompass, D. & Stamatakis, A. Pruning rogue taxa improves phylogenetic accuracy: an efficient algorithm and webservice. Syst. Biol. 62, 162â166 (2013).
Mendes, F. K. & Hahn, M. W. Why concatenation fails near the anomaly zone. Syst. Biol. 67, 158â169 (2018).
Zhang, C., Rabiee, M., Sayyari, E. & Mirarab, S. ASTRAL-III: polynomial time species tree reconstruction from partially resolved gene trees. BMC Bioinformatics 19, 153 (2018).
Mai, U. & Mirarab, S. TreeShrink: fast and accurate detection of outlier long branches in collections of phylogenetic trees. BMC Genomics 19, 272 (2018).
Sayyari, E. & Mirarab, S. Fast coalescent-based computation of local branch support from quartet frequencies. Mol. Biol. Evol. 33, 1654â1668 (2016).
Revell, L. J. phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol. 3, 217â223 (2012).
R Core Team. R: a Language and Environment for Statistical Computing. (2020).
Shimodaira, H. An approximately unbiased test of phylogenetic tree selection. Syst. Biol. 51, 492â508 (2002).
Hellemans, S. et al. Data from: Using ultraconserved elements to reconstruct the termite tree of life. Dryad, Dataset https://doi.org/10.5061/dryad.x0k6djhn0 (2022).
Arora, J. et al. Evidence of cospeciation between termites and their gut bacteria on a geological time scale. Dryad, Dataset https://doi.org/10.5061/dryad.tmpg4f53w (2023).
Hellemans, S., Wang, M., Kaymak, E. & Bourguignon, T. Data from: Genomic data provide insights into the classification of extant termites. Dryad, Dataset https://doi.org/10.5061/dryad.02v6wwqbm (2024).
Hagen, H. A. Hr. Peters Berichtete über die von ihm gesammelten und von Hrn. Dr. Hermann Hagen bearbeiteten Neuropteren aus Mossambique. Ber. Akad. Wiss. Berlin 18, 479â484 (1853).
Desneux, J. Isoptera Fam. Termitidæ. in Genera Insectorum (ed. Wytsman, P.) 25, 1â52+2pls (P. Wytsman, 1904).
Emerson, A. E. The relations of a relict South African Termite (Isoptera, Hodotermitidae, Stolotermes). Am. Museum Novit. 1â12 (1942).
Holmgren, K. & Holmgren, N. Report on a collection of termites from India. Mem. Dep. Agric. India 5, 135â171 (1917).
Engel, M. S. & Krishna, K. Family-group names for termites (Isoptera). Am. Mus. Novit. 3432, 1â9 (2004).
Holmgren, N. Termitenstudien. 1. Anatomische Untersuchungen. K. Sven. Vetensk. Akad. Handl. 44, 1â215 (1909).
Grassé, P.-P. & Noirot, C. Apicotermes arquieri (Isoptère): ses constructions, sa biologie. Considérations générales sur la sous-famille des Apicotermitinae nov. Ann. des Sci. Nat., Zool. (11e Série) 16, 345â388 (1955).
Dudley, P. H. Termites of the isthmus of Panama.âPart II. Trans. N. Y. Acad. Sci. 9, 157â180 (1890).
Wasmann, E. Beispiele rezenter Artenbildung bei Ameisengästen und Termitengästen. Biol. Cent. 26, 565â580 (1906).
Acknowledgements
We thank our esteemed colleagues P. Ban, J. Cabral, J. Chase, B. Forschler, M. Hasson, R. Hickman, N. Kanzaki, P. Kasangij, J. KÅeÄek, E. Kuswanto, B. Leru, J. Mangold, A. Mullins, T. Nishimura, R.G. Santos, R. Setter, and S. Traoré, for the collection of samples. We thank the Sequencing Section (SQC) and the Scientific Computing & Data Analysis Section (SCDA) of OIST for assistance with sequencing and providing access to the OIST computing cluster, respectively. This work was supported by the Czech Science Foundation (project No. 20-20548S to T.B.), as well as subsidiary funding from OIST to T.B., including funding for a workshop held at OIST in early December 2022 on the phylogeny and classification of termites. S.H. acknowledges the Fonds de la Recherche Scientifique â FNRS. The participation of M.S.E. was funded by CONCYTEC through the PROCIENCIA programme within the framework of the call âInterinstitutional Alliances for Doctorate Programmesâ, according to contract PE501084299-2023-PROCIENCIA-BM.
Author information
Authors and Affiliations
Contributions
S.H. and T.B. conceptualised the experiments. S.H. and T.B. wrote the original draft manuscript, with significant input from M.M.R., E.C., and M.S.E. S.H., G.J., Y.R., E.C., M.R., J.Å ., R.H.S., and T.B. collected the samples. M.W. and E.K. performed lab experiments and generated data. S.H. analysed the data. All authors (S.H., M.M.R., M.W., J.R.A., D.K.A., A.G.B., A.B., T.F.C., T.C., C.C., J.P.C., R.C., F.D., J.D., P.E., T.A.E., R.H., M.C.H., M.H., G.J., C.J., C.M.K., E.K., J.K., C.Y.L., F.L., H.F.L., N.L., T.L., K.e.M., K.i.M., D.M., N.M., D.E.O., M.P., D.S.D., N.Y.S., G.T., E.V., J.L.W., J.S., R.H.S., E.C., Y.R., M.S.E., and T.B.) edited and accepted the final version of this manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks Diego Souza and the other anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the articleâs Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
About this article
Cite this article
Hellemans, S., Rocha, M.M., Wang, M. et al. Genomic data provide insights into the classification of extant termites. Nat Commun 15, 6724 (2024). https://doi.org/10.1038/s41467-024-51028-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-024-51028-y