|Home | About | Journals | Submit | Contact Us | Français|
Aedes aegypti is the principal vector of dengue and yellow fever viruses. The availability of the sequenced and annotated genome enables genome-wide analyses of gene expression in this mosquito. The large amount of data resulting from these analyses requires efficient cataloguing before it becomes useful as the basis for new insights into gene expression patterns and studies of the underlying molecular mechanisms for generating these patterns.
We provide a publicly-accessible database and data-mining tool, aeGEPUCI, that integrates 1) microarray analyses of sex- and stage-specific gene expression in Ae. aegypti, 2) functional gene annotation, 3) genomic sequence data, and 4) computational sequence analysis tools. The database can be used to identify genes expressed in particular stages and patterns of interest, and to analyze putative cis-regulatory elements (CREs) that may play a role in coordinating these patterns. The database is accessible from the address http://www.aegep.bio.uci.edu.
The combination of gene expression, function and sequence data coupled with integrated sequence analysis tools allows for identification of expression patterns and streamlines the development of CRE predictions and experiments to assess how patterns of expression are coordinated at the molecular level.
The completed sequence of the Ae. aegypti genome  has enhanced the development of novel methods of manipulating vector populations to effect control of disease transmission . In order to further the prospects of such endeavours, we generated and organized data using gene expression microarrays to quantify genome-wide transcription in adult males and females in different developmental stages. Adult male and female mosquitoes feed on sugar obtained mostly from nectar of flowers and honeydew to meet the energy demands of basal metabolism and flight. In addition, female mosquitoes also feed on blood for egg development. Since this behaviour is associated with reproduction and disease transmission, our study explored changes in gene expression following a blood meal. Arrangement of these data into a searchable format has streamlined the elucidation of those genes that are expressed in a stage- and sex-specific/enhanced manner. In addition, by integrating DNA sequence comparison tools with a pattern-finding interface, analyses of putative cis-regulatory elements (CREs) can be performed on sets of genes that share similar patterns of expression. Building upon our foundation of the Anopheles gambiae Gene Expression Profile at UC Irvine, http://www.angaged.bio.uci.edu, [3-5], a study of an African vector of human malaria, we provide here a public database and web-based data-mining tool that combine staged expression microarray data, functional annotation, genomic sequence data, and integrated DNA sequence comparison algorithms to gain insight into gene expression and regulation in Ae. aegypti.
Stage-specific transcriptional signal values were imported from genome-wide microarray analyses of four-day old adult males, four-day old nonblood-fed adult females, and blood-fed adult females at 3, 12, 24, 48, 72, and 96 hours after a bloodmeal. RNA was extracted from whole mosquitoes and hybridized according to standard protocols to custom-designed microarrays that survey 16,222 Ae. aegypti transcripts (Platform GPL10542). Expression signals were normalized for background within chips with the Agilent spatial correction algorithm (gprocessed signals) and analyzed with JMP Genomics software http://www.jmp.com/software/genomics/ and Cyber-T http://cybert.microarray.ics.uci.edu/. Microarray data and detailed experimental protocols have been submitted to the Gene Expression Omnibus under the accession series: GSE22339. A total of 5,081 (32%) transcripts were identified as having sex-specific or preferential expression (p-values < 0.001). Of those genes exhibiting sex-differential expression, 2,557 accumulated transcripts at higher levels in males, while the remaining 2,524 were found preferentially or specifically in females (Figure (Figure1A).1A). Extensive variation in gene expression was observed in blood-fed females. A total of 4,773 transcripts were found to vary (p-values < 0.001) in accumulation during oogenesis in at least one of the analysed experimental time points when compared with nonblood-fed females (Figure (Figure1B).1B). Several transcripts, such as those corresponding to members of the D7-related group (AAEL006424-RA, AAEL007394-RB) and an apyrase (AAEL006347-RA), are expressed preferentially in adult female salivary glands and assist blood ingestion . Accordingly, these were detected as female-enhanced in our dataset. Similarly, transcripts characterized previously as accumulating in female mosquitoes following a blood meal, such as those involved with digestion (Late Trypsin, AAEL013284-RA)  or reproduction (Vitellogenin, AAEL006126-RA) , also were determined in our dataset to be upregulated after blood ingestion, supporting the legitimacy of the dataset presented here.
Functional gene annotation was imported from the AegyXcel database http://exon.niaid.nih.gov/transcriptome.html#aegyxcel to populate aeGEPUCI with keywords and annotation from the ENSEMBL, NCBI non-redundant, GO, PFAM, and SMART databases. Putative promoter sequences were selected as regions 2.0 kilobases (kb) in length immediately adjacent to the 5'-ends of the annotated 5' untranslated regions (UTRs) where available, or 5'-ends of coding sequences using genomic data from Vectorbase.org (Assembly: AaegL1, Oct 2005; Genebuild: VectorBase, Aug 2006; Database version: 55.1d).
The data have been stored as a MySQL relational database that is accessible directly through an Apache web server http://www.apache.org. A web-based data-mining interface is used to manage queries to identify genes that meet specific expression, keyword, and sequence criteria (Figure (Figure2).2). Sequence comparison and motif-finding tools, including MEME  and AlignACE [10,11], are accessible directly through the data-mining interface for analysis of CREs corresponding to a selected gene set.
The data-mining interface is accessible from the main page of the site and allows users to focus on specific genes that satisfy desired criteria based on: 1) stage- and sex-specific expression, 2) annotated keywords, or 3) DNA sequences in promoter, 3' UTR, 5' UTR, or cDNA (Figure (Figure2).2). The data-mining tool allows users to arrive at an increasingly specific gene set by imposing in a stepwise fashion each new filter criterion upon all preceding criteria. After selecting a set of genes of interest, users can access the analysis menu to perform CRE analyses, view expression profiles in batch, or export associated DNA sequences (Figure (Figure3).3). Detailed annotation and expression data for each gene also can be viewed by selecting the gene-identifier link to open the description of a gene entry.
Each gene in the Ae. aegypti genome has a corresponding data page that may be accessed by selecting the gene-identifier link during data retrieval. Gene entry pages display microarray expression values and functional annotation as gathered by AegyXcel from ENSEMBL, NCBI non-redundant, GO, PFAM, and SMART databases (Figure (Figure4).4). A link to Vectorbase http://www.vectorbase.org on each entry page provides access to additional, centralized gene data. Furthermore, any corresponding links to orthologous genes in the Anopheles gambiae Gene Expression Profile at UC Irvine  are provided as a bridge between these two mosquito species.
After identifying gene sets that show similar patterns of expression, users can use the menu to directly execute analyses to identify corresponding putative CREs in the set (Figure (Figure3).3). The interface allows searching for conserved motifs in regulatory domains within UTRs and protein-encoding sequences. Programs currently available for sequence analysis include MEME and AlignACE, which are run using default parameters. MEME allows users to specify the number of motifs for which to search. AlignACE does not require this parameter. Users are allowed to enter their e-mail addresses to receive an e-mail notification when their analysis job has completed and is available for viewing.
The analysis menu also allows users to view the transcription profiles of all genes in a gene set in batch. Resulting graphs print gene-expression profiles for all genes in the set according to developmental stage: sugar-fed four-day old adult males, four-day old nonblood-fed adult females, and blood-fed adult females at 3, 12, 24, 48, 72, and 96 hours after a bloodmeal (Figure (Figure55).
The import gene set link can be used to load a set of gene or transcript identifiers into the data-mining interface for viewing and analysis in batch. This feature allows the application of sequence analysis tools to gene sets built in contexts and with concepts different from the ones utilized in this database. Similarly, gene sets also can be exported from the data-mining interface by using the analysis menu, easing the transfer of the selected genes or sequences for analysis by additional methods.
The aeGEPUCI database has the capacity to house, integrate, and display additional microarray studies that examine gene expression in Ae. aegypti. The Submit Study link provides a form for uploading microarray data and specifications for review and possible integration into the database.
The aeGEPUCI database identifies genes co-expressed in similar patterns and incorporates keyword searching and sequence analysis into one unified data-mining tool. A case study best illustrates the utility of this integration. In this example, we identify genes linked to the complex regulation of phenoloxidase, an enzyme involved in the melanization of invading parasites and micro-organisms as part of invertebrate innate immunity . Specifically, we search for pro-phenoloxidase genes that are highly expressed 24 hours after bloodfeeding. Two filters are used to complete this inquiry (Figure (Figure2).2). The first filter selects genes that contain the keyword "prophenoloxidase" in their functional annotation. Seventy-five of the 16,221 transcripts in the Ae. aegypti genome contain this keyword. Second, a stage-specific filter identifies 3 of these 75 transcripts that show 5-fold up-regulated expression 24 hours after bloodfeeding (BF24h) as compared to nonblood-fed mosquitoes (NBF).
The analysis menu can be used to search the 5'-flanking putative promoter regions of genes in this gene set for conserved DNA sequence motifs. Analysis of the promoter regions of the three prophenoloxidase-related genes using MEME shows the occurrence of multiple conserved DNA sequence motifs organized in similar order in putative promoter regions (Figure (Figure66).
In addition to its current microarray data, aeGEPUCI has been built with the goal of expanding its scope to house, integrate, and display information from additional gene expression studies of Ae. aegypti. This flexibility assures that aeGEPUCI is capable of growing alongside the increasing quantity of data being produced from other studies. By working closely with Vectorbase http://www.vectorbase.org and other laboratories in this way, it is hoped that aeGEPUCI will act as a catalyst in accelerating the study and understanding of gene expression and regulation in this important vector of disease.
The Aedes aegypti Gene Expression Profile at UCI is publicly accessible from the URL: http://www.aegep.bio.uci.edu. Questions and comments are welcomed through the site.
The authors declare that they have no competing interests.
SND designed and implemented the website, database, and wrote the manuscript. OM designed, conducted and analyzed the microarray experiments and wrote the manuscript. MHW analyzed and normalized the microarray data for entry into the database. WAD analyzed data and edited the manuscript. JMCR constructed the AegyXcel database and edited the manuscript. AAJ and GY assisted in the editing of the manuscript. All authors read and approved the final manuscript.
The authors thank Lynn Olson for help in preparing the manuscript. This work was supported by a grant from the National Institutes of Health (AI29746). Because J.M.C.R. is a government employees and this is a government work, the work is in the public domain in the United States. Notwithstanding any other agreements, the NIH reserves the right to provide the work to PubMedCentral for display and use by the public, and PubMedCentral may tag or modify the work consistent with its customary practices. You can establish rights outside of the U.S. subject to a government use license.