- Research article
- Open Access
Metatranscriptomic profiles of Eastern subterranean termites, Reticulitermes flavipes (Kollar) fed on second generation feedstocks
© Rajarapu et al.; licensee BioMed Central. 2015
- Received: 23 February 2015
- Accepted: 27 March 2015
- Published: 22 April 2015
Second generation lignocellulosic feedstocks are being considered as an alternative to first generation biofuels that are derived from grain starches and sugars. However, the current pre-treatment methods for second generation biofuel production are inefficient and expensive due to the recalcitrant nature of lignocellulose. In this study, we used the lower termite Reticulitermes flavipes (Kollar), as a model to identify potential pretreatment genes/enzymes specifically adapted for use against agricultural feedstocks.
Metatranscriptomic profiling was performed on worker termite guts after feeding on corn stover (CS), soybean residue (SR), or 98% pure cellulose (paper) to identify (i) microbial community, (ii) pathway level and (iii) gene-level responses. Microbial community profiles after CS and SR feeding were different from the paper feeding profile, and protist symbiont abundance decreased significantly in termites feeding on SR and CS relative to paper. Functional profiles after CS feeding were similar to paper and SR; whereas paper and SR showed different profiles. Amino acid and carbohydrate metabolism pathways were downregulated in termites feeding on SR relative to paper and CS. Gene expression analyses showed more significant down regulation of genes after SR feeding relative to paper and CS. Stereotypical lignocellulase genes/enzymes were not differentially expressed, but rather were among the most abundant/constitutively-expressed genes.
These results suggest that the effect of CS and SR feeding on termite gut lignocellulase composition is minimal and thus, the most abundantly expressed enzymes appear to encode the best candidate catalysts for use in saccharification of these and related second-generation feedstocks. Further, based on these findings we hypothesize that the most abundantly expressed lignocellulases, rather than those that are differentially expressed have the best potential as pretreatment enzymes for CS and SR feedstocks.
- Corn stover
- Soybean residue
Second generation biofuels have gained interest in the field of biofuel production due to their several advantages over the first generation biofuels. Second generation biofuels are produced from various kinds of plant biomass including agricultural/forest by-products, organic wastes and other dedicated feedstocks . Lignocellulose, a complex of lignin, cellulose and hemicellulose, is the major component of plant biomass and second generation feedstocks. Breakdown of the lignocellulose matrix is required for the release of free sugars for biofuel production, and is performed in the three major steps of pretreatment, hydrolysis and fermentation. Pretreatment is the limiting step in biofuel production and involves the disintegration of lignocellulose polymers to increase the availability of cellulose and hemicellulose for downstream biochemical processes . There are several physical and chemical pretreatment methods available (e.g., steam, liquid hot water, acid treatment) for lignocellulosic biomass which have been extensively reviewed [2-5]. However, the current pretreatment methods have disadvantages wherein the major concern is the inhibition of cellulolytic enzymes by lignin, mono-lignols and other phenolic metabolites in downstream processes [6,7]. To date, one of the most significant challenges to biofuel production is the development of economically efficient pretreatment processes.
Termites are efficient lignocellulose decomposers with the unique ability to utilize most of the cellulose from their diet [8-11]. The lower termite, Reticulitermes flavipes (Kollar), is used as a model in this study to investigate potential mechanisms of lignocellulose feedstock digestion. R. flavipes, along with its microbial symbionts, produces a bouquet of enzymes that participate in lignocellulose degradation [12,13]. Cellulases, hemicellulases and candidate ligninases were previously identified in R. flavipes [13-15]. A recombinant enzyme cocktail of R. flavipes enzymes was shown to result in the release of free sugars from wood substrate in vitro . This evidence supports the further investigation and application of R. flavipes digestion mechanism(s) in biofuel production from lignocellulose biomasses.
In this study we evaluated changes in the metatranscriptomic profiles of the R. flavipes gut microenvironment in response to feeding on two second generation feedstocks (corn stover [CS] and soybean hull residue [SR]). Our hypothesis was that genes encoding key lignocellulose processing enzymes would be differentially expressed among these treatments. To test this hypothesis, we studied the gut metatranscriptomes to identify potential lignocellulase genes of R. flavipes workers fed on different feedstocks via the Illumina Hi-Seq platform. Our findings show a clear ability of R. flavipes to feed and survive on two very different plant feedstocks (CS and SR); however, transcripts encoding wood-associated lignocellulases identified through previous work were not affected by the diets. Along with these results our Illumina deep sequencing results significantly contribute to the molecular resources available for R. flavipes, and also shed light on novel feedstock-associated genes from R. flavipes.
Feedstock composition analysis
Lignin content in the feedstocks was estimated using the acetyl bromide method of Fukushima and Kereley [17,18]. Briefly, cell walls of the feedstocks were extracted using the following method. Finely milled feedstocks (60 mm particle size) were incubated in sodium phosphate buffer (0.1 M, pH 7) for 30 minutes at 50°C followed by centrifugation at 3000 g for 10 min. Supernatant was discarded and the pellet was washed with 70% ethanol at 70°C for 5 to 6 times. One volume of acetone was added and centrifuged at 3000 g for 10 minutes. Supernatant was discarded and pellet was dried in oven at 45 ± 5°C for overnight. Approximately 5 mg of dried cell wall was dissolved in 2.5 ml of 25% acetyl bromide v/v in glacial acetic acid at 70°C for 30 min in a glass tube closed with a Teflon lined cap. Acetyl bromide with dissolved cell wall was then transferred to a 50 ml volumetric flask containing 12 ml glacial acetic acid and 2.5 ml of 2 N sodium hydroxide followed by the addition of 0.5 ml of 7.5 M hydroxylamine hydrochloride. The final volume was made to 50 ml and absorbance was measured at 280 nm in a spectrophotometer (BioRad SmartSpec™ 3000). Two feedstocks including CS and SR were run in triplicates; blanks were included to correct for background absorbance [17,18].
Carbohydrate analysis was performed following the National Renewable Energy Lab protocol with slight modifications . Briefly, 300 mg of feedstock including CS and SR was hydrolyzed in 3 ml of concentrated sulfuric acid in pressure tubes at 30°C for 1 hr with occasional stirring. Hydrolysis was stopped by adding 84 ml of water and further hydrolyzed at 121°C in an autoclave for 1 hour. The resultant liquid was aliquoted in 1 ml glass vials for quantification of carbohydrates using a Waters e2695 separation module HPLC system at the Analytical chemistry center, Laboratory of Renewable Resources Engineering, Purdue University.
Bioassays, gut removal and RNA isolation
Feeding assays were performed using R. flavipes workers collected from West Lafayette, Indiana (single colony) and Florida (2 colonies) which served as three independent biological replicates. Collected termites were maintained under laboratory conditions (22°C with 70% relative humidity and a 0:24 light:dark photoperiod) on wood shims and paper towels before the feeding assays. Diets used in this study were Whatman no.1 filter paper (Maidstone, UK), CS (Specialty 3557) and SR (Williams 82). The CS was donated by Dr. Nathan Mosier of Purdue University Agricultural and Biological Engineering department and SR was donated by Dr. Karen Hudson of Purdue University USDA-Agronomy department. Diets were ground into a fine powder using a DCG-20 coffee grinder (Cuisinart; Stamford, CT) and made into “cookies” using nanopure water. Cookies were then dried at 50 ± 5°C for 48 hr and weighed before feeding. Feeding assays were performed in 35 mm Petri dishes with 30 termites per diet treatment. Three independent replicates were performed per treatment. After 7 days, the weight of termites and any remaining diet were recorded after the assay period to account for diet consumption. Whole termite guts including symbionts were then removed and pooled in RNA isolation buffer, and total RNA extracted with the SV Total RNA Isolation System (Promega; Madison, WI). Total extracted RNA (>10 μg) was assessed for quality using a Nanodrop spectrophotometer (Thermo Scientific) and submitted to the Purdue University Genomics Core Facility (PGCF) for cDNA library synthesis and sequencing.
Illumina sequencing and bioinformatics analysis
Two microliters of the total RNA was further analyzed for quality using a Bioanalyzer (Agilent Inc.). One microgram each of total RNA was enriched for mRNA using polyT hybridization and cDNA libraries were bar-coded for all the 9 replicate samples by the PGCF using TruSeq™ RNA sample preparation kit (Illumina). The bar-coded libraries were paired end sequenced with a read length of 150 per read using the Illumina Hiseq2500 platform. Adapters were removed using Trimmomatic and pre-processed for downstream analysis. Low quality bases (Phred33 score of < 30) were trimmed using FASTX toolkit (v 0.0.13)  after accessing the quality using FastQC (v 0.10.0) . Annotation of raw reads was performed in a hierarchy for maximum annotation. The first annotation method used was MetaCV (v 2.2.9) analysis using a database of protist and termite protein sequences from NCBI . Previous cDNA sequencing of R. flavipes gut tissue by Tartar et al.  has identified a significant number of bacterial sequences in a symbiont library and thus the remaining reads were processed with MetaCV using a second database of all bacterial protein sequences extracted from NCBI. The non-annotated reads were joined using FastqJoin (v 1.1.2-537) and in-house perl scripts. The entire set of joined Fastq reads was further analyzed using BLAST against termite and protist database (v 2.2.28) to further classify the joined reads that failed to annotate during each round of MetaCV . The remaining joined reads that could not be annotated from BLAST were used as input for RAPSearch2 (v 2.12) . Each annotation step resulted in GIs (GenInfo Identifiers), which were used to extract genus and species names using the BLAST tool blastdbcmd. An abundance matrix containing counts was generated for all unique genera. Differential analysis of genera abundance for pairwise comparisons of all 3 treatments were conducted using R (version 3.0.3) and the edgeR (v 3.4.2) package with default parameters [26,27].
In a separate pipeline, total reads from all the libraries were joined via FastqJoin (v 1.1.2-537) and were uploaded into MG-RAST (Metagenomics – Rapid Annotation using Subsystem Technology ) through a private repository. The data were compared to M5NR database using a maximum e-value of 1e-5, a minimum identity of 60%, and a minimum alignment length of 15 measured in amino acids for protein and base pairs for RNA databases. The functionally annotated reads obtained from MG-RAST were used to generate a pathway and KO (KEGG orthology) term count matrix to evaluate differential functional expression using edgeR (v. 3.4.2). A Principal Co-ordinate Analysis (PCoA) was performed for genera and functional abundance profiles using Bray-Curtis distance measures for genera abundance and Euclidean distance measures for functional abundance using MG-RAST tool. Bacterial abundances were analyzed using best hit classification method in MG-RAST. Analyses of significant differences between bacterial phyla were carried out using two non-parametric tests ANOSIM  and Adonis . Bacterial abundance values were converted to a Bray-Curtis distance matrix using the ‘vegdist’ function from vegan package (Version 2.0.10) (R version 3.1.1). This distance matrix was used as an input for ANOSIM (Analysis of Similarity) and Adonis test. Proc GLM (General Linear Model) Analysis of variance (ANOVA) with “feedstocks” as the main effect was used to analyze the differences between bacterial abundance among the diets (Paper, CS and SR) in SAS (version 9.2). Data were transformed when required to follow normal distribution.
A post-hoc feeding assay was set up with 10 termites per treatment per colony as described earlier, using five different lab-maintained termite colonies which were in the lab for a similar time period as the sequenced colonies. Protists were counted from the hindguts of termites using a Sedgewick rafter cell counter. Hindguts of termites were collected in 1X phosphate buffered saline (100 μl/gut). The resultant stock solution was diluted appropriately to 1 ml to achieve a countable dilution of cells per square of the cell counter. Cells were counted using a phase contrast microscope under the 20X objective. The final calculation for cell count in the stock solution was carried out using the following formula: cells/ml = (dilution factor) (cells counted) (calibration constant)/number of squares containing cells counted . Mean counts were compared statistically using non-parametric Kruskal-Wallis tests followed by a Tukey’s test for pairwise comparisons (Sigma plot v11.0).
Differentially expressed genes along with the abundantly expressed putative lignocellulases including the Cell −1 (endoglucanase), GHF-7 (Cellobiohydrolase) were validated using RT-qPCR. Primers were designed using Primer3  (Additional file 1: Table S1). Feeding bioassays were repeated with the same experimental design as mentioned earlier and three different R. flavipes colonies collected from West Lafayette, Indiana were used, which were also assessed for protist counts. Validations were performed on the sequenced samples and also the repeated experimental replicates (i.e., new biological samples). Total RNA isolated from each experimental replicate was used to synthesize first strand cDNA using the iScript cDNA synthesis kit (BioRad). RT-qPCR reactions were performed in a 10 μl reaction mixture containing 1 μl of cDNA template, 0.5 μl of each forward and reverse primers (10 μM), 5 μl SYBR green mastermix (2X, BioRad) and 3 μl of nuclease free water. All the primers designed for validation were tested a priori for specificity using PCR and for efficiency using standard curves of diluted template R. flavipes specific translation elongation factor delta (TEF-δ) was tested for stability across the treatments and used as the reference gene in this study. Thermocycler conditions for RT-qPCR were: 95°C for 3 min, 40 cycles of 95°C for 30 sec, 60°C for 30 sec, and 72°C for 30 sec. Relative fold changes were calculated using delta Ct method . Fold changes from RNA-Seq and RT-qPCR are compared using a paired T-test (Sigma plot v11.0). A Spearmann correlation was performed for the fold changes of RNA-seq and RT-qPCR. Correlation was also tested for RT-qPCR fold changes from Illumina sequenced samples and experimentally replicated samples (i.e., new biological samples).
Illumina sequencing resulted in approximately 800 million reads (400 million paired-end reads) in total from all treatments and replicates (Additional file 1: Table S3). An informatics pipeline involving several steps and analysis platforms was used for analysis (Additional file 2: Figure S1). Out of the total reads obtained from all the treatments, 98% of raw reads were retained after the removal of low scored reads (phred33 < 30). MetaCV searches of total reads against termite and protist databases assigned taxa identities to 17% of all reads with an e-value of 10−3. Annotation with MetaCV using a bacterial database annotated 2% of the total reads. Further, joining of overlapping reads that were initially non-annotated made it possible to annotate an additional 1-2% of the total reads. Specifically, BLAST searches of these joined reads against termite and protist EST databases annotated 38% of these reads with an e-value of 10−3. A RAPSearch2 identified 7% of the total joined reads with an e-value of 10−3 against a non-redundant database (Additional file 1: Table S4A).
Analysis of joined reads obtained per treatment in MG-RAST removed an average of 46% of the joined reads due to quality control screening across the treatments. Another 46% of reads per treatment library that passed quality control were putatively annotated as containing ORFs. Nine percent of the predicted ORFs per library had homology to proteins and were assigned putative protein IDs. The final step in the MG-RAST pipeline assigned 5% of total reads on average per treatment to pathways. In total, approximately 2.5 M reads per replicate per feedstock treatment were functionally annotated (Additional file 1: Table S4B). Rarefaction curves show the intensive sampling of the nine metatranscriptome libraries (Additional file 2: Figure S2). Also, a majority of the annotated sequences across the treatment libraries had similarity to phylum Arthropoda and other eukaryotes (Additional file 2: Figure S3).
Community and functional profiles
Community level differentials
uncultured protist Unknown
Functional level differentials
Cysteine and methionine metabolism
D-Glutamine and D-glutamate metabolism
Glycine, serine and threonine metabolism
Valine, leucine and isoleucine biosynthesis
Fructose and mannose metabolism
Hematopoietic cell lineage
HIF-1 signaling pathway
Pathogenic Escherichia coli infection
Proximal tubule bicarbonate reclamation
Vasopressin-regulated water reabsorption
Secondary metabolite production
Cofactor and Vitamin metabolism
Meiosis - yeast
Here we compared gut metatranscriptomes of R. flavipes workers feeding on different second generation feedstocks. Our goal was to identify the lignocellulase genes responding to these diets for their potential use in the pretreatment, and potentially saccharification processes of biofuel production. While our results suggest that CS and SR diets have little impact on termite digestome composition, our findings also suggest that the necessary lignocellulases required for digestion are constitutively expressed. However, we also found that the microbial gut community profiles of termites fed CS and SR diets shifted relative to those fed a paper diet. The CS metatranscriptome was functionally similar to that of the paper diet; whereas the SR metatranscriptome differed from the paper diet. It is notable that SR-fed termites were affected at the microbial community, pathway and functional levels. Protein and carbohydrate pathways were down-regulated in termites fed on SR; whereas, CS-fed termites exhibited minimal effects at the pathway and gene levels.
Interestingly, in the seven-day feeding bioassay, R. flavipes workers consumed more CS and SR than pure-cellulose paper (Figure 1A). Relative to paper, CS and SR have lesser carbohydrate and higher lignin content which might have resulted in higher consumption by termites to sustain themselves on these diets. The metatranscriptome profiles of both CS and SR-fed termites differed at the symbiont community level. A significant decrease in protist abundance in SR-fed termites relative to paper and CS-fed termites might be the effect of secondary compounds in SR. Isoflavones are the major secondary chemicals in SR, and in CS it is DIMBOA (hydroxamic acid, 2,4-dihydroxy-7-methoxy-1,4-benzoxazin-3-one). Several secondary plant metabolites are known to possess antimicrobial and anti-protozoal properties . A study in ruminants has similarly shown a decrease in protozoan abundance upon feeding on certain plant materials . A single study with the lower termite, Kalotermes flavicollis, showed a decrease in flagellate protists upon feeding on plant based essential oils . However, there is no direct evidence suggesting antiprotozoal activity of soybean and corn extracts. Among the bacterial sequences obtained in this study there were few groups of bacteria that were significantly affected by the tested feedstocks. The effects of diet on the bacterial gut community cannot be reliably concluded as relatively few bacterial sequences with a potential polyA tail might have been sequenced. However, earlier findings by Boucias et al.  showed that feeding different lignocellulose diets for the same amount of time as the present study had minimal effects on the termite gut bacterial symbiont community .
At the pathway level (Table 2), the downregulation of protein and carbohydrate metabolism pathways in SR-fed termites might have been related to decreases in protist abundance. Several protists and their prokaryotic symbionts are well known to participate in protein and carbohydrate metabolism in termites [11,38]. A decline in the number of protists might have affected protein and carbohydrate metabolism in termites fed the SR diet. Pathway-level expression patterns were further supported by gene or function-level expression patterns. Genes participating in protein, carbohydrate and energy metabolism were downregulated in SR-fed termites relative to paper. Protein and carbohydrate metabolism pathways along with the genes participating in these pathways were more highly expressed in CS-fed termites relative to SR-fed termites. This further suggests that CS has less of a deleterious effect on the digestion mechanisms of R. flavipes.
There was one up-regulated gene in each treatment; 5-carboxymethyl-2-hydroxymuconate isomerase was up-regulated in CS-fed termites, which acts on ring-opened metabolites of lignin phenolics  and participates in phenyl acetate catabolism for carbon in E. coli [40,41]. Upregulation of this gene either indicates the requirement of carbon by gut bacteria or a mechanism to cope with excess phenolic substrates encountered in the gut. MAPK41, upregulated in SR-fed termites relative to paper-fed, belongs to the superfamily of mitogen-activated protein kinases which are vital players in signal transduction. MAPK41 along with other proteins regulate several cellular activities including transcription, cell division, metabolism, movement and apoptosis . Condensin complex subunit was also more highly expressed in SR-fed termites compared to CS-fed termites. Condensin complex subunit is important for the organization of chromosomes during cell division and is ubiquitously distributed in archaea, bacteria and eukaryotes . The higher expression of MAPK41 and Condensin complex subunit is consistent with the possibility that there is a high turnover of cells in SR-fed termites.
Among the KO terms that had higher transcript levels in CS-fed relative to SR-fed termite guts, ~50% are terms related to amino acid metabolism and malate dehydrogenase isoforms. Higher levels of amino acid metabolism terms in CS-fed termites relative to SR-fed further supports the detrimental effects of SR compared to CS on termites and symbionts. Malate dehydrogenases are ubiquitously distributed and participate in a variety of metabolic pathways . In prokaryotes, malate dehydrogenase plays a role in carbon fixation . Higher expression of malate dehydrogenases in CS-fed termites compared to SR-fed suggests higher metabolism in termites or their symbionts with CS feeding.
Expression patterns similar to RNA-seq were observed for the tested KO term transcripts using RT-qPCR. However, EdgeR predicted differentials were not statistically significant in RT-qPCR, which is likely a result of differences in the normalization for the two methods. Further, the lack of correlation observed is due to differences in the magnitude of fold changes between the two techniques. However, a significant correlation was observed when using RT-qPCR to compare Illumina sequenced samples with independent biological samples (Additional file 2: Figure S6), which along with other results, provides reasonable validation of RNA-Seq results.
In conclusion, this study provides a novel glimpse, using termites as a model system, into possible mechanisms of lignocellulosic feedstock digestion. The sequences obtained from this study also tremendously enrich existing R. flavipes sequence databases [16,23,45]. Our findings here were unexpected in that they did not support our initial hypothesis (i.e., that genes encoding stereotypical digestive and other biomass-active enzymes would be induced by feeding termites agricultural feedstocks). The metatranscriptomic profiles of CS and SR-fed termites were more similar than expected to those of paper-fed termites in terms of lignocellulase gene expression. Also, although SR effects were more severe, both feedstocks had negative impacts on protist abundance, as well as probable functions associated with protists. However, the CS and SR diets had no effect on the expression levels of stereotypical lignocellulases which were constitutively expressed. For example, a host endoglucanase identified through previous work, Cell-1 , was constitutively expressed at high levels irrespective of the feedstock fed to termites. Prior in vitro analyses demonstrated the central role of Cell-1 in efficiently releasing glucose from wood substrates independent of symbiont action . Also, termites are wood-feeding specialists and the genes that are used for wood digestion might be genetically determined to express constitutively and may therefore be non-responsive to CS and SR regulation. Thus, we hypothesize that constitutively expressed lignocellulases (Cell-1, aldo-keto reductases, laccases, catalases, etc. ) may serve as biocatalysts for both pretreatment and saccharification of CS and SR feedstocks. This hypothesis is being further tested at the translational level using proteomics, along with other integrative approaches, and will be the focus of a forthcoming report.
We thank Amit Sethi for conducting feeding assays and sample preparation, the Purdue University Genomics Core Facility for sequencing and first-pass informatics, Jesse Hoteling for assistance with RT-qPCR; Iman Beheshti-Tabar for assistance with carbohydrate analysis; Brittany Peterson and Ameya Gondhalekar for helpful manuscript comments. This work was supported by National Science Foundation grant No. CBET 1233484 awarded to MES, and the O.W. Rollins/Orkin Endowment in the Department of Entomology at Purdue University.
- Naik SN, Goud VV, Rout PK, Dalai AK. Production of first and second generation biofuels: a comprehensive review. Renew Sust Energ Rev. 2010;14(2):578–97.View ArticleGoogle Scholar
- Mosier N, Wyman C, Dale B, Elander R, Lee YY, Holtzapple M, et al. Features of promising technologies for pretreatment of lignocellulosic biomass. Bioresour Technol. 2005;96(6):673–86.View ArticlePubMedGoogle Scholar
- Hendriks ATWM, Zeeman G. Pretreatments to enhance the digestibility of lignocellulosic biomass. Bioresour Technol. 2009;100(1):10–8.View ArticlePubMedGoogle Scholar
- Wyman CE, Dale BE, Elander RT, Holtzapple M, Ladisch MR, Lee Y. Coordinated development of leading biomass pretreatment technologies. Bioresour Technol. 2005;96(18):1959–66.View ArticlePubMedGoogle Scholar
- Kumar P, Barrett DM, Delwiche MJ, Stroeve P. Methods for pretreatment of lignocellulosic biomass for efficient hydrolysis and biofuel production. Ind Eng Chem Res. 2009;48(8):3713–29.View ArticleGoogle Scholar
- Nakagame S, Chandra RP, Kadla JF, Saddler JN. The isolation, characterization and effect of lignin isolated from steam pretreated Douglas-fir on the enzymatic hydrolysis of cellulose. Bioresour Technol. 2011;102(6):4507–17.View ArticlePubMedGoogle Scholar
- Li H, Pu Y, Kumar R, Ragauskas AJ, Wyman CE. Investigation of lignin deposition on cellulose during hydrothermal pretreatment, its effect on cellulose hydrolysis, and underlying mechanisms. Biotechnol Bioeng. 2014;111(3):485–92.View ArticlePubMedGoogle Scholar
- Geib SM, Filley TR, Hatcher PG, Hoover K, Carlson JE, Jimenez-Gasco MM, et al. Lignin degradation in wood-feeding insects. Proc Natl Acad Sci U S A. 2008;105(35):12932–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Scharf ME, Tartar A. Termite digestomes as sources for novel lignocellulases. Biofuels Bioproducts & Biorefining-Biofpr. 2008;2(6):540–52.View ArticleGoogle Scholar
- Ohkuma M. Termite symbiotic systems: efficient bio-recycling of lignocellulose. Appl Microbiol Biotechnol. 2003;61(1):1–9.View ArticlePubMedGoogle Scholar
- Ohkuma M. Symbioses of flagellates and prokaryotes in the gut of lower termites. Trends Microbiol. 2008;16(7):345–52.View ArticlePubMedGoogle Scholar
- Ni JF, Tokuda G. Lignocellulose-degrading enzymes from termites and their symbiotic microbiota. Biotechnol Adv. 2013;31(6):838–50.View ArticlePubMedGoogle Scholar
- Scharf ME, Karl ZJ, Sethi A, Boucias DG. Multiple levels of synergistic collaboration in Termite lignocellulose digestion. PLoS ONE. 2011;6(7):e21709.View ArticlePubMed CentralPubMedGoogle Scholar
- Scharf ME, Kovaleva ES, Jadhao S, Campbell JH, Buchman GW, Boucias DG. Functional and translational analyses of a beta-glucosidase gene (glycosyl hydrolase family 1) isolated from the gut of the lower termite Reticulitermes flavipes. Insect Biochem Mol Biol. 2010;40(8):611–20.View ArticlePubMedGoogle Scholar
- Sethi A, Kovaleva ES, Slack JM, Brown S, Buchman GW, Scharf ME. A GHF7 Cellulase from the protist symbiont community of Reticulitermes flavipes enables more efficeint lignocellulose processing by host enzymes. Arch Insect Biochem Physiol. 2013;84(4):175–93.View ArticlePubMedGoogle Scholar
- Sethi A, Slack JM, Kovaleva ES, Buchman GW, Scharf ME. Lignin-associated metagene expression in a lignocellulose-digesting termite. Insect Biochem Mol Biol. 2013;43(1):91–101.View ArticlePubMedGoogle Scholar
- Fukushima RS, Hatfield RD. Extraction and isolation of lignin for utilization as a standard to determine lignin concentration using the acetyl bromide spectrophotometric method. J Agric Food Chem. 2001;49(7):3133–9.View ArticlePubMedGoogle Scholar
- Fukushima RS, Hatfield RD. Comparison of the acetyl bromide spectrophotometric method with other analytical lignin methods for determining lignin concentration in forage samples. J Agric Food Chem. 2004;52(12):3713–20.View ArticlePubMedGoogle Scholar
- Sluiter A, Hames B, Ruiz R, Scarlata C, Sluiter J, Templeton D, et al. NREL Biomass Program: Determination of structural carbohydrates and lignin in biomass. In: Biomass Analysis Technology Team, Laboratory AnalyticalProcedure. 2006.Google Scholar
- FastxToolkit available at http://hannonlab.cshl.edu/fastx_toolkit/index.html.
- A S. FASTQC : A quality control tool for high throughput sequence data. 2010.Google Scholar
- Liu J, Wang H, Yang H, Zhang Y, Wang J, Zhao F, et al. Composition-based classification of short metagenomic sequences elucidates the landscapes of taxonomic and functional enrichment of microorganisms. Nucleic Acids Res. 2012;41:e3.View ArticlePubMed CentralPubMedGoogle Scholar
- Tartar A, Wheeler MM, Zhou X, Coy MR, Boucias DG, Scharf ME. Parallel metatranscriptome analyses of host and symbiont gene expression in the gut of the termite Reticulitermes flavipes. Biotechnol Biofuels. 2009;2(1):25.View ArticlePubMed CentralPubMedGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.View ArticlePubMedGoogle Scholar
- Zhao Y, Tang H, Ye Y. RAPSearch2: a fast and memory-efficient protein similarity search tool for next-generation sequencing data. Bioinformatics. 2012;28(1):125–6.View ArticlePubMed CentralPubMedGoogle Scholar
- R Development Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria; 2013. [http://www.R-project.org/ webcite] [ISBN 3-900051-07-0]
- Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.View ArticlePubMed CentralPubMedGoogle Scholar
- Meyer F, Paarmann D, D'Souza M, Olson R, Glass EM, Kubal M, et al. The metagenomics RAST server–a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC bioinf. 2008;9(1):386.View ArticleGoogle Scholar
- CLARKE KR. Non‐parametric multivariate analyses of changes in community structure. Aust J Ecol. 1993;18(1):117–43.View ArticleGoogle Scholar
- Anderson MJ. A new method for non‐parametric multivariate analysis of variance. Austral ecol. 2001;26(1):32–46.Google Scholar
- Cell counting with the Sedgewick-Rafter chamber and Whipple Micrometer Disc [http://www.protocol-online.org/prot/Protocols/Cell-Counting-with-the-Sedgewick-Rafter-Chamber-and-Whipple-Micrometer-Disc-4315.html]
- Rozen S, Skaletsky H. Primer3 on the WWW for general users and for biologist programmers. In: Bioinformatics methods and protocols. Springer; 1999: 365–386Google Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using Real-Time Quantitative PCR and the 2−ΔΔCT method. Methods. 2001;25(4):402–8.View ArticlePubMedGoogle Scholar
- Wallace RJ. Antimicrobial properties of plant secondary metabolites. Proc Nutr Soc. 2004;63(4):621–9.View ArticlePubMedGoogle Scholar
- Santra A, Banerjee A, Das SK, Chatterjee A. Effect of plants containing secondary metabolites on ruminal fermentation and methanogenesis in vitro. Indian J Anim Sci. 2012;82(2):194–9.Google Scholar
- Alfazairy AAM. Antimicrobial activity of certain essential oils against hindgut symbionts of the dry wood termite Kalotermes flavicollis Fabr. and prevalent fungi on termite-infested wood. Egypt J Biol Pest Control. 2005;15(1/2):27–33.Google Scholar
- Boucias DG, Cai Y, Sun Y, Lietze V-U, Sen R, Raychoudhury R, et al. The hindgut lumen prokaryotic microbiota of the termite Reticulitermes flavipes and its responses to dietary lignocellulose composition. Mol Ecol. 2013;22(7):1836–53.View ArticlePubMedGoogle Scholar
- Ohkuma M, Brune A. Diversity, structure, and evolution of the termite gut microbial community. In: Bignell DE, Roisin Y, Lo N, editors. Biology of termites: a modern synthesis. Dordrecht, The Netherlands: Springer; 2011. p. 413–38.Google Scholar
- Lucey KS, Leadbetter JR. Catechol 2, 3-dioxygenase and other meta-cleavage catabolic pathway genes in the ‘anaerobic’termite gut spirochete Treponema primitia. Mol Ecol. 2014;23(6):1531–43.View ArticlePubMedGoogle Scholar
- Johnson Jr WH, Hajipour G, Whitman CP. Stereochemical Studies of 5-(Carboxymethyl)-2-hydroxymuconate Isomerase and 5-(Carboxymethyl)-2-oxo-3-hexene-1, 6-dioate Decarboxylase from Escherichia coli C: Mechanistic and Evolutionary Implications. J Am Chem Soc. 1995;117(34):8719–26.View ArticleGoogle Scholar
- Hajipour G, Johnson Jr WH, Dauben PD, Stolowich NJ, Whitman CP. Chemical and enzymic ketonization of 5-(carboxymethyl)-2-hydroxymuconate. J Am Chem Soc. 1993;115(9):3533–42.View ArticleGoogle Scholar
- Johnson GL, Lapadat R. Mitogen-Activated Protein Kinase pathways mediated by ERK, JNK, and p38 Protein Kinases. Science. 2002;298(5600):1911–2.View ArticlePubMedGoogle Scholar
- Hirano T. Condensins: universal organizers of chromosomes with diverse functions. Genes Dev. 2012;26(15):1659–78.View ArticlePubMed CentralPubMedGoogle Scholar
- Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.View ArticlePubMed CentralPubMedGoogle Scholar
- Raychoudhury R, Sen R, Cai Y, Sun Y, Lietze VU, Boucias DG, et al. Comparative metatranscriptomic signatures of wood and paper feeding in the gut of the termite Reticulitermes flavipes (Isoptera: Rhinotermitidae). Insect Mol Biol. 2013;22(2):155–71.View ArticlePubMedGoogle Scholar
- Zhou X, Smith JA, Oi FM, Koehler PG, Bennett GW, Scharf ME. Correlation of cellulase gene expression and cellulolytic activity throughout the gut of the termite Reticulitermes flavipes. Gene. 2007;395(1–2):29–39.View ArticlePubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.