- Short report
- Open Access
Expression and accumulation of the two-domain odorant-binding protein AaegOBP45 in the ovaries of blood-fed Aedes aegypti
Parasites & Vectors volume 6, Article number: 364 (2013)
Aedes aegypti mosquitoes are the main vectors of dengue viruses. Despite global efforts to reduce the prevalence of dengue using integrated vector management strategies, innovative alternatives are necessary to help prevent virus transmission. Detailed characterizations of Ae. aegypti genes and their products provide information about the biology of mosquitoes and may serve as foundations for the design of new vector control methods.
We studied the Ae. aegypti gene, AAEL010714, that encodes a two-domain odorant-binding protein, AaegOBP45. The predicted gene structure and sequence were validated, although single nucleotide polymorphisms were observed. Transcriptional and translational products accumulate in the ovaries of blood fed females and are not detected or are at low abundance in other tissues.
We validated the Ae. aegypti AAEL010714 gene sequence and characterized the expression profile of a two-domain OBP expressed in ovaries. We propose that AaegOBP45 function as a component of the mosquito eggshell.
Aedes aegypti is the main vector of dengue, viruses that causes 390 million infections per year worldwide . Novel methods for the control of this vector to impair dengue transmission are needed urgently. Mosquito genomics and post-genomic research provide new opportunities to explore the biology of mosquitoes and support innovative strategies in vector control .
Odorant-binding proteins (OBPs) are present in a wide range of insect species and experimental evidence shows their involvement in chemoreception [3–6]. However, additional functions have been proposed for these proteins [7–9]. Subgroups of OBPs have been defined according to their amino acid sequences. Atypical OBPs were identified first in Anopheles gambiae, and subsequently in Ae. albopictus, Ae. aegypti and Culex quinquefasciatus. Atypical OBPs contain two amino acid sequence domains similar in sequence to Classic OBPs . Although it was proposed that atypical OBPs are unique to mosquito species [8, 10, 11, 13, 14], structural similarities support their inclusion in a family of proteins called Dimer OBPs, described previously in Drosophila[15, 16]. Accordingly, atypical OBPs have been renamed as two-domain OBP proteins . Aedes aegypti OBP45 (AaegOBP45) is encoded by the gene AAEL010714 and belongs to the matype4 protein cluster . AAEL010714 is up-regulated by blood feeding, attaining a maximum transcript accumulation at 48 hours post blood meal (hPBM), a pattern similar to other Ae. aegypti two-domain OBP genes . Moreover, four OBPs from An. gambiae (OBP11, OBP1, OB44 and OBP13) have a high abundance of their corresponding transcripts also at 48 hPBM and were characterized as mosquito eggshell components . Here we describe the accumulation and distribution of a two-domain odorant-binding protein, AaegOBP45, and its corresponding mRNA.
Aedes aegypti (Higgs white-eye strain) larvae, pupae and adults were reared using standard laboratory procedures .
DNA and RNA extractions and first-strand cDNA synthesis
Genomic DNA extraction was carried out on a pool of 5 mosquitoes using DNeasy Blood & Tissue kit (Qiagen). Whole-body and tissue total RNA preparations were made using TRIZOL (Invitrogen). Samples were extracted from 4th instar larvae, early pupae (0–4 h post-pupation), late pupae (24 h post-pupation), adult males and females (1-day old) and 5-day old females (sugar-fed or 24, 48, 60, 72, 84 and 96 hPBM). Additional samples were derived from 48 hPBM females tissues (fat bodies, ovaries, midguts and heads) dissected in 1x phosphate-buffered saline (PBS). Samples were resuspended in DEPC-treated water and stored at −70°C until required. Total RNA (2.2 μg) was treated with DNaseI (Invitrogen). Approximately 0.2 μg of each treated RNA sample was used as the template for gene amplification with β-actin oligonucleotide primers to confirm the absence of genomic DNA contamination. First-strand cDNA synthesis was done using 2.0 μg of treated total RNA (except from head samples, for which 0.7 μg was used). Oligo(dT)12-18 primer or gene-specific reverse primers and SuperScript III (Invitrogen) were used for RT-PCR.
Oligonucleotide primers designed for AAEL010714 genomic amplification of the corresponding coding region as well as primers for RT-PCR and qRT-PCR experiments, were designed using Primer3 software (http://bioinfo.ut.ee/primer3-0.4.0/) (Additional file 1). The internal forward primer, In-F, was used in combination with the reverse primer to generate a 150 base-pairs (bp) fragment for tissue-specific and qRT-PCR analysis of AAEL010714 mRNA accumulation.
RT-PCRs were performed in final reaction volumes of 25 μL, containing 2.0 mM Mg++, 20.0 mM Tris–HCl (pH 8,4), 50 mM KCl, 0.2 μM of each respective primer, 0.2 mM each dNTP, 2.5 U of Taq DNA polymerase (Invitrogen) and 2 μL of cDNA. The thermocycler program was configured for an initial denaturation step at 94°C for 2 m, followed by 30 cycles (except for 18S ribosomal RNA – 20 cycles) at 94°C for 1 m, 60°C for 30 s and 72°C for 1 m, and a final elongation step at 72°C for 7 m. The annealing temperature was 58°C for the AAEL010714-RA transcript template.
qRT-PCR assay was performed in Mastercycler Realplex 2 thermocycler (Eppendorf) with ABSOLUTE™ QPCR SYBR® Green Mix (ABgene). Amplification conditions were 95°C for 15 m, followed by 45 cycles of 94°C for 30 s, 60°C for 30 s and 72°C for 30 s. Serial dilutions of target plasmids containing the 150 bp cloned cDNA fragments in concentrations ranging from 10-9 g to 10-16 g were used to generate standard curves to obtain the number of transcript copies. Each pool of 10 vitellogenic female samples (36, 48, 60, 72, 84, 96 hPBM) was analyzed in triplicate and the assay was repeated three times with independent biological samples, except for 84 hPBM, which was measured in duplicate. The copy number of β-actin transcripts was used to normalize variation in total cDNA concentration as an endogenous control (Additional file 2) [20, 21]. Statistical analyses were performed using GraphPad Prism® (version 5.00) for Windows (GraphPad Software, San Diego, CA, USA). One-way analyses of variance following Tukey’s Multiple Comparison test were used to estimate statistical significance.
Amplification and cloning procedures for isolation of the AAEL010714 putative gene promoter region and transcription factor binding site prediction
A 2.5 kilobase (kb) region at the 5′-end of the putative translation initiation site of the AAEL010714 gene was identified based on the nucleotide sequence of Supercontig 1.500 (obtained at Vectorbase.org). A primer pair was designed using Primer3 software and included restriction endonuclease cleavage sites for Bam HI (Additional file 1). PCR was performed using genomic DNA as the template and TAQ platinum DNA polymerase (Invitrogen) with the amplification cycles 2 m at 94°C, followed by 35 cycles of 1 m at 94°C, 30 s at 60°C, 2 m at 72°C and a final step of 10 m at 72°C. Amplified products were cloned and sequenced using M13For and M13Rev primers and three internal gene promoter primers (Additional file 1). The consensus DNA sequence was submitted to transcription factor binding sites prediction analyses using the JASPAR database  based on Drosophila melanogaster JASPAR matrix models with a defined 95% profile score threshold. Manual searches were performed for other motifs.
Alignments, commercial synthesis of peptide and antisera production
Nucleotide and protein alignments were performed using BioEdit software . The peptide GNQFSSSDIDGL was used to generate an AaegOBP45 antiserum. The peptide sequence was subjected to a Blastp search against the NCBI database to ensure that it represents a unique sequence (<10% overlap). Specifically, the peptide sequence was selected to prevent cross-reactions of the antiserum with a related protein, AaegOBP44, encoded by AAEL01718. Polyclonal antiserum to AaegOBP45 was generated in rabbits using a standard protocol (Proteimax Biotechnology Ltda., Cotia, São Paulo-Brazil).
Proteins extracted from ovaries dissected at 24, 48 and 72 hPBM (50 μg for each time point) were fractionated by SDS-PAGE on a 12% polyacrylamide gel. Duplicate gels either were stained with 0.1% Coomassie Blue G-250 solution or electroblotted onto Trans- Blot transfer medium- supported nitrocellulose membrane (Bio-Rad). Membrane blocking, antisera incubations and detection procedures were performed following instructions described in the ECL Plus product booklet and detection of immune reacted moieties was carried out using ECL Plus Western blotting reagent (GE Healthcare). Primary antiserum was used at a dilution of 1:500 and HRP-conjugated secondary antibodies at 1:2500.
AAEL01714 gene sequence
Nucleotide polymorphisms are anticipated between mosquito strains  and therefore 3,716 bp of genomic DNA corresponding to the AAEL010714 coding sequence and putative promoter from the Higgs strain were sequenced and compared with the Liverpool reference genome (Vectorbase.org). The genomic sequences showed 99.8% identity and nucleotide polymorphisms were observed only in putative promoter regions.
The putative 5′-end regulatory region of AAEL010714 contains a canonical TATA box and an arthropod initiator motif (Inr). Additionally, potential binding sites for transcription factors involved in the hierarchical control of genes induced by ecdysone during reproduction were identified  (Additional file 3).
AAEL010714 has a paralogue, AAEL010718, with 95% identity in the coding sequences. Although both transcripts are annotated in Vectorbase, RT-PCR experiments using primer pairs F and RT-R or IN-F and RT-R (Additional file 4) resulted in 687 bp and 156 bp amplicons, respectively. The sequences of these amplicons contain a 7 nucleotide sequence that is unique to AAEL010714-RA (Additional file 4), supporting that only this target is amplified. The sequenced cDNAs revealed single nucleotide polymorphisms within the Higgs strain and also between Higgs and Liverpool-derived AAEL010714-RA sequences. Blastx analysis of the conceptual translation product of the 687 bp cDNA sequence using the Vectorbase databank showed the highest similarity with AAEL010714-PA followed by AAEL010718-PA.
Temporal- and tissue-specific expression of AAEL010714
AAEL010714 transcripts are observed only in females and are accumulated highly at 48 hPBM (Figure 1A). The transcripts are most abundant in ovaries with comparatively lower levels in the fat body (Figure 1B). No transcripts are detected in female midguts. Despite the AEL010714 translation product being annotated as an OBP, no corresponding transcripts were detected in head samples (Figure 1C).
A one-way ANOVA statistical analysis performed with normalized values of transcript copy numbers confirmed that AAEL010714 transcripts are most abundant in vitellogenic females at 48 hPBM. The measured transcript levels differ significantly (p < 0.0001) between females at 48 hPBM and all other sample points (Figure 1D).
Detection of AaegOBP45 in mosquito ovaries
Polyclonal antisera against the AaegOBP45-GNQFSSSDIDGL (Figure 2) peptide reacts with a protein present in 48 and 72 hPBM ovaries (Figure 3). The protein is not detected at 24 hPBM, consistent with the gene transcriptional profile (Figure 3). The predicted molecular mass of the mature AaegOBP45 is 26.8 kDa. However, the protein reacting with the antibodies displays a higher apparent molecular mass. Two potential sites for N-linked glycosylation (Asn44 and Asn147) exist in AaegOBP45; therefore, post-translational modifications may account for the difference between the predicted and experimentally-derived apparent mass.
Only two publications to date report the expression profiles of atypical OBPs in Culicidae. These atypical OBPs in An. gambiae are transcribed in chemosensory organs, as well as other tissues in adult and immature aquatic stages . Two of the atypical OBPs in Ae. albopictus are expressed both in preadult stages and in adult hemolymph, and a third is expressed exclusively in antennae and legs . In addition to chemosensory functions, transport of hydrophobic ligands in the hemolymph has been proposed as a role for the atypical Ae. albopictus OBPs.
Although 15 atypical OBPs were annotated in the Ae. aegypti genome , their transcriptional profiles have not been studied in detail. Our work validated the expression profile previously described , and further defined the accumulation of AAEL010714 transcripts in vitellogenic ovaries. Also, we demonstrated that the protein product AaegOBP45 is accumulated in the ovaries leading us to propose a new possible function involved with the reproductive mechanisms responsible for oocyte maturation. It has been shown that seven An. gambiae odorant-binding proteins are associated with mosquito eggshell, and four of them show restricted expression 48 hPBM . Therefore, although functional validation of the atypical/two-domain AaegOBP45 is required, we hypothesize its involvement with eggshell formation.
ALCS and MLC conceived the study. ALCS performed the molecular experiments as RNA extractions, RT-PCR, qRT-PCR, SDS-PAGE and Western blot assays. The author conducted all other analyses of the work as alignments and the statistical analysis and drafted the manuscript. BBK performed DNA genomic amplifications and sequencing, and participated in the design of the study. MLC, OM and AAJ participated in study design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, Drake JM, Brownstein JS, Hoen AG, Sankoh O, Myers MF, George DB, Jaenisch T, Wint GRW, Simmons CP, Scott TW, Farrar JJ, Hay SI: The global distribution and burden of dengue. Nature. 2013, 496: 504-507. 10.1038/nature12060.
Terenius O, Marinotti O, Sieglaff D, James AA: Molecular genetic manipulation of vector mosquitoes. Cell Host Microbe. 2008, 4: 417-423. 10.1016/j.chom.2008.09.002.
Xu P, Atkinson R, Jones DNM, Smith DP: Drosophila OBP LUSH is required for activity of pheromone-sensitive neurons. Neuron. 2005, 45: 193-200. 10.1016/j.neuron.2004.12.031.
Matsuo T, Sugaya S, Yasukawa J, Aigaki T, Fuyama Y: Odorant-binding proteins OBP57d and OBP57e affect taste perception and host-plant preference in Drosophila sechellia. PLoS Biol. 2007, 5: e118-10.1371/journal.pbio.0050118.
Deng Y, Yan H, Gu J, Xu J, Wu K, Tu Z, James AA, Chen X: Molecular and functional characterization of odorant-binding protein genes in an invasive vector mosquito, Aedes albopictus. PLoS One. 2013, 8: e68836-10.1371/journal.pone.0068836.
Hoffman SA, Aravind L, Velmurugan S: Female Anopheles gambiae antennae: increased transcript accumulation of the mosquito-specific odorant-binding-protein OBP2. Parasit Vector. 2012, 5: 27-10.1186/1756-3305-5-27.
Li S, Picimbon JF, Ji S, Kan Y, Chuanling Q, Zhou JJ, Pelosi P: Multiple functions of an odorant-binding protein in the mosquito Aedes aegypti. Biochem Biophys Res Commun. 2008, 372: 464-468. 10.1016/j.bbrc.2008.05.064.
Armbruster P, White S, Dzundza J, Crawford J, Zhao X: Identification of genes encoding atypical odorant-binding proteins in Aedes albopictus (Diptera: Culicidae). J Med Entomol. 2009, 46: 271-280. 10.1603/033.046.0211.
Pelletier J, Leal WS: Genome analysis and expression patterns of odorant-binding proteins from the southern house mosquito Culex pipiens quinquefasciatus. PLoS One. 2009, 4: e6237-10.1371/journal.pone.0006237.
Xu PX, Zwiebel LJ, Smith DP: Identification of a distinct family of genes encoding atypical odorant-binding proteins in the malaria vector mosquito, Anopheles gambiae. Insect Mol Biol. 2003, 12: 549-560. 10.1046/j.1365-2583.2003.00440.x.
Manoharan M, Ng Fuk Chong M, Vaïtinadapoule A, Frumence E, Sowdhamini R, Offmann B: Comparative genomics of odorant binding proteins in Anopheles gambiae, Aedes aegypti and Culex quinquefasciatus. Genome Biol Evol. 2013, 5: 163-180. 10.1093/gbe/evs131.
Vieira FG, Rozas J: Comparative genomics of the odorant-binding and chemosensory protein gene families across the Arthropoda: origin and evolutionary history of the chemosensory system. Genome Biol Evol. 2011, 3: 476-490. 10.1093/gbe/evr033.
Li ZX, Pickett JA, Field LM, Zhou JJ: Identification and expression of odorant-binding proteins of the malaria-carrying mosquitoes Anopheles gambiae and Anopheles arabiensis. Arch Insect Biochem Physiol. 2005, 58: 175-189. 10.1002/arch.20047.
Zhou JJ, He XL, Pickett JA, Field LM: Identification of odorant-binding proteins of the yellow fever mosquito Aedes aegypti: genome annotation and comparative analyses. Insect Mol Biol. 2008, 17: 147-163. 10.1111/j.1365-2583.2007.00789.x.
Hekmat-Scafe DS, Scafe CR, Mckinney AJ, Mark A: Genome-wide analysis of the odorant-binding protein gene family in Drosophila melanogaster. Genome Res. 2002, 12: 1357-1369. 10.1101/gr.239402.
Zhou JJ, Huang W, Zhang GA, Pickett JA, Field LM: “Plus-C” odorant-binding protein genes in two Drosophila species and the malaria mosquito Anopheles gambiae. Gene. 2004, 327: 117-129. 10.1016/j.gene.2003.11.007.
Dissanayake SN, Ribeiro JM, Wang MH, Dunn WA, Yan G, James AA, Marinotti O: aeGEPUCI: a database of gene expression in the dengue vector mosquito, Aedes aegypti. BMC Res Notes. 2010, 3: 248-10.1186/1756-0500-3-248.
Amenya DA, Chou W, Li J, Yan G, Gershon PD, James AA, Marinotti O: Proteomics reveals novel components of the Anopheles gambiae eggshell. J Insect Physiol. 2010, 56: 1414-1419. 10.1016/j.jinsphys.2010.04.013.
Costa-da-Silva AL, Navarrete FR, Salvador FS, Karina-Costa M, Ioshino RS, Azevedo DS, Rocha DR, Romano CM, Capurro ML: Glytube: a conical tube and parafilm M-based method as a simplified device to artificially blood-feed the dengue vector mosquito, Aedes aegypti. PLoS One. 2013, 8: e53816-10.1371/journal.pone.0053816.
Chen L, Zhu J, Sun G, Raikhel AS: The early gene broad is involved in the ecdysteroid hierarchy governing vitellogenesis of the mosquito Aedes aegypti. J Mol Endocrinol. 2004, 33: 743-761. 10.1677/jme.1.01531.
Zhu J, Chen L, Raikhel AS: Distinct roles of broad isoforms in regulation of the 20-hydroxyecdysone effector gene, Vitellogenin, in the mosquito Aedes aegypti. Mol Cell Endocrinol. 2007, 267: 97-105. 10.1016/j.mce.2007.01.006.
Sandelin A, Alkema W, Engström P, Wasserman WW, Lenhard B: JASPAR: an open-access database for eukaryotic transcription factor binding profiles. Nucleic Acids Res. 2004, 32 (Database issue): D91-D94.
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999, 41: 95-98.
Bonizzoni M, Britton M, Marinotti O, Dunn WA, Fass J, James AA: Probing functional polymorphisms in the dengue vector, Aedes aegypti. BMC Genomics. 2013, 14: 739-10.1186/1471-2164-14-739.
Sanders HR, Evans AM, Ross LS, Gill SS: Blood meal induces global changes in midgut gene expression in the disease vector, Aedes aegypti. Insect Biochem Mol Biol. 2003, 33: 1105-1122. 10.1016/S0965-1748(03)00124-3.
We dedicate this paper to the memory of Alexandre A. Peixoto.
We thank Ediane Saraiva Fernandes and Neuza Saraiva Fernandes for the technical assistance. This work was funded by grant #2006/00331-3 from São Paulo Research Foundation (FAPESP). ALCS was supported by CNPq fellowship #143113/2006-2. OM and AAJ were supported in part by a grant from the NIH NIAID (AI29746).
The authors declare that they have no competing interests.
Electronic supplementary material
Additional file 2: Quantification of β-actin accumulated transcripts by qRT-PCR in adult females 36, 48, 60, 72, 84, 96 hPBM. RNA was extracted from three independent experiments, except for 84 hPBM (two independent), and each time point were measured three times. mRNA copy numbers were determined by comparison with known concentrations of a standard plasmid. No significant differences (p > 0.05) were observed between all the samples analyzed. All data represent mean values and error bars indicate standard error of the mean. (PPTX 130 KB)
Additional file 3: Nucleotide sequence of the AAEL010714 gene from Ae. aegypti generated by manual sequencing. Pentamer nucleotides under-lined with dots are a consensus putative arthropod initiator cap-site and the bold under-lined nucleotide is the putative start of transcription (+1). Solid and dashed underlined sequences represent exons (1 and 2) and introns, respectively. The putative TATA sequence is boxed. One putative hunchback-binding site is in a dashed box and potential BR-C gene family zinc-finger binding-sites are double underlined (non-aligned lines) highlighting two overlapped motifs. Nucleotide numbers are presented to the right and position was defined based on the putative first transcribed nucleotide. (PPTX 182 KB)
Additional file 4: Alignment of the transcript nucleotide sequences AAEL010714-RA and AAEL010718-RA from Vectorbase with nucleotide sequences obtained experimentally. BioEdit was used to align the sequences. Regions where primers were designed are boxed in gray line. Primer F and Primer RT-R were used for cDNA amplification. Primer IN-F and RT-R were used for cDNA amplification in tissue RT-PCR and Real-time RT-PCR experiments. Black boxed gap indicates that 7 nucleotides are missing in AAEL010718-RA (718-RA) which are present in AAEL010714-RA (714-RA) and in sequences obtained by RT-PCR (cDNA_partial/1-687 and cDNA partial/1-156). Dots represent gaps. Polymorphisms between all sequences are shaded gray. (PPTX 114 KB)
About this article
Cite this article
Costa-da-Silva, A.L., Kojin, B.B., Marinotti, O. et al. Expression and accumulation of the two-domain odorant-binding protein AaegOBP45 in the ovaries of blood-fed Aedes aegypti. Parasites Vectors 6, 364 (2013) doi:10.1186/1756-3305-6-364
- Aedes aegypti
- Odorant-binding protein