- Research article
- Open Access
Functional genomics of mountain pine beetle (Dendroctonus ponderosae) midguts and fat bodies
© Aw et al; licensee BioMed Central Ltd. 2010
- Received: 18 December 2009
- Accepted: 30 March 2010
- Published: 30 March 2010
The mountain pine beetle (Dendroctonus ponderosae) is a significant coniferous forest pest in western North America. It relies on aggregation pheromones to colonize hosts. Its three major pheromone components, trans-verbenol, exo-brevicomin, and frontalin, are thought to arise via different metabolic pathways, but the enzymes involved have not been identified or characterized. We produced ESTs from male and female midguts and associated fat bodies and used custom oligonucleotide microarrays to study gene expression patterns and thereby made preliminary identification of pheromone-biosynthetic genes.
Clones from two un-normalized cDNA libraries were directionally sequenced from the 5' end to yield 11,775 ESTs following sequence cleansing. The average read length was 550 nt. The ESTs clustered into 1,201 contigs and 2,833 singlets (4,034 tentative unique genes). The ESTs are broadly distributed among GO functional groups, suggesting they reflect a broad spectrum of the transcriptome. Among the most represented genes are representatives of sugar-digesting enzymes and members of an apparently Scolytid-specific gene family of unknown function. Custom NimbleGen 4-plex arrays representing the 4,034 tentative unique genes were queried with RNA from eleven different biological states representing larvae, pupae, and midguts and associated fat bodies of unfed or fed adults. Quantitative (Real-Time) RT-PCR (qRT-PCR) experiments confirmed that the microarray data accurately reflect expression levels in the different samples. Candidate genes encoding enzymes involved in terminal steps of biosynthetic pathways for exo-brevicomin and frontalin were tentatively identified.
These EST and microarray data are the first publicly-available functional genomics resources for this devastating forestry pest.
- Juvenile Hormone
- Pheromone Component
- Aggregation Pheromone
Understanding the metabolic pathways, their regulation, and the enzymes responsible for creating these chemicals is important in developing potentially targeted management strategies for this devastating forest pest and other destructive Dendroctonus spp. including the spruce beetle (D. rufipennis) and the southern pine beetle (D. frontalis). The three chemicals appear to arise via different metabolic pathways: trans-verbenol is likely the cytochrome P450-mediated hydroxylation product of ingested host tree α-pinene [5, 6]; exo-brevicomin is produced from a fatty acyl-derived precursor [7, 8]; and frontalin is synthesized de novo via the mevalonate pathway . While common enzymes synthesizing their precursors are relatively easily identified, those creating the final pheromone components are more difficult to recognize.
Fortunately, advances in genomic and proteomic technologies allow creation of sequence and/or gene expression databases to investigate economically important but genetically inconvenient "non-model" species. For example, functional genomics studies of the pine engraver beetle, Ips pini (Coleoptera, Scolytidae) [10–12], assisted the identification of important pheromone-biosynthetic genes that otherwise would have been difficult to discover [13–15] and led to subsequent understanding of their regulation and biochemical roles [[13, 16, 17], Figueroa Teran et al., unpublished]. Similarly, a recent effort to characterize the proteome of D. frontalis pronota catalogued differences between males and females that may lead to new control measures . We wished to create functional genomics tools for D. ponderosae to facilitate the development of new strategies to mitigate its economic impact. We produced ESTs representing over 4,000 tentative unique genes, and used custom oligonucleotide arrays to study their regulation in eleven different biological states. We have used these data to make preliminary identification of putative pheromone-biosynthetic genes. Furthermore, the ESTs and microarrays are resources for forest insect researchers working to understand bark beetles.
Two cDNA libraries, "MPB" and "DPG," were prepared in pDONR222 (Invitrogen). The MPB library includes cDNAs from midguts and fat bodies of juvenile hormone (JH) III-treated and acetone-treated (control) insects because of the known role for JH III in stimulating pheromone biosynthesis in the Coleoptera , while the insects used for the DPG library were fed or unfed, but not stimulated with applied hormone and thus can be expected to have more normal biological expression patterns. The MPB library was constructed from 4.9 μg poly(A)+RNA, yielding a primary titre of 1.4 × 106 colony forming units (cfu). The average insert size, not including flanking vector sequences, was 1.3 ± 0.1 kb. The DPG library was constructed from 1.6 μg poly(A)+RNA, yielding a primary titre of 1.27 × 107 cfu and an average insert size of 1.2 ± 0.2 kb. In total, 2,867 and 9,594 (12,461 total) templates were sequenced from the MPB and DPG libraries respectively, yielding 12,119 total ESTs that passed the sequencer's quality control threshold (2,776 from MPB and 9,297 from DPG). The average read length for both libraries was 550 nt.
Further cleansing with EGAssembler  removed 344 additional ESTs. The remaining 11,775 ESTs were masked to remove vector and repeat sequences and clustered with CAP3 to yield 1,201 contigs and 2,833 singlets representing 4,034 tentative unique genes (TUGs). Of these, 2,040 returned at least one hit when compared against GenBank nr using BlastX and an E-value minimum cutoff of 10-5. The ESTs and assembled contigs have been deposited in GenBank [GO484341-GO495894] and [EZ114957-EZ116155], respectively.
Genes with the highest representation by number of ESTs
BlastX hit or identity
A major goal of this study was to combine sequence and microarray expression data to help identify pheromone biosynthetic and potential resin detoxification genes. To this end, we chose 11 biological states that spanned the beetle's life history and pheromone component profile (Figure 1B) to query custom oligonucleotide microarrays representing the 4,034 TUGs. Four biological replicate pools were generated for each of the 11 states. Visual inspection of all 44 hybridizations showed no gross spatial variation due to fibers or bubbles. The normalized distributions of expression values for all 44 arrays were very similar, with no apparent outlying arrays. Digestion curves suggested that trends in RNA degradation between 5' and 3' ends of each probe set were very similar. Thus we concluded that the quality of the arrays and the RNA used to query them was high. The array data have been deposited in NCBI [GEO: GPL9118, Sample ID numbers GSM446276 - GSM446319].
Quantitative (Real-Time) RT-PCR information
Primer pair (Forward/Reverse)
Amplicon Length (bp)
Amplification Efficiency (%)
Genes active in the terminal steps of ipsdienol biosynthesis in I. pini share three common expression characteristics: they are induced in pheromone-biosynthetic tissues by juvenile hormone, their basal expression levels are higher in males than in females, and they are coordinately regulated with mevalonate pathway genes that function earlier in the pathway [11, 12]. Combining these data allowed the identification of previously unknown pheromone-biosynthetic genes, including the dual-function geranyldiphosphate synthase/myrcene synthase [15, 16], myrcene hydroxylase , and a novel oxidoreductase (Figueroa Teran et al., unpublished data). Therefore, we clustered the current microarray data to assist in the identification of D. ponderosae pheromone-biosynthetic genes.
Putative pheromone-biosynthetic genes
The considerations noted above support the assertion by Keeling et al.  that microarray data alone are not reliable indicators of a gene's potential role. While the relatively small EST sample size confers a risk that some pheromone-biosynthetic genes were missed, the available data were still useful for preliminary identification of pheromone-biosynthetic genes, particularly those involved in later steps. For example, exo-brevicomin biosynthesis requires the fatty acyl-derived precursor, (Z)-6-nonen-2-one, to be epoxidized to 6,7-epoxynonan-2-one, a reaction that is likely catalyzed by a cytochrome P450 . One cluster of three genes with a "male-enriched" expression profile contains a cytochrome P450 (CYP6CR1), a putative dehydrogenase, and an unknown protein (Figure 6A). Further qRT-PCR analyses indicate that the P450 gene has expression characteristics consistent with exo-brevicomin biosynthesis, suggesting that it may carry out the epoxidation step, while the dehydrogenase has characteristics implying a role in fatty-acyl-CoA metabolism upstream of epoxidation (G. Song et al., unpublished data). Similarly, frontalin biosynthesis requires carbon to be shunted from the mevalonate pathway either by GPPS (a short-chain isoprenyl diphosphate transferase) , or by a dioxygenase active on GPP or longer-chain isoprenoids, or by another, as yet unknown pathway. One gene, represented by contig126, has a best blastx hit to geranylgeranyl diphosphate synthase (GGPPS, a long chain isoprenyl diphosphate transferase) and an expression profile congruent with HMGR and HMGS (Figure 6B, Figure 7) and with frontalin production: it is induced by feeding, with most of the mRNA localizing to the midgut (Figure 7), where frontalin is produced . While the expression profile is consistent with one expected of a putative GPPS, the activity of this enzyme must be determined. The lack of sequence similarity between this gene and I. pini GPPS may perhaps reflect both broad evolutionary divergence and potentially different enzyme activities. I. pini GPPS also functions as myrcene synthase , an activity that one would not expect in D. ponderosae because myrcene is not a predicted intermediate in frontalin biosynthesis. Both CYP6CR1 (putative epoxidase) and contig126 (putative GGPPS) are current subjects of post-genomic experiments to confirm their functions.
We have generated EST and microarray tools for D. ponderosae. The ESTs and microarray data will be useful for researchers working to develop control strategies for this important pest insect. The microarray data suggest genes encoding enzymes in early steps of metabolic pathways implicated in pheromone component biosynthesis are not clearly coordinately regulated. However, candidate enzymes catalyzing late steps of exo-brevicomin and frontalin biosynthesis were identified, and are subjects of ongoing post-genomic research.
Insects and treatments
Dendroctonus ponderosae were obtained from infested lodgepole pine (Pinus contorta) in Little Valley, NV (39°15'00" N, 119°52'30"W). Infested trees were cut into 1-1.5 m long bolts and placed in rearing boxes as per . Adults that emerged from the bolts were sexed and stored at 4°C on moist paper towels for up to two weeks. Larvae, pupae, and teneral adults were removed from beneath the bark and used immediately. For juvenile hormone (JH) treatment, 10 μg (±)JH III (Sigma) in 0.5 μl acetone was applied to the ventral abdomens of females or males (30 insects/group), followed by incubation in the dark at room temperature in 60 ml plastic containers for 8 or 16 hours. Control insects were treated with 0.5 μl acetone and incubated similarly. For fed samples, females or males were placed in ~4 mm holes drilled through the outer bark of P. contorta bolts and held there with wire mesh for the desired feeding period, after which the beetles were removed from under the bark. Unfed females or males were incubated in plastic cups in the dark for the same amount of time.
Beetles were chilled on ice and then dissected in distilled H2O to remove the midguts and associated fat bodies: in our hands it was virtually impossible to isolate midguts without some fat body remaining attached. The midgut/fat body tissues were pooled according to treatment group and stored at -84°C until use.
cDNA library construction
Two cDNA libraries were prepared in pDONR222 (Invitrogen). The first, "MPB," was prepared using RNA isolated from midguts and fat bodies of eight treatment groups: females or males that had been unfed for 24 h, fed (each sex alone) for 24 h, or treated with JH III and incubated for 8 or 16 h. There were five sets of each treatment group, with 30 beetles/set. The second library, "DPG," was prepared from RNA isolated from midguts and fat bodies of six treatment groups: unfed females or males, females that had been fed alone for 6 or 24 h, females that had been fed alone for 24 h and then with males for 24 h ("fed with mates"), and males that were fed 24 h in the presence of females that had been feeding for 24 h. There were four replicates of 20 beetles/replicate for each condition. For both libraries, total RNA was isolated from each sample using RNeasy Plant Mini Kit (Qiagen). Poly(A)+RNA was isolated by twice passing the pooled total RNA over oligo-dT latex beads (Nucleotrap, Clontech) according to the manufacturer's protocol. cDNAs were prepared and directionally cloned into pDONR222 using the CloneMiner cDNA library construction kit (Invitrogen) and stored at -84°C until use. Insert sizes were determined by PCR amplification with vector primers across the inserts of 20 random clones for each library. Primary titres were determined by growth assays of serially-diluted library stocks.
Sequencing and bioinformatics
A fraction of the MPB library was plated on LB agar containing 50 μg/ml kanamycin and 2,880 colonies were manually picked into 100 μl of TB-kanamycin broth in wells of 30 96-well plates and provided to the Nevada Genomics Center for incubation, plasmid preparation, and single-pass sequencing from the 5' ends using the M13 (Forward) vector primer. Sequencing reactions were performed using an ABI Prism 3730 DNA analyzer. For the DPG library, Agencourt Biosciences was contracted to prepare templates and Sanger sequence single-pass 5' reads of 10,000 individual colonies. The Sanger method was chosen for sequencing because it was the most effective technology available; pyrosequencing technology was not developed enough for our purposes at the time this study was done. The sequences from both libraries were trimmed of low confidence signals using the sequencers' software. The trimmed sequences were combined and further cleansed to remove vector sequences and short reads (<100 nt) using the EGAssembler website . The cleansed ESTs were clustered using CAP3  set with the following parameters: o = 40, p = 95, s = 401, e = 12. Contigs and singlets were annotated by BLASTX alignments to GenBank nr using Blast2GO freeware . Multiple sequence alignments were done with the ClustalW2  server at EBI using default settings.
Roche NimbleGen was contracted to fabricate and hybridize 4-plex microarrays based on 4,034 D. ponderosae sequences from EST assembly and four negative control sequences: pDONR222, pBluescript, Hevea brasiliensis 2-C-methyl-D-erythritol 2,4-cyclodiphosphate synthase [GenBank: AB294705], and Pinus halepensis partial mRNA for rubisco large subunit [GenBank: AJ271897]. Each feature was represented by six 60-mer probes distributed through the sequence, replicated in three non-contiguous blocks for each array.
The arrays were queried with cDNA prepared from 20-30 pooled insects representing 11 different biological states (Figure 1B): larvae (Larv), pupae (Pupa), teneral ("pre-emerged") females (Fpe) or males (Mpe), females that were unfed for 24 h (F24u) or 48 h (F48u), females that fed on P. contorta phloem for 24 h (F24f), females fed for 24 h alone and then with males for an additional 24 h (F48f), females that were fed with males for 24 h without "pre-feeding" alone (F24fM), males that were unfed for 24 h (M24u), and males that were allowed to feed for 24 h in the presence of females (M24fF). Total RNA was prepared from whole bodies of immature stages or from isolated midguts and fat bodies of adults as described above, treated with DNase, and precipitated at -84°C until use. The quality of RNA was confirmed with an Agilent 2100 Bioanalyzer and quantity was determined with a NanoDrop spectrometer (Thermo Scientific). Cy3-labeled cDNA was prepared from the RNA and hybridized to the arrays by Roche NimbleGen. Four biological replicates were prepared for each biological state, for a total of 44 hybridizations.
Microarray data processing and analysis
All NimbleGen custom oligonucleotide array images were examined visually for gross spatial variation due to fibers or bubbles. All array data were processed and normalized first by Robust Multi-Array Average (RMA)  using the R package affy . Specifically, expression values were computed by applying the RMA model of probe-specific correction of perfect match probes. The processed probe values were then normalized via quantile normalization, and a median polish was applied to compute one expression measure from all probe values. A visual inspection showed that the normalized distributions of expression values of all 44 arrays were very similar, with no apparent outlying arrays. Curves representing the trends in RNA degradation between the 5' end and the 3' end in each probeset were generated and inspected and all proved very similar.
To ensure strict reproducibility standards, the quadruplicated expression measurements for each biological state were inspected individually. Any set of quadruplicates in which one of the measures exhibited a standard deviation of more than 1.25 (the maximum possible standard deviation for four measures is 1.499), and a coefficient of variation of greater than 0.75 for the quadruplicate set was scrutinized. If one single measure was near 1.5, this indicated that the remaining three measures were very similar, and that the fourth replicate was at its maximum outlying capacity, and thus this one quadruplicate value was removed. This procedure left three replicates within the set of which the mean was used for subsequent analyses. Only 504 replicate sets (1.2% of all replicated sets) had one measure excluded by this rule. Additionally, any replicate sets that continued to exhibit a coefficient of variation of greater than 0.75 were removed entirely. This included only 235 sets of triplicated measures (0.5% of all replicate measurements), and reduced the mean coefficient of variation of all quadruplicates to 0.19. We found that these thresholds allowed us to identify gross outlying individual measurements within a quadruplicate set .
A simple 1-way ANOVA was performed on the normalized data to determine which features on the array were differentially expressed across the eleven states. A multiple testing correction was applied to the p-values of the ANOVA , and any feature with a significant adjusted p-value P < 0.05 was examined further with a Tukey's post-hoc test. Specifically, any feature that showed a significant difference in means across one of pairwise comparisons: M24u-M24fF, Mpe-M24u, Larv-F24u, M24u-Larv, M24fF-F24fM, F24fM-F24f, M24u-F24f, F24u-F24f, M24fF-Larv, Mpe-M24fF, Mpe-Fpe, F48f-F24f, F48u-F24f, Larv-F24fM, M24u-F24u was retained for further analysis.
The retained features (3,573 total) were then subjected to hierarchical clustering using Pearson's correlation coefficient as distance metric and the average agglomeration method. Clustering dendrograms were examined below the 0.15 height threshold, allowing a close inspection of genes clustered at or above a cluster-average correlation coefficient of 0.85.
Nine features for which preliminary qRT-PCR experiments showed clean amplification efficiencies and melting curves were used to confirm the microarray data. First strand cDNA was prepared and primer pairs were selected as described below. For each of the nine features, log2-transformed ratios of expression relative to unfed females (F24u, arbitrarily chosen as a reference sample) were determined for the qRT-PCR data, normalized to the endogenous control gene, Ubiquitin. The ten qRT-PCR ratios were compared to the analogous log-transformed ratios of the microarray data of the other ten biological states using Pearson's correlation coefficient. A linear regression of 10 pairs of log-transformed qRT-PCR/microarray expression ratios was performed for all nine genes (90 pairs, total). A hypothesis test was performed to evaluate the association between PCR and microarray data, using a test statistic based on Pearson's product moment correlation coefficient, and resulted in a p-value P << 0.001. Steps for performing these verifications were identical to those in .
To confirm the microarray data, first strand cDNA template was prepared from RNA isolated from combined midguts and fat bodies of adults, or whole bodies of immature stages using the RNeasy Plant Mini Kit (Qiagen) followed by RNase-free DNase digestion and further purification using the MasterPure RNA Purification kit (Epicentre). There were 20 insects/replicate, except for the pupal samples, for which six individuals were used. For the tissue distribution assays, beetles were dissected in water into head, anterior midgut, posterior midgut, hindgut, fat body (not associated with the anterior midgut), and carcass with the aid of a dissecting microscope essentially as per . Tissues were pooled from five to eight beetles/sample, frozen in N2(l) and stored at -84°C prior to RNA extraction and purification. RNA pellets were washed with 70% ethanol, resuspended in TE, and checked for quantity and integrity as described above. First strand template cDNA was produced from 500 ng RNA for each sample using random primers (Invitrogen) and SuperScript III reverse transcriptase (Invitrogen) according to the manufacturer's protocol.
Primers for qRT-PCR were designed from the EST sequences using Primer Express v 2.0 software (Applied Biosystems). Selected primers (Table 2) were screened for potential primer-dimers and hairpin loop formation using Vector NTI Advance 9 (Invitrogen), and selected primer sets were tested for non-specific amplification by visual inspection of melting curves, and their amplification efficiencies were determined using a relative standard curve method. A survey of candidate normalizing genes using qBase  indicated tubulin and Ubiquitin [GenBank: EZ115790 and EZ115624, respectively] had the most stable expression across all samples (not shown). Relative expression values for all genes were determined using the ΔΔCT method  normalized to tubulin and Ubiquitin. There were four replicates of each biological sample, each containing material pooled from five to eight insects. There were two template preparations per sample with three technical replicates for each PCR reaction.
We thank members of the laboratories of GJB and CT for assistance with collecting beetles, dissections and helpful advice. The US Forest Service and managers of Whittell Forest gave permission to collect beetle-infested trees. This work was supported by the USDA-NRI (2005-35604-16727), NSF (IOS 0642182), and NIH (P20 RR-016464 from the INBRE Program of the National Center for Research Resources), and HATCH (NEV00339) and McIntyre-Stennis (NEV0369) grants from the Nevada Agriculture Experiment Station (NAES).
- Mock KE, Bentz BJ, O'Neill EM, Chong JP, Orwin J, Pfrender ME: Landscape-scale genetic variation in a forest outbreak species, the mountain pine beetle (Dendroctonus ponderosae). Molecular Ecology. 2007, 16: 553-568. 10.1111/j.1365-294X.2006.03158.x.PubMedView ArticleGoogle Scholar
- Keeling CI, Bohlmann J: Genes, enzymes, and chemicals of terpenoid diversity in the constitutive and induced defence of conifers against insects and pathogens. New Phytologist. 2006, 170: 657-675. 10.1111/j.1469-8137.2006.01716.x.PubMedView ArticleGoogle Scholar
- Borden JH, Ryker LC, Chong LJ: Response of the mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Scolytidae), to five semiochemicals in British Columbia lodgepole pine forests. Canadian Journal of Forest Research. 1986, 17: 118-6 128.Google Scholar
- Pureswaran DS, Gries R, Borden JH, Pierce JHD: Dynamics of pheromone production and communication in the mountain pine beetle, Dendroctonus ponderosae Hopkins, and the pine engraver, Ips pini (Say) (Coleoptera: Scolytidae). Chemoecology. 2000, 10: 153-168. 10.1007/PL00001818.View ArticleGoogle Scholar
- Pitman GB, Vite JP, Kinzer GW, Fentiman AF: Specificity of population-aggregating pheromones in Dendroctonus. Journal of Insect Physiology. 1969, 15: 363-366. 10.1016/0022-1910(69)90282-0.View ArticleGoogle Scholar
- Francke W, Bartels J, Meyer H, Schroder F, Kohnle U, Baade E, Vite JP: Semiochemicals From Bark Beetles: New Results, Remarks, and Reflections. Journal of Chemical Ecology. 1995, 21: 1043-1063. 10.1007/BF02033807.PubMedView ArticleGoogle Scholar
- Francke W, Schroder F, Philipp P, Meyer H, Sinnwell V, Gries G: Identification and synthesis of new bicyclic acetals from the mountain pine beetle, Dendroctonus ponderosae Hopkins (Col.: Scol.). Bioorganic and Medicinal Chemistry. 1996, 4: 363-19. 10.1016/0968-0896(96)00013-2.PubMedView ArticleGoogle Scholar
- Vanderwel D, Gries G, Singh SM, Borden JH, Oehlschlager AC: (E)- and (Z)-6-nonen-2-one: biosynthetic precursors of endo - and exo-brevicomin in two bark beetles (Coleoptera: Scolytidae). Journal of Chemical Ecology. 1992, 18: 1389-1404. 10.1007/BF00994364.PubMedView ArticleGoogle Scholar
- Barkawi LS, Francke W, Blomquist GJ, Seybold SJ: Frontalin: de novo biosynthesis of an aggregation pheromone component by Dendroctonus spp. bark beetles (Coleoptera: Scolytidae). Insect Biochemistry and Molecular Biology. 2003, 33: 773-788. 10.1016/S0965-1748(03)00069-9.PubMedView ArticleGoogle Scholar
- Eigenheer AL, Keeling CI, Young S, Tittiger C: Comparison of gene representation in midguts from two phytophagous insects, Bombyx mori and Ips pini, using expressed sequence tags. Gene. 2003, 316C: 127-136. 10.1016/S0378-1119(03)00749-2.View ArticleGoogle Scholar
- Keeling CI, Bearfield JC, Young S, Blomquist GJ, Tittiger C: Effects of juvenile hormone on gene expression in the pheromone-producing midgut of the pine engraver beetle, Ips pini. Insect Molecular Biology. 2006, 15: 207-216. 10.1111/j.1365-2583.2006.00629.x.PubMedView ArticleGoogle Scholar
- Keeling CI, Blomquist GJ, Tittiger C: Coordinated gene expression for pheromone biosynthesis in the pine engraver beetle, Ips pini (Coleoptera: Scolytidae). Naturwissenschaften. 2004, 91: 324-328. 10.1007/s00114-004-0523-y.PubMedGoogle Scholar
- Sandstrom P, Ginzel MD, Bearfield JC, Welch WH, Blomquist GJ, Tittiger C: Myrcene hydroxylases do not determine enantiomeric composition of pheromonal ipsdienol in Ips spp. Journal of Chemical Ecology. 2008, 34: 1584-1592. 10.1007/s10886-008-9563-7.PubMedView ArticleGoogle Scholar
- Sandstrom P, Welch WH, Blomquist GJ, Tittiger C: Functional expression of a bark beetle cytochrome P450 that hydroxylates myrcene to ipsdienol. Insect Biochemistry. 2006, 36: 835-845. 10.1016/j.ibmb.2006.08.004.View ArticleGoogle Scholar
- Gilg AB, Bearfield JC, Tittiger C, Welch WH, Blomquist GJ: Isolation and functional expression of the first animal geranyl diphosphate synthase and its role in bark beetle pheromone biosynthesis. Proceedings of the National Academy of Sciences of the. 2005, 102: 9760-9765. 10.1073/pnas.0503277102.View ArticleGoogle Scholar
- Gilg AB, Tittiger C, Blomquist GJ: Unique animal prenyltransferase with monoterpene synthase activity. Naturwissenschaften. 2009, 96: 731-735. 10.1007/s00114-009-0521-1.PubMedView ArticleGoogle Scholar
- Bearfield JC, Henry AG, Tittiger C, Blomquist GJ, Ginzel MD: Two regulatory mechanisms of monoterpenoid pheromone production in Ips spp. of bark beetles. Journal of Chemical Ecology. 2009, 35: 689-697. 10.1007/s10886-009-9652-2.PubMedView ArticleGoogle Scholar
- Pechanova O, Stone WD, Monroe W, Nebeker TE, Klepzig KD, Yuceer C: Global and comparative protein profiles of the pronotum of the southern pine beetle, Dendroctonus frontalis. Insect Molecular Biology. 2008, 17: 261-277. 10.1111/j.1365-2583.2008.00801.x.PubMedView ArticleGoogle Scholar
- Tillman JA, Seybold SJ, Jurenka RA, Blomquist GJ: Insect pheromones--an overview of biosynthesis and endocrine regulation. Insect Biochemistry and Molecular Biology. 1999, 29: 481-514. 10.1016/S0965-1748(99)00016-8.PubMedView ArticleGoogle Scholar
- Masoudi-Nejad A, Tonomura K, Kawashima S, Moriya Y, Suzuki M, Itoh M, Kanehisa M, Endo T, Goto S: EGAssembler: online bioinformatics service for large-scale processing, clustering and assembling ESTs and genomic DNA fragments. Nucleic Acids Research. 2006, 34: W459-462. 10.1093/nar/gkl066.PubMed CentralPubMedView ArticleGoogle Scholar
- Girard C, Jouanin L: Molecular cloning of cDNAs encoding a range of digestive enzymes from a phytophagous beetle, Phaedon cochleariae. Insect Biochemistry and. 1999, 29: 1129-1142. 10.1016/S0965-1748(99)00104-6.View ArticleGoogle Scholar
- Pauchet Y, Wilkinson P, van Munster M, Augustin S, Pauron D, ffrench-Constant RH: Pyrosequencing of the midgut transcriptome of the poplar leaf beetle Chrysomela tremulae reveals new gene families in Coleoptera. Insect Biochemistry and Molecular Biology. 2009, 39: 403-413. 10.1016/j.ibmb.2009.04.001.PubMedView ArticleGoogle Scholar
- Casals N, Buesa C, Piulachs MD, Cabano J, Marrero PF, Belles X, Hegardt FG: Coordinated expression and activity of 3-hydroxy-3-methylglutaryl coenzyme A synthase and reductase in the fat body of Blattella germanica (L.) during vitellogenesis. Insect Biochemistry and Molecular Biology. 1996, 26: 837-843. 10.1016/S0965-1748(96)00044-6.PubMedView ArticleGoogle Scholar
- Browne LE: An emergence cage and refrigerated collector for wood-boring insects and their associates. Journal of Economic Entomology. 1972, 65: 1499-1501.View ArticleGoogle Scholar
- Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Research. 1999, 9: 868-877. 10.1101/gr.9.9.868.PubMed CentralPubMedView ArticleGoogle Scholar
- Gotz S, Garcia-Gomez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talon M, Dopazo J, Conesa A: High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Research. 2008, 36: 3420-3435. 10.1093/nar/gkn176.PubMed CentralPubMedView ArticleGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404.PubMedView ArticleGoogle Scholar
- Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics (Oxford, England). 2003, 4: 249-264.View ArticleGoogle Scholar
- Gautier L, Cope L, Bolstad BM, Irizarry RA: affy-analysis of Affymetrix GeneChip data at the probe level. Bioinformatics (Oxford, England). 2004, 20: 307-315. 10.1093/bioinformatics/btg405.View ArticleGoogle Scholar
- Miller G, Schlauch K, Tam R, Cortes D, Torres MA, Shulaev V, Dangl JL, Mittler R: The plant NADPH oxidase RBOHD mediates rapid systemic signaling in response to diverse stimuli. Science Signaling. 2009, 2: ra45-10.1126/scisignal.2000448.PubMedGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B (Methodological). 1995, 57: 289-300.Google Scholar
- Cushman JC, Tillett RL, Wood JA, Branco JM, Schlauch KA: Large-scale mRNA expression profiling in the common ice plant, Mesembryanthemum crystallinum, performing C3 photosynthesis and Crassulacean acid metabolism (CAM). Journal of Experimental Botany. 2008, 59: 1875-1894. 10.1093/jxb/ern008.PubMedView ArticleGoogle Scholar
- Hellmans J, Mortier G, De Paepe A, Speleman F, Vandersomple J: qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. BMC Genomics. 2007, 8: R19-10.1186/1471-2164-8-19.View ArticleGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.