RNA-Seq and molecular docking reveal multi-level pesticide resistance in the bed bug
© Mamidala et al; licensee BioMed Central Ltd. 2012
Received: 28 September 2011
Accepted: 6 January 2012
Published: 6 January 2012
Bed bugs (Cimex lectularius) are hematophagous nocturnal parasites of humans that have attained high impact status due to their worldwide resurgence. The sudden and rampant resurgence of C. lectularius has been attributed to numerous factors including frequent international travel, narrower pest management practices, and insecticide resistance.
We performed a next-generation RNA sequencing (RNA-Seq) experiment to find differentially expressed genes between pesticide-resistant (PR) and pesticide-susceptible (PS) strains of C. lectularius. A reference transcriptome database of 51,492 expressed sequence tags (ESTs) was created by combining the databases derived from de novo assembled mRNA-Seq tags (30,404 ESTs) and our previous 454 pyrosequenced database (21,088 ESTs). The two-way GLMseq analysis revealed ~15,000 highly significant differentially expressed ESTs between the PR and PS strains. Among the top 5,000 differentially expressed ESTs, 109 putative defense genes (cuticular proteins, cytochrome P450s, antioxidant genes, ABC transporters, glutathione S-transferases, carboxylesterases and acetyl cholinesterase) involved in penetration resistance and metabolic resistance were identified. Tissue and development-specific expression of P450 CYP3 clan members showed high mRNA levels in the cuticle, Malpighian tubules, and midgut; and in early instar nymphs, respectively. Lastly, molecular modeling and docking of a candidate cytochrome P450 (CYP397A1V2) revealed the flexibility of the deduced protein to metabolize a broad range of insecticide substrates including DDT, deltamethrin, permethrin, and imidacloprid.
We developed significant molecular resources for C. lectularius putatively involved in metabolic resistance as well as those participating in other modes of insecticide resistance. RNA-Seq profiles of PR strains combined with tissue-specific profiles and molecular docking revealed multi-level insecticide resistance in C. lectularius. Future research that is targeted towards RNA interference (RNAi) on the identified metabolic targets such as cytochrome P450s and cuticular proteins could lay the foundation for a better understanding of the genetic basis of insecticide resistance in C. lectularius.
Cimex lectularius (the bed bug), a hematophagous ectoparasite of humans, is now spreading at alarming rates across the globe [1–6]. The bites of these nocturnal blood feeders result in cutaneous manifestations, urticarial reactions, and occasionally anaphylaxis . Further, scratching of bite sites promotes secondary bacterial infections . Though the role of C. lectularius in disease transmission remains unclear, bed bugs are known to carry methicillin-resistant Staphylococcus aureus (MRSA) and vancomycin-resistant Enterococcus faecium (VRE) [9, 10]. C. lectularius associated with MRSA and VRE strains possibly amplify infections in impoverished urban communities .
It has been more than half a century since C. lectularius first showed resistance to DDT [11, 12]. The sudden resurgence of C. lectularius is purportedly due to increased resistance to broad-spectrum insecticides, changes in pest management practices, frequent international travel and passive dispersal (clothing, luggage and second-hand furniture) [13, 14]. In the recent past, C. lectularius severely affected the hospitality industry, wherein some hotels were closed due to heavy infestations . The lack of effective pest management tools for this blood-feeding insect has led to its successful establishment around the globe .
Insects develop resistance to insecticides through four modes: penetration resistance (thicker cuticle for decreased entry of insecticides), behavioral resistance (avoidance of the toxic compounds), target site resistance (knockdown resistance, kdr), and metabolic resistance (detoxification primarily through the action of cytochrome P450s, glutathione S-transferases [GSTs], and carboxylesterases) . Among these four modes, target site resistance (kdr) has been well characterized in C. lectularius [17, 18], but little to no knowledge exists regarding the other modes of resistance. Recent transcriptomic studies of C. lectularius have set the foundation for understanding the potential contribution of cytochrome P450s (metabolic resistance) in pesticide resistance [19, 20]. However, a more comprehensive and/or global understanding of the genes involved in pesticide resistance and their "modus operandi" in C. lectularius are necessary for understanding the genetic factors establishing resistance in C. lectularius and for improving existing control strategies or devising new ones. The next-generation sequencing (NGS) methods via different platforms (Illumina Genome Analyzer, Applied Biosystems SOLiD, Helicos Biosciences Heliscope, Roche 454 Life Sciences) have revolutionized functional genomics research in non-model organisms [21–24]. Illumina deep sequencing (RNA sequencing, RNA-Seq) has emerged as a powerful tool for simultaneous transcriptome characterization and differential gene expression (DGE) analysis to better understand eco-physiological adaptations of insects [24–28].
In this study, we used RNA-Seq to compare and contrast pesticide-resistant (PR) and pesticide-susceptible (PS) strains of C. lectularius. As a result, we have significantly enriched the existing transcriptomic database of C. lectularius  and subsequently identified candidate genes putatively involved in insecticide resistance. Lastly, molecular docking was performed for a candidate CYP3 clan member (CYP397A1V2) to further determine the likely contribution of P450 mediated insecticide resistance in C. lectularius.
Results and Discussion
De novo assembly
Summary statistics of Illumina and Roche 454 GS FLX reads of Cimex lectularius
Total number of reads
Number of reads after trimming for quality a
Number of reads used for assembly
Number of contigs
Longest contig (bp)
Number of reads mapped
Percentage of reads mapped to the FA
The read variation number observed in the current study has been commonly associated with the Illumina platform as different lanes can produce different number of counts [25, 27]. Other possible reasons for read variation include different percentages of poly A+ in different samples, differential efficiency of conversion of RNA to cDNA, size or GC composition variation among libraries (smaller and/or lower GC content amplicons amplify more efficiently during PCR for library quantitation). All these factors may contribute to slight variations in the estimation of the library concentration that is applied onto the flow-cell, finally resulting in differences in cluster densities. The coverage and number of reads (i.e., depth of sequencing) obtained in this study far exceeds that of another recently published article on the C. lectularius transcriptome involving differential expression between susceptible and resistant populations .
To determine the completeness of our final assembly, each of the derived assembled sequences was compared with its putative ortholog from Acyrthosiphon pisum, Anopheles gambiae, Drosophila melanogaster, and Pediculus humanus corporis. Nearly, 33-35% of the assembled sequences exceeded a 0.8 ortholog hit ratio and 48-50% had more than a 0.5 ortholog hit ratio (Figure 2B). An ortholog ratio of one indicates a fully assembled transcriptome, while a value close to zero indicates a poor assembly . The current transcriptome characterization resulted in a 23-fold enrichment of the existing EST database of C. lectularius (1,971 ESTs http://www.ncbi.nlm.nih.gov accessed on 25th August 2011) which allowed fine characterization of the C. lectularius transcriptome as has been observed in many non-model insect species [24–28].
A total of 15,540 out of the 51,492 ESTs (reference sequences) were annotated and Gene Ontology (GO) terms were found to be distributed in a wide variety of functional categories (Additional file 3). Annotation with the D. melanogaster genome revealed a similar distribution of GO categories with C. lectularius showing no notable bias towards any category. The majority of the ESTs annotated with biological processes were involved in development and biological regulation, while the ESTs annotated with molecular function revealed high catalytic and binding activity (Additional file 3). A number of putative viral sequences were represented in the GO terms for biological processes (viral life cycle) and cellular component (viroids) (Additional file 3 and Additional file 4). Though viral transmission by C. lectularius is still not clear, the Hepatitis B virus (HBV) has frequently been detected in wild C. lectularius . Future studies are required to both validate these viral sequences and determine their biological relevance within the C. lectularius genome. The assigned metabolic pathways revealed a high number of sequences (1,557/8,159) to be involved in synthesis of carbohydrates, proteins, lipids and nucleotides (Additional file 5). These predicted pathways together with gene annotation will better help in revealing gene function in C. lectularius .
Differential gene expression
RNA-Seq libraries in the current study were constructed from paired biological replicates of resistant and susceptible RNA samples of C. lectularius. The four RNA samples were processed in two batches: Batch-1 (Resistant-1, Susceptible-1) and Batch-2 (Resistant-2, Susceptible -2). Hence, we applied GLMseq (a custom R-script for fitting a two-way Generalized Linear Model) analysis to the four un-normalized RNA-Seq counts at each EST. To find differentially expressed ESTs, Illumina counts of all four populations (PS1, PS2, PR1and PR2, refer to the methods section for more details) were aligned to the assembled reference database which resulted in ~15,000 highly significant differentially expressed ESTs (P < 0.005). All differentially expressed ESTs, along with their Absolute Log2 Fold Change Adjusted (ALFCA) and their description based on BLAST hits, are presented in Additional file 6.
The primary objective of GLMseq analysis was to determine the amount and significance of the main resistance effects, setting aside any secondary batch effects due to transcript-specific interactions between batch effects and resistance effects. Currently available methods for RNA-Seq data analysis, such as DESeq, are limited to one-way analysis of replicate groups and consequently cannot correctly handle two-way interactions. Therefore, a custom R-script (GLMseq) was used to fit a two-way generalized linear model consisting of a main resistance effect and a secondary interaction effect.
GO analysis of differentially expressed ESTs
Summary of differentially expressed genes involved in insecticide resistance of Cimex lectularius
Mode of resistance
Number in the top 5,000 differentials
BB_Contig_103, 18015, 19601, 22399
BB_Contig_1762,1766, 17694, 21630, 48951
The enriched GO terms for actin binding, structural constitution of cuticle, and structural constituent of muscle in PR strains of C. lectularius could imply the strengthening of structural components such as the cuticle and midgut-associated structures. Indeed among the up-regulated cluster of ESTs, we found a high number (46) of transcripts encoding for cuticular proteins (Table 2). Cuticular proteins are major components of insect cuticle. The constitutive or induced expression of such proteins potentially establishes penetration resistance in C. lectularius, as observed in other insect systems . It is interesting to observe that only 8 transcripts encoding cuticular proteins were found in the down-regulated differential cluster in PR strains of C. lectularius (Table 2).
For the down-regulated cluster of differentially expressed ESTs, the most enriched GO terms included metalloexopeptidase, aminopeptidase activities, manganese ion binding, lipid transporter activity, etc. (Figure 4B). These trends could imply the transcriptomic adjustments within C. lectularius upon encountering insecticides. While Adelman et al.  attribute multi-level insecticide resistance of bed bugs primarily toward the activity of metabolic genes (cytochrome P450, carboxylesterases and GSTs) plus kdr mutations, their analysis failed to identify major players involved in penetration resistance (e.g., cuticular proteins). Further, in addition to the metabolic genes identified in their study, our unbiased/global approach identified other key players of metabolic resistance such as ABC transporters and antioxidant genes (Table 2).
We also screened the ESTs that were uniquely expressed in susceptible and resistant populations of C. lectularius, wherein a high number of ESTs (77) were found in PS strains and a few (7 ESTs) in PR strains (Additional File 9). In both cases, a majority of the ESTs were TUFs (85.71% in PR and 74.02% in PS strains). At the current time, it is too speculative to explain the uniquely expressed ESTs encoding for cathepsin in PR strains and RP45 and keratin-associated proteins among PS strains (Additional file 9).
RNA-Seq validation of candidate genes
To validate the expression profiles obtained through GLMseq analysis, quantitative real time RT-PCR (qRT-PCR) was performed on 12 selected candidate genes (cuticular proteins, CYP3 and mitochondrial P450 clan members, ABC transporter, superoxide dismutase, and acetylcholinesterase) belonging to the top differential cluster (Additional file 10). Except for one of the candidate genes profiled (mitochondrial CYP member, CYP 301A2), the qRT-PCR results correlated with the GLMseq profiles (Additional file 10). The expression patterns of CYP301A2 of C. lectularius are in corroboration with DDT-resistant and DDT-susceptible fruit flies (D. melanogaster), wherein no appreciable differences of CYP301A1 transcript levels were reported between the susceptible and resistant flies. Further, DDT treatment did not induce the expression levels of CYP301A1 . The mitochondrial P450 clan members are exclusively found in animals with conserved and diversified groups potentially involved in physiological and detoxification processes in insects [41, 42]. Though these clan members have been well documented to detoxify insecticides and plant allelochemicals [42–44], the CYP301A2 of C. lectularius might plausibly also participate in the synthesis of 20-hydroxyecdysone (20-HE) as observed in other insect studies [45, 46].
The PR strains of C. lectularius displayed higher transcript levels for all of the assayed cuticular proteins (larval cuticle protein [LCP], pupal cuticle protein [PCP], chitin synthase [CHS], chitin deacetylase [CDA] and cuticular protein analogous to peritrophin [CPAP]) and are consistent with the RNA-Seq profiles (Additional file 10). The higher expression of LCP and PCP in PR strains could contribute to increased pyrethroid resistance . In addition to their role in other physiological processes, insects commonly up-regulate transcript levels of cuticular proteins to reduce the penetration rate of insecticides [39, 47, 48]. Higher transcript levels for CDA and CHS could contribute altered chitin and chitosan ratios in cuticle, which influences insect survival . Co-expression of the CHS and CDA transcripts might potentially reduce the entry of xenobiotic compounds into the insect body. The higher expression of CPAP in PR strains (Additional file 10) indicates a possible role of the peritrophic membrane (gut) in sequestering, detoxifying ingested xenobiotics in addition to its role in peritrophic membrane formation and protection from invasive parasites [50, 51].
All three cytochrome P450s of CYP3 clan members (CYP397A1V2, putative CYP6A2 and CYP6A13) revealed higher mRNA levels in the PR strains (Additional file 10). This is in agreement with other insect systems showing resistance to DDT or pyrethroids . The transcript levels of a superoxide dismutase and an ABC transporter were also found to be higher in the PR strains (Additional file 10), which suggests that these proteins are potentially involved in the elimination and efflux of intracellular toxins thereby reducing their interaction with intracellular targets [35, 53]. The transcripts encoding acetylcholinesterase (insecticide target protein) showed significantly higher expression in the PR strains (Additional file 10) possibly resulting in increased production of acetylcholinesterase to minimize the toxic effects of insecticides . The current qRT-PCR profiles of three P450s and an ABC transporter together with our previous findings of a Phase II protein (GST)  suggest a phase-wise detoxification of xenobiotic compounds in C. lectularius. Further, the qRT-PCR profiles also validate our GLMseq analysis employed for generating the differentially expressed ESTs.
Tissue and development-specific expression of CYP397A1V2 and CYP6A13
Given our interest in the role of P450s in metabolic resistance, we further profiled the transcript levels of CYP397A1V2 and CYP6A13 in different tissues (cuticle, midgut and Malpighian tubules) of adults and different developmental stages (early and late instar nymphs and adults) of the PR and PS strains. Both genes showed lesser transcript levels in PS strains, therefore the PS strains were used as a calibrator to determine fold change in tissues and developmental stages of the PR samples.
The qRT-PCR analysis of CYP397A1V2 and CYP6A13 showed higher mRNA levels in all developmental stages of PR compared to PS strains, specifically peak levels in early instar nymphs for CYP397A1V2 and late instar nymphs for CYP6A13 (Figure 5B). Our findings correlate with those in other insect systems and therefore suggest the up-regulation of P450s in C. lectularius to efficiently detoxify a broad range of toxic substrates at different tissue interfaces and during development [19, 42, 56–58].
The expression of CYP397A1V2 transcripts in adults are in agreement with recent transcriptomic studies of C. lectularius . However, our results clearly demonstrate the broader appearance of key metabolic players of both penetration resistance as well as metabolic resistance compared to Adelman et al. .
DDT and pyrethroids (permethrin and deltamethrin) are preferential sodium channel modulators in insects . The frequent exposures to DDT, an early generation insecticide during the last century may have altered the development of resistance in C. lectularius to the next generation insecticides through a phenomenon called cross resistance or multi-level resistance . The latter scenario might be the likely phenomenon associated with modern C. lectularius. Our docking effort on CYP397A1V2 supports this hypothesis. Similar observations were found in Papilio glaucus, wherein Pc CYP6B4 had the ability to metabolize a broad range of substrates . The unlikely binding of diazinon with CYP397A1V2 suggests potential involvement of other P450s or metabolic proteins that interact with polar compounds.
Validation for kdr mutation
We also examined the PR and PS strains for mutations in the voltage-gated sodium ion channels, as this is a well-characterized kdr mechanism found in C. lectularius . All three PR strains showed both mutations (V419L and L925I) (Additional file 15). Although the PS1 and PS2 strains did not carry any of these mutations, the PS3 strain showed a single mutation at V419L (Additional file 15). These results further suggest that resistance in the PR strains could be attributed to multiple mechanisms of resistance including mutations as well as penetration and metabolic resistance.
In summary, we have developed significant molecular resources for C. lectularius and have identified several candidate genes potentially involved in different phases of insecticide metabolism as well as those participating in other modes of insecticide resistance in this species. Our GLMseq approach further revealed significant differentially expressed ESTs across biologically variant samples. Specifically, the high occurrence of up-regulated cuticular proteins and the expression patterns of P450s in cuticular tissue might represent unique sites for penetration resistance as well as metabolic resistance in C. lectularius. Molecular modeling and docking studies revealed the potential of P450s to metabolize multiple insecticides in C. lectularius. Future functional studies (RNA interference) on cuticular and P450 proteins could lay the foundation for identifying hot-spots for insecticide resistance in C. lectularius, which could provide the basis for developing effective management strategies.
Insect material for RNA-Seq
Six populations (3 pesticide-susceptible strains--Harlan 1, Harlan 2, and FV strain hereinafter referred to as PS1, PS2 and PS3; 3 pesticide-resistant strains collected during 2010 from three residences in different zip codes in Columbus, OH, hereinafter referred to as PR1, PR2 and PR3) of C. lectularius were used in the current study. Bed bug populations were maintained under ambient conditions in the laboratory (~22 ± 1.22°C and 35 ± 6% RH). For each bed bug population, all stages and representatives of multiple generations were housed together on filter paper strips contained in a glass jar (13 cm high by 7 cm diameter) (narrow-mouth Mason pint jar, Ball Corp., Broomfield, CO) with an organza fabric and filter paper covering held in place using a screw-on ring. Multiple jars were needed for large populations. Each bed bug population was fed in situ approximately every 2 weeks on heparinized chicken blood using the Hemotek 5W1 membrane feeding system for blood-sucking insects (Discovery Workshops, Accrington, United Kingdom) with Parafilm® as the membrane.
The resistance status of bed bug populations was assessed using a discriminating dose adapted from Romero et al. . The Harlan strain was included as a standard susceptible strain given its long-term laboratory rearing status (since 1973). For this test, 10 adult bed bugs from each test population were placed on filter paper discs treated at a rate of 0.13 mg/cm2 with technical grade deltamethrin (99% purity, Chem Service, Westchester, PA) dissolved in acetone. Acetone-treated filter paper was used as a control. Each population was replicated in triplicate. The condition of the bugs assessed after 24 h exposure to the treated filter paper. At 24 h, 100% of susceptible bed bugs (Harlan strains [PS1, PS2]) were either moribund or dead after exposure to technical deltamethrin, whereas resistant bed bugs (PR1, PR2, PR3) did not show any signs of intoxication. Preliminary resistance testing with PS3 indicated that despite having a known kdr mutation (V419L) (Additional file 15), it was succeptible to dried residues of both lambda-cyhalothrin and deltamethrin at or near the label rate.
RNA isolation and Illumina Paired End (PE) library preparation
Total RNA was extracted separately from each individual (8 adults per population of PS1, PS2, PR1 and PR2) using TRIzol (Invitrogen) and pooled before conducting a quality check (Thermo Scientific Nanodrop 2000 and Agilent Bioanalyzer, Ohio State University, Columbus, OH, USA). Illumina high-throughput sequencing was done using the GAII platform at the Molecular and Cellular Imaging Center, Ohio Agriculture Research and Development Center, Wooster, OH, USA. Samples were prepared for Illumina GA sequencing using the PE library preparation kit (Catalogue Number PE-102-1001, Illumina, San Diego, CA, USA) and TruSeq RNA Sample Prep Kits (Catalogue Number FC-122-1001, Illumina, San Diego, CA, USA) per manufacturer instructions. Briefly, Sera-mag Magnetic Oligo (dT) beads were used for the poly (A) RNA enrichment and divalent cations were used to fragment the purified mRNA (100-400 bp) by heating the reaction mixture at 94°C for 5 minutes. The fragmented RNA was used for double strand synthesis followed by cDNA synthesis. The DNA was end-repaired and phosphorylated using T4 DNA polymerase, Klenow DNA polymerase and T4 PNK. These fragments were 3' adenylated using Klenow Exo- (3' to 5' exo minus) and Illumina PE adapters were ligated using DNA Ligase. To select appropriate size (for PE library preparation kit) and to eliminate unligated adapters, the adapter ligated products were purified on a 2% TAE-agarose gel (Certified Low-Range Ultra Agarose, Cat. No., 161-3107, Life Science Research, Hercules, CA, USA). cDNA fragments with 200 ± 25 bp were cut from the gel and purified using Qiagen gel purification kit (Cat. No. 28704, Qiagen, Valencia, CA, USA). For the Truseq RNA sample prep kit, size selection and purification from unligated adapter were done using Agencourt AMPure XP beads (Catalogue Number: A63881 Beckman Coulter Genomics Danvers, MA, USA). To enrich adapter ligated fragments and to add additional sequences necessary for flowcell binding, 15 rounds of PCR were performed using Illumina PE 1.0 and PE 2.0 primers (for the PE library preparation kit). The libraries were validated on an Agilent Technologies 2100 Bioanalyzer using the Agilent DNA 1000 chip kit and quantified using quantitative PCR with PhiX sequencing control as a standard (Cat. No. CT-901-2001 Illumina, San Diego, CA, USA). These libraries were sequenced using an Illumina GA II sequencer.
Raw data analysis and de novo assembly
To eliminate low quality nucleotides, raw Illumina reads were trimmed using a custom Perl script with windowed adaptive trimming (Phred quality threshold of 20 and minimum read length of 20 nt). To do the de novo assembly, the trimmed reads from PR sample 1 and PS sample 1 were separately fed into Rnnotator, automated de novo RNAseq assembly pipeline , to remove duplicates and erroneous sequences. Multiple rounds on velvet assemblies then were performed using different hash values to account for the different sequencing depth for different transcripts . Resulting contigs were merged using Minumus2 from the AMOS package [62, 63]. These were combined with previously published C. lectularius 454 sequences  for better coverage and quality using Minumus2 from the AMOS package to make the final reference contigs.
Validation and annotation of assembled contigs
To validate the unigenes, we used three criteria: BLASTx searches (E value < 10-3) between unigenes and the NCBI non-redundant protein (nr) database http://www.ncbi.nlm.nih.gov; BLASTx searches (E value < 10-3) between unigenes and several insect species with genome sequence information (A. pisum, A. gambiae, D. melanogaster, P. humanus); and prediction of putative open reading frames (ORFs) using http://www.scbi.uma.es/cgi-bin/full-lengther/full-lengther_login.cgi. The ortholog hit ratio was calculated using custom python script . To annotate the unigenes, BLASTx searches (E value < 10-3) were performed between nr and unigenes. BLAST results were imported to Blast2GO program for further annotation of the unigenes . After the mapping step, those gene ontology (GO) terms with E value < 1E-3, annotation cut-off > 45, and GO weight > 5 were used for annotation. To find the pathways in which putative peptides of the unigenes are involved, analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) was performed using Blast2GO . To find enriched GO terms, enrichment analysis (Fisher's Exact Test) tool in Blast2GO software was used with term filter value below 0.05, term filter mode "FDR" and two-tailed test. To categorize the GO terms into different GO categories, a web-based tool, CateGOrizer with "Aqua" tool was used and these categories were compared to a precomputed GO terms from Drosophila melanogaster [67, 68].
Differential gene expression analyses
Differential gene expression analysis among PR and PS samples was performed using two-way GLMseq, a custom R-script for fitting a two-way Generalized Linear Model to the four un-normalized RNA-Seq counts at each contig. The GLM function in R was used to estimate and test both the replicate effect and the resistance effect for each quartet of Poisson counts, with offsets determined by the total mapped reads for each sample. This model-free method makes no assumptions about the four unknown Poisson parameters and requires no pseudo-reference for the calculation of p-values.
The model treated each count as an independently measured Poisson variate with offset determined by library size. The examples below illustrate the effectiveness of GLMseq for setting aside secondary batch effects: In example A (Contig_18664), the raw counts for Batch-1 (4, 4) and Batch-2 (75359, 22) clearly suggest a large interaction effect (Additional file 16). In this case the GLMseq estimate of the main Resistant/Susceptible fold-change was 1.14, with a two-sided p-value of 0.85, indicating the absence of a significant resistance effect. In example B (Contig_27586), the raw counts for Batch-1 (18, 4) and Batch-2 (15041, 527) clearly suggest a positive resistance effect (Additional file 16). In this case the GLMseq estimate of the main Resistant/Susceptible fold-change was 5.15, with a two-sided p-value of 0.003, indicating the presence of a significant positive resistance effect. In both examples, naïve analysis of the pooled Resistant and Susceptible counts would have mistakenly assigned much lower p-values to the estimated fold-changes.
For finding uniquely expressed genes in PR and PS, the read counts were normalized using EDASeq method, which performs normalization for sequence length and GC content for intra and inter samples . We considered ESTs with less than 10 short reads aligned in both biological replicates as not reliably expressed; for the expressed genes, we used more than 100 short reads aligned in both biological replicates as an arbitrary cut off.
cDNA preparation for quantitative PCR analysis
Individuals from the six strains (3 PR and 3 PS) of C. lectularius were categorized into various development stages (early instar nymphs, late instar nymphs and adults) as previously described [19, 70]. The tissues (cuticle, midgut and Malphigian tubules) were dissected as per Mamidala et al. . Total RNA was extracted using TriZOL and evaluated with Nanodrop. Further, RNA was treated with TURBO DNase™ (AMBION, Inc., Austin, TX) to remove any genomic contamination, and stored (at -80°C) until further use. First strand cDNA was synthesized in 20 μL reactions using ~0.5 μg of RNA, oligo dT primer and Super Script III First-Strand Synthesis Super Mix (Invitrogen) and the resultant cDNA was diluted to 20 ng/μl for qRT-PCR studies.
Gene characterization of CYP 397A1V2
To develop full length CYP397A1V2, we performed 3' RACE (Rapid Amplification of cDNA Ends) using an oligo dT primer (Invitrogen) and gene specific primers (Additional file 17). The amplicon was sent for sequencing at the Molecular and Cellular Imaging Center, Wooster, OH, USA. Two additional gene specific primers were designed to the ends of the full length gene to confirm the sequence (Additional file 17). The identification and annotation of the CYP3 clan member was performed using BLASTx search against nr database at NCBI http://www.ncbi.nlm.nih.gov/. For the phylogenetic analysis, an unrooted neighbor-joining tree was constructed with 500 bootstrap replicates and excluding positions with gaps using MEGA version 5 .
Primer designing and qRT-PCR
Twelve potential genes involved in insecticide resistance that showed differential expression in RNA-Seq were profiled including 4 cytochrome P450s (CYP9, CYP6A2, CYP6A13 and CYP301A2), 5 cuticular proteins (larval cuticle protein, pupal cuticle protein, cuticle protein analogous to peritrophin, chitin deacetylase and chitin synthase), an antioxidant gene (superoxide dismutase), an acetylcholinesterase gene, and an ABC transporter. Primers were designed using Beacon Designer 7 software (Additional file 17). Prior to qRT-PCR, standard PCR was performed for all primer pairs and their products were run on agarose gel electrophoresis to ensure single bands. qRT-PCR reactions were performed in 96-well plates using a BioRad thermocycler (CFX-96) as per Bai et al. . The analysis included three biological replicates and three technical replicates (within each biological replicate) for RNA-Seq validation and for developmental stages whereas two biological replicates (PS1, PS2, PR1 and PR2) and two technical replicates (within each biological replicate) were included for tissues. Relative gene expression among PR and PS adult populations was analyzed as per Bai et al. . For tissues and developmental stages, the fold change in gene expression between PR and PS strains of C. lectularius were derived by the 2-ΔΔCT method  using ribosomal protein (RPL18) of C. lectularius as the internal control gene .
The Wilcoxon rank sum test revealed statistically significant differences (P < 0.0001) between the adult samples for PR and PS strains . The transformed CT values (2-CT) were used for statistical analysis using SAS (SAS/STAT User's Guide, Version 9.1, SAS Institute Inc.) with a significance level (α) of 0.05. A PROC MIXED analysis of variance (ANOVA) was performed for tissues and developmental stages of PR and PS strains.
Construction of the CYP397A1V2 model
The MODELLER 9v8 program  was used to construct the CYPP397A1V2 structure. MODELLER is a general program that implements comparative protein structure modeling by satisfying spatial restraints in terms of probability density functions. To this end, the crystal structure of Human Cytochrome P450 3A4 (PDB code 1W0E) was used as a 3D template (Additional file 12). The MODELLER program was applied to generate 40 satisfactory models, including the iron-oxo group of CYP397A1V2. The model with the lowest energy and the lowest restraint violation was selected. The initial model was improved by energy minimization. First, energy minimization of 1,000 steps of steepest descent (SD) followed by 2,000 steps of conjugate gradient (CG) was carried out in order to release the conflicting contacts among residues. After performing 2,000 steps of conjugate gradient (CG) minimization, MD simulation was carried out to examine the quality of the model structure, by checking the stability via performing 2000 ps simulations at a constant temperature 300 K. An explicit solvent model SPC216 water was used . All energy minimization and MD simulations were accomplished by GROMACS4.0 software package [77, 78] using the GROMOS96 43a1 force field [79, 80].
Molecular docking, preparation of the protein and the ligand
The molecular docking program Autodock, which uses a powerful Lamarkian genetic algorithm (LGA), was used to dock the ligands to the protein active site . Insecticide substrate models (DDT, diazinon, deltamethrin, imidacloprid and permethrin) were built, and their geometry was optimized through discover3 in the InsightII/Builder program http://www.accelrys.com. To recognize the binding sites in CYPP397A1V2, blind docking was carried out with the grid size set to 126, 126 and 126 along the X, Y and Z axes with 0.397 Å grid spacing. The AutoDocking parameters used were GA population size: 150 and maximum number of energy evolutions: 250,000. During docking, a maximum number of 10 conformers was considered, and the root-mean-square (rms) cluster tolerance set to 1.5 Å. One of the lowest energy conformations was considered for further analysis.
DNA extraction and gene sequencing for kdr mutations
For each of the bed bug populations (PS1, PS2, PS3, PR1, PR2, PR3), five whole adults were homogenized and genetic material was extracted using the EZNA DNA extraction kit (Omega Bio-Tek). Gene segments containing known mutation sites in the voltage-gated sodium ion channel were amplified using PCR as described in Zhu et al.  (Additional file 17). PCR products were subsequently purified and sequenced by Functional Biosciences (Madison, WI, USA) http://www.functionalbio.com/web/.
All the contigs and CYP397A1V2 were deposited in Genbank under accession numbers GSE31823 and JN624742, respectively.
We thank Joshua Bryant and George Keeney for their help in collection and/or maintenance of bed bugs. We gratefully acknowledge Harold Harlan's contribution of pesticide-susceptible bed bugs. Thanks to Swapna Priya Rajarapu, Binny Bhandary and John Shukle for their help in dissections and qRT-PCR analysis and to Dr. Peter Piermarini for his insightful comments on the manuscript. We also thank Dr. David R Nelson, University of Tennessee, Memphis, TN for his assistance with the nomenclature for CYP397A1V2. This research was supported by State and Federal funds appropriated to the Ohio Agricultural Research and Development Center, The Ohio State University.
- Doggett SL, Geary MJ, Russel R: The resurgence of bedbugs in Australia: with notes on their ecology and control. Env Health. 2004, 4: 30-38.Google Scholar
- Hwang SW, Svoboda TJ, De Jong LJ, Kabasele KJ, Gogosis E: Bed bug infestations in an urban environment. Emerg Infect Dis. 2005, 11: 533-538.PubMed CentralView ArticlePubMedGoogle Scholar
- Harlan JH: Bed bug control: challenging and still evolving. Outlooks on Pest Management. 2007, 18: 57-61. 10.1564/18apr04.View ArticleGoogle Scholar
- Masetti M, Bruschi F: Bedbug infestations recorded in central Italy. Parasitol Int. 2007, 56: 81-83. 10.1016/j.parint.2006.12.002.View ArticlePubMedGoogle Scholar
- Reinhardt K, Harder A, Holland S, Hooper J, Leake-Lyall C: Who knows the bed bug? Knowledge of adult bed bug appearance increases with people's age in three counties of Great Britain. J Med Entomol. 2008, 45: 956-958. 10.1603/0022-2585(2008)45[956:WKTBBK]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Lee IY, Ree HI, An SJ, Linton JA, Yong TS: Reemergence of the bedbug Cimex lectularius in Seoul, Korea. Korean J Parasitol. 2008, 46: 269-271. 10.3347/kjp.2008.46.4.269.PubMed CentralView ArticlePubMedGoogle Scholar
- Delaunay P, Blanc V, Del Giudice P, Levy-Bencheton A, Chosidow O, Marty P, Brouqui P: Bedbugs and infectious diseases. Clin Infect Dis. 2011, 52: 200-210. 10.1093/cid/ciq102.PubMed CentralView ArticlePubMedGoogle Scholar
- CDC: 2011, [http://www.cdc.gov/parasites/bedbugs/health_professionals/index.html]
- Goddard J, de Shazo RD: Bedbugs (Cimex lectularius) and clinical consequences of their bites. JAMA. 2009, 301: 1358-1366. 10.1001/jama.2009.405.View ArticlePubMedGoogle Scholar
- Lowe CF, Romney MG: Bedbugs as vectors for drug-resistant bacteria. Emer Infect Dis. 2011, 17: 1132-1134. 10.3201/eid1706.101978.View ArticleGoogle Scholar
- Busvine JR: Insecticide-resistant strains of insects of public health importance. Trans R Soc Trop Med Hyg. 1957, 51: 11-31. 10.1016/0035-9203(57)90002-0.View ArticlePubMedGoogle Scholar
- Lofgren CS, Keller JC, Burden GS: Resistance tests with the bed bug and evaluation of insecticides for its control. J Econ Entomol. 1958, 51: 241-244.View ArticleGoogle Scholar
- Romero A, Potter MF, Potter DA, Haynes KF: Insecticide resistance in the bed bug: a factor in the pest's sudden resurgence?. J Med Entomol. 2007, 44: 175-178. 10.1603/0022-2585(2007)44[175:IRITBB]2.0.CO;2.PubMedGoogle Scholar
- Boase C: Bed bugs: research and resurgence. Emerging pests and vector-borne diseases in Europe. Edited by: Takken W, Knols BGJ. 2007, Wageningen Academic Publishers, Wageningen, 1: 261-280.Google Scholar
- Anderson A: The decade of bedbugs and fear. Environ Health Insights. 2011, 5: 53-54.PubMed CentralView ArticlePubMedGoogle Scholar
- Mamidala P, Jones SC, Mittapalli O: Metabolic resistance in bed bugs. Insects. 2011, 2: 36-48. 10.3390/insects2010036.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhu F, Wigginton J, Romero A, Moore A, Ferguson K, Palli R, Potter MF, Haynes KF, Palli SR: Widespread distribution of knockdown resistance mutations in the bed bug, Cimex lectularius (Hemiptera: Cimicidae), populations in the United States. Arch Insect Biochem Physiol. 2010, 73: 245-257.PubMedGoogle Scholar
- Yoon KS, Kwon DH, Strycharz JP, Hollingsworth CS, Lee SH, Clark JM: Biochemical and molecular analysis of deltamethrin resistance in the common bed bug (Hemiptera: Cimicidae). J Med Entomol. 2008, 45: 1092-1101. 10.1603/0022-2585(2008)45[1092:BAMAOD]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Bai X, Mamidala P, Rajarapu SP, Jones SC, Mittapalli O: Transcriptomics of the bed bug (Cimex lectularius). PLoS ONE. 2011, 6: e16336-10.1371/journal.pone.0016336.PubMed CentralView ArticlePubMedGoogle Scholar
- Adelman ZN, Kilcullen KA, Koganemaru R, Anderson MAE, Anderson TD, Miller DM: Deep sequencing of pyrethroid-resistant bed bugs reveals multiple mechanisms of resistance within a single population. PLoS ONE. 2011, 6: e26228-10.1371/journal.pone.0026228.PubMed CentralView ArticlePubMedGoogle Scholar
- Li P, Ponnala L, Gandotra N, Wang L, Si Y, Tausta SL, Kebrom TH, Provart N, Patel R, Myers CR, et al: The developmental dynamics of the maize leaf transcriptome. Nat Genet. 2010, 42: 1060-1067. 10.1038/ng.703.View ArticlePubMedGoogle Scholar
- Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009, 10: 57-63. 10.1038/nrg2484.PubMed CentralView ArticlePubMedGoogle Scholar
- Siebert S, Robinson MD, Tintori SC, Goetz F, Helm RR, Smith SA, Shaner N, Haddock SH, Dunn CW: Differential gene expression in the siphonophore Nanomia bijuga (Cnidaria) assessed with multiple next-generation sequencing workflows. PLoS ONE. 2011, 6: e22953-10.1371/journal.pone.0022953.PubMed CentralView ArticlePubMedGoogle Scholar
- Crawford JE, Guelbeogo WM, Sanou A, Traoré A, Vernick KD, Sagnon N, Lazzaro BP: De Novo Transcriptome Sequencing in Anopheles funestus Using Illumina RNA-Seq technology. PLoS ONE. 2010, 5: e14202-10.1371/journal.pone.0014202.PubMed CentralView ArticlePubMedGoogle Scholar
- Bonizzoni M, Dunn WA, Campbell CL, Olson KE, Dimon MT, Marinotti O, James AA: RNA-seq analyses of blood-induced changes in gene expression in the mosquito vector species, Aedes aegypti. BMC Genomics. 2011, 12: 82-10.1186/1471-2164-12-82.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu F, Li W, Li Z, Zhang S, Chen S, Su S: High-abundance mRNAs in Apis mellifera: Comparison between nurses and foragers. J Insect Physiol. 2011, 57: 274-279. 10.1016/j.jinsphys.2010.11.015.View ArticlePubMedGoogle Scholar
- Pitts RJ, Rinker DC, Jones PL, Rokas A, Zwiebel LJ: Transcriptome profiling of chemosensory appendages in the malaria vector Anopheles gambiae reveals tissue- and sex-specific signatures of odor coding. BMC Genomics. 2011, 27: 271-View ArticleGoogle Scholar
- DiGuistini S, Wang Y, Liao NY, Taylor G, Tanguay P, Feau N, Henrissat B, Chan SK, Hesse-Orce U, Alamouti SM, et al: Genome and transcriptome analyses of the mountain pine beetle-fungal symbiont Grosmannia clavigera, a lodgepole pine pathogen. Proc Natl Acad Sci USA. 2011, 108: 2504-2509. 10.1073/pnas.1011289108.PubMed CentralView ArticlePubMedGoogle Scholar
- O'Neil ST, Dzurisin JD, Carmichael RD, Lobo NF, Emrich SJ, Hellmann JJ: Population-level transcriptome sequencing of nonmodel organisms Erynnis propertius and Papilio zelicaon. BMC Genomics. 2010, 11: 310-10.1186/1471-2164-11-310.PubMed CentralView ArticlePubMedGoogle 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 Biochem Mol Biol. 2009, 39: 403-13. 10.1016/j.ibmb.2009.04.001.View ArticlePubMedGoogle Scholar
- Karatolos N, Pauchet Y, Wilkinson P, Chauhan R, Denholm I, Gorman K, Nelson DR, Bass C, et al: Pyrosequencing the transcriptome of the greenhouse whitefly, Trialeurodes vaporariorum reveals multiple transcripts encoding insecticide targets and detoxifying enzymes. BMC Genomics. 2011, 12: 56-10.1186/1471-2164-12-56.PubMed CentralView ArticlePubMedGoogle Scholar
- Pedra JHF, McIntyre LM, Scharf ME, Pittendrigh BR: Genome-wide transcription profile of field- and laboratory-selected dichlorodiphenyltrichloroethane (DDT)-resistant Drosophila. Proc Natl Acad Sci USA. 2004, 101: 7034-7039. 10.1073/pnas.0400580101.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu N, Liu H, Zhu F, Zhang L: Differential expression of genes in pyrethroid resistant and susceptible mosquitoes, Culex quinquefasciatus (S.). Gene. 2007, 394: 61-68. 10.1016/j.gene.2007.01.032.View ArticlePubMedGoogle Scholar
- Labbe R, Caveney S, Donly C: Genetic analysis of the xenobiotic resistance associated ABC gene subfamilies of the Lepidoptera. Insect Mol Biol. 2011, 20: 243-256. 10.1111/j.1365-2583.2010.01064.x.View ArticlePubMedGoogle Scholar
- Mittapalli O, Neal JJ, Shukle RH: Antioxidant defense response in a galling insect. Proc Natl Acad Sci USA. 2007, 104: 1889-1894. 10.1073/pnas.0604722104.PubMed CentralView ArticlePubMedGoogle Scholar
- Scott JG: Cytochromes P450 and insecticide resistance. Insect Biochem Mol Biol. 1999, 29: 757-777. 10.1016/S0965-1748(99)00038-7.View ArticlePubMedGoogle Scholar
- Strode C, Wondji CS, David JP, Hawkes NJ, Lumjuan N, Nelson DR, Drane DR, Karunaratne SH, Hemingway J, Black WC, Ranson H: Genomic analysis of detoxification genes in the mosquito Aedes aegypti. Insect Biochem Mol Biol. 2008, 38: 113-123. 10.1016/j.ibmb.2007.09.007.View ArticlePubMedGoogle Scholar
- Brun-Barale A, Héma O, Martin T, Suraporn S, Audant P, Sezutsu H, Feyereisen R: Multiple P450 genes overexpressed in deltamethrin-resistant strains of Helicoverpa armigera. Pest Manag Sci. 2010, 66: 900-909.PubMedGoogle Scholar
- Awolola TS, Oduola OA, Strode C, Koekemoer LL, Brooke B, Ranson H: Evidence of multiple pyrethroid resistance mechanisms in the malaria vector Anopheles gambiae sensu stricto from Nigeria. Trans R Soc Trop Med Hyg. 2009, 103: 1139-1145. 10.1016/j.trstmh.2008.08.021.View ArticlePubMedGoogle Scholar
- Brandt A, Scharf M, Pedra JHF, Holmes G, Dean A, Kreitman M, Pittendrigh BR: Differential expression and induction of two Drosophila cytochrome P450 genes near the Rst (2) DDT locus. Insect Mol Biol. 2002, 11: 337-341. 10.1046/j.1365-2583.2002.00344.x.View ArticlePubMedGoogle Scholar
- Feyereisen R: Evolution of insect P450. Biochem Soc Trans. 2006, 34: 1252-1255.View ArticlePubMedGoogle Scholar
- Mittapalli O, Bai X, Mamidala P, Rajarapu SP, Bonello P, Herms DA: Tissue specific transcriptomics of the exotic invasive insect pest emerald ash borer. PLoS ONE. 2010, 5: e13708-10.1371/journal.pone.0013708.PubMed CentralView ArticlePubMedGoogle Scholar
- Bogwitz MR, Chung H, Magoc L, Rigby S, Wong W, O'Keefe M, McKenzie JA, Batterham P, et al: CYP12A4 confers lufenuron resistance in a natural population of Drosophila melanogaster. Proc Natl Acad Sci USA. 2005, 102: 12807-12812. 10.1073/pnas.0503709102.PubMed CentralView ArticlePubMedGoogle Scholar
- Guzov VM, Unnithan GC, Chernogolov AA, Feyereisen R: CYP12A1, a mitochondrial cytochrome P450 from the house fly. Arch Biochem Biophys. 1998, 359: 231-240. 10.1006/abbi.1998.0901.View ArticlePubMedGoogle Scholar
- Rewitz KF, Rybczynski R, Warren JT, Gilbert LI: The Halloween genes code for cytochrome P450 enzymes mediating synthesis of the insect moulting hormone. Biochem Soc Trans. 2006, 34: 1256-1260. 10.1042/BST0341256.View ArticlePubMedGoogle Scholar
- Rewitz KF, Gilbert LI: Daphnia Halloween genes that encode cytochrome P450s mediating the synthesis of the arthropod molting hormone: evolutionary implications. BMC Evol Biol. 2008, 8: 60-10.1186/1471-2148-8-60.PubMed CentralView ArticlePubMedGoogle Scholar
- Wood OR, Hanrahan S: Cuticle thickening associated with pyrethroid resistance in the major malaria vector Anopheles funestus. Parasites and Vectors. 2010, 3: 67-10.1186/1756-3305-3-67.PubMed CentralView ArticlePubMedGoogle Scholar
- Wigglesworth B: The physiology of insect cuticle. Ann Rev Entomol. 1957, 2: 37-54. 10.1146/annurev.en.02.010157.000345.View ArticleGoogle Scholar
- Arakane Y, Dixit R, Begum K, Park Y, Specht CA, Merzendorfer H, Kramer KJ, Muthukrishnan S, Beeman RW: Analysis of functions of the chitin deacetylase gene family in Tribolium castaneum. Insect Biochem Mol Biol. 2009, 39: 355-365. 10.1016/j.ibmb.2009.02.002.View ArticlePubMedGoogle Scholar
- Abedi ZH, Brown AWA: Peritrophic membrane as vehicle for DDT and DDE excretion in Aedes aegypti larvae. Ann Entomol Soc Am. 1961, 54: 539-542.View ArticleGoogle Scholar
- Elvin CM, Vuocolo T, Pearson RD, East IJ, Riding GA, Eisemann CH, Tellam RL: Characterization of a major peritrophic membrane protein, Peritrophin-44, from the larvae of Lucilia cuprina cDNA and deduced amino acid sequences. J Biol Chem. 1996, 271: 8925-35. 10.1074/jbc.271.15.8925.View ArticlePubMedGoogle Scholar
- Vontas J, Blass C, Koutsos AC, David JP, Kafatos FC, Louis C, Hemingway J, Christophides GK, Ranson H: Gene expression in insecticide resistant and susceptible Anopheles gambiae strains constitutively or after insecticide exposure. Insect Mol Biol. 2005, 14: 509-521. 10.1111/j.1365-2583.2005.00582.x.View ArticlePubMedGoogle Scholar
- Labbe R, Caveney S, Donly C: Genetic analysis of the xenobiotic resistance-associated ABC gene subfamilies of the Lepidoptera. Insect Mol Biol. 2010, 20: 243-256.View ArticlePubMedGoogle Scholar
- Revuelta L, Piulachs MD, Bellés X, Castañera P, Ortego F, Díaz-Ruíz JR, Hernández-Crespo P, Tenllado F: RNAi of ace1 and ace2 in Blattella germanica reveals their differential contribution to acetylcholinesterase activity and sensitivity to insecticides. Insect Biochem Mol Biol. 2009, 39: 913-919. 10.1016/j.ibmb.2009.11.001.View ArticlePubMedGoogle Scholar
- Chung H, Sztal T, Pasricha S, Sridhar M, Batterham P, Daborn PJ: Characterization of Drosophila melanogaster cytochrome P450 genes. Proc Natl Acad Sci USA. 2009, 106: 5731-5736. 10.1073/pnas.0812141106.PubMed CentralView ArticlePubMedGoogle Scholar
- Aw T, Schlauch K, Keeling CI, Young S, Bearfield JC, Blomquist GJ, Tittiger C: Functional genomics of mountain pine beetle (Dendroctonus ponderosae) midguts and fat bodies. BMC Genomics. 2010, 11: 215-PubMed CentralView ArticlePubMedGoogle Scholar
- Zhu F, Parthasarathy R, Bai H, Woithe K, Kaussmann M, Nauen R, Harrison DA, Palli SR: A brain-specific cytochrome P450 responsible for the majority of deltamethrin resistance in the QTC279 strain of Tribolium castaneum. Proc Natl Acad Sci USA. 2010, 107: 8557-8562. 10.1073/pnas.1000059107.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhou X, Ma C, Li M, Sheng C, Liu H, Qiu X: CYP9A12 and CYP9A17 in the cotton bollworm, Helicoverpa armigera: sequence similarity, expression profile and xenobiotic response. Pest Manag Sci. 2010, 66: 65-73. 10.1002/ps.1832.View ArticlePubMedGoogle Scholar
- Davies TGE, Field LM, Usherwood PNR, Williamson MS: DDT, pyrethrins, pyrethroids, and insect sodium channels. IUBMB Life. 2007, 59: 151-162. 10.1080/15216540701352042.View ArticlePubMedGoogle Scholar
- Romero A, Potter MF, Haynes KF: Insecticide resistant bed bugs: implications for the industry. Pest Control Technol. 2007, 35: 42-50.Google Scholar
- Li W, Schuler MA, Berenbaum MR: Diversification of furanocoumarin metabolizing cytochrome P450 monooxygenases in two papilionids: specificity and substrate encounter rate. Proc Natl Acad Sci USA. 2003, 100: 14593-14598. 10.1073/pnas.1934643100.PubMed CentralView ArticlePubMedGoogle Scholar
- Martin J, Bruno VM, Fang Z, Meng X, Blow M, Zhang T, Sherlock G, Snyder M, Wang Z: Rnnotator: an automated de novo transcriptome assembly pipeline from stranded RNA-Seq reads. BMC Genomics. 2010, 11: 663-10.1186/1471-2164-11-663.PubMed CentralView ArticlePubMedGoogle Scholar
- Sommer DD, Delcher AL, Salzberg SL, Pop M: Minimus: a fast, lightweight genome assembler. BMC Bioinformatics. 2007, 8: 64-10.1186/1471-2105-8-64.PubMed CentralView ArticlePubMedGoogle Scholar
- Ewen-Campen B, Shaner N, Panfilio KA, Suzuki Y, Roth S, Extavour CG: The maternal and early embryonic transcriptome of the milkweed bug Oncopeltus fasciatus. BMC Genomics. 2011, 12: 61-10.1186/1471-2164-12-61.PubMed CentralView ArticlePubMedGoogle Scholar
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.View ArticlePubMedGoogle Scholar
- Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M: KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 1999, 27: 29-34. 10.1093/nar/27.1.29.PubMed CentralView ArticlePubMedGoogle Scholar
- B2G-FAR: A Species Centered GO Annotation Repository. http://bioinfo.cipf.es/b2gfar/showspecies?species=7227
- Hu Z, Jie B, Reecy J: CateGOrizer: A Web-Based Program to Batch Analyze Gene Ontology Classification Categories. Online Journal of Bioinform. 2008, 9: 108-112.Google Scholar
- Risso D, Schwartz K, Sherlock G, Dudoit S: GC-Content Normalization for RNA-Seq Data. 2011, University of California, Berkeley, Division of Biostatistics, Technical report #291, [http://www.bepress.com/ucbbiostat/paper291/]Google Scholar
- Usinger R: Monograph of Cimicidae. Thomas Say Foundation, Vol. 7. College Park, MD: Entomological Society of America. 1966Google Scholar
- Mamidala P, Rajarapu SP, Jones SC, Mittapalli O: Identification and validation of reference genes for quantitative real-time PCR in the bed bug. J Med Entomol. 2011, 48: 947-951. 10.1603/ME10262.View ArticlePubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, doi:10.1093/molbev/msr121,Google Scholar
- Schmittgen TD, Livak KJ: Analyzing real-time PCR data by the comparative C(T) method. Nat Protocols. 2008, 3: 1101-1118. 10.1038/nprot.2008.73.View ArticlePubMedGoogle Scholar
- R Development Core Team: R: A language and environment for statistical computing. R Foundation for Statistical Computing. 2010, Vienna, Austria, ISBN 3-900051-07-0, URL http://www.R-project.orgGoogle Scholar
- Sali A, Blundell TL: Comparative protein modeling by satisfaction of spatial restraints. J Mol Biol. 1993, 234: 779-815. 10.1006/jmbi.1993.1626.View ArticlePubMedGoogle Scholar
- Berendsen HJC, Postma JPM, van Gunsteren WF, Hermans J: Interaction models for water in relation to protein hydration. Intermolecular forces. Edited by: Pullman B. 1981, 331-342.View ArticleGoogle Scholar
- Berendsen HJC, van der Spoel D, van Drunen R: GROMACS: a message-passing parallel molecular dynamics implementation. Comp Phys Comm. 1995, 91: 43-56. 10.1016/0010-4655(95)00042-E.View ArticleGoogle Scholar
- Lindah E, Hess B, van der Spoel D: GROMACS 3.0: a package for molecular simulation and trajectory analysis. J Mol Model. 2001, 7: 306-317.Google Scholar
- Van Gunsteren WF, Billeter SR, Eising AA, Hünenberger PH, Krüger PKHC, Mark AE, Scott WRP, Tironi IG: Biomolecular simulation: The GROMOS96 manual and user guide. 1996, Vdf Hochschulverlag AG, ZürichGoogle Scholar
- Van Gunsteren WF, Daura X, Mark AE: The GROMOS force field. Encycloped Computat Chem. Edited by: Von Rague Schleyer P. 1998, 1211-1216.Google Scholar
- Morris GM, Goodsell DS, Halliday RS, Huey R, Hart WE, Belew RK, Olson AJ: Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function. J Comput Chem. 1998, 19: 1639-1662. 10.1002/(SICI)1096-987X(19981115)19:14<1639::AID-JCC10>3.0.CO;2-B.View ArticleGoogle Scholar