Expression and accumulation of the two-domain odorant-binding protein AaegOBP45 in the ovaries of blood-fed Aedes aegypti
© Costa-da-Silva et al.; licensee BioMed Central Ltd. 2013
Received: 15 October 2013
Accepted: 10 December 2013
Published: 24 December 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
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
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.
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).
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.