BACKGROUND
class hexapoda is one of the most ancient and diverse animal lineages on the planet  <cit> , containing organisms capable of occupying many recalcitrant niches and persisting under intense environmental conditions including extreme temperatures, periods of desiccation, and exposure to toxins  <cit> .
many of its members are capable of thriving on suboptimal, nutritionally-deficient substrates  <cit> , including wood-boring beetles belonging to family cerambycidae that feed exclusively on woody tissue.
specifically, beetles in the genus anoplophora are wood-boring insects of great interest because some of its members preferentially target healthy host trees and have relatively broad host ranges  <cit> .
for example, the asian longhorned beetle  was introduced from china into the united states, canada, and several countries in europe and has been documented to complete development in approximately  <dig> deciduous tree species worldwide, including several genera commonly planted as feedstock   <cit> .
acer spp.
are the predominant hosts in the introduced range  <cit> .
this beetle poses a significant threat to urban streetscapes, has the potential to destroy up to 35% of the urban tree canopy in its introduced range, and has already caused millions of dollars in damage to urban landscapes  <cit> .
wood-borers, like a. glabripennis, are especially challenging to control in both their natural and invasive ranges because the larvae spend 1–2 years living deep inside their host trees  <cit> .
natural enemies are rare and treatment of host trees with systemic insecticides is costly and has variable efficacy against a. glabripennis larvae  <cit> .
the most effective method for eradication is destruction of infested and nearby host trees and implementation of strict quarantine measures to contain the infestation.
therefore, understanding the digestive physiology of this cerambycid at the genetic level is paramount to devising novel control strategies.
because a. glabripennis spends the majority of its lifecycle in the larval stage and feeds primarily in the heartwood of a broad range of healthy deciduous trees, it must overcome challenges of digesting intractable woody tissue in order to acquire sufficient nutrients to complete development  <cit> .
glucose is a predominant wood sugar, but it is present in the form of complex polysaccharides, including cellulose, hemicellulose, callose, and pectin, which are inherently difficult to digest and require a complex of hydrolytic enzymes for efficient degradation and liberation of sugar monomers  <cit> .
extensive hydrogen bonding coupled with linear configurations increases the crystallinity of these cell wall polysaccharides and decreases their permeability, further hindering the activity of hydrolytic enzymes.
plant cell wall polysaccharides are further protected from hydrolytic enzymes by lignin, a biopolymer containing over  <dig> types of chemical bonds that is extensively cross-linked to both cellulose and hemicellulose  <cit> , shielding them from digestion.
due to the random, heterogeneous nature of these cross-linkages and the high resilience of carbon-carbon and β-aryl ether linkages that dominate this macromolecule, lignin polymers can only be efficiently degraded through oxidative depolymerization, a process that has only been conclusively documented to be catalyzed by enzymes produced by a small number of wood degrading fungi  <cit> .
nitrogen is also extremely limited in woody tissues  <cit>  and plant cell wall proteins are intricately cross-linked with lignin and cellulose, making them difficult to access  <cit> .
other essential nutrients, including fatty acids, sterols, and vitamins are present in low concentrations or are completely absent  <cit> .
lastly, wood-feeding insects must overcome plant secondary metabolites that often accumulate to high concentrations in the heartwood through detoxification or sequestration processes  <cit> .
many wood-feeding beetles cultivate extracellular symbiotic fungi to facilitate digestion of woody tissue and nutrient acquisition, which are carried in mycangia or other specialized structures on their body  <cit> .
for example, bark beetles utilize a mass attack strategy, in which a mycangial fungus is directly inoculated into a host tree during oviposition to facilitate pre-digestion of woody tissue and mitigation of host tree defenses.
an alternative strategy is to preferentially colonize stressed trees  <cit>  whose woody components have already been pre-digested by wood-rotting microbes.
however, a. glabripennis is distinct from many other wood-feeding beetles in the sense that a single larvae can successfully develop in a healthy tree without requiring mass attack and the majority of the challenging reactions, including digestion of lignocellulose and hemicellulose and detoxification of plant metabolites, can occur within the gut itself .
while the midgut community associated with a. glabripennis has the metabolic potential to overcome many of the challenges associated with feeding in woody tissue, including degradation of lignin, cellulose, and hemicellulose and acquisition of nitrogen and other essential nutrients  <cit> , the contributions of insect-derived digestive and nutrient acquiring enzymes cannot be ignored since insects themselves can produce a diverse array of digestive enzymes, including cellulases, hemicellulases, pectinases, and enzymes that enhance lignin degradation .
insects have also evolved other sophisticated abilities to evade host plant defenses and often possess extensive suites of enzymes involved in detoxification of plant metabolites and phytohormones, digestive proteinase inhibitors  <cit> , and cyanates and cyanoamino acids  <cit> , as well as enzymes capable of disrupting jasmonic acid signaling pathways  <cit> .
furthermore, insects produce many cytochrome p450s  <cit> , which are integrally involved in xenobiotic metabolic processes that ultimately lead to oxidative destruction of toxic compounds, including plant derived secondary metabolites and pesticides  <cit> .
the primary goals of this study were to survey the endogenous digestive and physiological capabilities of larval a. glabripennis through shotgun sequencing of midgut derived messenger rna and to identify insect-derived genes that are highly expressed in the midgut while actively feeding in wood.
the a. glabripennis midgut transcriptome library was also compared to all publically available transcriptome libraries sampled from other plant feeding insects to identify core groups of genes that are associated with digestive processes that could facilitate nutrient recovery from woody tissue regardless of insect taxa.
this study represents an important addendum to the growing database of genomic and transcriptomic resources available for coleopterans and fills an important gap, representing the first transcriptome sampled from a wood-feeding cerambycid and the first comprehensive analysis of endogenous genes associated with wood-feeding in insects.
these findings offer unique opportunities to bioprospect for enzymes that could be exploited for cellulosic biofuel production or other industrial processes, and to develop novel control methods for this destructive wood-boring pest and other wood-feeding insects.
RESULTS
454- and illumina-based transcriptome sequencing
to develop a comprehensive profile of the endogenous digestive and physiological capabilities of a. glabripennis, mrna was collected from the midguts of third instar larvae feeding in the heartwood of a preferred host  and was sequenced using both roche  <dig> pyrosequencing and illumina technologies.
in total,  <dig>  shotgun sequence reads were produced using the roche  <dig> flx platform using two separate runs.
<dig>  reads , ranging in length from  <dig> to  <dig> nt , were generated on a half plate and  <dig>  reads  ranging from  <dig> to  <dig> nt , were generated on a quarter plate.
these runs correspond to e4gebh <dig> sff and e5ty7pb <dig> sff from sra , respectively.
reads from both runs were pooled and were quality filtered and assembled together.
approximately  <dig>   of the total  <dig> flx reads passed quality filtering and were utilized in the assembly.
to enhance sequencing depth and acquire a more complete inventory of the endogenous digestive and metabolic capabilities of a. glabripennis,  <dig> million 100 bp paired illumina reads  with a library insert size of  <dig> nucleotides  were generated on a single lane using the illumina hiseq  <dig> .
after quality filtering and adapter removal, over  <dig> million read pairs  remained and were utilized in downstream processing and analyses.
digital k-mer normalization reduced the number of illumina read pairs to  <dig> , <dig>  which were ultimately used for co-assembly with the  <dig> flx reads.
assembly and annotation statistics
<dig> assembly and annotation statistics for comparative transcriptomics
to facilitate comparisons to transcriptome libraries prepared from the guts of other herbivorous insects, which were derived solely from  <dig> reads, the  <dig> reads were first assembled and analyzed without the illumina reads.
of the  <dig>  shotgun reads generated through  <dig> pyrosequencing , approximately  <dig>  reads assembled into  <dig>  contigs , ranging in length from  <dig> nt to  <dig>  nt with an n <dig> contig length of  <dig> nt .
assembled contigs that shared common reads were placed into isogroups.
these contigs are often broken at branch points between exon boundaries in multiple transcript isoforms from the same unigene.
contig branch structures within each isogroup were then traversed to create  <dig>  isotigs , which represent unique assembled transcripts or transcript fragments.
the n <dig> isotig length was  <dig>  nt and isotigs were grouped into  <dig>  isogroups, representing a gene locus or unigene.
of these isogroups,  <dig>  were comprised of a single transcript isoform and the average number of isotigs within an isogroup was  <dig> .
the maximum number of isotigs classified to the same isogroup was  <dig>  for downstream comparative analyses, isogroups were treated as unigenes and isotigs associated with the same isogroup were treated as transcript isoforms.
roughly  <dig>  reads  were singletons and were not included in the assembly.
of the singletons, approximately  <dig>  reads  were flagged as high quality and, to increase the amount of information present in the transcriptome dataset, these singleton reads were concatenated to the assembly and the pooled dataset was utilized in downstream transcriptome comparisons.
assembly metrics from the 454-based assembly are presented in table  <dig>  after clustering the isotigs and high quality singletons with cd-hit-est using a sequence similarity threshold of  <dig>  to group transcripts that likely represented allelic variants of the same gene, the total number of isotigs and singletons was reduced to around  <dig> .
seventy-eight of these isotigs and reads were classified as ribosomal rnas, while none were classified as trnas.
roughly  <dig>  isotigs and singletons had blastx alignments to protein sequences housed in the non-redundant protein database at an e-value threshold of 1e- <dig> or lower.
of the isotigs and singletons that had blastx alignments,  <dig>  were classified to class hexapoda .
annotation statistics for this assembly are summarized in table  <dig> figure  <dig>
histogram of isotig lengths generated from  <dig> flx reads.
<dig>  shotgun reads were assembled into  <dig>  isotigs  using newbler, which represent assembled transcripts.
the n <dig> isotig length was  <dig>  nt.
assembly metrics for
a. glabripennis
midgut transcriptome assembly generated using  <dig> reads only
a. glabripennis
midgut transcriptome  <dig> assembly
hybrid illumina/ <dig> transcriptome assembly
co-assembly with illumina paired-end sequences using trinity substantially improved the assembly metrics, resulting in the assembly of more full-length transcripts.
for this reason, discussion of the digestive and metabolic capabilities of a. glabripennis are focused mainly on genes and transcripts detected in the co-assembly and the 454-only assembly is used strictly for comparisons to other herbivorous insect gut transcriptomes.
the final 454/illumina co-assembly contained  <dig>  transcripts  ranging in length from  <dig> to  <dig>  nt with an n <dig> transcript length of  <dig> nt .
approximately  <dig>  transcripts had predicted protein coding regions and, of these, over  <dig>  transcripts contained full length open reading frames  with discernible start and stop codons.
these transcripts were classified to  <dig>  unigenes, bringing the average number of transcript isoforms per locus to  <dig> .
the highest number of isoforms detected for an individual gene/locus was  <dig> and transcripts assigned to this unigene were predicted to encode tropomyosin.
full assembly and annotation metrics for the 454-illumina hybrid assembly are presented in table  <dig>  of the unigenes predicted to contain full length or partial orfs,  <dig>   had blastp alignments at an e-value threshold of 1e- <dig> or lower, while  <dig> unigenes were predicted to encode rrnas and  <dig> transcripts were predicted to encode trnas.
approximately  <dig>   of the unigenes that had blastp alignments were classified to class hexapoda.
annotation metrics are presented in table  <dig>  to assess the potential completeness and quality of the larval midgut transcriptome assembly, several kegg metabolic pathways known to be conserved, functional, and complete in insects were examined to determine if all genes associated with these pathways were represented in the assembly.
full pathways for glycolysis and gluconeogenesis, pyrimidine metabolism, purine metabolism, pyruvate metabolism, the citric acid cycle, and phosphatidylinositol signaling systems were successfully constructed from protein coding transcripts in the assembly.figure  <dig>
histogram of transcript lengths generated from  <dig> flx and illumina paired end reads using trinity.
co-assembly of  <dig> shotgun reads and illumina paired-end sequences using trinity yielded  <dig>  transcripts  ranging in length from  <dig> to  <dig>  nt with an n <dig> transcript length of  <dig> nt.
this represented a substantial improvement over the  <dig> assembly and generated more full length transcripts in comparison.
assembly metrics for
a. glabripennis
midgut transcriptome illumina- <dig> co-assembly
a. glabripennis
midgut transcriptome illumina- <dig> hybrid assembly
overall, the most abundant pfam assignments detected in transcripts generated from the illumina/ <dig> co-assembly were primarily structural domains, including wd  <dig>  ankyrin, spectrin, and i-set, and domains associated with regulatory proteins, including reverse transcriptase, protein kinases, and zinc finger domain proteins.
the most dominant unigenes predicted to encode enzymes that were detected in this assembly were annotated as trypsins, dde superfamily endonucleases, carboxylesterases, cytochrome p450s, and glycoside hydrolase family one .
the majority of the unigenes detected in the midgut were assigned to the general functional prediction kog category, indicating that many of the unigenes detected in the midgut have not been definitively assigned to metabolic pathways and suggesting that they may be involved in novel or uncharacterized processes.
other highly abundant kog categories included signal transduction and carbohydrate transport and metabolism .
kog assignments of unigenes with putative signal peptides that could be involved in digestive processes were also conducted .figure  <dig>
relative abundance of the  <dig> most abundant pfam domain assignments in the
a. glabripennis
midgut transcriptome assembly.
unigenes generated from the 454/illumina assembly were scanned for pfam hmms.
over  <dig>  unigenes contained pfam domains.
kog assignments for midgut unigenes.
unigenes detected in the a. glabripennis midgut transcriptome were assigned to  <dig> different eukaryotic clusters of orthologous gene  categories.
overall, the majority of the unigenes were assigned to the signal transduction mechanisms and the general functional prediction only kog categories.
glycoside hydrolases and plant cell wall digesting enzymes
transcripts predicted to encode hemicellulases
over  <dig> different unigenes assigned to  <dig> gh families were identified, many of which have annotations consistent with involvement in plant cell wall degradation in the a. glabripennis midgut .
of particular interest are enzymes capable of degrading cellulose and hemicellulose, which are the two most predominant polysaccharides found in hardwoods.
few insect enzymes involved in large-scale degradation of xylan   <cit>  have been expressed and biochemically characterized in vitro.
through in-gel zymograms infused with birch xylan and madli-tof-based peptide sequencing, it was previously demonstrated that a. glabripennis was capable of producing at least one enzyme with hydrolytic activity directed at birch xylan, suggesting that the beetle has the endogenous capacity to degrade this hardwood polysaccharide  <cit> .
eight transcript isoforms of this gh  <dig> xylanase were detected in the transcriptome assembly, indicating that xylan degrading transcripts in a. glabripennis may be more numerous that previously reported.
the identification of these transcripts is significant and redefines the role of insects in processing xylan as it has generally been presumed that xylanases are only produced by microbial symbionts  <cit> .
it is possible that other gh transcripts detected in the a. glabripennis midgut may also encode xylanases or β-xylosidases.
for example, gh family  <dig> is predominately comprised of β-xylosidases  <cit>  and over  <dig> unigenes with gh  <dig> functional domains were detected in the a. glabripennis midgut transcriptome.
however, the ability to predict polysaccharide substrates and catalytic potentials of these enzymes was impeded by the lack of specific annotations in the databases because very few of the highest scoring blastp alignments have corresponding kegg e.c.
annotations.
more refined annotations would require in-depth functional genomics approaches.figure  <dig>
distribution of glycoside hydrolase families found in the
a. glabripennis
midgut transcriptome.
approximately  <dig> unique unigenes assigned to  <dig> different glycoside hydrolase families were detected.
the most dominant families were gh  <dig> and  <dig> while gh families  <dig>   <dig>   <dig>  and  <dig> were represented by a single unigene.
despite a. glabripennis’ endogenous ability to degrade long-chain xylan into shorter oligosaccharides, no insect-derived transcripts capable of releasing xylose monomers  from xylo-oligomers or converting xylose to atp or acetyl coa  were detected.
endogenous xylose utilization capabilities have not been described in cerambycids  <cit>  and it is generally hypothesized that these beetles depend on yeasts or other microbes in the gut to supply these enzymes  <cit> .
previous metagenomic profiling of the a. glabripennis midgut microbiota revealed that yeasts and lactic acid bacteria associated with the gut have the metabolic potential to ferment five-carbon sugars, converting them to ethanol and other compounds that could be used directly by a. glabripennis for energy and fatty acid production  <cit> .
furthermore, the presence of a large number of a. glabripennis-derived transcripts predicted to encode alcohol and aldehyde dehydrogenases could suggest a role in processing ethanol, acetate, and other metabolites generated through xylose fermentation by microbes colonizing the gut.
a. glabripennis also possesses full fatty acid biosynthetic pathways capable of incorporating acetate, acetyl coa, and microbial fermentation products into fatty acids.
other minor polysaccharides present in heartwood hemicellulose include glucuronoxylan, arabinoxylan, glucomannan, and xyloglucan, which are comprised of mannose, galactose, rhamnose, arabinose, glucuronic acid, and galacturonic acid residues residues  <cit> .
despite the fact that these polysaccharides make up a relatively minor component of plant cell walls in the heartwood of deciduous trees, many transcripts predicted to encode enzymes that release mannose and galactose residues from polysaccharides were detected in the a. glabripennis midgut.
for example,  <dig> unigenes predicted to encode gh  <dig> exo-β-galactosidases and β-galactosidases,  <dig> unigenes predicted to encode gh  <dig> α-mannosidases and mannosyl-oligosaccharide α- <dig> - <dig>  mannosidases, and  <dig> unigenes predicted to encode gh  <dig> α-mannosidases were detected and could be utilized to liberate mannose and galactose from the hemicellulose matrix.
other transcripts predicted to encode enzymes responsible for processing minor polysaccharides present in hemicellulose included β-mannosidase, lactase, β-thioglucosidase, and β-fucosidase  <cit> .
further, a. glabripennis actively expressed transcripts involved in processing and utilizing mannose and galactose sugars via glycolysis, suggesting that these sugars can be directly utilized for energy production.
transcripts predicted to encode cellulases and callases
like many other wood-feeding insects, a. glabripennis also produces a number of transcripts predicted to encode cellulases.
one of the most striking discoveries in the midgut transcriptome was the presence of six gh  <dig> cellulase unigenes, which all had highest scoring blastp alignments to gh  <dig> endo-β- <dig>  glucanases previously detected in the guts of other wood-feeding cerambycid beetles.
recombinant protein expression assays revealed that cellulases associated with other cerambycids  catalyzed the release of cello-oligomers from crystalline cellulose  <cit> .
these were not flagged as transcript isoforms by trinity, suggesting that genes encoding cellulases are represented in multiple copies in the a. glabripennis genome.
the purpose of this redundancy is unknown, but several other coleopterans harbor multiple copies of cellulases belonging to the same gh family  <cit> .
these enzymes could function under different physiological conditions, which is consistent with the contrasting ph and oxygen gradients that can be found in different regions of cerambycid midguts  <cit> .
alternatively, these unigenes could encode enzymes with slightly different catalytic capabilities that act on different cellulose macromolecule substructures  <cit> , target soluble or insoluble fractions of cellulose  <cit> , or process cello-oligomers into cellobiose .
in addition to gh  <dig> cellulases, a. glabripennis also expresses endo-β- <dig> -glucanases classified to gh families  <dig> and  <dig>  a. glabripennis also produced a large number of β-glucosidases, which hydrolyze cellobiose to release glucose.
the majority of these were classified to gh family  <dig>  which was the most abundant gh family detected in the a. glabripennis midgut transcriptome.
the overabundance of β-glucosidases relative to cellulases is common in many wood-feeding insects and wood-degrading microorganisms and is hypothesized to serve as a mechanism to indirectly enhance cellulase activity.
these β-glucosidases can often act quickly and efficiently to release glucose from cellobiose, reducing the impact of end product inhibition on cellulase activity  <cit> .
in addition to transcripts encoding enzymes predicted to disrupt major hardwood polysaccharides, several transcripts involved in degrading minor polysaccharides were detected.
for example, callose is a linear polysaccharide comprised of β- <dig>  and β- <dig>  linked glucose.
although callose is normally associated with the fleshy and metabolically active regions of plants, such as leaves and stems, it is also sporadically deposited in cell walls of secondary growth  <cit>  and represents suitable stores of glucose that could be liberated and assimilated by a. glabripennis.
several β- <dig>  and β- <dig>  glucanases detected in the midgut transcriptome could be involved in liberating glucose from this polysaccharide.
transcripts predicted to encode enzymes that contribute to lignin and phenylpropanoid degradation
while lignin is highly abundant in the heartwood of deciduous trees where the a. glabripennis larvae were collected for this study, no transcripts predicted to encode enzymes that are capable of yielding the types of lignin degradation products previously observed in a. glabripennis frass  <cit>  were detected.
a single laccase unigene with a signal peptide for extracellular targeting was detected in addition to several extracellular copper oxidase domain proteins, peroxidases, aldo-keto reductases, and alcohol dehydrogenases.
laccases are involved in lignin degradation in some white rot fungal taxa  <cit> , and an endogenous termite laccase capable of degrading lignin alkali and lignin phenolics was recently characterized  <cit> .
however, despite their reported ability to degrade lignin phenolics, many laccases require extracellular redox mediators to disrupt the non-phenolic β-aryl ether and c-c linkages that dominant hardwood lignins to yield the types of degradation products observed in a. glabripennis frass  <cit> .
while pathways for synthesis of these redox mediators have been identified in some white rot fungi, insects are unlikely to have the endogenous ability to synthesize them since all characterized laccase redox mediators are comprised of aromatic rings, which insects cannot inherently synthesize  <cit> .
based on these observations, we hypothesize that lignin degrading activities in the gut should be directly enhanced through interactions with microbial enzymes capable of synthesizing aromatic redox mediators or liberating aromatic compounds from lignin.
lignin metabolites released from the biopolymer can also be used as laccase mediators.
in addition to laccases,  <dig> unigenes predicted to encode aldo-keto reductases were detected in the a. glabripennis transcriptome.
in a recent study, expression levels of termite-produced aldo-keto reductases were correlated with feeding on wood and a recombinant aldo-keto reductase expressed in conjunction with other termite-derived cellulases enhanced sugar release from pine saw dust  <cit> , suggesting a role in enhancing lignocellulose digestion.
additionally, aldo-keto reductases have been shown to enhance xylose metabolism  <cit> , degrade xenobiotics and carbohydrates  <cit> , function as aryl alcohol dehydrogenases to facilitate the degradation of β-aryl ethers in lignin  <cit> , and are induced by exposure to phenolics and aromatic compounds in bacteria and yeasts  <cit> .
the abundance of these aldo-keto reductases in the midgut suggests that they could work in collaboration with other insect and microbial enzymes to facilitate penetration of lignin.
other enzymes encoded by the a. glabripennis transcriptome capable of disrupting bonds that cross-link hemicellulose to lignin  <cit>  included esterases, which liberate polysaccharide termini from the cell wall matrix, exposing them to hydrolytic enzymes and enhancing sugar release from this group of polysaccharides.
additionally,  <dig> unigenes predicted to encode alcohol dehydrogenases were detected in the midgut transcriptome; although these enzymes have not been shown to break linkages in polymeric lignin, they are hypothesized to enhance lignin oxidation in the guts of termites  <cit>  and they could serve similar roles in the a. glabripennis midgut.
finally, a number of extracellular peroxidases were also detected.
although the roles of insect-derived peroxidases in digestion and physiology are numerous and diverse  <cit> , direct roles for insect peroxidases in lignin degradation have not been explored.
lignin degradation releases phenylpropanoids , which are often toxic; however, a. glabripennis produces enzymes capable of degrading phenylpropanoid subunits, including epoxide hydrolases, which are often involved in polycyclic aromatic compound metabolism  <cit> .
other transcripts predicted to encode detoxification enzymes and antioxidants that could make contributions to degradation or inactivation of toxic lignin metabolites include alcohol dehydrogenases, aldehyde dehydrogenases, cytochrome p450s, glutathione s-transferases, catalases, carboxylesterases, enzymes involved in aromatic compound degradation, and glucuronosyl transferases.
additionally, aldo-keto reductases are capable of degrading phenolic compounds, including tannins and phenylpropanoids released from lignin degradation, and could be primed for detoxification roles.
transcripts predicted to encode detoxification enzymes
a. glabripennis eggs hatch directly beneath the bark of hardwood trees and first and second instars feed on primary phloem and xylem  <cit> , which serve as diffuse transport systems for toxic tree defensive compounds  <cit> , before tunneling into the heartwood as later instars.
though heartwood is not as metabolically active as the primary phloem and xylem, it accumulates potentially toxic secondary metabolites, including alkaloids, tannins, hydroxycinnamic acids, and phenolic glycosides, defending the plant from herbivory and protecting structural polysaccharides and biopolymers from biotic assaults  <cit> .
given that a. glabripennis completes development in over  <dig> different tree species  <cit>  and that it feeds in the phloem and xylem before eventually making its way into the heartwood, this insect must have mechanisms to detoxify or sequester the breadth of defensive plant secondary metabolites it encounters throughout its life cycle.
the gut represents the first line of defense against ingested host plant allelochemicals, pesticides, and other toxins and many transcripts predicted to encode detoxification enzymes were detected.
for example,  <dig> cytochrome p450-like unigenes were detected in the a. glabripennis midgut transcriptome.
these enzymes have versatile oxidoreductive properties, are highly involved in degrading lipophilic toxins  <cit> , and have been shown to confer resistance to pesticides  <cit>  as well as small aromatic toxins that can accumulate to high concentrations in the heartwood of trees .
the majority of the cytochrome p <dig> unigenes detected in the a. glabripennis midgut transcriptome had highest scoring blastp alignments to cytochrome p450s identified in the tribolium castaneum genome  <cit> , although the percent similarity at the amino acid level ranged from 34% to 66%, reflecting a relatively large degree of divergence from previously annotated cytochrome p450s .table  <dig>
cytochrome p <dig> annotations from
a. glabripennis
midgut transcriptome
t. castaneum
t. castaneum
l. decemlineata
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
l. decemlineata
l. decemlineata
l. decemlineata
l. decemlineata
l. decemlineata
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
l. decemlineata
t. castaneum
t. castaneum
t. castaneum
t. castaneum
t. castaneum
l. decemlineata
t. castaneum
t. castaneum
cytochrome p <dig> unigenes detected in the a. glabripennis midgut were putatively assigned to clans and families based on the annotations of the highest scoring blastp alignments  <cit> .
to conform to cytochrome p <dig> classification conventions, only alignments sharing ≥ 40% amino acid identity were used to annotate these unigenes.
while a handful of unigenes predicted to encode mitochondrial cytochrome p450s were identified, non-mitochondrial clans  <dig> and  <dig> were more highly represented in comparison.
non-mitochondrial cytochrome p450s were classified to five families, which were predominated by cyp <dig> , but also included cyp <dig> , cyp <dig> , cyp <dig> , and cyp <dig> .
cyp <dig> family genes are often present in multiple copies, occur in clusters in insect genomes  <cit> , have pivotal roles in the detoxification of host plant defensive chemicals such as xanthotoxins, gossypol, and chlorogenic acid   <cit> , and are often induced in herbivorous insects during periods of feeding.
unigenes predicted to encode carboxylesterases were the most dominant enzyme with putative involvement in detoxification processes detected in the midgut;  <dig> individual unigenes were predicted to encode these proteins.
although secreted carboxylesterases are generally involved in pheromone metabolism  <cit> , intracellular carboxylesterases are often implicated in pesticide and allelochemical metabolism and tolerance  <cit> .
over  <dig> of the carboxylesterase unigenes detected in the a. glabripennis midgut lacked secretory peptides and may be primed to serve detoxification roles.
for example, carboxylesterases are hypothesized to mediate resistance to phenolic glycosides in papilio canadensis <cit> , and are often found in high levels in the midgut  <cit> .
notably, trees in the family salicaceae, which include many of a. glabripennis’ preferred hosts , are notorious producers of phenolic glycosides   <cit>  and the abundance of carboxylesterases may promote colonization and survival in these hosts.
transcripts predicted to encode enzymes involved in conjugative deactivation of xenobiotic compounds were also detected from  <dig> and  <dig> unigenes predicted to encode udp-glucuronosyl transferases and glutathione s-transferases, respectively  <cit> .
these transferases can bind to xenobiotic compounds containing a diversity of functional groups, including oxygen, nitrogen, sulfur or carboxyl groups, enhancing their solubility and allowing them to be excreted or stored in the fat body for eventual elimination  <cit> .
they have been previously shown to detoxify cyanates and cinnamaldehydes  <cit> , which can be found in high concentrations in heartwood.
further, they can also conjugate and eliminate aromatic compounds, including tannins and toxic aromatic compounds stored in the heartwood or released from lignin degradation  <cit> .
transcripts predicted to encode enzymes involved in nitrogen acquisition
although nitrogen is scarce in the woody tissue of host trees, like other insects, a. glabripennis larvae have high demands for nitrogen during growth and development.
although nitrogen sources are present in very low abundance in woody tissues, microbes associated with the midgut have the metabolic capacity to synthesize all  <dig> amino acids, which could be assimilated and stored by the insect in the form of arylphorin and hexameric storage proteins encoded by a. glabripennis <cit> .
however, the insect possesses endogenous pathways for amino acid synthesis that may be complemented or augmented by microbial pathways, including complete pathways for the synthesis of alanine, aspartic acid, asparagine, proline, cysteine, glycine, and serine and for the synthesis of tyrosine from phenylalanine.
in addition, nearly complete pathways for the synthesis of arginine, glutamic acid and selenocysteine were detected, but argininosuccinate lyase, glutamate formiminotransferase, and selenocysteine synthases transcripts were absent.
these pathways may be incomplete because transcripts encoding these enzymes are simply not expressed in the midgut, they were expressed at low levels and were not detected at the sequencing depth obtained, or because they may be complemented by microbial enzymes that catalyze these reactions.
the phloem tissue where early instars feed is rich in amino acids relative to the heartwood  <cit>  where older larvae grow and develop.
therefore, recycling waste products of amino acid and nucleotide deamination reactions back into functional amino acids, nucleotides, and other nitrogen-containing compounds may be important to the nitrogen economy in a. glabripennis .
a number of transcripts with highest scoring blastp alignments to enzymes predicted to catalyze deamination reactions and liberating ammonia from a variety of nitrogen-containing compounds were detected in the midgut transcriptome.
these included adenine deaminases, cytosine deaminases, nitrilases, amidohydrolases, and chitin-degrading enzymes.
ammonia liberated from nitrogen-containing compounds could be directly converted to glutamine by glutamine synthetase and aspartic acid by aspartate ammonia ligase and indirectly incorporated into the synthesis of purines  <cit> , glutamate, alanine, asparagine, and proline by enzymes encoded by the insect.
furthermore, ammonia or amino acids constructed from recycled ammonia could be shuttled to microbes housed in the midgut to synthesize nonessential or essential amino acids, augmenting or complementing a. glabripennis’ physiological capabilities.
despite these potential contributions to nitrogen economy, the mechanisms of essential amino acid synthesis and recycling are not clear as the abundance of essential amino acids in woody tissue varies depending on tree species, but are significantly lower than the abundance of nonessential amino acids  <cit> .
as with other insects studied, no pathways for the synthesis and metabolism of essential amino acids  were detected in the midgut transcriptome, although these could be expressed elsewhere.
however, full pathways for the synthesis of nine essential amino acids were detected in the a. glabripennis midgut metagenome, including histidine, isoleucine, leucine, lysine, methionine, phenylalanine, threonine, tryptophan, and valine, and could serve as key sources of essential amino acids in this insect  <cit> .
furthermore, the gut community also contained numerous uricase and urease genes, which could be involved in recycling nitrogenous waste products produced by a. glabripennis or its gut microbes  <cit> .
although the waste products from a. glabripennis have not been biochemically characterized, all enzymes associated with the urea cycle were detected in the a. glabripennis midgut, including several arginase transcripts that catalyze the conversion of arginine to urea, suggesting that urea could be produced in this insect.
this urea pathway is also functional in the guts of several other insects  <cit> .figure  <dig>
proposed mechanisms for direct utilization of ammonia detected in the
a. glabripennis
midgut.
several methods for re-incorporating ammonia produced through nucleotide or amino acid deamination reactions into amino acid and nucleosides were reconstructed based on the presence of unigenes in the a. glabripennis midgut transcriptome.
ammonia can be directly integrated into glutamine by glutamine synthetase, which can subsequently be converted into glutamate by microbial-derived glutamate synthase detected in the a. glabripennis midgut metagenome.
glutamate could subsequently be converted into either proline by a series of reactions catalyzed by γ-glutamyl kinase, dehydrogenase, aminotransferase, and pyrroline-5-carboxylate or aspartate via transaminase, respectively.
furthermore, aspartate can be converted into inosine, which can be directly used for synthesis of purine nucleotides or asparagine via asparagine synthetase.
*microbial in origin.
1succinylaminoimidazole carboxamide ribotide synthase.
2aminoimidazole carboxamide ribotide transformylase.
transcripts predicted to encode several types of digestive proteinases, including serine and cysteine proteinases, were detected, which aid in protein acquisition from glycoproteins cross-linked in the cell wall matrix  <cit>  or from microbes housed in the midgut.
by far, serine proteinase unigenes were more expansive than cysteine proteinase unigenes and cystatins  were co-expressed and likely repress cysteine proteinase activities  <cit>  in the midgut.
these results are consistent with a previous study that reported high serine protease activities and the absence of cysteine and aspartic peptidase activities in the a. glabripennis midgut  <cit> .
however, the expression of both cysteine and aspartic proteinases in the midgut suggests that this insect still has the genetic capacity to produce these proteinases under certain conditions and these genes may be maintained in the genome as a strategy to combat digestive proteinase inhibitors produced by host plants.
analysis of transcripts involved in converting compounds in woody tissue into fatty acids and sterols was also conducted .
transcripts involved in facilitating interactions with gut microbes
although there is debate about how microbes associated with cerambycid guts contribute to digestive physiology  <cit> , a number of transcripts with putative involvement in mediating interactions with microbes were detected in the a. glabripennis midgut transcriptome.
several of these transcripts are likely involved in maintaining host-microbe homeostasis, including transcripts predicted to encode both antifungal and antibacterial proteins, dual oxidases  <cit> , mucin, which forms a protective barrier to protect the midgut from microbial invasion  <cit> , mpa <dig> allergen proteins, which have antimicrobial properties and are often upregulated during periods of stress  <cit> , and several encapsulation proteins involved in activating innate immune pathways.
seven unigenes predicted to encode hemocyanins were detected .
while they primarily function as oxygen carriers in crustaceans, hemocyanins are rare, but not completely absent in insect genomes and their physiological functions are not well characterized  <cit> .
despite their functional obscurity in insects, they can function as pro-phenol oxidases under certain circumstances, activating innate immune pathways and mediating insect-microbe interactions in the midgut  <cit> .
hemocyanins have also been hypothesized to serve roles in the degradation of lignin since transcripts encoding hemocyanins are highly expressed in a symbiont-free, wood-feeding marine isopod, but no direct involvement in this process has been demonstrated  <cit> .
additionally, signal peptides were not observed in any of the hemocyanin unigenes detected in a. glabripennis which, if not an assembly artifact, would preclude their involvement in extracellular digestive processes .figure  <dig>
structure of hemocyanin enzymes detected in the
a. glabripennis
midgut transcriptome.
seven unigenes predicted to encode hemocyanins were detected.
unlike many insect hemocyanins, the hemocyanins produced by a. glabripennis had intact and presumably functional copper binding domains.
these unigenes were predicted to have involvement in mediating interactions with gut microbes.
identification of highly expressed genes
transcripts originating from genomic and mitochondrial ribosomal rrnas were omitted from this analysis and the unigenes with the top  <dig> fpkm values  containing predicted coding regions were identified as highly expressed .
many of the highly expressed genes identified in the midgut have predicted involvement in stress and immune modulation and included unigenes predicted to encode several mpa <dig> allergen domain proteins, five carboxylesterases, two cathespins, two encapsulation related proteins, two mucin proteins, a lysosomal lipase, a cytochrome p <dig>  a thaumatin domain protein, and a lectin domain protein.
while carboxylesterases and cytochrome p450s have key involvements in detoxification and sterol acquisition, cathespins, encapsulation proteins, mucin proteins, lipases, thaumatin domain proteins and lectins are hypothesized to play fundamental roles in mediating host-microbe interactions.
not surprisingly, gh  <dig> and gh  <dig> cellulases and gh  <dig> β-glucosidases were also highly expressed in the midgut, reflecting the nutritional importance of cellulose to a. glabripennis.
gh  <dig> and gh  <dig> β-galactosidases were also highly expressed, suggesting that galactan polymers present in hardwood hemicellulose or on the cell surface of microbes are also crucial sources of sugar for this insect.
chitin deacetylase unigenes were highly abundant; these enzymes can liberate acetate from insect or fungal chitin, which can be recycled for energy or fatty acid production  <cit> .
several different types of digestive proteinases were also highly expressed and included m <dig> peptidases, m <dig> carboxypeptidases, serine proteinases, and cysteine proteinases, which likely serve key roles in nitrogen extraction from plant or microbial cell wall proteins.
in addition to these digestive proteins, several unigenes predicted to encode hypothetical proteins with unknown functions were abundant, suggesting that a. glabripennis encodes suites of novel proteins that could be relevant for digestive physiology and development.
glycoside hydrolase profile comparisons
through comparisons of transcriptome libraries sampled from a variety of herbivorous insects, no major trends were detected with regard to gh profiles and feeding habitats.
euclidean distances between insects that fed on similar substrates were large in many cases and reflected strong differences in gh compositions.
thus, the gut transcriptome libraries did not show any significant clustering trends by food source .
for example, agrilus planipennis and dendroctonus ponderosae both feed in phloem and were found in separate planes of the pca ordination, indicating that there were large differences in gh family composition between these two insects.
likewise, a. glabripennis, coptotermes formosanus, and reticulitermes flavipes all feed in wood and were also found in opposite quadrants in the pca ordination.
although these insects are all capable of producing endogenous cellulases, a. glabripennis produces different types of cellulases than the wood- and phloem-feeding insects compared in this study.
for example, the two termite species included in this analysis predominately produce gh  <dig> cellulases, while a. glabripennis produces gh  <dig>  gh  <dig>  and gh  <dig> cellulase transcripts and dendroctonous ponderosae produces gh  <dig> and gh  <dig> cellulase transcripts.figure  <dig>
multivariate comparison of glycoside hydrolase families detected in the gut transcriptomes of herbivorous insects.
principal components analysis was conducted to plot the abundance of each glycoside hydrolase family detected in the gut transcriptomes of wood-feeding insects.
pca values were plotted in sample space and variables are displayed as vectors on the pca biplot.
pca axis  <dig> and pca axis  <dig> explain  <dig> % and  <dig> % of the variation in the data, respectively.
monte carlo permutation procedure : p <  <dig>  for pca  <dig> and pca <dig>  no major clustering of insects feeding on similar food sources was detected.
despite the lack of clustering by feeding niche, there appeared to be some clustering by phylogenetic relatedness.
for example, most cerambycid and chrysomelid beetles were positioned along the positive x axis and, like a. glabripennis, some of these insects produce transcripts predicted to encode gh  <dig>   <dig>  and  <dig> cellulases, although they feed on very different parts of their host plants.
furthermore, it is interesting to note that gh  <dig> cellulases have not yet been found in any insect outside the order coleoptera, but the number of gh  <dig> cellulases unigenes detected in insect species from this order varied tremendously.
while gh  <dig> transcripts were not detected in association with many coleopterans, the chrysomelids gastrophysa viridula and callosobruchus maculatus encode one and four gh  <dig> unigenes, respectively.
phylogenetic analysis of translated proteinsequences revealed that, although chyrsomelid gh  <dig> cellulases and cerambycid  gh  <dig> cellulases share a common ancestor, chrysomelid cellulases have rapidly diverged from cerambycid cellulases .
in contrast, gh  <dig> cellulases within the cerambycidae seem to have multiplied and diversified through gene conversion or gene duplication events and are possibly more adapted to digesting highly insoluble cellulose associated with woody plants.figure  <dig>
phylogenetic analysis of gh  <dig> cellulases detected in the coleoptera.
unrooted maximum likelihood tree for insect, bacterial, and nematode derived gh  <dig> family proteins was constructed with garli  using the wag  <cit>  + i + f + g evolutionary model .
fully resolved bootstrap consensus trees were compiled using sum trees  and bootstrap values >  <dig> are displayed.
nodes are annotated with ncbi protein database accession numbers.
scale bar represents branch lengths.
furthering the hypothesis that the pca ordination was primarily driven by phylogenetic relatedness, family-specific trends in abundances of gh families were observed within the coleoptera.
in contrast to gh  <dig> cellulases, which seem to have multiplied in some cerambycid beetles, gh  <dig> and gh  <dig> cellulases were expressed as single copy genes in a. glabripennis .
in contrast, members of these gh families have multiplied and diversified in the chrysomelids and curculionid lineages, suggesting that coleopterans have undergone lineage specific adaptations to overcome challenges associated with different feeding regimes.
for example, the results of the gh  <dig> maximum likelihood analysis  suggest that gh  <dig> enzymes were likely encoded in the genome of the last common ancestor of coleopterans and that they underwent family-specific adaptations.
this scenario is supported since gh  <dig> proteins in each insect associated family formed their own supported clusters in the maximum likelihood tree.
in particular, genes encoding gh  <dig> enzymes were likely duplicated in the chrysomelidae.
all members of this family encode at least two gh  <dig> proteins and the branching topology suggests that the second gh  <dig> gene originated directly from the first.
likewise, gh  <dig> genes have also duplicated and proliferated throughout the chrysomelid and curcurlionid lineages, but the dynamics driving the evolution of this gh family seem to be more complex in comparison to the gh  <dig> family .
in some species, gh  <dig> genes have rapidly propagated and diversified , while in other cases, the insect expressed only a single copy of this gene .figure  <dig>
phylogenetic analysis of gh  <dig> cellulases detected in the coleoptera.
unrooted maximum likelihood tree for insect and bacterial derived gh  <dig> family proteins was constructed with garli  using the lg  <cit>  + i + f + g evolutionary model .
fully resolved bootstrap consensus trees were compiled using sum trees  and bootstrap values >  <dig> are displayed.
nodes are annotated with ncbi protein database accession numbers.
scale bar represents branch lengths.
phylogenetic analysis of gh  <dig> cellulases detected in the coleoptera.
unrooted maximum likelihood tree for coleopteran-derived gh  <dig> family proteins with garli  using the wag  <cit>  + g evolutionary model .
fully resolved bootstrap consensus trees were compiled using sum  and bootstrap values >  <dig> are displayed.
nodes are annotated with ncbi protein database accession numbers.
scale bar represents branch lengths.
the hypothesis that a. glabripennis benefits from microbial enzymes to facilitate nutrient acquisition is supported through the above comparisons.
for example, transcripts predicted to encode gh family  <dig>   <dig>  and  <dig> chitinases, β-hexosamidases, and glucosylceramidases are strongly associated with a. glabripennis.
although these chitin-degrading genes may be important for remodeling the gut peritrophic matrix, which is predominately composed of chitin, they could also play key roles in modulating interactions with fungal taxa associated with the midgut, including yeasts and fusarium solani, a soft rot fungal symbiont of a. glabripennis <cit> .
we hypothesize that the predominance of chitinases allows a. glabripennis to derive a portion of its carbohydrate and nitrogen resources from fungal chitinous cell walls, which are composed of polymers of n-acetylglucosamine.
non-entomopathogenic fungi associated with wood-feeding insects have been previously hypothesized to concentrate and/or recycle nitrogen  <cit>  and thus we also hypothesize that the ascomycota fungal strains found in association with the a. glabripennis midgut serve these same purposes.
multivariate transcriptome comparisons of gene ontology annotations
like the gh analysis, multivariate comparisons of level four go categories in midgut transcriptome libraries sampled from a variety of herbivorous insects revealed no significant clustering of transcriptome libraries by feeding habitat.
however, phylogenetic relatedness alone did not explain the observed pattern of clustering achieved for the transcriptome comparisons .
in addition, subsets of different go categories were enriched in each transcriptome library included in the comparison, while many of the go categories were present in approximately the same abundances in each library.
together, these findings suggest that most insects possess similar repertoires of gene families and that these genes have adapted in lineage specific manners optimal for overcoming digestive and nutritional challenges associated with specific feeding habitats and ecological niches.
for example, although most insects produce a similar number of gh unigenes , the gh family-level comparison uggested that each insect produced its own unique gh profile.
other go categories that are present in similar abundances in all insects included in this comparison include 4-α-glucanotransferases, heme binding and transporting proteins, and regulatory genes .figure  <dig>
two-way cluster analysis of level four gene ontology terms from herbivorous insect gut transcriptomes.
two-way cluster analysis of level four gene ontology terms was performed to identify potential correlations between abundances of level four go terms and feeding niche.
clustering did not appear to be driven by food source or phylogenetic relatedness.
however, the a. glabripennis midgut transcriptome is notably distinct from the other gut transcriptomes included in this comparison.
despite the lack of clustering by food source or phylogenetic relatedness, several trends were detected that distinguish a. glabripennis from the rest of the insects included in this comparison that could be pivotal to its ability to digest lignocellulose and other wood polysaccharides and extract nutrients from a broad range of deciduous host trees.
for example, in the midgut of a. glabripennis, more unigenes and transcript isoforms were produced with predicted monooxygenase and oxidoreductase activities relative to other insects included in this comparison, which could be relevant to its ability to detoxify allelochemicals from its broad range of host plants.
further examination of pfam domain abundances in each library revealed that unigenes and transcript isoforms predicted to encode carboxylesterases and cytochrome p450s were more abundant in the a. glabripennis midgut than many of the other insect libraries sampled.
a. glabripennis also has the broadest host range of any insect included in this comparison, suggesting that it needs to encode a broader arsenal of detoxification enzymes relative to other insects included in this comparison.
several unigenes predicted to encode digestive peptidases, ligases , and protein transporters were also overrepresented relative to other insect transcriptome libraries, which could be relevant for digesting and assimilating proteins produced by microbes associated with the midgut or from plant cell walls  <cit> .
these digestive peptidases are also overrepresented in g. viridula, c. formosanus, and m. sexta transcriptome libraries.
unigenes associated with hydrolase activity  were also highly abundant in a. glabripennis, many of which were predicted to encode atpases and other nucleosidases, dna binding proteins, rna binding proteins, nucleotide binding proteins, and transferases involved in transferring phosphorous containing groups.
the high abundance of unigenes for these nucleotide binding proteins and nucleosidases is likely associated with the high numbers of unigenes predicted to encode reverse transcriptases, transposases, and integrases that were detected in the midgut.
finally, unigenes predicted to encode proteins with substrate-specific and active transmembrane transporter activities, including major facilitator family  transporters, were also highly abundant in the a. glabripennis midgut.
mfs transporters are a diverse group of carriers involved in the absorption of small solutes, including sugar, aromatic amino acids, and other small compounds, which may be involved in assimilation and utilization of small microbial metabolites and/or small metabolites released from the degradation of woody tissue.
taken together, the differences in gh family and level four go compositions among insects with similar feeding regimes suggest that the ability to degrade polysaccharides found in woody tissue evolved through lineage specific adaptations rather than through convergent evolutionary processes.
insights into interactions with microbes
in a recent publication documenting the metabolic potential of the microbiota associated with the a. glabripennis midgut, we discovered a taxonomically diverse assemblage of bacteria and fungi primed to make key contributions to digestive physiology, wood digestion, and nutrient acquisition in this system  <cit> .
however, the potential contributions of the beetle to digestive processes were unknown as no comprehensive transcriptomic or genomic resources are currently available for cerambycids; thus, it is difficult to determine precisely how the gut microbes augment or complement physiological processes in the gut.
here, we demonstrate that a. glabripennis also has a relatively versatile metabolic potential, producing a number of its own key digestive enzymes.
in particular, a. glabripennis endogenously produces a diverse suite of detoxification enzymes, including cyp 450s and carboxylesterases, many of which act promiscuously on a diverse array of toxins and are hypothesized to facilitate the degradation and inactivation of host tree defensive compounds.
in contrast, comparatively few cyp 450s and carboxylesterases were detected in the microbiome of the beetle gut.
instead, the community predominantly produces enzymes that target specific classes of toxins, including salicylates and other compounds involved in salicylic acid mediated defenses, as well as arsenic, cyanide, and cyanoamino acids.
these could enhance endogenous detoxification processes driven by beetle-derived enzymes and could also serve to detoxify many of the defensive chemicals that act directly on the microbial community , thereby having indirect deleterious impacts on the beetle.
furthermore, a. glabripennis produces several different types of plant cell wall degrading enzymes, including cellulases, β-glucosidases, and xylanases that likely serve a major role in breaking down these prominent wood polysaccharides.
however, the community also produces a diverse array of plant cell wall degrading enzymes, which could not only serve to augment the cell wall degrading capacities of a. glabripennis, but could also complement these endogenous processes.
for example, while the beetle produces enzymes predicted to target β- <dig>   <dig> linkages in xylan, several genes predicted to encode enzymes with activity towards ester linkages in xylan were detected in association with the gut community and included acetyl xylan and feruloyl esterases.
the presence of these esterases could expedite degradation of xylan in the midgut by cleaving ester linkages that cross-link adjacent hemicellulose chains, exposing sugar residues to beetle-derived hydrolytic enzymes and enhancing digestibility.
in addition, although xylose is the second most abundant sugar in deciduous trees, no putative β-xylosidases were detected in the a. glabripennis midgut transcriptome.
however, gene tags predicted to encode β-xylosidases and enzymes associated with the pentose phosphate pathway were highly abundant in the midgut community, which could possibly serve to convert xylose sugars into compounds that can be directly used by a. glabripennis for energy and fatty acid production.
in addition, the potential for cooperation between a. glabripennis and its gut microbes was also noted with regards to lignin degradation and nitrogen recycling.
for example, a. glabripennis produces a small pool of transcripts that could facilitate degradation of lignin, including aldo-keto reductases, laccases, and peroxidases.
alone, these enzymes can only facilitate small-scale degradation of the phenolic linkages that comprise a relatively small percentage of the lignin biopolymer.
however, larger scale lignin degrading reactions could be accomplished through interactions with enzymes produced by the gut community.
for example, the microbial community contained a high abundance of gene tags involved in synthesizing aromatic redox mediators, which could work in tandem with insect-derived laccases to facilitate degradation of the β-aryl ether linkages that predominantly comprise lignin, and genes tags predicted to directly degrade β-aryl ethers and other abundant linkages in lignin.
the beetle requires essential nutrients to complete its development, which are lacking in woody tissue.
complete pathways for the synthesis of all  <dig> major amino acids and several essential vitamins were detected in the gut community, which could augment the production of non-essential amino acids endogenously synthesized by the beetle and contribute to the production of essential amino acids and other essential nutrients, which the beetle cannot synthesize.
although a. glabripennis has abilities to scavenge ammonia, reincorporating it into nonessential amino acids and nucleotides, the community has an expanded capacity to recycle nitrogenous waste products including urea, uric acid, xanthine, and arginine, to potentially reincorporate ammonia into both essential and non-essential amino acids, nucleotides, and other nitrogen-containing compounds.
thus, the community could serve as an additional source of nitrogen and non-esseential amino acids in this high c:n environment.
the microbiota also has the capacity to fix atmospheric nitrogen, providing additional sources of nitrogen to both the beetle and members of the gut community.
while this suggests that the beetle collaborates with its gut microbes to facilitate survival in woody tissue, it is unknown which of these microbial pathways are metabolically active in the a. glabripennis midgut.
because the transcriptome library was sampled primarily from midgut tissue, few microbial transcripts were detected in this dataset.
to gain further insight into transcriptional activity of microbes in the gut, we are currently sequencing rna collected from the midgut contents to prepare a sample that is more enriched in microbial rnas.
this was the same approach that was used to sequence the metagenome of the a. glabripennis midgut and we anticipate that this will expand our abilities to more conclusively model the interactions between a. glabripennis and its gut microbes that enhance fitness and/or are required for survival in woody tissue.
CONCLUSIONS
the a. glabripennis midgut transcriptome provides the first comprehensive insight into the endogenous digestive capabilities of wood-boring cerambycid larvae as they feed in a highly lignified and nutritionally deficient environment.
comparative transcriptome analysis clearly distinguished the a. glabripennis midgut transcriptome from the gut transcriptome libraries of other herbivorous insects that have been previously sampled for sequencing, which may contribute to its long life cycle and ability to feed and develop in a highly lignified food source.
our results highlighted gene categories that were enriched in the a. glabripennis midgut transcriptome, which were hypothesized to make key contributions to this insect’s lifestyle, including its ability to colonize a broad range of living host trees.
for example, unigenes predicted to encode monooxygenases, carboxylesterases, heat shock proteins, and other detoxification enzymes were highly abundant in the midgut relative to unigenes expressed in the gut transcriptomes of many other herbivorous insects included in this comparison.
furthermore, a. glabripennis expressed its own unique profile of gh unigenes for liberating sugars from woody tissue.
our results also highlight deficiencies in endogenous digestive and metabolic pathways that could be supplied by microbes associated with the gut, including enzymes for xylose and pentose sugar utilization, enzymes that facilitate the degradation of lignin, pathways for the synthesis of essential amino acids and nutrients, and pathways for nutrient scavenging.
for these reasons, we hypothesize that enzymes derived from midgut-associated microbes complement the expression of insect genes and serve vital roles in the digestive physiology of a. glabripennis.
this study also substantially expands the genetic resources available for coleopterans by providing transcriptome data for an insect that feeds in the heartwood of healthy host trees.
the results presented here may serve as a source for bioprospecting of novel enzymes to enhance industrial biofuels productions and for the pursuit of novel targets for controlling this serious pest.

 class hexapoda is one of the most ancient and diverse animal lineages on the planet, containing organisms capable of occupying many recalcitrant niches and persisting under intense environmental conditions including extreme temperatures, periods of desiccation, and exposure to toxins. many of its members are capable of thriving on suboptimal, nutritionally-deficient substrates, including wood-boring beetles belonging to family cerambycidae that feed exclusively on woody tissue. specifically, beetles in the genus anoplophora are wood-boring insects of great interest because some of its members preferentially target healthy host trees and have relatively broad host ranges. for example, the asian longhorned beetle was introduced from china into the united states, canada, and several countries in europe and has been documented to complete development in approximately deciduous tree species worldwide, including several genera commonly planted as feedstock. acer spp. are the predominant hosts in the introduced range. this beetle poses a significant threat to urban streetscapes, has the potential to destroy up to 35% of the urban tree canopy in its introduced range, and has already caused millions of dollars in damage to urban landscapes. wood-borers, like a. glabripennis, are especially challenging to control in both their natural and invasive ranges because the larvae spend 1–2 years living deep inside their host trees. natural enemies are rare and treatment of host trees with systemic insecticides is costly and has variable efficacy against a. glabripennis larvae. the most effective method for eradication is destruction of infested and nearby host trees and implementation of strict quarantine measures to contain the infestation. therefore, understanding the digestive physiology of this cerambycid at the genetic level is paramount to devising novel control strategies. because a. glabripennis spends the majority of its lifecycle in the larval stage and feeds primarily in the heartwood of a broad range of healthy deciduous trees, it must overcome challenges of digesting intractable woody tissue in order to acquire sufficient nutrients to complete development. glucose is a predominant wood sugar, but it is present in the form of complex polysaccharides, including cellulose, hemicellulose, callose, and pectin, which are inherently difficult to digest and require a complex of hydrolytic enzymes for efficient degradation and liberation of sugar monomers. extensive hydrogen bonding coupled with linear configurations increases the crystallinity of these cell wall polysaccharides and decreases their permeability, further hindering the activity of hydrolytic enzymes. plant cell wall polysaccharides are further protected from hydrolytic enzymes by lignin, a biopolymer containing over types of chemical bonds that is extensively cross-linked to both cellulose and hemicellulose, shielding them from digestion. due to the random, heterogeneous nature of these cross-linkages and the high resilience of carbon-carbon and β-aryl ether linkages that dominate this macromolecule, lignin polymers can only be efficiently degraded through oxidative depolymerization, a process that has only been conclusively documented to be catalyzed by enzymes produced by a small number of wood degrading fungi. nitrogen is also extremely limited in woody tissues and plant cell wall proteins are intricately cross-linked with lignin and cellulose, making them difficult to access. other essential nutrients, including fatty acids, sterols, and vitamins are present in low concentrations or are completely absent. lastly, wood-feeding insects must overcome plant secondary metabolites that often accumulate to high concentrations in the heartwood through detoxification or sequestration processes. many wood-feeding beetles cultivate extracellular symbiotic fungi to facilitate digestion of woody tissue and nutrient acquisition, which are carried in mycangia or other specialized structures on their body. for example, bark beetles utilize a mass attack strategy, in which a mycangial fungus is directly inoculated into a host tree during oviposition to facilitate pre-digestion of woody tissue and mitigation of host tree defenses. an alternative strategy is to preferentially colonize stressed trees whose woody components have already been pre-digested by wood-rotting microbes. however, a. glabripennis is distinct from many other wood-feeding beetles in the sense that a single larvae can successfully develop in a healthy tree without requiring mass attack and the majority of the challenging reactions, including digestion of lignocellulose and hemicellulose and detoxification of plant metabolites, can occur within the gut itself. while the midgut community associated with a. glabripennis has the metabolic potential to overcome many of the challenges associated with feeding in woody tissue, including degradation of lignin, cellulose, and hemicellulose and acquisition of nitrogen and other essential nutrients, the contributions of insect-derived digestive and nutrient acquiring enzymes cannot be ignored since insects themselves can produce a diverse array of digestive enzymes, including cellulases, hemicellulases, pectinases, and enzymes that enhance lignin degradation. insects have also evolved other sophisticated abilities to evade host plant defenses and often possess extensive suites of enzymes involved in detoxification of plant metabolites and phytohormones, digestive proteinase inhibitors, and cyanates and cyanoamino acids, as well as enzymes capable of disrupting jasmonic acid signaling pathways. furthermore, insects produce many cytochrome p450s, which are integrally involved in xenobiotic metabolic processes that ultimately lead to oxidative destruction of toxic compounds, including plant derived secondary metabolites and pesticides. the primary goals of this study were to survey the endogenous digestive and physiological capabilities of larval a. glabripennis through shotgun sequencing of midgut derived messenger rna and to identify insect-derived genes that are highly expressed in the midgut while actively feeding in wood. the a. glabripennis midgut transcriptome library was also compared to all publically available transcriptome libraries sampled from other plant feeding insects to identify core groups of genes that are associated with digestive processes that could facilitate nutrient recovery from woody tissue regardless of insect taxa. this study represents an important addendum to the growing database of genomic and transcriptomic resources available for coleopterans and fills an important gap, representing the first transcriptome sampled from a wood-feeding cerambycid and the first comprehensive analysis of endogenous genes associated with wood-feeding in insects. these findings offer unique opportunities to bioprospect for enzymes that could be exploited for cellulosic biofuel production or other industrial processes, and to develop novel control methods for this destructive wood-boring pest and other wood-feeding insects.  454- and illumina-based transcriptome sequencing to develop a comprehensive profile of the endogenous digestive and physiological capabilities of a. glabripennis, mrna was collected from the midguts of third instar larvae feeding in the heartwood of a preferred host and was sequenced using both roche pyrosequencing and illumina technologies. in total, shotgun sequence reads were produced using the roche flx platform using two separate runs. reads, ranging in length from to nt, were generated on a half plate and reads ranging from to nt, were generated on a quarter plate. these runs correspond to e4gebh sff and e5ty7pb sff from sra, respectively. reads from both runs were pooled and were quality filtered and assembled together. approximately of the total flx reads passed quality filtering and were utilized in the assembly. to enhance sequencing depth and acquire a more complete inventory of the endogenous digestive and metabolic capabilities of a. glabripennis, million 100 bp paired illumina reads with a library insert size of nucleotides were generated on a single lane using the illumina hiseq. after quality filtering and adapter removal, over million read pairs remained and were utilized in downstream processing and analyses. digital k-mer normalization reduced the number of illumina read pairs to, which were ultimately used for co-assembly with the flx reads. assembly and annotation statistics assembly and annotation statistics for comparative transcriptomics to facilitate comparisons to transcriptome libraries prepared from the guts of other herbivorous insects, which were derived solely from reads, the reads were first assembled and analyzed without the illumina reads. of the shotgun reads generated through pyrosequencing, approximately reads assembled into contigs, ranging in length from nt to nt with an n contig length of nt. assembled contigs that shared common reads were placed into isogroups. these contigs are often broken at branch points between exon boundaries in multiple transcript isoforms from the same unigene. contig branch structures within each isogroup were then traversed to create isotigs, which represent unique assembled transcripts or transcript fragments. the n isotig length was nt and isotigs were grouped into isogroups, representing a gene locus or unigene. of these isogroups, were comprised of a single transcript isoform and the average number of isotigs within an isogroup was. the maximum number of isotigs classified to the same isogroup was for downstream comparative analyses, isogroups were treated as unigenes and isotigs associated with the same isogroup were treated as transcript isoforms. roughly reads were singletons and were not included in the assembly. of the singletons, approximately reads were flagged as high quality and, to increase the amount of information present in the transcriptome dataset, these singleton reads were concatenated to the assembly and the pooled dataset was utilized in downstream transcriptome comparisons. assembly metrics from the 454-based assembly are presented in table  after clustering the isotigs and high quality singletons with cd-hit-est using a sequence similarity threshold of to group transcripts that likely represented allelic variants of the same gene, the total number of isotigs and singletons was reduced to around. seventy-eight of these isotigs and reads were classified as ribosomal rnas, while none were classified as trnas. roughly isotigs and singletons had blastx alignments to protein sequences housed in the non-redundant protein database at an e-value threshold of 1e- or lower. of the isotigs and singletons that had blastx alignments, were classified to class hexapoda. annotation statistics for this assembly are summarized in table  figure histogram of isotig lengths generated from flx reads.  shotgun reads were assembled into isotigs using newbler, which represent assembled transcripts. the n isotig length was nt. assembly metrics for a. glabripennis midgut transcriptome assembly generated using reads only a. glabripennis midgut transcriptome assembly hybrid illumina/ transcriptome assembly co-assembly with illumina paired-end sequences using trinity substantially improved the assembly metrics, resulting in the assembly of more full-length transcripts. for this reason, discussion of the digestive and metabolic capabilities of a. glabripennis are focused mainly on genes and transcripts detected in the co-assembly and the 454-only assembly is used strictly for comparisons to other herbivorous insect gut transcriptomes. the final 454/illumina co-assembly contained transcripts ranging in length from to nt with an n transcript length of nt. approximately transcripts had predicted protein coding regions and, of these, over transcripts contained full length open reading frames with discernible start and stop codons. these transcripts were classified to unigenes, bringing the average number of transcript isoforms per locus to. the highest number of isoforms detected for an individual gene/locus was and transcripts assigned to this unigene were predicted to encode tropomyosin. full assembly and annotation metrics for the 454-illumina hybrid assembly are presented in table  of the unigenes predicted to contain full length or partial orfs, had blastp alignments at an e-value threshold of 1e- or lower, while unigenes were predicted to encode rrnas and transcripts were predicted to encode trnas. approximately of the unigenes that had blastp alignments were classified to class hexapoda. annotation metrics are presented in table  to assess the potential completeness and quality of the larval midgut transcriptome assembly, several kegg metabolic pathways known to be conserved, functional, and complete in insects were examined to determine if all genes associated with these pathways were represented in the assembly. full pathways for glycolysis and gluconeogenesis, pyrimidine metabolism, purine metabolism, pyruvate metabolism, the citric acid cycle, and phosphatidylinositol signaling systems were successfully constructed from protein coding transcripts in the assembly.figure histogram of transcript lengths generated from flx and illumina paired end reads using trinity. co-assembly of shotgun reads and illumina paired-end sequences using trinity yielded transcripts ranging in length from to nt with an n transcript length of nt. this represented a substantial improvement over the assembly and generated more full length transcripts in comparison. assembly metrics for a. glabripennis midgut transcriptome illumina- co-assembly a. glabripennis midgut transcriptome illumina- hybrid assembly overall, the most abundant pfam assignments detected in transcripts generated from the illumina/ co-assembly were primarily structural domains, including wd ankyrin, spectrin, and i-set, and domains associated with regulatory proteins, including reverse transcriptase, protein kinases, and zinc finger domain proteins. the most dominant unigenes predicted to encode enzymes that were detected in this assembly were annotated as trypsins, dde superfamily endonucleases, carboxylesterases, cytochrome p450s, and glycoside hydrolase family one. the majority of the unigenes detected in the midgut were assigned to the general functional prediction kog category, indicating that many of the unigenes detected in the midgut have not been definitively assigned to metabolic pathways and suggesting that they may be involved in novel or uncharacterized processes. other highly abundant kog categories included signal transduction and carbohydrate transport and metabolism. kog assignments of unigenes with putative signal peptides that could be involved in digestive processes were also conducted.figure relative abundance of the most abundant pfam domain assignments in the a. glabripennis midgut transcriptome assembly. unigenes generated from the 454/illumina assembly were scanned for pfam hmms. over unigenes contained pfam domains. kog assignments for midgut unigenes. unigenes detected in the a. glabripennis midgut transcriptome were assigned to different eukaryotic clusters of orthologous gene categories. overall, the majority of the unigenes were assigned to the signal transduction mechanisms and the general functional prediction only kog categories. glycoside hydrolases and plant cell wall digesting enzymes transcripts predicted to encode hemicellulases over different unigenes assigned to gh families were identified, many of which have annotations consistent with involvement in plant cell wall degradation in the a. glabripennis midgut. of particular interest are enzymes capable of degrading cellulose and hemicellulose, which are the two most predominant polysaccharides found in hardwoods. few insect enzymes involved in large-scale degradation of xylan have been expressed and biochemically characterized in vitro. through in-gel zymograms infused with birch xylan and madli-tof-based peptide sequencing, it was previously demonstrated that a. glabripennis was capable of producing at least one enzyme with hydrolytic activity directed at birch xylan, suggesting that the beetle has the endogenous capacity to degrade this hardwood polysaccharide. eight transcript isoforms of this gh xylanase were detected in the transcriptome assembly, indicating that xylan degrading transcripts in a. glabripennis may be more numerous that previously reported. the identification of these transcripts is significant and redefines the role of insects in processing xylan as it has generally been presumed that xylanases are only produced by microbial symbionts. it is possible that other gh transcripts detected in the a. glabripennis midgut may also encode xylanases or β-xylosidases. for example, gh family is predominately comprised of β-xylosidases and over unigenes with gh functional domains were detected in the a. glabripennis midgut transcriptome. however, the ability to predict polysaccharide substrates and catalytic potentials of these enzymes was impeded by the lack of specific annotations in the databases because very few of the highest scoring blastp alignments have corresponding kegg e.c. annotations. more refined annotations would require in-depth functional genomics approaches.figure distribution of glycoside hydrolase families found in the a. glabripennis midgut transcriptome. approximately unique unigenes assigned to different glycoside hydrolase families were detected. the most dominant families were gh and while gh families and were represented by a single unigene. despite a. glabripennis’ endogenous ability to degrade long-chain xylan into shorter oligosaccharides, no insect-derived transcripts capable of releasing xylose monomers from xylo-oligomers or converting xylose to atp or acetyl coa were detected. endogenous xylose utilization capabilities have not been described in cerambycids and it is generally hypothesized that these beetles depend on yeasts or other microbes in the gut to supply these enzymes. previous metagenomic profiling of the a. glabripennis midgut microbiota revealed that yeasts and lactic acid bacteria associated with the gut have the metabolic potential to ferment five-carbon sugars, converting them to ethanol and other compounds that could be used directly by a. glabripennis for energy and fatty acid production. furthermore, the presence of a large number of a. glabripennis-derived transcripts predicted to encode alcohol and aldehyde dehydrogenases could suggest a role in processing ethanol, acetate, and other metabolites generated through xylose fermentation by microbes colonizing the gut. a. glabripennis also possesses full fatty acid biosynthetic pathways capable of incorporating acetate, acetyl coa, and microbial fermentation products into fatty acids. other minor polysaccharides present in heartwood hemicellulose include glucuronoxylan, arabinoxylan, glucomannan, and xyloglucan, which are comprised of mannose, galactose, rhamnose, arabinose, glucuronic acid, and galacturonic acid residues residues. despite the fact that these polysaccharides make up a relatively minor component of plant cell walls in the heartwood of deciduous trees, many transcripts predicted to encode enzymes that release mannose and galactose residues from polysaccharides were detected in the a. glabripennis midgut. for example, unigenes predicted to encode gh exo-β-galactosidases and β-galactosidases, unigenes predicted to encode gh α-mannosidases and mannosyl-oligosaccharide α- - mannosidases, and unigenes predicted to encode gh α-mannosidases were detected and could be utilized to liberate mannose and galactose from the hemicellulose matrix. other transcripts predicted to encode enzymes responsible for processing minor polysaccharides present in hemicellulose included β-mannosidase, lactase, β-thioglucosidase, and β-fucosidase. further, a. glabripennis actively expressed transcripts involved in processing and utilizing mannose and galactose sugars via glycolysis, suggesting that these sugars can be directly utilized for energy production. transcripts predicted to encode cellulases and callases like many other wood-feeding insects, a. glabripennis also produces a number of transcripts predicted to encode cellulases. one of the most striking discoveries in the midgut transcriptome was the presence of six gh cellulase unigenes, which all had highest scoring blastp alignments to gh endo-β- glucanases previously detected in the guts of other wood-feeding cerambycid beetles. recombinant protein expression assays revealed that cellulases associated with other cerambycids catalyzed the release of cello-oligomers from crystalline cellulose. these were not flagged as transcript isoforms by trinity, suggesting that genes encoding cellulases are represented in multiple copies in the a. glabripennis genome. the purpose of this redundancy is unknown, but several other coleopterans harbor multiple copies of cellulases belonging to the same gh family. these enzymes could function under different physiological conditions, which is consistent with the contrasting ph and oxygen gradients that can be found in different regions of cerambycid midguts. alternatively, these unigenes could encode enzymes with slightly different catalytic capabilities that act on different cellulose macromolecule substructures, target soluble or insoluble fractions of cellulose, or process cello-oligomers into cellobiose. in addition to gh cellulases, a. glabripennis also expresses endo-β- -glucanases classified to gh families and a. glabripennis also produced a large number of β-glucosidases, which hydrolyze cellobiose to release glucose. the majority of these were classified to gh family which was the most abundant gh family detected in the a. glabripennis midgut transcriptome. the overabundance of β-glucosidases relative to cellulases is common in many wood-feeding insects and wood-degrading microorganisms and is hypothesized to serve as a mechanism to indirectly enhance cellulase activity. these β-glucosidases can often act quickly and efficiently to release glucose from cellobiose, reducing the impact of end product inhibition on cellulase activity. in addition to transcripts encoding enzymes predicted to disrupt major hardwood polysaccharides, several transcripts involved in degrading minor polysaccharides were detected. for example, callose is a linear polysaccharide comprised of β- and β- linked glucose. although callose is normally associated with the fleshy and metabolically active regions of plants, such as leaves and stems, it is also sporadically deposited in cell walls of secondary growth and represents suitable stores of glucose that could be liberated and assimilated by a. glabripennis. several β- and β- glucanases detected in the midgut transcriptome could be involved in liberating glucose from this polysaccharide. transcripts predicted to encode enzymes that contribute to lignin and phenylpropanoid degradation while lignin is highly abundant in the heartwood of deciduous trees where the a. glabripennis larvae were collected for this study, no transcripts predicted to encode enzymes that are capable of yielding the types of lignin degradation products previously observed in a. glabripennis frass were detected. a single laccase unigene with a signal peptide for extracellular targeting was detected in addition to several extracellular copper oxidase domain proteins, peroxidases, aldo-keto reductases, and alcohol dehydrogenases. laccases are involved in lignin degradation in some white rot fungal taxa, and an endogenous termite laccase capable of degrading lignin alkali and lignin phenolics was recently characterized. however, despite their reported ability to degrade lignin phenolics, many laccases require extracellular redox mediators to disrupt the non-phenolic β-aryl ether and c-c linkages that dominant hardwood lignins to yield the types of degradation products observed in a. glabripennis frass. while pathways for synthesis of these redox mediators have been identified in some white rot fungi, insects are unlikely to have the endogenous ability to synthesize them since all characterized laccase redox mediators are comprised of aromatic rings, which insects cannot inherently synthesize. based on these observations, we hypothesize that lignin degrading activities in the gut should be directly enhanced through interactions with microbial enzymes capable of synthesizing aromatic redox mediators or liberating aromatic compounds from lignin. lignin metabolites released from the biopolymer can also be used as laccase mediators. in addition to laccases, unigenes predicted to encode aldo-keto reductases were detected in the a. glabripennis transcriptome. in a recent study, expression levels of termite-produced aldo-keto reductases were correlated with feeding on wood and a recombinant aldo-keto reductase expressed in conjunction with other termite-derived cellulases enhanced sugar release from pine saw dust, suggesting a role in enhancing lignocellulose digestion. additionally, aldo-keto reductases have been shown to enhance xylose metabolism, degrade xenobiotics and carbohydrates, function as aryl alcohol dehydrogenases to facilitate the degradation of β-aryl ethers in lignin, and are induced by exposure to phenolics and aromatic compounds in bacteria and yeasts. the abundance of these aldo-keto reductases in the midgut suggests that they could work in collaboration with other insect and microbial enzymes to facilitate penetration of lignin. other enzymes encoded by the a. glabripennis transcriptome capable of disrupting bonds that cross-link hemicellulose to lignin included esterases, which liberate polysaccharide termini from the cell wall matrix, exposing them to hydrolytic enzymes and enhancing sugar release from this group of polysaccharides. additionally, unigenes predicted to encode alcohol dehydrogenases were detected in the midgut transcriptome; although these enzymes have not been shown to break linkages in polymeric lignin, they are hypothesized to enhance lignin oxidation in the guts of termites and they could serve similar roles in the a. glabripennis midgut. finally, a number of extracellular peroxidases were also detected. although the roles of insect-derived peroxidases in digestion and physiology are numerous and diverse, direct roles for insect peroxidases in lignin degradation have not been explored. lignin degradation releases phenylpropanoids, which are often toxic; however, a. glabripennis produces enzymes capable of degrading phenylpropanoid subunits, including epoxide hydrolases, which are often involved in polycyclic aromatic compound metabolism. other transcripts predicted to encode detoxification enzymes and antioxidants that could make contributions to degradation or inactivation of toxic lignin metabolites include alcohol dehydrogenases, aldehyde dehydrogenases, cytochrome p450s, glutathione s-transferases, catalases, carboxylesterases, enzymes involved in aromatic compound degradation, and glucuronosyl transferases. additionally, aldo-keto reductases are capable of degrading phenolic compounds, including tannins and phenylpropanoids released from lignin degradation, and could be primed for detoxification roles. transcripts predicted to encode detoxification enzymes a. glabripennis eggs hatch directly beneath the bark of hardwood trees and first and second instars feed on primary phloem and xylem, which serve as diffuse transport systems for toxic tree defensive compounds, before tunneling into the heartwood as later instars. though heartwood is not as metabolically active as the primary phloem and xylem, it accumulates potentially toxic secondary metabolites, including alkaloids, tannins, hydroxycinnamic acids, and phenolic glycosides, defending the plant from herbivory and protecting structural polysaccharides and biopolymers from biotic assaults. given that a. glabripennis completes development in over different tree species and that it feeds in the phloem and xylem before eventually making its way into the heartwood, this insect must have mechanisms to detoxify or sequester the breadth of defensive plant secondary metabolites it encounters throughout its life cycle. the gut represents the first line of defense against ingested host plant allelochemicals, pesticides, and other toxins and many transcripts predicted to encode detoxification enzymes were detected. for example, cytochrome p450-like unigenes were detected in the a. glabripennis midgut transcriptome. these enzymes have versatile oxidoreductive properties, are highly involved in degrading lipophilic toxins, and have been shown to confer resistance to pesticides as well as small aromatic toxins that can accumulate to high concentrations in the heartwood of trees. the majority of the cytochrome p unigenes detected in the a. glabripennis midgut transcriptome had highest scoring blastp alignments to cytochrome p450s identified in the tribolium castaneum genome, although the percent similarity at the amino acid level ranged from 34% to 66%, reflecting a relatively large degree of divergence from previously annotated cytochrome p450s.table cytochrome p annotations from a. glabripennis midgut transcriptome t. castaneum t. castaneum l. decemlineata t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum l. decemlineata l. decemlineata l. decemlineata l. decemlineata l. decemlineata t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum l. decemlineata t. castaneum t. castaneum t. castaneum t. castaneum t. castaneum l. decemlineata t. castaneum t. castaneum cytochrome p unigenes detected in the a. glabripennis midgut were putatively assigned to clans and families based on the annotations of the highest scoring blastp alignments. to conform to cytochrome p classification conventions, only alignments sharing ≥ 40% amino acid identity were used to annotate these unigenes. while a handful of unigenes predicted to encode mitochondrial cytochrome p450s were identified, non-mitochondrial clans and were more highly represented in comparison. non-mitochondrial cytochrome p450s were classified to five families, which were predominated by cyp, but also included cyp, cyp, cyp, and cyp. cyp family genes are often present in multiple copies, occur in clusters in insect genomes, have pivotal roles in the detoxification of host plant defensive chemicals such as xanthotoxins, gossypol, and chlorogenic acid, and are often induced in herbivorous insects during periods of feeding. unigenes predicted to encode carboxylesterases were the most dominant enzyme with putative involvement in detoxification processes detected in the midgut; individual unigenes were predicted to encode these proteins. although secreted carboxylesterases are generally involved in pheromone metabolism, intracellular carboxylesterases are often implicated in pesticide and allelochemical metabolism and tolerance. over of the carboxylesterase unigenes detected in the a. glabripennis midgut lacked secretory peptides and may be primed to serve detoxification roles. for example, carboxylesterases are hypothesized to mediate resistance to phenolic glycosides in papilio canadensis, and are often found in high levels in the midgut. notably, trees in the family salicaceae, which include many of a. glabripennis’ preferred hosts, are notorious producers of phenolic glycosides and the abundance of carboxylesterases may promote colonization and survival in these hosts. transcripts predicted to encode enzymes involved in conjugative deactivation of xenobiotic compounds were also detected from and unigenes predicted to encode udp-glucuronosyl transferases and glutathione s-transferases, respectively. these transferases can bind to xenobiotic compounds containing a diversity of functional groups, including oxygen, nitrogen, sulfur or carboxyl groups, enhancing their solubility and allowing them to be excreted or stored in the fat body for eventual elimination. they have been previously shown to detoxify cyanates and cinnamaldehydes, which can be found in high concentrations in heartwood. further, they can also conjugate and eliminate aromatic compounds, including tannins and toxic aromatic compounds stored in the heartwood or released from lignin degradation. transcripts predicted to encode enzymes involved in nitrogen acquisition although nitrogen is scarce in the woody tissue of host trees, like other insects, a. glabripennis larvae have high demands for nitrogen during growth and development. although nitrogen sources are present in very low abundance in woody tissues, microbes associated with the midgut have the metabolic capacity to synthesize all amino acids, which could be assimilated and stored by the insect in the form of arylphorin and hexameric storage proteins encoded by a. glabripennis. however, the insect possesses endogenous pathways for amino acid synthesis that may be complemented or augmented by microbial pathways, including complete pathways for the synthesis of alanine, aspartic acid, asparagine, proline, cysteine, glycine, and serine and for the synthesis of tyrosine from phenylalanine. in addition, nearly complete pathways for the synthesis of arginine, glutamic acid and selenocysteine were detected, but argininosuccinate lyase, glutamate formiminotransferase, and selenocysteine synthases transcripts were absent. these pathways may be incomplete because transcripts encoding these enzymes are simply not expressed in the midgut, they were expressed at low levels and were not detected at the sequencing depth obtained, or because they may be complemented by microbial enzymes that catalyze these reactions. the phloem tissue where early instars feed is rich in amino acids relative to the heartwood where older larvae grow and develop. therefore, recycling waste products of amino acid and nucleotide deamination reactions back into functional amino acids, nucleotides, and other nitrogen-containing compounds may be important to the nitrogen economy in a. glabripennis. a number of transcripts with highest scoring blastp alignments to enzymes predicted to catalyze deamination reactions and liberating ammonia from a variety of nitrogen-containing compounds were detected in the midgut transcriptome. these included adenine deaminases, cytosine deaminases, nitrilases, amidohydrolases, and chitin-degrading enzymes. ammonia liberated from nitrogen-containing compounds could be directly converted to glutamine by glutamine synthetase and aspartic acid by aspartate ammonia ligase and indirectly incorporated into the synthesis of purines, glutamate, alanine, asparagine, and proline by enzymes encoded by the insect. furthermore, ammonia or amino acids constructed from recycled ammonia could be shuttled to microbes housed in the midgut to synthesize nonessential or essential amino acids, augmenting or complementing a. glabripennis’ physiological capabilities. despite these potential contributions to nitrogen economy, the mechanisms of essential amino acid synthesis and recycling are not clear as the abundance of essential amino acids in woody tissue varies depending on tree species, but are significantly lower than the abundance of nonessential amino acids. as with other insects studied, no pathways for the synthesis and metabolism of essential amino acids were detected in the midgut transcriptome, although these could be expressed elsewhere. however, full pathways for the synthesis of nine essential amino acids were detected in the a. glabripennis midgut metagenome, including histidine, isoleucine, leucine, lysine, methionine, phenylalanine, threonine, tryptophan, and valine, and could serve as key sources of essential amino acids in this insect. furthermore, the gut community also contained numerous uricase and urease genes, which could be involved in recycling nitrogenous waste products produced by a. glabripennis or its gut microbes. although the waste products from a. glabripennis have not been biochemically characterized, all enzymes associated with the urea cycle were detected in the a. glabripennis midgut, including several arginase transcripts that catalyze the conversion of arginine to urea, suggesting that urea could be produced in this insect. this urea pathway is also functional in the guts of several other insects.figure proposed mechanisms for direct utilization of ammonia detected in the a. glabripennis midgut. several methods for re-incorporating ammonia produced through nucleotide or amino acid deamination reactions into amino acid and nucleosides were reconstructed based on the presence of unigenes in the a. glabripennis midgut transcriptome. ammonia can be directly integrated into glutamine by glutamine synthetase, which can subsequently be converted into glutamate by microbial-derived glutamate synthase detected in the a. glabripennis midgut metagenome. glutamate could subsequently be converted into either proline by a series of reactions catalyzed by γ-glutamyl kinase, dehydrogenase, aminotransferase, and pyrroline-5-carboxylate or aspartate via transaminase, respectively. furthermore, aspartate can be converted into inosine, which can be directly used for synthesis of purine nucleotides or asparagine via asparagine synthetase. *microbial in origin. 1succinylaminoimidazole carboxamide ribotide synthase. 2aminoimidazole carboxamide ribotide transformylase. transcripts predicted to encode several types of digestive proteinases, including serine and cysteine proteinases, were detected, which aid in protein acquisition from glycoproteins cross-linked in the cell wall matrix or from microbes housed in the midgut. by far, serine proteinase unigenes were more expansive than cysteine proteinase unigenes and cystatins were co-expressed and likely repress cysteine proteinase activities in the midgut. these results are consistent with a previous study that reported high serine protease activities and the absence of cysteine and aspartic peptidase activities in the a. glabripennis midgut. however, the expression of both cysteine and aspartic proteinases in the midgut suggests that this insect still has the genetic capacity to produce these proteinases under certain conditions and these genes may be maintained in the genome as a strategy to combat digestive proteinase inhibitors produced by host plants. analysis of transcripts involved in converting compounds in woody tissue into fatty acids and sterols was also conducted. transcripts involved in facilitating interactions with gut microbes although there is debate about how microbes associated with cerambycid guts contribute to digestive physiology, a number of transcripts with putative involvement in mediating interactions with microbes were detected in the a. glabripennis midgut transcriptome. several of these transcripts are likely involved in maintaining host-microbe homeostasis, including transcripts predicted to encode both antifungal and antibacterial proteins, dual oxidases, mucin, which forms a protective barrier to protect the midgut from microbial invasion, mpa allergen proteins, which have antimicrobial properties and are often upregulated during periods of stress, and several encapsulation proteins involved in activating innate immune pathways. seven unigenes predicted to encode hemocyanins were detected. while they primarily function as oxygen carriers in crustaceans, hemocyanins are rare, but not completely absent in insect genomes and their physiological functions are not well characterized. despite their functional obscurity in insects, they can function as pro-phenol oxidases under certain circumstances, activating innate immune pathways and mediating insect-microbe interactions in the midgut. hemocyanins have also been hypothesized to serve roles in the degradation of lignin since transcripts encoding hemocyanins are highly expressed in a symbiont-free, wood-feeding marine isopod, but no direct involvement in this process has been demonstrated. additionally, signal peptides were not observed in any of the hemocyanin unigenes detected in a. glabripennis which, if not an assembly artifact, would preclude their involvement in extracellular digestive processes.figure structure of hemocyanin enzymes detected in the a. glabripennis midgut transcriptome. seven unigenes predicted to encode hemocyanins were detected. unlike many insect hemocyanins, the hemocyanins produced by a. glabripennis had intact and presumably functional copper binding domains. these unigenes were predicted to have involvement in mediating interactions with gut microbes. identification of highly expressed genes transcripts originating from genomic and mitochondrial ribosomal rrnas were omitted from this analysis and the unigenes with the top fpkm values containing predicted coding regions were identified as highly expressed. many of the highly expressed genes identified in the midgut have predicted involvement in stress and immune modulation and included unigenes predicted to encode several mpa allergen domain proteins, five carboxylesterases, two cathespins, two encapsulation related proteins, two mucin proteins, a lysosomal lipase, a cytochrome p a thaumatin domain protein, and a lectin domain protein. while carboxylesterases and cytochrome p450s have key involvements in detoxification and sterol acquisition, cathespins, encapsulation proteins, mucin proteins, lipases, thaumatin domain proteins and lectins are hypothesized to play fundamental roles in mediating host-microbe interactions. not surprisingly, gh and gh cellulases and gh β-glucosidases were also highly expressed in the midgut, reflecting the nutritional importance of cellulose to a. glabripennis. gh and gh β-galactosidases were also highly expressed, suggesting that galactan polymers present in hardwood hemicellulose or on the cell surface of microbes are also crucial sources of sugar for this insect. chitin deacetylase unigenes were highly abundant; these enzymes can liberate acetate from insect or fungal chitin, which can be recycled for energy or fatty acid production. several different types of digestive proteinases were also highly expressed and included m peptidases, m carboxypeptidases, serine proteinases, and cysteine proteinases, which likely serve key roles in nitrogen extraction from plant or microbial cell wall proteins. in addition to these digestive proteins, several unigenes predicted to encode hypothetical proteins with unknown functions were abundant, suggesting that a. glabripennis encodes suites of novel proteins that could be relevant for digestive physiology and development. glycoside hydrolase profile comparisons through comparisons of transcriptome libraries sampled from a variety of herbivorous insects, no major trends were detected with regard to gh profiles and feeding habitats. euclidean distances between insects that fed on similar substrates were large in many cases and reflected strong differences in gh compositions. thus, the gut transcriptome libraries did not show any significant clustering trends by food source. for example, agrilus planipennis and dendroctonus ponderosae both feed in phloem and were found in separate planes of the pca ordination, indicating that there were large differences in gh family composition between these two insects. likewise, a. glabripennis, coptotermes formosanus, and reticulitermes flavipes all feed in wood and were also found in opposite quadrants in the pca ordination. although these insects are all capable of producing endogenous cellulases, a. glabripennis produces different types of cellulases than the wood- and phloem-feeding insects compared in this study. for example, the two termite species included in this analysis predominately produce gh cellulases, while a. glabripennis produces gh gh and gh cellulase transcripts and dendroctonous ponderosae produces gh and gh cellulase transcripts.figure multivariate comparison of glycoside hydrolase families detected in the gut transcriptomes of herbivorous insects. principal components analysis was conducted to plot the abundance of each glycoside hydrolase family detected in the gut transcriptomes of wood-feeding insects. pca values were plotted in sample space and variables are displayed as vectors on the pca biplot. pca axis and pca axis explain % and % of the variation in the data, respectively. monte carlo permutation procedure: p <  for pca and pca no major clustering of insects feeding on similar food sources was detected. despite the lack of clustering by feeding niche, there appeared to be some clustering by phylogenetic relatedness. for example, most cerambycid and chrysomelid beetles were positioned along the positive x axis and, like a. glabripennis, some of these insects produce transcripts predicted to encode gh and cellulases, although they feed on very different parts of their host plants. furthermore, it is interesting to note that gh cellulases have not yet been found in any insect outside the order coleoptera, but the number of gh cellulases unigenes detected in insect species from this order varied tremendously. while gh transcripts were not detected in association with many coleopterans, the chrysomelids gastrophysa viridula and callosobruchus maculatus encode one and four gh unigenes, respectively. phylogenetic analysis of translated proteinsequences revealed that, although chyrsomelid gh cellulases and cerambycid gh cellulases share a common ancestor, chrysomelid cellulases have rapidly diverged from cerambycid cellulases. in contrast, gh cellulases within the cerambycidae seem to have multiplied and diversified through gene conversion or gene duplication events and are possibly more adapted to digesting highly insoluble cellulose associated with woody plants.figure phylogenetic analysis of gh cellulases detected in the coleoptera. unrooted maximum likelihood tree for insect, bacterial, and nematode derived gh family proteins was constructed with garli using the wag  + i + f + g evolutionary model. fully resolved bootstrap consensus trees were compiled using sum trees and bootstrap values >  are displayed. nodes are annotated with ncbi protein database accession numbers. scale bar represents branch lengths. furthering the hypothesis that the pca ordination was primarily driven by phylogenetic relatedness, family-specific trends in abundances of gh families were observed within the coleoptera. in contrast to gh cellulases, which seem to have multiplied in some cerambycid beetles, gh and gh cellulases were expressed as single copy genes in a. glabripennis. in contrast, members of these gh families have multiplied and diversified in the chrysomelids and curculionid lineages, suggesting that coleopterans have undergone lineage specific adaptations to overcome challenges associated with different feeding regimes. for example, the results of the gh maximum likelihood analysis suggest that gh enzymes were likely encoded in the genome of the last common ancestor of coleopterans and that they underwent family-specific adaptations. this scenario is supported since gh proteins in each insect associated family formed their own supported clusters in the maximum likelihood tree. in particular, genes encoding gh enzymes were likely duplicated in the chrysomelidae. all members of this family encode at least two gh proteins and the branching topology suggests that the second gh gene originated directly from the first. likewise, gh genes have also duplicated and proliferated throughout the chrysomelid and curcurlionid lineages, but the dynamics driving the evolution of this gh family seem to be more complex in comparison to the gh family. in some species, gh genes have rapidly propagated and diversified, while in other cases, the insect expressed only a single copy of this gene.figure phylogenetic analysis of gh cellulases detected in the coleoptera. unrooted maximum likelihood tree for insect and bacterial derived gh family proteins was constructed with garli using the lg  + i + f + g evolutionary model. fully resolved bootstrap consensus trees were compiled using sum trees and bootstrap values >  are displayed. nodes are annotated with ncbi protein database accession numbers. scale bar represents branch lengths. phylogenetic analysis of gh cellulases detected in the coleoptera. unrooted maximum likelihood tree for coleopteran-derived gh family proteins with garli using the wag  + g evolutionary model. fully resolved bootstrap consensus trees were compiled using sum and bootstrap values >  are displayed. nodes are annotated with ncbi protein database accession numbers. scale bar represents branch lengths. the hypothesis that a. glabripennis benefits from microbial enzymes to facilitate nutrient acquisition is supported through the above comparisons. for example, transcripts predicted to encode gh family and chitinases, β-hexosamidases, and glucosylceramidases are strongly associated with a. glabripennis. although these chitin-degrading genes may be important for remodeling the gut peritrophic matrix, which is predominately composed of chitin, they could also play key roles in modulating interactions with fungal taxa associated with the midgut, including yeasts and fusarium solani, a soft rot fungal symbiont of a. glabripennis. we hypothesize that the predominance of chitinases allows a. glabripennis to derive a portion of its carbohydrate and nitrogen resources from fungal chitinous cell walls, which are composed of polymers of n-acetylglucosamine. non-entomopathogenic fungi associated with wood-feeding insects have been previously hypothesized to concentrate and/or recycle nitrogen and thus we also hypothesize that the ascomycota fungal strains found in association with the a. glabripennis midgut serve these same purposes. multivariate transcriptome comparisons of gene ontology annotations like the gh analysis, multivariate comparisons of level four go categories in midgut transcriptome libraries sampled from a variety of herbivorous insects revealed no significant clustering of transcriptome libraries by feeding habitat. however, phylogenetic relatedness alone did not explain the observed pattern of clustering achieved for the transcriptome comparisons. in addition, subsets of different go categories were enriched in each transcriptome library included in the comparison, while many of the go categories were present in approximately the same abundances in each library. together, these findings suggest that most insects possess similar repertoires of gene families and that these genes have adapted in lineage specific manners optimal for overcoming digestive and nutritional challenges associated with specific feeding habitats and ecological niches. for example, although most insects produce a similar number of gh unigenes, the gh family-level comparison uggested that each insect produced its own unique gh profile. other go categories that are present in similar abundances in all insects included in this comparison include 4-α-glucanotransferases, heme binding and transporting proteins, and regulatory genes.figure two-way cluster analysis of level four gene ontology terms from herbivorous insect gut transcriptomes. two-way cluster analysis of level four gene ontology terms was performed to identify potential correlations between abundances of level four go terms and feeding niche. clustering did not appear to be driven by food source or phylogenetic relatedness. however, the a. glabripennis midgut transcriptome is notably distinct from the other gut transcriptomes included in this comparison. despite the lack of clustering by food source or phylogenetic relatedness, several trends were detected that distinguish a. glabripennis from the rest of the insects included in this comparison that could be pivotal to its ability to digest lignocellulose and other wood polysaccharides and extract nutrients from a broad range of deciduous host trees. for example, in the midgut of a. glabripennis, more unigenes and transcript isoforms were produced with predicted monooxygenase and oxidoreductase activities relative to other insects included in this comparison, which could be relevant to its ability to detoxify allelochemicals from its broad range of host plants. further examination of pfam domain abundances in each library revealed that unigenes and transcript isoforms predicted to encode carboxylesterases and cytochrome p450s were more abundant in the a. glabripennis midgut than many of the other insect libraries sampled. a. glabripennis also has the broadest host range of any insect included in this comparison, suggesting that it needs to encode a broader arsenal of detoxification enzymes relative to other insects included in this comparison. several unigenes predicted to encode digestive peptidases, ligases, and protein transporters were also overrepresented relative to other insect transcriptome libraries, which could be relevant for digesting and assimilating proteins produced by microbes associated with the midgut or from plant cell walls. these digestive peptidases are also overrepresented in g. viridula, c. formosanus, and m. sexta transcriptome libraries. unigenes associated with hydrolase activity were also highly abundant in a. glabripennis, many of which were predicted to encode atpases and other nucleosidases, dna binding proteins, rna binding proteins, nucleotide binding proteins, and transferases involved in transferring phosphorous containing groups. the high abundance of unigenes for these nucleotide binding proteins and nucleosidases is likely associated with the high numbers of unigenes predicted to encode reverse transcriptases, transposases, and integrases that were detected in the midgut. finally, unigenes predicted to encode proteins with substrate-specific and active transmembrane transporter activities, including major facilitator family transporters, were also highly abundant in the a. glabripennis midgut. mfs transporters are a diverse group of carriers involved in the absorption of small solutes, including sugar, aromatic amino acids, and other small compounds, which may be involved in assimilation and utilization of small microbial metabolites and/or small metabolites released from the degradation of woody tissue. taken together, the differences in gh family and level four go compositions among insects with similar feeding regimes suggest that the ability to degrade polysaccharides found in woody tissue evolved through lineage specific adaptations rather than through convergent evolutionary processes. insights into interactions with microbes in a recent publication documenting the metabolic potential of the microbiota associated with the a. glabripennis midgut, we discovered a taxonomically diverse assemblage of bacteria and fungi primed to make key contributions to digestive physiology, wood digestion, and nutrient acquisition in this system. however, the potential contributions of the beetle to digestive processes were unknown as no comprehensive transcriptomic or genomic resources are currently available for cerambycids; thus, it is difficult to determine precisely how the gut microbes augment or complement physiological processes in the gut. here, we demonstrate that a. glabripennis also has a relatively versatile metabolic potential, producing a number of its own key digestive enzymes. in particular, a. glabripennis endogenously produces a diverse suite of detoxification enzymes, including cyp 450s and carboxylesterases, many of which act promiscuously on a diverse array of toxins and are hypothesized to facilitate the degradation and inactivation of host tree defensive compounds. in contrast, comparatively few cyp 450s and carboxylesterases were detected in the microbiome of the beetle gut. instead, the community predominantly produces enzymes that target specific classes of toxins, including salicylates and other compounds involved in salicylic acid mediated defenses, as well as arsenic, cyanide, and cyanoamino acids. these could enhance endogenous detoxification processes driven by beetle-derived enzymes and could also serve to detoxify many of the defensive chemicals that act directly on the microbial community, thereby having indirect deleterious impacts on the beetle. furthermore, a. glabripennis produces several different types of plant cell wall degrading enzymes, including cellulases, β-glucosidases, and xylanases that likely serve a major role in breaking down these prominent wood polysaccharides. however, the community also produces a diverse array of plant cell wall degrading enzymes, which could not only serve to augment the cell wall degrading capacities of a. glabripennis, but could also complement these endogenous processes. for example, while the beetle produces enzymes predicted to target β- linkages in xylan, several genes predicted to encode enzymes with activity towards ester linkages in xylan were detected in association with the gut community and included acetyl xylan and feruloyl esterases. the presence of these esterases could expedite degradation of xylan in the midgut by cleaving ester linkages that cross-link adjacent hemicellulose chains, exposing sugar residues to beetle-derived hydrolytic enzymes and enhancing digestibility. in addition, although xylose is the second most abundant sugar in deciduous trees, no putative β-xylosidases were detected in the a. glabripennis midgut transcriptome. however, gene tags predicted to encode β-xylosidases and enzymes associated with the pentose phosphate pathway were highly abundant in the midgut community, which could possibly serve to convert xylose sugars into compounds that can be directly used by a. glabripennis for energy and fatty acid production. in addition, the potential for cooperation between a. glabripennis and its gut microbes was also noted with regards to lignin degradation and nitrogen recycling. for example, a. glabripennis produces a small pool of transcripts that could facilitate degradation of lignin, including aldo-keto reductases, laccases, and peroxidases. alone, these enzymes can only facilitate small-scale degradation of the phenolic linkages that comprise a relatively small percentage of the lignin biopolymer. however, larger scale lignin degrading reactions could be accomplished through interactions with enzymes produced by the gut community. for example, the microbial community contained a high abundance of gene tags involved in synthesizing aromatic redox mediators, which could work in tandem with insect-derived laccases to facilitate degradation of the β-aryl ether linkages that predominantly comprise lignin, and genes tags predicted to directly degrade β-aryl ethers and other abundant linkages in lignin. the beetle requires essential nutrients to complete its development, which are lacking in woody tissue. complete pathways for the synthesis of all major amino acids and several essential vitamins were detected in the gut community, which could augment the production of non-essential amino acids endogenously synthesized by the beetle and contribute to the production of essential amino acids and other essential nutrients, which the beetle cannot synthesize. although a. glabripennis has abilities to scavenge ammonia, reincorporating it into nonessential amino acids and nucleotides, the community has an expanded capacity to recycle nitrogenous waste products including urea, uric acid, xanthine, and arginine, to potentially reincorporate ammonia into both essential and non-essential amino acids, nucleotides, and other nitrogen-containing compounds. thus, the community could serve as an additional source of nitrogen and non-esseential amino acids in this high c:n environment. the microbiota also has the capacity to fix atmospheric nitrogen, providing additional sources of nitrogen to both the beetle and members of the gut community. while this suggests that the beetle collaborates with its gut microbes to facilitate survival in woody tissue, it is unknown which of these microbial pathways are metabolically active in the a. glabripennis midgut. because the transcriptome library was sampled primarily from midgut tissue, few microbial transcripts were detected in this dataset. to gain further insight into transcriptional activity of microbes in the gut, we are currently sequencing rna collected from the midgut contents to prepare a sample that is more enriched in microbial rnas. this was the same approach that was used to sequence the metagenome of the a. glabripennis midgut and we anticipate that this will expand our abilities to more conclusively model the interactions between a. glabripennis and its gut microbes that enhance fitness and/or are required for survival in woody tissue.  the a. glabripennis midgut transcriptome provides the first comprehensive insight into the endogenous digestive capabilities of wood-boring cerambycid larvae as they feed in a highly lignified and nutritionally deficient environment. comparative transcriptome analysis clearly distinguished the a. glabripennis midgut transcriptome from the gut transcriptome libraries of other herbivorous insects that have been previously sampled for sequencing, which may contribute to its long life cycle and ability to feed and develop in a highly lignified food source. our results highlighted gene categories that were enriched in the a. glabripennis midgut transcriptome, which were hypothesized to make key contributions to this insect’s lifestyle, including its ability to colonize a broad range of living host trees. for example, unigenes predicted to encode monooxygenases, carboxylesterases, heat shock proteins, and other detoxification enzymes were highly abundant in the midgut relative to unigenes expressed in the gut transcriptomes of many other herbivorous insects included in this comparison. furthermore, a. glabripennis expressed its own unique profile of gh unigenes for liberating sugars from woody tissue. our results also highlight deficiencies in endogenous digestive and metabolic pathways that could be supplied by microbes associated with the gut, including enzymes for xylose and pentose sugar utilization, enzymes that facilitate the degradation of lignin, pathways for the synthesis of essential amino acids and nutrients, and pathways for nutrient scavenging. for these reasons, we hypothesize that enzymes derived from midgut-associated microbes complement the expression of insect genes and serve vital roles in the digestive physiology of a. glabripennis. this study also substantially expands the genetic resources available for coleopterans by providing transcriptome data for an insect that feeds in the heartwood of healthy host trees. the results presented here may serve as a source for bioprospecting of novel enzymes to enhance industrial biofuels productions and for the pursuit of novel targets for controlling this serious pest. 
a. glabripennis produces a number of enzymes with putative roles in cell wall digestion, detoxification, and nutrient extraction, which likely contribute to its ability to thrive in a broad range of host trees. while a large diversity of insect-derived transcripts predicted to encode digestive and detoxification enzymes were detected, few transcripts predicted to encode enzymes required for lignin degradation or synthesis of essential nutrients were identified, suggesting that collaboration with microbial enzymes may be required for survival in woody tissue. in addition, transcripts predicted to encode enzymes involved in detoxification were detected, including a substantial number of unigenes classified as cytochrome p450s and carboxylesterases, which are hypothesized to play pivotal roles in detoxifying host tree defensive chemicals and could make important contributions to a. glabripennis’ expansive host range. the most highly-expressed genes predicted to encode digestive-related enzymes were trypsins, carboxylesterases, β-glucosidases, and cytochrome p450s. transcripts from putative protein coding regions of over insect-derived genes were identified in the a. glabripennis midgut transcriptome using a combination of shotgun and illumina paired-end reads. here, we present the midgut transcriptome of larval anoplophora glabripennis, a wood-boring beetle with documented lignin-, cellulose-, and hemicellulose- degrading capabilities, which provides valuable insights into how this insect overcomes challenges associated with feeding in woody tissue. wood-feeding insects often work in collaboration with microbial symbionts to degrade lignin biopolymers and release glucose and other fermentable sugars from recalcitrant plant cell wall carbohydrates, including cellulose and hemicellulose.
 wood-feeding insects often work in collaboration with microbial symbionts to degrade lignin biopolymers and release glucose and other fermentable sugars from recalcitrant plant cell wall carbohydrates, including cellulose and hemicellulose. here, we present the midgut transcriptome of larval anoplophora glabripennis, a wood-boring beetle with documented lignin-, cellulose-, and hemicellulose- degrading capabilities, which provides valuable insights into how this insect overcomes challenges associated with feeding in woody tissue.  transcripts from putative protein coding regions of over insect-derived genes were identified in the a. glabripennis midgut transcriptome using a combination of shotgun and illumina paired-end reads. the most highly-expressed genes predicted to encode digestive-related enzymes were trypsins, carboxylesterases, β-glucosidases, and cytochrome p450s. furthermore, unigenes predicted to encode glycoside hydrolases were identified and included several gh and cellulases, gh xylanases, and gh β-glucosidases. in addition, transcripts predicted to encode enzymes involved in detoxification were detected, including a substantial number of unigenes classified as cytochrome p450s and carboxylesterases, which are hypothesized to play pivotal roles in detoxifying host tree defensive chemicals and could make important contributions to a. glabripennis’ expansive host range. while a large diversity of insect-derived transcripts predicted to encode digestive and detoxification enzymes were detected, few transcripts predicted to encode enzymes required for lignin degradation or synthesis of essential nutrients were identified, suggesting that collaboration with microbial enzymes may be required for survival in woody tissue.  a. glabripennis produces a number of enzymes with putative roles in cell wall digestion, detoxification, and nutrient extraction, which likely contribute to its ability to thrive in a broad range of host trees. this system is quite different from the previously characterized termite fermentation system and provides new opportunities to discover enzymes that could be exploited for cellulosic ethanol biofuel production or the development of novel methods to control wood-boring pests. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords cellulasecarboxylesterasecytochrome p450xylanasecomparative transcriptomicshemicellulosebiofuelsbeta-glucosidaseissue-copyright-statement© biomed central ltd 2013 