An indirect test of inbreeding depression in the termites Reticulitermes flavipes and Reticulitermes virginicus

Pervasive relaxed selection in termite genomes

Genetic changes that enabled the evolution of eusociality have long captivated biologists. More recently, attention has focussed on the consequences of eusociality on genome evolution. Studies have reported higher molecular evolutionary rates in eusocial hymenopteran insects compared with their solitary relatives. To investigate the genomic consequences of eusociality in termites, we analysed nine genomes, including newly sequenced genomes from three non-eusocial cockroaches. Using a phylogenomic approach, we found that termite genomes have experienced lower rates of synonymous substitutions than those of cockroaches, possibly as a result of longer generation times. We identified higher rates of non-synonymous substitutions in termite genomes than in cockroach genomes, and identified pervasive relaxed selection in the former (24-31% of the genes analysed) compared with the latter (2-4%). We infer that this is due to reductions in effective population size, rather than gene-specific effects (e.g. indirect selection of caste-biased genes). We found no obvious signature of increased genetic load in termites, and postulate efficient purging of deleterious alleles at the colony level. Additionally, we identified genomic adaptations that may underpin caste differentiation, such as genes involved in post-translational modifications. Our results provide insights into the evolution of termites and the genomic consequences of eusociality more broadly.

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.

Short and long-term costs of inbreeding in the lifelong-partnership in a termite

Social life and lifelong partner commitments are expected to favor thorough partner choice, as an ill-suited partnership may have long-term consequences, adversely affecting the parents and spanning several cohorts of offspring. Here, we used ~1400 termite incipient colonies to estimate the short- and long-term costs of inbreeding upon the survival of the parents over a 15-month period, their productivity, and the resistance of their offspring toward pathogen pressure. We observed that foundation success was not influenced by the relatedness of partners, but by their levels of microbial load. We showed faster growth in inbred colonies with low levels of microbial load, revealing a potential tradeoff between pathogen defense and offspring production. Yet, inbreeding takes its toll later in colony development when offspring from incipient colonies face pathogen pressure. Although the success of a lifetime partnership is initially determined by the partner’s health, the cost of inbreeding in incipient colonies favors outbred colonies reaching maturity.

Studies of termite colonies over 15 months show that inbred colonies exhibit faster initial growth with low levels of microbial load, but higher mortality toward pathogens later in colony development.

Extensive human-mediated jump dispersal within and across the native and introduced ranges of the invasive termite Reticulitermes flavipes

As native ranges are often geographically structured, invasive species originating from a single source population only carry a fraction of the genetic diversity present in their native range. The invasion process is thus often associated with a drastic loss of genetic diversity resulting from a founder event. However, the fraction of diversity brought to the invasive range may vary under different invasion histories, increasing with the size of the propagule, the number of reintroduction events, and/or the total genetic diversity represented by the various source populations in a multiple-introduction scenario. In this study, we generated a SNP data set for the invasive termite Reticulitermes flavipes from 23 native populations in the eastern United States and six introduced populations throughout the world. Using population genetic analyses and approximate Bayesian computation random forest, we investigated its worldwide invasion history. We found a complex invasion pathway with multiple events out of the native range and bridgehead introductions from the introduced population in France. Our data suggest that extensive long-distance jump dispersal appears common in both the native and introduced ranges of this species, probably through human transportation. Overall, our results show that similar to multiple introduction events into the invasive range, admixture in the native range prior to invasion can potentially favour invasion success by increasing the genetic diversity that is later transferred to the introduced range.

Experimental Warming Reduces Survival, Cold Tolerance, and Gut Prokaryotic Diversity of the Eastern Subterranean Termite, Reticulitermes flavipes (Kollar)

Understanding the effects of environmental disturbances on insects is crucial in predicting the impact of climate change on their distribution, abundance, and ecology. As microbial symbionts are known to play an integral role in a diversity of functions within the insect host, research examining how organisms adapt to environmental fluctuations should include their associated microbiota. In this study, subterranean termites [Reticulitermes flavipes (Kollar)] were exposed to three different temperature treatments characterized as low (15°C), medium (27°C), and high (35°C). Results suggested that pre-exposure to cold allowed termites to stay active longer in decreasing temperatures but caused termites to freeze at higher temperatures. High temperature exposure had the most deleterious effects on termites with a significant reduction in termite survival as well as reduced ability to withstand cold stress. The microbial community of high temperature exposed termites also showed a reduction in bacterial richness and decreased relative abundance of Spirochaetes, Elusimicrobia, and methanogenic Euryarchaeota. Our results indicate a potential link between gut bacterial symbionts and termite's physiological response to environmental changes and highlight the need to consider microbial symbionts in studies relating to insect thermosensitivity.

Female-biased sex allocation and lack of inbreeding avoidance in Cubitermes termites

Sexually reproducing organisms face a strong selective pressure to find a mate and ensure reproduction. An important criterion during mate-selection is to avoid closely related individuals and subsequent potential fitness costs of resulting inbred offspring. Inbreeding avoidance can be active through kin recognition during mate choice, or passive through differential male and female-biased sex ratios, which effectively prevents sib-mating. In addition, sex allocation, or the resources allotted to male and female offspring, can impact mating and reproductive success. Here, we investigate mate choice, sex ratios, and sex allocation in dispersing reproductives (alates) from colonies of the termite Cubitermes tenuiceps. Termites have a short time to select a mate for life, which should intensify any fitness consequences of inbreeding. However, alates did not actively avoid inbreeding through mate choice via kin recognition based on genetic or environmental cues. Furthermore, the majority of colonies exhibited a female-biased sex ratio, and none exhibited a male-bias, indicating that differential bias does not reduce inbreeding. Sex allocation was generally female-biased, as females also were heavier, but the potential fitness effect of this costly strategy remains unclear. The bacterium Wolbachia, known in other insects to parasitically distort sex allocation toward females, was present within all alates. While Wolbachia is commonly associated with termites, parasitism has yet to be demonstrated, warranting further study of the nature of the symbiosis. Both the apparent lack of inbreeding avoidance and potential maladaptive sex allocation implies possible negative effects on mating and fitness.

Natural variation in colony inbreeding does not influence susceptibility to a fungal pathogen in a termite

Reduced genetic diversity through inbreeding can negatively affect pathogen resistance. This relationship becomes more complicated in social species, such as social insects, since the chance of disease transmission increases with the frequency of interactions among individuals. However, social insects may benefit from social immunity, whereby individual physiological defenses may be bolstered by collective-level immune responses, such as grooming or sharing of antimicrobial substance through trophallaxis. We set out to determine whether differences in genetic diversity between colonies of the subterranean termite, Reticulitermes flavipes, accounts for colony survival against pathogens. We sampled colonies throughout the United States (Texas, North Carolina, Maryland, and Massachusetts) and determined the level of inbreeding of each colony. To assess whether genetically diverse colonies were better able to survive exposure to diverse pathogens, we challenged groups of termite workers with two strains of a pathogenic fungus, one local strain present in the soil surrounding sampled colonies and another naïve strain, collected outside the range of this species. We found natural variation in the level of inbreeding between colonies, but this variation did not explain differences in susceptibility to either pathogen. Although the naïve strain was found to be more hazardous than the local strain, colony resistance was correlated between two strains, meaning that colonies had either relatively high or low susceptibility to both strains regardless of their inbreeding coefficient. Overall, our findings may reflect differential virulence between the strains, immune priming of the colonies via prior exposure to the local strain, or a coevolved resistance toward this strain. They also suggest that colony survival may rely more upon additional factors, such as different behavioral response thresholds or the influence of a specific genetic background, rather than the overall genetic diversity of the colony.

Area-Wide Elimination of Subterranean Termite Colonies Using a Novaluron Bait

We investigated the use of termite baiting, a proven system of targeted colony elimination, in an overall area-wide control strategy against subterranean termites. At two field sites, we used microsatellite markers to estimate the total number of Reticulitermes colonies, their spatial partitioning, and breeding structure. Termite pressure was recorded for two years before and after the introduction of Trelona® (active ingredient novaluron) to a large area of one of the sites. Roughly 70% of the colonies in the treatment site that were present at the time of baiting were not found in the site within two months after the introduction of novaluron. Feeding activity of the remaining colonies subsequently ceased over time and new invading colonies were unable to establish within this site. Our study provides novel field data on the efficacy of novaluron in colony elimination of Reticulitermes flavipes, as well as evidence that an area-wide baiting program is feasible to maintain a termite-free area within its native range.

Increased genetic diversity from colony merging in termites does not improve survival against a fungal pathogen

In some species of social insects the increased genetic diversity from having multiple breeders in a colony has been shown to improve pathogen resistance. Termite species typically found colonies from single mated pairs and therefore may lack the flexibility to buffer pathogen pressure with increased genetic diversity by varying the initial number of reproductives. However, they can later increase group diversity through colony merging, resulting in a genetically diverse, yet cohesive, workforce. In this study, we investigate whether the increased group diversity from colony fusion benefits social immunity in the subterranean termite Reticulitermes flavipes. We confirm previous findings that colonies of R. flavipes will readily merge and we show that workers will equally groom nestmates and non-nestmates after merging. Despite this, the survival of these merged colonies was not improved after exposure to a fungal pathogen, but instead leveled to that of the more susceptible or the more resistant colony. Our study brings little support to the hypothesis that colony fusion may improve immunity through an increase of genetic diversity in R. flavipes. Instead, we find that following exposure to a lethal pathogen, one colony is heavily influential to the entire group's survival after merging.

Genetic Analysis of Termite Colonies in Wisconsin

The objective of this study was to document current areas of subterranean termite activity in Wisconsin and to evaluate genetic characteristics of these northern, peripheral colonies. Here, amplified fragment-length polymorphism was used to characterize levels of inbreeding, expected heterozygosity, and percent polymorphism within colonies as well as genetic structure among populations sampled. Genetic analysis revealed two species of termites occur in Wisconsin, Reticulitermes flavipes (Kollar) and Reticulitermes tibialis Banks, both found in the southern half of the state. Colonies of both species in Wisconsin are thought to represent the northern boundary of their current distributions. Measurements of within colony genetic variation showed the proportion of polymorphic loci to be between 52.9-63.9% and expected heterozygosity to range from 0.122-0.189. Consistent with geographical isolation, strong intercolony genetic differences were observed, with over 50% of FST values above 0.25 and the remaining showing moderate levels of genetic differentiation. Combined with low levels of inbreeding in most collection locations (FIS 0.042-0.123), we hypothesize termites were introduced numerous times in the state, likely by anthropogenic means. We discuss the potential effects of these genetic characteristics on successful colony establishment of termites along the northern boundary compared with termites in the core region of their distribution.

Relationship between invasion success and colony breeding structure in a subterranean termite

Factors promoting the establishment and colonization success of introduced populations in new environments constitute an important issue in biological invasions. In this context, the respective role of pre-adaptation and evolutionary changes during the invasion process is a key question that requires particular attention. This study compared the colony breeding structure (i.e. number and relatedness among reproductives within colonies) in native and introduced populations of the subterranean pest termite, Reticulitermes flavipes. We generated and analysed a data set of both microsatellite and mtDNA loci on termite samples collected in three introduced populations, one in France and two in Chile, and in the putative source population of French and Chilean infestations that has recently been identified in New Orleans, LA. We also provided a synthesis combining our results with those of previous studies to obtain a global picture of the variation in breeding structure in this species. Whereas most native US populations are mainly composed of colonies headed by monogamous pairs of primary reproductives, all introduced populations exhibit a particular colony breeding structure that is characterized by hundreds of inbreeding reproductives (neotenics) and by a propensity of colonies to fuse, a pattern shared uniquely with the population of New Orleans. These characteristics are comparable to those of many invasive ants and are discussed to play an important role during the invasion process. Our finding that the New Orleans population exhibits the same breeding structure as its related introduced populations suggests that this native population is pre-adapted to invade new ranges.

Reference-free population genomics from next-generation transcriptome data and the vertebrate-invertebrate gap

In animals, the population genomic literature is dominated by two taxa, namely mammals and drosophilids, in which fully sequenced, well-annotated genomes have been available for years. Data from other metazoan phyla are scarce, probably because the vast majority of living species still lack a closely related reference genome. Here we achieve de novo, reference-free population genomic analysis from wild samples in five non-model animal species, based on next-generation sequencing transcriptome data. We introduce a pipe-line for cDNA assembly, read mapping, SNP/genotype calling, and data cleaning, with specific focus on the issue of hidden paralogy detection. In two species for which a reference genome is available, similar results were obtained whether the reference was used or not, demonstrating the robustness of our de novo inferences. The population genomic profile of a hare, a turtle, an oyster, a tunicate, and a termite were found to be intermediate between those of human and Drosophila, indicating that the discordant genomic diversity patterns that have been reported between these two species do not reflect a generalized vertebrate versus invertebrate gap. The genomic average diversity was generally higher in invertebrates than in vertebrates (with the notable exception of termite), in agreement with the notion that population size tends to be larger in the former than in the latter. The non-synonymous to synonymous ratio, however, did not differ significantly between vertebrates and invertebrates, even though it was negatively correlated with genetic diversity within each of the two groups. This study opens promising perspective regarding genome-wide population analyses of non-model organisms and the influence of population size on non-synonymous versus synonymous diversity.

Preadaptation for parthenogenetic colony foundation in subterranean termites Reticulitermes spp. (Isoptera: Rhinotermitidae)

Thelytokous (all-female producing) parthenogenesis, in some cases, involves reproductive advantages against obligate sexual reproduction. However, the completion of parthenogenesis takes multiple steps without the help of males, and thus preadaptation that meets those requirements will be an important factor for the evolution of parthenogenesis. The Japanese subterranean termite, Reticulitermes speratus, is known to have the ability of parthenogenetic colony foundation, where females that failed to mate with males found colonies cooperatively with partner females and reproduce by parthenogenesis. In this study, we compared the parthenogenetic ability and the colony initiation behavior among six Reticulitermes species in Japan. All species other than R. speratus were not able to reproduce parthenogenetically. Nevertheless, females of these species without the parthenogenetic ability performed homosexual female–female colony initiation and produced eggs without fertilization. In addition, in one species without parthenogenetic reproduction, R. kanmonensis, female–female pair initiated founding behavior as quickly as a heterosexual pair. These results suggest that female–female colony initiation and virgin egg-laying are predominant characters among the genus Reticulitermes and provide a preadaptive condition for parthenogenetic colony foundation in R. speratus.

Asexual queen succession in the subterranean termite Reticulitermes virginicus

Termite colonies are founded by a pair of primary reproductives. In many species, including subterranean termites (family Rhinotermitidae), the primary king and queen can be succeeded by neotenic reproductives that are produced from workers or nymphs within the colony. It is generally believed that these neotenics inbreed within the colony, sometimes for many generations. Here, we show that primary queens of the North American subterranean termite, Reticulitermes virginicus, are replaced by numerous parthenogenetically produced female neotenics. We collected functional female neotenics from five colonies of R. virginicus in North Carolina and Texas, USA. Genetic analysis at eight microsatellite loci showed that 91-100% of the neotenics present within a colony were homozygous at all loci, indicating that they were produced through automictic parthenogenesis with terminal fusion. In contrast, workers, soldiers and alates were almost exclusively sexually produced by mating between the female neotenics and a single king. This is the second termite species shown to undergo asexual queen succession, a system first described in the Japanese species, Reticulitermes speratus. Thus, the conditional use of sexual and asexual reproduction to produce members of different castes may be widespread within Reticulitermes and possibly other subterranean termites.

Colony breeding system influences cuticular bacterial load of Formosan subterranean termite (Isoptera: Rhinotermitidae) workers

The goal of this study was to test whether the breeding system and/or the degree of inbreeding of field colonies of the Formosan subterranean termite, Coptotermes formosanus, Shiraki (Isoptera: Rhinotermitidae) influences bacterial load on the cuticle of foraging workers. We enumerated bacterial load on the cuticle of groups of workers foraging in 20 inground monitoring stations surrounding the French Market in New Orleans, LA, and identified bacteria species using 16S rRNA gene sequencing. We used microsatellite genotyping to assign the 20 worker groups to seven simple family colonies (headed by a single pair of reproductives) and four extended family colonies (headed by multiple inbreeding reproductives) with a wide range of degrees of inbreeding. Workers from extended family colonies had a higher bacterial load than those from simple family colonies; however, bacterial load was not significantly correlated to the degree of inbreeding, possibly because of confounding factors in colony life history, such as age and/or size of colonies. Colonies with high bacterial load did not have a higher proportion of entomopathogens, and thus, bacterial load is not necessarily an indicator for disease risk. The majority of bacteria cultured from the cuticle of termites were soil bacteria with no known pathology against termites.

Biology of subterranean termites: insights from molecular studies of Reticulitermes and Coptotermes

Molecular genetic techniques have made contributions to studies on subterranean termites at all levels of biological organization. Most of this work has focused on Reticulitermes and Coptotermes, two ecologically and economically important genera. DNA sequence data have significantly improved our understanding of the systematics and taxonomy of these genera. Techniques of molecular biology have provided important new insights into the process of caste differentiation. Population genetic markers, primarily microsatellites, have furthered our understanding of the life history, population biology, community ecology, and invasion biology of subterranean termites. Recent results on the behavioral ecology of subterranean termites reveal a picture different from long-held views, especially those concerning colony breeding structures and foraging ranges. As additional molecular tools and genomic resources become available, and as more subterranean termite researchers incorporate molecular techniques into their approaches, we can expect accelerating advances in all aspects of the biology of this group.

Size and heterozygosity influence partner selection in the Formosan subterranean termite

In monogamous species that exhibit extensive biparental investment, such as termites, both sexes are predicted to be selective when choosing a mate. Size-related traits are expected to be important in partner selection because the fat reserves of the colony founders sustain the incipient colony. Partner relatedness and heterozygosity determine the degree of inbreeding and genetic diversity within the colony and may thus also influence partner selection. To test these predictions, we investigated whether phenotypic and genetic traits influence mate choice and/or competitive advantage during pair formation of Formosan subterranean termites, Coptotermes formosanus (Isoptera: Rhinotermitidae). Pair formation in termites normally occurs within a short period after swarming when alates form tandem pairs on the ground. Alates were collected from 5 light trap samples in the French Quarter of New Orleans, LA. From each sample, both tandem pairs and single individuals were collected and their sex, body weights, and head widths were recorded. Pairwise relatedness and individual levels of heterozygosity were determined by microsatellite genotyping. Males in tandem pairs with females had a significantly larger head width than males that did not form tandem pairs. Weights as well as head widths of tandem running partners were positively correlated. For the majority of the samples, relatedness between tandem partners did not differ from the relatedness to members of the other tandem pairs. Thus, no kin discrimination occurred during tandem running. However, females engaged in tandem running had a higher degree of heterozygosity than females that remained single. These findings suggest partner selection and/or competitive advantage based on size-related phenotypic parameters and genetic diversity. The pairing advantage of heterozygous females might explain previous findings of sex-biased alate production depending on the degree of inbreeding in colonies of several species of the genus Coptotermes.

Molecular-genetic analyses of dispersal and breeding behaviour in the Australian termite Coptotermes lacteus: evidence for non-random mating in a swarm-dispersal mating system

We used microsatellite DNA markers to infer the dispersal and breeding behaviour of Coptotermes lacteus, a termite whose large mounds are a conspicuous feature of Australia’s central east coast. We genotyped a subsample of neuter offspring for each of 38 colonies sampled over two spatially separated populations, one in a natural forest and the other in an exotic radiata pine plantation. All colonies showed offspring genotype frequencies consistent with a single reproductive pair. This result confirms that stable monogamy is the normal breeding arrangement for this species and that multi-reproductive colonies are rare. The two study populations were significantly differentiated and the distance separating them (~150 km) is therefore an effective constraint on gene flow. The populations themselves, however, were not noticeably subdivided above the level of colony. This lack of within-population viscosity is unexpected for weakly dispersing species and suggests that local gamete dispersal is in fact quite effective in C. lacteus. Nonetheless, dispersing sexuals do not appear to mate randomly. Instead, all four microsatellite loci are deficient in heterozygotes, indicating that populations are substantially inbred, irrespective of habitat. Evidence from hierarchical F-statistics, spatial genetic autocorrelation and relatedness calculations suggests that deviations from Hardy–Weinberg equilibrium may result from either a preference for non-sibling relatives over totally unrelated mates, or from random mating with viscosity – though evidence for the latter hypothesis was not detected. These findings suggest that swarm-dispersal mating systems, usually considered to produce outbreeding and panmixia, can instead involve a notable degree of non-random mating.

Inbreeding and disease resistance in a social insect: effects of heterozygosity on immunocompetence in the termite Zootermopsis angusticollis

Recent research has shown that low genetic variation in individuals can increase susceptibility to infection and group living may exacerbate pathogen transmission. In the eusocial diploid termites, cycles of outbreeding and inbreeding characterizing basal species can reduce genetic variation within nestmates during the life of a colony, but the relationship of genetic heterogeneity to disease resistance is poorly understood. Here we show that, one generation of inbreeding differentially affects the survivorship of isolated and grouped termites (Zootermopsis angusticollis) depending on the nature of immune challenge and treatment. Inbred and outbred isolated and grouped termites inoculated with a bacterial pathogen, exposed to a low dose of fungal pathogen or challenged with an implanted nylon monofilament had similar levels of immune defence. However, inbred grouped termites exposed to a relatively high concentration of fungal conidia had significantly greater mortality than outbred grouped termites. Inbred termites also had significantly higher cuticular microbial loads, presumably due to less effective grooming by nestmates. Genetic analyses showed that inbreeding significantly reduced heterozygosity and allelic diversity. Decreased heterozygosity thus appeared to increase disease susceptibility by affecting social behaviour or some other group-level process influencing infection control rather than affecting individual immune physiology.