Carbohydrate-active enzymes revealed in Coptotermes formosanus (Isoptera: Rhinotermitidae) transcriptome

Structural and biochemical analysis of family 92 carbohydrate-binding modules uncovers multivalent binding to β-glucans

Carbohydrate-binding modules (CBMs) are non-catalytic proteins found appended to carbohydrate-active enzymes. Soil and marine bacteria secrete such enzymes to scavenge nutrition, and they often use CBMs to improve reaction rates and retention of released sugars. Here we present a structural and functional analysis of the recently established CBM family 92. All proteins analysed bind preferentially to β-1,6-glucans. This contrasts with the diversity of predicted substrates among the enzymes attached to CBM92 domains. We present crystal structures for two proteins, and confirm by mutagenesis that tryptophan residues permit ligand binding at three distinct functional binding sites on each protein. Multivalent CBM families are uncommon, so the establishment and structural characterisation of CBM92 enriches the classification database and will facilitate functional prediction in future projects. We propose that CBM92 proteins may cross-link polysaccharides in nature, and might have use in novel strategies for enzyme immobilisation.

Uncovering the Chemosensory System of a Subterranean Termite, Odontotermes formosanus (Shiraki) (Isoptera: Termitidae): Revealing the Chemosensory Genes and Gene Expression Patterns

Termites are eusocial insects. Chemical signals between colony members are crucial to the smooth running of colony operations, but little is known about their olfactory system and the roles played by various chemosensory genes in this process. Chemosensory genes are involved in basic olfactory perception in insects. Odontotermes formosanus (Shiraki) is one of the most damaging pests to agricultural crops, forests, and human-made structures. To better understand the olfactory system and the genes involved in olfactory processing in O. formosanus, we produced a transcriptome of worker termites. In this study, we identified 13 OforOBPs, 1 OforCSP, 15 OforORs, 9 OforGRs, and 4 OforSNMPs. Multiple sequence alignments were used in the phylogenetic study, which included data from other termite species and a wide variety of insect species. Moreover, we also investigated the mRNA expression levels using qRT-PCR. The significantly high expression levels of OforCSP1, OforOBP2, OforOR1, and OforSNMP1 suggest that these genes may play important roles in olfactory processing in termite social behavior, including caste differentiation, nestmate and non-nestmate discrimination, and the performance of colony operations among members. Our research establishes a foundation for future molecular-level functional studies of chemosensory genes in O. formosanus, which might lead to the identification of novel targets for termite integrated pest management.

Trichoderma-derived pentapeptides from the infected nest mycobiome of the subterranean termite Coptotermes testaceus

Termites live in a dynamic environment where colony health is strongly influenced by surrounding microbes. However, little is known about the mycobiomes of lower termites and their nests, and how these change in response to disease. Here we compared the individual and nest mycobiomes of a healthy subterranean termite colony (Coptotermes testaceus) to one infected and ultimately eradicated by a fungal pathogen. We identified Trichoderma species in the materials of both nests, but they were also abundant in the infected termites. Methanolic extracts of Trichoderma sp. FHG000531, isolated from the infected nest, were screened for secondary metabolites by UHPLC‐HR MS/MS‐guided molecular networking. We identified many bioactive compounds with potential roles in the eradication of the infected colony, as well as a cluster of six unknown peptides. The novel peptide FE011 was isolated and characterized by NMR spectroscopy. The function of this novel peptide family as well as the role of Trichoderma species in dying termite colonies therefore requires further investigation.

Three endogenous cellulases from termite, Reticulitermes speratus KMT001

Among termites, lower termites need symbiotic microorganisms in the digestive tract for digestion and cellulose metabolism. In this symbiotic relationship, the decomposition of cellulose is initiated by endoglucanase in termite salivary glands and completed by β-glycosidase of symbiotic microorganisms in the hindgut. The expression of β-glycosidase in lower termites has been reported in recent studies. The expression of two endoglucanases and one β-glycosidase gene related to cellulose degradation was identified in Reticulitermes speratus, a lower termite, through transcriptomic analysis. The proposed enzyme activities of three identified cellulose degradation genes were confirmed by heterologous expression in Escherichia coli. In addition to the endoglucanase expressed in the salivary gland, additional endoglucanase and β-glycosidase genes suggest that R. speratus performs the overall cellulose digestion using its own enzymes at all stages.

A Century of Synergy in Termite Symbiosis Research: Linking the Past with New Genomic Insights

Termites have long been studied for their symbiotic associations with gut microbes. In the late nineteenth century, this relationship was poorly understood and captured the interest of parasitologists such as Joseph Leidy; this research led to that of twentieth-century biologists and entomologists including Cleveland, Hungate, Trager, and Lüscher. Early insights came via microscopy, organismal, and defaunation studies, which led to descriptions of microbes present, descriptions of the roles of symbionts in lignocellulose digestion, and early insights into energy gas utilization by the host termite. Focus then progressed to culture-dependent microbiology and biochemical studies of host-symbiont complementarity, which revealed specific microhabitat requirements for symbionts and noncellulosic mechanisms of symbiosis (e.g., N₂ fixation). Today, knowledge on termite symbiosis has accrued exponentially thanks to omic technologies that reveal symbiont identities, functions, and interdependence, as well as intricacies of host-symbiont complementarity. Moving forward, the merging of classical twentieth-century approaches with evolving omic tools should provide even deeper insights into host-symbiont interplay.

Insights into the Gryllus bimaculatus Immune-Related Transcriptomic Profiling to Combat Naturally Invading Pathogens

Natural pathogen pressure is an important factor that shapes the host immune defense mechanism. The current study primarily aimed to explore the molecular basis of the natural immune defense mechanism of a sporadic pest, Gryllus bimaculatus, during swarming by constructing cDNA libraries of the female mid-gut, male mid-gut, testes, and ovaries. The Illumina HiSeq platform generated an average of 7.9 G, 11.77 G, 10.07 G, and 10.07 G bases of outputs from the male mid-gut, female mid-gut, testes, and ovaries and libraries, respectively. The transcriptome of two-spotted field crickets was assembled into 233,172 UniGenes, which yielded approximately 163.58 million reads. On the other hand, there were 43,055 genes in common that were shared among all the biological samples. Gene Ontology analysis successfully annotated 492 immune-related genes, which comprised mainly Pattern Recognition Receptors (62 genes), Signal modulators (57 genes), Signal transduction (214 genes), Effectors (36 genes), and another immune-related 123 genes. In summary, the identified wide range of immune-related genes from G. bimaculatus indicates the existence of a sophisticated and specialized broad spectrum immune mechanism against invading pathogens, which provides, for the first time, insights into the molecular mechanism of disease resistance among two-spotted field crickets.

Neotropical termite microbiomes as sources of novel plant cell wall degrading enzymes

In this study, we used shotgun metagenomic sequencing to characterise the microbial metabolic potential for lignocellulose transformation in the gut of two colonies of Argentine higher termite species with different feeding habits, Cortaritermes fulviceps and Nasutitermes aquilinus. Our goal was to assess the microbial community compositions and metabolic capacity, and to identify genes involved in lignocellulose degradation. Individuals from both termite species contained the same five dominant bacterial phyla (Spirochaetes, Firmicutes, Proteobacteria, Fibrobacteres and Bacteroidetes) although with different relative abundances. However, detected functional capacity varied, with C. fulviceps (a grass-wood-feeder) gut microbiome samples containing more genes related to amino acid metabolism, whereas N. aquilinus (a wood-feeder) gut microbiome samples were enriched in genes involved in carbohydrate metabolism and cellulose degradation. The C. fulviceps gut microbiome was enriched specifically in genes coding for debranching- and oligosaccharide-degrading enzymes. These findings suggest an association between the primary food source and the predicted categories of the enzymes present in the gut microbiomes of each species. To further investigate the termite microbiomes as sources of biotechnologically relevant glycosyl hydrolases, a putative GH10 endo-β-1,4-xylanase, Xyl10E, was cloned and expressed in Escherichia coli. Functional analysis of the recombinant metagenome-derived enzyme showed high specificity towards beechwood xylan (288.1 IU/mg), with the optimum activity at 50 °C and a pH-activity range from 5 to 10. These characteristics suggest that Xy110E may be a promising candidate for further development in lignocellulose deconstruction applications.

RNA interference of endoglucanases in the formosan subterranean termite Coptotermes formosanus shiraki (Blattodea: Rhinotermitidae) by dsRNA injection or ingestion

Termites obtain energy and nutrition from wood and wood-related materials by utilizing endogenous and symbiotic cellulases. Endoglucanase is one of the key cellulases in cellulose digestion. Previous studies have shown that the inhibition of the cellulase enzyme system would be a plausible approach for termite control. In the present study, we studied the effect of RNAi on termites by targeting a conserved region of five endoglucanase genes from Coptotermes formosanus (CfEGs). Both dsRNA injection and oral delivery resulted in significant gene silencing of CfEGs and consequently led to mortality, reduced enzyme activity, and reduced weight compared to control worker termites. An injection dose of 150 ng and a feeding dose of 2 μg/cm² provided for the best RNAi efficiency. dsCfEG was further combined with flufenoxuron, an insect growth regulator used to manage/suppress subterranean termites, and when fed to workers, caused a lower enzyme activity compared to the dsCfEG- or flufenoxuron-only treatment. The weight loss (∼0.598 mg) and mortality (∼28%) observed in the combined dsCfEG and flufenoxuron treatment differed significantly from those observed in the flufenoxuron-only treatment (∼0.208 mg and ∼16%, respectively). Although the effects of these dsCfEG treatments on mortality were insufficient to serve as termiticides, dsCfEGs could be used in combination with other treatments to increase efficacy. This study provides a research basis for the use of RNAi in termiticides.

Functional metagenomics reveals abundant polysaccharide-degrading gene clusters and cellobiose utilization pathways within gut microbiota of a wood-feeding higher termite

Plant cell-wall polysaccharides constitute the most abundant but recalcitrant organic carbon source in nature. Microbes residing in the digestive tract of herbivorous bilaterians are particularly efficient at depolymerizing polysaccharides into fermentable sugars and play a significant support role towards their host's lifestyle. Here, we combine large-scale functional screening of fosmid libraries, shotgun sequencing, and biochemical assays to interrogate the gut microbiota of the wood-feeding "higher" termite Globitermes brachycerastes. A number of putative polysaccharide utilization gene clusters were identified with multiple fibrolytic genes. Our large-scale functional screening of 50,000 fosmid clones resulted in 464 clones demonstrating plant polysaccharide-degrading activities, including 267 endoglucanase-, 24 exoglucanase-, 72 β-glucosidase-, and 101 endoxylanase-positive clones. We sequenced 173 functionally active clones and identified ~219 genes encoding putative carbohydrate-active enzymes (CAZymes) targeting cellulose, hemicellulose and pectin. Further analyses revealed that 68 of 154 contigs encode one or more CAZyme, which includes 35 examples of putative saccharolytic operons, suggesting that clustering of CAZymes is common in termite gut microbial inhabitants. Biochemical characterization of a representative xylanase cluster demonstrated that constituent enzymes exhibited complementary physicochemical properties and saccharolytic capabilities. Furthermore, diverse cellobiose-metabolizing enzymes include β-glucosidases, cellobiose phosphorylases, and phopho-6-β-glucosidases were identified and functionally verified, indicating that the termite gut micro-ecosystem utilizes diverse metabolic pathways to interconnect hydrolysis and central metabolism. Collectively, these results provide an in-depth view of the adaptation and digestive strategies employed by gut microbiota within this tiny-yet-efficient host-associated ecosystem.

De novo transcriptome assembly of the bamboo snout beetle Cyrtotrachelus buqueti reveals ability to degrade lignocellulose of bamboo feedstock

BACKGROUND

The bamboo weevil Cyrtotrachelus buqueti, which is considered a pest species, damages bamboo shoots via its piercing-sucking mode of feeding. C. buqueti is well known for its ability to transform bamboo shoot biomass into nutrients and energy for growth, development and reproduction with high specificity and efficacy of bioconversion. Woody bamboo is a perennial grass that is a potential feedstock for lignocellulosic biomass because of its high growth rate and lignocellulose content. To verify our hypothesis that C. buqueti efficiently degrades bamboo lignocellulose, we assessed the bamboo lignocellulose-degrading ability of this insect through RNA sequencing for identifying a potential route for utilisation of bamboo biomass.

RESULTS

Analysis of carbohydrate-active enzyme (CAZyme) family genes in the developmental transcriptome of C. buqueti revealed 1082 unigenes, including 55 glycoside hydrolases (GH) families containing 309 GHs, 51 glycosyltransferases (GT) families containing 329 GTs, 8 carbohydrate esterases (CE) families containing 174 CEs, 6 polysaccharide lyases (PL) families containing 11 PLs, 8 auxiliary activities (AA) families containing 131 enzymes with AAs and 17 carbohydrate-binding modules (CBM) families containing 128 CBMs. We used weighted gene co-expression network analysis to analyse developmental RNA sequencing data, and 19 unique modules were identified in the analysis. Of these modules, the expression of MEyellow module genes was unique and the module included numerous CAZyme family genes. CAZyme genes in this module were divided into two groups depending on whether gene expression was higher in the adult/larval stages or in the egg/pupal stages. Enzyme assays revealed that cellulase activity was highest in the midgut whereas lignin-degrading enzyme activity was highest in the hindgut, consistent with findings from intestinal gene expression studies. We also analysed the expression of CAZyme genes in the transcriptome of C. buqueti from two cities and found that several genes were also assigned to CAZyme families. The insect had genes and enzymes associated with lignocellulose degradation, the expression of which differed with developmental stage and intestinal region.

CONCLUSION

Cyrtotrachelus buqueti exhibits lignocellulose degradation-related enzymes and genes, most notably CAZyme family genes. CAZyme family genes showed differences in expression at different developmental stages, with adults being more effective at cellulose degradation and larvae at lignin degradation, as well as at different regions of the intestine, with the midgut being more cellulolytic than the hindgut. This degradative system could be utilised for the bioconversion of bamboo lignocellulosic biomass.

Molecular Characterization and Potential Synthetic Applications of GH1 β-Glucosidase from Higher Termite Microcerotermes annandalei

A novel β-glucosidase from higher termite Microcerotermes annandalei (MaBG) was obtained via a screening method targeting β-glucosidases with increased activities in the presence of glucose. The purified natural MaBG showed a subunit molecular weight of 55 kDa and existed in a native form as a dimer without any glycosylation. Gene-specific primers designed from its partial amino acid sequences were used to amplify the corresponding 1,419-bp coding sequence of MaBG which encodes a 472-amino acid glycoside hydrolase family 1 (GH1) β-glucosidase. When expressed in Komagataella pastoris, the recombinant MaBG appeared as a ~ 55-kDa protein without glycosylation modifications. Kinetic parameters as well as the lack of secretion signal suggested that MaBG is an intracellular enzyme and not involved in cellulolysis. The hydrolytic activities of MaBG were enhanced in the presence of up to 3.5-4.5 M glucose, partly due to its strong transglucosylation activity, which suggests its applicability in biosynthetic processes. The potential synthetic activities of the recombinant MaBG were demonstrated in the synthesis of para-nitrophenyl-β-D-gentiobioside via transglucosylation and octyl glucoside via reverse hydrolysis. The information obtained from this study has broadened our insight into the functional characteristics of this variant of termite GH1 β-glucosidase and its applications in bioconversion and biotechnology.

dsRNA Injection Successfully Inhibited Two Endogenous β-Glucosidases in Coptotermes formosanus (Isoptera: Rhinotermitidae)

Cellulose digestion is an essential process of termites, and it is accomplished by three types of cellulases. β-Glucosidase (BG), one of the critical cellulases responsible for cellulose degradation and glucose production, has been considered as a potential target for pest management strategies. Previous experiments identified two new endogenous BG homologs, CfBG-Ia and CfBG-Ib, in the digestive system of Coptotermes formosanus Shiraki (Isoptera: Rhinotermitidae). The objectives of this study were to assess the impact of RNA interference on CfBG-Ia and CfBG-Ib expression and on termite survival. We tested the expression profiles of worker termites which were injected with gene-specific double-stranded RNA (dsRNA, targeting one gene at a time) and a dsRNA cocktail (targeting CfBG-Ia and CfBG-Ib simultaneously). The expression of CfBG-Ib showed a sharp decline in both dsCfBG-Ib and dsRNA cocktail treatments. The expression of CfBG-Ia reduced quickly and significantly in the dsRNA cocktail treatment; while in dsCfBG-Ia treatment, it decreased on the fifth day. Results showed that treatment with the dsRNA cocktail caused greater inhibition of the transcript expression and a shorter response time. However, the expression of nontarget BG homologs was increased as the target BG homologs were being repressed during the testing period in dsRNA cocktail treatment. These results demonstrate that targeting cellulase-coding genes may be a potential strategy to inhibit termite digestion process, or at least dsRNA cocktails serve as a means for identifying the most susceptible target gene families or biological processes.

Characterization of a laccase from a wood-feeding termite, Coptotermes formosanus

Coptotermes formosanus Shiraki is a wood-feeding termite which secretes a series of lignolytic and cellulolytic enzymes for woody biomass degradation. However, the lignin modification mechanism in the termite is largely elusive, and the characteristics of most lignolytic enzymes in termites remain unknown. In this study, a laccase gene lac1 from C. formosanus was heterogeneously expressed in insect Sf9 cells. The purified Lac1 showed strong activities toward hydroquinone (305 mU/mg) and 2,6-dimethoxyphenol (2.9 mU/mg) with low K_m values, but not veratryl alcohol or 2,2'-azino-bis (3-ethylbenzothiazoline-6-sulfonic acid). Lac1 could function well from pH 4.5 to 7.5, and its activity was significantly inhibited by H₂ O₂ at above 4.85 mmol/L (P < 0.01). In addition, the lac1 gene was found to be mainly expressed in the salivary glands and foregut of C. formosanus, and seldom in the midgut or hindgut. These findings suggested that Lac1 is a phenol-oxidizing laccase like RflacA and RflacB from termite Reticulitermes flavipes, except that Lac1 was found to be more efficient in phenol oxidation, and it did not require H₂ O₂ for its function. It is suspected that this kind of termite laccase might only be able to directly oxidize low redox-potential substrates, and the high redox-potential groups in lignin might be oxidized by other enzymes in the termite or by using the Fenton reaction.

Cross-reaction between Formosan termite (Coptotermes formosanus) proteins and cockroach allergens

Cockroach allergens can lead to serious allergy and asthma symptoms. Termites are evolutionarily related to cockroaches, cohabitate in human dwellings, and represent an increasing pest problem in the United States. The Formosan subterranean termite (Coptotermes formosanus) is one of the most common species in the southern United States. Several assays were used to determine if C. formosanus termite proteins cross-react with cockroach allergens. Expressed sequence tag and genomic sequencing results were searched for homology to cockroach allergens using BLAST 2.2.21 software. Whole termite extracts were analyzed by mass-spectrometry, immunoassay with IgG and scFv antibodies to cockroach allergens, and human IgE from serum samples of cockroach allergic patients. Expressed sequence tag and genomic sequencing results indicate greater than 60% similarity between predicted termite proteins and German and American cockroach allergens, including Bla g 2/Per a 2, Bla g 3/Per a 3, Bla g 5, Bla g 6/Per a 6, Bla g 7/Per a 7, Bla g 8, Per a 9, and Per a 10. Peptides from whole termite extract were matched to those of the tropomyosin (Bla g 7), arginine kinase (Per a 9), and myosin (Bla g 8) cockroach allergens by mass-spectrometry. Immunoblot and ELISA testing revealed cross-reaction between several proteins with IgG and IgE antibodies to cockroach allergens. Several termite proteins, including the hemocyanin and tropomyosin orthologs of Blag 3 and Bla g 7, were shown to crossreact with cockroach allergens. This work presents support for the hypothesis that termite proteins may act as allergens and the findings could be applied to future allergen characterization, epitope analysis, and clinical studies.

Polysaccharide Utilization Loci: Fueling Microbial Communities

The complex carbohydrates of terrestrial and marine biomass represent a rich nutrient source for free-living and mutualistic microbes alike. The enzymatic saccharification of these diverse substrates is of critical importance for fueling a variety of complex microbial communities, including marine, soil, ruminant, and monogastric microbiota. Consequently, highly specific carbohydrate-active enzymes, recognition proteins, and transporters are enriched in the genomes of certain species and are of critical importance in competitive environments. In Bacteroidetes bacteria, these systems are organized as polysaccharide utilization loci (PULs), which are strictly regulated, colocalized gene clusters that encode enzyme and protein ensembles required for the saccharification of complex carbohydrates. This review provides historical perspectives and summarizes key findings in the study of these systems, highlighting a critical shift from sequence-based PUL discovery to systems-based analyses combining reverse genetics, biochemistry, enzymology, and structural biology to precisely illuminate the molecular mechanisms underpinning PUL function. The ecological implications of dynamic PUL deployment by key species in the human gastrointestinal tract are explored, as well as the wider distribution of these systems in other gut, terrestrial, and marine environments.

Transcriptome response to elevated atmospheric CO2 concentration in the Formosan subterranean termite, Coptotermes formosanus Shiraki (Isoptera: Rhinotermitidae)

Background

Carbon dioxide (CO₂) is a pervasive chemical stimulus that plays a critical role in insect life, eliciting behavioral and physiological responses across different species. High CO₂ concentration is a major feature of termite nests, which may be used as a cue for locating their nests. Termites also survive under an elevated CO₂ concentration. However, the mechanism by which elevated CO₂ concentration influences gene expression in termites is poorly understood.

Methods

To gain a better understanding of the molecular basis involved in the adaptation to CO₂ concentration, a transcriptome of Coptotermes formosanus Shiraki was constructed to assemble the reference genes, followed by comparative transcriptomic analyses across different CO₂ concentration (0.04%, 0.4%, 4% and 40%) treatments.

Results

(1) Based on a high throughput sequencing platform, we obtained approximately 20 GB of clean data and revealed 189,421 unigenes, with a mean length and an N50 length of 629 bp and 974 bp, respectively. (2) The transcriptomic response of C. formosanus to elevated CO₂ levels presented discontinuous changes. Comparative analysis of the transcriptomes revealed 2,936 genes regulated among 0.04%, 0.4%, 4% and 40% CO₂ concentration treatments, 909 genes derived from termites and 2,027 from gut symbionts. Genes derived from termites appears selectively activated under 4% CO₂ level. In 40% CO₂ level, most of the down-regulated genes were derived from symbionts. (3) Through similarity searches to data from other species, a number of protein sequences putatively involved in chemosensory reception were identified and characterized in C. formosanus, including odorant receptors, gustatory receptors, ionotropic receptors, odorant binding proteins, and chemosensory proteins.

Discussion

We found that most genes associated with carbohydrate metabolism, energy metabolism, and genetic information processing were regulated under different CO₂ concentrations. Results suggested that termites adapt to ∼4% CO₂ level and their gut symbionts may be killed under high CO₂ level. We anticipate that our findings provide insights into the transcriptome dynamics of CO₂ responses in termites and form the basis to gain a better understanding of regulatory networks.

Comparative Analysis of Transcriptomes from Secondary Reproductives of Three Reticulitermes Termite Species

Termites are eusocial insects related to cockroaches that feed on lignocellulose. These insects are key species in ecosystems since they recycle a large amount of nutrients but also are pests, exerting major economic impacts. Knowledge on the molecular pathways underlying reproduction, caste differentiation or lignocellulose digestion would largely benefit from additional transcriptomic data. This study focused on transcriptomes of secondary reproductive females (nymphoid neotenics). Thirteen transcriptomes were used: 10 of Reticulitermes flavipes and R. grassei sequenced from a previous study, and two transcriptomes of R. lucifugus sequenced for the present study. After transcriptome assembly and read mapping, we examined interspecific variations of genes expressed by termites or gut microorganisms. A total of 18,323 orthologous gene clusters were detected. Functional annotation and taxonomic assignment were performed on a total of 41,287 predicted contigs in the three termite species. Between the termite species studied, functional categories of genes were comparable. Gene ontology (GO) terms analysis allowed the discovery of 9 cellulases and a total of 79 contigs potentially involved in 11 enzymatic activities used in wood metabolism. Altogether, results of this study illustrate the strong potential for the use of comparative interspecific transcriptomes, representing a complete resource for future studies including differentially expressed genes between castes or SNP analysis for population genetics.

Identification and Characterization of Two Endogenous β-Glucosidases from the Termite Coptotermes formosanus

Coptotermes formosanus is a well-known wood-feeding termite that can degrade lignocellulose polysaccharides efficiently with its unique multi-enzyme catalysis system. β-glucosidase (BG) is one of the important cellulases in its enzyme system. However, there may present multiple endogenous BGs in termite digestive system for various properties and functions. This study aims to characterize two BG homologs and reveal their potential coordinative effect. In this study, two endogenous BG homologs (CfGlu1B and CfGlu1C) from C. formosanus showed different substrate specificity. CfGlu1B favors cellobiose while CfGlu1C favors sucrose. Besides, CfGlu1C exhibited higher alkali resistance than CfGlu1B. Kinetic analysis revealed that CfGlu1B enzyme's activity toward p-NP-β-D-glucopyranoside (p-NPG) was higher than that of CfGlu1C, and the difference mainly attributes to the turnover number (K cat). In addition, the activity assay showed significant synergistic action of CfGlu1B and CfGlu1C in degrading D-lactosum. For effect of metals, Cu(2+) inhibited both enzymes and Ca(2+) increased the activity of CfGlu1C but not CfGlu1B. Site-directed mutagenesis analysis indicated that both enzymes lost activities when residues E190 of CfGlu1B and E168 of CfGlu1C were mutated to alanine, respectively, which were essential active centers of the GHF1 enzymes. Moreover, mutation H252N increased the activity of enzyme CfGlu1C by 2.1-fold. This study implies interesting possibilities for better practical biotechnological use in green energy production.

How insects overcome two-component plant chemical defence: plant β-glucosidases as the main target for herbivore adaptation

Insect herbivory is often restricted by glucosylated plant chemical defence compounds that are activated by plant β-glucosidases to release toxic aglucones upon plant tissue damage. Such two-component plant defences are widespread in the plant kingdom and examples of these classes of compounds are alkaloid, benzoxazinoid, cyanogenic and iridoid glucosides as well as glucosinolates and salicinoids. Conversely, many insects have evolved a diversity of counteradaptations to overcome this type of constitutive chemical defence. Here we discuss that such counter-adaptations occur at different time points, before and during feeding as well as during digestion, and at several levels such as the insects’ feeding behaviour, physiology and metabolism. Insect adaptations frequently circumvent or counteract the activity of the plant β-glucosidases, bioactivating enzymes that are a key element in the plant’s two-component chemical defence. These adaptations include host plant choice, non-disruptive feeding guilds and various physiological adaptations as well as metabolic enzymatic strategies of the insect’s digestive system. Furthermore, insect adaptations often act in combination, may exist in both generalists and specialists, and can act on different classes of defence compounds. We discuss how generalist and specialist insects appear to differ in their ability to use these different types of adaptations: in generalists, adaptations are often inducible, whereas in specialists they are often constitutive. Future studies are suggested to investigate in detail how insect adaptations act in combination to overcome plant chemical defences and to allow ecologically relevant conclusions.

Omic research in termites: an overview and a roadmap

Many recent breakthroughs in our understanding of termite biology have been facilitated by "omics" research. Omic science seeks to collectively catalog, quantify, and characterize pools of biological molecules that translate into structure, function, and life processes of an organism. Biological molecules in this context include genomic DNA, messenger RNA, proteins, and other biochemicals. Other permutations of omics that apply to termites include sociogenomics, which seeks to define social life in molecular terms (e.g., behavior, sociality, physiology, symbiosis, etc.) and digestomics, which seeks to define the collective pool of host and symbiont genes that collaborate to achieve high-efficiency lignocellulose digestion in the termite gut. This review covers a wide spectrum of termite omic studies from the past 15 years. Topics covered include a summary of terminology, the various kinds of omic efforts that have been undertaken, what has been revealed, and to a degree, what the results mean. Although recent omic efforts have contributed to a better understanding of many facets of termite and symbiont biology, and have created important new resources for many species, significant knowledge gaps still remain. Crossing these gaps can best be done by applying new omic resources within multi-dimensional (i.e., functional, translational, and applied) research programs.

Gene expression and molecular phylogenetic analyses of beta-glucosidase in the termite Reticulitermes speratus (Isoptera: Rhinotermitidae)

Beta-glucosidase (BG) is known as a multifunctional enzyme for social maintenance in terms of both cellulose digestion and social communication in termites. However, the expression profiles of each BG gene and their evolutionary history are not well understood. First, we cloned two types of BG homologs (RsBGI and RsBGII) from the termite Reticulitermes speratus (Kolbe). Gene expression analyses showed that RsBGI expression levels of primary queens and kings from 30 to 100 days after colony foundation were high, but those of reproductives dropped after day 400. Extremely low gene expression levels of RsBGI were observed in eggs, whereas workers had significantly higher expression levels than those of soldiers and other colony members. Consequently, RsBGI gene expression levels changed among each developmental stage, and RsBGI was shown to be involved in cellulose digestion. On the other hand, the RsBGII gene was consistently expressed in all castes and developmental stages examined, and notable expression changes were not observed among them, including in eggs. It was indicated that RsBGII is a main component involved in social communication, for example, the egg-recognition pheromone shown in this species previously. Finally, we obtained partial gene homologs from other termite and cockroach species, including the woodroach (genus Cryptocercus), which is the sister group to termites, and performed molecular phylogenetic analyses. The results showed that the origin of the BG gene homologs preceded the divergence of termites and cockroaches, suggesting that the acquisition of multifunctionality of the BG gene also occurred in cockroach lineages.

A GHF7 cellulase from the protist symbiont community of Reticulitermes flavipes enables more efficient lignocellulose processing by host enzymes

Termites and their gut microbial symbionts efficiently degrade lignocellulose into fermentable monosaccharides. This study examined three glycosyl hydrolase family 7 (GHF7) cellulases from protist symbionts of the termite Reticulitermes flavipes. We tested the hypotheses that three GHF7 cellulases (GHF7-3, GHF7-5, and GHF7-6) can function synergistically with three host digestive enzymes and a fungal cellulase preparation. Full-length cDNA sequences of the three GHF7s were assembled and their protist origins confirmed through a combination of quantitative PCR and cellobiohydrolase (CBH) activity assays. Recombinant versions of the three GHF7s were generated using a baculovirus-insect expression system and their activity toward several model substrates compared with and without metallic cofactors. GHF7-3 was the most active of the three cellulases; it exhibited a combination of CBH, endoglucanase (EGase), and β-glucosidase activities that were optimal around pH 7 and 30°C, and enhanced by calcium chloride and zinc sulfate. Lignocellulose saccharification assays were then done using various combinations of the three GHF7s along with a host EGase (Cell-1), beta-glucosidase (β-glu), and laccase (LacA). GHF7-3 was the only GHF7 to enhance glucose release by Cell-1 and β-glu. Finally, GHF7-3, Cell-1, and β-glu were individually tested with a commercial fungal cellulase preparation in lignocellulose saccharification assays, but only β-glu appreciably enhanced glucose release. Our hypothesis that protist GHF7 cellulases are capable of synergistic interactions with host termite digestive enzymes is supported only in the case of GHF7-3. These findings suggest that not all protist cellulases will enhance saccharification by cocktails of other termite or fungal lignocellulases.

Immune-related transcriptome of Coptotermes formosanus Shiraki workers: the defense mechanism

Formosan subterranean termites, Coptotermes formosanus Shiraki, live socially in microbial-rich habitats. To understand the molecular mechanism by which termites combat pathogenic microbes, a full-length normalized cDNA library and four Suppression Subtractive Hybridization (SSH) libraries were constructed from termite workers infected with entomopathogenic fungi (Metarhizium anisopliae and Beauveria bassiana), Gram-positive Bacillus thuringiensis and Gram-negative Escherichia coli, and the libraries were analyzed. From the high quality normalized cDNA library, 439 immune-related sequences were identified. These sequences were categorized as pattern recognition receptors (47 sequences), signal modulators (52 sequences), signal transducers (137 sequences), effectors (39 sequences) and others (164 sequences). From the SSH libraries, 27, 17, 22 and 15 immune-related genes were identified from each SSH library treated with M. anisopliae, B. bassiana, B. thuringiensis and E. coli, respectively. When the normalized cDNA library was compared with the SSH libraries, 37 immune-related clusters were found in common; 56 clusters were identified in the SSH libraries, and 259 were identified in the normalized cDNA library. The immune-related gene expression pattern was further investigated using quantitative real time PCR (qPCR). Important immune-related genes were characterized, and their potential functions were discussed based on the integrated analysis of the results. We suggest that normalized cDNA and SSH libraries enable us to discover functional genes transcriptome. The results remarkably expand our knowledge about immune-inducible genes in C. formosanus Shiraki and enable the future development of novel control strategies for the management of Formosan subterranean termites.

Termite hindguts and the ecology of microbial communities in the sequencing age

Advances in high-throughput nucleic acid sequencing have improved our understanding of microbial communities in a number of ways. Deeper sequence coverage provides the means to assess diversity at the resolution necessary to recover ecological and biogeographic patterns, and at the same time single-cell genomics provides detailed information about the interactions between members of a microbial community. Given the vastness and complexity of microbial ecosystems, such analyses remain challenging for most environments, so greater insight can also be drawn from analysing less dynamic ecosystems. Here, we outline the advantages of one such environment, the wood-digesting hindgut communities of termites and cockroaches, and how it is a model to examine and compare both protist and bacterial communities. Beyond the analysis of diversity, our understanding of protist community ecology will depend on using statistically sound sampling regimes at biologically relevant scales, transitioning from discovery-based to experimental ecology, incorporating single-cell microbiology and other data sources, and continued development of analytical tools.

Evidence of a conserved functional role for DNA methylation in termites

Many organisms are capable of developing distinct phenotypes from the same genotype. This developmental plasticity is particularly prevalent in insects, which can produce alternate adaptive forms in response to distinct environmental cues. The ability to develop divergent phenotypes from the same genotype often relies on epigenetic information, which affects gene function and is transmitted through cell divisions. One of the most important epigenetic marks, DNA methylation, has been lost in several insect lineages, yet its taxonomic distribution and functional conservation remain uninvestigated in many taxa. In the present study, we demonstrate that the signature of high levels of DNA methylation exists in the expressed genes of two termites, Reticulitermes flavipes and Coptotermes formosanus. Further, we show that DNA methylation is preferentially targeted to genes with ubiquitous expression among morphs. Functional associations of DNA methylation are also similar to those observed in other invertebrate taxa with functional DNA methylation systems. Finally, we demonstrate an association between DNA methylation and the long-term evolutionary conservation of genes. Overall, our findings strongly suggest DNA methylation is present at particularly high levels in termites and may play similar roles to those found in other insects.

Construction and characterization of normalized cDNA libraries by 454 pyrosequencing and estimation of DNA methylation levels in three distantly related termite species

In termites, division of labor among castes, categories of individuals that perform specialized tasks, increases colony-level productivity and is the key to their ecological success. Although molecular studies on caste polymorphism have been performed in termites, we are far from a comprehensive understanding of the molecular basis of this phenomenon. To facilitate future molecular studies, we aimed to construct expressed sequence tag (EST) libraries covering wide ranges of gene repertoires in three representative termite species, Hodotermopsis sjostedti, Reticulitermes speratus and Nasutitermes takasagoensis. We generated normalized cDNA libraries from whole bodies, except for guts containing microbes, of almost all castes, sexes and developmental stages and sequenced them with the 454 GS FLX titanium system. We obtained >1.2 million quality-filtered reads yielding >400 million bases for each of the three species. Isotigs, which are analogous to individual transcripts, and singletons were produced by assembling the reads and annotated using public databases. Genes related to juvenile hormone, which plays crucial roles in caste differentiation of termites, were identified from the EST libraries by BLAST search. To explore the potential for DNA methylation, which plays an important role in caste differentiation of honeybees, tBLASTn searches for DNA methyltransferases (dnmt1, dnmt2 and dnmt3) and methyl-CpG binding domain (mbd) were performed against the EST libraries. All four of these genes were found in the H. sjostedti library, while all except dnmt3 were found in R. speratus and N. takasagoensis. The ratio of the observed to the expected CpG content (CpG O/E), which is a proxy for DNA methylation level, was calculated for the coding sequences predicted from the isotigs and singletons. In all of the three species, the majority of coding sequences showed depletion of CpG O/E (less than 1), and the distributions of CpG O/E were bimodal, suggesting the presence of DNA methylation.

Characterization of head transcriptome and analysis of gene expression involved in caste differentiation and aggression in Odontotermes formosanus (Shiraki)

BACKGROUND

The subterranean termite Odontotermes formosanus (Shiraki) is a serious insect pest of trees and dams in China. To date, very little is known about genomic or transcriptomic data for caste differentiation and aggression in O. formosanus. Hence, studies on transcriptome and gene expression profiling are helpful to better understand molecular basis underlying caste differentiation and aggressive behavior in O. formosanus.

METHODOLOGY AND PRINCIPAL FINDINGS

Using the Illumina sequencing, we obtained more than 57 million sequencing reads derived from the heads of O. formosanus. These reads were assembled into 116,885 unique sequences (mean size  =  536 bp). Of the unigenes, 30,646 (26.22%) had significant similarity with proteins in the NCBI nonredundant protein database and Swiss-Prot database (E-value<10(-5)). Of these annotated unigenes, 10,409 and 9,009 unigenes were assigned to gene ontology categories and clusters of orthologous groups, respectively. In total, 19,611 (25.52%) unigenes were mapped onto 242 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG). A total of 11,661 simple sequence repeats (SSRs) were predicted from the current transcriptome database. Moreover, we detected seven putative genes involved in caste differentiation and six putative genes involved in aggression. The qPCR analysis showed that there were significant differences in the expression levels of the three putative genes hexamerin 2, β-glycosidase and bicaudal D involved in caste differentiation and one putative gene Cyp6a20 involved in aggression among workers, soldiers and larvae of O. formosanus.

CONCLUSIONS

To our knowledge, this is the first study to characterize the complete head transcriptome of a higher fungus-cultivating termite using high-throughput sequencing. Our study has provided the comprehensive sequence resources available for elucidating molecular basis underlying caste differentiation and aggressive behavior in O. formosanus.

ASGARD: an open-access database of annotated transcriptomes for emerging model arthropod species

The increased throughput and decreased cost of next-generation sequencing (NGS) have shifted the bottleneck genomic research from sequencing to annotation, analysis and accessibility. This is particularly challenging for research communities working on organisms that lack the basic infrastructure of a sequenced genome, or an efficient way to utilize whatever sequence data may be available. Here we present a new database, the Assembled Searchable Giant Arthropod Read Database (ASGARD). This database is a repository and search engine for transcriptomic data from arthropods that are of high interest to multiple research communities but currently lack sequenced genomes. We demonstrate the functionality and utility of ASGARD using de novo assembled transcriptomes from the milkweed bug Oncopeltus fasciatus, the cricket Gryllus bimaculatus and the amphipod crustacean Parhyale hawaiensis. We have annotated these transcriptomes to assign putative orthology, coding region determination, protein domain identification and Gene Ontology (GO) term annotation to all possible assembly products. ASGARD allows users to search all assemblies by orthology annotation, GO term annotation or Basic Local Alignment Search Tool. User-friendly features of ASGARD include search term auto-completion suggestions based on database content, the ability to download assembly product sequences in FASTA format, direct links to NCBI data for predicted orthologs and graphical representation of the location of protein domains and matches to similar sequences from the NCBI non-redundant database. ASGARD will be a useful repository for transcriptome data from future NGS studies on these and other emerging model arthropods, regardless of sequencing platform, assembly or annotation status. This database thus provides easy, one-stop access to multi-species annotated transcriptome information. We anticipate that this database will be useful for members of multiple research communities, including developmental biology, physiology, evolutionary biology, ecology, comparative genomics and phylogenomics.

Database URL:asgard.rc.fas.harvard.edu