- Research article
- Open Access
Effects of cold-acclimation on gene expression in Fall field cricket (Gryllus pennsylvanicus) ionoregulatory tissues
BMC Genomics volume 18, Article number: 357 (2017)
Cold tolerance is a key determinant of temperate insect distribution and performance. Chill-susceptible insects lose ion and water homeostasis during cold exposure, but prior cold acclimation improves both cold tolerance and defense of homeostasis. The mechanisms underlying these processes are mostly unknown; cold acclimation is thought to enhance ion transport in the cold and/or prevent leak of water and ions. To identify candidate mechanisms of cold tolerance plasticity we generated transcriptomes of ionoregulatory tissues (hindgut and Malpighian tubules) from Gryllus pennsylvanicus crickets and compared gene expression in warm- and cold-acclimated individuals.
We assembled a G. pennsylvanicus transcriptome de novo from 286 million 50-bp reads, yielding 70,037 contigs (~44% of which had putative BLAST identities). We compared the transcriptomes of warm- and cold-acclimated hindguts and Malpighian tubules. Cold acclimation led to a ≥ 2-fold change in the expression of 1493 hindgut genes (733 downregulated, 760 upregulated) and 2008 Malpighian tubule genes (1009 downregulated, 999 upregulated). Cold-acclimated crickets had altered expression of genes putatively associated with ion and water balance, including: a downregulation of V-ATPase and carbonic anhydrase in the Malpighian tubules and an upregulation of Na+-K+ ATPase in the hindgut. We also observed acclimation-related shifts in the expression of cytoskeletal genes in the hindgut, including actin and actin-anchoring/stabilizing proteins, tubulin, α-actinin, and genes involved in adherens junctions organization. In both tissues, cold acclimation led to differential expression of genes encoding cytochrome P450s, glutathione-S-transferases, apoptosis factors, DNA repair, and heat shock proteins.
This is the first G. pennsylvanicus transcriptome, and our tissue-specific approach yielded new candidate mechanisms of cold tolerance plasticity. Cold acclimation may reduce loss of hemolymph volume in the cold by 1) decreasing primary urine production via reduced expression of carbonic anhydrase and V-ATPase in the Malpighian tubules and 2) by increasing Na+ (and therefore water) reabsorption across the hindgut via increase in Na+-K+ ATPase expression. Cold acclimation may reduce chilling injury by remodeling and stabilizing the hindgut epithelial cytoskeleton and cell-to-cell junctions, and by increasing the expression of genes involved in DNA repair, detoxification, and protein chaperones.
Most insects are chill-susceptible, such that their thermal performance and survival are limited in the cold at temperatures well above the freezing point . Although ice formation causes direct injury, the mechanisms of chilling injury (i.e those not associated with ice formation) are not well-understood. Direct cold shock probably causes immediate damage to cells, for example via phase transition in membranes [2,3,4], disruption of the cytoskeleton [4,5,6,7,8], or induction of apoptosis . Indirect chilling injury accumulates over time, most likely resulting from a loss of ion and water balance in the cold [10,11,12,13,14] (although there is also evidence of roles for oxidative damage and disruption of signalling pathways [5, 15, 16]). Cold-acclimated insects maintain water and ion homeostasis and avoid chilling injury and mortality to lower temperatures than warm-acclimated conspecifics [17,18,19,20,21,22], although the underlying mechanisms are not completely understood.
Chill-susceptible insects lose hemolymph water and Na+ to the gut lumen during cold exposure [13, 23,24,25] and re-establish this homeostasis during recovery upon rewarming . Ion and water homeostasis in insects is primarily regulated at the Malpighian tubules and hindgut . The distal Malpighian tubules actively transport ions across a leaky epithelium to drive secretion of water, metabolic wastes, and other ions into the tubule lumen. This primary urine – which is isosmotic to the hemolymph  – is partially modified at the proximal tubule (a tight epithelium) prior to entering the gut lumen [28,29,30]. Water and ions are then selectively reabsorbed from the gut lumen by the hindgut (particularly at the rectum) . Na+-K+ ATPase (NKA) maintains high paracellular [Na+] in the rectal epithelium, driving paracellular migration of water (and concurrent reabsorption of some Na+ and Cl-) from the gut lumen to the hemolymph. Secretion and reabsorption are regulated by a suite of diuretic and antidiuretic peptides (see ), and these peptides may be important for recovery from cold stress in Drosophila . Loss and recovery of ion and water balance in the cold is likely dependent upon processes at the Malpighian tubule and hindgut epithelia; specifically, enhanced ion pumping rate at low temperatures could maintain ion balance for longer, and decreased epithelial permeability could reduce the rate of ion leakage.
Transport enzyme function is temperature-dependent , thus cold exposure should limit ion pumping rates while the rate of passive leak should remain relatively unchanged . To compensate, insects could increase transport capacity by expressing or mobilizing more transport enzymes [2, 36,37,38]. An obvious hypothesis, then, is that expression or activity of NKA and other transport enzymes in the Malpighian tubules and hindgut may be targeted for modification during cold acclimation. Water and ion leak across the ionoregulatory epithelia during cold exposure will depend on tissue permeability [13, 35, 39], and paracellular pathways may be especially important as these shunts are the primary pathway for water and anion movement across both Malpighian tubule and rectal epithelia [28, 31, 40, 41]. Because paracellular shunts are plastic and may be temperature-sensitive [42,43,44,45], they are likely targets for modification during cold acclimation to reduce epithelial permeability.
Surprisingly few genes have been directly associated with cold tolerance in insects [2, 46], but transcriptomics approaches have revealed many candidates. For example, modified expression of cuticular genes appears to underlie cold tolerance variation among New Zealand stick insects [47, 48], and the gene Frost has been associated with recovery from cold exposure in Drosophila [49,50,51,52,53], although its function remains elusive . Cold shock recovery in Sarcophaga bullata flesh flies changes the expression of genes related to the membrane and cytoskeletal structure, apoptosis, protein folding, oxidative stress, and signaling , and many of these genes (in addition to those involving autophagy and ion transport) are also modified with cold acclimation and rapid cold-hardening in D. melanogaster . Most transcriptomic studies have explored responses to acute cold exposure (e.g. [52, 57]), or compared natural variation among populations and species whose underlying population differences may render the specific differences that drive cold-related phenotypes difficult to detect (e.g. ). Acclimation of a single population is therefore a useful approach to identify candidate mechanisms associated specifically with plasticity of cold tolerance.
A few single - population transcriptomic studies have identified transcriptomic changes associated with cold acclimation, e.g. MacMillan et al.  and Gerken et al.  recently identified key pathways and 1000s of genes associated with cold acclimation in Drosophila melanogaster. In these Drosophila datasets, modification of ion transport (particularly differential expression of Na+ transporters ) and cellular adhesion is consistent with our expectation that modulating epithelial transport is associated with cold acclimation. However, these and other related studies have examined the transcriptomes of either entire animals (e.g. [56, 58]), or mixed tissues (e.g. the entire excised heads of stick insects [47, 48]). Because the Malpighian tubules and hindgut effectively work antagonistically in insect ion and water balance, transcriptomic shifts in these epithelia are likely to be masked in whole-animal homogenates. Thus, a tissue-specific approach to transcriptomics is urgently needed to more precisely determine the cellular- and tissue-level changes underlying cold acclimation in chill-susceptible insects.
The fall field cricket, Gryllus pennsylvanicus (Orthoptera: Gryllidae) is a generalist omnivore of grassy habitats across the Eastern North American temperate zone . The species is univoltine and overwinters in diapause in the soil as an egg . Adult G. pennsylvanicus are chill-susceptible: they develop chilling injuries in as little as 12 h at 0 °C, and are killed by 3–5 d at this temperature . Gryllus pennsylvanicus exhibits plasticity in cold tolerance, and has emerged as a model for understanding cold-induced loss of ion and water balance [13, 39, 61]. Briefly, when these crickets are exposed to cold, Na+, Ca2+, Mg2+, and water migrate from the hemolymph to the gut, hemolymph [K+] rises, and muscle equilibrium potential is lost [13, 39, 61]. Ion and water balance are actively re-established during recovery from cold . Cold-acclimated G. pennsylvanicus have improved defense of ion and water homeostasis in the cold, a lowered critical thermal minimum, faster chill coma recovery time, and suffer lower rates of injury and mortality following cold shock . Thus, the G. pennsylvanicus system is well-suited for exploring mechanisms of cold tolerance plasticity.
Here we took a tissue-specific comparative gene expression approach to understand the processes of cold acclimation in the transporting epithelia of chill-susceptible insects. We assembled a transcriptome for G. pennsylvanicus and compared the expression of Malpighian tubule and hindgut genes between warm- and cold-acclimated adults (with a focus on genes involved in ion and water homeostasis and cellular and junctional integrity). We aimed to generate mechanistic hypotheses about specific molecular underpinnings of cold acclimation, and provide insights into the causes of water and ion disruption during cold exposure.
Our colony of G. pennsylvanicus originated from individuals collected in 2004 from the University of Toronto at Mississauga campus, Ontario and was reared under constant summer-like conditions (25 °C, 14 L:10 D photoperiod, 70% RH), as described previously [13, 39]. At approximately 8 weeks post-hatch, and prior to sexual maturation, female crickets were separated from males to prevent mating. Adult females at approximately 11 weeks post-hatch were used for all experiments.
Acclimation and dissection
Crickets were first isolated individually in common summer-like conditions (25 °C, 14 L:10D photoperiod) in mesh-covered 177 mL transparent cups (Polar Plastics, Summit Food Distributors, London, ON, Canada) containing egg carton shelters, rabbit food, and water. Isolation prevented cannibalism and lasted 1 week. We then haphazardly assigned crickets into cold- and warm-acclimations (n = 30 per treatment). For warm acclimation, crickets remained in the rearing growth chamber under constant summer-like conditions. We cold-acclimated crickets in a Sanyo MIR 154 incubator (Sanyo Scientific, Bensenville, Illinois) at 10 L:14 D with temperature decreasing from 25 to 12 °C over 7 days followed by constant 12 °C for 3 weeks. This regime lowers the critical thermal minimum c. 2 °C, reduces chill coma recovery time by c. 65%, increases survival following 5 d at 0 °C by nearly 80%, and enhances maintenance of ion and water homeostasis in the cold . While the two acclimation temperatures are likely to affect physiological ageing, the 4 weeks of acclimation represent approximately 20% of the adult cricket lifespan. Therefore we assume that the effect of physiological age on gene expression will be less than the effect of acclimation temperature .
Cricket hindguts were dissected as described previously  from live crickets immediately following the 4 weeks of warm- or cold-acclimation. Under Ringer’s solution in a Petri dish the hindgut (rectum, colon, and ileum) was cut away from the gastrointestinal tract and flushed of fecal material with approximately 3 mL of Ringer’s using a syringe (this procedure took < 3 min). Malpighian tubules were removed as a single bunch by detaching the ureter with forceps, rinsing briefly in Ringer’s, and blotting on a tissue. Malpighian tubules and hindguts were flash-frozen in liquid nitrogen. Three biological replicates for sequencing of hindgut and Malpighian tubule transcriptomes were each comprised of pooled tissues from ten individuals. Following dissection, crickets were killed by placing them in a freezer at -20 °C. To maximize transcript representation for the de novo assembly, warm- and cold-acclimated whole male and female adult crickets, eggs, and warm-acclimated juveniles were pooled and added to an additional 1.5 mL microcentrifuge tube and flash frozen in liquid nitrogen. All samples were stored at -80 °C until RNA extraction.
RNA extraction & cDNA library preparation
We homogenized thawed tissues with a plastic micropestle (ThermoFisher Scientific, Ottawa ON, Canada) in 1.1 mL TRIzol (Invitrogen, Burlington ON, Canada), and extracted RNA according to manufacturer’s instructions. We purified RNA extracts using the RNeasy Mini kit (Qiagen, Mississauga ON, Canada) according to manufacturer’s instructions, measured absorbance at 260 nm to determine RNA concentrations, and checked for quality with an Agilent Bioanalyzer. cDNA library production and sequencing were performed by the Donnelly Sequencing Center (Toronto ON, Canada). At 13 samples per lane, each cDNA library was sequenced twice using the Illumina HiSeq2500 platform (Illumina, San Diego, CA) with single-end, 50-bp reads.
De novo transcriptome assembly and annotation
We removed Illumina adapter sequences and discarded sequences shorter than 15 nucleotides or containing unknown bases using the Galaxy web service . Sequenced libraries were then grouped and assembled de novo using Trinity release 2012-10-25 [64, 65] on the SHARCNET computing cluster (https://www.sharcnet.ca), with 1 GB Jellyfish Memory and a minimum contig length criterion of 100 nucleotides. We analyzed transcriptome assembly “completeness” compared to a database of arthropod Benchmark Universal Single Copy Orthologs (BUSCO) using BUSCO v1.22 . We compared the contigs from the Trinity assembly to the National Centre for Biotechnology Information (NCBI) non-redundant (nr) protein database (September 2013) by BLASTx (e-value threshold = 1 × 10-3). Gene Ontology (GO) annotation (e-value threshold = 1 × 10-6) was based on SwissProt BLAST matches using Blast2GO version 2.7.2 . To filter out transcriptional artifacts, misassembled transcripts, and poorly supported transcripts, we mapped the original cleaned sequence reads back onto the Trinity-assembly using Bowtie2 version 2.1.0 [68, 69] and reassembled them with the Cufflinks package version 2.1.1 . We used Blast2GO  and the NCBI database to obtain putative identities and GO annotation for mapped transcripts. We accepted one hit for each transcript at an e-value threshold of 1 × 10-3.
Gene expression analyses
We used normalized read counts of genes in warm- and cold-acclimated hindgut and Malpighian tubule libraries for differential gene expression analyses using the edgeR Bioconductor package  in R (v3.2.2, R Development Core Team, 2015 ). Because each biological replicate was sequenced twice (two technical sequencing replicates), read counts from these technical replicate libraries were summed for each gene. For analyses we retained only those genes with at least 10 counts per million in three of the six libraries being compared (warm- vs cold-acclimated hindguts each had three biological replicates) . Filtering yielded 11,140 and 11,066 contigs for differential gene expression analyses of hindgut and Malpighian tubules, respectively. We compared gene expression profiles within tissues (i.e. warm- vs cold-acclimated), and also compared the hindgut and Malpighian tubules for genes that were uniquely up- or downregulated between those tissues with cold acclimation. Individual genes were considered differentially expressed if the absolute fold change was ≥ 2 and if the P-value adjusted for false discovery rate (FDR) was < 0.05.
We used contigs that met our criteria for inclusion (fold change ≥ 2, FDR-adjusted P-value < 0.05) to identify the GO terms associated with the responses to cold acclimation in each tissue (note that we did not formally compare GO terms among tissues or treatments). Differentially-expressed pathways in warm- and cold-acclimated tissues were analyzed using the Kyoto Encyclopedia of Genes and Genomes (KEGG ). KEGG identities were assigned to contigs by the KEGG Automatic Annotation Server , and differential expression analyses of pathway components were performed using the Generally Applicable Gene-set Enrichment (GAGE) and Pathview Bioconductor packages [75, 76] in R. These packages identify coordinated differential expression in gene sets (pre-defined, functionally-related groups of genes) . We accepted pathways as differentially-expressed if the FDR-adjusted P-value was < 0.1.
Sequencing of 26 libraries yielded 286 million 50-bp reads, which we assembled into 70,037 contigs (Additional file 1: Table S1). Our transcriptome included 1808 (67.6%) complete and 415 (15.5%) fragmented arthropod BUSCOs; which is similar to other recent arthropod transcriptome assemblies (e.g. [77, 78]), and comparable to, or better than, the transcriptomes referred to by Simão et al. . Approximately 44% of these contigs in our transcriptome had putative identities by BLAST (Additional file 2: Spreadsheet S1), and of these approximately 36% aligned to genes of the termite Zootermopsis nevadensis. Cold acclimation was associated with a two-fold or greater change in 1493 genes in the hindgut and 2008 genes in the Malpighian tubules (Fig. 1). Within a given tissue, the number of genes up- and down-regulated with cold acclimation were approximately equal. Approximately 52% of all upregulated genes and 60% of all downregulated genes exhibited unique differential expression across the two tissues. Eighty-one genes that appear to be important for cold acclimation (those with a 10-fold or greater change in expression) were unidentifiable by BLAST (Additional file 2: Spreadsheet S1). These represented 22 upregulated and 11 downregulated genes in the hindgut, and 26 upregulated and 22 downregulated genes in the Malpighian tubules.
The GO domain profiles that changed with cold acclimation were similar across the two tissues; of all up- or downregulated GO terms, just over half related to ‘molecular function’ (i.e. transport, binding, enzyme and receptor activities), over one third related to ‘biological processes’ (series of molecular events with a defined beginning and end), and roughly one tenth were ‘cellular components’ (i.e. specific locations of subcellular structures and macromolecular complexes; Additional file 3: Spreadsheet S2). Within the molecular function GO domain, genes involved in nucleotide, protein, metal, and ion binding accounted for over 50% of the upregulated transcripts and approximately 30% of downregulated transcripts in both tissues. Most of the cellular components that were differentially-expressed with cold acclimation involved the cell membrane and extracellular region. In the Malpighian tubules, genes pertaining to the cytoskeleton and cellular junctions accounted for 5% of all upregulated cellular component transcripts, while V-ATPase accounted for 2% of the downregulated transcripts. Metabolic genes accounted for much of the downregulated biological process transcript counts in both tissues. However, some unique differences in biological process profiles between the two tissues were apparent; approximately 17% of downregulated transcripts in the Malpighian tubules were transport-related (2% of which were ion transport-specific), compared to just over 8% of downregulated transcripts in the hindgut. Of the upregulated biological processes over 30% of transcripts in the hindgut involved stress response, protein folding, and repair, while over 10% of transcripts in the Malpighian tubules related to transport.
Cold-acclimated G. pennsylvanicus had altered expression of putative gene orthologs related to apoptosis, the cytoskeleton, detoxification and repair, ion transport and pH regulation, phosphorylation, protein folding, and signal transduction in the hindgut (Tables 1 and 2). A complete list of differentially-expressed genes in the hindgut is provided in the Additional file 2: Spreadsheet S1. Upregulated genes involved in water and ion balance included those encoding atrial natriuretic peptide-converting enzyme, NKA α-subunit, and a Ca2+ release-activated Ca2+ channel protein, while downregulated genes included those putatively encoding bumetanide-sensitive Na+-Cl+ channel (the Na+-K+-2Cl- cotransporter, or NKCC), carbonic anhydrase (CA) isozymes 1 and 9, and a mitochondrial Na+-H+ exchanger (NHA). A gene encoding the putative microtubule-associated protein Jupiter increased during cold acclimation, while a putative enzyme involved in homeoviscous adaptation – ∆9 desaturase 1 – was downregulated 3.9-fold. The genes most differentially-expressed in cold-acclimated crickets related to repair and oxidative damage included those encoding cytochrome P450 (26-fold increase) and a putative cytochrome P450 cyp44 (5.5-fold decrease), glutathione-S-transferase (6.8-fold decrease), and vitellogenin (40-fold decrease). Cold acclimation corresponded with altered expression of hindgut heat shock proteins (hsp 70 and hsp 90 were upregulated, while hsp 67B and hsp β11 were downregulated), and altered expression of some apoptosis genes. Cold-acclimated crickets also exhibited upregulation of a number of protein kinases, phosphodiesterases, and adenylate cyclase in the hindgut. Circadian genes per, clock, and nocturnin were upregulated with cold acclimation, while timeless was downregulated.
More KEGG pathways in the hindgut were downregulated with cold acclimation than were upregulated (Fig. 2). Among 25 upregulated pathways, ‘adherens junction’ (Fig. 3) and ‘gap junction’ are likely to be directly relevant to ion and water balance. Actin regulation within the ‘adherens junction’ pathway was modified; some genes putatively encoding actin-associated proteins (FRG and α-actinin) were upregulated while others were downregulated (β/γ actin, vinculin, and α-catenin). The putative proteins vascular endothelial protein tyrosine phosphatase (VE-PTP), transforming growth factor β2 (TGFβ2), and partitioning defective protein 3 (PAR3) were also upregulated. Upregulation of the ‘gap junction’ pathway was driven by increased expression of a gene putatively encoding tubulin α (TUBA), and to some degree epidermal growth factor receptor (EGFR, or ErbB-1,1 listed as the receptor tyrosine kinase, RTK), while the gene encoding protein kinase C α (PKC-α) was downregulated. Many of the 47 downregulated KEGG pathways in the hindgut were related to metabolism, but also included ‘cardiac muscle contraction’ (Fig. 4) and ‘synaptic vesicle cycle’. Downregulation of the ‘cardiac muscle contraction pathway’ was driven by a decrease in expression of the gene encoding cytochrome c reductase, however there were also significant increases in the expression of genes encoding the NKA α subunit, tropomyosin 1, and myosin heavy chain 6/7. Downregulation of the ‘synaptic vesicle cycle’ was driven by a reduction in the expression of the putative proton pump (V-ATPase), however the gene encoding the dynamin GTPase increased.
As in the hindgut, the Malpighian tubules of cold-acclimated crickets had altered expression of genes related to apoptosis and autophagy, the cytoskeleton, detoxification and repair, ion transport, pH regulation, phosphorylation and signal transduction, and protein folding (Tables 3 and 4). A complete list of differentially-expressed genes in the Malpighian tubules is provided in Additional file 2: Spreadsheet S1. Differentially-expressed genes involved in ion balance and pH regulation included a putative V-ATPase (downregulated 2-fold), Ca2+ and anion transporters and, as in the hindgut, a downregulation of both CA 1 and 9. Cold acclimation led to variable expression of multiple cytoskeletal genes, increased expression of two apoptosis genes, and decreased expression of one gene involved in autophagosome formation. Similar to the hindgut, cold-acclimated Malpighian tubules also exhibited increased expression of hsp 70 and downregulation of hsp 67B, and both up- and downregulation of multiple repair and antioxidant genes (including those encoding putative cytochrome P450s and glutathione-S-transferases). Multiple kinase genes were upregulated in cold-acclimated Malpighian tubules (similar to the hindgut) in addition to a relatively large decrease (6.9-fold) in putative inositol monophosphatase expression. Altered expression of circadian genes following cold acclimation were also similar to that of the hindgut, and juvenile hormone expression was reduced nearly 11-fold.
Similar to patterns in the hindgut, more KEGG pathways were downregulated than were upregulated in cold-acclimated Malpighian tubules (Fig. 2). Many of the 20 upregulated pathways were involved in signaling, and most of the 47 downregulated pathways related to metabolism. The ‘cardiac muscle contraction’ pathway (appropriate to insect striated muscle, ) was downregulated based on reduced expression of a cytochrome c reductase gene. Unlike in the hindgut, NKA, tropomyosin, or myosin heavy chain components of this pathway were not upregulated in Malpighian tubules. The ‘synaptic vesicle cycle’ pathway exhibited downregulation overall (driven by downregulation of V-ATPase), however a number of genes involved in endocytosis and vesicle-membrane fusion were upregulated. These upregulated genes include those encoding putative N-ethylmaleimide-sensitive factor (NSF, an ATPase), dynamin, AP2 complex α (a protein associated with endocytosis of clathrin-coated vesicles), and syntaxin 1A (a protein involved in vesicle fusion for exocytosis).
Cold-acclimated Gryllus pennsylvanicus exhibited modified expression of a range of genes, the functions of which were broadly consistent with differentially-regulated genes associated with cold acclimation and rapid cold-hardening in Drosophila [56, 58, 80]. In crickets, genes involved in stress response, protein folding, and repair were prominent in cold acclimated hindguts, while cold acclimation in the Malpighian tubules was associated with shifts in transport-related genes. In both tissues, cold acclimation was accompanied by altered expression of genes encoding components of the membrane and extracellular space.
Only one gene with known function in insect water homeostasis – that encoding atrial natriuretic peptide-converting enzyme – was upregulated in the hindgut following cold acclimation. The mosquito homologs of this enzyme stimulate primary urine production by the Malpighian tubules by increasing secretion of Na+ and Cl- [32, 81]; the role of this peptide in the insect hindgut is less certain, but increased Na+ and Cl- transport could enhance water reabsorption  and help to defend hemolymph volume during cold exposure. Although some aquaporins have been associated with insect freeze tolerance [83,84,85], their role in cold acclimation among chill-susceptible insects is unknown, and none of the water-transporting insect aquaporins  were differentially-expressed in the hindgut or Malpighian tubules with cold acclimation.
Cold acclimation corresponded with altered expression of putative NKA, NKCC, CA, NHA, and V-ATPase–encoding genes, which are typically enriched in insect transporting epithelia . Most of these gene expression changes were observed in the hindgut. Although all of these transport enzymes contribute to primary urine production by the Malpighian tubules , cold-acclimated Malpighian tubules only exhibited downregulation of genes encoding V-ATPase and CAs 1 and 9.
CA catalyzes the hydration of CO2 to produce H+ + HCO3 -, a source of protons for export by apical V-ATPase by the Malpighian tubules [87, 89, 90], which (in exchange for Na+ or K+ by NHA and K+-H+ antiporters [87, 91]) drives passive excretion of water and anions [28, 92]. Downregulation of CA9 (membrane-bound), CA1 (cytosolic), and V-ATPase in the Malpighian tubules during cold acclimation could therefore have an antidiuretic effect, perhaps defending hemolymph volume in the cold. Indeed, cold-acclimated G. pennsylvanicus Malpighian tubules produce primary urine more slowly (Des Marteaux, Khazraeenia, and Sinclair, unpublished observations).
Paracellular Na+ gradients across the rectal pads drive passive reabsorption of water against osmotic gradients [31, 93, 94]. Failure of NKA to maintain these Na+ gradients during cold exposure could account for leak of Na+, and consequently water, to the gut. We might therefore expect cold-acclimated crickets to compensate for slower enzyme pumping rates at low temperatures by increasing NKA protein abundance; indeed, expression of the α (catalytic) subunit of NKA increased nearly 3-fold in the cricket hindgut after cold acclimation. In D. melanogaster, cold acclimation decreases whole-body NKA activity , suggesting a decoupling of Na+ gradients across the gut rather than compensation for slowed NKA activity in the cold . However, it is unclear how NKA activity changes in D. melanogaster ionoregulatory tissues specifically, so we cannot unequivocally suggest that Diptera and Orthoptera use different acclimation strategies with regards to Na+ balance. NKCC (which imports Na+ , K+, and Cl- basally, energizing apical ion exchange [87, 95]) was downregulated in cold-acclimated cricket hindguts, but the role of NKCC in insect hindgut transport has received little attention, making it difficult to evaluate how NKCC regulation might modify homeostasis.
Cell junctions and structure
Cold acclimation was associated with altered expression of cell and tissue structure-related genes (e.g. hindgut genes involved in cell growth, differentiation, and adhesion included endothelial growth factor [96,97,98]), which could indicate a modified rectal pad epithelium, and which we hypothesise would enhance cold tolerance by minimizing Na+ and water leak in the cold . In addition to possible roles in providing protons for ion transport, CA9 may have roles in cellular adhesion , proliferation, and differentiation (at least in mammals ). Thus, CA9 downregulation in hindgut and Malpighian tubules during cold acclimation may affect epithelial transport by increasing cellular adhesion and epithelial tightness. Tissue-specific post-translational modification of CA9 is an important means of regulating CA9 activity , which would not be captured in a transcriptome comparison such as the present study.
Cold-acclimated crickets had altered expression of hindgut cellular adhesion components associated with both adherens and septate junctions, which comprise the bulk of paracellular connections in the rectal pad [31, 102,103,104], and are closely related . Cold acclimation was also associated with differential expression of genes encoding actin-membrane anchors, which could influence cell junction characteristics or cell shape. We hypothesise that altered actin-membrane anchoring could reduce tension, shearing damage, or unwanted stretch-activation of membrane-bound ion channels [4, 106] when the gut is distended by water migration during cold exposure.
Some gap and tight junction components were also altered during cold acclimation. These will likely modify ion and water recycling between the cytoplasm and paracellular channels [107, 108], and selectivity of absorption [109, 110]. Upregulation of PAR3  and downregulation of PKC-α  suggests increased tight junction formation and therefore increased tightness of the hindgut epithelium during acclimation; we suggest that changes in tight junction morphology in the rectal pads following acclimation could be confirmed by transmission electron microscopy or immunostaining [112, 113]. After cold acclimation, we also observed shifts in the expression of multiple Malpighian tubule genes involving the cytoskeleton and cell junctions (e.g. protein shroom) [114, 115]. Whether these structural changes affect ion and water balance requires some assessment of Malpighian tubule permeability following cold acclimation.
Cold-attributed oxidative stress, disruption of homeostasis and signaling, protein mis-folding, and loss of membrane and cytoskeletal integrity may all contribute to chilling injury and mortality in chill-susceptible insects [12, 15, 116,117,118,119,120]. Cold acclimation was associated with upregulation of putative apoptosis genes (e.g. those encoding the apoptosis initiator caspase 6 and the apoptosis effector caspase 8 ), as well as shifts in autophagy-related gene expression (e.g. upregulation of Ras and ubiquitin signaling KEGG pathways). We hypothesize that the ability to clear cold-damaged cells or cell components is likely increased in cold-acclimated crickets. Polymorphisms or shifts in the expression of genes associated with apoptosis and autophagy appear to be common to the rapid cold-hardening process , and response to dehydration [122, 123] in other insects.
The cytoskeleton depolymerises at low temperatures in fish, mammals, and insects [6, 7, 119, 124], and damage to the cytoskeleton could well be associated with chilling injury in insects. Water loss – which occurs during cold exposure in chill-susceptible species – appears to drive shifts in cytoskeletal gene expression in other insects . Cold-acclimated crickets had altered expression of cytoskeletal branching and stabilizing genes, such as α-catenin, ARP2/3 (a nucleation site for actin polymerization ), and (in the hindgut) tropomyosin, microtubule protein Jupiter and MAP1A/B. Together, this suggests enhanced polymerisation and stabilization of actin and microtubules. Cytoskeleton-related genes were also upregulated in cold-acclimated Culex pipiens , Delia antiqua , alfalfa (Medicago sativa) , and D. melanogaster [58, 126].
Cold exposure appears to cause oxidative stress [16, 118, 127], and cold acclimation is associated with increased activity or expression of antioxidants (e.g. glutathione-S-transferase, catalase, and superoxide dismutase) in a number of insects [16, 46, 58, 128,129,130]. We observed increased expression of some putative DNA repair, and glutathione-S-transferase genes in both the hindgut and Malpighian tubules after cold acclimation (and decreased expression in some of these genes, which we attribute to decreased reactive oxygen species production in the cold [16, 131]). Insects often upregulate antioxidant expression during rewarming [16, 46, 128], so it is possible that a larger suite of antioxidants may become relevant if expression is also modified after cold exposure.
Cold acclimation was associated with increased expression of hsp 70 in both tissues (as well as hsp 90 in the hindgut). Heat shock proteins have a range of cellular-protective roles including as protein chaperones, and could therefore protect against a many aspects of thermal stress and maintain cellular integrity (and therefore perhaps epithelial function) during and after cold exposure. In cold-acclimated D. melanogaster, hsp 70 is one of the few upregulated genes that also increases in protein abundance . Genes encoding several hsps were downregulated, including the less-characterized hsp β11 (involved in the vertebrate heat stress response ) and hsp 67B2 (hsp 67Bb), which may protect against heat and oxidative stress in Drosophila [133, 134].
Other candidate genes
Membrane remodelling is likely an important aspect of cold acclimation [2, 135, 136], and a large proportion of differentially-expressed GO cellular components in the cold-acclimated Malpighian tubules and hindgut were membrane-associated. In the hindgut, a gene encoding ∆9 desaturase (which is involved in homeoviscous adaptation) was downregulated with cold acclimation. By contrast, cold acclimation in the flies Belgica antarctica and Delia radicum is associated with upregulation of ∆9 desaturase [36, 137, 138]. Changes in membrane composition and fluidity with cold acclimation have not been explored in G. pennsylvanicus.
Loss of signal transduction in the cold is a proposed mechanism of chill coma, and recent evidence suggests that paralysis during cold exposure results from direct inhibitory effects of low temperature on neuromuscular function [24, 25, 61, 139, 140]. The expression of some signal transduction and neurotransmission genes (e.g. cAMP, G-proteins, PKC, and a GABA transporter and Ca2+ channel) was altered with cold acclimation in crickets, however it is difficult to predict how these changes would enhance cold tolerance. Perhaps of greater interest were vesicle localization and fusion genes such as dynamin (which mediates membrane-vesicle fusion), NSF (a vesicle-fusing ATPase), AP2 (involved in vesicle endocytosis), and syntaxin 1A (which promotes vesicle-membrane docking), all of which were upregulated with cold acclimation. Vesicle-membrane localization is important for both neurotransmission and recruitment of ion transporters (the latter of which could directly affect ion homeostasis in the cold ), and is likely slowed in the cold.
Cold acclimation appears to affect the expression of some genes associated with circadian rhythm, storage and metabolism, development, and phosphorylation. Differential expression of circadian genes was likely a result of the cold (e.g. low temperatures decrease the expression of timeless and increase the expression of per in flesh flies ). The expression of hindgut genes encoding vitellogenin (a yolk protein precursor) was drastically lower with cold acclimation, and this also occurs in cold-acclimated D. melanogaster ; in this case, we suspect that vitellogenin mRNA came from traces of fat body remaining on the gut during dissections .
We have assembled the first transcriptome of Gryllus pennsylvanicus and our tissue-specific comparative approach yielded precise mechanistic hypotheses about the cold acclimation process. Cold acclimation appears to involve a modification of both ion transport function and cellular/junctional integrity (summarized in Fig. 5). The ion transport modifications likely defend hemolymph volume in the cold; decreased Malpighian tubule V-ATPase and CA expression should slow primary urine production, while upregulation of hindgut NKA should increase Na+ and water reabsorption. Remodeling of the cytoskeleton and adherens junctions (and potentially tight junctions) may mitigate paracellular leak of water and ions in both the hindgut and Malpighian tubules. Cold-acclimated crickets may prevent direct chilling injury by stabilizing the actin cytoskeleton and by changing the way in which actin anchors to the membrane. To repair chilling injuries, cold-acclimated insects may increase expression of antioxidant, DNA repair, apoptosis, autophagy, and chaperone genes. Upregulation of NKA (which should enhance Na+ and water reabsorption across the rectum) may also account for faster chill coma recovery in cold-acclimated insects.
Adaptor-related protein complex 2 subunit α
Actin-related proteins 2 and 3
Carbonic anhydrase isoform 1
Carbonic anhydrase isoform 9
- CTmin :
Critical thermal minimum
Epidermal growth factor receptor
Heat shock protein
Microtubule-associated proteins 1A/1B
Bumetanide-sensitive Na+-K+-2Cl- cotransporter
Partitioning defective 3
Protein kinase C iota
Protein kinase C α
Reactive oxygen species
Receptor tyrosine kinase
Transforming growth factor β2
Vacuolar-type H+ ATPase
Vascular endothelial protein tyrosine phosphatase
Bale JS. Classes of insect cold hardiness. Funct Ecol. 1993;7(6):751–3.
Clark MS, Worland MR. How insects survive the cold: molecular mechanisms—a review. J Comp Physiol B. 2008;178(8):917–33.
Overgaard J, Sørensen JG, Petersen SO, Loeschcke V, Holmstrup M. Changes in membrane lipid composition following rapid cold hardening in Drosophila melanogaster. J Insect Physiol. 2005;51(11):1173–82.
Örvar BL, Sangwan V, Omann F, Dhindsa RS. Early steps in cold sensing by plant cells: the role of actin cytoskeleton and membrane fluidity. Plant J. 2000;23(6):785–94.
Michaud MR, Denlinger DL. Molecular modalities of insect cold survival: current understanding and future trends. Int Congr Ser. 2004;1275:32–46.
Belous A. The role of regulatory systems modulating the state of cytoskeletal proteins under the temperature and osmotic effects. Problems Cryobiol. 1992;4:3–14.
Kayukawa T, Ishikawa Y. Chaperonin contributes to cold hardiness of the onion maggot Delia antiqua through repression of depolymerization of actin at low temperatures. PLoS One. 2009;4(12):e8277.
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(11):1226–33.
Yi S-X, Moore CW, Lee Jr RE. Rapid cold-hardening protects Drosophila melanogaster from cold-induced apoptosis. Apoptosis. 2007;12(7):1183–93.
MacMillan HA, Andersen JL, Davies SA, Overgaard J. The capacity to maintain ion and water homeostasis underlies interspecific variation in Drosophila cold tolerance. Sci Rep. 2015;5:18607.
MacMillan HA, Andersen JL, Loeschcke V, Overgaard J. Sodium distribution predicts the chill tolerance of Drosophila melanogaster raised in different thermal conditions. Am J Physiol. 2015;308(10):R823–31.
MacMillan HA, Baatrup E, Overgaard J. Concurrent effects of cold and hyperkalaemia cause insect chilling injury. Proc R Soc B. 2015;282(1817):20151483.
MacMillan HA, Sinclair BJ. The role of the gut in insect chilling injury: cold-induced disruption of osmoregulation in the fall field cricket, Gryllus pennsylvanicus. J Exp Biol. 2011;214(5):726–34.
MacMillan HA, Williams CM, Staples JF, Sinclair BJ. Reestablishment of ion homeostasis during chill-coma recovery in the cricket Gryllus pennsylvanicus. Proc Natl Acad Sci U S A. 2012;109(50):20750–5.
Teets NM, Yi S-X, Lee RE, Denlinger DL. Calcium signaling mediates cold sensing in insect tissues. Proc Natl Acad Sci U S A. 2013;110(22):9154–9.
Lalouette L, Williams C, Hervant F, Sinclair BJ, Renault D. Metabolic rate and oxidative stress in insects exposed to low temperature thermal fluctuations. Comp Biochem Phys A. 2011;158(2):229–34.
Armstrong GAB, Rodríguez EC, Meldrum Robertson R. Cold hardening modulates K+ homeostasis in the brain of Drosophila melanogaster during chill coma. J Insect Physiol. 2012;58(11):1511–6.
Koštál V, Renault D, Mehrabianová A, Bastl J. Insect cold tolerance and repair of chill-injury at fluctuating thermal regimes: role of ion homeostasis. Comp Biochem Phys A. 2007;147(1):231–8.
Chen C-P, Denlinger DL, Lee RE. Cold-shock injury and rapid cold hardening in the flesh fly Sarcophaga crassipalpis. Physiol Zool. 1987;60(3):297–304.
Coello Alvarado LE. Ion homeostasis and variation in low temperature performance in the fall and spring field cricket (Orthoptera: Gryllidae). London: University of Western Ontario; 2012.
Findsen A, Andersen JL, Calderon S, Overgaard J. Rapid cold hardening improves recovery of ion homeostasis and chill coma recovery time in the migratory locust, Locusta migratoria. J Exp Biol. 2013;216(9):1630–7.
MacMillan HA, Ferguson LV, Nicolai A, Donini A, Staples JF, Sinclair BJ. Parallel ionoregulatory adjustments underlie phenotypic plasticity and evolution of Drosophila cold tolerance. J Exp Biol. 2015;112:2882–7.
Koštál V, Vambera J, Bastl J. On the nature of pre-freeze mortality in insects: water balance, ion homeostasis and energy charge in the adults of Pyrrhocoris apterus. J Exp Biol. 2004;207(9):1509–21.
Koštál V, Yanagimoto M, Bastl J. Chilling-injury and disturbance of ion homeostasis in the coxal muscle of the tropical cockroach (Nauphoeta cinerea). Comp Biochem Phys B. 2006;143(2):171–9.
Findsen A, Pedersen TH, Petersen AG, Nielsen OB, Overgaard J. Why do insects enter and recover from chill coma? Low temperature and high extracellular potassium compromise muscle function in Locusta migratoria. J Exp Biol. 2014;217(8):1297–306.
Phillips J. Comparative physiology of insect renal function. Am J Physiol. 1981;241(5):R241–57.
Ramsay JA. Active transport of water by the Malpighian tubules of the stick insect, Dixippus morosus (Orthoptera, Phasmidae). J Exp Biol. 1954;31(1):104–13.
Beyenbach KW. Transport mechanisms of diuresis in Malpighian tubules of insects. J Exp Biol. 2003;206(21):3845–56.
Beyenbach KW. Mechanism and regulation of electrolyte transport in Malpighian tubules. J Insect Physiol. 1995;41(3):197–207.
Coast GM. The regulation of primary urine production in insects. In: Society for Experimental Biology: 1998. Cambridge: Cambridge University Press; 1998. pp. 189–209.
Phillips JE, Hanrahan J, Chamberlin M, Thomson B. Mechanisms and control of reabsorption in insect hindgut. Adv Insect Physiol. 1987;19:329–422.
Schooley D, Horodyski F, Coast G. Hormones controlling homeostasis in insects. In: Insect endocrinology. Edited by Gilbert L. Waltham: Academic Press; 2012. pp. 366–429.
Terhzaz S, Teets NM, Cabrero P, Henderson L, Ritchie MG, Nachman RJ, Dow JA, Denlinger DL, Davies S-A. Insect capa neuropeptides impact desiccation and cold tolerance. Proc Natl Acad Sci U S A. 2015;112(9):2882–7.
Wolfenden R, Snider M, Ridgway C, Miller B. The temperature dependence of enzyme rate enhancements. J Am Chem Soc. 1999;121(32):7419–20.
MacMillan HA, Sinclair BJ. Mechanisms underlying insect chill-coma. J Insect Physiol. 2011;57(1):12–20.
Kayukawa T, Chen B, Hoshizaki S, Ishikawa Y. Upregulation of a desaturase is associated with the enhancement of cold hardiness in the onion maggot, Delia antiqua. Insect Biochem Mol Biol. 2007;37(11):1160–7.
Storey KB, Storey JM. Biochemical strategies of overwintering in the gall gly larva, Eurosta solidaginis: Effect of low temperature acclimation on the activities of enzymes of intermediary metabolism. J Comp Physiol. 1981;144(2):191–9.
Fujiwara Y, Denlinger DL. p38 MAPK is a likely component of the signal transduction pathway triggering rapid cold hardening in the flesh fly Sarcophaga crassipalpis. J Exp Biol. 2007;210(18):3295–300.
Coello Alvarado LE, MacMillan HA, Sinclair BJ. Chill-tolerant Gryllus crickets maintain ion balance at low temperatures. J Insect Physiol. 2015;77:15–25.
O’Donnell MJ, Maddrell SH. Paracellular and transcellular routes for water and solute movements across insect epithelia. J Exp Biol. 1983;106(1):231–53.
O’Donnell MJ, Dow JAT, Huesmann GR, Tublitz NJ, Maddrell SH. Separate control of anion and cation transport in malpighian tubules of Drosophila melanogaster. J Exp Biol. 1996;199(5):1163–75.
Armitage WJ, Juss BK, Easty DL. Response of epithelial (MDCK) cell junctions to calcium removal and osmotic stress is influenced by temperature. Cryobiology. 1994;31(5):453–60.
Behrens J, Vakaet L, Friis R, Winterhager E, Van Roy F, Mareel MM, Birchmeier W. Loss of epithelial differentiation and gain of invasiveness correlates with tyrosine phosphorylation of the E-cadherin/beta-catenin complex in cells transformed with a temperature-sensitive v-SRC gene. J Cell Biol. 1993;120(3):757–66.
Gonzalez-Mariscal L, De Ramirez BC, Cereijido M. Effect of temperature on the occluding junctions of monolayers of epithelioid cells (MDCK). J Membr Biol. 1984;79(2):175–84.
Turner JR, Rill BK, Carlson SL, Carnes D, Kerner R, Mrsny RJ, Madara JL. Physiological regulation of epithelial tight junctions is associated with myosin light-chain phosphorylation. Am J Physiol. 1997;273(4):C1378–85.
Storey KB, Storey JM. Insect cold hardiness: metabolic, gene, and protein adaptation. Can J Zool. 2012;90(4):456–75.
Dunning LT, Dennis AB, Sinclair BJ, Newcomb RD, Buckley TR. Divergent transcriptional responses to low temperature among populations of alpine and lowland species of New Zealand stick insects (Micrarchus). Mol Ecol. 2014;23(11):2712–26.
Dennis AB, Dunning LT, Sinclair BJ, Buckley TR. Parallel molecular routes to cold adaptation in eight genera of New Zealand stick insects. Sci Rep. 2015;5:13965.
Sinclair BJ, Gibbs AG, Roberts SP. Gene transcription during exposure to, and recovery from, cold and desiccation stress in Drosophila melanogaster. Insect Mol Biol. 2007;16(4):435–43.
Colinet H, Lee SF, Hoffmann AA. Functional characterization of the Frost gene in Drosophila melanogaster: importance for recovery from chill coma. PLoS One. 2010;5(6):e10925.
Hoffmann AA, Blacket MJ, McKechnie SW, Rako L, Schiffer M, Rane RV, Good RT, Robin C, Lee SF. A proline repeat polymorphism of the Frost gene of Drosophila melanogaster showing clinal variation but not associated with cold resistance. Insect Mol Biol. 2012;21(4):437–45.
Qin W, Neal SJ, Robertson MR, Westwood JT, Walker VK. Cold hardening and transcriptional change in Drosophila melanogaster. Insect Mol Biol. 2005;14(6):607–13.
Bing X, Zhang J, Sinclair BJ. A comparison of Frost expression among species and life stages of Drosophila. Insect Mol Biol. 2012;21(1):31–9.
Udaka H, Percival‐Smith AJ, Sinclair BJ. Increased abundance of Frost mRNA during recovery from cold stress is not essential for cold tolerance in adult Drosophila melanogaster. Insect Mol Biol. 2013;22(5):541–50.
Teets NM, Peyton JT, Ragland GJ, Colinet H, Renault D, Hahn DA, Denlinger DL. Combined transcriptomic and metabolomic approach uncovers molecular mechanisms of cold tolerance in a temperate flesh fly. Physiol Genomics. 2012;44(15):764–77.
Gerken AR, Eller OC, Hahn DA, Morgan TJ. Constraints, independence, and evolution of thermal plasticity: probing genetic architecture of long- and short-term thermal acclimation. Proc Natl Acad Sci U S A. 2015;112(14):4399–404.
Zhang J, Marshall KE, Westwood JT, Clark MS, Sinclair BJ. Divergent transcriptomic responses to repeated and single cold exposures in Drosophila melanogaster. J Exp Biol. 2011;214(23):4021–9.
MacMillan HA, Knee JM, Dennis AB, Udaka H, Marshall KE, Merritt TJ, Sinclair BJ. Cold acclimation wholly reorganizes the Drosophila melanogaster transcriptome and metabolome. Sci Rep. 2016;6:28999.
Criddle N. Field crickets in Manitoba. Can Entomol. 1925;57(04):79–84.
Carrière Y, Simons AM, Roff DA. The effect of the timing of post-diapause egg development on survival, growth, and body size in Gryllus pennsylvanicus. Oikos. 1996;75(3):463–70.
Des Marteaux LE, Sinclair BJ. Ion and water balance in Gryllus crickets during the first twelve hours of cold exposure. J Insect Physiol. 2016;89:19–27.
Lai C-Q, Parnell LD, Lyman RF, Ordovas JM, Mackay TF. Candidate genes affecting Drosophila life span identified by integrating microarray gene expression analysis and QTL mapping. Mech Ageing Dev. 2007;128(3):237–49.
Goecks J, Nekrutenko A, Taylor J. Galaxy: a comprehensive approach for supporting accessible, reproducible, and transparent computational research in the life sciences. Genome Biol. 2010;11(8):R86.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.
Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, Couger MB, Eccles D, Li B, Lieber M. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 2013;8(8):1494–512.
Simão FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015;31:3210–2.
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(18):3674–6.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Meth. 2012;9(4):357–9.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7(3):562–78.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
Risso D, Ngai J, Speed TP, Dudoit S. Normalization of RNA-seq data using factor analysis of control genes or samples. Nat Biotechnol. 2014;32(9):896–902.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M. KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007;35 Suppl 2:W182–5.
Luo W, Brouwer C. Pathview: an R/Bioconductor package for pathway-based data integration and visualization. Bioinformatics. 2013;29(14):1830–1.
Luo W, Friedman MS, Shedden K, Hankenson KD, Woolf PJ. GAGE: generally applicable gene set enrichment for pathway analysis. BMC Bioinformatics. 2009;10(1):1.
Tassone EE, Geib SM, Hall B, Fabrick JA, Brent CS, Hull JJ. De novo construction of an expanded transcriptome assembly for the western tarnished plant bug, Lygus hesperus. GigaScience. 2016;5(6):1–5.
Theissinger K, Falckenhayn C, Blande D, Toljamo A, Gutekunst J, Makkonen J, Jussila J, Lyko F, Schrimpf A, Schulz R: De Novo assembly and annotation of the freshwater crayfish Astacus astacus transcriptome. Mar Genomics 2016, In press.
Piek T, Njio KD. Morphology and electrochemistry of insect muscle fibre membrane. Adv Insect Physiol. 1979;14:185–250.
Colinet H, Overgaard J, Com E, Sørensen JG. Proteomic profiling of thermal acclimation in Drosophila melanogaster. Insect Biochem Mol Biol. 2013;43(4):352–65.
Petzel DH, Hagedorn HH, Beyenbach KW. Preliminary isolation of mosquito natriuretic factor. Am J Physiol. 1985;249(4):R379–86.
Phillips JE, Audsley N. Neuropeptide control of ion and fluid transport across locust hindgut. Am Zool. 1995;35(6):503–14.
Philip BN, Yi S-X, Elnitsky MA, Lee RE. Aquaporins play a role in desiccation and freeze tolerance in larvae of the goldenrod gall fly, Eurosta solidaginis. J Exp Biol. 2008;211(7):1114–9.
Goto SG, Philip BN, Teets NM, Kawarasaki Y, Lee RE, Denlinger DL. Functional characterization of an aquaporin in the Antarctic midge Belgica antarctica. J Insect Physiol. 2011;57(8):1106–14.
Philip BN, Lee RE. Changes in abundance of aquaporin-like proteins occurs concomitantly with seasonal acquisition of freeze tolerance in the goldenrod gall fly, Eurosta solidaginis. J Insect Physiol. 2010;56(7):679–85.
Spring JH, Robichaux SR, Hamlin JA. The role of aquaporins in excretion in insects. J Exp Biol. 2009;212(3):358–62.
Chintapalli VR, Wang J, Herzyk P, Davies SA, Dow JA. Data-mining the FlyAtlas online resource to identify core functional motifs across transporting epithelia. BMC Genomics. 2013;14(1):1.
Ianowski JP, O’Donnell MJ. Basolateral ion transport mechanisms during fluid secretion by Drosophila Malpighian tubules: Na+ recycling, Na+:K+:2Cl- cotransport and Cl- conductance. J Exp Biol. 2004;207(15):2599–609.
Maddrell SH, O’Donnell MJ. Insect Malpighian tubules: V-ATPase action in ion and fluid transport. J Exp Biol. 1992;172(1):417–29.
Wessing A, Zierold K, Bertram G. Carbonic anhydrase supports electrolyte transport in Drosophila Malpighian tubules. Evidence by X-ray microanalysis of cryosections. J Insect Physiol. 1997;43(1):17–28.
Ianowski JP, O’Donnell MJ. Electrochemical gradients for Na+, K+, Cl - and H+ across the apical membrane in Malpighian (renal) tubule cells of Rhodnius prolixus. J Exp Biol. 2006;209(10):1964–75.
Nicolson S. The ionic basis of fluid secretion in insect Malpighian tubules: advances in the last ten years. J Insect Physiol. 1993;39(6):451–8.
Hanrahan J, Phillips J. Electrogenic, K+-dependent chloride transport in locust hindgut. Phil Trans R Soc B. 1982;299(1097):585–95.
Chapman RF. Regulation: digestion, nutrition, excretion. In: Comprehensive Insect Physiology, Biochemistry, and Pharmacology. Edited by Kerkut GA, Gilbert L, Vol. 4. Oxford: Pergamon; 1985.
Zeuthen T, MacAulay N. Cotransport of water by Na+–K+–2Cl− cotransporters expressed in Xenopus oocytes: NKCC1 versus NKCC2. J Physiol. 2012;590(5):1139–54.
Boyer AS, Ayerinskas II, Vincent EB, McKinney LA, Weeks DL, Runyan RB. TGFβ2 and TGFβ3 have separate and sequential activities during epithelial–mesenchymal cell transformation in the embryonic heart. Dev Biol. 1999;208(2):530–45.
Ebnet K, Aurrand-Lions M, Kuhn A, Kiefer F, Butz S, Zander K, zu Brickwedde M-KM, Suzuki A, Imhof BA, Vestweber D. The junctional adhesion molecule (JAM) family members JAM-2 and JAM-3 associate with the cell polarity protein PAR-3: a possible role for JAMs in endothelial cell polarity. J Cell Sci. 2003;116(19):3879–91.
Broermann A, Winderlich M, Block H, Frye M, Rossaint J, Zarbock A, Cagna G, Linnepe R, Schulte D, Nottebaum AF, et al. Dissociation of VE-PTP from VE-cadherin is required for leukocyte extravasation and for VEGF-induced vascular permeability in vivo. J Exp Med. 2011;208(12):2393–401.
Švastová E, Žilka N, Zat’ovičová M, Gibadulinová A, Čiampor F, Pastorek J, Pastoreková S. Carbonic anhydrase IX reduces E-cadherin-mediated adhesion of MDCK cells via interaction with β-catenin. Exp Cell Res. 2003;290(2):332–45.
Gut MO, Parkkila S, Vernerová Z, Rohde E, Závada J, Höcker M, Pastorek J, Karttunen T, Gibadulinová A, Závadová Z. Gastric hyperplasia in mice with targeted disruption of the carbonic anhydrase gene Car9. Gastroenterology. 2002;123(6):1889–903.
Hilvo M, Rafajová M, Pastoreková S, Pastorek J, Parkkila S. Expression of carbonic anhydrase IX in mouse tissues. J Histochem Cytochem. 2004;52(10):1313–21.
Satir P, Gilula NB. The fine structure of membranes and intercellular communication in insects. Annu Rev Entomol. 1973;18(1):143–66.
Tepass U, Tanentzapf G, Ward R, Fehon R. Epithelial cell polarity and cell junctions in Drosophila. Annu Rev Genet. 2001;35(1):747–84.
Adam AP. Regulation of endothelial adherens junctions by tyrosine phosphorylation. Mediat Inflamm. 2015;2015:24.
Matter K, Balda MS. Signalling to and from tight junctions. Nat Rev Mol Cell Biol. 2003;4(3):225–37.
Lecuit T, Lenne P-F. Cell surface mechanics and the control of cell shape, tissue patterns and morphogenesis. Nat Rev Mol Cell Biol. 2007;8(8):633–44.
Gupta BL, Wall BJ, Oschman JL, Hall T. Direct microprobe evidence of local concentration gradients and recycling of electrolytes during fluid absorption in the rectal papillae of Calliphora. J Exp Biol. 1980;88(1):21–48.
Haraguchi Y, Shimizu T, Yamato M, Kikuchi A, Okano T. Electrical coupling of cardiomyocyte sheets occurs rapidly via functional gap junction formation. Biomaterials. 2006;27(27):4765–74.
Reuss L. Tight junction permeability to ions and water. In: Cereijido M, Anderson JM, editors. Tight Junctions, vol. 2. 2001. p. 61–88.
Hartsock A, Nelson WJ. Adherens and tight junctions: structure, function and connections to the actin cytoskeleton. Biochim Biophys Acta. 2008;1778(3):660–9.
Rosson D, O’Brien TG, Kampherstein JA, Szallasi Z, Bogi K, Blumberg PM, Mullin JM. Protein kinase C-α activity modulates transepithelial permeability and cell junctions in the LLC-PK1 epithelial cell line. J Biol Chem. 1997;272(23):14950–3.
Dokladny K, Moseley PL, Ma TY. Physiologically relevant increase in temperature causes an increase in intestinal epithelial tight junction permeability. Am J Physiol. 2006;290(2):G204–12.
Berridge M. Transporting epithelia. New York: Academic; 1972.
Hildebrand JD. Shroom regulates epithelial cell shape via the apical positioning of an actomyosin network. J Cell Sci. 2005;118(22):5191–203.
Hildebrand JD, Soriano P. Shroom, a PDZ domain–containing actin-binding protein, is required for neural tube morphogenesis in mice. Cell. 1999;99(5):485–97.
Yu SP, Canzoniero LM, Choi DW. Ion homeostasis and apoptosis. Curr Opin Cell Biol. 2001;13(4):405–11.
Heimlich G, Bortner CD, Cidlowski JA. Apoptosis and cell volume regulation - The importance of ions and ion channels. In: Lauf PK, Adragna NC, editors. Cell Volume and Signaling, vol. 559. New York: Springer; 2004. p. 189–203.
Rojas RR, Leopold RA. Chilling injury in the housefly: evidence for the role of oxidative stress between pupariation and emergence. Cryobiology. 1996;33(4):447–58.
Kim M, Denlinger DL. Decrease in expression of beta‐tubulin and microtubule abundance in flight muscles during diapause in adults of Culex pipiens. Insect Mol Biol. 2009;18(3):295–302.
Štětina T, Koštál V, Korbelová J. The role of inducible Hsp70, and other heat shock proteins, in adaptive complex of cold tolerance of the fruit fly (Drosophila melanogaster). PLoS One. 2015;10(6):e0128976.
Elmore S. Apoptosis: a review of programmed cell death. Toxicol Pathol. 2007;35(4):495–516.
Lee RE, Denlinger DL. Stress tolerance in a polyextremophile: the southernmost insect. Can J Zool. 2014;93(9):679–86.
Teets NM, Denlinger DL. Autophagy in Antarctica: combating dehydration stress in the world’s southernmost insect. Autophagy. 2013;9(4):629–31.
Madara JL, Barenberg D, Carlson S. Effects of cytochalasin D on occluding junctions of intestinal absorptive cells: further evidence that the cytoskeleton may influence paracellular permeability and junctional charge selectivity. J Cell Biol. 1986;102(6):2125–36.
Drees F, Pokutta S, Yamada S, Nelson WJ, Weis WI. α-catenin is a molecular switch that binds E-cadherin-β-catenin and regulates actin-filament assembly. Cell. 2005;123(5):903–15.
Cottam DM, Tucker JB, Rogers‐Bald MM, Mackie JB, Macintyre J, Scarborough JA, Ohkura H, Milner MJ. Non-centrosomal microtubule-organising centres in cold-treated cultured Drosophila cells. Cell Motil Cytoskelet. 2006;63(2):88–100.
Huot J, Houle F, Spitz DR, Landry J. HSP27 phosphorylation-mediated resistance against actin fragmentation and cell death induced by oxidative stress. Cancer Res. 1996;56(2):273–9.
Joanisse D, Storey K. Oxidative stress and antioxidants in overwintering larvae of cold-hardy goldenrod gall insects. J Exp Biol. 1996;199(7):1483–91.
Joanisse DR, Storey KB. Oxidative stress and antioxidants in stress and recovery of cold-hardy insects. Insect Biochem Mol Biol. 1998;28(1):23–30.
Grubor‐Lajsic G, Block W, Telesmanic M, Jovanovic A, Stevanovic D, Baca F. Effect of cold acclimation on the antioxidant defense system of two larval Lepidoptera (Noctuidae). Arch Insect Biochem Physiol. 1997;36(1):1–10.
Niehaus AC, Wilson RS, Storm JJ, Angilletta Jr MJ. Fall field crickets did not acclimate to simulated seasonal changes in temperature. J Comp Physiol B. 2012;182(2):199–207.
Quinn NL, McGowan CR, Cooper GA, Koop BF, Davidson WS. Ribosomal genes and heat shock proteins as putative markers for chronic, sublethal heat stress in Arctic charr: applications for aquaculture and wild fish. Physiol Genomics. 2011;43(18):1056–64.
Pauli D, Tonka C-H, Ayme-Southgate A. An unusual split Drosophila heat shock gene expressed during embryogenesis, pupation and in testis. J Mol Biol. 1988;200(1):47–53.
Moskalev A, Shaposhnikov M, Turysheva E. Life span alteration after irradiation in Drosophila melanogaster strains with mutations of Hsf and Hsps. Biogerontology. 2008;10(1):3–11.
MacMillan HA, Guglielmo CG, Sinclair BJ. Membrane remodeling and glucose in Drosophila melanogaster: a test of rapid cold-hardening and chilling tolerance hypotheses. J Insect Physiol. 2009;55(3):243–9.
Lopez-Martinez G, Benoit JB, Rinehart JP, Elnitsky MA, Lee RE, Denlinger DL. Dehydration, rehydration, and overhydration alter patterns of gene expression in the Antarctic midge, Belgica antarctica. J Comp Physiol B. 2009;179(4):481–91.
Zerai DB, Fitzsimmons KM, Collier RJ. Transcriptional response of delta-9-desaturase gene to acute and chronic cold stress in Nile tilapia, Oreochromis niloticus. J World Aquacult Soc. 2010;41(5):800–6.
Tiku P, Gracey A, Macarteny A, Beynon R, Cossins A. Cold-induced expression of delta9-desaturase in carp by transcriptional and posttranslational mechanisms. Science. 1996;271(5250):815.
Andersen JL, MacMillan HA, Overgaard J. Muscle membrane potential and insect chill coma. J Exp Biol. 2015;218(16):2492–5.
MacMillan HA, Findsen A, Pedersen TH, Overgaard J. Cold-induced depolarization of insect muscle: differing roles of extracellular K+ during acute and chronic chilling. J Exp Biol. 2014;217(16):2930–8.
Bezzerides VJ, Ramsey IS, Kotecha S, Greka A, Clapham DE. Rapid vesicular translocation and insertion of TRP channels. Nat Cell Biol. 2004;6(8):709–20.
Goto SG, Denlinger DL. Short-day and long-day expression patterns of genes involved in the flesh fly clock mechanism: period, timeless, cycle and cryptochrome. J Insect Physiol. 2002;48(8):803–16.
Bownes M. Expression of the genes coding for vitellogenin (yolk protein). Annu Rev Entomol. 1986;31(1):507–31.
Thanks to Laura Ferguson for laboratory assistance, Alice Dennis and Morag Dick for bioinformatics guidance. Steve Xia, Tiffany Ng, Zayn Khamis, Jennifer Ho, and Kevin Lee assisted with insect rearing.
This research was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC) via Canada Graduate Scholarships to LED and JT, and a Discovery Grant to BJS, by and Ontario Graduate Scholarship to AHM, and a Japan Society for the Promotion of Science postdoctoral fellowship to HU. The funding bodies had no role in the design or implementation of the experiment nor in the collection, analysis, or interpretation of the data or the writing of the manuscript.
Availability of data and materials
The sequencing data sets supporting the results of this article are available via the NCBI Sequence Read Archive (accession number: SRP079970).
LED and BJS designed the study and drafted the manuscript, LED performed lab work, data analysis, LED, AHM, JT, and HU conceived and completed the bioinformatics pipeline. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
Summary of G. pennsylvanicus transcriptome de novo assembly. (DOCX 12 kb)
Additional file 2:
Spreadsheet S1. Complete list of genes differentially regulated in the hindgut and Malpighian tubules (MTs) following cold acclimation. Gene identities (Contig description) and GO term IDs were assigned by Blast2GO based on the top BLAST hit (GenBank accession number) of each contig. P-values were adjusted for false discovery rate (FDR). CPM, counts per million; FC, fold-change. (XLSX 356 kb)
Additional file 3:
Spreadsheet S2. List of GO terms that were differentially represented in the hindgut and Malpighian tubules (MTs) following cold acclimation. Read counts of contigs associated with each GO term were summed for cold- and warm-acclimated crickets. Positive Δread counts (cold-acclimated cricket read counts minus warm-acclimated cricket read counts) indicate upregulation of a GO term in cold-acclimated crickets, while negative Δread counts indicate downregulation in cold-acclimated crickets. (XLSX 82 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Des Marteaux, L.E., McKinnon, A.H., Udaka, H. et al. Effects of cold-acclimation on gene expression in Fall field cricket (Gryllus pennsylvanicus) ionoregulatory tissues. BMC Genomics 18, 357 (2017). https://doi.org/10.1186/s12864-017-3711-9
- Cold tolerance
- Ion pump
- Malpighian tubules
- Phenotypic plasticity