Genome wide gene-expression analysis of facultative reproductive diapause in the two-spotted spider mite Tetranychus urticae
BMC Genomics volume 14, Article number: 815 (2013)
Diapause or developmental arrest, is one of the major adaptations that allows mites and insects to survive unfavorable conditions. Diapause evokes a number of physiological, morphological and molecular modifications. In general, diapause is characterized by a suppression of the metabolism, change in behavior, increased stress tolerance and often by the synthesis of cryoprotectants. At the molecular level, diapause is less studied but characterized by a complex and regulated change in gene-expression. The spider mite Tetranychus urticae is a serious polyphagous pest that exhibits a reproductive facultative diapause, which allows it to survive winter conditions. Diapausing mites turn deeply orange in color, stop feeding and do not lay eggs.
We investigated essential physiological processes in diapausing mites by studying genome-wide expression changes, using a custom built microarray. Analysis of this dataset showed that a remarkable number, 11% of the total number of predicted T. urticae genes, were differentially expressed. Gene Ontology analysis revealed that many metabolic pathways were affected in diapausing females. Genes related to digestion and detoxification, cryoprotection, carotenoid synthesis and the organization of the cytoskeleton were profoundly influenced by the state of diapause. Furthermore, we identified and analyzed an unique class of putative antifreeze proteins that were highly upregulated in diapausing females. We also further confirmed the involvement of horizontally transferred carotenoid synthesis genes in diapause and different color morphs of T. urticae.
This study offers the first in-depth analysis of genome-wide gene-expression patterns related to diapause in a member of the Chelicerata, and further adds to our understanding of the overall strategies of diapause in arthropods.
The seasonal pattern in temperate regions poses a major challenge for arthropod development and survival, and as a result the reproduction and growth of poikilothermic organisms are largely limited to the warmer part of the year. In addition, specific mechanisms have evolved to survive the extreme climatic conditions of winter. One of those mechanisms is a developmental arrest, often considered as the most distinctive characteristic of diapause in insects and mites. In many species, diapause is limited to one particular life stage, and examples have been documented varying from the embryonic stage to larval, pupal and adult stages . Besides developmental arrest, diapause is characterized by a suppression of metabolism , altered behavior , increased stress tolerance  and increased energy reserves . To overcome injury caused by low temperatures, diapausing species can also synthesize cryoprotectants like polyols [6, 7], heat shock proteins [4, 8], and thermal hysteresis proteins [8, 9]. At the molecular level, it is clear that diapause is not simply characterized by a decreased expression of a large number of genes, but that the climate-induced expression of a specific set of genes is crucial to regulate and accomplish the physiological adaptations mentioned above .
In spider mites (Acari: Tetranychidae), diapause is limited to the egg or female adult stage, depending on the species . Two varieties of true diapause, i.e. aestival and hibernal diapause have been described, of which the latter is by far the best studied. In the two-spotted spider mite Tetranychus urticae, a hibernal facultative reproductive diapause allows adult female mites to survive winter . Since this spider mite species is a highly polyphagous and difficult to control plant pest [11, 12], winter survival has important consequences for crop protection in temperate climates. Therefore, within the group of Acari, diapause is best studied in this species, mainly on the level of behavior and to some extent on the level of biochemistry [3, 13]. Adult female spider mites initiate diapause in response to the reduction of day length during juvenile development. The diapausing females of the green color morphs are characterized by a deeply orange color and suppression of the ovarian development . Next to photoperiod, other factors influence the onset of diapause such as temperature, food quality, food quantity  and predation risk [15, 16]. The incidence of diapause can vary greatly between the predominant natural color morphs, red and green, and different T. urticae populations. In some populations of the red color morphs, a high percentage of the population enters into diapause after stimulus, while others lack the ability . For the green form of T. urticae, this variability is profoundly examined and appears to be dependent on light:day regime [18, 19], the presence of host plants during hibernation period  and genetic variability [18, 21]. A recent classical genetic study reveals that in some populations, a single recessive allele is responsible for the non-diapause phenotype . Next to the change in color and developmental arrest, diapause in T. urticae is characterized by various biological and physiological changes. A morphological difference between diapause and non-diapause forms is found in the structure of the integumentary striae on the dorsal surface of the cuticle , that significantly reduce the evaporative surface and thus the rate of evaporation . At the physiological level, it is known that diapausing mites deposit fat in the body tissues  and further limit transpiration through the cuticle by closing the stigmata . Recently, it was experimentally confirmed that T. urticae diapausing females show an increased cold hardiness, as the supercooling point of acclimated diapausing female mites is 6°C lower than non-diapausing females .
Diapause in T. urticae is also accompanied by a number of behavioral changes: mites almost completely cease feeding, leave the host plant in search of hibernation sites and mate immediately after the last molting stage (teliochrysalis), but do not lay eggs . At the time of color change to the typical orange forms, the hindgut is emptied and mites become positively geotactic  and negatively phototactic to find shelter [29, 30].
Despite the fact that a number of studies have addressed the changes in morphology, physiology and behavior of diapausing T. urticae, the knowledge about the molecular modifications underlying and associated with diapause in spider mites is very limited. Recently, the draft genome of T. urticae was reported and a high quality annotation of genes is available . Exploiting this annotated genome sequence and a previously developed whole genome gene expression micro-array platform , we studied gene expression changes during diapause in T. urticae. Such genome wide transcriptional changes associated with diapause were previously studied in insect species such as the flesh fly Sarcophaga crassipalpis and compared with the nematode Caenorhabditis elegans and the fruit fly Drosophila melanogaster. This comparison reveals that there may be diverse molecular mechanisms, but also that a number of adaptations related to diapause are broadly conserved. In this study, we extend the knowledge on molecular mechanisms underlying diapause in arthropods to the spider mite T. urticae, a member of the Chelicerata that diverged more than 450 Mya from other arthropod lineages such as insects and crustaceans .
The green strain LS-VL of T. urticae was originally collected in October 2000 near Ghent  and was ever since maintained on potted kidney bean plants Phaseolus vulgaris L. var. Prelude in controlled conditions at 24 ± 0.5°C and 60% RH with a 16:8h light:dark photoperiod (standard incubation conditions). For induction of diapause, 500 LS-VL adult females were transferred to a single bean plant to allow the deposition of eggs and were kept in standard conditions for 4 days until the larvae hatched. Subsequently, this bean plant and 7 other uninfected plants were held under diapause inducing conditions at 17 ± 0.5°C, 80% RH with 8:16h light:dark photoperiod. After 3 weeks, diapausing females with a distinguishing orange color (Figure 1C) and non-diapausing females (Figure 1A) were collected for RNA extraction. Other strains (London, MR-VL and Tu-SB9) used in this study for qPCR experiments were kept in the same standard incubation conditions. The London strain is a green morph of T. urticae that was originally collected in the Vineland region, Ontario, Canada, and originates from the culture used in the T. urticae genome project . MR-VL is a red morph and a well characterized multi-resistant strain . Tu-SB9 was collected in Crete, Greece in 2006 and is a red morph (Figure 1B) of T. urticae.
RNA preparation for microarray experiments
In order to extract total RNA, four replicates of 250 adult female mites with distinguishing characteristics of diapause (deep orange color, Figure 1C) were collected and homogenized according to the protocol of RNeasy mini kit (Qiagen). The 4-fold replicated control groups consisted of 100 adult females that did not enter diapause under similar diapause inducing conditions (Figure 1A). An additional washing step with one volume of chloroform:iso-amylalcohol (24:1) was performed to decrease the amount of pigments in the supernatant . The quality and quantity of the total RNA was analyzed by a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies) and by running an aliquot on a 1% agarose gel.
RNA preparation for qPCR
For qPCR experiments, total RNA was extracted from 3 replicates of 100 adult females of the strains LS-VL, London, MR-VL and Tu-SB9. After homogenization, RNA was extracted according to the protocol of the RNeasy mini kit (Qiagen). The quality and quantity of total RNA was analyzed as abovementioned.
A custom Sureprint genome wide G3 Gene Expression 8x60K microarray (GEO Platform accession GPL16890) was designed using the Agilent eArray platform (Agilent Technologies) based on the T. urticae gene annotation file frozen in April 2011, including coding sequences of 18,217 predicted unigenes, as previously described [12, 38]. The slide layout consisted of eight arrays per slide.
Microarray preparation, hybridization and analysis
One hundred nanograms of RNA of non-diapausing (4 replicates) and diapausing mites (4 replicates) were used to generate Cy3- and Cy5-labeled cRNA respectively, using the Agilent Low Input Quick Amplification labeling kit (version 6.5, Agilent Technologies). RNA spike-in controls (Agilent Technologies) were added to diluted aliquots of RNA samples before cRNA synthesis. The labeled cRNA was purified with the RNeasy mini kit (Qiagen). The dye content and concentration of cRNA was measured by NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies). Cy3- and Cy5-labeled cRNAs were pooled and hybridized using the Gene Expression Hybridization Kit (Agilent Technologies) for 17 h in a rotating hybridization oven at 10 rpm and 65°C. After hybridization, slides were washed using the Gene Expression Wash Buffer kit (Agilent Technologies), treated with Stabilization and Drying solution (Agilent Technologies), protected by an Ozone-Barrier cover (Agilent Technologies) until scanned by an Agilent Microarray High Resolution Scanner with default settings for 8 × 60K G3 microarrays. Data were normalized by the Agilent Feature Extraction software version 10.5 (Agilent Technologies) with default parameter settings for gene expression two-color microarrays (protocol GE2_107_SEP09) and transferred to GeneSpring GX 11.0 software (Agilent Technologies) for further statistical evaluation. Next, probes were flag filtered (only probes that had flag-value ‘detected’ in 50% of the replicates were retained) and linked to the most recent annotation file (September 2011, ) using the ‘Create New Gene-Level Experiment”-option. Genes with a Benjamini-Hochberg false discovery rate (FDR) corrected p-value < 0.05 and with an absolute fold change (FC) ≥ 2 were considered as differentially expressed. The microarray data reported in this paper have been deposited in the Gene Expression Omnibus (GEO) (accession number: GSE48858).
Differentially expressed genes were exported to Blast2GO software v.2.6.3 . This tool enables homology searches based on Protein Basic Local Alignment Search Tool (BLASTp) against the NCBI non-redundant protein database using an E-value cut-off of 1e-15. Subsequently, Blast2GO mapping performs different steps to link all best BLAST hits to information stored in the Gene Ontology (GO) database. All these mapping results are associated to an Evidence Code which provides information about the quality of mapping. After mapping, the results were subjected to GO annotation whereas GO terms were selected from the GO pool with a threshold of 1e-15 and assigned to the sequences. Further annotation was done using InterPro annotation in Blast2GO and corresponding GO terms were transferred to the sequences and merged with already existent GO terms. For the complete genome, the functional analysis stopped here. But for the genes that were differentially expressed in the microarray, GO terms were modulated using the Augment annotation tool by Annex  followed by GO-Slim which is a reduced version of the Gene Ontology database that contains a selected number of relevant nodes.
Gene ontology terms enrichment analysis
Within Blast2GO, an analysis was performed on gene function information, by the statistical assessment of GO term enrichment in a group of genes of interest when compared with a reference group i.e. to assess the functional differences between two sets of functional annotations. This was integrated in Blast2GO by Gossip  that performs a Fisher’s Exact test in combination with a Benjamini-Hochberg False discovery rate (FDR) correction for multiple testing. The option “Two-Tailed” was selected to perform a test which considers both directions as extreme and this translates to over- and under-represented Gene Ontology functions in the test-set compared to a reference set. The result of this test returns a list of significant GO terms ranked by p-values. GO terms with a Fisher Exact test FDR corrected p-values < 0.01 were considered as statistically significant and in order to reduce the size of the result-set to existing, statistically significant, child GO terms, the function “reduce to most specific terms” was used. This function enables to simplify the outcome of GSEA in case of a very large list of enriched GO terms by removing parent terms of already existing, statistically significant, child GO terms that represent the same functional concept but at different levels of specificity.
To identify statistically over-represented functions in our dataset, a Fisher’s exact test was carried out with the selected differentially expressed genes compared to the completely functional annotated genome of T. urticae. Between the sets of up-and down-regulated genes, a second enrichment test was performed to detect different biological roles of these sets during diapause.
Microarray confirmation by qPCR
In order to validate the microarray results, gene specific primers were designed for 13 differentially expressed T. urticae genes (7 up- and 6 down-regulated genes) using Primer 3 v0.4.0. . T. urticae genes coding for actin, ubiquitin and RP49 were used as housekeeping genes after being tested for their suitability as reference genes by comparing expression patterns between diapausing and non-diapausing forms (Additional file 1). Total RNA was extracted as described above and cDNA was synthesized with 2 μg of total RNA using the Maxima First Strand cDNA synthesis kit for RT-PCR (Fermentas Life Sciences). Three biological and two technical replicates were used to conduct these experiments and no-template-controls were added to exclude sample contamination. All qPCR reactions were carried out with the thermal cycler Mx3005P (Stratagene). Reactions were prepared with Maxima SYBR Green qPCR/Master Mix following the manufacturer’s instructions (Fermentas Life Sciences). The reactions were run with the following protocol: initial denaturation at 95°C for 10s followed by 35 cycles of 95°C for 15s, 55°C for 30s, 72°C for 30s. At the end of these cycles, a melting curve (from 55°C to 95°C, 1°C per 2s) was generated to confirm the absence of non-specific amplification. Standard curves were constructed for every primer pair using different cDNA dilutions to calculate the primer specific amplification efficiency. These efficiencies were incorporated in calculations of the expression values. The housekeeping genes with the highest efficiency and constant Ct values were selected and the obtained Ct values of these genes were used for normalization. Analysis of qPCR results was performed according to Pfaffl , producing relative expression values of the target gene. Significant differences in gene expression of the target gene were tested with pairwise fixed reallocation randomization .
T. urticae antifreeze proteins (AFPs) analysis
Gene expression analysis revealed the upregulation of 14 genes belonging to the same hypothetical protein family. A tBLASTn analysis (E-value cutoff = 1e-5) using these hypothetical proteins as queries was conducted against the T. urticae genome (http://bioinformatics.psb.ugent.be/orcae/overview/Tetur). Gene models were refined or created on the basis of homology, RNA-seq and/or EST support and RT-PCR. A BLASTp search, using these hypothetical proteins as queries, was performed against the NCBI nr protein sequence database to detect homologues in other organisms. Identity and similarity matrices were calculated using MatGAT 2.0  while an alignment was created using MUSCLE . T. urticae AFP structures were predicted using the Phyre2 server , in order to create a 3D model for the query sequence, and further edited in Swiss-Pdb viewer . The top match produced by Phyre2 is the protein sequence that showed the highest raw alignment score with the query sequence and is based on the number of aligned residues and the quality of alignment. This first match was selected to fit our protein to a 3D model. SignalP 4.0 was used for prediction of signal peptides of AFP protein sequences using default cutoff values for eukaryotes .
The LS-VL strain was used for the additional qPCR experiments assessing the effect of temperature on expression levels of T. urticae AFP genes. For the samples of 24°C, adult female mites were collected from the culture maintained at standard conditions. The group which experienced a cold shock consisted of female mites that were maintained on a leaf disk at 5°C during 7 days. Subsequently, RNA extraction and cDNA preparation were carried out as abovementioned.
Results and discussion
Micro-array analysis of diapausing T. urticae females
A laboratory strain (LS-VL) of T. urticae was grown under environmental conditions that induce diapause (light:dark regime of 8:16h at a temperature of 17°C). Under these conditions, approximately 10–30% of the population entered diapause. This made it possible to select and collect both diapausing and non-diapausing mites with a common genetic background that developed under identical environmental conditions (Figure 1A and C). The gene expression differences between these two distinctive groups should be mainly linked to diapause and were the foundation of this analysis.
We compared gene expression levels between diapausing and non-diapausing mites by analyzing genome-wide microarray experiments, including probes for a total number of 18,217 predicted unigenes . Analysis revealed that 1,994 genes (11% of all predicted genes) were differentially expressed (absolute FC ≥ 2, FDR <0.05) (Figure 2), of which 1078 were downregulated and 916 were upregulated in diapausing females. The expression level of the 20 most upregulated genes varied from a FC of 45 to 407 and included many genes with an unknown function (Additional file 2). Genes with known homologues included a T. urticae homologue of 5’ nucleotidase (tetur16g01680, FC 67) and a gene encoding a protein with a thioredoxin domain (tetur46g00020, FC 46). Among the 20 most downregulated genes in diapausing mites, genes with known functions and homologues included six vitellogenin genes (tetur39g00700, tetur39g00720, tetur39g00810, tetur43g00010, tetur43g00020, tetur516g00020, FC between 172 and 511), one lipase A gene (tetur09g06700 with a FC of 593), a gene belonging to the peroxidase family (tetur13g03760, FC 124), a midA gene (tetur09g06670, FC 116), a phosphatidylserine decarboxylase (tetur12g01500, FC 102), a glucose dehydrogenase (tetur03g09330, FC 99) and a serine protease (tetur19g00740, FC 95) (Additional file 2). Noteworthy, the downregulation of vitellogenin genes in diapausing T. uticae females was considered by Kawakami et al. (2011)  as a validation of true diapause.
A confirmation of the microarray experiment was conducted by qPCR. A similar trend was found between fold changes extracted from the complete genome microarray and the mean fold change data obtained by qPCR (Additional file 3).
Functional analysis and Gene Set Enrichment Analysis (GSEA) of differentially expressed genes in diapausing T. urticae females
A functional analysis of the differentially expressed genes revealed that 947 genes, 48% of total differentially expressed genes, could be assigned to a GO category. When the complete genome was categorized in a similar way, 6,105 genes, 33% of total number of T. urticae genes, were classified in GO terms. Details of the GO classification can be found in Additional file 4. Based on this GO classification, a Gene Set Enrichment Analysis (GSEA) was carried out to screen for an enrichment of the differentially expressed genes compared to the complete genome of T. urticae (Fisher Exact test FDR corrected p-values < 0.01, reduced to most specific terms) (Table 1). Notable was the enrichment of gene sets related to oxidoreductase activity (FDR = 3.46e-05), negative regulation of multicellular organism growth (FDR = 5.77e-05), carboxylic ester hydrolase activity (FDR = 5.19e-04) and fatty acid synthase activity (FDR = 2e-03). These GO terms are indicative of several expected adaptations on the molecular level further discussed below. Other diapause enriched gene sets were the categories of response to acid (FDR = 1e-03) and response to wounding (FDR = 2e-03) which are both related to stimuli responses. Underrepresented gene sets of the biological processes of translation (FDR = 6.3e-04), cell cycle (FDR = 2e-03) and RNA processing (FDR = 4e-03), indicate the arrest of cell growth and a decrease in protein synthesis. Subsequently, a GSEA was performed to retrieve information on the over- or underrepresentation of the upregulated genes compared to the downregulated genes (Figure 3). Clear differences were observed for a number of categories. For signal transduction, upregulated genes (n = 60) were more represented than downregulated genes (n = 13) (FDR = 1.86e-11). The cytoskeleton and actin binding related GO terms were also more abundant in diapausing mites (FDR = 2.13e-04 and FDR = 2e-03, respectively). On the other hand, the GO category of peptidase activity was overrepresented in the downregulated genes (FDR = 4.13e-08) as well as the GO terms of catabolic processes, carbohydrate binding, carbohydrate metabolic process and the cellular component of lysosome (FDR = 4.97e-04, 7e-03, 1e-04 and 3.2e-04 respectively) indicating a shutdown of anabolic and catabolic mechanisms in diapausing females. Combining this analysis with available literature and the draft genome sequence of T. urticae, we focused our following analysis on specific processes related to diapause in more detail below.
Adaptations in feeding and detoxification
A noticeable behavioral change indicating that T. urticae is preparing for a facultative reproductive diapause, is the fact that it stops feeding and empties the stomach, before migrating to hibernation shelters . This phenomenon was clearly reflected in the GO analysis. Many genes coding for enzymes involved in feeding, like digestive and detoxifying enzymes, were differentially expressed during diapause. About 50% of the total number of T. urticae cysteine peptidases (44 out of 87), showed a significant downregulation. Cysteine peptidases are involved in hydrolysis of peptide bonds using a catalytic cysteine and at present, 72 families of cysteine peptidases are identified . In diapausing mites, the downregulated genes belong to the papain-like (C1A) family (32 genes), the legumain (C13) family (2 genes) and the calpain (C2) family (10 genes)  (Additional file 5). Particularly the C1A family is known to be involved in the digestion of food [52, 53] but this family is also implicated in embryogenesis and metamorphosis [54, 55]. In the insects Praon volucre and D. melanogaster, the overexpression of these cysteine peptidases is associated with the termination of diapause to allow to resume digestion .
The genome of T. urticae harbors a remarkable number of genes encoding for ABC transporters and well-known detoxifying enzymes such as cytochrome P450 monooxygenases (CYPs), glutathione-S-transferases (GSTs), carboxyl/choline esterases (CCEs). This high number of detoxifying enzymes and ABC transporters in T. urticae is thought to be the result of lineage-specific expansions and has recently been shown to be involved in feeding, pesticide resistance and host plant adaptation [12, 31, 38]. The effects of diapause were also reflected in the regulation of these specialized genes. Of all 82 CYPs annotated in the genome, 28 were differentially expressed. Interestingly, the majority of differentially expressed CYP genes are upregulated in diapause (18 CYP genes belonging to the CYP2, CYP3 and CYP4 clan, fold change between 2.0-4.1) (Additional file 6). Next to detoxification, CYPs contribute to several vital processes including biosynthesis and degradation of ecdysteroids and juvenile hormones and the metabolism of plant allelochemicals and pesticides . In the silkmoth, Antheraea yamamai, a CYP gene (CYP4G15) was highly expressed during diapause in the pharate first instar larvae and was suggested to be involved in the diapause termination pathway . Another dormancy involved CYP is the daf-9 gene of the nematode C. elegans, belonging to the CYP2 clan and resembles steroidogenic and fatty acid hydroxylases as well as xenobiotic detoxifying genes. This gene serves as a central point of developmental control, produces hormonal signals and influences the decision of the nematode to choose between a developmental arrest as a third instar dauer larva or the completion to a reproductive adult . In addition, 10 CYPs were downregulated in T. urticae (FC between 2.0-23.2), of which tetur05g04000 was the most downregulated T. urticae CYP gene. The majority of these downregulated CYPs belong to the 389 subfamily that is one of the spider-mite-specific expansions and the regulation of this group appears to be greatly affected by the host plant range , suggesting their involvement in xenobiotic detoxification.
Next to CYPs, the expression of T. urticae carboxyl/choline esterases (CCEs) was also influenced by the diapausing status of T. urticae: 13 genes were downregulated and 18 genes showed an upregulation (Additional file 7). All differentially expressed CCEs in diapausing spider mites, except tetur19g00850, belonged to two new clades under the neurodevelopmental root of CCE phylogeny . Tetur19g00850 encodes the enzyme acetylcholinesterase, one of the best characterized CCEs, which is responsible for the hydrolysis of acetylcholine in the nervous system. This gene was two-fold upregulated in diapausing mites, and expression also appears to fluctuate in the moth Celerio euphorbiae during diapausing pupal development during winter . Also in other organisms, cholinesterases have been reported to be involved in diapause. In the brain of giant silkmoths, they have been associated with a reduction of electrical activity during pupal diapause . However, other studies observed changes in the cholinergic system that are more influenced by the growth and development of the nervous system than by the mechanism of diapause induction and termination [63, 64]. The downregulation of 13 CCE genes probably also reflects the decrease in food uptake during diapause.
Twenty-five percent of the 32 T. urticae GSTs  were downregulated (FC 2–5) while only one (tetur26g01490) was upregulated under diapausing conditions (Additional file 8). Three differentially expressed genes were classified into the Delta-class GSTs and 5 belonged to the Mu-class which was previously believed to be vertebrate-specific . One of the downregulated genes, tetur01g02320 (FC 2), belonged to the Omega-class. Recently, it was suggested that an Apis cerana cerana gene belonging to this class displayed protective effects against oxidative stress . In general, GSTs are important in Phase II of the detoxification process, acting by conjugating glutathione to xenobiotics or their derivatives. The lower expression of Delta- and Mu-class GSTs, can be associated with the lack of feeding and the reduced intake of potential xenobiotics. In addition, glutathione is often associated with the removal of reactive oxygen species (ROS) and results in a change of the ratio of reduced glutathione (GSH) over oxidized glutathione (GSSG). This ratio serves as an indicator of oxidative stress and is also influenced by GSTs, which are capable of catalyzing the oxidation of GSH via conjugation. It has been shown for the European corn borer Ostrinia nubilalis that GST activity in mitochondria was lower under diapause than under non-diapause conditions indicating a lower GSH/GSSG ratio creating oxidizing conditions [66, 67]. Similar to insects, the downregulation of GST genes in the two-spotted spider mite could indicate the transition of a reducing to an oxidative gut environment under diapause conditions.
The expression of ABC-transporters was also affected and of the 103 ABC genes annotated in the T. urticae genome , 20 genes displayed a differential expression (Additional file 9). The differentially expressed genes coding for ABC-transporters of class C were all downregulated (FC 2.1-7.3) except for one gene (FC 2). This class of ATP-binding cassette transporters constitutes efflux pumps, also named multidrug resistance associated proteins (MRP), and many proteins of this class have been implicated in reducing the cellular concentration of toxic compounds [68, 69]. The fact that diapausing spider mites do not feed, implies that also the intake of toxic compounds will decrease and downregulation of the ABCC genes could be a consequence of this. On the other hand, in the nematode C. elegans, an MRP-1, was detected to mediate the regulation of dauer larvae formation , and the precise role of ABCC genes in diapause needs to be further investigated. Four ABC-transporters class G were also differentially expressed. One gene, tetur17g02510 was upregulated (FC 3.5) and previously identified as a clear orthologue of D. melanogaster CG3327, also known as Early gene 23 (E23) . This E23 is a 20-OH ecdysone-induced ABC transporter that is capable of regulating 20E responses during metamorphosis, probably by removing 20E from cells . In Drosophila flies, it is also believed that E23 controls the circadian clock in adult flies through ecdysone-mediated expression of the clock gene vrille. Although T. urticae uses a different molting hormone than flies do (ponasterone A instead of 20E) , the upregulation of tetur17g02510 during diapause and its function at this life stage could be crucial and merits further investigation. In the ABC-transporters class H, six genes showed a differential expression ranging from a downregulation of 5.6 to an upregulation of 4.5. This transporters class was not associated with feeding, but recent studies suggest an involvement in cold hardening in D. melanogaster, pesticide resistance in Plutella xylostella pesticide and lipid transport to the cuticle in Tribolium castaneum.
Finally, genes encoding for intradiol ring-cleavage dioxygenases (ID-RCDs) were greatly influenced by the condition of diapause. Out of 17 ID-RCD genes, 11 were differentially expressed and all except one, showed a significant downregulation ranging from 2 to 37 times (Additional file 10). These ID-RCDs are among the most compelling cases of horizontal gene transfer in the genome of T. urticae, as they are specific for bacteria and fungi, with no known homologues in Metazoa [12, 31]. They seem to have an important role in host plant adaptation, and are also constitutely overexpressed in several highly multi resistant spider mite strains . In bacteria and fungi, they cleave aromatic rings of catecholic substrates, but their precise function in T. urticae remains until now unclear . Among detoxifying families, ID-RCDs showed one of the most prominent responses to diapause strongly suggesting an important role in digestion or detoxification.
Cryoprotection mechanisms: a potential new type of antifreeze proteins
Evolutionary adaptations among arthropods to enter diapause are not limited to those that downregulate the metabolism and energy production, but also include metabolic changes that allow to survive the direct effect of harsh conditions such as life-threatening temperatures. Already in 1957, it was reported that diapausing eggs of the silkworm are accumulating sorbitol and glycerol that serve as cryoprotectants . Thereupon, more polyols and sugars were discovered in various arthropods including mannitol, glucose, ribitol, arabinitol and threitol as well as trehalose and the amino acid alanine [6, 77, 78]. Synthesis of polyols is not exclusively assigned to diapause but is also often attributed to rapid cold-hardening [6, 79]. A biochemical study of diapause in T. urticae previously revealed an enrichment of potential cryoprotectants . It was found that during diapause, the level of glucose, maltose, inositol and ribitol increased significantly compared to non-diapausing mites and lowering the temperature to 5°C caused a shift in accumulation to mannitol and sorbitol. Worth mentioning is that the concentration of alanine, together with glutamate, did not decrease during diapause . Alanine is also found to serve as a cryoprotectant in Gynaephora groenlandicae[80, 81] and associated with freeze survival in Enchytraeus albidus and a freeze tolerant frog Rana sylvatica.
The results of our gene expression analysis partly support these metabolic changes at the molecular level. Inositol or 1D-myo-inositol are intermediates of the inositol phosphate metabolism and some direct key enzymes were differentially expressed (Figure 4, Additional file 11). Inositol monophosphatase (IMPA) dephosphorylates several molecules into 1D-myo-inositol and T. urticae orthologues of this enzyme, tetur32g00440 and tetur32g01900, were upregulated (FC 2.0 and 2.1 respectively) in diapausing spider mites. Another important enzyme is inositol oxygenase (MIOX) which converts 1D-myo-inositol into D-glucuronate and plays a key role in the ascorbate metabolism. Its orthologue, tetur19g00780, was 4.3 fold downregulated in diapausing spider mites, possibly to maintain high levels of these cryoprotectants. The enzyme, phosphatidylinositol phospholipase C, PLC, converts 1-phosphatidyl-1D-myo-inositol-4,5P2 to 1D-myo-inositol-1,4,5-P3 and its orthologue tetur05g05350 was upregulated (FC 2.2) during diapause. Subsequently, 1D-myo-inositol-1,4,5-P3 is dephosphorylated by inositol-1,4,5-trisphosphate 5-phosphatase, INPP5, to 1D-myo-inositol-1,4-P2. Tetur08g02900, one of the genes encoding for this INPP5 enzyme, was upregulated 2.5 times. In conclusion, the differential expression of these enzymes belonging to the inositol phosphate pathway endorses the hypothesis that T. urticae synthesizes inositol which is involved in cryoprotection. It was first discovered that inositol accumulates in hibernating adults of beetles [84, 85] and afterwards associated with cold hardening in several insects [85–87]. Also the hibernating adult house spider, which belongs together with T. urticae to the Arachnida, was shown to protect against cold-shock damage by the synthesis of inositol . In hydrophilic soil arthropods, the accumulation of myo-inositol is essential to increase the osmotic pressure and to maintain their body fluids in a hyperosmotic state. This causes net water uptake from the atmosphere by passive diffusion . Next to inositol related pathways, our analysis also revealed that a maltase gene (tetur06g01620) is 4.4 times downregulated, and this could indicate an increase of maltose, which was previously observed by direct measurements in T. urticae. Maltose is a disaccharide formed from two units of glucose and its hydrolysis is catalyzed by the enzyme maltase which regenerates glucose levels. This maltose sugar is also enriched in several species during diapause and after cold hardening [90, 91] and is already used as a supplementary cryoprotectant for algae, protozoa  and in pretreatment of embryogenic cultures of the maritime pine tree .
One of the most striking differentially expressed genes in our analysis consisted of a small family of ‘hypothetical proteins’ (Additional file 2). This family was manually annotated in the genome database of T. urticae, and consisted of a set of 20 genes and 2 pseudogenes (tetur22g03073 and tetur63g00100). Of those, only 16 genes were represented by probes on the microarray and 14 of these were upregulated in diapausing spider mites, with fold changes varying between 4 and 164 (Additional file 12). Gene loci are distributed over three genomic scaffolds: 22 (15), 63 (4) and 283 (1) and all genes except tetur22g03033 are intronless. Genes on scaffold 22 are tightly clustered within a 90 kb region, but the proteins showed only moderate sequence similarity, suggesting their proliferation in the T. urticae genome is the result of “ancient” tandem duplication events. Tetur283g00030, tetur63g00030, tetur63g00050, tetur63g00070 and tetur63g00090 showed very high identity values (96.7% - 100%) with tetur22g02550, tetur22g02640, tetur22g02670, tetur22g02730 and tetur22g02690, respectively (Additional file 13). This might suggest a recent duplication event, although the possibility that the genes on scaffold 22 and 63/283 actually represent the same gene due to assembly issues caused by allelic variants cannot be excluded at this stage.
All members of this hypothetical protein family, except tetur22g02690, tetur22g03063 and tetur63g00090, have a predicted signal peptide  and hence are probably secreted by cells. Strikingly, an InterPro-scan revealed that the majority (16/20) of these hypothetical proteins contained the “Insect Antifreeze Protein” motif (IPR016133 with an E-value between 1e-04 and 1e-06). A BLASTp search against the non-redundant protein database of NCBI also hit to antifreeze proteins of Coleoptera with a low to moderate E-value (between 1e-02 and 1e-06). Furthermore, all members from this family were predicted as an insect AFP by AFP-Pred, a recently developed software tool using a “random forest” approach for the prediction of antifreeze proteins .
Accordingly, we aligned and compared the T. urticae mite sequences with well-studied insect AFP sequences of two beetles: Dendroides canadiensis (AAF86362 and AAB94303) and Tenebrio molitor (1L1I_A) (Additional file 14). These insect AFPs consist of 78 to 148 amino acids and their cysteine (Cys) content ranges from 15 to 19%, whereas the spider mite sequences measure between 92 and 210 amino acids and consist for 19 to 25% of Cys. In addition, typical insect AFPs are characterized by 7 repeats of 12-or 13-mer repeats (Thr-Cys-Thr-X-Ser-X-X-Cys-X-X-Ala-X) with at least every sixth residue a cysteine [95, 96]. The Cys residues can form stable disulfide bridges throughout the protein and are flanked by Thr residues, which are thought to be responsible for the ice binding sites . In T. urticae, the predicted AFPs comprise a continuous repeated sequence (Asn-Cys-Thr-X-Cys-X-X-Cys-X-Asn-Cys-X). It is clear that two of the Cys (at the sixth residue, involved in disulfide bridging) are conserved, but two other Cys are found at positions where in beetles a conserved Ser and Ala residue is found. Furthermore, the conserved Cys are not only flanked by Thr, but a variety of Asp, Asn and Thr. These hydrophilic residues are capable of ice-binding and inhibition of ice crystal growth [98, 99]. The specific arrangement and position of these residues is however fundamental for the activity of AFPs in order to bind to the ice surface with a particular orientation. Hence, all 20 protein sequences were submitted to Phyre2 web service for structure prediction. Six out of 20 query sequences (tetur22g02640, tetur22g02690, 22g03033, tetur22g03063, tetur63g00030 and tetur63g00090) had a top scoring match with an insect AFP template structure (pdb: d1ezga) in the Phyre2 library. Subsequently, 2D and 3D models for these six T. urticae proteins were returned based on this match. For example, the 3D model of tetur22g03033 (53% of the tetur22g02640 protein sequence has been modeled with 97.1% confidence, Additional file 15) showed the typical configuration of cysteine rich insect AFP consisting of six β-strands containing the typical Cys-residues (Figure 5) (at six AA distance) producing disulfide bridges. The important ice-binding amino acids (Asn, Asp and Thr) are located on the outside of the protein. Interestingly, the extra Cys-residues (at three AA distance) specific for mite AFP sequences, are pointing inwards and may form alternative cysteine bridges in vivo. The typical ice-binding surface [100–102] was not apparent in the 3D homology model of tetur22g03033, but these outfacing hydrophilic residues might still be involved in ice binding. Recently, divergent tertiary structures were also reported for Rhagium inquisitor and snow fleas . In conclusion, sequence alignments, and in silico predictions including a homology model strongly suggest that this small family of proteins might represent a new class of arthropod AFPs. Nevertheless, functional expression combined with activity tests using nanoliter osmometry should provide formal evidence of the activity and ice binding properties of these proteins [105–107].
The production of antifreeze proteins is a protective mechanism against the effects of lowered environmental temperatures and freezing . Typically, cysteine-rich insect antifreeze proteins adsorb to seeded ice crystals and inhibit enlargment of these crystals  which causes a difference between the freezing temperature and the melting point, a process that is known as thermal hysteresis. Thermal hysteresis antifreeze proteins have also been found in vertebrates, invertebrates, fungi, bacteria and plants . At least fifty species of insects and many terrestrial arthropods [107, 109] are known to produce AFPs and although many have low levels of thermal hysteresis, some insects produce hyperactive antifreeze proteins. The common yellow mealworm beetle, T. molitor, contains antifreeze proteins that can account a thermal hysteresis of 5.5°C at a concentration of 1 mg/ml . Comparative results are found in Hypogastrura harveyi[110, 111] and C. fumiferana who display similar thermal hysteresis properties . The magnitude of thermal hysteresis effect occurring in a species depends not only on the structural type of proteins, but can differ greatly depending on the moment of sampling or the level of acclimatization of insects .
We performed additional qPCR experiments to investigate the association of AFPs with diapause and/or cold stress of 5 typical AFP genes (tetur22g02690, tetur22g02640, tetur22g02730, tetur22g02790 and tetur22g02670). Expression levels were compared between non-diapausing mites at 24°C (standard rearing conditions), and mites that were submitted to a cold shock of 7 days, together with diapausing and non-diapausing mites at 17°C (the conditions of the microarray experiments). Results showed that the expression of these genes was overall more affected by the diapause condition, than the reduction in temperature in T. urticae, indicating that for some genes the expression is probably regulated by physiological, rather than environmental changes (Figure 6). All genes were significantly high upregulated in diapausing mites at 17°C when compared to the expression of these genes in active mites at 17°C, confirming microarray results. Only one gene, tetur22g02690, was equally high expressed both in diapause mites at 17°C and in mites treated by cold shock. Of the 4 other genes tested, two did respond to cold stress (tetur22g02670, tetur22g02730) but with small changes in expression, while tetur22g02790 and tetur22g02640 were only upregulated to high levels in diapausing mites.
It was previously shown that the production of antifreeze proteins is not exclusively associated with diapause. In T. molitor, D. canadensis, R. inquisitor and H. harveyi, antifreeze proteins could be collected from insects that were exposed to cold. In the spruce budworm, transcripts of antifreeze proteins were found in both the first larval instar and in the second diapausing instar . This indicates that both cold hardening and diapause can trigger the production of AFPs in insects, similar to what was shown for T. urticae in this study. It was demonstrated previously in T. urticae that diapausing forms had a lower supercooling point than non-diapausing forms, even if those active forms were cold acclimated for 10 days at 5 or 0°C . Furthermore, the 5°C cold acclimated active mites showed a lower supercooling point than the not acclimated active mites, indicating cryoprotection caused by cold and diapause, but the involvement of mite AFPs identified in this study was not assessed and should be further studied. Of particular note, homologues of spider mite AFPs were not found in other mite and tick species for which a draft genome or transcriptome is available, such as Varroa destructor, Metaseiulus occidentalis, Ixodes scapularis and Panonychus citri.
The cytoskeleton plays a critical role in the process of cold-acclimatization in arthropods. For Culex pipiens, it was reported that two actin genes were highly overexpressed in early diapause [114, 115]. In addition, the upregulation of actin genes was also associated with a redistribution of polymerized actin during exposure to cold in non-diapausing and diapausing mosquitoes. The functional analysis of our enriched GO terms showed that more genes associated with the cellular component of the cytoskeleton and the molecular function of actin binding were upregulated than downregulated during diapause (Figure 3). Of the 15 genes associated with the cytoskeleton, five showed high identity with the myosin heavy chain protein (a BLASTp hit (E-value < 1e-144) of Camponotus floridanus). Myosin heavy chain (MHC) isoforms are complex, multifunctional contractile proteins with an ATPase activity and are associated with other MHCs and myosin light chains to form the sarcomeric thick filament [116, 117]. Many organisms have several MHC isoforms that fulfill various functions and differ in their expression patterns. They were found to be involved in translocation of actin filaments  and are selectively repressed during dissolution of insect skeletal muscles at the end of the metamorphosis of Manduca sexta. In crustaceans, it has been shown that the increased expression of MHC was temperature related and could not specifically be ascribed to diapause . More recently, a sequence that matched with a Drosophila MHC was found in the diapause-destined Artemia fransicana embryo . In T. urticae, two genes, tetur20g00190 and tetur20g00160, that are coding for muscle-specific proteins 300 (MSP300) are also upregulated and are associated with the actin microfilament system . A similar muscular protein, MSP20, was also found to be upregulated in the diapausing P. volucre. The gene tetur23g00300 is related to a formin homology 2 domain containing protein (a BLASTp hit (E-value = 8.59e-133) with formin homology 2 domain containing 3 of Homo sapiens) and is known as an actin-organizing protein. Tetur23g00300 was upregulated and the GO terms of actin binding and cytoskeleton components were assigned to this gene. In the GO-term pool of actin binding, two genes, tetur11g04370 and tetur11g04320, with a LIM domain were upregulated. These LIM domains are present in many proteins with diverse functions. One of these functions involves a direct or indirect interaction with the actin cytoskeleton. Besides these genes that came up in the GSEA, other cytoskeleton-related genes also showed a differential expression. Four ankyrin related genes, tetur01g08790, tetur11g02270, tetur15g02730 and tetur19g00130, were upregulated, while a cuticular protein, tetur06g01680, was downregulated. This is exactly the opposite of what was observed in the proteomic profiling of the aphid parasitoid P. volucre. Ankyrin proteins are implicated in coupling integral membrane proteins to the cytoskeleton network which can cause dynamic interactions among cytoskeletal filaments. Moreover, the upregulation of ankyrin-like proteins is linked to the recovery of insects from cold stress . The differential gene expression of genes coding for cytoskeleton-related proteins suggests that a reorganization of these structures might take place during diapause and that the cytoskeleton is considerably distinct from that of non-diapausing spider mites.
Already in 1974, Veerman  detected that diapausing females of T. urticae showed a more than two-fold increase in total keto-carotenoid pigment content than the green summer females. Further research demonstrated that carotenoids not only cause the striking color difference between diapausing and non-diapausing females, but also that carotenoids are functionally involved in photoperiodism [124, 125]. The perception of this photoperiod by mites requires the presence of a photoreceptor and was long unknown. Rearing albino strains on carotenoid-containing and carotenoid-free diets demonstrated the potential of vitamin A as photoreceptor . When wild-type females fed on the carotenoid free diets, they could still enter diapause which according to Bosse et al. (1996) could be allocated to the transfer of maternal carotenoids to the egg . In this study it was also assumed that the amount of carotenoids present in the egg was not only enough for diapause induction but also for the body coloration of the diapausing female. Recently, Grbić et al. (2011) demonstrated that T. urticae possesses two horizontally transferred fused carotenoid cyclase/synthase genes and three carotenoid desaturase genes . These genes most likely originate from a Mucorales-related fungal donor and are similar to the pea aphid, Acyrthosiphon pisum, carotenoid genes . Varying numbers of diversified carotenoid desaturases have been found in 34 aphid species and the related insect group Adelgidae . Most recently, Cobbs et al. (2013)  reported the third case of fungal carotenoid biosynthesis gene homologues in gall midges (Diptera: Cecidomyiidae). Several other cases of horizontal gene transfer have been described in arthropods and most of them have their origin in bacteria, particularly in the Wolbachia genus. Only few of those transferred genes to arthropods have been proven to be functional [130, 131].
In our microarray experiment, four genes involved in carotenoid synthesis were clearly differentially expressed in diapausing females (Additional file 16). Tetur01g11260 and tetur11g04840 encoding for the fused lycopene cyclase/ phytoene synthase genes showed an upregulation (FC 6.2) and a downregulation (FC 1.8) respectively. Two out of three genes of the phytoene desaturase genes showed the same trend of regulation. One gene,tetur01g11270, displayed an upregulation with fold change of 5.8 and another phytoene desaturase gene, tetur11g04820, was downregulated (FC 24.5). Subsequently, we analyzed the expression of this complete set of carotenoid genes by qPCR in two green morphs (LS-VL and London), two red morphs (MR-VL and Tu-SB9) and one diapausing condition of one of the green morphs (LS-VL, 17°C and 8:16 L:D) (Figure 1). All results were compared to one of the green morphs (LS-VL) cultured at 24°C and 16:8 L:D (Figure 7). The expression of the fused lycopene cyclase/ phytoene synthase genes varied depending on the color and physiological state (diapause or not). The cyclase-synthase gene tetur01g11260 showed a consistent higher expression in diapausing and red morphs, while the cyclase-synthase gene tetur11g04840 was downregulated in all conditions compared to non-diapausing LS-VL. The expression of the three phytoene desaturase genes resembles the trend of the cyclase/synthase genes, as one gene,tetur01g11270, is upregulated in both red morphs and diapausing mites and the other two genes, tetur11g04810 and tetur11g04820, are downregulated in red and diapausing mites. Moreover, the downregulation of tetur11g04810 was more pronounced: fold changes for the MR-VL red morph and diapausing mites were 14.3 and 6.3 respectively. In the green London morph, this gene was upregulated with a fold change of 4.7 when compared to the green LS-VL morph. From this experiment, it seems that two out of five carotenoid genes are actively involved in the carotenoid synthesis during diapause conditions and in the active red morphs. This data strengthens the general expectation that these genes will enable mites, aphids and gall midges to synthesize their own carotenes, however the functionality of these genes has yet to be determined. Besides this group of horizontally transferred carotenoid genes, genes encoding for GGPP synthase and carotenoid oxygenases were also detected in T. urticae. One of these carotenoid oxygenase genes, tetur06g06440 showed a downregulation in diapausing spider mites which could indicate a shift in beta-carotene molecules.
Previously, the gene regulation of one carotenoid desaturase (tetur01g11270) and one fused carotenoid cyclase-synthase (tetur01g11260), two of the whole set of five genes, were examined in green morphs, red morphs and diapausing mites by Altincicek et al. (2011) . Their expression was subsequently investigated and showed that both genes were upregulated in red morphs in comparison to the green morphs. Diapause conditions had a great influence on the expression, however at the final point of sampling, the green diapause morphs did not display the typical orange color indicative of diapause. Hence, possibly mites in this study had not truly entered diapause, as many populations do not enter it, and results reflect that the conditions to induce diapause are sufficient to trigger the differential expression. Alternatively, pigmentation in the adult morphs might only slowly develop in these populations after the last molt, which has been previously documented .
It was suggested that carotenoids and their derivates fulfill many important roles in insect physiology and life history . According to Bosse and Veerman, [126, 134], the photoperiodic induction of diapause requires carotenoids in moths, butterflies, wasps and spider mites. Moreover, the receptor responsible for the photoperiod measurement appears to be an opsin photoreceptor that needs vitamin A, requiring carotenoids for its synthesis . Results suggest that the horizontal transfer of carotenoid synthesis genes may allow spider mites to enter diapause without depending on dietary intake of carotenoids. Another function of carotenoids is to protect arthropods from ultraviolet radiation and oxidative stress . In T. urticae, it was confirmed that diapausing females exhibit a lower mortality at different doses of UV than the green summer females . Both these threats are known to seriously damage the cellular components of arthropods. In insects, diapause is often accompanied by oxidative stress [32, 136–138] and besides the protective function of GSTs, an upregulation of carotenoids during diapause in T. urticae could increase the chance to survive these harsh conditions.
Signal transduction and G-protein coupled receptors
The token stimuli for the two-spotted spider mite to initiate diapause are a decrease in temperature or in the length of the photoperiod. These signals are most likely transmitted through messengers like hormones or other chemical substances in order to initiate the structural and functional changes related to diapause. In this experiment, four genes that were predicted to operate as G-protein coupled receptors (GPCRs) displayed a differential expression. Tetur14g00960, tetur05g00270 and tetur03g05860 show high similarity with the neuropeptide Y receptor, a tachykinin receptor and neuropeptide F, respectively, and were all upregulated (FC ranging from 2.2 to 7.6). Also one relaxin-like GPCR, tetur15g00700, was downregulated (FC 5.1) . Furthermore, the GO enrichment analysis revealed a high number of genes associated with signal transduction (GO:0007165), including a number of kinases/phosphatases and hormone responsive genes (Additional file 17). Interestingly, a number of nuclear receptors were downregulated, including the retinoid X receptor beta, tetur01g09220, which has been previously implicated in diapause . In addition, one of the 8 paralogs of Drosophila DHR96, tetur36g00260, was downregulated in diapausing mites. An increased expression of this paralogue has been previously associated with the response to xenobiotics in T. urticae. The downregulation of DHR96 in diapausing mites fits to the overall pattern of a decreased expression of genes associated with a reduced intake of xenobiotics. In Bombyx mori, a diapause hormone was the first chemical substance that was identified as a maternal control factor to block progeny development . Subsequently, the GPCR for the insect diapause hormone was also discovered and turned out to serve as an inducer for embryonic diapause . Lately, also other GPCRs belonging to a pyrokinin/pheromone biosynthesis activating neuropeptide family have been found in the European corn borer and were downregulated during larval diapause . The differential expression of the GPCRs found in T. urticae might indicate a role in the mechanisms underlying triggering and regulation of diapause.
Comparison of diapause responses
The gene regulation and gene expression patterns related to diapause appear to vary considerably among different organisms. Ragland et al. (2009) compared arrays of dormant and non-dormant phenotypes of S. crassipalpis, D. melanogaster and C. elegans and deduced that only a small set of 10 differentially expressed genes were regulated similarly . Even among dipterans, patterns of gene expression were not more similar to each other than to C. elegans. This finding could imply that the regulation of diapause is not very conserved. On the other hand, two common physiological responses, i.e. those related to metabolism and stress resistance, were affected similarly during diapause across these three species. In order to investigate to which extent regulation of diapause by T. urticae aligns with that of these three species, we compared the pattern of the 10 across-species markers genes with that of their homologs in spider mites (Table 2). For 5 of these genes, the expression patterns were similar across the four species. One of those genes, pyruvate carboxylase, is associated with both glycolysis as gluconeogenesis. The other irreversible gluconeogenesis enzyme phosphoenolpyruvate carboxykinase showed a slight downregulation where it displayed an upregulation in the dipterans and C. elegans. This variation in gene expression during dormancy i.e. dauer and diapause most likely reflects insects, nematodes and Acari having highly diverged life histories and acquired markedly different life styles during the course of their evolution. Furthermore, the life stage and profundity of diapause varies between these species and is likely reflected in diapause gene regulation.
In this study, a genome-wide microarray system was used to determine patterns of differential gene expression between diapausing and non-diapausing females of the two-spotted spider mite T. urticae. Spider mites are chelicerates, a lineage that diverged more than 450 Mya from other arthropod lineages. In these experiments, RNA was collected from a single population in which approximately 30% of the individuals had entered diapause under the applied inducing conditions. This allowed assessing gene expression differences between samples with a highly similar genetic background and which were reared under identical environmental conditions. GO analysis of differentially expressed genes revealed that many metabolic pathways are affected in diapausing females, especially those related to digestion and detoxification, cryoprotection, carotenoid synthesis and the organization of the cytoskeleton. Some of these transcriptional responses confirmed earlier metabolite studies in spider mites on the molecular level. However, we also discovered previously unknown adaptations, such as the proliferation in T. urticae of an unique class of putative antifreeze proteins. These proteins were among the genes most strongly upregulated in diapausing females. Although further validation of their activity is needed, these proteins might be suitable for biotechnological applications related to cryoprotection. Furthermore, we documented differential expression of genes involved in signaling and signal transduction such as G-coupled receptor genes, possibly being the starting point towards understanding the regulation of diapause in this species. Comparison with other genome-wide diapause expression studies suggests that some fundamental changes are conserved, but that overall specific strategies have evolved in different species.
ATP-binding cassette (transporter)
Cytochrome P450 monooxygenase
False Discovery Rate
Gene Set Enrichment Analysis
G-protein coupled receptors
intradiol ring-cleavage dioxygenases
Myosin Heavy Chain
Reactive Oxygen Species.
Denlinger D: Regulation of diapause. Annu Rev Entomol. 2002, 47: 93-122. 10.1146/annurev.ento.47.091201.145137.
Denlinger D: Induction and termination of pupal diapause in Sarcophaga (Diptera: Sarcophagidae). Biol Bull. 1972, 142: 11-24. 10.2307/1540242.
Veerman A: Diapause. Spider mites their Biol Nat enemies Control. Edited by: Sabelis MW, Helle W. 1985, Amsterdam: Elsevier, 279-317.
Rinehart JP, Li A, Yocum GD, Robich RM, Hayward S, Denlinger D: Up-regulation of heat shock proteins is essential for cold survival during insect diapause. Proc Natl Acad Sci U S A. 2007, 104: 11130-11137. 10.1073/pnas.0703538104.
Hahn DA, Denlinger DL: Meeting the energetic demands of insect diapause: nutrient storage and utilization. J Insect Physiol. 2007, 53: 760-773. 10.1016/j.jinsphys.2007.03.018.
Kostal V, Tamura M, Tollarova M, Zahradnickova H: Enzymatic capacity for accumulation of polyol cryoprotectants changes during diapause development in the adult red firebug, Pyrrhocoris apterus. Physiol Entomol. 2004, 29: 344-355. 10.1111/j.0307-6962.2004.00396.x.
MacRae TH: Gene expression, metabolic regulation and stress tolerance during diapause. Cell Mol Life Sci. 2010, 67: 2405-2424. 10.1007/s00018-010-0311-0.
Clark MS, Worland MR: How insects survive the cold: molecular mechanisms-a review. J Comp Physiol B. 2008, 178: 917-933. 10.1007/s00360-008-0286-4.
Doucet D, Walker VK, Qin W: The bugs that came in from the cold: molecular adaptations to low temperatures in insects. Cell Mol Life Sci. 2009, 66: 1404-1418. 10.1007/s00018-009-8320-6.
Veerman A: Physiological aspects of diapause in plant-inhabiting mites. Acari Reprod Dev life-history Strateg. Edited by: Schuster R, Murphy PW. 1991, NY: Chapman and Hall, 245-265. 1
Van Leeuwen T, Vontas J, Tsagkarakou A, Dermauw W, Tirry L: Acaricide resistance mechanisms in the two-spotted spider mite Tetranychus urticae and other important Acari: a review. Insect Biochem Mol Biol. 2010, 40: 563-572. 10.1016/j.ibmb.2010.05.008.
Dermauw W, Wybouw N, Rombauts S, Menten B, Vontas J, Grbic M, Clark RM, Feyereisen R, Van Leeuwen T: A link between host plant adaptation and pesticide resistance in the polyphagous spider mite Tetranychus urticae. Proc Natl Acad Sci U S A. 2013, 110: E113-E122. 10.1073/pnas.1213214110.
Khodayari S, Moharramipour S, Larvor V, Hidalgo K, Renault D: Deciphering the metabolic changes associated with diapause syndrome and cold acclimation in the two-spotted spider mite Tetranychus urticae. PLoS One. 2013, 8: e54025-10.1371/journal.pone.0054025.
Kawakami Y, Goto SG, Ito K, Numata H: Suppression of ovarian development and vitellogenin gene expression in the adult diapause of the two-spotted spider mite Tetranychus urticae. J Insect Physiol. 2009, 55: 70-77. 10.1016/j.jinsphys.2008.10.007.
Kroon A, Veenendaal RL, Bruin J, Egas M, Sabelis MW: Predation risk affects diapause induction in the spider mite Tetranychus urticae. Exp Appl Acarol. 2004, 34: 307-314.
Kroon A, Veenendaal RL, Egas M, Bruin J, Sabelis MW: Diapause incidence in the two-spotted spider mite increases due to predator presence, not due to selective predation. Exp Appl Acarol. 2005, 35: 73-81. 10.1007/s10493-004-1980-x.
Gotoh T, Tokioka T: Genetic compatibility among diapausing red, non-diapausing red and diapausing green forms of the two-spotted spider mite, Tetranychus urticae Koch (Acari: Tetranychidae). Jpn J Ent. 1996, 64: 215-225.
Helle W: Genetic variability of photoperiodic response to an arrhentokous mite (Tetranychus urticae). Entomol Exp Appl. 1968, 11: 101-113. 10.1111/j.1570-7458.1968.tb00078.x.
Veerman A: Aspects of the induction of diapause in a laboratory strain of the mite Tetranychus urticae. J Insect Physiol. 1977, 33: 703-711.
Goka K, Takafuji A: Genetical studies on the diapause of the two-spotted spider mite Tetranychus urticae (Koch) (2). Appl Entomol Zool. 1991, 26: 77-84.
Ignatowicz S, Helle W: Genetics of diapause supression in the two-spotted spider mite, Tetranychus urticae Koch. Exp Appl Acarol. 1986, 2: 161-172. 10.1007/BF01213759.
Kawakami Y, Numata H, Ito K, Goto SG: Dominant and recessive inheritance patterns of diapause in the two-spotted spider mite Tetranychus urticae. J Hered. 2010, 101: 20-25. 10.1093/jhered/esp085.
Pritchard AE, Baker EW: A guide to the spider mites of deciduous fruit trees. Hilgardia. 1952, 21: 253-287.
Boudreaux HB: The effect of relative humidity on egg-laying, hatching, and survival in various spider mites. J Insect Physiol. 1958, 2: 65-72. 10.1016/0022-1910(58)90029-5.
Boudreaux HB: Biological aspects of som phytophagous mites. Annu Rev Entomol. 1963, 8: 137-154. 10.1146/annurev.en.08.010163.001033.
McEnroe WD: The Control of Water Loss by the Two-Spotted Spider Mite (Tetranychus telarius). Ann Entomol Soc Am. 1961, 54: 883-887.
Khodayari S, Moharramipour S, Kamali K, Jalali Javaran M, Renault D: Effects of acclimation and diapause on the thermal tolerance of the two-spotted spider mite Tetranychus urticae. J Therm Biol. 2012, 37: 419-423. 10.1016/j.jtherbio.2012.04.005.
Foott WH: eotactic response of the two-spotted spider mite, Tetranychus urticae Koch (Acarina, Tetranychidae). Proc Entomol Soc OntarioOntario, Volume95. 1965, 106-108.
Parr WJ, Hussey NW: Diapause in the glasshouse red spider mite (Tetranychus urticae Koch): a synthesis of present knowledge. Hortic Res. 1966, 6: 1-21.
Suzuki T, Kojima T, Takeda M, Sakuma M: Photo-orientation regulates seasonal habitat selection in the two-spotted spider mite, Tetranychus urticae. J Exp Biol. 2013, 216 (Pt 6): 977-983.
Grbić M, Van Leeuwen T, Clark RM, Rombauts S, Rouzé P, Grbić V, Osborne EJ, Dermauw W, Ngoc PCT, Ortego F, Hernández-Crespo P, Diaz I, Martinez M, Navajas M, Sucena É, Magalhães S, Nagy L, Pace RM, Djuranović S, Smagghe G, Iga M, Christiaens O, Veenstra J a, Ewer J, Villalobos RM, Hutter JL, Hudson SD, Velez M, Yi S, Zeng J: The genome of Tetranychus urticae reveals herbivorous pest adaptations. Nature. 2011, 479: 487-492. 10.1038/nature10640.
Ragland GJ, Denlinger DL, Hahn DA: Mechanisms of suspended animation are revealed by transcript profiling of diapause in the flesh fly. Proc Natl Acad Sci USA. 2010, 107: 14909-14914. 10.1073/pnas.1007075107.
Dunlop JA: Geological history and phylogeny of Chelicerata. Arthropod Struct Dev. 2010, 39: 124-142. 10.1016/j.asd.2010.01.003.
Van Leeuwen T, Stillatus V, Tirry L: Genetic analysis and cross-resistance spectrum of a laboratory-selected chlorfenapyr resistant strain of two-spotted spider mite (Acari: Tetranychidae). Exp Appl Acarol. 2004, 32: 249-261.
Van Leeuwen T, Van Pottelberge S, Tirry L: Comparative acaricide susceptibility and detoxifying enzyme activities in field-collected resistant and susceptible strains of Tetranychus urticae. Pest Manag Sci. 2005, 61: 499-507. 10.1002/ps.1001.
Van Leeuwen T, Demaeght P, Osborne EJ, Dermauw W, Gohlke S, Nauen R, Grbic M, Tirry L, Merzendorfer H, Clark RM: Population bulk segregant mapping uncovers resistance mutations and the mode of action of a chitin synthesis inhibitor in arthropods. Proc Natl Acad Sci U S A. 2012, 109: 4407-4412. 10.1073/pnas.1200068109.
Kim S-H, Hamada T: Rapid and reliable method of extracting DNA and RNA from sweetpotato, Ipomoea batatas (L). Lam. Biotechnol Lett. 2005, 27: 1841-1845. 10.1007/s10529-005-3891-2.
Demaeght P, Dermauw W, Tsakireli D, Khajehali J, Nauen R, Tirry L, Vontas J, Lümmen P, Van Leeuwen T: Molecular analysis of resistance to acaricidal spirocyclic tetronic acids in Tetranychus urticae: CYP392E10 metabolizes spirodiclofen, but not its corresponding enol. Insect Biochem Mol Biol. 2013, 43: 544-554. 10.1016/j.ibmb.2013.03.007.
ORCAE: online resource for community annotation of eukaryotes.http://bioinformatics.psb.ugent.be/orcae/overview/Tetur,
Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.
Myhre S, Tveit H, Mollestad T, Laegreid A: Additional gene ontology structure for improved biological reasoning. Bioinformatics. 2006, 22: 2020-2027. 10.1093/bioinformatics/btl334.
Blüthgen N, Kiełbasa SM, Herzel H: Inferring combinatorial regulation of transcription in silico. Nucleic Acids Res. 2005, 33: 272-279. 10.1093/nar/gki167.
Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Bioinforma methods Protoc, Volume 132. Edited by: Misener S, Krawetz SA. 1999, Totowa NJ: Humana Press Inc., 365-386.
Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: 2002-2007.
Pfaffl MW, Horgan GW, Dempfle L: Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 2002, 30: e36-10.1093/nar/30.9.e36.
Campanella JJ, Bitincka L, Smalley J: MatGAT: an application that generates similarity/identity matrices using protein or DNA sequences. BMC Bioinformatics. 2003, 4: 29-10.1186/1471-2105-4-29.
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.
Kelley LA, Sternberg MJE: Protein structure prediction on the Web: a case study using the Phyre server. Nat Protoc. 2009, 4: 363-371. 10.1038/nprot.2009.2.
Guex N, Peitsch MC: SWISS-MODEL and the Swiss-Pdb Viewer: An environment for comparative protein modeling. Electrophoresis. 1997, 18: 2714-2723. 10.1002/elps.1150181505.
Petersen TN, Brunak S, von Heijne G, Nielsen H: SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011, 8: 785-786. 10.1038/nmeth.1701.
Santamaría ME, Hernández-Crespo P, Ortego F, Grbic V, Grbic M, Diaz I, Martinez M: Cysteine peptidases and their inhibitors in Tetranychus urticae: a comparative genomic approach. BMC Genomics. 2012, 13: 307-10.1186/1471-2164-13-307.
Cristofoletti PT, Ribeiro AF, Terra WR: The cathepsin L-like proteinases from the midgut of Tenebrio molitor larvae: sequence, properties, immunocytochemical localization and function. Insect Biochem Mol Biol. 2005, 35: 883-901. 10.1016/j.ibmb.2005.03.006.
Soares-Costa A, Dias AB, Dellamano M, de Paula FFP, Carmona AK, Terra WR, Henrique-Silva F: Digestive physiology and characterization of digestive cathepsin L-like proteinase from the sugarcane weevil Sphenophorus levis. J Insect Physiol. 2011, 57: 462-468. 10.1016/j.jinsphys.2011.01.006.
Cho W-L: Mosquito cathepsin B-like protease involved in embryonic degradation of vitellin is produced as a latent extraovarian precursor. J Biol Chem. 1999, 274: 13311-13321. 10.1074/jbc.274.19.13311.
Liu J, Shi G-P, Zhang W-Q, Zhang G-R, Xu W-H: Cathepsin L function in insect moulting: molecular cloning and functional analysis in cotton bollworm, Helicoverpa armigera. Insect Mol Biol. 2006, 15: 823-834. 10.1111/j.1365-2583.2006.00686.x.
Colinet H, Renault D, Charoy-Guével B, Com E: Metabolic and proteomic profiling of diapause in the aphid parasitoid Praon volucre. PLoS One. 2012, 7: e32606-10.1371/journal.pone.0032606.
Baker DA, Russell S: Gene expression during Drosophila melanogaster egg development before and after reproductive diapause. BMC Genomics. 2009, 10: 242-10.1186/1471-2164-10-242.
Feyereisen R: Insect P450 enzymes. Annu Rev Entomol. 1999, 44: 507-533. 10.1146/annurev.ento.44.1.507.
Yang P, Tanaka H, Kuwano E, Suzuki K: A novel cytochrome P450 gene (CYP4G25) of the silkmoth Antheraea yamamai: cloning and expression pattern in pharate first instar larvae in relation to diapause. J Insect Physiol. 2008, 54: 636-643. 10.1016/j.jinsphys.2008.01.001.
Gerisch B, Antebi A: Hormonal signals produced by DAF-9/cytochrome P450 regulate C. elegans dauer diapause in response to environmental cues. Development. 2004, 131: 1765-1776. 10.1242/dev.01068.
Grzelak K, Lassota Z, Wroniszewska A: Acetylcholinesterase in Celerio euphorbiae during pupal development with diapause. J Insect Physiol. 1971, 17: 1961-1968. 10.1016/0022-1910(71)90140-5.
Van Der Kloot W: The control of neurosecretion and diapause by physiological changes in the brain of the cecropia silkworm. Biol Bull. 1955, 109: 276-294. 10.2307/1538728.
Mansingh A, Smallman BN: The cholinergic system in insect diapause. J Insect Physiol. 1967, 13: 447-467. 10.1016/0022-1910(67)90084-4.
Shappirio D, Eichenbaum D, Locke B: Cholinesterase in the brain of the cecropia silkmoth during metamorphosis and pupal diapause. Biol Bull. 1967, 132: 108-125. 10.2307/1539881.
Zhang Y, Yan H, Lu W, Li Y, Guo X, Xu B: A novel Omega-class glutathione S-transferase gene in Apis cerana cerana: molecular characterisation of GSTO2 and its protective effects in oxidative stress. Cell Stress Chaperones. 2013, 18: 503-516. 10.1007/s12192-013-0406-2.
Jovanović-Galović A, Blagojevic D, Grubor-Lajsic G, Worland M, Spasic M: Antioxidant defense in mitochondria during diapause and postdiapause development of European corn borer (Ostrinia nubilalis, Hubn.). Arch Insect Biochem Physiol. 2007, 64: 111-119. 10.1002/arch.20160.
Jovanović-Galović A, Blagojević DP, Grubor-Lajsić G, Worland R, Spasić MB: Role of antioxidant defense during different stages of preadult life cycle in European corn borer (Ostrinia nubilalis, Hubn.): Diapause and metamorphosis. Arch Insect Biochem Physiol. 2004, 55: 79-89. 10.1002/arch.10126.
Dermauw W, Osborne EJ, Clark RM, Grbi M, Tirry L, Van Leeuwen T: A burst of ABC genes in the genome of the polyphagous spider mite Tetranychus urticae. BMC Genomics. 2013, 14: 317-10.1186/1471-2164-14-317.
Pohl PC, Klafke GM, Carvalho DD, Martins JR, Daffre S, da Silva Vaz I, Masuda A: ABC transporter efflux pumps: a defense mechanism against ivermectin in Rhipicephalus (Boophilus) microplus. Int J Parasitol. 2011, 41: 1323-1333. 10.1016/j.ijpara.2011.08.004.
Yabe T, Suzuki N, Furukawa T, Ishihara T, Katsura I: Multidrug resistance-associated protein MRP-1 regulates dauer diapause by its export activity in Caenorhabditis elegans. Development. 2005, 132: 3197-3207. 10.1242/dev.01909.
Hock T, Cottrill T, Keegan J, Garza D: The E23 early gene of Drosophila encodes an ecdysone-inducible ATP-binding cassette transporter capable of repressing ecdysone-mediated gene activation. Proc Natl Acad Sci USA. 2000, 97: 9519-9524. 10.1073/pnas.160271797.
Itoh TQ, Tanimura T, Matsumoto A: Membrane-bound transporter controls the circadian transcription of clock genes in Drosophila. Genes Cells. 2011, 16: 1159-1167. 10.1111/j.1365-2443.2011.01559.x.
Qin W, Neal SJ, Robertson RM, Westwood JT, Walker VK: Cold hardening and transcriptional change in Drosophila melanogaster. Insect Mol Biol. 2005, 14: 607-613. 10.1111/j.1365-2583.2005.00589.x.
You M, Yue Z, He W, Yang X, Yang G, Xie M, Zhan D, Baxter SW, Vasseur L, Gurr GM, Douglas CJ, Bai J, Wang P, Cui K, Huang S, Li X, Zhou Q, Wu Z, Chen Q, Liu C, Wang B, Li X, Xu X, Lu C, Hu M, Davey JW, Smith SM, Chen M, Xia X, Tang W: A heterozygous moth genome provides insights into herbivory and detoxification. Nat Genet. 2013, 45: 220-225. 10.1038/ng.2524.
Broehan G, Kroeger T, Lorenzen M, Merzendorfer H: Functional analysis of the ATP-binding cassette (ABC) transporter gene family of Tribolium castaneum. BMC Genomics. 2013, 14: 6-10.1186/1471-2164-14-6.
Chino C: Conversion of glycogen to sorbitol and glycerol in the diapause egg of the Bombyx silkworm. Nature. 1957, 180: 606-607.
Sømme L: Supercooling and winter survival in terrestrial arthropods. Comp Biochem Physiol Part A Physiol. 1982, 73A: 519-543.
Miller L, Smith JS: Production of threitol and sorbitol by an adult insect: association with freezing tolerance. Nature. 1975, 258: 519-520. 10.1038/258519a0.
Michaud MR, Denlinger DL: Shifts in the carbohydrate, polyol, and amino acid pools during rapid cold-hardening and diapause-associated cold-hardening in flesh flies (Sarcophaga crassipalpis): a metabolomic comparison. J Comp Physiol B. 2007, 177: 753-763. 10.1007/s00360-007-0172-5.
Kukal O, Duman JG, Serianni AS: Glycerol metabolism in a freeze-tolerant arctic insect: An in vivo 13-C NMR study. J Comp Physiol B. 1988, 158: 175-183. 10.1007/BF01075831.
Kukal O, Ayres MP, Scriber JM: Cold tolerance of the pupae in relation to the distribution of swallowtail butterflies. Can J Zool. 1991, 69: 3028-3037. 10.1139/z91-427.
Slotsbo S, Maraldo K, Malmendal A, Nielsen NC, Holmstrup M: Freeze tolerance and accumulation of cryoprotectants in the enchytraeid Enchytraeus albidus (Oligochaeta) from Greenland and Europe. Cryobiology. 2008, 57: 286-291. 10.1016/j.cryobiol.2008.09.010.
Storey K, Storey JM: Freeze tolerant frogs: cryo-protectants and tissue metabolism during freeze-thaw cycles. Can J Zool. 1986, 64: 49-56. 10.1139/z86-008.
Hoshikawa K: Accumulation of inositol by hibernating adults of Coccinellid and Chrysomelid beetles. Low Temp Sci. 1981, 39: 45-48.
Kostál V, Zahradnícková H, Simek P, Zelený J: Multiple component system of sugars and polyols in the overwintering spruce bark beetle, Ips typographus. J Insect Physiol. 2007, 53: 580-586. 10.1016/j.jinsphys.2007.02.009.
Saeidi F, Moharramipour S, Barzegar M: Seasonal patterns of cold hardiness and cryoprotectant profiles in Brevicoryne brassicae (Hemiptera: Aphididae). Environ Entomol. 2012, 41: 1638-1643. 10.1603/EN12127.
Wang H-S, Kang L: Effect of cooling rates on the cold hardiness and cryoprotectant profiles of locust eggs. Cryobiology. 2005, 51: 220-229. 10.1016/j.cryobiol.2005.07.003.
Tanaka K: Seasonal change in glycogen and inositol/sorbitol contents of the house spider, Achaearanea tepidariorum (Araneae: Theridiidae). Comp Biochem Physiol Part B Biochem Mol Biol. 1995, 110: 539-545. 10.1016/0305-0491(94)00183-U.
Bayley M: Water vapor absorption in Arthropods by accumulation of myoinositol and glucose. Science (80- ). 1999, 285: 1909-1911. 10.1126/science.285.5435.1909.
Colinet H, Nguyen TTA, Cloutier C, Michaud D, Hance T: Proteomic profiling of a parasitic wasp exposed to constant and fluctuating cold exposure. Insect Biochem Mol Biol. 2007, 37: 1177-1188. 10.1016/j.ibmb.2007.07.004.
Overgaard J, Malmendal A, Sørensen JG, Bundy JG, Loeschcke V, Nielsen NC, Holmstrup M: Metabolomic profiling of rapid cold hardening and cold shock in Drosophila melanogaster. J Insect Physiol. 2007, 53: 1218-1232. 10.1016/j.jinsphys.2007.06.012.
Hubálek Z: Protectants used in the cryopreservation of microorganisms. Cryobiology. 2003, 46: 205-229. 10.1016/S0011-2240(03)00046-4.
Marum L, Estêvão C, Oliveira MM, Amâncio S, Rodrigues L, Miguel C: Recovery of cryopreserved embryogenic cultures of maritime pine - effect of cryoprotectant and suspension density. Cryo Letters. 2004, 25: 363-374.
Kandaswamy KK, Chou K-C, Martinetz T, Möller S, Suganthan PN, Sridharan S, Pugalenthi G: AFP-Pred: A random forest approach for predicting antifreeze proteins from sequence-derived properties. J Theor Biol. 2011, 270: 56-62. 10.1016/j.jtbi.2010.10.037.
Graham LA, Liou YC, Walker VK, Davies PL: Hyperactive antifreeze protein from beetles. Nature. 1997, 388: 727-728. 10.1038/41908.
Liou YC, Tocilj A, Davies PL, Jia Z: Mimicry of ice structure by surface hydroxyls and water of a beta-helix antifreeze protein. Nature. 2000, 406: 322-324. 10.1038/35018604.
Graether SP, Kuiper MJ, Gagné SM, Walker VK, Jia Z, Sykes BD, Davies PL: β-Helix structure and ice-binding properties of a hyperactive antifreeze protein from an insect. Nature. 2000, 406: 325-328. 10.1038/35018610.
Garnham CP, Campbell RL, Walker VK, Davies PL: Novel dimeric β-helical model of an ice nucleation protein with bridged active sites. BMC Struct Biol. 2011, 11: 36-10.1186/1472-6807-11-36.
Zhang D-Q, Liu B, Feng D-R, He Y-M, Wang S-Q, Wang H-B, Wang J-F: Significance of conservative asparagine residues in the thermal hysteresis activity of carrot antifreeze protein. Biochem J. 2004, 377 (Pt 3): 589-595.
Marshall C: Identification of the ice-binding face of antifreeze protein from Tenebrio molitor. FEBS Lett. 2002, 529: 261-267. 10.1016/S0014-5793(02)03355-0.
Scotter AJ, Marshall CB, Graham L a, Gilbert J a, Garnham CP, Davies PL: The basis for hyperactivity of antifreeze proteins. Cryobiology. 2006, 53: 229-239. 10.1016/j.cryobiol.2006.06.006.
Graham LA, Qin W, Lougheed SC, Davies PL, Walker VK: Evolution of hyperactive, repetitive antifreeze proteins in beetles. J Mol Evol. 2007, 64: 387-398. 10.1007/s00239-005-0256-3.
Hakim A, Nguyen JB, Basu K, Zhu DF, Thakral D, Davies PL, Isaacs FJ, Modis Y, Meng W: Crystal structure of an insect antifreeze protein and its implications for ice binding. J Biol Chem. 2013, 288: 12295-12304. 10.1074/jbc.M113.450973.
Lin F-H, Graham L a, Campbell RL, Davies PL: Structural modeling of snow flea antifreeze protein. Biophys J. 2007, 92: 1717-1723. 10.1529/biophysj.106.093435.
Nicodemus J, O’tousa JE, Duman JG: Expression of a beetle, Dendroides canadensis, antifreeze protein in Drosophila melanogaster. J Insect Physiol. 2006, 52: 888-896. 10.1016/j.jinsphys.2006.05.009.
Lin X, Wisniewski ME, Duman JG: Expression of two self-enhancing Antifreeze Proteins from the beetle Dendroides canadensis in Arabidopsis thaliana. Plant Mol Biol Report. 2011, 29: 802-813. 10.1007/s11105-011-0287-4.
Venketesh S, Dayananda C: Properties, potentials, and prospects of antifreeze proteins. Crit Rev Biotechnol. 2008, 28: 57-82. 10.1080/07388550801891152.
Barrett J: Thermal hysteresis proteins. Int J Biochem Cell Biol. 2001, 33: 105-117. 10.1016/S1357-2725(00)00083-2.
Duman JG: Antifreeze and ice nucleator proteins in terrestrial arthropods. Annu Rev Physiol. 2001, 63: 327-357. 10.1146/annurev.physiol.63.1.327.
Graham L, Davies P: Glycine-rich antifreeze proteins from snow fleas. Science. 2005, 310: 461-10.1126/science.1115145.
Qin W, Doucet D, Tyshenko MG, Walker VK: Transcription of antifreeze protein genes in Choristoneura fumiferana. Insect Mol Biol. 2007, 16: 423-434. 10.1111/j.1365-2583.2007.00743.x.
Duman JG, Li N, Verleye D, Goetz FW, Wu DW, Andorfer C a, Benjamin T, Parmelee DC: Molecular characterization and sequencing of antifreeze proteins from larvae of the beetle Dendroides canadensis. J Comp Physiol B. 1998, 168: 225-232. 10.1007/s003600050140.
Kristiansen E, Pedersen S, Ramløv H, Zachariassen KE: Antifreeze activity in the cerambycid beetle Rhagium inquisitor. J Comp Physiol B. 1999, 169: 55-60. 10.1007/s003600050193.
Kim M, Robich RM, Rinehart JP, Denlinger DL: Upregulation of two actin genes and redistribution of actin during diapause and cold stress in the northern house mosquito, Culex pipiens. J Insect Physiol. 2006, 52: 1226-1233. 10.1016/j.jinsphys.2006.09.007.
Robich RM, Rinehart JP, Kitchen LJ, Denlinger DL: Diapause-specific gene expression in the northern house mosquito, Culex pipiens L., identified by suppressive subtractive hybridization. J Insect Physiol. 2007, 53: 235-245. 10.1016/j.jinsphys.2006.08.008.
George E, Ober M, Emerson C: Functional domains of the Drosophila melanogaster muscle myosin heavy-chain gene are encoded by alternatively spliced exons. Mol Cell Biol. 1989, 9: 2957-2974.
Wells L, Edwards K, Bernstein S: Myosin heavy chain isoforms regulate muscle function but not myofibril assembly. EMBO J. 1996, 15: 4454-4459.
Lowey S, Waller G, Trybus K: Function of skeletal muscle myosin heavy and light chain isoforms by an in vitro motility assay. J Biol Chem. 1993, 268: 20414-20416.
Schwartz L, Jones M, Kosz L, Kuah K: Selective repression of actin and myosin heavy chain expression during the programmed death of insect skeletal muscle. Dev Biol. 1993, 158: 448-455. 10.1006/dbio.1993.1202.
Whiteley N, Taylor E, Haj AE: Seasonal and latitudinal adaptation to temperature in crustaceans. J Therm Biol. 1997, 22: 419-427. 10.1016/S0306-4565(97)00061-2.
Qiu Z, Tsoi SCM, MacRae TH: Gene expression in diapause-destined embryos of the crustacean, Artemia franciscana. Mech Dev. 2007, 124: 856-867. 10.1016/j.mod.2007.09.001.
Rosenberg-Hasson Y, Renert-Pasca M, Volk T: A Drosophila dystrophin-related protein, MSP-300, is required for embryonic muscle morphogenesis. Mech Dev. 1996, 60: 83-94. 10.1016/S0925-4773(96)00602-8.
Veerman A: Carotenoid metabolism in Tetranychus urticae koch (Acari: Tetranychidae). Comp Biochem Physiol Part B Comp Biochem. 1974, 47: 101-116. 10.1016/0305-0491(74)90095-9.
Veerman A, Helle W: Evidence for the functional involvement of carotenoids in the photoperiodic reaction of spider mites. Nature. 1978, 275: 234-
Veerman A: Functional involvement of carotenoids in photoperiodic induction of diapause in the spider mite, Tetranychus urticae. Physiol Entomol. 1980, 5: 291-300. 10.1111/j.1365-3032.1980.tb00237.x.
Bosse TC, Veerman A: Involvement of vitamin A in the photoperiodic induction of diapause in the spider mite Tetranychus urticae is demonstrated by rearing an albino mutant on a semi-synthetic diet with and without β-carotene or vitamin A. Physiol Entomol. 1996, 21: 188-192. 10.1111/j.1365-3032.1996.tb00854.x.
Moran N, Jarvik T: Lateral transfer of genes from fungi underlies carotenoid production in aphids. Science. 2010, 328: 624-627. 10.1126/science.1187113.
Nováková E, Moran N: Diversification of genes for carotenoid biosynthesis in aphids following an ancient transfer from a fungus. Mol Biol Evol. 2012, 29: 313-323. 10.1093/molbev/msr206.
Cobbs C, Heath J, Stireman JO, Abbot P: Carotenoids in unexpected places: Gall midges, lateral gene transfer, and carotenoid biosynthesis in animals. Mol Phylogenet Evol. 2013, 68: 221-228. 10.1016/j.ympev.2013.03.012.
Acuña R, Padilla BE, Flórez-Ramos CP, Rubio JD, Herrera JC, Benavides P, Lee S-J, Yeats TH, Egan AN, Doyle JJ, Rose JKC: Adaptive horizontal transfer of a bacterial gene to an invasive insect pest of coffee. Proc Natl Acad Sci U S A. 2012, 109: 4197-4202. 10.1073/pnas.1121190109.
Wybouw N, Balabanidou V, Ballhorn DJ, Dermauw W, Grbić M, Vontas J, Van Leeuwen T: A horizontally transferred cyanase gene in the spider mite Tetranychus urticae is involved in cyanate metabolism and is differentially expressed upon host plant change. Insect Biochem Mol Biol. 2012, 42: 881-889. 10.1016/j.ibmb.2012.08.002.
Altincicek B, Kovacs JL, Gerardo NM:Horizontally transferred fungal carotenoid genes in the two-spotted spider miteTetranychus urticae.Biol Lett. 2012, 8: 253-257. 10.1098/rsbl.2011.0704.
Heath JJ, Cipollini DF, Stireman JO: The role of carotenoids and their derivatives in mediating interactions between insects and their environment. Arthropod Plant Interact. 2012, 7: 1-20.
Veerman A: Photoperiodic time measurement in insects and mites: a critical evaluation of the oscillator-clock hypothesis. J Insect Physiol. 2001, 47: 1097-1109. 10.1016/S0022-1910(01)00106-8.
Suzuki T, Watanabe M, Takeda M: UV tolerance in the two-spotted spider mite, Tetranychus urticae. J Insect Physiol. 2009, 55: 649-654. 10.1016/j.jinsphys.2009.04.005.
Rinehart JP, Robich RM, Denlinger DL:Isolation of diapause-regulated genes from the flesh fly,Sarcophaga crassipalpisby suppressive subtractive hybridization. J Insect Physiol. 2009, 56: 603-609.
Sim C, Denlinger DL: Catalase and superoxide dismutase-2 enhance survival and protect ovaries during overwintering diapause in the mosquito Culex pipiens. J Insect Physiol. 2011, 57: 628-634. 10.1016/j.jinsphys.2011.01.012.
Bao B, Xu W-H: Identification of gene expression changes associated with the initiation of diapause in the brain of the cotton bollworm, Helicoverpa armigera. BMC Genomics. 2011, 12: 224-10.1186/1471-2164-12-224.
Veenstra JA, Rombauts S, Grbić M: In silico cloning of genes encoding neuropeptides, neurohormones and their putative G-protein coupled receptors in a spider mite. Insect Biochem Mol Biol. 2012, 42: 277-295. 10.1016/j.ibmb.2011.12.009.
Clark K a J, Brierley AS, Pond DW, Smith VJ: Changes in seasonal expression patterns of ecdysone receptor, retinoid X receptor and an A-type allatostatin in the copepod, Calanus finmarchicus, in a sea loch environment: an investigation of possible mediators of diapause. Gen Comp Endocrinol. 2013, 189: 66-73.
Yamashita O: Diapause hormone of the silkworm, Bombyx mori: structure, gene expression and function. J Insect Physiol. 1996, 42 (I): 669-679.
Homma T, Watanabe K, Tsurumaru S, Kataoka H, Imai K, Kamba M, Niimi T, Yamashita O, Yaginuma T: G protein-coupled receptor for diapause hormone, an inducer of Bombyx embryonic diapause. Biochem Biophys Res Commun. 2006, 344: 386-393. 10.1016/j.bbrc.2006.03.085.
Nusawardani T, Kroemer J, Choi M-Y, Jurenka R: Identification and characterization of the pyrokinin/pheromone biosynthesis activating neuropeptide family of G protein-coupled receptors from Ostrinia nubilalis. Insect Mol Biol. 2013, 22: 331-340. 10.1111/imb.12025.
Hall T: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999, 41: 95-98.
TVL and WD are post-doctoral fellows of the Fund for Scientific Research Flanders (FWO). This work was supported by FWO grant 3G061011 (to TVL) and 3G009312 (to TVL and LT), a Ghent University Special Research Fund grant 01J13711 (to TVL and LT), the Government of Canada through Genome Canada and the Ontario Genomics Institute OGI-046. Nicky Wybouw is supported by the Institute for the Promotion of Innovation by Science and Technology in Flanders (IWT, grant IWT/SB/101451). We thank dr. ir. Bartel Vanholme for his valuable advice on interpreting and editing the 3D model, Prof. dr. ir. Björn Menten for providing the microarray facilities and ir. Jeroen Crappé for running the program AFP-Pred.
The authors declare that they have no competing interests.
AB and TVL designed research. AB, WD and NW performed the experiments. AB analyzed data with input from WD, TVL and NW. AB and TVL wrote the manuscript, with input from WD, NW and LT. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 3:Microarray validation by qPCR. qPCR validation was performed for 13 genes that were differentially expressed in our microarray experiment (tetur23g00860 (C1A cysteine peptidase), tetur03g09330 (glucose dehydrogenase), tetur26g00570 (low-density lipoprotein receptor), tetur17g03230 (lipase A), tetur22g02640, tetur22g02670, tetur22g02690, tetur22g02730, tetur22g02790 (antifreeze proteins), tetur01g11260, tetur11g04820 (carotenoid desaturases), tetur01g11270 and tetur11g04840 (carotenoid synthases). Error bars represent the standard error of the calculated mean based on three biological replicates. Asterisks indicate significantly differentially expressed genes (random reallocation test) compared to the reference condition (green, non-diapausing LS-VL strain at 17°C). Microarray expression data from this selection of genes are shown next to their qPCR expression data. qPCR primers can be found in Additional file 1. (TIFF 2 MB)
Additional file 4:Blast2GO data distribution of differentially expressed genes in diapausing T. urticae females and of all protein coding genes in the T. urticae genome.(DOCX 14 KB)
Additional file 6:Differentially expressed cytochrome P450 monooxygenases (CYPs) in diapausing T. urticae females.(DOCX 22 KB)
Additional file 7:Differentially expressed carboxyl/cholinesterases (CCEs) in diapausing T. urticae females.(DOCX 21 KB)
Additional file 8:Differentially expressed glutathione S-transferases (GSTs) in diapausing T. urticae females.(DOCX 18 KB)
Additional file 10:Differentially expressed intradiol ring-cleavage dioxygenases (ID-RCDs) in diapausing T. urticae females.(DOCX 18 KB)
Additional file 14:Alignment of T. urticae antifreeze proteins (AFPs) with those of T. molitor (pdb: 1L1I_A) and D. canadiensis (AFP1A and AFP7 (GenBank acc. numbers AAB94303.1 and AAF86362.1, respectively). AFP protein sequences, without their predicted signal peptides, were aligned using MUSCLE  and the resulting alignment was graphically edited using BioEdit version 7.2 . (PDF 78 KB)
Additional file 15:Predicted 3D structure of a putative T. urticae AFP, tetur22g03033. (best viewed with Swiss-Pdb-Viewer (http://spdbv.vital-it.ch).(PDB 38 KB)
Additional file 16:Differentially expressed genes involved in carotenoid synthesis in diapausing T. urticae females.(DOCX 16 KB)
Additional file 17:Differentially expressed genes involved in signal transduction in diapausing T. urticae females.(DOCX 23 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Bryon, A., Wybouw, N., Dermauw, W. et al. Genome wide gene-expression analysis of facultative reproductive diapause in the two-spotted spider mite Tetranychus urticae. BMC Genomics 14, 815 (2013). https://doi.org/10.1186/1471-2164-14-815
- Thermal hysteresis
- Ice binding