Efficient but occasionally imperfect vertical transmission of gut mutualistic protists in a wood‐feeding termite

Termites and subsocial roaches inherited many bacterial-borne carbohydrate-active enzymes (CAZymes) from their common ancestor

Termites digest wood using Carbohydrate-Active Enzymes (CAZymes) produced by gut bacteria with whom they have cospeciated at geological timescales. Whether CAZymes were encoded in the genomes of their ancestor's gut bacteria and transmitted to modern termites or acquired more recently from bacteria not associated with termites is unclear. We used gut metagenomes from 195 termites and one Cryptocercus, the sister group of termites, to investigate the evolution of termite gut bacterial CAZymes. We found 420 termite-specific clusters in 81 bacterial CAZyme gene trees, including 404 clusters showing strong cophylogenetic patterns with termites. Of the 420 clusters, 131 included at least one bacterial CAZyme sequence associated with Cryptocercus or Mastotermes, the sister group of all other termites. Our results suggest many bacterial CAZymes have been encoded in the genomes of termite gut bacteria since termite origin, indicating termites rely upon many bacterial CAZymes endemic to their guts to digest wood.

Updated classification of the phylum Parabasalia

The phylum Parabasalia includes very diverse single-cell organisms that nevertheless share a distinctive set of morphological traits. Most are harmless or beneficial gut symbionts of animals, but some have turned into parasites in other body compartments, the most notorious example being Trichomonas vaginalis in humans. Parabasalians have garnered attention for their nutritional symbioses with termites, their modified anaerobic mitochondria (hydrogenosomes), their character evolution, and the wholly unique features of some species. The molecular revolution confirmed the monophyly of Parabasalia, but considerably changed our view of their internal relationships, prompting a comprehensive reclassification 14 years ago. This classification has remained authoritative for many subgroups despite a greatly expanded pool of available data, but the large number of species and sequences that have since come out allow for taxonomic refinements in certain lineages, which we undertake here. We aimed to introduce as little disruption as possible but at the same time ensure that most taxa are truly monophyletic, and that the larger clades are subdivided into meaningful units. In doing so, we also highlighted correlations between the phylogeny of parabasalians and that of their hosts.

Transmission dynamics of symbiotic protist communities in the termite gut: association with host adult eclosion and dispersal

The fidelity of vertical transmission is a critical factor in maintaining mutualistic associations with microorganisms. The obligate mutualism between termites and intestinal protist communities has been maintained for over 130 million years, suggesting the faithful transmission of diverse protist species across host generations. Although a severe bottleneck can occur when alates disperse with gut protists, how protist communities are maintained during this process remains largely unknown. In this study, we examined the dynamics of intestinal protist communities during adult eclosion and alate dispersal in the termite Reticulitermes speratus. We found that the protist community structure in last-instar nymphs differed significantly from that in workers and persisted intact during adult eclosion, whereas all protists disappeared from the gut during moults between worker stages. The number of protists in nymphs and alates was substantially lower than in workers, whereas the proportion of protist species exhibiting low abundance in workers was higher in nymphs and alates. Using a simulation-based approach, we demonstrate that such changes in the protist community composition of nymphs and alates improve the transmission efficiency of whole protist species communities. This study thus provides novel insights into how termites have maintained mutualistic relationships with diverse gut microbiota for generations.

Protist symbionts of termites: diversity, distribution, and coevolution

The symbiosis between termites and their hindgut protists is mutually obligate and vertically inherited. It was established by the late Jurassic in the cockroach ancestors of termites as they transitioned to wood feeding. Since then, protist symbionts have been transmitted from host generation to host generation by proctodeal trophallaxis (anal feeding). The protists belong to multiple lineages within the eukaryotic superphylum Metamonada. Most of these lineages have evolved large cells with complex morphology, unlike the non-termite-associated Metamonada. The species richness and taxonomic composition of symbiotic protist communities varies widely across termite lineages, especially within the deep-branching clade Teletisoptera. In general, closely related termites tend to harbour closely related protists, and deep-branching termites tend to harbour deep-branching protists, reflecting their broad-scale co-diversification. A closer view, however, reveals a complex distribution of protist lineages across hosts. Some protist taxa are common, some are rare, some are widespread, and some are restricted to a single host family or genus. Some protist taxa can be found in only a few, distantly related, host species. Thus, the long history of co-diversification in this symbiosis has been complicated by lineage-specific loss of symbionts, transfer of symbionts from one host lineage to another, and by independent diversification of the symbionts relative to their hosts. This review aims to introduce the biology of this important symbiosis and serve as a gateway to the diversity and systematics literature for both termites and protists. A searchable database with all termite-protist occurrence records and taxonomic references is provided as a supplementary file to encourage and facilitate new research in this field.

Cellulolytic and Ethanologenic Evaluation of Heterotermes indicola's Gut-Associated Bacterial Isolates

Cellulose is the basic component of lignocellulosic biomass (LCB) making it a suitable substrate for bioethanol fermentation. Cellulolytic and ethanologenic bacteria possess cellulases that convert cellulose to glucose, which in turn yields ethanol subsequently. Heterotermes indicola is a subterranean termite that causes destructive damage by consuming wooden structures of infrastructure, LCB products, etc. Prospectively, the study envisioned the screening of cellulolytic and ethanologenic bacteria from the termite gut. Twenty six bacterial strains (H1-H26) based on varied colonial morphologies were isolated. Bacterial cellulolytic activity was tested biochemically. Marked gas production in the form of bubbles (0.1-4 cm) in Durham tubes was observed in H3, H7, H13, H15, H17, H21, and H22. Sugar degradation of all isolates was indicated by pink to maroon color development with the tetrazolium salt. Hallow zones (0.42-11 mm) by Congo red staining was exhibited by all strains except H2, H7, H8, and H19. Among the 26 bacterial isolates, 12 strains were identified as efficient cellulolytic bacteria. CMCase activity and ethanol titer of all isolates varied from 1.30 ± 0.03 (H13) to 1.83 ± 0.01 (H21) umol/mL/min and 2.36 ± 0.01 (H25) to 7.00 ± 0.01 (H21) g/L, respectively. Likewise, isolate H21 exhibited an ethanol yield of 0.40 ± 0.10 g/g with 78.38 ± 2.05% fermentation efficiency. Molecular characterization of four strains, Staphylococcus sp. H13, Acinetobacter baumanni H17, Acinetobacter sp. H21, and Acinetobacter nosocomialis H22, were based on the maximum cellulolytic index and the ethanol yield. H. indicola harbor promising and novel bacteria with a natural cellulolytic tendency for efficient bioconversion of LCB to value-added products. Hence, the selected cellulolytic bacteria can become an excellent addition for use in enzyme purification, composting, and production of biofuel at large.

Weight and protozoa number but not bacteria diversity are associated with successful pair formation of dealates in the Formosan subterranean termite, Coptotermes formosanus

New colonies of Formosan subterranean termites are founded by monogamous pairs. During swarming season, alates (winged reproductives) leave their parental colony. After swarming, they drop to the ground, shed their wings, and male and female dealates find suitable nesting sites where they mate and become kings and queens of new colonies. The first generation of offspring is entirely dependent on the nutritional resources of the founder pair consisting of the fat and protein reserves of the dealates and their microbiota, which include the cellulose-digesting protozoa and diverse bacteria. Since termite kings and queens can live for decades, mate for life and colony success is linked to those initial resources, we hypothesized that gut microbiota of founders affect pair formation. To test this hypothesis, we collected pairs found in nest chambers and single male and female dealates from four swarm populations. The association of three factors (pairing status, sex of the dealates and population) with dealate weights, total protozoa, and protozoa Pseudotrichonympha grassii numbers in dealate hindguts was determined. In addition, Illumina 16S rRNA gene sequencing and the QIIME2 pipeline were used to determine the impact of those three factors on gut bacteria diversity of dealates. Here we report that pairing status was significantly affected by weight and total protozoa numbers, but not by P. grassii numbers and bacteria diversity. Weight and total protozoa numbers were higher in paired compared to single dealates. Males contained significantly higher P. grassii numbers and bacteria richness and marginally higher phylogenetic diversity despite having lower weights than females. In conclusion, this study showed that dealates with high body weight and protozoa numbers are more likely to pair and become colony founders, probably because of competitive advantage. The combined nutritional resources provided by body weight and protozoa symbionts of the parents are important for successful colony foundation and development.

Selective enrichment of founding reproductive microbiomes allows extensive vertical transmission in a fungus-farming termite

Mutualistic coevolution can be mediated by vertical transmission of symbionts between host generations. Termites host complex gut bacterial communities with evolutionary histories indicative of mixed-mode transmission. Here, we document that vertical transmission of gut bacterial strains is congruent across parent to offspring colonies in four pedigrees of the fungus-farming termite Macrotermes natalensis. We show that 44% of the offspring colony microbiome, including more than 80 bacterial genera and pedigree-specific strains, are consistently inherited. We go on to demonstrate that this is achieved because colony-founding reproductives are selectively enriched with a set of non-random, environmentally sensitive and termite-specific gut microbes from their colonies of origin. These symbionts transfer to offspring colony workers with high fidelity, after which priority effects appear to influence the composition of the establishing microbiome. Termite reproductives thus secure transmission of complex communities of specific, co-evolved microbes that are critical to their offspring colonies. Extensive yet imperfect inheritance implies that the maturing colony benefits from acquiring environmental microbes to complement combinations of termite, fungus and vertically transmitted microbes; a mode of transmission that is emerging as a prevailing strategy for hosts to assemble complex adaptive microbiomes.

Understanding the Role of Free-Living Bacteria in the Gut of the Lower Termite Coptotermes gestroi Based on Metagenomic DNA Analysis

Termites' digestive systems, particularly in lower termites with the presence of protozoa, are unique ecological niches that shelter a diverse microbiota with a variety of functions for the host and the environment. In 2012, the metagenomic DNA (5.4 Gb) of the prokaryotes that freely live in the gut of the lower termite Coptotermes gestroi were sequenced. A total of 125,431 genes were predicted and analyzed in order to mine lignocellulolytic genes. however, the overall picture of the structure, diversity, and function of the prokaryotic gut microbiota was not investigated. In the present study, these 125,431 genes were taxonomically classified by MEGAN and functionally annotated by the Kyoto Encyclopedia of Genes and Genomes (KEGG) and by the Carbohydrate-Active enZYmes (CAZy) and HMMER databases. As a result, 95,751 bacterial genes were classified into 35 phyla. The structure of the bacteria, typified by a high ratio of Firmicutes to Bacterioidetes, was distinct from the structure of the entirety of the bacteria in the lower or higher termites' guts. The archaea (533 genes) were distributed into 4 phyla, 10 classes, 15 orders, 21 families, 47 genera, and 61 species. Although freely living in the guts, the prokaryotic community was formed, developed, and adapted to exhibit unique interactions in order to perform mutual roles of benefit to their hosts. Methanobacteriales, accounting for 61% of the archaea symbionts, seem to play an important role in methanogenesis. Concomitantly, bacterial methanotrophs in the gut utilize methane and combine with other bacterial groups, including potential lignocellulolytic degraders, acetogens, sulfur bacteria, and nitrogen-recycling bacteria, to efficiently convert wood with little nitrogen into acetates via certain pathway modules specified by prokaryotes that freely live in the gut. This forms an important energy source for the termites. Furthermore, bacteria carry 2223 genes involved in the biosynthesis of 17 antibiotic groups. The gut bacteria also possess genes for the degradation of 18 toxic aromatic compounds, of which four are commercial pesticides against termites commonly used for the preservation of wooden constructions. Eight of the eighteen pathways were the first to be reported from the termite gut. Overall, this study sheds light on the roles of the freely living bacteria and archaea in the C. gestroi gut, providing evidence that the gut microbiome acts as the second host genome, contributing both nutrients and immunity to support the host's existence, growth, and development.

New Parabasalia symbionts Snyderella spp. and Daimonympha gen. nov. from South American Rugitermes termites and the parallel evolution of a cell with a rotating "head"

Most Parabasalia are symbionts in the hindgut of "lower" (non-Termitidae) termites, where they widely vary in morphology and degree of morphological complexity. Large and complex cells in the class Cristamonadea evolved by replicating a fundamental unit, the karyomastigont, in various ways. We describe here four new species of Calonymphidae (Cristamonadea) from Rugitermes hosts, assigned to the genus Snyderella based on diagnostic features (including the karyomastigont pattern) and molecular phylogeny. We also report a new genus of Calonymphidae, Daimonympha, from Rugitermes laticollis. Daimonympha's morphology does not match that of any known Parabasalia, and its SSU rRNA gene sequence corroborates this distinction. Daimonympha does however share a puzzling feature with a few previously described, but distantly related, Cristamonadea: a rapid, smooth, and continuous rotation of the anterior end of the cell, including the many karyomastigont nuclei. The function of this rotatory movement, the cellular mechanisms enabling it, and the way the cell deals with the consequent cell membrane shear, are all unknown. "Rotating wheel" structures are famously rare in biology, with prokaryotic flagella being the main exception; these mysterious spinning cells found only among Parabasalia are another, far less understood, example.

Molecular phylogeny of Spirotrichonymphea (Parabasalia) with emphasis on Spironympha, Spirotrichonympha, and three new genera Pseudospironympha, Nanospironympha, and Brugerollina

Spirotrichonymphea, one of the six classes of phylum Parabasalia, are characterized by bearing many flagella in spiral rows, and they occur exclusively in the guts of termites. Phylogenetic relationships among the 13 described genera are not well understood due to complex morphological evolution and a paucity of molecular data. One such understudied genus is Spironympha. It has been variously considered a valid genus, a subgenus of Spirotrichonympha, or an "immature" life cycle stage of Spirotrichonympha. To clarify this, we sequenced the small subunit rRNA gene sequences of Spironympha and Spirotrichonympha cells isolated from the hindguts of Reticulitermes species and Hodotermopsis sjostedti and confirmed the molecular identity of H. sjostedti symbionts using fluorescence in situ hybridization. Spironympha as currently circumscribed is polyphyletic, with both H. sjostedti symbiont species branching separately from the "true" Spironympha from Reticulitermes. Similarly, the Spirotrichonympha symbiont of H. sjostedti branches separately from the "true" Spirotrichonympha found in Reticulitermes. Our data support Spironympha from Reticulitermes as a valid genus most closely related to Spirotrichonympha, though its monophyly and interspecific relationships are not resolved in our molecular phylogenetic analysis. We propose three new genera to accommodate the H. sjostedti symbionts and two new species of Spirotrichonympha from Reticulitermes.

Updated Virophage Taxonomy and Distinction from Polinton-like Viruses

Virophages are small dsDNA viruses that hijack the machinery of giant viruses during the co-infection of a protist (i.e., microeukaryotic) host and represent an exceptional case of "hyperparasitism" in the viral world. While only a handful of virophages have been isolated, a vast diversity of virophage-like sequences have been uncovered from diverse metagenomes. Their wide ecological distribution, idiosyncratic infection and replication strategy, ability to integrate into protist and giant virus genomes and potential role in antiviral defense have made virophages a topic of broad interest. However, one limitation for further studies is the lack of clarity regarding the nomenclature and taxonomy of this group of viruses. Specifically, virophages have been linked in the literature to other "virophage-like" mobile genetic elements and viruses, including polinton-like viruses (PLVs), but there are no formal demarcation criteria and proper nomenclature for either group, i.e., virophage or PLVs. Here, as part of the ICTV Virophage Study Group, we leverage a large set of genomes gathered from published datasets as well as newly generated protist genomes to propose delineation criteria and classification methods at multiple taxonomic ranks for virophages 'sensu stricto', i.e., genomes related to the prototype isolates Sputnik and mavirus. Based on a combination of comparative genomics and phylogenetic analyses, we show that this group of virophages forms a cohesive taxon that we propose to establish at the class level and suggest a subdivision into four orders and seven families with distinctive ecogenomic features. Finally, to illustrate how the proposed delineation criteria and classification method would be used, we apply these to two recently published datasets, which we show include both virophages and other virophage-related elements. Overall, we see this proposed classification as a necessary first step to provide a robust taxonomic framework in this area of the virosphere, which will need to be expanded in the future to cover other virophage-related viruses such as PLVs.

Succession of the microbiota in the gut of reproductives of Macrotermes subhyalinus (Termitidae) at colony foundation gives insights into symbionts transmission

Termites have co-evolved with a complex gut microbiota consisting mostly of exclusive resident taxa, but key forces sustaining this exclusive partnership are still poorly understood. The potential for primary reproductives to vertically transmit their gut microbiota (mycobiome and bacteriome) to offspring was investigated using colony foundations from field-derived swarming alates of Macrotermes subhyalinus. Metabarcoding based on the fungal internal transcribed spacer (ITS) region and the bacterial 16S rRNA gene was used to characterize the reproductives mycobiome and bacteriome over the colony foundation time. The mycobiome of swarming alates differed from that of workers of Macrotermitinae and changed randomly within and between sampling time points, highlighting no close link with the gut habitat. The fungal ectosymbiont Termitomyces was lost early from the gut of reproductives, confirming the absence of vertical transmission to offspring. Unlike fungi, the bacteriome of alates mirrored that of workers of Macroterminae. Key genera and core OTUs inherited from the mother colony mostly persisted in the gut of reproductive until the emergence of workers, enabling their vertical transmission and explaining why they were found in offspring workers. These findings demonstrate that the parental transmission may greatly contribute to the maintenance of the bacteriome and its co-evolution with termite hosts at short time scales.

The functional evolution of termite gut microbiota

BACKGROUND

Termites primarily feed on lignocellulose or soil in association with specific gut microbes. The functioning of the termite gut microbiota is partly understood in a handful of wood-feeding pest species but remains largely unknown in other taxa. We intend to fill this gap and provide a global understanding of the functional evolution of termite gut microbiota.

RESULTS

We sequenced the gut metagenomes of 145 samples representative of the termite diversity. We show that the prokaryotic fraction of the gut microbiota of all termites possesses similar genes for carbohydrate and nitrogen metabolisms, in proportions varying with termite phylogenetic position and diet. The presence of a conserved set of gut prokaryotic genes implies that essential nutritional functions were present in the ancestor of modern termites. Furthermore, the abundance of these genes largely correlated with the host phylogeny. Finally, we found that the adaptation to a diet of soil by some termite lineages was accompanied by a change in the stoichiometry of genes involved in important nutritional functions rather than by the acquisition of new genes and pathways.

CONCLUSIONS

Our results reveal that the composition and function of termite gut prokaryotic communities have been remarkably conserved since termites first appeared ~ 150 million years ago. Therefore, the "world's smallest bioreactor" has been operating as a multipartite symbiosis composed of termites, archaea, bacteria, and cellulolytic flagellates since its inception. Video Abstract.

Fine-scale genetic diversity and putative ecotypes of oxymonad protists coinhabiting the hindgut of Reticulitermes speratus

The hindgut of lower termites is generally coinhabited by multiple morphologically identifiable protist species. However, it is unclear how many protist species truly coexist in this miniaturized environment, and moreover, it is difficult to define the fundamental unit of protist diversity. Species delineation of termite gut protists has therefore been guided without a theory-based concept of species. Here, we focused on the hindgut of the termite Reticulitermes speratus, where 10 or 11 morphologically distinct oxymonad cell types, that is, morphospecies, coexist. We elucidated the phylogenetic structure of all co-occurring oxymonads and addressed whether their diversity can be explained by the "ecotype" hypothesis. Oxymonad-specific 18S rRNA gene amplicon sequencing analyses of whole-gut samples, combined with single-cell 18S rRNA sequencing of the oxymonad morphospecies, identified 210 one-nucleotide-level variants. The phylogenetic analysis of these variants revealed the presence of microdiverse clusters typically within 1% sequence divergence. Each known oxymonad morphospecies comprised one to several monophyletic or paraphyletic microdiverse clusters. Using these sequence data sets, we conducted computational simulation to predict the rates of ecotype formation and periodic selection, and to demarcate putative ecotypes. Our simulations suggested that the oxymonad genetic divergence is constrained primarily by strong selection, in spite of limited population size and possible bottlenecks during intergenerational transmission. A total of 33 oxymonad ecotypes were predicted, and most of the putative ecotypes were consistently detected among different colonies and host individuals. These findings provide a possible theoretical basis for species diversity and underlying mechanisms of coexistence of termite gut protists.

pr2-primers: An 18S rRNA primer database for protists

Metabarcoding of microbial eukaryotes (collectively known as protists) has developed tremendously in the last decade, almost solely relying on the 18S rRNA gene. As microbial eukaryotes are extremely diverse, many primers and primer pairs have been developed. To cover a relevant and representative fraction of the protist community in a given study system, an informed primer choice is necessary, as no primer pair can target all protists equally well. As such, a smart primer choice is very difficult even for experts and there are very few online resources available to list existing primers. We built a database listing 285 primers and 83 unique primer pairs that have been used for eukaryotic 18S rRNA gene metabarcoding. In silico performance of primer pairs was tested against two sequence databases: PR² version 4.12.0 for eukaryotes and a subset of silva version 132 for bacteria and archaea. We developed an R-based web application enabling browsing of the database, visualization of the taxonomic distribution of the amplified sequences with the number of mismatches, and testing any user-defined primer or primer set (https://app.pr2-primers.org). Taxonomic specificity of primer pairs, amplicon size and location of mismatches can also be determined. We identified universal primer sets that matched the largest number of sequences and analysed the specificity of some primer sets designed to target certain groups. This tool enables guided primer choices that will help a wide range of researchers to include protists as part of their investigations.

Alteration of gut microbiota with a broad-spectrum antibiotic does not impair maternal care in the European earwig

The microbes residing within the gut of an animal host often increase their own fitness by modifying their host's physiological, reproductive and behavioural functions. Whereas recent studies suggest that they may also shape host sociality and therefore have critical effects on animal social evolution, the impact of the gut microbiota on maternal care remains unexplored. This is surprising, as this behaviour is widespread among animals, often determines the fitness of both juveniles and parents, and is essential in the evolution of complex animal societies. Here, we tested whether life-long alterations of the gut microbiota with rifampicin-a broad-spectrum antibiotic-impair pre- and post-hatching maternal care in the European earwig. Our results first confirm that rifampicin altered the mothers' gut microbial communities and indicate that the composition of the gut microbiota differs before and after egg care. Contrary to our predictions, however, the rifampicin-induced alterations of the gut microbiota did not modify pre- or post-hatching care. Independent of maternal care, rifampicin increased the females' faeces production and resulted in lighter eggs and juveniles. By contrast, rifampicin altered none of the other 21 physiological, reproductive and longevity traits measured over the 300 days of a female's lifetime. Overall, these findings reveal that altering the gut microbiota with a large spectrum antibiotic such as rifampicin does not necessarily affect host sociality. They also emphasize that not all animals have evolved a co-dependence with their microbiota and call for caution when generalizing the central role of gut microbes in host biology.

Biogeography and Independent Diversification in the Protist Symbiont Community of Heterotermes tenuis

The eukaryotic microbiome of “lower” termites is highly stable and host-specific. This is due to the mutually obligate nature of the symbiosis and the direct inheritance of protists by proctodeal trophallaxis. However, vertical transmission is occasionally imperfect, resulting in daughter colonies that lack one or more of the expected protist species. This phenomenon could conceivably lead to regional differences in protist community composition within a host species. Here, we have characterized the protist symbiont community of Heterotermes tenuis (Hagen) (Blattodea: Rhinotermitidae) from samples spanning South and Central America. Using light microscopy, single cell isolation, and amplicon sequencing, we report eight species-level protist phylotypes belonging to four genera in the phylum Parabasalia. The diversity and distribution of each phylotype’s 18S rRNA amplicon sequence variants (ASVs) mostly did not correlate with geographical or host genetic distances according to Mantel tests, consistent with the lack of correlation we observed between host genetic and geographical distances. However, the ASV distances of Holomastigotoides Ht3 were significantly correlated with geography while those of Holomastigotoides Ht1 were significantly correlated with host phylogeny. These results suggest mechanisms by which termite-associated protist species may diversify independently of each other and of their hosts, shedding light on the coevolutionary dynamics of this important symbiosis.

The Complete Protist Symbiont Communities of Coptotermes formosanus and Coptotermes gestroi: Morphological and Molecular Characterization of Five New Species

Coptotermes formosanus Shiraki and Coptotermes gestroi (Wasmann) (Blattoidea: Rhinotermitidae) are invasive subterranean termite pest species with a major global economic impact. However, the descriptions of the mutualistic protist communities harbored in their respective hindguts remain fragmentary. The C. formosanus hindgut has long been considered to harbor three protist species, Pseudotrichonympha grassii (Trichonymphida), Holomastigotoides hartmanni, and Cononympha (Spirotrichonympha) leidyi (Spirotrichonymphida), but molecular data have suggested that the diversity may be higher. Meanwhile, the C. gestroi community remains undescribed except for Pseudotrichonympha leei. To complete the characterization of these communities, hindguts of workers from both termite species were investigated using single-cell PCR, microscopy, cell counts, and 18S rRNA amplicon sequencing. The two hosts were found to harbor intriguingly parallel protist communities, each consisting of one Pseudotrichonympha species, two Holomastigotoides species, and two Cononympha species. All protist species were unique to their respective hosts, which last shared a common ancestor ~18 MYA. The relative abundances of protist species in each hindgut differed remarkably between cell count data and 18S rRNA profiles, calling for caution in interpreting species abundances from amplicon data. This study will enable future research in C. formosanus and C. gestroi hybrids, which provide a unique opportunity to study protist community inheritance, compatibility, and potential contribution to hybrid vigor.