SUMMARY
Dormancy is a crucial adaptation allowing insects to withstand harsh environmental conditions. The pre-programmed developmental arrest of diapause is a form of dormancy that is distinct from quiescence, in which development arrests in immediate response to hardship. Much progress has been made in understanding the environmental and hormonal controls of diapause. However, studies identifying transcriptional changes unique to diapause, rather than quiescence, are lacking, making it difficult to disentangle the transcriptional profiles of diapause from dormancy in general. The Asian tiger mosquito, Aedes albopictus, presents an ideal model for such a study, as diapausing and quiescent eggs can be staged and collected for global gene expression profiling using a newly developed transcriptome. Here, we use RNA-Seq to contrast gene expression during diapause with quiescence to identify transcriptional changes specific to the diapause response. We identify global trends in gene expression that show gradual convergence of diapause gene expression upon gene expression during quiescence. Functionally, early diapause A. albopictus show strong expression differences of genes involved in metabolism, which diminish over time. Of these, only expression of lipid metabolism genes remained distinct in late diapause. We identify several genes putatively related to hormonal control of development that are persistently differentially expressed throughout diapause, suggesting these might be involved in the maintenance of diapause. Our results identify key biological differences between diapausing and quiescent pharate larvae, and suggest candidate pathways for studying metabolism and the hormonal control of development during diapause in other species.
INTRODUCTION
Seasonal transitions require insects to respond to harsh environmental changes in order to survive. Diapause is an alternative developmental program that is initiated in response to a token stimulus, often photoperiod, that occurs well in advance of physiologically limiting environmental factors. Physiological changes during diapause result in developmental arrest, metabolic restructuring and stress tolerance, which allows insects to withstand seasonally occurring environmental insults, such as the harsh conditions of winter (Tauber and Tauber, 1976). For many insect species, developmental suppression continues after the physiological limitation has been lifted, until specific environmental changes or endogenous processes lead to diapause termination (Koštál, 2006).
Because of the paramount adaptive importance of diapause for insect survival during seasonal change, there has been a sustained interest in diapause physiology, which has led to the discovery of many of its environmental and hormonal controls (Denlinger, 2002; Denlinger et al., 2005). However, knowledge of the molecular regulation of diapause is only beginning to be thoroughly explored, in part because the stage in which diapause is expressed varies among insects, complicating efforts to identify common mechanisms of regulation (Denlinger, 2002). Additionally, the traditional model insect of choice, Drosophila melanogaster, has only a weak diapause phenotype (Emerson et al., 2009; Schmidt et al., 2005; but see Williams et al., 2010), minimizing the utility of this classical system for research on the molecular basis of diapause. Global gene expression profiling using microarrays or RNA sequencing (RNA-Seq) is increasingly being applied to ‘non-model’ insects – that are nonetheless excellent experimental systems for diapause – and has enabled substantial progress in documenting important transcriptional changes throughout diapause (Bao and Xu, 2011; Emerson et al., 2010; Poelchau et al., 2013b; Ragland et al., 2010; Ragland et al., 2011).
In contrast to diapause, quiescence is an alternative form of insect dormancy, in which physiological processes halt in immediate response to the reduction of an environmental, physiologically limiting factor (Hand and Podrabsky, 2000; Koštál, 2006). Once the environmental factor returns to non-limiting levels, normal activity is resumed. Both diapause and quiescence present important adaptations to avoid environmental exigencies, but there are important distinctions between these two forms of dormancy. While quiescence occurs in immediate response to an unpredictable environmental change, diapause is induced in advance of seasonally recurring changes, has an extended preparatory period, often over more than a generation, and the developmental arrest of diapause cannot be broken by an external stimulus until diapause has terminated. Metabolism is depressed in both types of dormancy, but the mode of depression can differ (Hand and Podrabsky, 2000). Diapausing and quiescent individuals can also have other important phenotypic differences that likely relate to different strategies of energy metabolism, such as differences in lipid content, the time interval required for achieving metabolic depression, and desiccation resistance (e.g. Hand and Podrabsky, 2000; Reynolds et al., 2012; Urbanski et al., 2010).
Gene expression profiling studies that have sought to identify transcriptional distinctions of diapause have usually characterized transitions between diapause stages (i.e. pre-diapause versus diapause versus post-diapause) (Bao and Xu, 2011; Emerson et al., 2010; Ragland et al., 2010; Ragland et al., 2011; but see Reynolds and Hand, 2009b). These studies have revealed important insights into the degree and mode of metabolic depression and developmental arrest, and accompanying physiological changes, in diapausing insects. However, an equally important question is how diapause is transcriptionally distinct from quiescence, and how this distinction changes throughout the dynamic progression of diapause (i.e. from early to late stages of developmental arrest). This comparison may yield alternative insights into the molecular distinctions of diapause that could be missed in other experimental designs, because two states of dormancy are being contrasted.
The Asian tiger mosquito, Aedes albopictus (Skuse), presents an excellent model system to identify molecular components of diapause. Aedes albopictus is a highly invasive vector species (Benedict et al., 2007) that enters diapause as a pharate larva within the chorion of the egg (Mori et al., 1981; Wang, 1966). Temperate populations of A. albopictus undergo a photoperiodic diapause in which a ‘short-day’ photoperiod experienced during the maternal pupal and adult stages stimulates the production of offspring destined for diapause. How diapause terminates in A. albopictus is not clear, but a certain period of time must elapse before diapause is broken [up to several months (Pumpuni, 1989), and in our laboratory, ca. 60 days], and this period can be influenced by temperature and photoperiod (Pumpuni, 1989). Like many other insect species (Hodek, 1996; Tauber et al., 1986), in A. albopictus diapause termination is followed by a period of post-diapause quiescence, which ends when environmental conditions that are favorable for direct development (i.e. immersion in water and high temperatures) stimulate direct development (i.e. hatching of the pharate larva from the egg). In contrast, maternal mosquitoes that experience a ‘long-day’ photoperiod oviposit eggs capable of quiescence: once embryonic development is complete, in the absence of a hatch stimulus, fully developed pharate larvae remain dormant within the eggs. This state of dormancy is distinct from diapause, because quiescent larvae will immediately hatch once the appropriate stimulus is received [e.g. flooding (reviewed in Hawley, 1988)]. Because the two types of dormancy are easily induced in A. albopictus, staged quiescent and diapause eggs can be easily matched and gene expression compared (Fig. 1). Additionally, there are substantial genomic resources available for expression studies in A. albopictus: while as yet there is no genome sequence available, the genome and accompanying annotations of the closely related Aedes aegypti (Nene et al., 2007) provide a powerful resource for global gene expression studies of A. albopictus (Poelchau et al., 2013a; Poelchau et al., 2013b; Poelchau et al., 2011).
Here, we identify global gene expression differences between diapausing and quiescent A. albopictus at three separate time points representing early, middle and late diapause. We find that gene expression patterns converge over time between the two states of dormancy to a quiescence expression profile. We identify key metabolic distinctions between diapause and quiescence that are important early in diapause, of which only differences in lipid metabolism remain throughout the course of diapause. Finally, we identify several genes with putative hormonal functions that are implicated throughout diapause, suggesting future avenues of investigation of the hormonal control of diapause maintenance.
MATERIALS AND METHODS
Experimental design
The experimental design (Fig. 1), insect rearing and RNA extraction have been described in a previous paper (Poelchau et al., 2013a). Tissue was generated from a laboratory F13 A. albopictus strain collected from Manassas, VA, USA. Larvae were reared at 21°C, ca. 80% relative humidity and a 16 h:8 h light:dark photoperiod until pupation (see Armbruster and Hutchinson, 2002; Armbruster and Conn, 2006). At pupation, mosquitoes were transferred to adult cages maintained under either diapause-inducing, short-day (D; 8 h:16 h light:dark) or non-diapause-inducing, long-day photoperiod treatments (ND; 16 h:8 h light:dark). We established four separate 9.5 liter adult cages (biological replicates) for each photoperiod treatment (D, ND) with ca. 100 mosquitoes per cage. Females were blood-fed on a human host 9–16 days after eclosion, and ca. 7 and 14 days thereafter. Egg collection to provide pharate larvae for RNA extraction and diapause measurements began 3 days after the first blood-feeding. Females were allowed to oviposit into a small brown jar lined with unbleached seed germination paper and half-filled with ca. 50 ml deionized water, which was placed into each cage 6–7 h after lights on. Egg papers were removed and replaced every 24 h for 26 days. Collected eggs were slowly air-dried on the papers 72 h after removal, and kept at 80% relative humidity, a 8 h:16 h light:dark cycle and 21°C until further use. We note that diapausing pharate larvae may be more metabolically active at 21°C than at lower temperatures, but this temperature is still ecologically relevant, because diapausing pharate larvae spend considerable time at higher temperatures during the fall before the onset of winter. Additionally, we are able to eliminate confounding effects of temperature by maintaining pharate larvae from all treatments at the same temperature. Eggs designated for RNA extraction were snap-frozen in liquid nitrogen and stored at −80°C at 11 (early diapause), 21 (mid-diapause) or 40 (late diapause or post-diapause quiescence) days post-oviposition (dpov; counted from the start of the oviposition period). Frozen eggs from each photoperiod, time point and replicate were ground in TRI Reagent (Sigma-Aldrich, St Louis, MO, USA), followed by RNA extraction according to the manufacturer's instructions. DNA was removed using Turbo-DNAfree (Applied Biosystems/Ambion, Austin, TX, USA). Three biological replicates from each photoperiod and development stage were chosen from the four available replicates based on RNA quality and quantity, as measured on an RNA chip (Bioanalyzer 2100, Agilent Technologies, Santa Clara, CA, USA). Only two biological replicates for 40 dpov pharate larvae reared on an ND photoperiod were chosen because of low RNA quality in the remaining replicates. Specific libraries chosen for each time period are listed in supplementary material Table S1. Incubator malfunction resulted in temperature irregularities for some 40 dpov eggs (ca. 4°C fluctuations on three consecutive days), but we discarded eggs scheduled for snap-freezing on these days. Furthermore, these temperature fluctuations should not result in systematic differences in gene expression between ND and D treatments, because ND and D eggs were stored together and experienced the same environmental conditions throughout the experiment.
For diapause incidence measurements, for each biological replicate, 14- to 28-day-old eggs were hatched, the number of hatched larvae recorded, and the egg papers with remaining, un-hatched eggs re-dried. This procedure was repeated twice, after 7 and 14 days. The remaining eggs were bleached (Trpis, 1970) to visualize and record the number of embryonated but unhatched (diapause) eggs. Diapause incidence was calculated as: (no. embryonated unhatched eggs)/(no. hatched eggs + no. embryonated unhatched eggs) (Urbanski et al., 2012). Percent embryonation was calculated as: (no. embryonated unhatched eggs + no. hatched eggs)/total no. eggs.
Sequence assembly and annotation
Sequencing, assembly and annotation are described in detail in Poelchau et al. (Poelchau et al., 2013a). Briefly, paired-end, barcoded Illumina mRNA-Seq libraries were constructed from each of the 17 RNA samples, and a proportion of each library was sequenced on three lanes on an Illumina HiSeq 2000 sequencer by the University of Maryland Genomics Institute. Cleaned reads were assembled into contigs after digital normalization (C. T. Brown, A. Howe, Q. Zhang, A. B. Pyrkosz and T. H. Brom, unpublished, arXiv:1203.4802) using Velvet (Zerbino and Birney, 2008) and Oases (Schulz et al., 2012). Contigs were then merged with two previous assemblies (Poelchau et al., 2011; Poelchau et al., 2013b) using a reference-based assembly approach outlined in Poelchau et al. (Poelchau et al., 2013a). Resulting contigs were annotated based on protein models from A. aegypti, Culex quinquefasciatus, Anopheles gambiae and D. melanogaster, and based on the A. aegypti genome sequence (Nene et al., 2007). Raw reads are available in NCBI's short read archive under accession number SRA063587, and the assembly can be downloaded at http://www.albopictusexpression.org/?q=data.
Gene expression analysis
Transcriptome assemblies without a genomic reference will generate redundant contigs for each identified gene model, in part because of allelic variation and/or alternative splicing. To account for this redundancy in our gene expression calculations, we used the program RSEM v.1.2.0 (Li and Dewey, 2011) to generate composite gene expression measures for each identified gene model. We mapped cleaned read pairs to the A. albopictus transcriptome using the program's default parameters.
Read counts were processed with the program edgeR (Robinson and Oshlack, 2010) in the R software environment (www.r-project.org). Only contigs with annotations to proteins or A. aegypti genome features based on gene set AaegL1.2 (as opposed to un-annotated A. aegypti genome sequence) were used in all subsequent analyses, as these required functional annotations. Additionally, gene models with fewer than two counts per million reads across all libraries were removed from the analysis, as these are not likely to show statistically significant differential expression (see Robinson et al., 2010). Read counts were TMM normalized (Robinson and Oshlack, 2010), which accounts for library size and expression bias, and log2 fold-change and its significance was calculated for each gene model between D and ND libraries for each time point (11, 21 and 40 dpov). We classified a gene as differentially expressed (DE) if its absolute log2 fold-change was greater than 0.5, with a Benjamini–Hochberg-corrected P<0.05. Previous RNA-Seq studies performed in our laboratory using the same A. albopictus strain, sequencing center and normalization methods show strong congruence with qRT-PCR results (Poelchau et al., 2013b), and RNA-Seq expression data have repeatedly been shown to produce accurate gene expression estimates, given proper normalization (Bullard et al., 2010; Feng et al., 2010; Fu et al., 2009; Poelchau et al., 2013b).
A distance matrix of gene expression patterns (R function dist) was summarized using multi-dimensional scaling (R function cmdscale) after transformation for linear modeling via the function voom in limma (Smyth, 2004; Smyth, 2005). Standardized expression patterns of all DE genes were also visualized as Z-scores in heat maps generated by hierarchical clustering (function hclust in R). Variability of gene expression within all D and ND DE genes, calculated as coefficients of variation (CV), was assessed via a Wilcoxon signed-rank test. All expression information is available at http://www.albopictusexpression.org/?q=data.
Gene Ontology and Kegg pathway enrichment analyses
Global gene expression data sets, such as those derived from RNA-Seq experiments, can provide insights into changes involving functionally related groups of genes that underlie specific physiological processes, e.g. Gene Ontology (GO) categories (Ashburner et al., 2000) or Kegg pathways (Kanehisa and Goto, 2000; Kanehisa et al., 2012). We asked whether these functional groups were over-represented among DE genes at each time period using the program GOseq, which corrects enrichment analyses for biases arising from variable transcript lengths in RNA-Seq data sets (Young et al., 2010). We also performed the same analysis for genes that were DE throughout all three time points in order to identify functional groups of genes that were differentially expressed throughout diapause. Generic GO Slim assignments for each gene model were downloaded from EnsemblMetazoa BioMart (Haider et al., 2009), and Kegg pathway assignments from http://www.genome.jp/kegg/. In addition to the suite of GO Slim categories and Kegg pathways, we manually composed gene lists representative of pathways or physiological processes with likely relevance for diapause in A. albopictus based on gene expression studies from other organisms (following Poelchau et al., 2013b): insulin signaling, which is instrumental for insect growth and metabolism (see Ragland et al., 2010; Wu and Brown, 2006); ecdysone signaling [from ‘molting’ genes in Brody (Brody, 1999)]; and heat shock proteins (Hsps), which are a subset of the gene ontology category ‘response to stress’ (GO:0006950). Uncorrected P-values from the GOseq analysis were Benjamini–Hochberg-corrected for multiple testing using the p.adjust function in limma (Smyth, 2004; Smyth, 2005). Functional groups with corrected P-values <0.05 and five or more DE genes were considered significantly enriched.
RESULTS
Diapause incidence
Diapause incidence of each biological replicate ranged from 87.5 to 100% in the diapause-inducing photoperiod treatments, and percent hatch ranged from 77.4 to 82.9% in the non-diapause treatment (supplementary material Table S1) (Poelchau et al., 2013a). Embryonation ranged from 82.9 to 98.9% across all replicates. Diapause incidence was not 100% for all replicates, indicating that a mixture of mostly diapause, but some quiescent, pharate larvae were sequenced in the D libraries at some time points. However, this is not likely to generate spurious results, but rather makes our analysis more conservative, as fewer genes are likely to be detected as DE.
Gene expression analysis
The number of genes that were DE between D and ND conditions decreased from early to middle to late diapause (Fig. 2, Table 1). Multi-dimensional scaling of the gene expression results clustered libraries generally by day post-oviposition and photoperiod treatment, although one D library at 40 dpov clustered with the two ND libraries (Fig. 3). The first MDS axis appeared to separate libraries by day post-oviposition (explaining 37.6% of the variation in the data), with the libraries clustering in chronological order, and tighter clustering occurring between the 21 and 40 dpov libraries. Also, the 11 dpov D libraries were located towards the ‘earlier’ side of the axis relative to the 11 dpov ND libraries. This is interesting, given previous observations from A. albopictus that suggest a developmental delay of embryos during diapause preparation (Poelchau et al., 2013b). The 21 and 40 dpov libraries, in contrast, show no such temporal separation between the D and ND treatments. The second axis, which explains 14.1% of the variation, roughly separated the libraries by photoperiod. Taken together, these results indicate that gene expression during A. albopictus quiescence and diapause converges over time. Normalized read counts, log-fold changes and their P-values, and descriptions of all genes in the data set are available in supplementary material Table S2.
We visualized gene expression of all DE genes to identify trends in their expression convergence over time. Heat maps of standardized expression scores suggested that change in expression over time is driven by change during diapause, not quiescence (Fig. 4). To determine whether the decrease in the number of DE genes over time was driven more by change in ND or D expression, we asked whether the coefficient of variation (CV) in gene expression across time periods differed between D and ND genes. D genes had higher CVs that differed significantly from ND genes (paired Wilcoxon-rank-sum test, P<2.2e–16; ND, mean CV=0.357, D, mean CV=0.453), indicating that gene expression change over time occurs more in D rather than in ND pharate larvae.
Gene Ontology and Kegg pathway enrichment analyses
The enrichment analyses overwhelmingly point to differential regulation of metabolic processes as a key distinction between diapause and quiescence in A. albopictus. The number of enriched processes decreased over time, reflecting the decreasing number of DE genes at each time point (Table 1). In addition, these categories converged in function: the category ‘lipid metabolism’ remained enriched across time periods, and was the only enriched category remaining at 40 dpov. We describe the categories at individual time points below.
A diversity of metabolic processes were enriched for DE genes at 11 dpov (Table 1). Several were related to amino acid metabolism, with most genes under-expressed in D conditions. In contrast, genes involved in carbohydrate metabolism, including pyruvate metabolism and glycolysis/gluconeogenesis, were primarily over-expressed. Genes involved in lipid metabolism had mostly higher expression under D conditions, although the direction of expression was mixed. The category ‘extracellular region’ was mainly comprised of lipid metabolism genes, and genes that contained chitin-binding domains. The categories ‘biological process’ and ‘metabolic pathways’ have a strong overlap in gene composition with other categories, and will not be discussed further.
At 21 dpov, the categories ‘lipid metabolism’ and ‘extracellular region’ were enriched for DE genes, where ‘extracellular region’ genes were primarily lipases and genes with chitin-binding domains. At 40 dpov, only ‘lipid metabolism’ remained enriched. Similar to 11 dpov, gene expression tended to be higher under D conditions for all of these categories at 21 and 40 dpov.
Twenty-five genes were DE throughout all three time points (Table 1), and all were expressed in the same direction continuously. The set of DE genes was enriched for genes in the category ‘extracellular region’, which contained two lipases, two triacyl-glycerol lipases and one conserved hypothetical protein containing conotoxin domains, indicating that most of the genes were involved in lipid metabolism. This is consistent with the enrichment of ‘lipid metabolism’ at all three time points (see above; Table 1). While not contributing to an enriched category, three genes related to juvenile hormone (JH) binding and metabolism were also DE: a putative JH-inducible protein, AAEL012680, under-expressed in D conditions; an ortholog of juvenile hormone esterase (JHE), AAEL005200, over-expressed; and a gene containing a JH binding protein domain, AAEL000500, over-expressed (supplementary material Table S3).
DISCUSSION
In this study, we contrast gene expression patterns throughout A. albopictus diapause and quiescence, two alternative developmental states for pharate first instar larvae, in order to identify transcriptional regulation of processes relevant to diapause. To our knowledge, this is the first study to identify global transcriptional components of diapause maintenance for pharate larval diapause, and the first RNA-Seq study to use staged comparisons of diapause versus quiescence to identify diapause-enriched transcripts. These comparisons allowed us to demonstrate on a global scale how gene expression patterns during A. albopictus diapause maintenance reflect biological distinctions of A. albopictus diapause from quiescence. Aedes albopictus eggs are refractory to hatching stimuli during diapause, yet return to a period of post-diapause quiescence. In turn, A. albopictus gene expression patterns demonstrate convergence of diapause upon quiescence over time (Figs 2, 4), rather than achieving a state entirely distinct from quiescence once diapause is broken. We note that because diapause was only measured once for each replicate between 14 and 28 dpov, the transcriptional changes between 11 and 40 dpov could be due to diapause termination leading to post-diapause quiescence in a proportion of individuals in the diapause group. Previous experiments suggest that most eggs terminate diapause between 30 and 60 dpov under the conditions utilized in this experiment (Pumpuni, 1989) (P.A.A., unpublished data). Thus, transcriptional convergence between the diapause and quiescence profiles could be driven by transcriptional change during diapause and/or by diapause loss itself. However, regardless of the mechanism, diminishing gene expression differences over time are likely to reflect a physiological convergence of diapause towards quiescence. This observation supports the model of diapause as a dynamic process, rather than a static condition (Denlinger, 2002; Koštál, 2006; Tauber et al., 1986).
Metabolic gene expression is the main transcriptional distinction between early diapause and quiescence (Table 1). This general result corroborates findings from many other gene expression studies of insect diapause, which have documented profound transcriptional changes related to metabolism in diapausing larvae, pupae and adult insects (e.g. Emerson et al., 2010; Ragland et al., 2010; Ragland et al., 2011; Reynolds et al., 2012), and many other physiological studies that have outlined various mechanisms of metabolic restructuring during diapause (reviewed in Hahn and Denlinger, 2007; Hahn and Denlinger, 2011; Hand et al., 2011). In general, diapausing animals survive extended periods of developmental arrest by increasing nutrient stores during diapause preparation, and reducing metabolism during developmental arrest. How these metabolic changes are achieved differs among species: insects vary in the degree that metabolism is suppressed during diapause (Chaplin and Wells, 1982; Denlinger et al., 1972; Reynolds and Hand, 2009a), and in the specific nutrient compositions that are stored and later utilized, although energy stores in form of triacylglycerides (Danks, 1987; Hahn and Denlinger, 2007), glycogen (Danks, 1987; Zhou and Miesfeld, 2009) and specialized storage proteins (Burmester, 1999; Denlinger et al., 2005) are common. Below, we discuss gene expression changes underlying the metabolism of different types of nutrient stores during diapause.
Carbohydrate metabolism
Glycolysis is the first step in generating metabolic energy from glucose. Gluconeogenesis reverses the glycolytic process, generating glucose from pyruvate, and uses many of the steps of glycolysis in reverse. Upregulation of the gluconeogenetic process has been implicated in previous studies of diapause gene expression (Baker and Russell, 2009; Emerson et al., 2010; Ragland et al., 2010; Ragland et al., 2011), and in a diapause context is considered consistent with reliance on anerobic metabolism (Hahn and Denlinger, 2011). At 11 dpov, both the glycolysis/gluconeogenesis pathway and the pyruvate metabolism pathway were enriched for DE genes, most of which had higher expression in diapause (Table 1). Pepck (phosphoenolpyruvate carboxykinase; AAEL000006, AAEL000080), which encodes a rate-limiting enzyme in gluconeogenesis, and a gapdh homolog (glyceraldehyde 3 phosphate dehydrogenase; AAEL016984) had higher expression under diapausing conditions. These results are consistent with upregulation of gluconeogenesis, which suggests a shift towards anerobic metabolism. In contrast, pyk (pyruvate kinase; AAEL012576, AAEL014913) also had higher expression, which should indicate reliance on glycolysis, because PyK converts phosphoenolpyruvate to the end-product of the glycolysis pathway, pyruvate. However, these results are not necessarily contradictory, as PyK can be inhibited via post-translational modification under fasting conditions (Feliú et al., 1976; Llorente et al., 1970).
In addition to A. albopictus, a diverse range of other organisms show upregulation of pepck during diapause in gene expression scans, such as Sarcophaga crassipalpis, Rhagoletis pomonella, Wyeomia smithii and Caenorhabditis elegans (Emerson et al., 2010; McElwee et al., 2006; Ragland et al., 2010; Ragland et al., 2011), suggesting that this enzyme may have a ubiquitous role in the metabolic restructuring of diapausing animals. In A. albopictus, pepck expression is high throughout diapause induction (Poelchau et al., 2011), preparation (Poelchau et al., 2013b) and early diapause. Collectively, our results on pepck suggest a reliance on anerobic metabolism in preparation for diapause and during early diapause that exceeds that of quiescent pharate larvae.
Amino acid metabolism
Amino acid metabolic pathways synthesize proteins, hormones and enzymes; they can also degrade amino acids to generate metabolic intermediates of glucose to be used in the citric acid cycle (Klowden, 2007). Amino acids are thought to mediate cold and desiccation resistance during diapause (e.g. Michaud and Denlinger, 2007) or to play a role in nutrient storage (Morgan and Chippendale, 1983). Several amino acid metabolic pathways were enriched at 11 dpov, mostly for genes with lower expression under D conditions (Table 1, supplementary material Table S4). Many of these genes are involved in glutamine, glycine and serine metabolism (supplementary material Table S4). These results suggest a downregulation of these pathways, which could result in (1) lower provisioning of the citric acid cycle with metabolic intermediates, consistent with a shift towards anaerobic metabolism, and (2) higher concentrations of amino acids because of decreased degradation. Consistent with this interpretation, preliminary metabolomics data from A. albopictus show higher levels of amino acids in diapausing versus non-diapausing eggs (leucine, serine, threonine, tyrosine, lysine and proline; data not shown). In general, these data point towards a key role of amino acids in early A. albopictus diapause that is consistent with increased anaerobic metabolism and increased cold and desiccation resistance.
Lipid metabolism
Lipids can serve as a fundamental energy source for diapausing insects, and are the primary fuel for embryonic development (Arrese and Soulages, 2010; Van Handel, 1993). Because of their high caloric content and water yield, they store energy more efficiently than carbohydrate-based sources (Hahn and Denlinger, 2011). Diapausing animals, which often do not feed and thus must rely on stored nutrients for survival, can be provisioned with higher lipid reserves, in particular triacylglycerides, than their non-diapause counterparts (Danks, 1987; Hahn and Denlinger, 2007; McElwee et al., 2006; Tauber et al., 1986). These stores can then be metabolized during diapause via lipases, which catalyze the hydrolysis of triacylglycerides, to generate energy. Diapausing 11 dpov pharate larvae were enriched for lipid metabolism genes (Table 1), especially genes involved in lipid store mobilization, such as lipases and hydrolases (supplementary material Table S4). Expression patterns of these genes were mixed; however, the majority of lipases, in particular the triacylglycerol lipases, were upregulated rather than downregulated, suggesting that diapausing A. albopictus pharate larvae metabolize lipid stores as an energy source at this stage. Reliance on lipid stores as an energy source during diapause is consistent with a previous physiological study of A. albopictus: 10- to 14-day-old A. albopictus eggs contained ~30% more total lipid than quiescent eggs, and pre-diapause embryos showed expression evidence of lipid storage relative to non-diapause embryos (Reynolds et al., 2012).
Lipid metabolism persisted as a distinct feature of all sampled diapause stages: lipid metabolism genes were enriched throughout all sampled time periods (Table 1), and remained primarily over-expressed under diapause conditions. This points to a consistent role of lipids as an energy store provisioning pharate larvae throughout diapause, relative to quiescence. Other studies in diapausing insects have indicated different temporal profiles of lipid metabolism across diapause: for example, in the adult diapause of the mosquito Culex pipiens, lipase expression is low in early diapause, then increases in late diapause (Sim and Denlinger, 2009). The cotton bollworm Helicoverpa armigera also downregulates lipase expression in early diapause, presumably to promote lipid storage for use as an energy source later in diapause (Bao and Xu, 2011). Our data indicate that at 11 dpov, lipids are already an important source of energy for diapausing A. albopictus, suggesting that the relative importance of this energy source is higher than for other species.
Hormone action during diapause maintenance
Insect hormones play a fundamental role in the control of development (Fraenkel, 1935; Riddiford, 1994; Wigglesworth, 1934) and diapause (Denlinger, 2002). The relative abundance of two major hormones, ecdysone and JH, during development dictates the developmental progression of the insect (Klowden, 2007; Riddiford, 1994). Changes in the relative and absolute abundance of JH or ecdysone are known to be important during the initiation, maintenance and termination of diapause in many insect species, but the nature of these changes will depend on the life-cycle stage of diapause developmental arrest (Denlinger, 2002). For example, ecdysteroids play a regulatory role in the pharate larval diapause of the gypsy moth Lymantria dispar (Lee and Denlinger, 1997; Lee et al., 1997), whereas JH mediates the hormonal control of adult diapause in C. pipiens (Readio et al., 1999; Spielman, 1974).
The mechanisms of hormonal control of pharate larval diapause in A. albopictus are unknown. Previous transcriptome analyses of A. albopictus suggested a role for ecdysteroid signaling during the preparatory stage of diapause (Poelchau et al., 2013b; Poelchau et al., 2011). In our analysis of A. albopictus developmental arrest, we identified three genes (out of a total of 25) with putative functions related to JH action that had consistent differential expression across all time points (supplementary material Table S3). This conspicuous pattern would suggest that JH, or its absence, has a role in the maintenance of A. albopictus diapause. Endogenous JH production begins in late embryonic development, and its presence is thought to be important for dorsal closure, first-instar larval cuticle formation, and differentiation of the midgut (reviewed in Riddiford, 1994). JH continues to be present during larval feeding, inter-molt and molting phases, and its titer rises before the molt to the next larval instar. JH levels can be influenced by the environment: for example, starvation can increase JH titers (Truman et al., 2006). The JH titer is a function of JH synthesis and degradation. The enzyme that degrades JH is JHE (Klowden, 2007). A JHE homolog (AAEL005200), which has been verified experimentally in A. aegypti (Bai et al., 2007), was over-expressed under diapausing conditions, which would suggest it functions to keep JH levels low throughout diapause. Consistent with this pattern, a putative JH-inducible protein (AAEL012680), which should increase in expression under higher JH levels, was under-expressed, suggesting lower JH levels in diapausing pharate larvae. In contrast, a gene containing multiple JH binding domains (AAEL000500), which generally function to transport JH and protect it from degradation by JHEs, was over-expressed. Therefore, it is difficult to deduce a mode of JH action in the maintenance of diapause in A. albopictus. Accordingly, an experiment using adult females reared under short day-lengths did not show conclusive effects of JH topical application on subsequent hatching rates (Pumpuni, 1989). However, our data from the pharate larval stage strongly suggest that further research into JH as a regulatory hormone of diapause maintenance is worthwhile.
An inspection of the ‘biological process’ enriched category revealed a list of 15 members of the cytochrome P450 family (supplementary material Table S4). One of the diverse functions of the cytochrome P450 family is steroid hormone biosynthesis (Miller, 1988). Four of the cytochrome P450s were DE at both 11 and 21 dpov; and one of these, cyp18a1, encodes an enzyme that inactivates steroid hormones in D. melanogaster; loss-of-function mutations in D. melanogaster cause an extended final larval instar and lethality during metamorphosis (Guittard et al., 2011). Drosophila melanogaster cyp18a1 is also homologous to C. elegans daf-9, which regulates dauer, larval growth and longevity (Gerisch et al., 2001). Cyp18a1 had lower expression in diapausing A. albopictus pharate larvae, which is intriguing, given that daf-9 loss-of-function mutants form constitutive dauer larvae (Gerisch et al., 2001; Jia et al., 2002). Because of its conspicuous expression pattern – lower expression during early and mid-diapause, and lack of differential expression late in diapause – and because of the documented function of related genes, this gene represents a promising candidate for future studies into the hormonal regulation of diapause maintenance in A. albopictus.
Stress resistance
Diapausing insects use various mechanisms to tolerate adverse environmental conditions, such as extreme cold, aridity and hypoxia (Denlinger, 2002; MacRae, 2010). Hsps, in particular Hsp70, are often, but not always, upregulated in diapausing insects as protection against cold injury (Hayward et al., 2005; Rinehart et al., 2007). For example, one of the few studies contrasting gene expression during diapause versus post-diapause quiescence found that Hsp70 and Hsp23 were upregulated during Sarcophaga crassipalpis diapause (Hayward et al., 2005). Interestingly, this study also found close parallels between diapause and post-diapause quiescence: expression of these genes continued at a high level after diapause was broken, but before adult development resumed. We did not find conspicuous evidence for diapause upregulation of Hsps in our analysis of A. albopictus. This result aligns with findings from other insect species (Rinehart et al., 2007), such as C. pipiens, where Hsp70 was not upregulated in diapausing adults, despite the fact that diapausing individuals were more cold tolerant (Rinehart et al., 2006; Rinehart et al., 2007). These results suggest that other protective measures against low temperatures, such as amino acid provisioning (see ‘Amino acid metabolism’, above) or synthesis of classic cryoprotectants, distinguish A. albopictus diapause from quiescence. Finally, several immune-related genes – for example, homologs of a putative cecropin anti-microbial peptide and of Gram-negative binding proteins 3 and 4 – were over-expressed at 11 dpov (supplementary material Table S4), suggesting that pathogen defense may be particularly important in early diapausing pharate larvae. Similarly, several immune-responsive genes had higher expression in diapausing S. crassipalpis (Ragland et al., 2010). This suggests that higher investment in pathogen defense may be a common strategy during insect diapause.
Conclusions
In this study, we use gene expression profiling to gain fundamental insights into the molecular and physiological distinctions between diapause and quiescence in A. albopictus. Very little is known about the molecular mechanisms of diapause maintenance in pharate larval diapause. Because of the increasing importance of A. albopictus as a disease vector (Benedict et al., 2007), understanding the molecular regulation of this crucial life history trait could potentially provide a platform for novel vector control strategies based on the genetic or chemical disruption of diapause. We find gradual convergence of global diapause gene expression patterns towards quiescence. Metabolic differences, which are the primary distinguishing factor between early diapause and quiescent gene expression, decline over time to only include small differences in lipid metabolism, likely the main source of energy for diapausing pharate larvae. The data also suggest a role for juvenile hormone, and a member of the cytochrome P450 family, in facilitating diapause maintenance. With our experimental design, we can effectively characterize diapause as a physiological state distinct from quiescence, and therefore identify more subtle and likely important components of diapause that would be missed if compared with actively developing first-instar larvae.
Acknowledgements
We thank the Armbruster, Elsik and Denlinger labs for helpful comments and suggestions on this work.
FOOTNOTES
FUNDING
This work was supported by the National Institutes of Health [grant no. 5R21AI081041-02 to P.A.A., C.G.E. and D.L.D.] and Georgetown University. Deposited in PMC for release after 12 months.
LIST OF SYMBOLS AND ABBREVIATIONS
REFERENCES
COMPETING INTERESTS
No competing interests declared.