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.

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

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

Enhanced Mutation Rate, Relaxed Selection, and the "Domino Effect" are associated with Gene Loss in Blattabacterium, A Cockroach Endosymbiont.

Intracellular endosymbionts have reduced genomes that progressively lose genes at a timescale of tens of million years. We previously reported that gene loss rate is linked to mutation rate in Blattabacterium, however, the mechanisms causing gene loss are not yet fully understood. Here, we carried out comparative genomic analyses on the complete genome sequences of a representative set of 67 Blattabacterium strains, with sizes ranging between 511 and 645 kb. We found that 200 of the 566 analyzed protein-coding genes were lost in at least one lineage of Blattabacterium, with the most extreme case being one gene that was lost independently in 24 lineages. We found evidence for three mechanisms influencing gene loss in Blattabacterium. First, gene loss rates were found to increase exponentially with the accumulation of substitutions. Second, genes involved in vitamin and amino acid metabolism experienced relaxed selection in Cryptocercus and Mastotermes, possibly triggered by their vertically inherited gut symbionts. Third, we found evidence of epistatic interactions among genes leading to a "domino effect" of gene loss within pathways. Our results highlight the complexity of the process of genome erosion in an endosymbiont.

Vicariance and dispersal events inferred from mitochondrial genomes and nuclear genes (18S, 28S) shaped global Cryptocercus distributions.

Cryptocercus Scudder, a genus of wingless, subsocial cockroaches, has low vagility but exhibits a disjunct distribution in eastern and western North America, and in China, South Korea and the Russian Far East. This distribution provides an ideal model for testing hypotheses of vicariance through plate tectonics or other natural barriers versus dispersal across oceans or other natural barriers. We sequenced 45 samples of Cryptocercus to resolve phylogenetic relationships among members of the genus worldwide. We identified four types of tRNA rearrangements among samples from the Qin-Daba Mountains. Our maximum-likelihood and Bayesian phylogenetic trees, based on mitochondrial genomes and nuclear genes (18S, 28S), strongly supported six major lineages of Cryptocercus, which displayed a clear geographical distribution pattern. We used Bayesian molecular dating to estimate the evolutionary timescale of the genus, and reconstructed Cryptocercus ancestral ranges using statistical dispersal-vicariance analysis (S-DIVA) in RASP. Two dispersal events and six vicariance events for Cryptocercus were inferred with high support. The initial vicariance event occurred between American and Asian lineages at 80.5 Ma (95% credibility interval: 60.0-104.7 Ma), followed by one vicariance event within the American lineage 43.8 Ma (95% CI: 32.0-57.5 Ma), and two dispersal 31.9 Ma (95% CI: 25.8-39.5 Ma), 21.7 Ma (95% CI: 17.3-27.1 Ma) plus four vicariance events c. 29.3 Ma, 27.2 Ma, 24.8 Ma and 16.7 Ma within the Asian lineage. Our analyses provide evidence that both vicariance and dispersal have played important roles in shaping the distribution and diversity of these woodroaches.

Shrinking in the dark: Parallel endosymbiont genome erosions are associated with repeated host transitions to an underground life.

Microbial symbioses have had profound impacts on the evolution of animals. Conversely, changes in host biology may impact the evolutionary trajectory of symbionts themselves. Blattabacterium cuenoti is present in almost all cockroach species and enables hosts to subsist on a nutrient-poor diet. To investigate if host biology has impacted Blattabacterium at the genomic level, we sequenced and analyzed 25 genomes from Australian soil-burrowing cockroaches (Blaberidae: Panesthiinae), which have undergone at least seven separate subterranean, subsocial transitions from above-ground, wood-feeding ancestors. We find at least three independent instances of genome erosion have occurred in Blattabacterium strains exclusive to Australian soil-burrowing cockroaches. These shrinkages have involved the repeated inactivation of genes involved in amino acid biosynthesis and nitrogen recycling, the core role of Blattabacterium in the host-symbiont relationship. The most drastic of these erosions have occurred in hosts thought to have transitioned underground the earliest relative to other lineages, further suggestive of a link between gene loss in Blattabacterium and the burrowing behavior of hosts. As Blattabacterium is unable to fulfill its core function in certain host lineages, these findings suggest soil-burrowing cockroaches must acquire these nutrients from novel sources. Our study represents one of the first cases, to our knowledge, of parallel host adaptations leading to concomitant parallelism in their mutualistic symbionts, further underscoring the intimate relationship between these two partners.

Maintenance of essential amino acid synthesis pathways in the Blattabacterium cuenoti symbiont of a wood-feeding cockroach.

In addition to harbouring intestinal symbionts, some animal species also possess intracellular symbiotic microbes. The relative contributions of gut-resident and intracellular symbionts to host metabolism, and how they coevolve are not well understood. Cockroaches and the termite Mastotermes darwiniensis present a unique opportunity to examine the evolution of spatially separated symbionts, as they harbour gut symbionts and the intracellular symbiont Blattabacterium cuenoti. The genomes of B. cuenoti from M. darwiniensis and the social wood-feeding cockroach Cryptocercus punctulatus are each missing most of the pathways for the synthesis of essential amino acids found in the genomes of relatives from non-wood-feeding hosts. Hypotheses to explain this pathway degradation include: (i) feeding on microbes present in rotting wood by ancestral hosts; (ii) the evolution of high-fidelity transfer of gut microbes via social behaviour. To test these hypotheses, we sequenced the B. cuenoti genome of a third wood-feeding species, the phylogenetically distant and non-social Panesthia angustipennis. We show that host wood-feeding does not necessarily lead to degradation of essential amino acid synthesis pathways in B. cuenoti, and argue that ancestral high-fidelity transfer of gut microbes best explains their loss in strains from M. darwiniensis and C. punctulatus.

Coevolution of Metabolic Pathways in Blattodea and Their Blattabacterium Endosymbionts, and Comparisons with Other Insect-Bacteria Symbioses.

Many insects harbor bacterial endosymbionts that supply essential nutrients and enable their hosts to thrive on a nutritionally unbalanced diet. Comparisons of the genomes of endosymbionts and their insect hosts have revealed multiple cases of mutually-dependent metabolic pathways that require enzymes encoded in 2 genomes. Complementation of metabolic reactions at the pathway level has been described for hosts feeding on unbalanced diets, such as plant sap. However, the level of collaboration between symbionts and hosts that feed on more variable diets is largely unknown. In this study, we investigated amino acid and vitamin/cofactor biosynthetic pathways in Blattodea, which comprises cockroaches and termites, and their obligate endosymbiont Blattabacterium cuenoti (hereafter Blattabacterium). In contrast to other obligate symbiotic systems, we found no clear evidence of "collaborative pathways" for amino acid biosynthesis in the genomes of these taxa, with the exception of collaborative arginine biosynthesis in 2 taxa, Cryptocercus punctulatus and Mastotermes darwiniensis. Nevertheless, we found that several gaps specific to Blattabacterium in the folate biosynthetic pathway are likely to be complemented by their host. Comparisons with other insects revealed that, with the exception of the arginine biosynthetic pathway, collaborative pathways for essential amino acids are only observed in phloem-sap feeders. These results suggest that the host diet is an important driving factor of metabolic pathway evolution in obligate symbiotic systems. IMPORTANCE The long-term coevolution between insects and their obligate endosymbionts is accompanied by increasing levels of genome integration, sometimes to the point that metabolic pathways require enzymes encoded in two genomes, which we refer to as "collaborative pathways". To date, collaborative pathways have only been reported from sap-feeding insects. Here, we examined metabolic interactions between cockroaches, a group of detritivorous insects, and their obligate endosymbiont, Blattabacterium, and only found evidence of collaborative pathways for arginine biosynthesis. The rarity of collaborative pathways in cockroaches and Blattabacterium contrasts with their prevalence in insect hosts feeding on phloem-sap. Our results suggest that host diet is a factor affecting metabolic integration in obligate symbiotic systems.

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.

Increased Mutation Rate Is Linked to Genome Reduction in Prokaryotes.

The evolutionary processes that drive variation in genome size across the tree of life remain unresolved. Effective population size (Ne) is thought to play an important role in shaping genome size [1-3]-a key example being the reduced genomes of insect endosymbionts, which undergo population bottlenecks during transmission [4]. However, the existence of reduced genomes in marine and terrestrial prokaryote species with large Ne indicate that genome reduction is influenced by multiple processes [3]. One candidate process is enhanced mutation rate, which can increase adaptive capacity but can also promote gene loss. To investigate evolutionary forces associated with prokaryotic genome reduction, we performed molecular evolutionary and phylogenomic analyses of nine lineages from five bacterial and archaeal phyla. We found that gene-loss rate strongly correlated with synonymous substitution rate (a proxy for mutation rate) in seven of the nine lineages. However, gene-loss rate showed weak or no correlation with the ratio of nonsynonymous/synonymous substitution rate (dN/dS). These results indicate that genome reduction is largely associated with increased mutation rate, while the association between gene loss and changes in Ne is less well defined. Lineages with relatively high dS and dN, as well as smaller genomes, lacked multiple DNA repair genes, providing a proximate cause for increased mutation rates. Our findings suggest that similar mechanisms drive genome reduction in both intracellular and free-living prokaryotes, with implications for developing a comprehensive theory of prokaryote genome size evolution.

Parallel and Gradual Genome Erosion in the Blattabacterium Endosymbionts of Mastotermes darwiniensis and Cryptocercus Wood Roaches.

Almost all examined cockroaches harbor an obligate intracellular endosymbiont, Blattabacterium cuenoti. On the basis of genome content, Blattabacterium has been inferred to recycle nitrogen wastes and provide amino acids and cofactors for its hosts. Most Blattabacterium strains sequenced to date harbor a genome of ∼630 kbp, with the exception of the termite Mastotermes darwiniensis (∼590 kbp) and Cryptocercus punctulatus (∼614 kbp), a representative of the sister group of termites. Such genome reduction may have led to the ultimate loss of Blattabacterium in all termites other than Mastotermes. In this study, we sequenced 11 new Blattabacterium genomes from three species of Cryptocercus in order to shed light on the genomic evolution of Blattabacterium in termites and Cryptocercus. All genomes of Cryptocercus-derived Blattabacterium genomes were reduced (∼614 kbp), except for that associated with Cryptocercus kyebangensis, which comprised 637 kbp. Phylogenetic analysis of these genomes and their content indicates that Blattabacterium experienced parallel genome reduction in Mastotermes and Cryptocercus, possibly due to similar selective forces. We found evidence of ongoing genome reduction in Blattabacterium from three lineages of the C. punctulatus species complex, which independently lost one cysteine biosynthetic gene. We also sequenced the genome of the Blattabacterium associated with Salganea taiwanensis, a subsocial xylophagous cockroach that does not vertically transmit gut symbionts via proctodeal trophallaxis. This genome was 632 kbp, typical of that of nonsubsocial cockroaches. Overall, our results show that genome reduction occurred on multiple occasions in Blattabacterium, and is still ongoing, possibly because of new associations with gut symbionts in some lineages.

An Efficient Strategy Developed for Next-Generation Sequencing of Endosymbiont Genomes Performed Using Crude DNA Isolated from Host Tissues: A Case Study of Blattabacterium cuenoti Inhabiting the Fat Bodies of Cockroaches.

Whole-genome sequencing has emerged as one of the most effective means to elucidate the biological roles and molecular features of obligate intracellular symbionts (endosymbionts). However, the de novo assembly of an endosymbiont genome remains a challenge when host and/or mitochondrial DNA sequences are present in a dataset and hinder the assembly of the genome. By focusing on the traits of genome evolution in endosymbionts, we herein developed and investigated a genome-assembly strategy that consisted of two consecutive procedures: the selection of endosymbiont contigs from an output obtained from a de novo assembly performed using a TBLASTX search against a reference genome, named TBLASTX Contig Selection and Filtering (TCSF), and the iterative reassembling of the genome from reads mapped on the selected contigs, named Iterative Mapping and ReAssembling (IMRA), to merge the contigs. In order to validate this approach, we sequenced two strains of the cockroach endosymbiont Blattabacterium cuenoti and applied this strategy to the datasets. TCSF was determined to be highly accurate and sensitive in contig selection even when the genome of a distantly related free-living bacterium was used as a reference genome. Furthermore, the use of IMRA markedly improved sequence assemblies: the genomic sequence of an endosymbiont was almost completed from a dataset containing only 3% of the sequences of the endosymbiont's genome. The efficiency of our strategy may facilitate further studies on endosymbionts.