Open Access

Effects of cold-acclimation on gene expression in Fall field cricket (Gryllus pennsylvanicus) ionoregulatory tissues

  • Lauren E. Des Marteaux1Email author,
  • Alexander H McKinnon1, 2,
  • Hiroko Udaka1, 3,
  • Jantina Toxopeus1 and
  • Brent J. Sinclair1
BMC Genomics201718:357

https://doi.org/10.1186/s12864-017-3711-9

Received: 18 August 2016

Accepted: 20 April 2017

Published: 8 May 2017

Abstract

Background

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.

Results

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.

Conclusions

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.

Keywords

Insect Cold tolerance RNA-Seq Acclimation Ion pump Malpighian tubules Rectum Phenotypic plasticity Gryllus

Background

Most insects are chill-susceptible, such that their thermal performance and survival are limited in the cold at temperatures well above the freezing point [1]. 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 [24], disruption of the cytoskeleton [48], or induction of apoptosis [9]. Indirect chilling injury accumulates over time, most likely resulting from a loss of ion and water balance in the cold [1014] (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 [1722], although the underlying mechanisms are not completely understood.

Chill-susceptible insects lose hemolymph water and Na+ to the gut lumen during cold exposure [13, 2325] and re-establish this homeostasis during recovery upon rewarming [14]. Ion and water homeostasis in insects is primarily regulated at the Malpighian tubules and hindgut [26]. 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 [27] – is partially modified at the proximal tubule (a tight epithelium) prior to entering the gut lumen [2830]. Water and ions are then selectively reabsorbed from the gut lumen by the hindgut (particularly at the rectum) [31]. 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 [32]), and these peptides may be important for recovery from cold stress in Drosophila [33]. 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 [34], thus cold exposure should limit ion pumping rates while the rate of passive leak should remain relatively unchanged [35]. To compensate, insects could increase transport capacity by expressing or mobilizing more transport enzymes [2, 3638]. 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 [4245], 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 [4953], although its function remains elusive [54]. 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 [55], 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 [56]. 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. [47]). 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. [58] and Gerken et al. [56] 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 [58]) 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 [59]. The species is univoltine and overwinters in diapause in the soil as an egg [60]. 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 [13]. 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 [14]. 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 [39]. 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.

Methods

Insect rearing

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 [39]. 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 [62].

Cricket hindguts were dissected as described previously [61] 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 [63]. 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 [66]. 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 [67]. 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 [70]. We used Blast2GO [67] 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 [71] in R (v3.2.2, R Development Core Team, 2015 [72]). 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) [71]. 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 [73]). KEGG identities were assigned to contigs by the KEGG Automatic Annotation Server [74], 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) [76]. We accepted pathways as differentially-expressed if the FDR-adjusted P-value was < 0.1.

Results

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. [66]. 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.
Fig. 1

Number of genes upregulated, down-regulated, or unchanged in the hindgut and Malpighian tubules following cold acclimation. Differentially-expressed genes are those with an FDR alpha < 0.05 and a fold-change > 2. Note that due to some overlap in contigs the sum of genes up- or down-regulated across both tissues is less than the sum of genes up- or down-regulated in separate tissues (1439 and 1424 unique genes were up- and downregulated across both tissues, respectively)

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.

Hindgut

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.
Table 1

Selected genes upregulated in the hindgut following cold acclimation whose putative function in relation to cold tolerance is discussed in the text

Function

Description

Fold change

P-value

Species

Apoptosis

Caspase-6

4.0

1.8E-31

Gs

Caspase-8

4.0

9.1E-39

Zn

Circadian

Clock

2.7

1.7E-17

Gb

Nocturnin

2.4

1.8E-16

Zn

Period

7.4

1.6E-50

Gb

Cytoskeleton

Microtubule-Associated Protein Jupiter

4.2

5.3E-31

Zn

Diuresis

Atrial Natriuretic Peptide-Converting Enzyme

3.5

3.7E-15

Zn

Ion transport

Ca2+ Release-Activated Ca2+ Channel Protein 1

2.3

4.3E-17

Zn

Na+-K+ ATPase Alpha Subunit

2.8

2.7E-21

Ll

Neurotransmission

Na+- and Cl--Dependent GABA Transporter Ine

2.1

8.8E-10

Md

Phosphorylation

Dual Specificity Tyrosine-Phosphorylation-Regulated Kinase 2

2.7

2.2E-04

Zn

Serine Threonine-Protein Kinase Rio3

2.9

3.7E-07

Zn

Serine Threonine-Protein Kinase Sik3-Like Isoform X3

2.9

8.7E-16

Ap

Protein folding/chaperone

Heat Shock Protein 90

3.0

1.7E-28

Gf

Hsp70 Family Member

2.1

8.0E-14

Lm

Repair/detoxicant/antioxidant

Cytochrome P450

26.0

7.7E-108

Aa

Cytochrome P450 4C1

2.0

7.2E-05

Zn

DNA Mismatch Repair Protein Mlh1

3.5

1.2E-22

Xt

DNA Repair Protein Complementing Xp-G Cells

2.7

4.8E-20

Zn

Glutathione S-Transferase D7

2.1

2.6E-05

Zn

Signal transduction

cAMP-Specific 3’,5’-Cyclic Phosphodiesterase, Isoform F Isoform X2

2.1

3.3E-12

Tc

Dual 3’,5’ Cyclic-AMP and -GMP Phosphodiesterase 11

2.5

1.5E-17

Zn

G Kinase-Anchoring Protein 1

3.8

6.5E-34

Zn

G-Protein Coupled Receptor Mth2-Like

2.2

4.5E-14

Ap

Protein Kinase C Iota (partial)

2.4

1.5E-16

Zn

Signalling/gut contraction

Adenylate Cyclase Type (partial)

2.5

1.2E-13

Zn

For a full list of the 760 upregulated hindgut genes, see Additional file 2: Spreadsheet S1. P-values were adjusted for false discovery rate (FDR). For each gene, the species with the highest sequence similarity via BLAST is given. Species codes: Aa (Aedes aegypti), Ap (Acyrthosiphon pisum), Gb (Gryllus bimaculatus), Gf (Gryllus firmus), Gm (Galleria mellonella), Ll (Lutzomyia longipalpis), Lm (Locusta migratoria), Md (Microplitis demolitor), Ms (Modicogryllus siamensis), Tc (Tribolium castaneum), Xt (Xenopus tropicalis), Zn (Zootermopsis nevadensis)

Table 2

Selected genes downregulated in the hindgut following cold acclimation whose putative function in relation to cold tolerance is discussed in the text

Function

Description

Fold change

P-value

Species

Apoptosis

Apoptosis-Inducing Factor 3-Like

-4.4

6.1E-29

Nv

Circadian

Timeless

-4.3

2.3E-36

Gb

Ion transport

Bumetanide-Sensitive Na+-Cl- (partial) (NKCC)

-2.9

6.9E-25

Zn

Na+-Independent Sulfate Anion Transporter-Like

-2.8

1.9E-12

Mr

Organic Cation Transporter Protein

-2.6

4.2E-19

Zn

Ion transport/pH regulation

Carbonic Anhydrase 1

-4.6

4.9E-49

Tc

Carbonic Anhydrase 9

-3.9

6.5E-40

Zn

Mitochondrial Na+-H+ Exchanger NHA2

-2.9

7.0E-18

Zn

Phospholipid biochemistry

∆9 Desaturase 1

-3.9

1.4E-10

Ad

Protein folding/chaperone

Heat Shock Protein 67B2

-2.1

3.3E-09

Zn

Heat Shock Protein β-11

-4.5

5.1E-29

Zn

Repair/detoxicant/antioxidant

Antioxidant Enzyme

-2.0

1.6E-12

Go

Cytochrome P450 4C1

-4.1

4.0E-42

Zn

Cytochrome P450 6A14

-3.8

1.3E-12

Zn

Cytochrome P450 9E1

-2.5

9.5E-19

Dp

Cytochrome P450 Cyp44

-5.5

1.9E-43

Zn

Epsilon Glutathione S-Transferase

-2.1

7.8E-11

Lm

Glutathione S-Transferase

-6.8

1.2E-35

Bt

Glutathione S-Transferase-Like

-2.5

2.3E-15

Md

Glutathione S-Transferase Sigma 1

-2.5

2.1E-10

Sg

Glutathione S-Transferase Sigma 7

-4.8

1.5E-26

Lm

Glutathione S-Transferase Theta 1

-2.8

2.5E-21

Lm

Peroxiredoxin

-2.6

4.3E-17

Go

Vitellogenin

-40.0

1.7E-47

Ar

Vitellogenin-2

-18.4

1.2E-50

Ps

For a full list of the 733 downregulated hindgut genes, see Additional file 2: Spreadsheet S1. P-values were adjusted for false discovery rate (FDR). For each gene, the species with the highest sequence similarity via BLAST is given. Species codes: Ac (Acheta domesticus), Ap (Aphis gossypii), Ar (Athalia rosae), Bt (Bemisia tabaci), Dp (Diploptera punctata), Gb (Gryllus bimaculatus), Go (Gryllus orientalis), Lm (Locusta migratoria), Md (Microplitis demolitor), Mr (Megachile rotundata), Nv (Nasonia vitripennis), Ps (Plautia stali), Sg (Schistocerca gregaria), Tc (Tribolium castaneum), Zn (Zootermopsis nevadensis)

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.
Fig. 2

Heat map of differentially-expressed KEGG pathways in the hindgut and Malpighian tubules following cold acclimation. Upregulated pathways are given in orange and downregulated pathways are given in blue. Each heat map contains three column indicating three cold-acclimated biological replicates (CA 1-3) each compared to the mean expression among warm-acclimated replicates. For a complete description of each pathway, see the KEGG online resource (http://www.genome.jp/kegg/)

Fig. 3

Shifts in the expression of ‘adherens junction’ KEGG pathway components in the cricket hindgut following cold acclimation, as an example of a pathway that was significantly differentially regulated. Each pathway component contains three color bars indicating three cold-acclimated biological replicates each compared to the mean expression among warm-acclimated replicates. For cold-acclimated crickets relative to warm-acclimated crickets, shifts in expression are either upregulated (orange), downregulated (blue), or unchanged (grey). For a complete description of each pathway component, see the KEGG ‘adherens junction’ reference pathway (http://www.genome.jp/kegg-bin/show_pathway?ko04520)

Fig. 4

Shifts in expression of the ‘cardiac muscle contraction’ KEGG pathway components in the cricket hindgut following cold acclimation, as an example of a pathway that was significantly differentially regulated. Each pathway component contains three color bars indicating three cold-acclimated biological replicates each compared to the mean expression among warm-acclimated replicates. For cold-acclimated crickets relative to warm-acclimated crickets, shifts in expression are either upregulated (orange), downregulated (blue), or unchanged (grey). ATP - Na+-K+ ATPase α subunit, Cyto - cytochrome c reductase iron-sulfur subunit, TPM - tropomyosin 1, Myosin - myosin heavy chain 6/7. For a complete description of each pathway component, see the KEGG ‘cardiac muscle contraction’ reference pathway (http://www.genome.jp/kegg-bin/show_pathway?ko04260)

Malpighian tubules

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.
Table 3

Selected upregulated genes in the Malpighian tubules following cold acclimation whose putative function in relation to cold tolerance is discussed in the text

Function

Description

Fold change

P-value

Species

Apoptosis

Caspase-8

2.1

1.7E-5

Zn

Programmed Cell Death Protein 2

2.4

7.0E-10

Nv

Circadian

Clock

2.5

6.4E-9

Gb

Nocturnin

3.5

1.2E-23

Zn

Period

7.6

2.5E-63

Ms

Cytoskeleton

Gamma-Tubulin Complex Component 6

2.2

8.6E-10

Zn

Kinesin-Like Protein Costa

2.8

2.6E-9

Zn

Protein Shroom

2.5

1.6E-9

Zn

Microtubule-Associated Protein Jupiter

2.2

1.8E-8

Zn

Ion transport

Ca2+-Transporting ATPase Type 2C Member 1

2.2

1.1E-8

Zn

Solute Carrier Organic Anion Transporter Family Member 5a1

2.4

3.0E-11

Hs

Phosphorylation

Inositol Polyphosphate Multikinase

2.2

2.6E-15

Zn

Serine Threonine-Protein Kinase rio3

3.4

5.4E-35

Zn

Serine Threonine-Protein Kinase pctaire-2

2.6

3.9E-14

Zn

Serine Threonine-Protein Kinase Tousled-Like 2

2.9

2.7E-11

Zn

Tyrosine-Protein Kinase Transmembrane Receptor ror1

2.5

6.6E-10

Phc

Tyrosine-Protein Phosphatase Non-Receptor Type 23

2.5

1.6E-15

Zn

Protein folding/chaperone

Hsp 70 Family Member

2.0

1.0E-10

Lm

Hsp 70-binding Protein 1

3.3

2.7E-22

Zn

Repair/detoxicant/antioxidant

Alkylated DNA Repair Protein Alkb-Like Protein 4

2.6

1.4E-12

Zn

Cytochrome P450

14.6

3.7E-58

Aa

Cytochrome P450 2j2

3.4

4.6E-7

Zn

Cytochrome P450 partial

6.1

9.1E-11

Zn

DNA Mismatch Repair Protein mlh1

4.2

3.9E-24

Xt

Glutathione S-Transferase Sigma 7

2.3

9.5E-9

Zn

Peroxiredoxin-6

2.1

4.0E-9

Lm

Probable Cytochrome P450 (partial)

6.2

2.6E-13

Tc

Signal transduction

cAMP-Specific 3’,5’-Cyclic Isoform F Isoform X2

3.2

2.4E-30

Tc

G Kinase-Anchoring Protein 1

3.7

1.2E-28

Zn

For a full list of the 999 upregulated Malpighian tubule genes, see Additional file 2: Spreadsheet S1. P-values were adjusted for false discovery rate (FDR). For each gene, the species with the highest sequence similarity via BLAST is given. Species codes: Aa (Aedes aegypti), Gb (Gryllus bimaculatus), Hs (Harpegnathos saltator), Lm (Locusta migratoria), Ms (Modicogryllus siamensis), Nv (Nasonia vitripennis), Phc (Pediculus humanus corporis), Tc (Tribolium castaneum), Xt (Xenopus tropicalis), Zn (Zootermopsis nevadensis)

Table 4

Selected downregulated genes in the Malpighian tubules following cold acclimation whose putative function in relation to cold tolerance is discussed in the text

Function

Description

Fold change

P-value

Species

Autophagy

Autophagy-Related Protein 2-Like Protein B

-2.1

1.5E-7

Zn

Circadian

Timeless

-3.0

3.3E-13

Gb

Cytoskeleton

Microtubule-Associated Proteins 1A 1B Light Chain 3C

-4.9

8.5E-23

Zn

Epidermal Growth Factor Receptor Kinase Substrate 8-Like Isoform X1

-2.4

1.9E-9

Md

Gamma-Tubulin Complex Component 3

-2.0

4.7E-9

Zn

Development

Juvenile Hormone-Inducible

-10.6

3.8E-50

Aa

Ion transport/pH regulation

Carbonic Anhydrase 9

-4.2

2.1E-33

Zn

Carbonic Anhydrase 1

-2.2

5.9E-15

Tc

V-ATPase Subunit D

-2.1

1.2E-14

Lm

Phosphorylation

Inositol Monophosphatase

-6.9

1.1E-58

Zn

Serine Threonine-Protein Phosphatase 2B Catalytic Subunit 2-Like Isoform X2

-2.1

9.0E-14

Tc

Protein folding/chaperone

Heat Shock Protein 67B2

-2.8

1.6E-16

Zn

Repair/detoxicant/antioxidant

Antioxidant Enzyme

-2.9

3.7E-23

Go

Cytochrome P450

-2.9

1.4E-11

Bt

Cytochrome P450 4C1

-2.7

3.2E-22

Zn

Cytochrome P450 6A14

-2.5

3.2E-13

Zn

Cytochrome P450 9E2

-2.3

1.4E-11

Zn

Cytochrome P450 (partial)

-2.3

3.6E-12

Zn

Cytochrome P450 9E1

-2.2

3.3E-14

Dp

Glutathione S-Transferase

-3.2

1.8E-23

Bt

Glutathione S-Transferase Theta 1

-3.1

1.9E-28

Lm

For a full list of the 1009 downregulated Malpighian tubule genes, see Additional file 2: Spreadsheet S1. P-values were adjusted for false discovery rate (FDR). For each gene, the species with the highest sequence similarity via BLAST is given. Species codes: Aa (Aedes aegypti), Bt (Bemisia tabaci), Dp (Diplotera punctata), Gb (Gryllus bimaculatus), Go (Gryllotalpa orientalis), Lm (Locusta migratoria), Md (Microplitis demolitor), Tc (Tribolium castaneum), Zn (Zootermopsis nevadensis)

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, [79]) 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).

Discussion

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.

Water balance

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 [82] and help to defend hemolymph volume during cold exposure. Although some aquaporins have been associated with insect freeze tolerance [8385], their role in cold acclimation among chill-susceptible insects is unknown, and none of the water-transporting insect aquaporins [86] were differentially-expressed in the hindgut or Malpighian tubules with cold acclimation.

Ion transport

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 [87]. 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 [88], 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 [22], suggesting a decoupling of Na+ gradients across the gut rather than compensation for slowed NKA activity in the cold [22]. 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 [9698]), 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 [13]. In addition to possible roles in providing protons for ion transport, CA9 may have roles in cellular adhesion [99], proliferation, and differentiation (at least in mammals [100]). 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 [101], 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, 102104], and are closely related [105]. 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 [97] and downregulation of PKC-α [111] 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.

Chilling injury

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, 116120]. 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 [121]), 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 [56], 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 [122]. Cold-acclimated crickets had altered expression of cytoskeletal branching and stabilizing genes, such as α-catenin, ARP2/3 (a nucleation site for actin polymerization [125]), 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 [8], Delia antiqua [7], alfalfa (Medicago sativa) [4], 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, 128130]. 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 [80]. Genes encoding several hsps were downregulated, including the less-characterized hsp β11 (involved in the vertebrate heat stress response [132]) 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 [141]), 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 [142]). 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 [80]; in this case, we suspect that vitellogenin mRNA came from traces of fat body remaining on the gut during dissections [143].

Conclusions

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.
Fig. 5

Candidate mechanisms of cold acclimation in Gryllus pennsylvanicus. a Increased expression of NKA in rectal pad epithelia should increase Na+ and water reabsorption; this may counteract leak of water and ions and aid in chill coma recovery). b Downregulation of CA and V-ATPase expression in the Malpighian tubules should slow primary urine production, thereby retaining hemolymph volume. c Cytoskeletal and junctional remodelling of the hindgut may mitigate water and ion leak during cold exposure

Abbreviations

AP2: 

Adaptor-related protein complex 2 subunit α

ARP2/3: 

Actin-related proteins 2 and 3

bp: 

Base pair

CA1: 

Carbonic anhydrase isoform 1

CA9: 

Carbonic anhydrase isoform 9

CTmin

Critical thermal minimum

EGFR/ErbB-1: 

Epidermal growth factor receptor

e-value: 

Expectation value

hsp: 

Heat shock protein

MAP1A/1B: 

Microtubule-associated proteins 1A/1B

NHA: 

Na+-H+ exchanger

NKA: 

Na+-K+ ATPase

NKCC: 

Bumetanide-sensitive Na+-K+-2Cl- cotransporter

NSF: 

N-ethylmaleimide-sensitive factor

PAR3: 

Partitioning defective 3

PKC-i: 

Protein kinase C iota

PKC-α: 

Protein kinase C α

ROS: 

Reactive oxygen species

RTK: 

Receptor tyrosine kinase

TGFβ2: 

Transforming growth factor β2

TUBA: 

Tubulin α

V-ATPase: 

Vacuolar-type H+ ATPase

VE-PTP: 

Vascular endothelial protein tyrosine phosphatase

Declarations

Acknowledgements

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.

Funding

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).

Authors’ contributions

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.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open AccessThis 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.

Authors’ Affiliations

(1)
Department of Biology, The University of Western Ontario
(2)
Present Address: Faculty of Medicine, University of Manitoba
(3)
Present Address: Graduate School of Science, Biological Sciences, Kyoto University

References

  1. Bale JS. Classes of insect cold hardiness. Funct Ecol. 1993;7(6):751–3.Google Scholar
  2. Clark MS, Worland MR. How insects survive the cold: molecular mechanisms—a review. J Comp Physiol B. 2008;178(8):917–33.PubMedView ArticleGoogle Scholar
  3. 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.PubMedView ArticleGoogle Scholar
  4. Ö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.PubMedView ArticleGoogle Scholar
  5. Michaud MR, Denlinger DL. Molecular modalities of insect cold survival: current understanding and future trends. Int Congr Ser. 2004;1275:32–46.View ArticleGoogle Scholar
  6. 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.Google Scholar
  7. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  8. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  9. Yi S-X, Moore CW, Lee Jr RE. Rapid cold-hardening protects Drosophila melanogaster from cold-induced apoptosis. Apoptosis. 2007;12(7):1183–93.PubMedView ArticleGoogle Scholar
  10. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  11. 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.Google Scholar
  12. MacMillan HA, Baatrup E, Overgaard J. Concurrent effects of cold and hyperkalaemia cause insect chilling injury. Proc R Soc B. 2015;282(1817):20151483.PubMedPubMed CentralView ArticleGoogle Scholar
  13. 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.PubMedView ArticleGoogle Scholar
  14. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  15. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  16. 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.View ArticleGoogle Scholar
  17. 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.PubMedView ArticleGoogle Scholar
  18. 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.View ArticleGoogle Scholar
  19. 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.View ArticleGoogle Scholar
  20. 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.Google Scholar
  21. 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.PubMedView ArticleGoogle Scholar
  22. 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.Google Scholar
  23. 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.PubMedView ArticleGoogle Scholar
  24. 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.View ArticleGoogle Scholar
  25. 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.PubMedView ArticleGoogle Scholar
  26. Phillips J. Comparative physiology of insect renal function. Am J Physiol. 1981;241(5):R241–57.PubMedGoogle Scholar
  27. 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.Google Scholar
  28. Beyenbach KW. Transport mechanisms of diuresis in Malpighian tubules of insects. J Exp Biol. 2003;206(21):3845–56.PubMedView ArticleGoogle Scholar
  29. Beyenbach KW. Mechanism and regulation of electrolyte transport in Malpighian tubules. J Insect Physiol. 1995;41(3):197–207.View ArticleGoogle Scholar
  30. Coast GM. The regulation of primary urine production in insects. In: Society for Experimental Biology: 1998. Cambridge: Cambridge University Press; 1998. pp. 189–209.Google Scholar
  31. Phillips JE, Hanrahan J, Chamberlin M, Thomson B. Mechanisms and control of reabsorption in insect hindgut. Adv Insect Physiol. 1987;19:329–422.View ArticleGoogle Scholar
  32. 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.Google Scholar
  33. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  34. Wolfenden R, Snider M, Ridgway C, Miller B. The temperature dependence of enzyme rate enhancements. J Am Chem Soc. 1999;121(32):7419–20.View ArticleGoogle Scholar
  35. MacMillan HA, Sinclair BJ. Mechanisms underlying insect chill-coma. J Insect Physiol. 2011;57(1):12–20.PubMedView ArticleGoogle Scholar
  36. 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.PubMedView ArticleGoogle Scholar
  37. 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.View ArticleGoogle Scholar
  38. 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.PubMedView ArticleGoogle Scholar
  39. Coello Alvarado LE, MacMillan HA, Sinclair BJ. Chill-tolerant Gryllus crickets maintain ion balance at low temperatures. J Insect Physiol. 2015;77:15–25.PubMedView ArticleGoogle Scholar
  40. 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.PubMedGoogle Scholar
  41. 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.PubMedGoogle Scholar
  42. 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.PubMedView ArticleGoogle Scholar
  43. 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.PubMedView ArticleGoogle Scholar
  44. 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.PubMedView ArticleGoogle Scholar
  45. 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.PubMedGoogle Scholar
  46. Storey KB, Storey JM. Insect cold hardiness: metabolic, gene, and protein adaptation. Can J Zool. 2012;90(4):456–75.Google Scholar
  47. 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.PubMedView ArticleGoogle Scholar
  48. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  49. 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.Google Scholar
  50. 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.Google Scholar
  51. 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.Google Scholar
  52. 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.Google Scholar
  53. 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.Google Scholar
  54. 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.Google Scholar
  55. 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.PubMedView ArticleGoogle Scholar
  56. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  57. 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.PubMedView ArticleGoogle Scholar
  58. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  59. Criddle N. Field crickets in Manitoba. Can Entomol. 1925;57(04):79–84.View ArticleGoogle Scholar
  60. 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.View ArticleGoogle Scholar
  61. 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.PubMedView ArticleGoogle Scholar
  62. 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.PubMedView ArticleGoogle Scholar
  63. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  64. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  65. 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.PubMedView ArticleGoogle Scholar
  66. 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.PubMedView ArticleGoogle Scholar
  67. 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.PubMedView ArticleGoogle Scholar
  68. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  69. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Meth. 2012;9(4):357–9.View ArticleGoogle Scholar
  70. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  71. 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.PubMedView ArticleGoogle Scholar
  72. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  73. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.PubMedPubMed CentralView ArticleGoogle Scholar
  74. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  75. Luo W, Brouwer C. Pathview: an R/Bioconductor package for pathway-based data integration and visualization. Bioinformatics. 2013;29(14):1830–1.PubMedPubMed CentralView ArticleGoogle Scholar
  76. 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.View ArticleGoogle Scholar
  77. 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.Google Scholar
  78. 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.Google Scholar
  79. Piek T, Njio KD. Morphology and electrochemistry of insect muscle fibre membrane. Adv Insect Physiol. 1979;14:185–250.View ArticleGoogle Scholar
  80. 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.PubMedView ArticleGoogle Scholar
  81. Petzel DH, Hagedorn HH, Beyenbach KW. Preliminary isolation of mosquito natriuretic factor. Am J Physiol. 1985;249(4):R379–86.PubMedGoogle Scholar
  82. Phillips JE, Audsley N. Neuropeptide control of ion and fluid transport across locust hindgut. Am Zool. 1995;35(6):503–14.View ArticleGoogle Scholar
  83. 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.PubMedView ArticleGoogle Scholar
  84. 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.PubMedView ArticleGoogle Scholar
  85. 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.PubMedView ArticleGoogle Scholar
  86. Spring JH, Robichaux SR, Hamlin JA. The role of aquaporins in excretion in insects. J Exp Biol. 2009;212(3):358–62.PubMedView ArticleGoogle Scholar
  87. 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.View ArticleGoogle Scholar
  88. 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.PubMedView ArticleGoogle Scholar
  89. Maddrell SH, O’Donnell MJ. Insect Malpighian tubules: V-ATPase action in ion and fluid transport. J Exp Biol. 1992;172(1):417–29.PubMedGoogle Scholar
  90. 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.PubMedView ArticleGoogle Scholar
  91. 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.PubMedView ArticleGoogle Scholar
  92. 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.View ArticleGoogle Scholar
  93. Hanrahan J, Phillips J. Electrogenic, K+-dependent chloride transport in locust hindgut. Phil Trans R Soc B. 1982;299(1097):585–95.PubMedView ArticleGoogle Scholar
  94. Chapman RF. Regulation: digestion, nutrition, excretion. In: Comprehensive Insect Physiology, Biochemistry, and Pharmacology. Edited by Kerkut GA, Gilbert L, Vol. 4. Oxford: Pergamon; 1985.Google Scholar
  95. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  96. 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.PubMedView ArticleGoogle Scholar
  97. 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.PubMedView ArticleGoogle Scholar
  98. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  99. Š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.PubMedView ArticleGoogle Scholar
  100. 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.PubMedView ArticleGoogle Scholar
  101. 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.PubMedView ArticleGoogle Scholar
  102. Satir P, Gilula NB. The fine structure of membranes and intercellular communication in insects. Annu Rev Entomol. 1973;18(1):143–66.PubMedView ArticleGoogle Scholar
  103. Tepass U, Tanentzapf G, Ward R, Fehon R. Epithelial cell polarity and cell junctions in Drosophila. Annu Rev Genet. 2001;35(1):747–84.PubMedView ArticleGoogle Scholar
  104. Adam AP. Regulation of endothelial adherens junctions by tyrosine phosphorylation. Mediat Inflamm. 2015;2015:24.View ArticleGoogle Scholar
  105. Matter K, Balda MS. Signalling to and from tight junctions. Nat Rev Mol Cell Biol. 2003;4(3):225–37.PubMedView ArticleGoogle Scholar
  106. 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.PubMedView ArticleGoogle Scholar
  107. 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.Google Scholar
  108. 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.PubMedView ArticleGoogle Scholar
  109. Reuss L. Tight junction permeability to ions and water. In: Cereijido M, Anderson JM, editors. Tight Junctions, vol. 2. 2001. p. 61–88.Google Scholar
  110. Hartsock A, Nelson WJ. Adherens and tight junctions: structure, function and connections to the actin cytoskeleton. Biochim Biophys Acta. 2008;1778(3):660–9.PubMedView ArticleGoogle Scholar
  111. 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.PubMedView ArticleGoogle Scholar
  112. 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.Google Scholar
  113. Berridge M. Transporting epithelia. New York: Academic; 1972.View ArticleGoogle Scholar
  114. Hildebrand JD. Shroom regulates epithelial cell shape via the apical positioning of an actomyosin network. J Cell Sci. 2005;118(22):5191–203.PubMedView ArticleGoogle Scholar
  115. 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.PubMedView ArticleGoogle Scholar
  116. Yu SP, Canzoniero LM, Choi DW. Ion homeostasis and apoptosis. Curr Opin Cell Biol. 2001;13(4):405–11.PubMedView ArticleGoogle Scholar
  117. 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.View ArticleGoogle Scholar
  118. 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.View ArticleGoogle Scholar
  119. 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.Google Scholar
  120. Š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.PubMedPubMed CentralView ArticleGoogle Scholar
  121. Elmore S. Apoptosis: a review of programmed cell death. Toxicol Pathol. 2007;35(4):495–516.PubMedPubMed CentralView ArticleGoogle Scholar
  122. Lee RE, Denlinger DL. Stress tolerance in a polyextremophile: the southernmost insect. Can J Zool. 2014;93(9):679–86.Google Scholar
  123. Teets NM, Denlinger DL. Autophagy in Antarctica: combating dehydration stress in the world’s southernmost insect. Autophagy. 2013;9(4):629–31.PubMedPubMed CentralView ArticleGoogle Scholar
  124. 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.PubMedView ArticleGoogle Scholar
  125. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  126. 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.View ArticleGoogle Scholar
  127. 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.PubMedGoogle Scholar
  128. 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.PubMedGoogle Scholar
  129. 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.View ArticleGoogle Scholar
  130. 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.View ArticleGoogle Scholar
  131. 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.PubMedView ArticleGoogle Scholar
  132. 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.PubMedView ArticleGoogle Scholar
  133. 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.PubMedView ArticleGoogle Scholar
  134. 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.PubMedView ArticleGoogle Scholar
  135. 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.PubMedView ArticleGoogle Scholar
  136. 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.PubMedView ArticleGoogle Scholar
  137. 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.View ArticleGoogle Scholar
  138. 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.PubMedView ArticleGoogle Scholar
  139. Andersen JL, MacMillan HA, Overgaard J. Muscle membrane potential and insect chill coma. J Exp Biol. 2015;218(16):2492–5.PubMedView ArticleGoogle Scholar
  140. 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.PubMedView ArticleGoogle Scholar
  141. 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.PubMedView ArticleGoogle Scholar
  142. 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.PubMedView ArticleGoogle Scholar
  143. Bownes M. Expression of the genes coding for vitellogenin (yolk protein). Annu Rev Entomol. 1986;31(1):507–31.View ArticleGoogle Scholar

Copyright

© The Author(s). 2017