- Research article
- Open Access
Integrated proteomic and transcriptomic analysis of the Aedes aegyptieggshell
BMC Developmental Biologyvolume 14, Article number: 15 (2014)
Mosquito eggshells show remarkable diversity in physical properties and structure consistent with adaptations to the wide variety of environments exploited by these insects. We applied proteomic, transcriptomic, and hybridization in situ techniques to identify gene products and pathways that participate in the assembly of the Aedes aegypti eggshell. Aedes aegypti population density is low during cold and dry seasons and increases immediately after rainfall. The survival of embryos through unfavorable periods is a key factor in the persistence of their populations. The work described here supports integrated vector control approaches that target eggshell formation and result in Ae. aegypti drought-intolerant phenotypes for public health initiatives directed to reduce mosquito-borne diseases.
A total of 130 proteins were identified from the combined mass spectrometric analyses of eggshell preparations.
Classification of proteins according to their known and putative functions revealed the complexity of the eggshell structure. Three novel Ae. aegypti vitelline membrane proteins were discovered. Odorant-binding and cysteine-rich proteins that may be structural components of the eggshell were identified. Enzymes with peroxidase, laccase and phenoloxidase activities also were identified, and their likely involvements in cross-linking reactions that stabilize the eggshell structure are discussed.
Insect eggshells provide embryos protection from environmental insults, both physical and biological [1–3]. Although the primary protective function of eggshells is conserved among insects, differences in ecology, in particular oviposition substrates, are strong forces driving the evolution of eggshell morphology, organization and composition [4–6]. Accordingly, mosquito (Diptera, Culicidae) eggshells show a remarkable diversity in physical properties and structure permitting adaptations to a wide variety of environments. For example, Aedes aegypti eggs are highly-resistant to desiccation allowing embryos to survive for months in dry conditions [7, 8], while the eggshells of other mosquito species, including anophelines, which comprise human malaria vectors, are more permeable, restricting their development and survival to humid environments [9, 10]. Greater knowledge of the proteins that comprise mosquito eggshells is required to understand these differences and to elucidate how they contribute to mosquito reproductive strategies.
Only a small number of Ae. aegypti eggshell components have been characterized [11–20]. Here we applied proteomic, transcriptomic, and hybridization in situ techniques to identify additional components and describe processes that participate in eggshell assembly. Understanding roles for the discovered proteins in long-term embryo survival may lead to novel approaches for intervening in the reproduction of these mosquitoes and consequently reduce vector-borne transmission of pathogens to humans.
Results and discussion
Aedes aegypti eggshells were isolated and their proteins resolved by SDS-PAGE (Figure 1). Several proteins are visibly enriched in the purified eggshell samples when compared with whole ovaries extracts prepared at 72 hours post blood meal (hPBM). Ten gel slices (A-J) containing the most abundant eggshell proteins identified by Coomassie blue staining were analyzed by mass spectrometry. Thirty-four distinct proteins were identified in this analysis (Additional file 1).
The buffer containing SDS/β-mercaptoethanol used for the SDS-PAGE sample preparations did not dissolve all eggshell components since a pellet was recovered by centrifugation of the samples prior to electrophoresis. Protein crosslinking reactions may form higher-order covalent complexes and prevent solubilization of some eggshell proteins, or prevent them from entering the polyacrylamide gel. A second, ‘shotgun’, approach not based on SDS-PAGE was used to mitigate the possible under-representation of proteins in these samples. Tryptic digestion of the whole eggshell preparations followed by mass spectrometry discovered 127 distinct proteins (Additional file 2). All but three of the proteins identified from the excised gel slices were present in the new list of proteins identified from the second experimental approach. Therefore, a total of 130 proteins were identified from the combined mass spectrometric analyses of eggshell preparations. Differences in the two proteomes may arise from various factors including (a) different eggshell preparations, (b) large (≥ 500 kDa) and small (≤ 5 kDa) proteins not resolved during electrophoresis, (c) proteins between gel slices that were not excised after SDS-PAGE, and (d) differences in thresholds parameters for database search results.
Temporal patterns of transcript accumulation – transcriptome and quantitative RT-PCR
An insoluble fraction was recovered from the eggshell preparation even after treatment with trypsin. Therefore, we undertook a complementary, non-protein based approach to evaluate the complexity of the eggshells. The transcriptome of ovaries at the latter stages of eggshell formation (60-72 hPBM) was compared with the published transcriptome of recently-laid Ae. aegypti eggs . Transcripts enriched in the late ovaries and either absent or present at negligible levels in embryos most likely correspond to genes that are expressed in the follicle cells. These are candidates for encoding components of the outermost eggshell layers. Transcripts enriched in the ovaries and maintained at high levels in embryos could be maternally-derived transcripts that serve functions in early embryonic development [22–24].
A total of 536 transcripts were found to be at least 10-fold more represented in ovaries than in embryos (p values < 5E−05) (Additional file 3). Publicly-available expression profiles of those transcripts were retrieved and analyzed [25, 26] (Additional file 4). We performed quantitative RT-PCR to examine accumulation profiles of selected transcripts to validate the high-throughput derived data (Additional file 5). Quantitative RT-PCR, microarray  and RNA-Seq  derived data were concordant for the majority (>90%), and evidenced genes with distinct expression patterns.
Not all proteins identified in the eggshell proteome described above were encoded by these transcripts and this is consistent with differential expression profiles of the eggshell components [25, 26] (Additional files 4 and 5). For example, vitelline membrane proteins (VPMs) encoding genes are expressed during early oogenesis and were not represented among the abundant transcripts during the late stages of egg formation. Also, a number of the transcripts present in ovaries and absent, or detected at low levels in embryos, encode predicted intracellular components that although displaying differential expression profiles are not likely components of the eggshell.
Spatial patterns of transcript accumulation - hybridizations in situ
Hybridizations in situ were performed to further validate expression of four selected genes in vitellogenic ovaries (Figure 2). Three odorant binding protein (OBP) encoding transcripts (AAEL01714-RA, AAEL006396-RA, AAEL006387-RA) accumulated in the follicle cells of ovaries dissected at 48 hPBM. The transcripts of a gene (AAEL007112) encoding a protein with a chitin-binding domain, also accumulated in follicle cells. Hybridizations in situ with ovaries dissected later than 48 hPBM to map accumulation of transcripts encoding cysteine-rich proteins (CRPs) failed. Previous attempts to perform hybridizations in situ at these late stages of oogenesis also failed . Structural and/or physiological changes in the ovaries during late egg formation could either degrade the probes or block their access to targets within the follicle cells.
Functional assignments of eggshell components
Grouping the transcripts and proteins identified in this study according to known and putative functions and expression profiles reveals multiple components participating in a sequential process of eggshell formation (Figure 3). Morphological aspects of the events resulting in mosquito eggshell formation have been described previously [27–29]. The proteins identified here provide molecular support to the earlier ultrastructural observations. Genes either expressed constitutively or specifically at early, intermediate or late stages of egg formation were grouped in six categories described below, five of which are listed in Figure 3.
Vitelline membrane proteins
The vitelline membrane is the inner layer of the eggshell and vitelline membrane proteins (VMPs) are its main components. Three Ae. aegypti VMPs, 15a-1, 15a-2 and 15a-3  were recovered from gel slices (A and B) that contain proteins with MW values between 10 and 15 kDa (Figure 1), sizes consistent with their predicted sequences. Gel slices C and D, containing proteins with MWs between 15 and 18 kDa, also contained VMPs 15a-1 and 15a-3, supporting the interpretation that alternative post-translational modification or other factors may result in their abnormal and distinctive SDS-PAGE migration.
The products of AAEL006670, identified in gel slice B, and AAEL017471 and AAEL0017501, identified in the whole eggshell proteome, are VMP paralogs and annotated here as novel Ae. aegypti VMPs. Two characteristics that define insect vitelline membrane proteins are essential hydrophobic domains that are highly-conserved in all members of the protein family and a short period of expression and accumulation at the onset of egg formation . The three novel Ae. aegypti VPMs display hydrophobic domains of 30-50 amino acids in length and their corresponding transcripts are up-regulated following a blood meal, patterns consistent with their proposed function (Figures 3 and 4; Additional files 4, 5 and 6). AAEL017471-RA is unique among the Ae. aegypti VMP-encoding transcripts in that its accumulation was maximum at 48 hPBM, while all other VMP transcripts attained their highest levels at 24-36 hPBM. A similar pattern was observed with one of the D. melanogaster VMPs, VM32E, which is expressed only at stage 10 of oogenesis and considered a ‘late’ vitelline membrane gene . VM32E is an integral component of both the vitelline membrane and the endochorion and it may serve as a bridge between the two eggshell layers . Therefore, VM32E and AAEL017471 could perform comparable functions.
Proteins with chitin-binding and chitinase-like domains
Chitin is present in several insect structures at various developmental stages. In addition to the exoskeleton, chitin and chitin-like materials are found in the peritrophic matrix (PM), embryos, ovaries, compound eyes and mosquito eggshells [32–34]. Proteins with chitin-binding domains and/or chitinase domains (CBPs), products of AAEL006328, AAEL002400, AAEL000144, and AAEL007112, are included in our datasets of eggshell components (Figure 3, Additional files 1 and 2). AAEL002400-PA and AAEL000144-PA are present in gel slices H (Molecular Weight, MW ~ 60 kDa) and J (250 kDa) (Figure 1), possibly a result of their association with higher-order covalent complexes.
Feeding the chitin synthesis inhibitor lufenuron to adult female mosquitoes results in inhibited oogenesis and egg hatching. Furthermore, eggs laid by lufenuron-treated mosquitoes are fragile  and this supports an important role for chitin in maintaining egg integrity. Our findings of chitin-binding proteins in Ae. aegypti eggshells support the hypothesis that an organized structure composed of chitin and chitin-binding proteins provides rigidity to the mosquito eggshell. Eggshell proteins with chitinase domains may play a role in embryo hatching [33, 35] and be stored as proenzymes that are activated immediately prior to hatching.
Odorant binding proteins
Gel slices E, F and G, containing proteins migrating with apparent MWs between 30 and 48 kDa (Figure 1) and the shotgun proteomic approach revealed a series of odorant binding proteins (OBPs) as eggshell components (Additional files 1 and 2). Insect OBPs are reported to be expressed in insect antennae and are proposed to bind small hydrophobic odorant molecules and carry and present them to the olfactory receptors . Structurally-related proteins also have been identified in non-sensory organs of insects, supporting the conclusion that OBPs and OBP-like proteins could have roles unrelated to olfaction. For example, D7-related (D7r) proteins similar in structure to OBPs are abundant components in the saliva of hematophagous insects. The An. gambiae D7r proteins have been shown to bind the biogenic amines serotonin, norepinephrine, and histamine with high affinity thereby reducing the concentrations of these effectors at the feeding site and resulting in an anti-hemostatic activity . The Tenebrio molitor THP12 gene encodes a member of the OBP family that is accumulated in the hemolymph and has the capacity to bind fatty acids with high affinity, a finding consistent with a role in the transport of small hydrophobic molecules .
All OBPs identified in the Ae. aegypti eggshell preparations, with the exception of AAEL011483-PA (AaegOBP 50c), belong to the atypical, two-domain subfamily or two-domain clade of OBPs [39, 40]. As the nomenclature indicates, proteins in this subfamily have two OBP domains that likely originated from the fusion of two Classic OBP genes. AAEL011483 encodes a protein that belongs to the “PlusC” OBP family, however its expression pattern indicates that it is not expressed in the ovaries (Additional file 4). AAEL000319-RA (AaegOBP93), although expressed in ovaries, showed enhanced expression in larvae and likely has an alternative function. Consistent with our findings, six of eight OBPs identified in the eggshell of An. gambiae also belong to the atypical, two-domain subfamily or two-domain clade. Although little is known about natural ligands of the OBPs expressed and accumulated in the ovaries, they could carry substrates for eggshell protein crosslinking and assembly reactions. Alternatively or additionally, odorant molecules may act as chemo-attractants for sperm , and the OBPs identified in this study could enrich specific attractant molecules at the surface of the untanned mosquito eggshell.
The accumulation of OBPs in vitellogenic mosquito ovaries evidenced by immunoblot analyses , their high abundance in the eggshell demonstrated by SDS-PAGE and proteomics in this study, and their transcript accumulation patterns support the hypothesis that they may be structural components of the chorion intermediate layer. OBPs delivering crosslinking substrates could become trapped within the eggshell structure and may have been co-opted through evolutionary selection processes for additional functions.
Cysteine-rich proteins - putative chorion proteins
A group of small proteins displaying high cysteine content were described as putative An. gambiae chorion proteins . Similarly, we found 10 cysteine-rich proteins (CRPs) as components of the Ae. aegypti eggshell. They compose a group of short, putatively secreted proteins that contain a trypsin inhibitor-like, cysteine-rich domain (IPR002919). The genes encoding these proteins are among those expressed most highly in ovaries at 60-72 hPBM and their transcripts are absent from embryos. This pattern supports their function as components of the most external layers of the chorion deposited during the last stages of oogenesis. An alignment of the eggshell CRPs sequences revealed a conserved motif 59 amino acids in length that defines a hydrophobic domain possibly analogous to that seen in the vitelline membrane proteins (Figure 5, Additional file 7).
Chorion proteins are proposed to have evolved rapidly, providing adaptation to stresses such as desiccation, plant toxins, pathogen invasion and predation . Culicidae (mosquito) and Drosophilidae (fruit fly) lineages are estimated to have diverged ∼ 250 million years ago  and their females lay eggs in ecologically-distinct sites. The strong selective pressures on chorion components resulting from their respective life histories provide opportunities for engagement of distinct genes and gene families in the process of eggshell assembly. We could not identify any mosquito eggshell component with amino acid sequence similarity to the well-characterized D. melanogaster chorion proteins . The mosquito CRPs possibly function as structural components of the external chorionic layer of the eggshell and their gene families have many members in both Ae. aegypti and An. gambiae. Interestingly, the genes encoding eggshell CRPs in Ae. aegypti are not annotated as orthologous to those described previously as encoding An. gambiae putative chorion proteins .
Peroxidases, phenoloxidases, laccases, dopachrome-conversion enzymes, and transglutaminases are enzymes present in the eggshell samples (Figure 3, Additional files 1, 2 and 3). Peroxidase-catalyzed chorion protein crosslinking is a phenomenon documented previously in insect eggshell formation  with AAEL004386-PA, HPX8C, identified in the proteome, having been characterized previously as an Ae. aegypti chorion peroxidase (CPO) . The peroxidases AAEL004390-PA (HPX8B) and AAEL000507-PA also were present in our datasets (gel slices H and I, respectively) likely representing new enzymes involved in eggshell formation. Laccase (AAEL016992) was identified as an Ae. aegypti eggshell component. Laccases have p-diphenol oxidase enzymatic activity and participate in cuticular tanning in insects. Pan et al.  reported that Culex pipiens pallens laccase 2 (CpLac2) is expressed abundantly in the egg developmental stages and postulated a role in eggshell tanning. The presence of the An. gambiae Laccase 2 transcripts (AGAP006176-RA and -RB) in ovaries and corresponding peptides in the eggshell provide additional support for this hypothesis . Pigmentation of Ae. albopictus eggshells was blocked when laccase 2 expression was knocked down by RNAi . Furthermore, eggs produced after laccase 2 knockdown showed incomplete sclerotization, were deformed and did not complete development .
Phenoloxidases are involved in chorion melanization . AePPO5 (AAEL013492-PA) is present in Ae. aegypti eggshells (Additional file 2), however, the accumulation profile of its corresponding transcript indicates that the enzyme is synthesized mainly in tissues other than ovaries  (Additional file 4). AAEL013492 is an ortholog of the An. gambiae prophenoloxidase 9 (AgPPO9; AGAP004978-PA) and the Cu. quinquefasciatus product of CPIJ016564 . AePPO5 and AgPPO9 lack secretory signal peptides, and mechanisms involved in the transport of these enzymes to the eggshell are unknown. A putative transglutaminase (AAEL007641-PA) identified as a component of the eggshell (Additional file 2) also has an expression pattern inconsistent with its synthesis in the ovaries  (Additional file 4). Transglutaminases participate in multiple biological processes by catalyzing protein crosslinking and playing a key role in cuticle morphogenesis and sclerotization in insects . Therefore, AAEL007641-PA may catalyze protein cross-linking during Ae. aegypti eggshell formation.
The translation products of AAEL006985, AAEL006830, AAEL007096, AAEL010848, and AAEL002333 were identified in gel slices F and H and contain motifs similar to those found in the Major Royal Jelly protein family. This protein family includes the products of the Drosophila yellow genes, which encode dopachrome-conversion enzymes that catalyze cuticle melanization. The yellow-g and yellow-g2 gene products participate in the cross-linking of fruit fly eggshell proteins . Those and our proteome and gene expression data support the hypothesis that mosquito Major Royal Jelly proteins play a catalytic role in the crosslinking of eggshell proteins. However, dopachrome conversion enzyme-encoding genes also are expressed in tissues other than ovaries and at various developmental stages, supporting a role for them in other mosquito metabolic functions. Our proteome data did not include AAEL014238-PA, a dopachrome-conversion enzyme described previously as participating in Ae. aegypti eggshell formation . Its absence indicates that the eggshell protein set reported here, although larger than those included in the previous reports, may be incomplete.
Other cellular and structural proteins
Proteins with similarity to cytoplasmic or cell membrane components were identified in eggshell preparations and could indicate novel functions for these molecules (Additional files 1 and 2). Actin (AAEL001673-PA), gamma glutamyl transpeptidase (AAEL004240-PA), transcriptional regulator ATRX (AAEL010502-PA), ATP-binding cassette transporter (AAEL012192-PA), vacuolar sorting protein (AAEL005861-PA), glycosyl transferase (AAEL001590-PA), coronin (AAEL009249-PA), alpha (1,3) fucosyltransferase (AAEL000244-PA), and a protein of unknown function (AAEL002382-PA) were identified by the proteomics approach. However, a number of non-eggshell proteins may co-purify with the eggshells with the methods used in this study. Other molecules with putative functions in eggshell formation identified by our differential transcriptomic methodology include serine-type endopeptidases (AAEL014188-PA, AAEL003060-PA), a putative 67B2 family heat shock protein (AAEL007339-PA) and myosin (AAEL005733-PA). Although actin, myosin, heat shock proteins and serine proteinases also were identified in D. melanogaster eggshell preparations , their functions in the mosquito eggshell remain to be validated.
Aedes aegypti is an important world-wide vector of dengue and yellow fever viruses. Therefore, the management of Ae. aegypti populations is of prime importance to public health and welfare. In general, Ae. aegypti population densities are low during cold and dry seasons and increase immediately after rainfall. This pattern is due to the hatching of a large number of estivated, drought-resistant larvae that hatch in response to water accumulation in natural and artificial breeding sites. Mosquito survival until the next favorable climatic period is then a key factor for the persistence of Ae. aegypti populations . Thus, targeting Ae. aegypti eggs to make them less resistant, or intolerant to desiccation would be beneficial.
Our data fill a gap in knowledge regarding the proteins and sequential processes involved in Ae. aegypti eggshell formation. Proteins involved in the hardening and tanning processes of the mosquito eggshell may represent targets for the development of enzyme-specific inhibitors that added to the egg-laying sites and/or fed to females via blood or nectar could result in abnormal egg formation or a defective chorion. For example, phoroglucinol, a natural phenolic peroxidase inhibitor, affects eggshell formation when fed to adult female insects . This substance or derivatives are potential candidates for controlling vector populations in integrated insect control programs. High-throughput searches for molecules with inhibitory activity on enzymes [53, 54] also can be applied to insect control. The discovery or development of highly specific, Ae. aegypti eggshell enzyme inhibitors would be advantageous and desirable to minimize unwanted environmental effects of insect control strategies on non-target organisms. Design of insecticides that target cuticle and chorion formation, such as lufenuron used to control fleas , are possible with an understanding of the biochemical and physiological processes involved in eggshell formation in mosquitoes.
The data documented here also support control strategies involving transgenic insects [56–59]. Transgenic mosquitoes with a drought-intolerant embryo phenotype could be engineered, easily raised under controlled humidity laboratory conditions, and upon release would display little or no fitness cost during rainy seasons. However, their eggs would not survive through dry seasons. The population suppression proposed here together with other efforts that result in local or regional vector elimination are anticipated to delay the onset of increased mosquito populations that follows the onset of a rainy season, result in reduced pathogen transmission and provide a consequent reduction in human morbidity and mortality.
The Ae. aegypti Liverpool strain was maintained at 26°C ± 1°C with 77% humidity and 12 h day/night, 30 min dusk/dawn lighting cycle. Larvae were fed a diet of powdered fish food (Tetramin). Adult males and females were kept together in a cage with unlimited access to water and sugar (10% sucrose) until blood feeding. Mosquitoes aged 3-5 days after eclosion were allowed to feed on anaesthetized mice. The study was carried out in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health.
Eggshell isolation was carried out as described  with modifications. Ovaries were dissected from females at 72 hPBM in 0.15 M NaCl and homogenized in a HEPES buffer, pH 8.0 (10 mM HEPES, 130 mM NaCl, 4.7 mM KCl, 0.5 mM phenyl thiocarbamide, 0.1 mM p-nitrophenyl-p’-guanidinobenzoate, 5 mM ethylene diamine tetraacetic acid, 1% Triton X-100, 0.1 mM diethyldithiocarbamic acid, 1 mM phenylmethylsulfonyl fluoride). Eggshells were collected by sedimentation for 15 min on ice and the supernatant decanted. The sediment was washed five times for 15 min each with homogenization buffer. Eggshells then were sonicated for 30 seconds and centrifuged (100 g, 10 min, 4°C). Five additional cycles of resuspending the pellet in homogenization buffer followed by centrifugation were performed to remove cellular contaminants.
Purified eggshell and whole ovary samples were dissolved in SDS/β-mercaptoethanol for 2 min at 100°C and then resolved by polyacrylamide gel electrophoresis. Proteins were stained with Coomassie Brilliant Blue R-250 and excised as slices from the gel using a sharp scalpel. Gel slices were destained in 100 μl of 100 mM ammonium bicarbonate and 50% acetonitrile for 10 min at room temperature (RT) while agitating at 700 rpm. Two identical destaining steps were performed. The gel slices were dehydrated in 30 μl of acetonitrile for 10 min at RT, while agitating at 700 rpm. The acetonitrile was removed and the dehydration step was repeated. After removal of the acetonitrile, gel slices were dried in a Savant ISS110 SpeedVac Concentrator for 10 min at 45°C. Disulfide bonds were disrupted in 50 μL of 10 mM dithiothreitol in 50 mM ammonium bicarbonate, incubating for 1 hr at 60°C. The dithiothreitol solution was removed and cysteines were alkylated with 50 mM iodoacetamide in 50 mM ammonium bicarbonate during incubation for 45 min at 45°C in the dark. Gel slices then were washed with 20 μl of 50 mM ammonium bicarbonate, followed by a wash with 50 μl of 100 mM ammonium bicarbonate. After two 10 min washes with 30 μl of acetronitrile, gel slices were dried for 5-10 min under vacuum. Each sample was incubated in trypsin (50 μL at 20 ng/mL) for 45 min at 4°C and then at 37°C, overnight. Tryptic peptides were extracted by adding 30 μL of 50% (v/v) acetonitrile/1% (v/v) formic acid and incubating for 10 min at RT. The extraction solution was transferred to another tube and the extraction step was performed three additional times. All extracts of each gel slice were combined and samples evaporated to dryness in a rotary vacuum and reconstituted in 20 μL of HPLC loading buffer consisting of 50% (v/v) acetonitrile/1.0% (v/v) formic acid. LC-MS was carried out for data acquisition using a Waters ACQUITY UPLC® system coupled to a SYNAPT®G2 mass spectrometer. Injected samples were washed on a trap column prior to being loaded on a separation column. The sample was separated by a linear gradient from 95% A (0.1% v/v formic acid in water) and 5% B (0.1% v/v formic acid in acetonitrile) to 65% A and 35% B over 90 min. Masslynx 3.5 software (Waters Corporation, Milford, MA, USA) was used for data acquisition, and processing. Peptides were identified using ProteinLynx Global Server Web (Waters) with SwissProt database for Ae. aegypti.
Whole eggshell nanoLC-MS/MS
Protein extraction and digestion for mass spectrometry was adapted from . Briefly, eggshells were dissolved in a buffered solution containing sodium dodecyl sulfate (SDS) and tris(2-carboxyethyl)phosphine (TCEP) (Pierce) and the resulting solution placed on a 30 kDa ultrafiltration unit (YM-30, Millipore) which was used for buffer exchange, carbamidomethylation of cysteine, and protein trypsinization. Eluted peptides were cleaned in pipette tips packed with strong cation-exchange (SCX) and C18 filters (3 M Inc.) . Peptides eluted from the C18 phase were transluted to the SCX phase to remove remaining contaminants. Peptides were eluted with ammonium acetate and subjected to nanoLC-MS/MS on an LTQ Velos Pro (ThermoFisher) with EASY-nLC 1000 running a 3 hr acetonitrile gradient at a flow rate of 0.2 nL/min. The Velos Pro instrument method fragmented the 20 most-abundant ions per precursor spectrum via collision induced dissociation, with dynamic exclusion. Searches using Mascot 2.4 were against an Ae. aegypti custom database (alongside a database of common contaminants), with trypsin specificity. The sole fixed and variable modifications were Carbamidomethyl (C) and Oxidation (M), respectively. Peptide and fragment mass tolerances were ± 2 Da and ± 0.8 Da, respectively, with a maximum of one missed cleavage.
Ilumina transcriptomes were generated as in . Data from two samples, 0-2 hour embryos  and ovaries dissected 60-72 hPBM (this work), were compared. The 0-2 hr embryo sample has 6,629,888 reads 33 bp in length, and the 60-72 hPBM ovary sample comprises 15,691,856 40 bp long reads. Reads were aligned by bowtie  to the 18,769 annotated Ae. aegypti transcripts (AaegL1.3 database, http://vectorbase.org/) to calculate relative transcript accumulation. The number of reads mapped to each transcript was counted and normalized, resulting in values expressed as reads per kilobase per million mapable reads (RPKM).
The DEGseq method in the DEGseq package was used to determine the significance and fold-changes of transcript abundance between samples. The MA-plot-based method with random sampling model (MARS) with the default p-value < 0.001 was used for our analysis. The input data were the values of normalized mapped reads for each transcriptome sample.
RNA samples were extracted using TRIZOL (Invitrogen) from pools of ovaries dissected from females kept exclusively on a sugar diet or collected at defined intervals after blood feeding. Following DNAse I (Invitrogen) treatment, a total of 10 μg of RNA were used for cDNA synthesis with superscript III (Invitrogen) and random primers. Real-time quantitative PCR reactions of 20 μl were performed in triplicate with SYBR Green Supermix (Biorad) and 0.3 μM of each primer (Additional file 8) on three sequential five-fold dilutions each of the original cDNA. Real-time quantitative PCR reactions were run on an iQ3 system (Biorad). No primer dimers were detected when inspecting the melting curves and primer pairs with at least 90% amplification efficiency were selected. Fold-changes in mRNA abundance were derived by the comparative Ct method  using a constitutively-expressed gene, rp49 (GenBank Acc. No.:AY539746; AAEL003396) as the reference transcript.
Hybridization in situ
Hybridizations in situ in whole-mount ovaries were conducted as described previously . A One-step RT-PCR kit (Qiagen, Valencia, CA) was used to generate the probes, for subsequent cDNA amplification reactions. Sequences of the oligonucleotide primers are listed in Additional file 8. Reaction mixtures were incubated for 30 min at 50°C and 15 min at 95°C. Amplification conditions were 3 min at 95°C followed by 30 cycles of 30 s at 95°C, 30 s at 60°C and 1 min at 72°C. RT-PCR products were cloned into the pCR®4-TOPO® cloning vector (Invitrogen) and sequenced to confirm their identity. Digoxygenin (DIG)-labeled sense and antisense RNA probes for each gene product was synthesized in vitro using T3 or T7 RNA polymerases (Ambion, Austin, TX).
Margaritis LH: Structure And Physiology Of The Eggshell. Comprehensive Insect Physiology, Biochemistry and Pharmacology. Volume I Embryogenesis and Reproduction. Edited by: Kerkut GA, Gilbert LI. 1985, Oxford: Pergamon Press, 153-230.
Waring GL: Morphogenesis of the eggshell in Drosophila. Int Rev Cytol. 2000, 198: 67-108.
Cavaliere V, Bernardi F, Romani P, Duchi S, Gargiulo G: Building up the Drosophila eggshell: first of all the eggshell genes must be transcribed. Dev Dyn. 2008, 237: 2061-2072. 10.1002/dvdy.21625.
Jagadeeshan S, Singh RS: Rapid evolution of outer egg membrane proteins in the Drosophila melanogaster subgroup: a case of ecologically driven evolution of female reproductive traits. Mol Biol Evol. 2007, 24: 929-938. 10.1093/molbev/msm009.
Suman DS, Shrivastava AR, Parashar BD, Pant SC, Agrawal OP, Prakash S: Variation in morphology and morphometrics of eggs of Culex quinquefasciatus mosquitoes from different ecological regions of India. J Vector Ecol. 2009, 34: 191-199.
Sahlén G: Eggshell ultrastructure in four mosquito genera (Diptera, culicidae). J Am Mosq Control Assoc. 1996, 12: 263-270.
Christophers SR: Aedes Aegypti. The Yellow Fever Mosquito. Its Life History, Bionomics And Structure. 1960, London: Cambridge University Press
Sota T, Mogi M: Survival time and resistance to desiccation of diapause and non-diapause eggs of temperate Aedes (Stegomyia) mosquitoes. Entomol Exp Appl. 1992, 63: 155-161. 10.1111/j.1570-7458.1992.tb01570.x.
Deane MP, Causey OR: Viability of Anopheles gambiae eggs and morphology of unusual types found in Brazil. Am J Trop Med Hyg. 1943, 23: 95-102.
Beier JC, Copeland R, Oyaro C, Masinya A, Odago WO, Oduor S, Koech DK, Roberts CR: Anopheles gambiae complex egg-stage survival in dry soil from larval development sites in western Kenya. J Am Mosq Control Assoc. 1990, 6: 105-109.
Lin Y, Hamblin MT, Edwards MJ, Barillas-Mury C, Kanost MR, Knipple DC, Wolfner MF, Hagedorn HH: Structure, expression, and hormonal control of genes from the mosquito, Aedes aegypti, which encode proteins similar to the vitelline membrane proteins of Drosophila melanogaster. Dev Biol. 1993, 155: 558-568. 10.1006/dbio.1993.1052.
Li J: Egg chorion tanning in Aedes aegypti mosquito. Comp Biochem Physiol. 1994, 109A: 835-843.
Ferdig MT, Li J, Severson DW, Christensen BM: Mosquito dopa decarboxylase cDNA characterization and blood-meal-induced ovarian expression. Insect Mol Biol. 1996, 5: 119-126. 10.1111/j.1365-2583.1996.tb00046.x.
Li J, Hodgeman BA, Christensen BM: Involvement of peroxidase in chorion hardening in Aedes aegypti. Insect Biochem Mol Biol. 1996, 26: 309-317. 10.1016/0965-1748(95)00099-2.
Edwards MJ, Severson DW, Hagedorn HH: Vitelline envelope genes of the yellow fever mosquito, Aedes aegypti. Insect Biochem Mol Biol. 1998, 28: 915-925. 10.1016/S0965-1748(98)00083-6.
Han Q, Li G, Li J: Purification and characterization of chorion peroxidase from Aedes aegypti eggs. Arch Biochem Biophys. 2000, 378: 107-115. 10.1006/abbi.2000.1821.
Yao R, Li J: Towards global analysis of mosquito chorion proteins through sequential extraction, two-dimensional electrophoresis and mass spectrometry. Proteomics. 2003, 3: 2036-2043. 10.1002/pmic.200300536.
Li J, Kim SR, Li J: Molecular characterization of a novel peroxidase involved in Aedes aegypti chorion protein crosslinking. Insect Biochem Mol Biol. 2004, 34: 1195-1203. 10.1016/j.ibmb.2004.08.001.
Kim SR, Yao Y, Han Q, Christensen BM, Li J: Identification and molecular characterization of a prophenoloxidase involved in Aedes aegypti chorion melanization. Insect Mol Biol. 2005, 14: 185-194. 10.1111/j.1365-2583.2004.00547.x.
Li JS, Li J: Major chorion proteins and their crosslinking during chorion hardening in Aedes aegypti mosquitoes. Insect Biochem Mol Biol. 2006, 36: 1195-1203.
Biedler JK, Hu W, Tae H, Tu Z: Identification of early zygotic genes in the yellow fever mosquito Aedes aegypti and discovery of a motif involved in early zygotic genome activation. PLoS One. 2012, 7: e33933-10.1371/journal.pone.0033933.
Jaeger J, Manu , Reinitz J: Drosophila blastoderm patterning. Curr Opin Genet Dev. 2012, 22: 533-541. 10.1016/j.gde.2012.10.005.
Calvo E, Walter M, Adelman ZN, Jimenez A, Onal S, Marinotti O, James AA: Nanos (nos) genes of the vector mosquitoes, Anopheles gambiae, Anopheles stephensi and Aedes aegypti. Insect Biochem Mol Biol. 2005, 35: 789-798. 10.1016/j.ibmb.2005.02.007.
Juhn J, James AA: Oskar gene expression in the vector mosquitoes, Anopheles gambiae and Aedes aegypti. Insect Molec Biol. 2006, 15: 363-372. 10.1111/j.1365-2583.2006.00655.x.
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.
Akbari OS, Antoshechkin I, Amrhein H, Williams B, Diloreto R, Sandler J, Hay BA: The developmental transcriptome of the mosquito Aedes aegypti, an invasive species and major arbovirus vector. G3 (Bethesda). 2013, 3: 1493-1509. 2013.
Anderson WA, Spielman A: Incorporation of RNA and protein precursors by ovarian follicles of Aedes aegypti mosquitoes. J Submicrosc Cytol. 1973, 5: 181-198.
Mathew G, Rai KS: The structure and formation of egg membranes in Aedes aegypti. Internat J Insect Morphol Embryol. 1975, 4: 369-380. 10.1016/0020-7322(75)90037-9.
Raikhel AS, Lea AO: Control of follicular epithelium development and vitelline envelope formation in the mosquito; role of juvenile hormone and 20-hydroxyecdysone. Tissue Cell. 1991, 23: 577-5791. 10.1016/0040-8166(91)90015-L.
Gigliotti S, Graziani F, De Ponti L, Rafti F, Manzi A, Lavorgna G, Gargiulo G, Malva C: Sex-, tissue-, and stage-specific expression of a vitelline membrane protein gene from region 32 of the second chromosome of Drosophila melanogaster. Dev Genet. 1989, 10: 33-41. 10.1002/dvg.1020100106.
Andrenacci D, Cernilogar FM, Taddel C, Rotoli D, Cavaliere V, Graziani F, Gargiulo G: Specific domains drive VM32E protein distribution and integration in Drosophila eggshell layers. J Cell Sci. 2001, 114: 2819-2829.
Merzendorfer H, Zimoch L: Chitin metabolism in insects: structure, function and regulation of chitin synthases and chitinases. J Exp Biol. 2003, 206: 4393-4412. 10.1242/jeb.00709.
Moreira MF, Dos Santos AS, Marotta HR, Mansur JF, Ramos IB, Machado EA, Souza GH, Eberlin MN, Kaiser CR, Kramer KJ, Muthukrishnan S, Vasconcellos AM: A chitin-like component in Aedes aegypti eggshells, eggs and ovaries. Insect Biochem Mol Biol. 2007, 37: 1249-61. 10.1016/j.ibmb.2007.07.017.
Cohen E: Chitin synthesis and inhibition: a revisit. Pest Manag Sci. 2001, 57: 946-950. 10.1002/ps.363.
Arnold K, Brydon LJ, Chappell LH, Gooday GW: Chitinolytic activities in Heligmosomoides polygyrus and their role in egg hatching. Mol Biochem Parasitol. 1993, 58: 317-323. 10.1016/0166-6851(93)90054-2.
Biessmann H, Walter MF, Dimitratos S, Woods D: Isolation of cDNA clones encoding putative odourant binding proteins from the antennae of the malaria-transmitting mosquito, Anopheles gambiae. Insect Mol Biol. 2002, 11: 123-132. 10.1046/j.1365-2583.2002.00316.x.
Calvo E, Mans BJ, Ribeiro JM, Andersen JF: Multifunctionality and mechanism of ligand binding in a mosquito anti-inflammatory protein. Proc Natl Acad Sci. 2009, 106: 3728-3733. 10.1073/pnas.0813190106.
Rothemund S, Liou YC, Davies PL, Krause E, Sönnichsen FD: A new class of hexahelical insect proteins revealed as putative carriers of small hydrophobic ligands. Structure. 1999, 15: 1325-1332.
Manoharan M, Ng Fuk Chong M, Vaïtinadapoulé 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-80. 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.
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.
Fukuda N, Yomogida K, Okabe M, Touhara K: Functional characterization of a mouse testicular olfactory receptor and its role in chemosensing and in regulation of sperm motility. J Cell Sci. 2004, 15: 5835-5845.
da Silva AL C, Kojin BB, Marinotti O, James AA, Capurro ML: Expression and accumulation of the two-domain odorant binding protein AaegOBP45 in the ovaries of blood-fed Aedes aegypti. Parasites and Vectors. 2013, 6: 364-10.1186/1756-3305-6-364.
Gaunt MW, Miles MA: An insect molecular clock dates the origin of the insects and accords with palaeontological and biogeographical landmarks. Mol Biol Evol. 2002, 19: 748-761. 10.1093/oxfordjournals.molbev.a004133.
Han Q, Li G, Li J: Chorion peroxidase-mediated NADH/O(2) oxidoreduction cooperated by chorion malate dehydrogenase-catalyzed NADH production: a feasible pathway leading to H(2)O(2) formation during chorion hardening in Aedes aegypti mosquitoes. Biochim Biophys Acta. 2000, 1523: 246-253. 10.1016/S0304-4165(00)00129-X.
Pan C, Zhou Y, Mo J: The clone of laccase gene and its potential function in cuticular penetration resistance of Culex pipiens pallens to fenvalerate. Pestic Biochem Physiol. 2009, 93: 105-111. 10.1016/j.pestbp.2008.12.003.
Wu X, Zhan X, Gan M, Zhang D, Zhang M, Zheng X, Wu Y, Li Z, He A: Laccase2 is required for sclerotization and pigmentation of Aedes albopictus eggshell. Parasitol Res. 2013, 112: 1929-1934. 10.1007/s00436-013-3349-8.
Shibata T, Ariki S, Shinzawa N, Miyaji R, Suyama H, Sako M, Inomata N, Koshiba T, Kanuka H, Kawabata S: Protein crosslinking by transglutaminase controls cuticle morphogenesis in Drosophila. PLoS One. 2010, 5: e13477-10.1371/journal.pone.0013477.
Claycomb JM, Benasutti M, Bosco G, Fenger DD, Orr-Weaver TL: Gene amplification as a developmental strategy: isolation of two developmental amplicons in Drosophila. Dev Cell. 2004, 6: 145-155. 10.1016/S1534-5807(03)00398-8.
Fakhouri M, Elalayli M, Sherling D, Hall JD, Miller E, Sun X, Wells L, LeMosy EK: Minor proteins and enzymes of the Drosophila eggshell matrix. Dev Biol. 2006, 293: 127-141. 10.1016/j.ydbio.2006.01.028.
Fischer S, Alem IS, De Majo MS, Campos RE, Schweigmann N: Cold season mortality and hatching behavior of Aedes aegypti L. (Diptera: Culicidae) eggs in Buenos Aires City, Argentina. J Vector Ecol. 2011, 36: 94-99. 10.1111/j.1948-7134.2011.00145.x.
Keramaris KE, Margaritis LH, Zografou EN, Tsiropoulos GJ: Egg laying suppression in Drosophila melanogaster (Diptera:Drosophilidae) and Dacus (Bractocera) oleae (Diptera: Tephritidae) by phloroglucinol, a peroxidase inhibitor. Bull Entomol Res. 1996, 86: 369-375. 10.1017/S0007485300034957.
Willardsen JA, Dudley DA, Cody WL, Chi L, McClanahan TB, Mertz TE, Potoczak RE, Narasimhan LS, Holland DR, Rapundalo ST, Edmunds JJ: Design, synthesis, and biological activity of potent and selective inhibitors of blood coagulation factor Xa. J Med Chem. 2004, 47: 4089-4099. 10.1021/jm0497491.
Huang P, Ramphal J, Wei J, Liang C, Jallal B, McMahon G, Tang C: Structure-based design and discovery of novel inhibitors of protein tyrosine phosphatases. Bioorg Med Chem. 2003, 11: 1835-1849. 10.1016/S0968-0896(03)00039-7.
Meola RW, Dean SR, Meola SM, Sittertz-Bhatkar H, Schenker R: Effect of lufenuron on chorionic and cuticular structure of unhatched larval Ctenocephalides felis (Siphonaptera: Pulicidae). J Med Entomol. 1999, 36: 92-100.
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.
O’Brochta DA, Handler AM: Perspectives on the state of insect transgenics. Adv Exp Med Biol. 2008, 627: 1-18. 10.1007/978-0-387-78225-6_1.
Fu G, Lees RS, Nimmo D, Aw D, Jin L, Gray P, Berendonk TU, White-Cooper H, Scaife S, Kim Phuc H, Marinotti O, Jasinskiene N, James AA, Alphey L: Female-specific flightless phenotype for mosquito control. Proc Natl Acad Sci U S A. 2010, 107: 4550-4554. 10.1073/pnas.1000251107.
Robert MA, Okamoto K, Lloyd AL, Gould F: A reduce and replace strategy for suppressing vector-borne diseases: insights from a deterministic model. PLoS One. 2013, 8: e73233-10.1371/journal.pone.0073233.
Wiśniewski JR, Zougman A, Nagaraj N, Mann M: Universal sample preparation method for proteome analysis. Nat Methods. 2009, 6: 359-362. 10.1038/nmeth.1322.
Rappsilber J, Mann M, Ishihama Y: Protocol for micro-purification, enrichment, pre-fractionation and storage of peptides for proteomics using StageTips. Nat Protoc. 2007, 2: 1896-1906. 10.1038/nprot.2007.261.
Bonizzoni M, Dunn WA, Campbell CL, Olson KE, Marinotti O, James AA: Strain variation in the transcriptome of the dengue fever vector, Aedes aegypti. G3 (Bethesda). 2012, 2: 103-114. 2012.
Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.
Schmittgen TD, Livak KJ: Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008, 3: 1101-1118. 10.1038/nprot.2008.73.
Bailey TL, Bodén M, Buske FA, Frith M, Grant CE, Clementi L, Ren J, Li WW, Noble WS: MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009, 37: W202-W208. 10.1093/nar/gkp335.
Simossis VA, Heringa J: PRALINE: a multiple sequence alignment toolbox that integrates homology-extended and secondary structure information. Nucleic Acids Res. 2005, 33: W289-W294. 10.1093/nar/gki390.
Gasteiger E, Hoogland C, Gattiker A, Duvaud S, Wilkins MR, Appel RD, Bairoch A: Protein Identification and Analysis Tools on the ExPASy Server. The Proteomics Protocols Handbook. Edited by: Walker JM. 2005, Totowa, New Jersey: Humana Press, 571-607.
Kyte J, Doolittle R: A simple method for displaying the hydropathic character of a protein. J Mol Biol. 1982, 157: 105-132. 10.1016/0022-2836(82)90515-0.
We thank Judy Coleman and Aniko Fazenas for assistance with mosquito rearing. The research was supported by grants from the NIH NIAID (AI29746).
We declare that we have no significant competing financial, professional or personal interests that might have influenced the performance or presentation of the work described in this manuscript.
OM, PDG and AAJ conceived and designed the experiments. OM, TN, BBK, SPC, BN, JJ, RCL, PM, MFW and PDG performed experiments. OM, TN, JJ, RCL, XJ, ZT, PDG and AAJ analyzed the data. TN, XJ, JJ, ZT and PDG contributed reagents/materials/analysis tools. OM, JJ, PDG, ZT and AAJ wrote the manuscript. All authors read and approved the final manuscript.