Proteomics of intracellular freezing survival

Panagrolaimus sp. DAW1, a nematode cultured from the Antarctic, has the extraordinary physiological ability to survive total intracellular freezing throughout all of its compartments. While a few other organisms, all nematodes, have subsequently also been found to survive freezing in this manner, P. sp. DAW1 has so far shown the highest survival rates. In addition, P. sp. DAW1 is also, depending on the rate or extent of freezing, able to undergo cryoprotective dehydration. In this study, the proteome of P. sp DAW1 is explored, highlighting a number of differentially expressed proteins and pathways that occur when the nematodes undergo intracellular freezing. Among the strongest signals after being frozen is an upregulation of proteases and the downregulation of cytoskeletal and antioxidant activity, the latter possibly accumulated before freezing much in the way the sugar trehalose has been shown to be stored during acclimation.

Genistein and estradiol have common and specific impacts on the sea bass (Dicentrarchus labrax) skin-scale barrier

Teleost fish scales play important roles in animal protection and homeostasis. They can be targeted by endogenous estrogens and by environmental estrogenic endocrine disruptors. The phytoestrogen genistein is ubiquitous in the environment and in aquaculture feeds and is a disruptor of estrogenic processes in vertebrates. To test genistein disrupting actions in teleost fish we used a minimally invasive approach by analysing scales plucked from the skin of sea bass (Dicentrarchus labrax). Genistein transactivated all three fish nuclear estrogen receptors and was most potent with the Esr2, had the highest efficacy with Esr1, but reached, in all cases, transactivation levels lower than those of estradiol. RNA-seq revealed 254 responsive genes in the sea bass scales transcriptome with an FDR < 0.05 and more than 2-fold change in expression, 1 or 5 days after acute exposure to estradiol or to genistein. 65 genes were specifically responsive to estradiol and 106 by genistein while 83 genes were responsive to both compounds. Estradiol specifically regulated genes of protein/matrix turnover and genistein affected sterol biosynthesis and regeneration, while innate immune responses were affected by both compounds. This comprehensive study revealed the impact on the fish scale transcriptome of estradiol and genistein, providing a solid background to further develop fish scales as a practical screening tool for endocrine disrupting chemicals in teleosts.

Signatures of the Evolution of Parthenogenesis and Cryptobiosis in the Genomes of Panagrolaimid Nematodes

Most animal species reproduce sexually and fully parthenogenetic lineages are usually short lived in evolution. Still, parthenogenesis may be advantageous as it avoids the cost of sex and permits colonization by single individuals. Panagrolaimid nematodes have colonized environments ranging from arid deserts to Arctic and Antarctic biomes. Many are obligatory meiotic parthenogens, and most have cryptobiotic abilities, being able to survive repeated cycles of complete desiccation and freezing. To identify systems that may contribute to these striking abilities, we sequenced and compared the genomes and transcriptomes of parthenogenetic and outcrossing panagrolaimid species, including cryptobionts and non-cryptobionts. The parthenogens are triploids, most likely originating through hybridization. Adaptation to cryptobiosis shaped the genomes of panagrolaimid nematodes and is associated with the expansion of gene families and signatures of selection on genes involved in cryptobiosis. All panagrolaimids have acquired genes through horizontal gene transfer, some of which are likely to contribute to cryptobiosis.

Estradiol and genistein effects on the sea bass (Dicentrarchus labrax) scales: Transcriptome dataset

Fish scales are mineralized structures that play important roles in protection and mineral homeostasis. This tissue expresses multiple estrogen receptor subtypes and can be targeted by estrogens or estrogenic endocrine-disrupting compounds, but their effects are poorly explored. The transcriptome data here presented support the findings reported in the research article “Genistein and estradiol have common and specific impacts on the sea bass (Dicentrarchus labrax) skin-scale barrier” [1]. Juvenile sea bass were exposed to estradiol and the phytoestrogen genistein for 1 and 5 days, by intraperitoneal injections, and the effects on scale transcript expression were analysed by RNA-seq using an Illumina Hi-seq 1500. The raw reads of the 30 libraries produced have been deposited in the NCBI-SRA database with the project accession number SRP102504. Mapping of RNA-seq reads against the sea bass reference genome using the Cufflinks/TopHat package identified 371 genes that had significant (FDR<0.05) differential expression with the estradiol or genistein treatments in relation to the control scales at each exposure time, 254 of which presented more than a 2-fold change in expression. The identity of the differentially expressed genes was obtained using both automatic and manual annotations against multiple public sequence databases and they were grouped according to their patterns of expression using hierarchical clustering and heat-maps. The biological processes and KEGG pathways most significantly affected by the estradiol and/or genistein treatments were identified using Cytoscape/ClueGO enrichment analyses.

Molecular mechanisms underpinning transgenerational plasticity in the green sea urchin Psammechinus miliaris

The pre-conditioning of adult marine invertebrates to altered conditions, such as low pH, can significantly impact offspring outcomes, a process which is often referred to as transgenerational plasticity (TGP). This study describes for the first time, the gene expression profiles associated with TGP in the green sea urchin Psammechinus miliaris and evaluates the transcriptional contribution to larval resilience. RNA-Seq was used to determine how the expression profiles of larvae spawned into low pH from pre-acclimated adults differed to those of larvae produced from adults cultured under ambient pH. The main findings demonstrated that adult conditioning to low pH critically pre-loads the embryonic transcriptional pool with antioxidants to prepare the larvae for the “new” conditions. In addition, the classic cellular stress response, measured via the production of heat shock proteins (the heat shock response (HSR)), was separately evaluated. None of the early stage larvae either spawned in low pH (produced from both ambient and pre-acclimated adults) or subjected to a separate heat shock experiment were able to activate the full HSR as measured in adults, but the capacity to mount an HSR increased as development proceeded. This compromised ability clearly contributes to the vulnerability of early stage larvae to acute environmental challenge.

Investigating trehalose synthesis genes after cold acclimation in the Antarctic nematode Panagrolaimus sp. DAW1

Panagrolaimus sp. DAW1 is a freeze-tolerant Antarctic nematode which survives extensive intracellular ice formation. The molecular mechanisms of this extreme adaptation are still poorly understood. We recently showed that desiccation-enhanced RNA interference (RNAi) soaking can be used in conjunction with quantitative polymerase chain reaction (qPCR) to screen for phenotypes associated with reduced expression of candidate genes in Panagrolaimus sp. DAW1. Here, we present the use of this approach to investigate the role of trehalose synthesis genes in this remarkable organism. Previous studies have shown that acclimating Panagrolaimus sp. DAW1 at 5°C before freezing or desiccation substantially enhances survival. In this study, the expression of tps-2 and other genes associated with trehalose metabolism, as well as lea-1, hsp-70 and gpx-1, in cold-acclimated and non-acclimated nematodes was analyzed using qPCR. Pd-tps-2 and Pd-lea-1 were significantly upregulated after cold acclimation, indicating an inducible expression in the cold adaptation of Panagrolaimus sp. DAW1. The role of trehalose synthesis genes in Panagrolaimus sp. DAW1 was further investigated by RNAi. Compared to the controls, Pd-tps-2a(RNAi)-treated and cold-acclimated nematodes showed a significant decrease in mRNA, but no change in trehalose content or freezing survival. The involvement of two other trehalose synthesis genes (tps-2b and gob-1) was also investigated. These findings provide the first functional genomic investigation of trehalose synthesis genes in the non-model organism Panagrolaimus sp. DAW1. The presence of several trehalose synthesis genes with different RNAi sensitivities suggests the existence of multiple backup systems in Panagrolaimus sp. DAW1, underlining the importance of this sugar in preparation for freezing.

Summary: Functional genomics was used to investigate trehalose synthesis genes after cold acclimation in Panagrolaimus sp. DAW1, an Antarctic nematode with the ability to survive intracellular freezing.

Molecular snapshot of an intracellular freezing event in an Antarctic nematode

The Antarctic nematode, Panagrolaimus sp. DAW1 (formerly called Panagrolaimus davidi), is the best documented example of an organism able to survive intracellular ice formation in all of its compartments. Not only is it able to survive such extreme physiological disruption, but it is able to produce progeny once thawed from such a state. In addition, under slower rates, or less extreme degrees, of cooling, its body remains unfrozen and the vapour pressure difference between the supercooled body fluids and the surrounding ice leads to a process termed cryoprotective dehydration. In contrast to a fairly large body of work in building up our molecular understanding of cryoprotective dehydration, no comparable work has been undertaken on intracellular freezing. This paper describes an experiment subjecting cultures of Panagrolaimus sp. DAW1 to a range of temperatures including a rapid descent to -10 °C, in a medium just prior to, and after, freezing. Through deep sequencing of RNA libraries we have gained a snapshot of which genes are highly abundant when P. sp. DAW1 is undergoing an intracellular freezing event. The onset of freezing correlated with a high production of genes involved in cuticle formation and subsequently, after 24 h in a frozen state, protease production. In addition to the mapping of RNA sequencing, we have focused on a select set of genes arising both from the expression profiles, as well as implicated from other cold tolerance studies, to undertake qPCR. Among the most abundantly represented transcripts in the RNA mapping is the zinc-metalloenzyme, neprilysin, which also shows a particularly strong upregulated signal through qPCR once the nematodes have frozen.

Biodiversity in marine invertebrate responses to acute warming revealed by a comparative multi-omics approach

Understanding species' responses to environmental change underpins our abilities to make predictions on future biodiversity under any range of scenarios. In spite of the huge biodiversity in most ecosystems, a model species approach is often taken in environmental studies. To date, we still do not know how many species we need to study to input into models and inform on ecosystem-level responses to change. In this study, we tested current paradigms on factors setting thermal limits by investigating the acute warming response of six Antarctic marine invertebrates: a crustacean Paraceradocus miersi, a brachiopod Liothyrella uva, two bivalve molluscs, Laternula elliptica, Aequiyoldia eightsii, a gastropod mollusc Marseniopsis mollis and an echinoderm Cucumaria georgiana. Each species was warmed at the rate of 1 °C h^-1 and taken to the same physiological end point (just prior to heat coma). Their molecular responses were evaluated using complementary metabolomics and transcriptomics approaches with the aim of discovering the underlying mechanisms of their resilience or sensitivity to warming. The responses were species-specific; only two showed accumulation of anaerobic end products and three exhibited the classical heat shock response with expression of HSP70 transcripts. These diverse cellular measures did not directly correlate with resilience to heat stress and suggested that each species may have a different critical point of failure. Thus, one unifying molecular mechanism underpinning response to warming could not be assigned, and no overarching paradigm was supported. This biodiversity in response makes future ecosystems predictions extremely challenging, as we clearly need to develop a macrophysiology-type approach to cellular evaluations of the environmental stress response, studying a range of well-rationalized members from different community levels and of different phylogenetic origins rather than extrapolating from one or two arbitrary model species.

Establishing RNAi in a Non-Model Organism: The Antarctic Nematode Panagrolaimus sp. DAW1

The Antarctic nematode Panagrolaimus sp. DAW1 is one of the only organisms known to survive extensive intracellular freezing throughout its tissues. Although the physiological mechanisms of this extreme adaptation are partly understood, the molecular mechanisms remain largely unknown. RNAi is a method that allows the examination of gene function in a direct, targeted manner, by knocking out specific mRNAs and revealing the effects on the phenotype. In this study we have explored the viability of RNAi in Panagrolaimus sp. DAW1. In the first trial, nematodes were fed E. coli expressing Panagrolaimus sp. DAW1 dsRNA of the embryonic lethal genes rps-2 and dhc, and the blister gene duox. Pd-rps-2(RNAi)-treated nematodes showed a significant decrease in larval hatching. However, qPCR showed no significant decrease in the amount of rps-2 mRNA in Pd-rps-2(RNAi)-treated animals. Several soaking protocols for dsRNA uptake were investigated using the fluorescent dye FITC. Desiccation-enhanced soaking showed the strongest uptake of FITC and resulted in a significant and consistent decrease of mRNA levels of two of the four tested genes (rps-2 and tps-2a), suggesting effective uptake of dsRNA-containing solution by the nematode. These findings suggest that RNAi by desiccation-enhanced soaking is viable in Panagrolaimus sp. DAW1 and provide the first functional genomic approach to investigate freezing tolerance in this non-model organism. RNAi, in conjunction with qPCR, can be used to screen for candidate genes involved in intracellular freezing tolerance in Panagrolaimus sp. DAW1.

Beyond connectedness: why pairwise metrics cannot capture community stability

The connectedness of species in a trophic web has long been a key structural characteristic for both theoreticians and empiricists in their understanding of community stability. In the past decades, there has been a shift from focussing on determining the number of interactions to taking into account their relative strengths. The question is: How do the strengths of the interactions determine the stability of a community? Recently, a metric has been proposed which compares the stability of observed communities in terms of the strength of three‐ and two‐link feedback loops (cycles of interaction strengths). However, it has also been suggested that we do not need to go beyond the pairwise structure of interactions to capture stability. Here, we directly compare the performance of the feedback and pairwise metrics. Using observed food‐web structures, we show that the pairwise metric does not work as a comparator of stability and is many orders of magnitude away from the actual stability values. We argue that metrics based on pairwise‐strength information cannot capture the complex organization of strong and weak links in a community, which is essential for system stability.

Transcriptomic SNP discovery for custom genotyping arrays: impacts of sequence data, SNP calling method and genotyping technology on the probability of validation success

BACKGROUND

Single nucleotide polymorphism (SNP) discovery is an important goal of many studies. However, the number of 'putative' SNPs discovered from a sequence resource may not provide a reliable indication of the number that will successfully validate with a given genotyping technology. For this it may be necessary to account for factors such as the method used for SNP discovery and the type of sequence data from which it originates, suitability of the SNP flanking sequences for probe design, and genomic context. To explore the relative importance of these and other factors, we used Illumina sequencing to augment an existing Roche 454 transcriptome assembly for the Antarctic fur seal (Arctocephalus gazella). We then mapped the raw Illumina reads to the new hybrid transcriptome using BWA and BOWTIE2 before calling SNPs with GATK. The resulting markers were pooled with two existing sets of SNPs called from the original 454 assembly using NEWBLER and SWAP454. Finally, we explored the extent to which SNPs discovered using these four methods overlapped and predicted the corresponding validation outcomes for both Illumina Infinium iSelect HD and Affymetrix Axiom arrays.

RESULTS

Collating markers across all discovery methods resulted in a global list of 34,718 SNPs. However, concordance between the methods was surprisingly poor, with only 51.0 % of SNPs being discovered by more than one method and 13.5 % being called from both the 454 and Illumina datasets. Using a predictive modeling approach, we could also show that SNPs called from the Illumina data were on average more likely to successfully validate, as were SNPs called by more than one method. Above and beyond this pattern, predicted validation outcomes were also consistently better for Affymetrix Axiom arrays.

CONCLUSIONS

Our results suggest that focusing on SNPs called by more than one method could potentially improve validation outcomes. They also highlight possible differences between alternative genotyping technologies that could be explored in future studies of non-model organisms.

The transcriptome of metamorphosing flatfish

BACKGROUND

Flatfish metamorphosis denotes the extraordinary transformation of a symmetric pelagic larva into an asymmetric benthic juvenile. Metamorphosis in vertebrates is driven by thyroid hormones (THs), but how they orchestrate the cellular, morphological and functional modifications associated with maturation to juvenile/adult states in flatfish is an enigma. Since THs act via thyroid receptors that are ligand activated transcription factors, we hypothesized that the maturation of tissues during metamorphosis should be preceded by significant modifications in the transcriptome. Targeting the unique metamorphosis of flatfish and taking advantage of the large size of Atlantic halibut (Hippoglossus hippoglossus) larvae, we determined the molecular basis of TH action using RNA sequencing.

RESULTS

De novo assembly of sequences for larval head, skin and gastrointestinal tract (GI-tract) yielded 90,676, 65,530 and 38,426 contigs, respectively. More than 57 % of the assembled sequences were successfully annotated using a multi-step Blast approach. A unique set of biological processes and candidate genes were identified specifically associated with changes in morphology and function of the head, skin and GI-tract. Transcriptome dynamics during metamorphosis were mapped with SOLiD sequencing of whole larvae and revealed greater than 8,000 differentially expressed (DE) genes significantly (p < 0.05) up- or down-regulated in comparison with the juvenile stage. Candidate transcripts quantified by SOLiD and qPCR analysis were significantly (r = 0.843; p < 0.05) correlated. The majority (98 %) of DE genes during metamorphosis were not TH-responsive. TH-responsive transcripts clustered into 6 groups based on their expression pattern during metamorphosis and the majority of the 145 DE TH-responsive genes were down-regulated.

CONCLUSIONS

A transcriptome resource has been generated for metamorphosing Atlantic halibut and over 8,000 DE transcripts per stage were identified. Unique sets of biological processes and candidate genes were associated with changes in the head, skin and GI-tract during metamorphosis. A small proportion of DE transcripts were TH-responsive, suggesting that they trigger gene networks, signalling cascades and transcription factors, leading to the overt changes in tissue occurring during metamorphosis.

Characterisation of the mantle transcriptome and biomineralisation genes in the blunt-gaper clam, Mya truncata

Members of the Myidae family are ecologically and economically important, but there is currently very little molecular data on these species. The present study sequenced and assembled the mantle transcriptome of Mya truncata from the North West coast of Scotland and identified candidate biomineralisation genes. RNA-Seq reads were assembled to create 20,106 contigs in a de novo transciptome, 18.81% of which were assigned putative functions using BLAST sequence similarity searching (cuttoff E-value 1E-10). The most highly expressed genes were compared to the Antarctic clam (Laternula elliptica) and showed that many of the dominant biological functions (muscle contraction, energy production, biomineralisation) in the mantle were conserved. There were however, differences in the constitutive expression of heat shock proteins, which were possibly due to the M. truncata sampling location being at a relatively low latitude, and hence relatively warm, in terms of the global distribution of the species. Phylogenetic analyses of the Tyrosinase proteins from M. truncata showed a gene expansion which was absent in L. elliptica. The tissue distribution expression patterns of putative biomineralisation genes were investigated using quantitative PCR, all genes showed a mantle specific expression pattern supporting their hypothesised role in shell secretion. The present study provides some preliminary insights into how clams from different environments - temperate versus polar - build their shells. In addition, the transcriptome data provides a valuable resource for future comparative studies investigating biomineralisation.

Extremophiles in an Antarctic Marine Ecosystem

Recent attempts to explore marine microbial diversity and the global marine microbiome have indicated a large proportion of previously unknown diversity. However, sequencing alone does not tell the whole story, as it relies heavily upon information that is already contained within sequence databases. In addition, microorganisms have been shown to present small-to-large scale biogeographical patterns worldwide, potentially making regional combinations of selection pressures unique. Here, we focus on the extremophile community in the boundary region located between the Polar Front and the Southern Antarctic Circumpolar Current in the Southern Ocean, to explore the potential of metagenomic approaches as a tool for bioprospecting in the search for novel functional activity based on targeted sampling efforts. We assessed the microbial composition and diversity from a region north of the current limit for winter sea ice, north of the Southern Antarctic Circumpolar Front (SACCF) but south of the Polar Front. Although, most of the more frequently encountered sequences  were derived from common marine microorganisms, within these dominant groups, we found a proportion of genes related to secondary metabolism of potential interest in bioprospecting. Extremophiles were rare by comparison but belonged to a range of genera. Hence, they represented interesting targets from which to identify rare or novel functions. Ultimately, future shifts in environmental conditions favoring more cosmopolitan groups could have an unpredictable effect on microbial diversity and function in the Southern Ocean, perhaps excluding the rarer extremophiles.

Age-related thermal response: the cellular resilience of juveniles

Understanding species' responses to environmental challenges is key to predicting future biodiversity. However, there is currently little data on how developmental stages affect responses and also whether universal gene biomarkers to environmental stress can be identified both within and between species. Using the Antarctic clam, Laternula elliptica, as a model species, we examined both the tissue-specific and age-related (juvenile versus mature adult) gene expression response to acute non-lethal warming (12 h at 3 °C). In general, there was a relatively muted response to this sub-lethal thermal challenge when the expression profiles of treated animals, of either age, were compared with those of 0 °C controls, with none of the "classical" stress response genes up-regulated. The expression profiles were very variable between the tissues of all animals, irrespective of age with no single transcript emerging as a universal biomarker of thermal stress. However, when the expression profiles of treated animals of the different age groups were directly compared, a very different pattern emerged. The profiles of the younger animals showed significant up-regulation of chaperone and antioxidant transcripts when compared with those of the older animals. Thus, the younger animals showed evidence of a more robust cellular response to warming. These data substantiate previous physiological analyses showing a more resilient juvenile population.

Deciphering the molecular adaptation of the king scallop (Pecten maximus) to heat stress using transcriptomics and proteomics

BACKGROUND

The capacity of marine species to survive chronic heat stress underpins their ability to survive warming oceans as a result of climate change. In this study RNA-Seq and 2-DE proteomics were employed to decipher the molecular response of the sub-tidal bivalve Pecten maximus, to elevated temperatures.

RESULTS

Individuals were maintained at three different temperatures (15, 21 and 25 °C) for 56 days, representing control conditions, maximum environmental temperature and extreme warming, with individuals sampled at seven time points. The scallops thrived at 21 °C, but suffered a reduction in condition at 25 °C. RNA-Seq analyses produced 26,064 assembled contigs, of which 531 were differentially expressed, with putative annotation assigned to 177 transcripts. The proteomic approach identified 24 differentially expressed proteins, with nine identified by mass spectrometry. Network analysis of these results indicated a pivotal role for GAPDH and AP-1 signalling pathways. Data also suggested a remodelling of the cell structure, as revealed by the differential expression of genes involved in the cytoskeleton and cell membrane and a reduction in DNA repair. They also indicated the diversion of energetic metabolism towards the mobilization of lipid energy reserves to fuel the increased metabolic rate at the higher temperature.

CONCLUSIONS

This work provides preliminary insights into the response of P. maximus to chronic heat stress and provides a basis for future studies examining the tipping points and energetic trade-offs of scallop culture in warming oceans.

Variability among individuals is generated at the gene expression level

Selection acts on individuals, specifically on their differences. To understand adaptation and responses to change therefore requires knowledge of how variation is generated and distributed across traits. Variation occurs on different biological scales, from genetic through physiological to morphological, yet it is unclear which of these carries the most variability. For example, if individual variation is mainly generated by differences in gene expression, variability should decrease progressively from coding genes to morphological traits, whereas if post-translational and epigenetic effects increase variation, the opposite should occur. To test these predictions, we compared levels of variation among individuals in various measures of gene expression, physiology (including activity), and morphology in two abundant and geographically widespread Antarctic molluscs, the clam Laternula elliptica and the limpet Nacella concinna. Direct comparisons among traits as diverse as heat shock protein QPCR assays, whole transcription profiles, respiration rates, burying rate, shell length, and ash-free dry mass were made possible through the novel application of an established metric, the Wentworth Scale. In principle, this approach could be extended to analyses of populations, communities, or even entire ecosystems. We found consistently greater variation in gene expression than morphology, with physiological measures falling in between. This suggests that variability is generated at the gene expression level. These findings have important implications for refining current biological models and predictions of how biodiversity may respond to climate change.

Key metabolic pathways involved in xenobiotic biotransformation and stress responses revealed by transcriptomics of the mangrove oyster Crassostrea brasiliana

The Brazilian oyster Crassostrea brasiliana was challenged to three common environmental contaminants: phenanthrene, diesel fuel water-accommodated fraction (WAF) and domestic sewage. Total RNA was extracted from the gill and digestive gland, and cDNA libraries were sequenced using the 454 FLX platform. The assembled transcriptome resulted in ̃20,000 contigs, which were annotated to produce the first de novo transcriptome for C. brasiliana. Sequences were screened to identify genes potentially involved in the biotransformation of xenobiotics and associated antioxidant defence mechanisms. These gene families included those of the cytochrome P450 (CYP450), 70kDa heat shock, antioxidants, such as glutathione S-transferase, superoxide dismutase, catalase and also multi-drug resistance proteins. Analysis showed that the massive expansion of the CYP450 and HSP70 family due to gene duplication identified in the Crassostrea gigas genome also occurred in C. brasiliana, suggesting these processes form the base of the Crassostrea lineage. Preliminary expression analyses revealed several candidates biomarker genes that were up-regulated during each of the three treatments, suggesting the potential for environmental monitoring.

Adult acclimation to combined temperature and pH stressors significantly enhances reproductive outcomes compared to short-term exposures

This study examined the effects of long-term culture under altered conditions on the Antarctic sea urchin, Sterechinus neumayeri. Sterechinus neumayeri was cultured under the combined environmental stressors of lowered pH (-0.3 and -0.5 pH units) and increased temperature (+2 °C) for 2 years. This time-scale covered two full reproductive cycles in this species and analyses included studies on both adult metabolism and larval development. Adults took at least 6-8 months to acclimate to the altered conditions, but beyond this, there was no detectable effect of temperature or pH. Animals were spawned after 6 and 17 months exposure to altered conditions, with markedly different outcomes. At 6 months, the percentage hatching and larval survival rates were greatest in the animals kept at 0 °C under current pH conditions, whilst those under lowered pH and +2 °C performed significantly less well. After 17 months, performance was not significantly different across treatments, including controls. However, under the altered conditions urchins produced larger eggs compared with control animals. These data show that under long-term culture adult S. neumayeri appear to acclimate their metabolic and reproductive physiology to the combined stressors of altered pH and increased temperature, with relatively little measureable effect. They also emphasize the importance of long-term studies in evaluating effects of altered pH, particularly in slow developing marine species with long gonad maturation times, as the effects of altered conditions cannot be accurately evaluated unless gonads have fully matured under the new conditions.

Molecular analysis of the cold tolerant Antarctic nematode, Panagrolaimus davidi

Isolated and established in culture from the Antarctic in 1988, the nematode Panagrolaimus davidi has proven to be an ideal model for the study of adaptation to the cold. Not only is it the best-documented example of an organism surviving intracellular freezing but it is also able to undergo cryoprotective dehydration. As part of an ongoing effort to develop a molecular understanding of this remarkable organism, we have assembled both a transcriptome and a set of genomic scaffolds. We provide an overview of the transcriptome and a survey of genes involved in temperature stress. We also explore, in silico, the possibility that P. davidi will be susceptible to an environmental RNAi response, important for further functional studies.

Transcriptome of the Atlantic halibut (Hippoglossus hippoglossus)

Although the Atlantic halibut (Hippoglossus hippoglossus) is an important commercial species, there is still a deficit with regard to the number of transcripts in the databases, which can be accessed and exploited for targeted candidate gene and pathway studies. In this study, the RNAs from head, skin and GI tract from different developmental stages were sequenced to generate 22,272 contigs of 500 base pairs or greater as a molecular resource for this species.

Deep sequencing of the mantle transcriptome of the great scallop Pecten maximus

RNA-Seq transcriptome data were generated from mantle tissue of the great scallop, Pecten maximus. The consensus data were produced from a time course series of animals subjected to a 56-day thermal challenge at 3 different temperatures. A total of 26,064 contigs were assembled de novo, providing a useful resource for both the aquaculture community and researchers with an interest in mollusc shell production.

Interaction strengths in balanced carbon cycles and the absence of a relation between ecosystem complexity and stability

The strength of interactions is crucial to the stability of ecological networks. However, the patterns of interaction strengths in mathematical models of ecosystems have not yet been based upon independent observations of balanced material fluxes. Here we analyse two Antarctic ecosystems for which the interaction strengths are obtained: (1) directly, from independently measured material fluxes, (2) for the complete ecosystem and (3) with a close match between species and 'trophic groups'. We analyse the role of recycling, predation and competition and find that ecosystem stability can be estimated by the strengths of the shortest positive and negative predator-prey feedbacks in the network. We show the generality of our explanation with another 21 observed food webs, comparing random-type parameterisations of interaction strengths with empirical ones. Our results show how functional relationships dominate over average-network topology. They make clear that the classic complexity-instability paradox is essentially an artificial interaction-strength result.

Comparative analysis of a teleost skeleton transcriptome provides insight into its regulation

An articulated endoskeleton that is calcified is a unifying innovation of the vertebrates, however the molecular basis of the structural divergence between terrestrial and aquatic vertebrates, such as teleost fish, has not been determined. In the present study long-read next generation sequencing (NGS, Roche 454 platform) was used to characterize acellular perichondral bone (vertebrae) and chondroid bone (gill arch) in the gilthead sea bream (Sparus auratus). A total of 15.97 and 14.53Mb were produced, respectively from vertebrae and gill arch cDNA libraries and yielded 32,374 and 28,371 contigs (consensus sequences) respectively. 10,455 contigs from vertebrae and 10,625 contigs from gill arches were annotated with gene ontology terms. Comparative analysis of the global transcriptome revealed 4249 unique transcripts in vertebrae, 4201 unique transcripts in the gill arches and 3700 common transcripts. Several core gene networks were conserved between the gilthead sea bream and mammalian skeleton. Transcripts for putative endocrine factors were identified in acellular gilthead sea bream bone suggesting that in common with mammalian bone it can act as an endocrine tissue. The acellular bone of the vertebra, in contrast to current opinion based on histological analysis, was responsive to a short fast and significant (p<0.05) down-regulation of several transcripts identified by NGS, osteonectin, osteocalcin, cathepsin K and IGFI occurred. In gill arches fasting caused a significant (p<0.05) down-regulation of osteocalcin and up-regulation of MMP9.

Transcriptome and peptidome characterisation of the main neuropeptides and peptidic hormones of a euphausiid: the Ice Krill, Euphausia crystallorophias

Background

The Ice krill, Euphausia crystallorophias is one of the species at the base of the Southern Ocean food chain. Given their significant contribution to the biomass of the Southern Ocean, it is vitally important to gain a better understanding of their physiology and, in particular, anticipate their responses to climate change effects in the warming seas around Antarctica.

Methodology/Principal Findings

Illumina sequencing was used to produce a transcriptome of the ice krill. Analysis of the assembled contigs via two different methods, produced 36 new pre-pro-peptides, coding for 61 neuropeptides or peptide hormones belonging to the following families: Allatostatins (A, B et C), Bursicon (α and β), Crustacean Hyperglycemic Hormones (CHH and MIH/VIHs), Crustacean Cardioactive Peptide (CCAP), Corazonin, Diuretic Hormones (DH), the Eclosion Hormone (EH), Neuroparsin, Neuropeptide F (NPF), small Neuropeptide F (sNPF), Pigment Dispersing Hormone (PDH), Red Pigment Concentrating Hormone (RPCH) and finally Tachykinin. LC/MS/MS proteomics was also carried out on eyestalk extracts, which are the major site of neuropeptide synthesis in decapod crustaceans. Results confirmed the presence of six neuropeptides and six precursor-related peptides previously identified in the transcriptome analyses.

Conclusions

This study represents the first comprehensive analysis of neuropeptide hormones in a Eucarida non-decapod Malacostraca, several of which are described for the first time in a non-decapod crustacean. Additionally, there is a potential expansion of PDH and Neuropeptide F family members, which may reflect certain life history traits such as circadian rhythms associated with diurnal migrations and also the confirmation via mass spectrometry of several novel pre-pro-peptides, of unknown function. Knowledge of these essential hormones provides a vital framework for understanding the physiological response of this key Southern Ocean species to climate change and provides a valuable resource for studies into the molecular phylogeny of these organisms and the evolution of neuropeptide hormones.

Identification of molecular and physiological responses to chronic environmental challenge in an invasive species: the Pacific oyster, Crassostrea gigas

Understanding the environmental responses of an invasive species is critical in predicting how ecosystem composition may be transformed in the future, especially under climate change. In this study, Crassostrea gigas, a species well adapted to the highly variable intertidal environment, was exposed to the chronic environmental challenges of temperature (19 and 24°C) and pH (ambient seawater and a reduction of 0.4 pH units) in an extended 3-month laboratory-based study. Physiological parameters were measured (condition index, shell growth, respiration, excretion rates, O:N ratios, and ability to repair shell damage) alongside molecular analyses. Temperature was by far the most important stressor, as demonstrated by reduced condition indexes and shell growth at 24°C, with relatively little effect detected for pH. Transcriptional profiling using candidate genes and SOLiD sequencing of mantle tissue revealed that classical “stress” genes, previously reported to be upregulated under acute temperature challenges, were not significantly expressed in any of the treatments, emphasizing the different response between acute and longer term chronic stress. The transcriptional profiling also elaborated on the cellular responses underpinning the physiological results, including the identification of the PI3K/AKT/mTOR pathway as a potentially novel marker for chronic environmental challenge. This study represents a first attempt to understand the energetic consequences of cumulative thermal stress on the intertidal C. gigas which could significantly impact on coastal ecosystem biodiversity and function in the future.

Cross-amplification and validation of SNPs conserved over 44 million years between seals and dogs

High-density SNP arrays developed for humans and their companion species provide a rapid and convenient tool for generating SNP data in closely-related non-model organisms, but have not yet been widely applied to phylogenetically divergent taxa. Consequently, we used the CanineHD BeadChip to genotype 24 Antarctic fur seal (Arctocephalus gazella) individuals. Despite seals and dogs having diverged around 44 million years ago, 33,324 out of 173,662 loci (19.2%) could be genotyped, of which 173 were polymorphic and clearly interpretable. Two SNPs were validated using KASP genotyping assays, with the resulting genotypes being 100% concordant with those obtained from the high-density array. Two loci were also confirmed through in silico visualisation after mapping them to the fur seal transcriptome. Polymorphic SNPs were distributed broadly throughout the dog genome and did not differ significantly in proximity to genes from either monomorphic SNPs or those that failed to cross-amplify in seals. However, the nearest genes to polymorphic SNPs were significantly enriched for functional annotations relating to energy metabolism, suggesting a possible bias towards conserved regions of the genome.

Hypoxia impacts large adults first: consequences in a warming world

Future oceans are predicted to contain less oxygen than at present. This is because oxygen is less soluble in warmer water and predicted stratification will reduce mixing. Hypoxia in marine environments is thus likely to become more widespread in marine environments and understanding species-responses is important to predicting future impacts on biodiversity. This study used a tractable model, the Antarctic clam, Laternula elliptica, which can live for 36 years, and has a well-characterized ecology and physiology to understand responses to hypoxia and how the effect varied with age. Younger animals had a higher condition index, higher adenylate energy charge and transcriptional profiling indicated that they were physically active in their response to hypoxia, whereas older animals were more sedentary, with higher levels of oxidative damage and apoptosis in the gills. These effects could be attributed, in part, to age-related tissue scaling; older animals had proportionally less contractile muscle mass and smaller gills and foot compared with younger animals, with consequential effects on the whole-animal physiological response. The data here emphasize the importance of including age effects, as large mature individuals appear to be less able to resist hypoxic conditions and this is the size range that is the major contributor to future generations. Thus, the increased prevalence of hypoxia in future oceans may have marked effects on benthic organisms' abilities to persist and this is especially so for long-lived species when predicting responses to environmental perturbation.

Transcriptome of the dead: characterisation of immune genes and marker development from necropsy samples in a free-ranging marine mammal

Background

Transcriptomes are powerful resources, providing a window on the expressed portion of the genome that can be generated rapidly and at low cost for virtually any organism. However, because many genes have tissue-specific expression patterns, developing a complete transcriptome usually requires a 'discovery pool' of individuals to be sacrificed in order to harvest mRNA from as many different types of tissue as possible. This hinders transcriptome development in large, charismatic and endangered species, many of which stand the most to gain from such approaches. To circumvent this problem in a model pinniped species, we 454 sequenced cDNA from testis, heart, spleen, intestine, kidney and lung tissues obtained from nine adult male Antarctic fur seals (Arctocephalus gazella) that died of natural causes at Bird Island, South Georgia.

Results

After applying stringent quality control criteria based on length and annotation, we obtained 12,397 contigs which, in combination with 454 data previously obtained from skin, gave a total of 23,096 unique contigs. Homology was found to 77.0% of dog (Canis lupus familiaris) transcripts, suggesting that the combined assembly represents a substantial proportion of this species' transcriptome. Moreover, only 0.5% of transcripts revealed sequence similarity to bacteria, implying minimal contamination, and the percentage of transcripts involved in cell death was low at 2.6%. Transcripts with immune-related annotations were almost five-fold enriched relative to skin and represented 13.2% of all spleen-specific contigs. By reference to the dog, we also identified transcripts revealing homology to five class I, ten class II and three class III genes of the Major Histocompatibility Complex and derived the putative genomic distribution of 17,121 contigs, 2,119 in silico mined microsatellites and 9,382 single nucleotide polymorphisms.

Conclusions

Our findings suggest that transcriptome development based on samples collected post mortem may greatly facilitate genomic studies, not only of marine mammals but also more generally of species that are of conservation concern.

Thermal reaction norms and the scale of temperature variation: latitudinal vulnerability of intertidal nacellid limpets to climate change

The thermal reaction norms of 4 closely related intertidal Nacellid limpets, Antarctic (Nacella concinna), New Zealand (Cellana ornata), Australia (C. tramoserica) and Singapore (C. radiata), were compared across environments with different temperature magnitude, variability and predictability, to test their relative vulnerability to different scales of climate warming. Lethal limits were measured alongside a newly developed metric of "duration tenacity", which was tested at different temperatures to calculate the thermal reaction norm of limpet adductor muscle fatigue. Except in C. tramoserica which had a wide optimum range with two break points, duration tenacity did not follow a typical aerobic capacity curve but was best described by a single break point at an optimum temperature. Thermal reaction norms were shifted to warmer temperatures in warmer environments; the optimum temperature for tenacity (T(opt)) increased from 1.0°C (N. concinna) to 14.3°C (C. ornata) to 18.0°C (an average for the optimum range of C. tramoserica) to 27.6°C (C. radiata). The temperature limits for duration tenacity of the 4 species were most consistently correlated with both maximum sea surface temperature and summer maximum in situ habitat logger temperature. Tropical C. radiata, which lives in the least variable and most predictable environment, generally had the lowest warming tolerance and thermal safety margin (WT and TSM; respectively the thermal buffer of CT(max) and T(opt) over habitat temperature). However, the two temperate species, C. ornata and C. tramoserica, which live in a variable and seasonally unpredictable microhabitat, had the lowest TSM relative to in situ logger temperature. N. concinna which lives in the most variable, but seasonally predictable microhabitat, generally had the highest TSMs. Intertidal animals live at the highly variable interface between terrestrial and marine biomes and even small changes in the magnitude and predictability of their environment could markedly influence their future distributions.

Metagenomic analysis of a southern maritime antarctic soil

Our current understanding of Antarctic soils is derived from direct culture on selective media, biodiversity studies based on clone library construction and analysis, quantitative PCR amplification of specific gene sequences and the application of generic microarrays for microbial community analysis. Here, we investigated the biodiversity and functional potential of a soil community at Mars Oasis on Alexander Island in the southern Maritime Antarctic, by applying 454 pyrosequencing technology to a metagenomic library constructed from soil genomic DNA. The results suggest that the commonly cited range of phylotypes used in clone library construction and analysis of 78–730 OTUs (de-replicated to 30–140) provides low coverage of the major groups present (∼5%). The vast majority of functional genes (>77%) were for structure, carbohydrate metabolism, and DNA/RNA processing and modification. This study suggests that prokaryotic diversity in Antarctic terrestrial environments appears to be limited at the generic level, with Proteobacteria, Actinobacteria being common. Cyanobacteria were surprisingly under-represented at 3.4% of sequences, although ∼1% of the genes identified were involved in CO₂ fixation. At the sequence level there appeared to be much greater heterogeneity, and this might be due to high divergence within the relatively restricted lineages which have successfully colonized Antarctic terrestrial environments.

Estimating long-term survival temperatures at the assemblage level in the marine environment: towards macrophysiology

Defining ecologically relevant upper temperature limits of species is important in the context of environmental change. The approach used in the present paper estimates the relationship between rates of temperature change and upper temperature limits for survival in order to evaluate the maximum long-term survival temperature (Ts). This new approach integrates both the exposure time and the exposure temperature in the evaluation of temperature limits. Using data previously published for different temperate and Antarctic marine environments, we calculated Ts in each environment, which allowed us to calculate a new index: the Warming Allowance (WA). This index is defined as the maximum environmental temperature increase which an ectotherm in a given environment can tolerate, possibly with a decrease in performance but without endangering survival over seasonal or lifetime time-scales. It is calculated as the difference between maximum long-term survival temperature (Ts) and mean maximum habitat temperature. It provides a measure of how close a species, assemblage or fauna are living to their temperature limits for long-term survival and hence their vulnerability to environmental warming. In contrast to data for terrestrial environments showing that warming tolerance increases with latitude, results here for marine environments show a less clear pattern as the smallest WA value was for the Peru upwelling system. The method applied here, relating upper temperature limits to rate of experimental warming, has potential for wide application in the identification of faunas with little capacity to survive environmental warming.

Physiological plasticity, long term resistance or acclimation to temperature, in the Antarctic bivalve, Laternula elliptica

To further investigate the previously reported limited acclimation capacities of Antarctic marine stenotherms, the Antarctic mud clam, Laternula elliptica (King and Broderip, 1830-1831), was incubated at 3.0°C for 89days. The thermal windows of a suite of biochemical and physiological metrics that characterise tissue aerobic status, were then measured in response to acute temperature elevation (2-2.5°C increase per week). To test if acclimation had occurred at the higher temperature, results were compared with published data, from the preceding year, for L. elliptica which had been incubated at ambient temperature (0.0°C) and then subjected to the same acute temperature treatments. Incubation to 3.0°C led to a temperature induced increase of tissue aerobic status (reduced intracellular cCO(2) with increased O(2) consumption, PLA (phospho-L-arginine) and ATP). At the highest acute temperature (7.5°C) the increase in anaerobic pathways (summed acetate/succinate and propionate) was less after 3.0°C than 0.0°C incubation. No other metric shifted its reaction norm in response to acute temperature elevation and so whole animal acclimation had not occurred, even after 3months at 3.0°C. Combined with the constant mortality throughout the 3.0°C incubation period, these data suggest that the recorded physiological changes were either the early stages of acclimation or, more likely, time limited resistance mechanisms.

Long-term survival of hydrated resting eggs from Brachionus plicatilis

Background

Several organisms display dormancy and developmental arrest at embryonic stages. Long-term survival in the dormant form is usually associated with desiccation, orthodox plant seeds and Artemia cysts being well documented examples. Several aquatic invertebrates display dormancy during embryonic development and survive for tens or even hundreds of years in a hydrated form, raising the question of whether survival in the non-desiccated form of embryonic development depends on pathways similar to those occurring in desiccation tolerant forms.

Methodology/Principal Findings

To address this question, Illumina short read sequencing was used to generate transcription profiles from the resting and amictic eggs of an aquatic invertebrate, the rotifer, Brachionus plicatilis. These two types of egg have very different life histories, with the dormant or diapausing resting eggs, the result of the sexual cycle and amictic eggs, the non-dormant products of the asexual cycle. Significant transcriptional differences were found between the two types of egg, with amictic eggs rich in genes involved in the morphological development into a juvenile rotifer. In contrast, representatives of classical “stress” proteins: a small heat shock protein, ferritin and Late Embryogenesis Abundant (LEA) proteins were identified in resting eggs. More importantly however, was the identification of transcripts for messenger ribonucleoprotein particles which stabilise RNA. These inhibit translation and provide a valuable source of useful RNAs which can be rapidly activated on the exit from dormancy. Apoptotic genes were also present. Although apoptosis is inconsistent with maintenance of prolonged dormancy, an altered apoptotic pathway has been proposed for Artemia, and this may be the case with the rotifer.

Conclusions

These data represent the first transcriptional profiling of molecular processes associated with dormancy in a non-desiccated form and indicate important similarities in the molecular pathways activated in resting eggs compared with desiccated dormant forms, specifically plant seeds and Artemia.

Skin healing and scale regeneration in fed and unfed sea bream, Sparus auratus

Background

Fish scales are an important reservoir of calcium and phosphorus and together with the skin function as an integrated barrier against environmental changes and external aggressors. Histological studies have revealed that the skin and scales regenerate rapidly in fish when they are lost or damaged. In the present manuscript the histological and molecular changes underlying skin and scale regeneration in fed and fasted sea bream (Sparus auratus) were studied using a microarray 3 and 7 days after scale removal to provide a comprehensive molecular understanding of the early stages of these processes.

Results

Histological analysis of skin/scales revealed 3 days after scale removal re-epithelisation and formation of the scale pocket had occurred and 53 and 109 genes showed significant up or down-regulation, respectively. Genes significantly up-regulated were involved in cell cycle regulation, cell proliferation and adhesion, immune response and antioxidant activities. 7 days after scale removal a thin regenerated scale was visible and only minor changes in gene expression occurred. In animals that were fasted to deplete mineral availability the expression profiles centred on maintaining energy homeostasis. The utilisation of fasting as a treatment emphasised the competing whole animal physiological requirements with regard to barrier repair, infection control and energy homeostasis.

Conclusions

The identification of numerous genes involved in the mitotic checkpoint and cell proliferation indicate that the experimental procedure may be useful for understanding cell proliferation and control in vertebrates within the context of the whole animal physiology. In response to skin damage genes of immune surveillance were up-regulated along with others involved in tissue regeneration required to rapidly re-establish barrier function. Additionally, candidate fish genes were identified that may be involved in cytoskeletal re-modelling, mineralization and stem cells, which are of potential use in aquaculture and fish husbandry, as they may impact on the ability of the fish to produce structural proteins, such as muscle, efficiently.

Proteomics of cryoprotective dehydration in Megaphorura arctica Tullberg 1876 (Onychiuridae: Collembola)

The Arctic springtail, Megaphorura arctica Tullberg 1876 (Onychiuridae: Collembola), is one of the few organisms known to survive the extreme stresses of its environment by using cryoprotective dehydration. We have undertaken a proteomics study comparing M. arctica, acclimated at -2°C, the temperature known to induce the production of the anhydroprotectant trehalose in this species, and -6°C, the temperature at which trehalose expression plateaus, against control animals acclimated at +5°C. Using difference gel electrophoresis, and liquid chromatography tandem mass spectrometry, we identified three categories of differentially expressed proteins with specific functions, up-regulated in both the -2°C and -6°C animals, that were involved in metabolism, membrane transport and protein folding. Proteins involved in cytoskeleton organisation were only up-regulated in the -6°C animals.

Antarctic krill 454 pyrosequencing reveals chaperone and stress transcriptome

BACKGROUND

The Antarctic krill Euphausia superba is a keystone species in the Antarctic food chain. Not only is it a significant grazer of phytoplankton, but it is also a major food item for charismatic megafauna such as whales and seals and an important Southern Ocean fisheries crop. Ecological data suggest that this species is being affected by climate change and this will have considerable consequences for the balance of the Southern Ocean ecosystem. Hence, understanding how this organism functions is a priority area and will provide fundamental data for life history studies, energy budget calculations and food web models.

METHODOLOGY/PRINCIPAL FINDINGS

The assembly of the 454 transcriptome of E. superba resulted in 22,177 contigs with an average size of 492bp (ranging between 137 and 8515bp). In depth analysis of the data revealed an extensive catalogue of the cellular chaperone systems and the major antioxidant proteins. Full length sequences were characterised for the chaperones HSP70, HSP90 and the super-oxide dismutase antioxidants, with the discovery of potentially novel duplications of these genes. The sequence data contained 41,470 microsatellites and 17,776 Single Nucleotide Polymorphisms (SNPs/INDELS), providing a resource for population and also gene function studies.

CONCLUSIONS

This paper details the first 454 generated data for a pelagic Antarctic species or any pelagic crustacean globally. The classical "stress proteins", such as HSP70, HSP90, ferritin and GST were all highly expressed. These genes were shown to be over expressed in the transcriptomes of Antarctic notothenioid fish and hypothesized as adaptations to living in the cold, with the associated problems of decreased protein folding efficiency and increased vulnerability to damage by reactive oxygen species. Hence, these data will provide a major resource for future physiological work on krill, but in particular a suite of "stress" genes for studies understanding marine ectotherms' capacities to cope with environmental change.

Transcription profiling of acute temperature stress in the Antarctic plunderfish Harpagifer antarcticus

Harpagifer antarcticus (the Antarctic plunderfish), a shallow-water benthic fish distributed around the Antarctic Peninsula, is a member of the notothenioid family, one of whose adaptations to the cold waters of Antarctica has been the loss of the classic heat shock response. In order to gain a more comprehensive understanding of the effects of temperature stress on H. antarcticus, we constructed a liver cDNA library and a 10,371 feature microarray. This was hybridized with material from a time course series of animals held at 6°C for 48h. The resulting expression profiles show that this fish displays the classical vertebrate acute inflammatory response. There was also a pronounced signal for increased energy requirements via up-regulation of genes involved in the β oxidation of fatty acids and also a strong signature of response to oxidative stress. Genes in the latter category did not include the "classic" antioxidants such as glutathione S-transferase, but genes involved in the production of reducing potential in the form of NADPH, peroxisome proliferation via peroxisomal acyl co-enzyme A oxidase 1 and genes known to be up-regulated by hypoxia-inducible factor 1 (HIF1). These identifications provide clear support for oxygen being the whole animal limiting factor at least in acute short-term temperature challenges. The classical heat shock proteins were not up-regulated during this trial, although numerous clones for each were present on the gene chip, confirming the lack of this response in this species. These data significantly increase our knowledge of the cellular stress response from animals in this unique environment.

Gene expression associated with changes in cold tolerance levels of the Antarctic springtail, Cryptopygus antarcticus

The ability of the Antarctic microarthropod Cryptopygus antarcticus (Collembola, Isotomidae) to survive low temperatures has been well studied at the physiological level, with recent investigations indicating the importance of the moulting process in conferring this ability. This study investigated gene expression in groups of C. antarcticus that have distinct differences in their ability to survive low temperatures. A microarray containing c. 5400 C. antarcticus expressed sequence tags was used to investigate gene expression differences between groups of animals with different supercooling points (SCP), and to low temperatures close to their SCP. By demonstrating the involvement of moult-related genes in the differential survival of two groups of C. antarcticus with distinct SCP profiles, the results of this investigation add support to the suggestion that moulting plays a role in conferring cold tolerance in C. antarcticus.

Cold hardening processes in the Antarctic springtail, Cryptopygus antarcticus: clues from a microarray

The physiology of the Antarctic microarthropod, Cryptopygus antarcticus, has been well studied, particularly with regard to its ability to withstand low winter temperatures. However, the molecular mechanisms underlying this phenomenon are still poorly understood. 1180 sequences (Expressed Sequence Tags or ESTs) were generated and analysed, from populations of C. antarcticus. This represents the first publicly available sequence data for this species. A sub-set (672 clones) were used to generate a small microarray to examine the differences in gene expression between summer acclimated cold tolerant and non-cold tolerant springtails. Although 60% of the clones showed no sequence similarity to annotated genes in the datasets, of those where putative function could be inferred via database homology, there was a clear pattern of up-regulation of structural proteins being associated with the cold tolerant group. These structural proteins mainly comprised cuticle proteins and provide support for the recent theory that summer SCP variation within Collembola species could be a consequence of moulting, with moulting population having lowered SCPs.