Transcriptomic and phylogenetic analysis of Culex pipiens quinquefasciatus for three detoxification gene families
© Yan et al.; licensee BioMed Central Ltd. 2012
Received: 29 June 2012
Accepted: 8 November 2012
Published: 10 November 2012
The genomes of three major mosquito vectors of human diseases, Anopheles gambiae, Aedes aegypti, and Culex pipiens quinquefasciatus, have been previously sequenced. C. p. quinquefasciatus has the largest number of predicted protein-coding genes, which partially results from the expansion of three detoxification gene families: cytochrome P450 monooxygenases (P450), glutathione S-transferases (GST), and carboxyl/cholinesterases (CCE). However, unlike An. gambiae and Ae. aegypti, which have large amounts of gene expression data, C. p. quinquefasciatus has limited transcriptomic resources. Knowledge of complete gene expression information is very important for the exploration of the functions of genes involved in specific biological processes. In the present study, the three detoxification gene families of C. p. quinquefasciatus were analyzed for phylogenetic classification and compared with those of three other dipteran insects. Gene expression during various developmental stages and the differential expression responsible for parathion resistance were profiled using the digital gene expression (DGE) technique.
A total of 302 detoxification genes were found in C. p. quinquefasciatus, including 71 CCE, 196 P450, and 35 cytosolic GST genes. Compared with three other dipteran species, gene expansion in Culex mainly occurred in the CCE and P450 families, where the genes of α-esterases, juvenile hormone esterases, and CYP325 of the CYP4 subfamily showed the most pronounced expansion on the genome. For the five DGE libraries, 3.5-3.8 million raw tags were generated and mapped to 13314 reference genes. Among 302 detoxification genes, 225 (75%) were detected for expression in at least one DGE library. One fourth of the CCE and P450 genes were detected uniquely in one stage, indicating potential developmentally regulated expression. A total of 1511 genes showed different expression levels between a parathion-resistant and a susceptible strain. Fifteen detoxification genes, including 2 CCEs, 6 GSTs, and 7 P450s, were expressed at higher levels in the resistant strain.
The results of the present study provide new insights into the functions and evolution of three detoxification gene families in mosquitoes and comprehensive transcriptomic resources for C. p. quinquefasciatus, which will facilitate the elucidation of molecular mechanisms underlying the different biological characteristics of the three major mosquito vectors.
KeywordsCarboxyl/cholinesterases Cytochrome P450 monooxygenases Glutathione S-transferases Insecticide resistance Gene expansion Gene expression
Mosquitoes are the most important vectors of human diseases. The Culex pipiens complex has a broad geographic distribution and is the vector of the West Nile virus and the Wuchereria bancrofti nematode, which causes filariasis. Over the last several decades, chemical insecticides have been intensively applied to control disease transmission. However, such control is undermined seriously by the increased insecticide resistance of vector mosquitoes. Three gene families are implicated in insecticide metabolism in mosquitoes: cytochrome P450 monooxygenases (P450s) are responsible for pyrethroid resistance, glutathione S-transferases (GSTs) are responsible for DDT resistance, and carboxyl/cholinesterases (CCEs) are responsible for organophosphate and carbamate resistance. Many insect species show rapid expansion and diversification of detoxification genes, as disclosed by their sequenced genomes. The expansion or restriction of detoxification genes likely helps insects adapt to their particular ecological niches and enable them to survive natural and man-made insecticide selection.
The genomes of three major taxonomic mosquitoes, including Anopheles gambiae, Aedes aegypti, and Culex pipiens quinquefasciatus, have been analyzed and released to the public[4–6]. Of the three, Ae. aegypti has the largest genome size (1376 Mb), while C. p. quinquefasciatus has the largest number of predicted protein-coding genes (18883), which is 22% larger than that of Ae. aegypti and 52% larger than that of An. gambiae. The extra number of protein-coding genes partially results from the expansion of its three detoxification gene families. However, unlike An. gambiae and Ae. aegypti, which have large amounts of gene expression data, such as various expressed sequence tag libraries and transcriptomes, C. p. quinquefasciatus has limited gene expression resources, with only several salivary gland transcriptomes currently reported[7, 8]. Knowledge of complete gene expression information is very important for the exploration of the functions of genes involved in specific biological processes and for the discovery of new candidate genes.
In the present study, the three detoxification gene families of C. p. quinquefasciatus were subjected to phylogenetic analysis and compared with those of three other dipteran insects. The CCE and P450 families were found to undergo large gene expansion. Digital gene expression tag profiling (DGE) technology was used to perform a deep transcriptome analysis of C. p. quinquefasciatus during development and in response to organophosphate insecticide selection. The gene expression profiles obtained provide an invaluable resource for the identification of genes involved in the development and insecticide resistance of C. p. quinquefasciatus.
Results and discussion
C. p. quinquefasciatus detoxification gene families
Classification of detoxification gene families in Drosophila melanogaster , Anopheles gambiae , Aedes aegypti , and Culex pipiens quinquefasciatus
C. p. quinquefasciatus
B class (α-esterases)
D class (integument esterases)
E class (β-esterases)
F class (dipteran JH esterases)
G class (lepidopteran JH esterases)
H class (glutactins)
I class (unknown)
J class (acetylcholinesterases)
K class (gliotactins)
L class (neuroligins)
M class (neurotactins)
CYP3 (include CYP6 and CYP9)
The number of CCEs in the neuro/developmental class was relatively conserved among the four dipteran insects. Similar conservation occurs in hymenopteran (Nasonia vitripennis, A. mellifera) and coleopteran (Tribolium castaneum) genomes, which reflects the relatively ancient origins of this class, where all members are catalytically inactive except for the acetylcholinesterases.
The secreted β-esterases (clade E) are comparatively conserved among the three mosquito species: from 2 to 4 β-esterase genes were found in mosquito genomes, largely different from the expansion (11 β-esterase genes) found in N. vitripennis genome. Some members of the β-esterases have well-described functions in other insects, such as E4 and FE4 esterases, which confer OP insecticide resistance in Myzus persicae, and the antennal Apo1PDE esterase of the silkworm Antheraea polyphemus, which degrades sex pheromones. The functions of the β-esterases in Culex need to be further investigated.
Why does C. p. quinquefasciatus have such an abundance of detoxification genes compared to other insect species? Several biological characteristics of mosquitoes may provide clues. The aquatic breeding sites of larvae and pupae contain numerous microorganisms, phenolic products of plant degradation, and pesticides. Adults feed on plant nectars and mammalian blood, which contain some harmful substances, such as heme and plant toxins. As viral pathogen vectors, mosquitoes have to deal with the generation of toxic endogenous compounds and reactive oxygen species during the immune response. But these cannot account for the gene expansion in C. p. quinquefasciatus compared to Anopheles and Aedes species. Perhaps its more polluted larval habitat and more diverse geographic range have exerted a greater selective pressure on C. p. quinquefasciatus so as to produce a larger repertoire of detoxification enzymes.
DGE library sequencing and mapping to genome
DGE sequencing statistics
Total raw tag
Total clean tag
Distinct clean tag
Distinct tag mapping to gene
Distinct unambiguous tag mapping to gene
Unambiguous tag-mapped genes
Distinct tag mapping to genome
Total unknown tag
Distinct unknown tag
GO and KEGG pathway classification of the genes expressed in C. p. quinquefasciatus
Life-stage specific detected genes
Developmental-stage specifically expressed genes of carboxylcholinesterase (CCE), glutathione-S-transferase (GST), and cytochrome P450 monooxygenase (P450)
For gene expansion clusters of detoxification genes, the expression profiles were different among the members. For example, among the six members of one expanded α-esterase cluster (Figure 2), CpipJ_CPIJ007825 was detected for expression in pupae and adults while the other five members were not detected in any stage. For the seven members of another expanded α-esterase cluster (Figure 2), CpipJ_CPIJ016025, CpipJ_CPIJ005694, and CpipJ_CPIJ008749 were not detected for expression; CpipJ_CPIJ000049 and CpipJ_CPIJ000051 were expressed in pupae and adults, CpipJ_CPIJ000050 in larvae and pupae, and CpipJ_CPIJ016026 only in eggs. Similar phenomenon was observed in the three expansion clusters of juvenile hormone esterases (Figure 3). For the five members of cluster, only CpipJ_CPIJ013027 was detected for expression and in eggs. No expression was found in the three members of cluster. For the four members of cluster, three were detected for expression: CpipJ_CPIJ016681 and CpipJ_CPIJ016682 in adults, and CpipJ_CPIJ017763 in larvae and adults. The different expression patterns of these duplicated detoxification genes are probably indicative of their subfunctionalization or retrogression as pseudogenes.
Differentially expressed genes between parathion resistant and susceptible larvae
When the third instar DGE library of the parathion-resistant strain SG was compared with the same stage in the susceptible strain S-lab, a total of 1511 genes showed different expression levels, among which 619 genes had up-regulated expression levels in the SG strain (Additional file8). The most prominent GO functions of these up-regulated genes were endopeptidase or serine-type peptidase activity, such as genes encoding trypsin, chymotrypsin, mast cell protease 2, urokinase-type plasminogen activator, and elastase. However, not all of the differentially expressed genes are responsible for parathion resistance because comparison strains were not selected from the same panmictic population such that genetic background differences could be ruled out.
Detoxification genes up-regulated in parathion resistant larvae of the SG strain
Egg, larva, adult
Larva, pupa, adult
C. p. quinquefasciatus is an important vector that transmits human diseases different from those by An. gambiae and Ae. aegypti. The lack of transcriptomic data available for this species has hampered characterization of the molecular mechanisms underlying the different biological characters of the three major mosquito vectors. The five DGE libraries described in the present study represent a dramatic expansion of the existing transcriptomic sequence available for C. p. quinquefasciatus. This expansion will facilitate the investigation of the fundamental biology of C. p. quinquefasciatus and its pathogenic interactions. In addition, the results of the present study provide new insights into the functions and evolution of the three detoxification gene families of mosquitoes. A larger number of detoxification genes were identified on the genome of C. p. quinquefasciatus compared with three other dipteran insect genomes, representing the widest gene expansion sequenced thus far. Comparative genomic analysis suggested that gene expansion mainly occurs in α-esterases, juvenile hormone esterases, and P450 CYP325. Some detoxification genes were expressed in all developmental stages, while some were developmentally regulated. The expression profiles were different among the members of gene expansion clusters, probably indicative of their subfunctionalization or retrogression as pseudogenes. Fifteen detoxification genes showed the potential to take part in the parathion resistance of Culex, including unexpected P450 genes.
Mosquito strains of C. p. quinquefasciatus used included S-lab, which was OP-susceptible and reared at the laboratory without any contact with insecticides for many years and Shengui (SG), a field population collected in Foshan, Guangdong Province, in 2007 and constantly treated with parathion at the laboratory. The parathion-resistance of SG was 115-fold that of S-lab before use in the DGE analysis. The mosquitoes were maintained at 26°C ± 1°C and a long-day photoperiod (14 h light/10 h darkness cycle). Fifty egg rafts, forty third instar larvae, forty pupae, and forty adults (twenty females and twenty males) of SG and forty third instar larvae of S-lab were collected and frozen at −80°C for further analysis.
Identification and phylogenetic classification of detoxification genes
Sequences encoding GSTs, P450s, and CCEs were identified from the protein set of the C. p. quinquefasciatus whole genome sequencing database at the Broad Institute (http://www.broadinstitute.org/annotation/genome/culex_pipiens) using the HMMER program (http://hmmer.janelia.org/) with the protein domains for CCEs (PF00135), GSTs (PF00043 and PF02798), and P450s (PF00067) as described in the Pfam database. A significance value of 1e-10 was used in the searches for CCEs and P450 and 2e-2 for GSTs. Community annotations and VectorBase were referred to verify the searches. Those candidate genes not supported by the community annotations as CCEs, P450 or GST were not accounted. P450s were named by the P450 nomenclature committee (http://drnelson.uthsc.edu/CytochromeP450.html). Known detoxification genes from An. gambiae, D. melanogaster, and Ae. aegypti were used as references for the phylogenetic classification of the detoxification genes from C. p. quinquefasciatus. Protein sequences were aligned with ClustalW2 at EMBL-EBI (http://www.ebi.ac.uk/Tools/msa/clustalw2). Unrooted distance neighbor-joining trees showing the phylogeny of detoxification gene families were constructed using the pairwise deletion and p-distance functions of Mega 4.0 software. Bootstrap analysis (1000 replicates) was applied to evaluate the internal support of the tree topology.
Pipeline of DGE
Six micrograms of total RNA from each of the above five mosquito samples were isolated using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. Tag library preparation was performed with an Illumina Gene Expression Sample Prep Kit. The raw data (tag sequences and counts) were deposited in the NCBI Sequence Read Archive (SRA) database under submission number SRA049959.
Pipeline of bioinformatics analysis on DGE
Sequencing-received raw image data were transformed by base calling into raw sequence data. Clean tags were obtained after raw sequences were filtered to remove adaptor sequences, empty tags, low quality tags, tags that were too long or too short, and tags with a copy number of 1. The distribution of clean tags was used to evaluate the normality of the whole data. Saturation analysis was performed to determine whether or not the number of detected genes continues to increase when the sequencing amount increases. Pearson correlation analysis of two parallel libraries was performed to evaluate the reliability and operational stability of the experimental results. All clean tags were mapped to C. p. quinquefasciatus whole genome reference sequences and allowed no more than 1 nucleotide mismatch. The number of unambiguous clean tags for each gene was calculated and then normalized to TPM (number of transcripts per million clean tags). When the expression of a gene was not detected, TPM was set to 0.01.
A rigorous custom written algorithm using the method described by Audic et al. was developed to identify differentially expressed genes between two samples. The p value corresponded to the differential gene expression test. False discovery rate (FDR) was used to determine the p value threshold in multiple tests and analyses. FDR ≤ 0.001 and the absolute value of log2Ratio ≥ 1 were used as thresholds to judge the significance of the gene expression difference.
Unigenes matched by clean tags were assigned to Gene Ontology (GO) terms using Blast2GO and canonical pathways in KEGG (Kyoto Encyclopedia of Genes and Genomes). GO or pathway enrichment analysis of the differentially expressed genes was performed based on the algorithm presented by GOstat. The difference between the differentially expressed gene group and the whole gene expression background was represented by a p value, which was approximated by a chi-square test. The Fisher exact test was used when any expected count value was below 5, which will result in inaccurate chi-square test results. Benjamini multiple-testing correction of the p value was done by FDR.
This study was financially supported by the Major State Basic Research Development Program of China (973 Program, Grant No. 2012CB114102), the Important National Science & Technology Specific Project (Grant No. 2012ZX10004219), and the Knowledge Innovation Program of Chinese Academy of Sciences (Grant No. KSCX3-IOZ-1006).
- Liu N, Xu Q, Zhu F, Zhang L: Pyrethroid resistance in mosquitoes. Insect Sci. 2006, 13: 159-166. 10.1111/j.1744-7917.2006.00078.x.View ArticleGoogle Scholar
- Grant DF, Dietze EC, Hammock BD: Glutathione Stransferase isozymes in Aedes aegypti: purification, characterization, and isozyme-specific regulation. Insect Biochem. 1991, 21: 421-433. 10.1016/0020-1790(91)90009-4.View ArticleGoogle Scholar
- Raymond M, Chevillon C, Guillemaud T, Lenormand T, Pasteur N: An overview of the evolution of overproduced esterase in mosquito Culex pipiens. Philos Trans R Soc Lond B. 1998, 353: 1701-1711. 10.1098/rstb.1998.0321.View ArticleGoogle Scholar
- Arensburger P, Megy K, Waterhouse RM, Abrudan J, Amedeo P, Antelo B, Bartholomay L, Bidwell S, Caler E, Camara F, Campbell CL, Campbell KS, Casola C, Castro MT, Chandramouliswaran I, Chapman SB, Christley S, Costas J, Eisenstadt E, Feschotte C, Fraser-Liggett C, Guigo R, Haas B, Hammond M, Hansson BS, Hemingway J, Hill SR, Howarth C, Ignell R, Kennedy RC, et al: Sequencing of Culex quinquefasciatus establishes a platform for mosquito comparative genomics. Science. 2010, 330: 86-88. 10.1126/science.1191864.PubMed CentralView ArticlePubMedGoogle Scholar
- Holt RA, Subramanian GM, Halpern A, Sutton GG, Charlab R, Nusskern DR, Wincker P, Clark AG, Ribeiro JMC, Wides R, Salzberg SL, Loftus B, Yandell M, Majoros WH, Rusch DB, Lai Z, Kraft CL, Abril JF, Anthouard V, Arensburger P, Atkinson PW, Baden H, de Berardinis V, Baldwin D, Benes V, Biedler J, Blass C, Bolanos R, Boscus D, Barnstead M, et al: The genome sequence of the malaria mosquito Anopheles gambiae. Science. 2002, 298: 129-149. 10.1126/science.1076181.View ArticlePubMedGoogle Scholar
- Nene V, Wortman JR, Lawson D, Haas B, Kodira C, Tu Z, Loftus B, Xi Z, Megy K, Grabherr M, Ren Q, Zdobnov EM, Lobo NF, Campbell KS, Brown SE, Bonaldo MF, Zhu J, Sinkins SP, Hogenkamp DG, Amedeo P, Arensburger P, Atkinson PW, Bidwell S, Biedler J, Birney E, Bruggner RV, Costas J, Coy MR, Crabtree J, Crawford M, et al: Genome sequence of Aedes aegypti, a major arbovirus vector. Science. 2007, 316: 1718-1723. 10.1126/science.1138878.View ArticlePubMedGoogle Scholar
- Girard YA, Mayhew GF, Fuchs JF, Li H, Schneider BS, McGee CE, Rocheleau TA, Helmy H, Christensen BM, Higgs S, Bartholomay LC: Transcriptome changes in Culex quinquefasciatus (Diptera: Culicidae) salivary glands during West Nile virus infection. J Med Entomol. 2010, 47: 421-435. 10.1603/ME09249.View ArticlePubMedGoogle Scholar
- Ribeiro JM, Charlab R, Pham VM, Garfield M, Valenzuela JG: An insight into the salivary transcriptome and proteome of the adult female mosquito Culex pipiens quinquefasciatus. Insect Biochem Mol Biol. 2004, 34: 543-563. 10.1016/j.ibmb.2004.02.008.View ArticlePubMedGoogle Scholar
- Oakeshott JG, Johnson RM, Berenbaum MR, Ranson H, Cristino AS, Claudianos C: Metabolic enzymes associated with xenobiotic and chemosensory responses in Nasonia vitripennis. Insect Mol Biol. 2010, 19 (Suppl. 1): 147-163.View ArticlePubMedGoogle Scholar
- Newcomb RD, Campbell PM, Ollis DL, Cheah E, Russell RJ, Oakeshott JG: A single amino acid substitution converts a carboxylesterase to an organophosphorus hydrolase and confers insecticide resistance on a blowfly. Proc Natl Acad Sci USA. 1997, 94: 7464-7468. 10.1073/pnas.94.14.7464.PubMed CentralView ArticlePubMedGoogle Scholar
- Korochkin LI, Ludwig MZ, Poliakova EV, Philinova MR: Some molecular genetic aspects of cellular differentiation in Drosophila. Soviet Sci Rev. 1987, 1: 411-466.Google Scholar
- Pasteur N, Iseki A, Georghiou GP: Genetic and biochemical studies of the highly active esterases A’ and B associated with organophosphate resistance in mosquitoes of the Culex pipiens complex. Biochem Genet. 1981, 19: 909-919. 10.1007/BF00504256.View ArticlePubMedGoogle Scholar
- Field LM, Williamson MS, Moores GD, Devonshire AL: Cloning and analysis of the esterase genes conferring insecticide resistance in the peach-potato aphid, Myzus persicae (Sulzer). Biochem J. 1993, 294: 569-574.PubMed CentralView ArticlePubMedGoogle Scholar
- Ishida Y, Leal WS: Rapid inactivation of a moth pheromone. Proc Natl Acad Sci USA. 2005, 102: 14075-14079. 10.1073/pnas.0505340102.PubMed CentralView ArticlePubMedGoogle Scholar
- David JP, Strode C, Vontas J, Nikou D, Vaughan A, Pignatelli PM, Louis C, Hemingway J, Ranson H: The Anopheles gambiae detoxification chip: a highly specific microarray to study metabolicbased insecticide resistance in malaria vectors. Proc Natl Acad Sci USA. 2005, 102: 4080-4084. 10.1073/pnas.0409348102.PubMed CentralView ArticlePubMedGoogle Scholar
- Shen B, Dong HQ, Tian HS, Ma L, Li XL, Wu GL: Cytochrome P450 genes expressed in the deltamethrin-susceptible and -resistant strains of Culex pipiens pallus. Pestic Biochem Physiol. 2003, 75: 19-26. 10.1016/S0048-3575(03)00014-2.View ArticleGoogle 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
- Li XC, Schuler MA, Berenbaum MR: Molecular mechanisms of metabolic resistance to synthetic and natural xenobiotics. Annu Rev Entomol. 2007, 52: 231-253. 10.1146/annurev.ento.51.110104.151104.View ArticlePubMedGoogle Scholar
- Puinean AM, Foster SP, Oliphant L, Denholm I, Field LM, Millar NS, Williamson MS, Bass C: Amplification of a cytochrome P450 gene is associated with resistance to neonicotinoid insecticides in the aphid Myzus persicae. PLoS Genet. 2010, 6: e1000999-10.1371/journal.pgen.1000999.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
- Pedra JH, 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
- Ranson H, Rossiter L, Ortelli F, Jensen B, Wang X, Roth CW, Collins FH, Hemingway J: Identification of a novel class of insect glutathione S-transferases involved in resistance to DDT in the malaria vector Anopheles gambiae. Biochem J. 2001, 359: 295-304. 10.1042/0264-6021:3590295.PubMed CentralView ArticlePubMedGoogle Scholar
- Motoyama N, Dauterman WC: The role of nonoxidative metabolism in organophosphorus resistance. J Agric Food Chem. 1974, 22: 350-356. 10.1021/jf60193a055.View ArticlePubMedGoogle Scholar
- Sabourault C, Guzov VM, Koener JF, Claudianos C, Plapp FW, Feyereisen R: Overproduction of a P450 that metabolizes diazinon is linked to a loss-of-function in the chromosome 2 ali-esterase (MdalphaE7) gene in resistant house flies. Insect Mol Biol. 2001, 10: 609-618. 10.1046/j.0962-1075.2001.00303.x.View ArticlePubMedGoogle Scholar
- Dunkov BC, Mocelin G, Shotkoski F, ffrench-Constant RH, Feyereisen R: The Drosophila cytochrome p450 gene Cyp6a2: structure, chromosomal localisation, heterologous expression and induction by Phenobarbital. DNA Cell Biol. 1997, 16: 1345-1356. 10.1089/dna.1997.16.1345.View ArticlePubMedGoogle Scholar
- Georghiou GP, Metcalf RL, Gidden FE: Carbamate resistance in mosquitoes: selection of Culex pipiens fatigans Wied. (= Culex quinquefasciatus) for resistance to Baygon. Bull World Health Org. 1966, 35: 691-708.PubMed CentralPubMedGoogle Scholar
- Lawson D, Arensburger P, Atkinson P, Besanksy NJ, Bruggner RV, Butler R, Campbell KS, Christophides GK, Christley S, Dialynas E, Hammond M, Hill CA, Konopinski N, Lobo NF, MacCallum RM, Madey G, Megy K, Meyer J, Redmond S, Severson DW, Stinson E, Topalis P, Birney E, Gelbart WM, Kafatos FC, Louis C, Collins FH: VectorBase: a data resource for invertebrate vector genomics. Nucleic Acids Res. 2009, 37: D583-D587. 10.1093/nar/gkn857.PubMed CentralView ArticlePubMedGoogle Scholar
- Ranson H, Claudianos C, Ortelli F, Abgrall C, Hemingway J, Sharakhova MV: Evolution of supergene families associated with insecticide resistance. Science. 2002, 298: 179-181. 10.1126/science.1076781.View ArticlePubMedGoogle Scholar
- Audic S, Claverie JM: The significance of digital gene expression profiles. Genome Res. 1997, 7: 986-995.PubMedGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol. 1995, 57: 289-300.Google Scholar
- Beissbarth T, Speed TP: GOstat: find statistically overrepresented gene ontologies within a group of genes. Bioinformatics. 2004, 20: 1464-1465. 10.1093/bioinformatics/bth088.View ArticlePubMedGoogle Scholar