- Research article
- Open Access
Desiccation survival in an Antarctic nematode: molecular analysis using expressed sequenced tags
© Adhikari et al; licensee BioMed Central Ltd. 2009
- Received: 15 September 2008
- Accepted: 09 February 2009
- Published: 09 February 2009
Nematodes are the dominant soil animals in Antarctic Dry Valleys and are capable of surviving desiccation and freezing in an anhydrobiotic state. Genes induced by desiccation stress have been successfully enumerated in nematodes; however we have little knowledge of gene regulation by Antarctic nematodes which can survive multiple environmental stresses. To address this problem we investigated the genetic responses of a nematode species, Plectus murrayi, that is capable of tolerating Antarctic environmental extremes, in particular desiccation and freezing. In this study, we provide the first insight into the desiccation induced transcriptome of an Antarctic nematode through cDNA library construction and suppressive subtractive hybridization.
We obtained 2,486 expressed sequence tags (ESTs) from 2,586 clones derived from the cDNA library of desiccated P. murrayi. The 2,486 ESTs formed 1,387 putative unique transcripts of which 523 (38%) had matches in the model-nematode Caenorhabditis elegans, 107 (7%) in nematodes other than C. elegans, 153 (11%) in non-nematode organisms and 605 (44%) had no significant match to any sequences in the current databases. The 1,387 unique transcripts were functionally classified by using Gene Ontology (GO) hierarchy and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The results indicate that the transcriptome contains a group of transcripts from diverse functional areas. The subtractive library of desiccated nematodes showed 80 transcripts differentially expressed during desiccation stress, of which 28% were metabolism related, 19% were involved in environmental information processing, 28% involved in genetic information processing and 21% were novel transcripts. Expression profiling of 14 selected genes by quantitative Real-time PCR showed 9 genes significantly up-regulated, 3 down-regulated and 2 continuously expressed in response to desiccation.
The establishment of a desiccation EST collection for Plectus murrayi, a useful model in assessing the structural, physiological, biochemical and genetic aspects of multiple stress tolerance, is an important step in understanding the genome level response of this nematode to desiccation stress. The type of transcript analysis performed in this study sets the foundation for more detailed functional and genome level analyses of the genes involved in desiccation tolerance in nematodes.
- Gene Ontology
- Differentially Express
- Late Embryogenesis Abundant
- Desiccation Tolerance
- Enzyme Commission
The Dry Valleys of Antarctica are one of the most extreme terrestrial environments on Earth . Soils in this cold desert ecosystem are subjected to freezing temperatures, desiccation and salt accumulation that affect biological water availability [2, 3]. Soil communities in Antarctic Dry Valleys are simple; primary production is largely limited to algae, and fauna are almost exclusively microbial grazers (mostly protozoa, rotifers, tardigrades and nematodes ). Nematodes are the dominant soil animals, present in 65% of the 415 soils sampled by Wall Freckman & Virginia  across four McMurdo Dry Valleys (MCM). Nematodes have been isolated from soil in an inactive coiled state called anhydrobiosis . Anhydrobiosis is a survival strategy employed by nematodes, rotifers, and tardigrades in response to desiccation . Nematodes in anhydrobiosis lose 95–99% of their body water content and can cease metabolic activity at any stage in their life cycle . While in an anhydrobiotic state, nematodes are capable of surviving desiccation  as well as extreme cold . Though Antarctic ecosystems are simple and have low species diversity compared to temperate ecosystems, nematodes are the most widely distributed and biologically diverse invertebrates in the Dry Valleys , with four species; Scottnema lindsayae, Eudorylaimus antarcticus, Plectus antarcticus, and Geomonhystera antarcticola . It has been suggested that specimens identified as P. antarcticus de Man 1904 in MCM are P. murrayi [11, 12] and we accept this nomenclature for the present paper.
Plectus murrayi, a bacteria feeding nematode , inhabits both semi-aquatic and terrestrial biotopes in the Dry Valleys, but is also reported from other parts of the Antarctica . Similar to S. lindsaye, another nematode endemic to the Southern continent , P. murrayi has a multiple year life cycle . The distribution of these nematodes in Antarctica is dependent on organic carbon and soil moisture  with high abundance in stream sediments . P. murrayi from the MCM are freeze tolerant, and can tolerate repeated freeze-thaw cycles in the laboratory (data not shown). Although adapted to the extreme desiccation and freezing encountered in its habitat , the biology and environmental tolerance of this nematode has not been well studied.
Despite recent work on behavioral, biochemical and molecular stress response mechanisms [17–19] the molecular mechanisms governing anhydrobiosis in nematodes are not fully understood. Anhydrobiosis in nematodes is reported to involve the biosynthesis of low molecular weight carbohydrates, proteins and glycerol [20, 21]. Recent research suggests anhydrobiotes synthesize many other compounds (primarily proteins) that are essential for survival [22–24]. Studies on these desiccation responsive compounds have resulted in the identification of many genes that play important roles in stress acclimation and survival. These responses include up-regulation of transcriptional regulators, molecular chaperones, antioxidants, hydrophilic proteins, and proteins involved in cell cycle regulation [19, 25–27]. The anhydrobiotic nematode Aphelenchus avenae synthesizes large amounts of trehalose in response to desiccation . However, it has become clear that such sugars are not sufficient for anhydrobiosis  and, indeed, that some anhydrobiotic organisms seem not to use them . As an effort to identify other adaptations required for anhydrobiosis Goyal et al.  characterised genes in the nematode A. avenae that requires a period of preconditioning to enter anhydrobiosis. During this preconditioning period, several genes, including trehalose synthase , hydrophilins (highly hydrophilic proteins), anhydrin, and a polypeptide, Aav-LEA-1, related to plant Group 3 late embryogenesis abundant (LEA) proteins were induced [23, 19]. Although similar gene classes were found to be associated with desiccation stress in many nematodes, none of the ESTs or proteins detected in these studies were encoded by the same gene  and their expression level was quite variable . To understand such molecular mechanisms activated during anhydrobiosis, a condition induced by slow dehydration, we identified gene expression patterns by gradually desiccating nematodes at relative humidity (RH) values reflective of the Antarctic environment.
To understand the mechanisms of desiccation survival we have initiated a genomic level analysis of gene expression during anhydrobiosis of P. murrayi. The first step in this process was to establish an EST collection that is representative of the desiccation induced transcripts and to identify the transcripts differentially expressed during desiccation stress. Here we present bioinformatics and molecular analysis of 2,486 ESTs from the gradually desiccated and anhydrobiotically induced nematode P. murrayi and 80 transcripts differentially expressed during the anhydrobiotic process. The bioinformatics approaches include EST cluster analyses, transcript abundancy estimations, and functional classifications based on Inter-Pro domains, Gene Ontology hierarchy, and KEGG biochemical classifications. The genetic information derived from P. murrayi informs the characterization of genes responding to desiccation stress, and is expected to further our understanding of the potential genetic determinants of desiccation tolerance in nematodes and perhaps other metazoans.
Sequencing and assembly of ESTs
Plectus murrayi EST summary
Total number of high quality sequences‡
Average length of sequences (bp)†
545 ± 156
Number of contigs§
Number of singletons
Number of putative unique transcripts¶
Unique transcripts with similarity to C. elegans database
Unique transcripts with similarity to other nematode database
Unique transcripts with similarity to other organisms
Total unique transcripts with significant similarity
Unique transcripts with no significant similarity
Comparison against public nematode ESTs
We used the 1,387 unique sequences to search a non-redundant protein data base using BLASTX [33, 34] (Table 1) and the Wormpep 190 database consisting of extensively curated C. elegans proteins from WormBase . A total of 782 unique sequences (56%) matched known proteins, including 523 unique sequences (38%) with significant match to C. elegans proteins at a cut-off expectation (E)-value of 10-5 or below. The remainder of the unique sequences (44%) had no meaningful matches (E > 10-5). We compared our unique sequences with the ESTs from other nematodes as well as non-nematodes using BLAST searches. Only 107 unique sequences (7%) matched other nematode ESTs and 153 unique sequences (11%) matched organisms other than nematodes at E < 10-5 (Table 1). Of 1,387 unique sequences, 36 had homologues in C. elegans which could be silenced by RNAi. The RNAi phenotypes (as described by WormBase) included mig-15 (migration), lin-8 (lineage), unc-16, 89 (uncoordinated), dpy-6 (dumpy), rde-1 (RNAi defective), drh-2 (dicer related), nhr-67 (nuclear hormone receptor) and ard-1 (alcohol/ribitol dehydrogenase).
Identification of differentially expressed genes
Listing of ESTs differentially expressed during desiccation of Plectus murrayi and their homologs in GenBank.
GenBank accession number
Annotation (Organism) number
Metabolism ɠ (22)
GPD family member [Caenorhabditis elegans]
ALDH family member [C. elegans]
Glutathione S-transferase [Haemonchus contortus]
Lipid Transfer protein family member [C. elegans]
Aspartyl protease protein 6 [C. elegans]
Plant LEA related family member [C. elegans]
TPS6 protein 1 [Brugia malayi]
FBA1, putative [B. malayi]
Thymidylate synthase [C. elegans]
Zinc finger protein family member [C. elegans]
Probable glycerol kinase [C. elegans]
MDH1, NAD (soluble) [Gallus gallus]
IDH [Magnetospirillum magnetotacticum]
ATP synthase subunit family member [C. elegans]
NADH-ubiquinone oxidoreductase [B. malayi]
Glycogen synthase family member [C. elegans]
ADP/ATP translocase [C. elegans]
Bi-functional glyoxylate cycle protein [C. elegans]
Putative glutamate dehydrogenase [H. contortus]
Thymidylate synthase [C. elegans]
Glutathione peroxidase [C. elegans]
ATP synthase sub unit family member [C. elegans]
Environmental information processing (15)
Cathepsin L cysteine protease [H. contortus]
JNK kinase family member jkk-1 [C. elegans]
Superoxide dismutase family member [C.elegans]
Type II antifreeze protein [Clupea harengus]
Heat shock protein 70 A [Heterodera glycines]
Small heat shock protein family member [C. elegans]
RAB family member [C. elegans]
Ras-related protein Rab-11B, putative [B. malayi]
Heat shock protein family member [C. elegans]
Inhibitor of apoptosis homolog [C. elegans]
Heat shock protein 90 protein [B. malayi]
DumPY: shorter than wild-type [C.elegans]
Collagen protein 170 [C. elegans]
Leucine rich repeat family protein [B. malayi]
Laminin receptor 1 [Xenopus laevis]
Genetic information processing (23)
AP inhibitor [Parelaphostrongylus tenuis]
Putative E2 enzyme [Oesophagostomum dentatum]
Histone H2B 2, putative [B. malayi]
Peroxiredoxin [H. contortus]
Ubiquitin specific protease 14 [D. rerio]
Ribosomal protein, LSU family member [C. elegans]
Ribosomal protein, SSU family member [C. elegans]
Ribosomal protein, SSU family member [C. elegans]
Ribosomal protein, SSU family member [C. elegans]
Ribosomal protein, LSU family member [C. elegans]
Ribosomal protein, LSU family member [C. elegans]
Ribosomal protein, SSU family member [C. elegans]
Ribosomal protein, LSU family member [C. elegans]
60S ribosomal protein L27a, putative [B. malayi]
60S ribosomal protein L39, putative [B. malayi]
40S ribosomal protein S6, putative [B. malayi]
EF family member [C. elegans]
EF1-alpha, putative [B. malayi]
EF family member [C. elegans]
EF1 beta, isoform A [Drosophila melanogaster]
EF family member [C. elegans]
NHR family member [C. elegans]
Transcription factor, putative [B. malayi]
Hypothetical proteins (3)
Hypothetical protein [C. briggsae]
Hypothetical protein [C. briggsae]
Predicted protein [X. laevis]
Novel proteins (17)
The most unexpected discovery among DE ESTs was Type II antifreeze protein (AFP) [GenBank: FK670242 which showed high similarity to that of Clupea harengus (Table 2), the Atlantic herring. This finding represents another case of an ice structuring protein from an Antarctic nematode, suggesting the possibility that Antarctic nematodes may use similar antifreeze proteins for stress adaptation heretofore observed only in one Antarctic nematode , some fishes , insects , plants  and fungi and bacteria .
Abundant transcripts expressed during desiccation
The most abundantly represented transcripts in the cDNA library
Number of ESTs
Cytochrome c oxidase subunit 2
Small heat shock protein family
Y-box family member
Zinc finger protein
Cu/Zn superoxide dismutase
Translation initiation factor
Chaperonin containing subunit
NADH-dehydrogenase subunit 1
Cathespin b -like cysteine proteinase
Elongation factor 1 alpha
Heat shock 70 kda protein
60S ribosomal protein
Elongation factor family member
Heat shock 90 kda protein
ATP synthase subunit
Nuclear hormone protein family
Functional classification based on gene ontology assignments
Functional classification based on KEGG analysis
KEGG biochemical mappings for Plectus murrayi clusters.
KEGG categories represented
Unique sequences (Number of enzymes)
Amino acid metabolism
Xenobiotics biodegradation and metabolism
Biosynthesis of secondary metabolites
Metabolism of other amino acids
Glyoxylate and dicarboxylate metabolism
Genetic information processing
Folding, sorting and degradation
Replication and repair
Environmental information processing
Signalling molecules and interaction
Cell growth and death
Refined gene-specific expression using quantitative real-time PCR
The expressed genome of P. murrayi showed that anhydrobiotic survival in nematodes involves a suite of genes from diverse functional areas, such as hormone signaling transduction, transcription regulation, ROS scavenging, reestablishment of homeostasis, molecular chaperoning and transcriptional regulation of ribosomal proteins and other genes. The 2,486 sequences comprised 1,387 putative unique sequences from P. murrayi and represent the first EST analysis from an Antarctic nematode. The number of ESTs analyzed in the study were relatively few; nonetheless, 44% of the unique transcripts (605) expressed by desiccated nematodes represent novel genes, which falls within the range reported from several other nematode EST analyses . Comparison of ESTs presented in this study with ESTs from GenBank showed that 56% of the unique transcripts (1,387) isolated from P. murrayi have been previously isolated from other organisms, including the model organism C. elegans (Table 1). Subtractive hybridization of cDNAs from desiccated and undesiccated nematodes used in library construction was done to enrich for the rare transcripts that are differentially expressed (DE). A total of 80 DE transcripts were identified, some of which are involved in genetic information processing followed by metabolism, and presumably are involved in stress survival of nematodes (Table 2). The expression level of 14 transcripts DE during desiccation survival using qRT-PCR indicates that these genes were DE between desiccated and undesiccated nematodes. The results from subtractive hybridization coupled with qRT-PCR showed that we appropriately enriched for the DE transcripts (Fig. 3). Our results showed that expressed sequence tag (EST) analysis coupled with subtractive hybridization is a powerful method to identify the genes involved in nematode desiccation stress. EST analysis is a commonly-used approach to identify genes involved in specific biological functions, especially in organisms where genomic data are not available .
Functional analysis of expressed genes
Gene ontology has been widely used to characterize gene function annotation and classification . GO describes gene function using controlled vocabulary and hierarchy, including molecular function, biological processes, and cellular components (Fig. 2; Additional file 1). A large number of the unique sequences from our study mapped to molecular functions and ligand binding, enzymes, molecular transducer and transporter subcategories. Each one of these subcategories represents catalytic activities that could be argued are important for a cell to survive major metabolic perturbation during desiccation. The second most represented category was biological process, which includes unique sequences associated with cell growth or maintenance, development and cellular communications. Within the cell growth or maintenance category most unique sequences were associated with metabolism, localization, transport and response to stress. This distribution is not surprising for ESTs (unique sequences) derived from an organism undergoing metabolic changes such as desiccation. The subcategories under the least represented category, cellular process, also seem to reflect the nature of the cellular disturbances that result from desiccation. Under this category there is significant representation under intracellular and membrane subcategories. Within these subcategories representation is most significant in mitochondria, cytosol, ribosome, nucleus and ribonucleoprotein complex. The importance of protein synthesis and cellular homeostastasis during desiccation may explain the preponderance of unique sequences associated with mitochondrial and ribosomal structural components.
As an alternative method of categorizing clusters by biochemical function, clusters were assigned to metabolic pathways using the KEGG database  (Table 4). Only 36% of the unique sequences mapped to the currently known KEGG pathways with 158 unique mappings. The paucity of EC assignments limits this aspect of analysis but the mapping of unique sequences to the KEGG metabolic and other pathways still presents some useful perspectives on the metabolic, protein folding and degradation and cellular repair emphasis of desiccating cells. Most of the unique transcripts belonged to genetic information processing (GIP) with protein folding, sorting and degradation. Transcripts encoding a cathespin L-like protease [GenBank: FK670238] enzyme found in this category have long been recognized for their role in intracellular and extracellular protein degradation in a range of cellular processes. In C. elegans cathespin L-protease (Ce-cpl) has been demonstrated to play a critical role during embryogenesis, larval development and moulting . Considering the putative homology of the transcript encoding P. murrayi cathespin L-protease (Pm-cpl) with C. elegans, it might have a similar role in nematode molting and thus its possible role during desiccation survival needs further investigation.
Metabolism was the second most represented category and pathways well-represented by the P. murrayi clusters were carbohydrate metabolism, amino acid metabolism, lipid metabolism, biosynthesis of secondary metabolites, glyoxylate, and decarboxylation metabolism (Table 4). The lipid metabolism pathway in anhydrobiotes is one of the most active pathways as lipids are the main reservoir of energy and the most likely source of carbon for the synthesis of trehalose . The GIP category included the pathways involved in protein folding, sorting and degradation, transcription, translation, replication and repair. Desiccation stress coupled with oxidative stress results in lipid and protein damage, leading to impaired cellular survival and functioning. Nematodes are reported to respond to these damages and repair them with increased expression of stress-protective genes and antioxidant enzymes . The cellular processes category has the least number of transcripts. This finding is suggestive of developmental arrest during anhydrobiosis. Indeed, anhydrobiosis may be an exaptation that allows the individual to survive unfavourable conditions by staying in an arrested state of development and reproduction (Table 4).
Abundantly expressed transcripts during anhydrobiosis
A high level of representation in a cDNA library generally correlates with high transcript abundance in the original biological sample , although artefacts of library construction can result in selection for or against representation of some transcripts. The genes with the most abundant transcripts were mostly involved in metabolism, molecular chaperones, reactive oxygen species scavenging and genetic information processing (Table 3). Ribosomal protein (RP) (GenBank: FG618944) was the most abundant transcript expressed during desiccation, suggesting the possible involvement of transcriptional regulation of ribosomal proteins during desiccation stress. Differential expression of plant ribosomal protein genes has been observed during development and following various stress or hormone treatments, including desiccation stress . The 'ribosomal filter hypothesis' of Mauro and Edelman  proposes that particular combinations of ribosomal proteins or rRNA could favour the translation of specific mRNAs, thereby providing a mechanism for translational control. Interestingly, many ribosomal proteins were DE and one of them was up-regulated (Pm-rpl-4), suggesting that RPs may provide a means for selectively translating specific mRNAs required for the desiccation response as suggested for the insect parasitic nematode Steinernema carpocapsae . Similarly, it may be part of an adaptive response to maintain ribosomal function in dehydrated cytoplasm.
A number of metabolism related transcripts encoding cytochrome c oxidase subunit 2 (GenBank: FG618989), glutamate synthase (GenBank: FG620419), NADH-dehydrogenase subunit 1 (GenBank: FG619004) and aldehyde dehydrogenase (GenBank: FG619103) were also abundantly expressed during desiccation stress. Several lines of evidence have suggested that nematodes activate metabolic pathways in response to desiccation shortly after exposure to dehydrating conditions .
Transcripts of a number of stress related genes were also abundantly expressed, including heat shock proteins, a chaperonin containing sub-unit, multiple stress-responsive zinc finger proteins and oxidative stress responsive genes. Desiccation stress in nematodes reportedly produces a large number of molecular chaperones to facilitate the synthesis, folding, assembly and intracellular transport of proteins, reduce protein denaturation and aggregation, and aid in protein renaturation [19, 52]. Transcripts of several oxidative stress related genes like glutathione S-transferase (GenBank: FG618951) and copper/zinc superoxide dismutase (GenBank: FG618942) were also abundantly expressed. Reactive oxygen species and other toxins produced by oxidative stress during desiccation of nematodes can damage membrane systems, proteins and nucleic acids. Therefore, the transcript abundance of several proteins that contribute to cellular survival after oxidative damage is not surprising.
Genes of general and secondary metabolism
Nematodes activate their metabolic pathways in response to desiccation shortly after exposure to dehydration . Several transcripts encoding metabolism related genes are differentially expressed by desiccation stress in P. murrayi. These genes include aldehyde dehydrogenase [GenBank: FK670237], trehalose-6-phosphate synthase [GenBank: FK670250], thymidylate synthase [GenBank: FK670263], glycerol kinase [GenBank: FK670258], glycogen synthase [GenBank: FK670264], ATP synthase [GenBank: FK670293], ADP/ATP translocase [GenBank: FK670265], and malate dehydrogenase [GenBank: FK670259]. Interestingly, a transcript encoding a bifunctional glyoxylate cycle protein (malate synthase), a distinct and anaplerotic variant of the tricarboxylic acid cycle, was also found to be highly expressed during desiccation (Table 2).
Nematodes are unique among animals in utilizing the glyoxylate cycle to generate carbohydrates from the beta-oxidation of fatty acids . The glyoxylate pathway, generally found in plants and micro-organisms, is similar to the citrate cycle, but relies on two critical enzymes, malate synthase and isocitrate lyase, to bypass two decarboxylation steps. Interestingly, the anhydrobiotic nematode A. avenae has been reported to use the glyoxylate cycle during induction of anhydrobiosis . One sequence unique to P. murrayi that mapped to a glyoxylate cycle protein [GenBank: FK670269] includes putative homologs of malate synthase (Pm-ms). The abundant expression of malate synthase transcripts in our EST collection and its up-regulation upon desiccation stress may provide experimental support for an active role of glyoxylate cycle proteins during induction of anhydrobiosis by P. murrayi.
Transcriptional regulation and signalling affected by desiccation
Transcriptional regulation and intracellular signalling cascades for nematode stress response in general and secondary metabolism in particular, are poorly understood. A number of desiccation responsive P. murrayi ESTs encode putative signalling molecules or transcription factors. One of the P. murrayi ESTs was most similar to the unc-16 gene of C. elegans, which encodes a c-Jun N-terminal kinase (JNK)-interacting protein , and one of the members of JNK kinase family (jkk-1) [GenBank: FK670239] was DE during desiccation. JNK (also known as stress-activated MAP kinases or SAPK) is a member of the mitogen-activated protein kinases (MAPKs) that regulate cellular responses to a variety of extracellular signals, including desiccation stress [56–58]. Transcripts homologous to genes encoding Zinc finger protein (GenBank: FG619080) were abundantly present in the expressed genome of P. murrayi. Zinc finger proteins are cellular proteins which play a major role in transcriptional regulation by binding with high affinity to specific regions of DNA. In conjunction with leucine zipper domains these proteins may form hetero- or homodimers and activate transcription, either constitutively, or in a regulatory manner, through post-translational modifications in response to external stimuli (although some may also be cell specific or developmentally regulated) . The abundance of genes encoding zinc finger protein in response to desiccation suggests that these proteins may regulate further events in the stress-response cascade of P. murrayi. Three transcripts from P. murrayi were most similar to the C. elegans gene encoding a predicted neurotransmitter gated ion-channel (GenBank: FG620960) protein. Neuronal signal transduction in response to desiccation stress would be required to initiate the coiling response [60, 70] of desiccating P. murrayi.
Stress response genes expressed during anhydrobiosis
Nematodes respond to desiccation stress by synthesizing a conserved set of proteins [23–25]. Our results demonstrate that for P. murrayi desiccation stress can significantly elevate stress related genes encoding trehalose 6-phosphate synthase, late embryogenesis abundant proteins, heat shock proteins, ubiquitin, c-type lectins, chaperone related proteins, and other stress responsive genes (Table 3). Three transcripts encoding trehalose 6-phosphate synthase, which synthesize the storage carbohydrate trehalose, were expressed during desiccation stress. A characteristic feature of anhydrobiotic organisms is their synthesis of high concentrations of non-reducing sugars during the induction of anhydrobiosis [29, 61]. Trehalose protects membranes and proteins from desiccation damage by replacing structural water , and contributes to the formation of an intracellular organic glass  which is thought to stabilize the cell's contents. The up-regulation of trehalose during desiccation stress of P. murrayi could be part of an adaptive response to desiccation (Fig. 3).
Protein aggregation during desiccation is likely to be a major potential hazard for anhydrobiotes; late embryogenesis abundant (LEA) proteins acting as molecular chaperones or molecular shields play an important role in prevention of this aggregation . Transcripts similar to plant LEA related family members [GenBank: FK670249] of C. elegans were up-regulated during desiccation stress of P. murrayi. An LEA group 3 gene Aav-lea-1 was strongly induced in A. avenae during the induction of anhydrobiosis . The C. elegans genome encodes three LEA genes  and silencing of the lea-1 gene by RNA interference (RNAi) caused a marked reduction in desiccation resistance in dauer larvae . We thus assume that LEA proteins contribute to protection and recovery from desiccation stress in anhydrobiotic nematodes.
Molecular chaperones, such as the Hsp70 family and the Hsp60 chaperonin complexes, are commonly perceived as heat shock proteins (Hsps), being up-regulated by stress. The P. murrayi transcripts of 70 [GenBank: FK670245] and 90 kda heat shock protein [GenBank: FK670262] and small heat shock proteins were abundantly expressed following desiccation stress. Heat shock proteins have been implicated in response to desiccation in many nematodes , but they appear to be constitutively expressed in P. murrayi (Additional file 2). Though Hsps may contribute to enhanced stress resistance overall, our results showed no evidence that expression levels of these Hsps were altered by desiccation. It has been shown that other Antarctic organisms constitutively express Hsp70 and Hsp90, showing no or modest up-regulation of this gene in response to thermal stress [65–67]. Characteristically, Hsps are expressed not at the animal's normal habitat condition, but only as part of the organism's stress response. In fact, the expression of these genes is thought to be incompatible with ongoing protein synthesis and the progression of development [68, 69]. It may be that P. murrayi has evolved a mechanism to maintain Hsp function without disrupting normal metabolism and growth that requires synthesis of other proteins. Consistent with the observed constitutive expression, it is possible that desiccation stress did not activate these genes, and the mild desiccation failed to boost Hsp expression. An alternative explanation for the continuous up-regulation of Hsps is that because Antarctic nematodes are frequently, although unpredictably, exposed to a variety of environmental stressors such as desiccation, high pH, extreme osmotic shock, freezing, and anoxia as well as temperature , their survival depends on maintaining continuous expression of molecular chaperones. Furthermore, because of the unpredictability and the potential rapidity of exposure to diverse environmental stresses, the continuous production of these molecular chaperones may be energetically justified.
Unexpectedly, one of the ESTs encoded a protein similar to the type II antifreeze protein (AFP) [GenBank: FK670242] of Atlantic herring (Clupea harengus). This finding is somewhat surprising since antifreeze proteins from phylogenetically divergent organisms (insects, fish, plants, etc.,) generally have very little similarity to each other . The transcript encoding antifreeze protein (AFP) (Pm-afp) was down-regulated during desiccation stress of P. murrayi. It is found that Pm-afp is up-regulated only during freezing and once nematodes are exposed to desiccation other stress response mechanisms are activated by down-regulating AFP (B.N. Adhikari, unpublished work). Many overwintering organisms, including insects, fish, bacteria, fungi, and plants, accumulate AFPs that bind to the faces of ice crystals during freezing and inhibit their growth . While the high degree of AFP similarity between P. murrayi and C. harengus raises the idea that it originated in either of the two organisms via horizontal gene transfer, we suggest that such speculation is premature, particularly given their geographic, ecological, and phylogenetic disparity, but also because such a small fraction of the phylogenetic diversity of nematode genomes has been explored. Further sequencing and characterization are currently underway to better understand the origin and evolution of this gene in P. murrayi.
Oxidative stress genes expressed during anhydrobiosis
Desiccation stress induces the generation of reactive oxygen species (ROS) in nematodes, and therefore it is important for nematodes to have effective ROS-scavenging mechanisms. A number of ESTs encoding proteins which detoxify reactive oxygen species like superoxide dismutase (SOD) [GenBank: FK670240], Ras-related protein [GenBank: FK670252], and glutathione S-transferase (GST) (GenBank: FK670241) were expressed in P. murrayi in response to desiccation (Table 2). The SOD enzymes are a family of metalloenzymes responsible for quenching the potentially deleterious effects of superoxide radicals. There was abundant expression of ESTs similar to C. elegans sod-1 that encode a copper/zinc superoxide dismutase (Table 3). Three transcripts encoding glutathione s-transferase 1 were expressed in the ESTs of P. murrayi and one of them was up-regulated following desiccation stress. GSTs are a diverse super-family of multifunctional proteins that play prominent roles in detoxification metabolism in nematodes . More than a dozen different GSTs have been isolated from C. elegans  and these detoxifying enzymes are reported to be involved in several functions, including xenobiotic detoxification and oxidative stress tolerance [72, 74]. Differential expression and up-regulation of detoxifying enzymes like Pm-sod-1 and Pm-gst-1 suggests that P. murrayi has efficient ROS scavenging mechanisms under desiccation stress.
Membrane and transport-related protein expressed during desiccation
In this study, we identified several ESTs encoding proteins involved in transport facilitation. A total of 15 genes encoding ion and water transporters included ABC transporter proteins (GenBank: GH229101), water channel proteins (GenBank: GH196899), ATPase (GenBank: FG618924) and lipid transfer proteins (LTP) [GenBank: FK670243]. Though the mechanism of coupling ion and water flow through membrane channels is not well studied, there is evidence that up-regulation of such genes is correlated with sensitivity to different types of stress in nematodes . Sixty ABC transporters have been identified and functionally characterized in C. elegans. Members of this protein family are responsible for resistance to heavy metals . In addition to short-term response and regulatory mechanisms, a functional system for reestablishing homeostasis is vital to desiccation tolerance. We found transcripts encoding proteins involved in ion homeostasis, such as vacuolar H+-ATPase and ATP synthase [GenBank: FK670298], showing that P. murrayi has efficient homeostatic pathways. Desiccation tolerance involves changes in the levels and composition of fatty acids of the major glycerolipids in nematodes . Our data showed that gene encoding non-specific LTP was differentially regulated following desiccation, suggesting an active lipid metabolism and desiccation resistance.
Refined gene-specific expression using quantitative real-time PCR
In general, we observed larger changes in gene expression using qRT-PCR, likely reflecting the greater dynamic range of detection and sensitivity of this method for gene expression profiling. Desiccation stress caused significant up-regulation of transcripts encoding stress related genes like Pm-lea (late embryogenesis abundant protein) and Pm-tps (trehalose-6-phosphate synthase) and down-regulation of Pm-afp (antifreeze protein). Similarly, there was significant up-regulation of the ROS scavenging enzyme Pm-gst-1(glutathione s-transferase 1) and metabolism related genes like Pm-gpx (glutathione peroxidise) and Pm-ms (malate synthase), and down-regulation of Pm-afp and Pm-gsy (glycogen synthase). Similar gene expression patterns, except for antifreeze protein, were reported for the free-living mycophagous nematode A. avenae  and insect parasitic nematode Steinernema feltiae  during desiccation acclimation and survival (Fig. 3). Although the observed changes in gene expression are common to many nematodes (except for the antifreeze protein and heat shock proteins), we observed significant variability in the magnitude of transcript abundance in P. murrayi. It is likely that the expression of these genes in response to desiccation stress is influenced by the thermal history of the organism, including seasonal variation experienced in the field, and by phylogenetic constraints.
Nematodes have an important yet poorly understood place in the study of stress response survival, and desiccation tolerance in particular. Nematodes are the most abundant metazoans on earth  and prominent drivers of Antarctic ecosystem functioning . Similar to the Antarctic nematode Panagrolaimus davidi , P. murrayi are desiccation as well as freeze tolerant, which establishes them as a useful model in assessing the structural, physiological, biochemical and genetic aspects of multiple stress tolerance, and the mechanisms by which organisms respond to and survive in extreme environments. To date very few studies have focused on the molecular aspects of desiccation tolerance in Antarctic nematodes and no genomic data is available in existing databases. The progression of Antarctic nematode genomics will be critical to the development of more useful and manipulable models for understanding desiccation tolerance and the nature of extremophiles. It is to these ends that we have initiated this study into the transcriptome of P. murrayi as they respond to a major stress event, in this case desiccation.
In the present study of 2,486 ESTs, a significant portion of transcripts had no known homologue in any nematode or other organisms for which sequence data are currently available in public databases. These molecules are particularly interesting as they may represent genes that are species/lineage specific or unique to Antarctic nematodes. There is considerable scope in exploring such molecules in the future and using a combination of genomic and proteomic approaches could provide valuable information on the survival strategy of extremophiles. Although the number of transcripts analyzed under this study is not large, we have identified and validated several genes differentially expressed during desiccation stress. It is crucial to emphasize that changes in mRNA accumulation may not necessarily correlate with protein/enzyme activity level . However, the expression profiles provide starting points for more in-depth studies on candidate genes using additional genetic approaches. The annotation of transcripts with significant fold change and detection of consistently DE transcripts by SSH and qRT-PCR strongly suggest that these putative genes have an important role in P. murrayi desiccation stress response. In combination with the bioinformatics analysis presented here, further functional genomics analyses are required to generate a more complete picture of the cellular response of a free-living metazoan to extreme abiotic stress and determine the biological roles of these genes within the larger context of its Antarctic ecosystem.
Live nematodes (P. murrayi) from Antarctica were reared in culture in order to produce sufficient quantities for experiments. Nematodes were cultured in sand agar plates containing Bold Modified Basal Freshwater Nutrient Media (BMBFN) (Sigma Aldrich Inc. USA) and stored at 15°C until used. Briefly, 20 ml/L BMBFN solution was mixed with 15 g/L Bacto-Agar in deionised water and autoclaved. The autoclaved Agar media was poured into Petri dishes and Standard Ottawa Sand (EMD Chemicals Inc., Gibbstown, NJ) was spread over each plate before the media solidified. About 5–10 nematodes were released into each plate, and a thin layer of water was maintained on the surface of the media for the duration of culture. The plates were incubated at 26°C for one to two weeks followed by 16°C for 3–5 weeks. In the laboratory culture at 16°C, nematodes completed their life cycle in about 6–8 weeks. Nematodes were also cultured in Nematode growth medium (NGM) plates with Escherichia coli OP50. About 15–20 nematodes were placed in one NGM plate and incubated at 16°C for about 5–8 weeks. Nematodes were sub-cultured every 3–4 weeks by chunking the agar with a sterile scalpel and moving a chunk of agar from an old plate to a fresh plate containing a lawn of E. coli OP50. To harvest the nematodes, the Agar media was cut into small pieces and poured into #40 standard testing sieve (Advantech Manufacturing, New Berlin, WI). Agar was removed by sieving and nematodes were washed twice with deionized water before collection by centrifugation.
Nematodes (about 3,000–5,000) were desiccated in glass desiccation chambers. The relative humidity (RH) was controlled by using saturated salt solutions as described by Winston & Bates  with some modifications. The requisite RH was maintained at 23°C ± 1°C in the desiccation chamber for 3 days prior to the addition of nematodes for equilibration. Relative humidity was maintained as 100% RH with distilled water vapour; 97% with K2SO4, 87% RH with HCl and 75% with NaCl. The treatments were an initial desiccation at 97% RH for 3 days followed by exposure to 87% for 2 days. All experiments were repeated under identical conditions using mixed stage populations of P. murrayi.
The desiccated and control nematodes were transferred to 10 volumes of Trizol Reagent (Molecular Research Center Inc., Cincinnati, OH) and exposed to freeze thaw cycles using liquid nitrogen and a 37°C water bath. The suspension was ground using mortar and pestle and vortexed. RNA was phase separated using chloroform, precipitated by isopropanol and pelleted. Total RNA was quantified and quality checked by spectrophotometer (NanoDrop, ND-1000 Spectrophotometer) and running agarose gel with RNA Century™ Plus Marker (Ambion Inc., Austin, TX).
cDNA Library Construction and Sequencing
Total RNA was extracted from the gradually desiccated nematodes and used as starting material for complementary deoxyribonucleic acid (cDNA) library construction. A full-length cDNA library was prepared using Creator™ SMART™ cDNA library construction kit (Clontech, Mountain View, CA) according to the manufacturer's protocol. Briefly, 2 μg of total RNA was reverse transcribed to cDNA using SMART™ IV Oligonucleotide and PowerScript Reverse Transcriptase and amplified using long distance (LD) PCR. The cycling parameters were 95°C, 20 s; two 20 cycles of 95°C, 5 s; and 68°C, 6 min. cDNA was analyzed using agarose gel electrophoresis, proteinase K digested and restriction digested by SfiI. cDNA was size fractionated using CHROMA SPIN-400 columns and analyzed via agarose gel electrophoresis. The size fractionated cDNA was directionally cloned into pDNR-LIB vector and incubated at 16°C for 15 h. Transformation was done at 37°C overnight, and colonies were picked and grown for 18 h. Template DNA was extracted using QIAprep Spin Miniprep Kit (Qiagen Inc., Valencia, CA), amplified, cleaned using ExoSAP-IT® (USB Corporation, Cleveland, Ohio) and sequenced using an ABI PRISM 377 automated DNA sequencer (Applied Biosystems, CA) at the DNA Sequencing Center (DNASC) at Brigham Young University (BYU). Sequencing reactions were performed with 2 μl of ABI PRISM BigDye Terminators v3.0 (Applied Biosystems) and 3.5 pmol of primer in 12 μl reaction volumes, followed by a sequencing reaction clean up to remove residual dye and enzyme. Sequencing was with M13 forward (GTAAAACGACGGCCAG) and reverse (CAGGAAACAGCTATGAC) primers for 25 cycles of 96°C, 10 s; 50°C, 5 s; and 60°C, 4 min. The sequenced products were sorted and analyzed with Sequencher™ 4.8 (Gene Codes Corporation, Ann Arbor, MI).
Subtractive hybridization was performed with the PCR-Select™ cDNA subtraction kit (Clontech, Palo Alto, CA). Random hexamer primers (Clontech) were used to convert mRNA that had been extracted from desiccated and control nematodes into cDNA. The cDNA under investigation was designated tester, and the reference cDNA was designated driver. In the forward run cDNA from desiccated samples was designated as tester, in a reverse run control nematode cDNA was used as tester to investigate both up- and down-regulation of mRNA molecules. Second strands of driver and tester DNA were synthesized using T4 DNA polymerase and restriction digested by RsaI to create blunt-ended double stranded cDNA for subtraction. The tester cDNAs were subdivided into two portions and each group ligated to different adaptor sequence and subjected to two levels of subtractive hybridizations as described in the manufacturer's protocol. The DE cDNAs were selectively amplified by two PCR reactions using 50× Advantage cDNA Polymerase mix. The first PCR amplified double stranded cDNA with different adaptor sequences and the second PCR further reduced background and enriched for DE sequences. The PCR-amplified cDNA fragments generated by SSH were then ligated into the TOPO® TA Cloning® vector (Invitrogen Corporation, USA) and sequenced using ABI 3730xl DNA Analyzer (Applied Biosystems, CA) at the BYU DNASC.
Sequence analysis and EST clustering
Base calling was performed using PHRED software (versions 0.000925.c) [82, 83] with the quality cut-off set at PHRED 20. Raw sequences were then imported into the Vector NTI Advance™ 10 software (Invitrogen Corporation) and subjected to trimming of vector sequences and 5' adapter sequences (for subtractive hybridization) using default settings. Sequences with less than 100 quality bases (PHRED 20 or better) after trimming were discarded. Sequences having polyA tails of 100 bases or more were eliminated from analysis. EST sequences representing contamination from bacterial, yeast or fungal sources were identified using BLASTN algorithm [34, 35] and removed from further analyses. ESTs were aligned and assembled into contigs using CAP3 software  when the criterion of a minimum identity of 95% over 50 bp was met. When an EST could not be assembled with others in a contig, it remained as a "singleton". The contigs and the singletons should thus correspond to sequences of unique genes. The consensus sequences of the contigs and the sequences of the singlets were compared to the sequences in GenBank's non-redundant (nr) and Uniprot database using the TBLASTX and the BLASTX [38, 39] algorithms and C. elegans Wormpep 190 database  using the BLASTX algorithm . The cut-off for sequence similarity was E-value < 10-5 for all analyses.
Functional analysis and pathway assignments
Gene ontology (GO) term annotation and function-based analysis  of unique sequences were performed using Blast2go (V 1.6.2) , a sequence-based tool to assign GO terms, extracting them for each BLAST hit obtained by mapping to extant annotation associations. GO terms for each of the three main categories (biological process, molecular function, and cellular component) was obtained from sequence similarity using the application default parameters. From these annotations, pie charts were made using 2nd level GO terms based on biological process, molecular function, and cellular component. Pathway assignments were carried out according to Kyoto encyclopedia of genes and genomes (KEGG) mapping . Enzyme commission (EC) numbers were assigned to unique sequences that had BLASTX scores with a cut-off value of E = 10-5 or less upon searching protein databases. The sequences were mapped to KEGG biochemical pathways according to the EC distribution in the pathway database.
All the primers used in quantitative real-time PCR (qRT-PCR) were designed with IDT SciTools (Integrated DNA Technologies, Coralville, IA) by aligning EST sequences with similar sequences from NCBI. All the primers used in this study were synthesized by Operon (Operon Biotechnologies Inc., Huntsville, AL). The sequences of the primers and product sizes are listed in Additional file 3.
Quantitative Real-time PCR
Total RNA extracted from dehydrated and control P. murrayi nematodes was reverse transcribed using ImProm-II™ reverse transcriptase (Promega corporation, Madison, WI) and subjected to qRT-PCR analysis using Light cycler 480 SYBER Green I mastermix and gene specific primers in a Light cycler 480 RT-PCR system (Roche Applied Science, Mannheim, Germany) equipped with light cycler 480 software. High-resolution gel electrophoresis was used to verify that the qRT-PCR amplification product from each examined gene was a single-band product. Thermal cycling was performed in accordance with the manufacturer's instructions for a total of 40 cycles at an annealing temperature of 58°C for each primer pair. Quantitative RT-PCR analysis was performed with Lightcycler 480 software, the threshold cycle was automatically calculated by the second-derivative maximum method, and the copy number of the specific mRNA in the experimental samples was calculated by extrapolation from the gene-specific standard curve.
The copy number of specific cDNA molecules present in the samples was determined by absolute quantification method of qPCR analysis . A range of six dilutions (107–102 copies) of the cDNA was made and a gene-specific external standard curve was generated by using cDNA standards that were run simultaneously with the experimental samples. Change in target gene expression was calculated using equation 2-ΔΔCT . The fold change in the target gene, normalized to 18S rRNA (Pm-18s) and relative to the expression of control, was calculated for each sample. A gene with a relative abundance of one is equal to the abundance of 18S rRNA in the same sample. An F-test at a significance level of P < 0.05 was used to compare the ratio of the mean gene expression of desiccated samples with that of the control. To minimize mRNA quantification errors, genomic DNA contamination biases and to correct for inter-sample variations, we used 18s ribosomal RNAs (rRNAs) of P. murrayi as an internal control.
This work was supported by the National Science Foundation McMurdo Long Term Ecological Research Project (OPP #98-10219) to DHW, the United States Department of Agriculture CSREES NRI (2005-00903) to BJA, and a Brigham Young University Mentored Environment Grant to BJA. We thank C-Y Lin, J. Griffitts, and G. F. Burton for their critical advice during the planning and execution of these experiments. We are grateful to J. A. Udall for bioinformatics assistance, and three anonymous reviewers whose constructive criticism significantly improved the paper. Funding for the Open Access publication charges for this article was provided by the Brigham Young University College of Life Sciences and Department of Biology.
- Priscu JC: Ecosystem Dynamics in Polar Desert: The McMurdo Dry Valleys of Antarctica. 1998, Washington DC: American Geophysical UnionView ArticleGoogle Scholar
- Campbell IB, Claridge GGC, Campbell DI, Balks MR: The soil environment of the McMurdo Dry Valleys, Antarctica. Ecosystem Dynamics in a Polar Desert: the McMurdo Dry Valleys, Antarctica. Edited by: Priscu JC. 1998, Washington DC: American Geophysical Union, 297-322.Google Scholar
- Wall Freckman D, Virginia RA: Soil biodiversity and community structure in the McMurdo Dry Valleys, Antarctica. Ecosystem Dynamics in a Polar Desert: the McMurdo Dry Valleys Antarctica. Edited by: Priscu JC. 1998, Washington DC: American Geophysical Union, 323-325.Google Scholar
- Treonis AM, Wall DH, Virginia RA: Invertebrate biodiversity in Antarctic Dry Valley soils and sediments. Ecosystems. 1999, 2: 482-92. 10.1007/s100219900096.View ArticleGoogle Scholar
- Treonis AM, Wall DH, Virginia RA: The use of anhydrobiosis by soil nematodes in the Antarctic Dry Valleys. Funct Ecol. 2000, 14: 460-467. 10.1046/j.1365-2435.2000.00442.x.View ArticleGoogle Scholar
- Crowe JH, Madin KAC: Anhydrobiosis in tardigrades and nematodes. Trans Am Microsc Soc. 1974, 93: 513-524. 10.2307/3225155.View ArticleGoogle Scholar
- Crowe JH, Madin KAC: Anhydrobiosis in nematodes: Evaporative water loss and survival. J Exp Zool. 1975, 193: 323-34. 10.1002/jez.1401930308.View ArticleGoogle Scholar
- Pickup J, Rothery P: Water-loss and anhydrobiotic survival in nematodes of Antarctic fellfields. Oikos. 1991, 61: 379-88. 10.2307/3545245.View ArticleGoogle Scholar
- Wharton DA, Block DW: Freezing tolerance of some Antarctic nematodes. Funct Ecol. 1993, 7: 578-584. 10.2307/2390134.View ArticleGoogle Scholar
- Adams BJ, Bardgett RD, Ayres E, Wall DH, Aislabie J, Bamforth S, Bargafli R, Cary C, Cavacini P, Connell L, Convey P, Fell JW, Frati F, Hogg ID, Newsham KK, O'Donnell A, Russel N, Seppelt RD, Stevens MI: Diversity and distribution of Victoria Land biota. Soil Biol Biochem. 2006, 38: 3003-3018. 10.1016/j.soilbio.2006.04.030.View ArticleGoogle Scholar
- Andrassy I: Nematodes in the sixth continent. J Nemat Morph Syst. 1998, 1: 107-187.Google Scholar
- Andrassy I: On the male of the Antarctic nematode species, Plectus murrayi Yeates, 1970. J Nemat Morph Syst. 2008, 11: 87-89.Google Scholar
- Freckman DW, Virginia RA: Low-diversity Antarctic soil nematode communities: distribution and response to disturbance. Ecology. 1997, 78: 363-369.View ArticleGoogle Scholar
- Wall DH, Virginia RA: Controls on soil biodiversity: insights from extreme environments. Appl Soil Ecol. 1999, 13: 137-150. 10.1016/S0929-1393(99)00029-3.View ArticleGoogle Scholar
- Overhoff A, Freckman DW, Virginia RA: Life cycle of the microbivorous Antarctic Dry Valley nematode Scottnema lindsayae (Timm 1971). Polar Biol. 1993, 13: 151-156. 10.1007/BF00238924.View ArticleGoogle Scholar
- Powers LE, Ho M, Freckman DW, Virginia RA: Distribution, community structure, and microhabitat of soil invertebrates along an elevational gradient in Taylor Valley, Anarctica. Arct Alpine Res. 1998, 30: 133-141. 10.2307/1552128.View ArticleGoogle Scholar
- Liu QZ, Piggott SJ, Solomon A, Glazer I: Physiological and biochemical changes in nematodes from the Genus Heterorhabditi s Following desiccation. Phytoparasitica. 2000, 30: 253-261.Google Scholar
- Gal T, Glazer I, Sherman A, Koltai H: Protein interaction of Nucleosome assembly Protein-1 and casein kinase 2 during desiccation response in the insect-killing nematode Steinernema feltiae IS-6. J Parasitol. 2005, 91: 691-693. 10.1645/GE-402R.View ArticlePubMedGoogle Scholar
- Browne JA, Dolan KM, Tyson T, Goyal K, Tunnacliffe A, Burnell AM: Dehydration-specific induction of hydrophilic protein genes in the anhydrobiotic nematode Aphelenchus avenae. Eukar Cell. 2004, 3: 966-975. 10.1128/EC.3.4.966-975.2004.View ArticleGoogle Scholar
- Crowe LM: Lessons from nature: the role of sugars in anhydrobiosis. Comp Biochem Physiol A. 2002, 131: 505-513. 10.1016/S1095-6433(01)00503-7.View ArticleGoogle Scholar
- Wharton DA: The environmental physiology of Antarctic terrestrial nematodes: a review. J Comp Physiol B. 2003, 173: 621-628. 10.1007/s00360-003-0378-0.View ArticlePubMedGoogle Scholar
- Solomon A, Solomon R, Paperna I, Glazer I: Desiccation stress of entomopathogenic nematodes induces the accumulation of a novel heat-stable protein. Parasitol. 2000, 121: 409-416. 10.1017/S0031182099006563.View ArticleGoogle Scholar
- Browne J, Tunnacliffe A, Burnell AM: Plant desiccation gene found in a nematode. Nature. 2002, 416: 38-10.1038/416038a.View ArticlePubMedGoogle Scholar
- Oliver AE, Hincha DH, Crowe JH: Looking beyond sugars: the role of amphiphilic solutes in preventing adventitious reactions in anhydrobiotes at low water contents. Comp Biochem Physiol A. 2002, 131: 515-525. 10.1016/S1095-6433(01)00514-1.View ArticleGoogle Scholar
- Gal TZ, Glazer I, Koltai H: Differential gene expression during desiccation stress in the insect-killing nematode Steinernema feltiae IS-6. J Parasitol. 2003, 89: 761-766. 10.1645/GE-3105.View ArticlePubMedGoogle Scholar
- Chen SB, Glazer I, Gollop N, Cash P, Argo E, Innes A, Stewart E, Davidson I, Wilson MJ: Proteomic analysis of the entomopathogenic nematode Steinernema feltiae IS-6 IJs under evaporative and osmotic stresses. Mol Biochem Parasitol. 2006, 145: 195-204. 10.1016/j.molbiopara.2005.10.003.View ArticlePubMedGoogle Scholar
- Tyson T, Reardon W, Browne JA, Burnell AM: Gene induction by desiccation stress in the entomopathogenic nematode Steinernema carpocapsae revelas parallels with drought tolerance mechanisms in plants. Int J Parasitol. 2007, 37: 763-776. 10.1016/j.ijpara.2006.12.015.View ArticlePubMedGoogle Scholar
- Madin KAC, Crowe JH: Anhydrobiosis in nematodes: carbohydrate and lipid metabolism during dehydration. J Exp Zool. 1975, 193: 335-342. 10.1002/jez.1401930309.View ArticleGoogle Scholar
- Hoekstra AF, Golorina EA, Buitink J: Mechanisms of plant desiccation tolerance. Trends Plant Sci. 2001, 6: 431-438. 10.1016/S1360-1385(01)02052-0.View ArticlePubMedGoogle Scholar
- Tunnacliffe A, Lapinski J: Resurrecting Van Leeuwenhoek's rotifers: a reappraisal of the role of disaccharides in anhydrobiosis. Philos Trans Roy Soc London B. 2003, 358: 1755-1771. 10.1098/rstb.2002.1214.View ArticleGoogle Scholar
- Goyal K, Walton LJ, Browne JA, Burnell AM, Tunnacliffe A: Molecular Anhydrobiology: identifying Molecules Implicated in Invertebrate Anhydrobiosis. Integr Comp Biol. 2005, 45: 702-709. 10.1093/icb/45.5.702.View ArticlePubMedGoogle Scholar
- Somvanshi VS, Koltai H, Glazer I: Expression of different desiccation tolerance related genes in various species of entomopathogenic nematodes. Mol Biochem Parasitol. 2008, 158: 65-71. 10.1016/j.molbiopara.2007.11.012.View ArticlePubMedGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.View ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- Rogers A, Antoshechkin I, Bieri T, Blasiar D, Bastiani C, Canaran P, Chan J, Chen WJ, Davis P, Fernandes J, Fiedler TJ, Han M, Harris TW, Kishore R, Lee R, McKay S, Muller HM, Nakamura C, Ozersky P, Petcherski A, Schindelman G, Schwarz EM, Spooner W, Tuli MA, Van Auken K, Wang D, Wang X, Williams G, Yook K, Durbin R, Stein LD, Spieth J, Sternberg PW: WormBase 2007. Nucleic Acids Res. 2008, D612-617. 36 DatabaseGoogle Scholar
- Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL: GenBank. Nucleic Acids Res. 2007, 35: D21-25. 10.1093/nar/gkl986.PubMed CentralView ArticlePubMedGoogle Scholar
- Wharton DA, Barrett J, Goodall G, Marshall CJ, Ramlov H: Ice-active proteins from the Antarctic nematode Palangrolaimus davidi. Cryobiology. 2005, 51: 198-207. 10.1016/j.cryobiol.2005.07.001.View ArticlePubMedGoogle Scholar
- Chen L, DeVries AL, Cheng CC: Evolution of antifreeze glycoprotein gene from a trypsinogen gene in Antarctic notothenioid fish. Proc Natl Acad Sci USA. 1997, 94: 3811-3816. 10.1073/pnas.94.8.3811.PubMed CentralView ArticlePubMedGoogle Scholar
- Graham LA, Liou YC, Walker VK, Davies PL: Hyperactive antifreeze protein from beetles. Nature. 388: 727-728. 10.1038/41908.Google Scholar
- Griffith M, Yaish M: Antifreeze proteins in overwintering plants: a tale of two activities. Trends Plant Sci. 2004, 9: 399-405. 10.1016/j.tplants.2004.06.007.View ArticlePubMedGoogle Scholar
- Duman JG, Olsen TM: Thermal hysteresis protein activity in bacteria, fungi and phylogenetically diverse plants. Cryobiology. 1993, 30: 322-328. 10.1006/cryo.1993.1031.View ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry MJ, Davis AP, Dilinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.PubMed CentralView ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S: KEGG: Kyoto Encyclopaedia of Genes and Genomes. Nucleic Acids Res. 2000, 28: 27-30. 10.1093/nar/28.1.27.PubMed CentralView ArticlePubMedGoogle Scholar
- Bustin SA: Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. J Mol Endocrinol. 2000, 25: 169-193. 10.1677/jme.0.0250169.View ArticlePubMedGoogle Scholar
- Wasmuth J, Schmid R, Hedley A, Blaxter M: On the extent and origins of genic novelty in the Phylum Nematoda. PLoS Negl Trop Dis. 2008, 2: e258-10.1371/journal.pntd.0000258.PubMed CentralView ArticlePubMedGoogle Scholar
- Nagaraj SH, Gasser RB, Ranganathan S: A hitchhiker's guide to expressed sequence tag (EST) analysis. Brief Bioinform. 2007, 8: 6-21. 10.1093/bib/bbl015.View ArticlePubMedGoogle Scholar
- Hashmi S, Tawe W, Lustigman S: Caenorhabditis elegans and study of gene functions in parasites. Trends Parasitol. 2001, 17: 387-393. 10.1016/S1471-4922(01)01986-9.View ArticlePubMedGoogle Scholar
- Lamitina T, Huang CG, Strange K: Genome-wide RNAi screening identifies protein damage as a regulator of osmoprotective gene expression. Proc Natl Acad Sci USA. 2006, 103: 12173-12178. 10.1073/pnas.0602987103.PubMed CentralView ArticlePubMedGoogle Scholar
- Audic S, Claverie JM: The significance of digital gene expression profiles. Genome Res. 1997, 7: 986-995.PubMedGoogle Scholar
- McIntosh KB, Bonham-Smith PC: Ribosomal protein gene regulation: what about plants?. Can J Bot. 2006, 84: 342-362. 10.1139/B06-014.View ArticleGoogle Scholar
- Mauro VP, Edelman GM: The ribosome filter hypothesis. Proc Natl Acad Sci USA. 2002, 99: 12031-12036. 10.1073/pnas.192442499.PubMed CentralView ArticlePubMedGoogle Scholar
- Feder ME, Hofmann GE: Heat-shock proteins, molecular chaperones, and the stress response: evolutionary and ecological physiology. Annu Rev Physiol. 1999, 61: 243-282. 10.1146/annurev.physiol.61.1.243.View ArticlePubMedGoogle Scholar
- Barrett J, Wright DJ: Intermediary metabolism. The Physiology and Biochemistry of Free-living and Plant-Parasitic Nematodes. Edited by: Perry RN, Wright DJ. 1998, Oxford UK: CABI Publishing, 331-353.Google Scholar
- Madin KAC, Loomis SH, Crowe JH: Anhydrobiosis in nematodes-control of carbon flow through the glyoxylate cycle. J Exp Zool. 1985, 234: 341-350. 10.1002/jez.1402340303.View ArticleGoogle Scholar
- Byrd DT, Kawasaki M, Walcoff M, Hisamoto N, Matsumoto K, Jin YS: UNC-16, a JNK-signaling scaffold protein, regulates vesicle transport in Caenorhabditis elegans. Neuron. 2001, 32: 787-800. 10.1016/S0896-6273(01)00532-3.View ArticlePubMedGoogle Scholar
- Davis RJ: Signal transduction by the JNK group of MAP Kinases. Cell. 2000, 103: 239-252. 10.1016/S0092-8674(00)00116-1.View ArticlePubMedGoogle Scholar
- Huang ZB, Tunnacliffe A: Response of human cells to desiccation: comparison with hyperosmotic stress response. J Physiol. 2004, 58: 181-191. 10.1113/jphysiol.2004.065540.View ArticleGoogle Scholar
- Nakagami H, Pitzschke A, Hirt H: Emerging MAP kinase pathways in plant stress signalling. Trends Plant Sci. 2005, 10: 339-346. 10.1016/j.tplants.2005.05.009.View ArticlePubMedGoogle Scholar
- Hurst LD: Embryonic growth and the evolution of the mammalian Y chromosome I. The Y as an attractor for selfish growth factors. Heredity. 1994, 73: 223-232. 10.1038/hdy.1994.127.View ArticlePubMedGoogle Scholar
- Freckman DW, Virginia RA: Extraction of nematodes from Dry Valley Antarctic soils. Polar Biol. 1993, 13: 483-487. 10.1007/BF00233139.View ArticleGoogle Scholar
- Crowe JH, Hoekstra FA, Crowe LM: Anhydrobiosis. Annu Rev Physiol. 1992, 54: 579-599. 10.1146/annurev.ph.54.030192.003051.View ArticlePubMedGoogle Scholar
- Crowe JH, Carpenter JF, Crowe LM: The role of vitrification in anhydrobiosis. Annu Rev Physiol. 1998, 60: 73-103. 10.1146/annurev.physiol.60.1.73.View ArticlePubMedGoogle Scholar
- Goyal K, Walton LJ, Tunnacliffe A: LEA proteins prevent protein aggregation due to water stress. Biochem J. 2005, 388: 151-157. 10.1042/BJ20041931.PubMed CentralView ArticlePubMedGoogle Scholar
- Gal T, Glazer I, Koltai H: A LEA3 family member is involved in survival of C. elegans during exposure to dehydration stress. FEBS Lett. 2004, 577: 21-26. 10.1016/j.febslet.2004.09.049.View ArticlePubMedGoogle Scholar
- Buckley BA, Place SP, Hofmann GE: Regulation of heatshock genes in isolated hepatocytes from an Antarctic fish, Trematomus bernacchii. J Exp Biol. 2004, 207: 3649-3656. 10.1242/jeb.01219.View ArticlePubMedGoogle Scholar
- La Terza A, Papa G, Miceli C, Luporini P: Divergencebetween two Antarctic species of the ciliate Euplotes, E. focardii and E. nobilii, in the expression of heat-shock protein 70 genes. Mol Ecol. 2001, 10: 1061-1067. 10.1046/j.1365-294X.2001.01242.x.View ArticlePubMedGoogle Scholar
- Rinehart JP, Hayward SAL, Elnitsky MA, Sandro LH, Lee RE, Denlinger DL: Continuous up-regulation of heat shock proteins in larvae, but not adults, of a polar insect. Proc Natl Acad Sci USA. 2006, 103: 14223-14227. 10.1073/pnas.0606840103.PubMed CentralView ArticlePubMedGoogle Scholar
- Feder JH, Rossi JM, Solomon J, Solomon N, Lindquist : The consequences of expressing hsp70 in Drosophila cells at normal temperatures. Gene Dev. 1992, 6: 1402-1413. 10.1101/gad.6.8.1402.View ArticlePubMedGoogle Scholar
- Krebs RA, Feder ME: Deleterious consequences of Hsp70 overexpression in Drosophila melanogaster larvae. Cell Stress Chaperones. 2: 60-71. 10.1379/1466-1268(1997)002<0060:DCOHOI>2.3.CO;2.Google Scholar
- Treonis AM, Wall DH: Soil nematodes and desiccation survival in the extreme arid environment of the Antarctic Dry Valleys. Integr Comp Biol. 2005, 45: 741-750. 10.1093/icb/45.5.741.View ArticlePubMedGoogle Scholar
- Duman JG: Antifreeze and ice nucleator proteins in terrestrial arthropods. Annu Rev Physiol. 2001, 63: 327-357. 10.1146/annurev.physiol.63.1.327.View ArticlePubMedGoogle Scholar
- Lindblom TH, Dodd AK: Xenobiotic detoxification in the nematode Caenorhabditis elegans. J Exp Zool. 2006, 305: 720-730. 10.1002/jez.a.324.View ArticleGoogle Scholar
- Van Rossum AJ, Jefferies JR, Young CJ, Barrett J, Tait A, Brophy PM: Glutathione transferase (GST) functional genomics: the role of Caenorhabditis elegans in investigating GST expression in parasitic nematodes. Chem-Biol Interact. 2001, 133: 274-277.Google Scholar
- Leiers B, Kampkötter A, Grevelding CG, Link CD, Johnson TE, Henkle-Dührsen K: A stress-responsive glutathione S-transferase confers resistance to oxidative stress in Caenorhabditis elegans. Free Radical Biol Med. 2003, 34: 1405-1415. 10.1016/S0891-5849(03)00102-3.View ArticleGoogle Scholar
- Sheps JA, Ralph S, Zhao Z, Baillie DL, Ling V: The ABC transporter gene family of Caenorhabditis elegans has implications for the evolutionary dynamics of multidrug resistance in eukaryotes. Genome Biol. 2004, 5: R15-10.1186/gb-2004-5-3-r15.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhao Z, Sheps JA, Ling V, Fang LL, Baillie DL: Expression analysis of ABC transporters reveals differential functions of tandemly duplicated genes in Caenorhabditis elegans. J Mol Biol. 2004, 344: 409-417. 10.1016/j.jmb.2004.09.052.View ArticlePubMedGoogle Scholar
- Bongers T: The Maturity index, the evolution of nematode life history traits, adaptive radiation and cp-scaling. Plant Soil. 1999, 212: 13-22. 10.1023/A:1004571900425.View ArticleGoogle Scholar
- Barrett JE, Virginia RA, Wall DH, Adams BJ: Decline in a dominant invertebrate species contributes to altered carbon cycling in a low-diversity soil ecosystem. Global Change Biol. 2008, 14: 1-11.Google Scholar
- Smith T, Wharton DA, Marshall C: Cold tolerance of an Antarctic nematode that survives intracellular freezing: comparisons with other nematode species. J Comp Physiol B. 2008, 178: 93-100. 10.1007/s00360-007-0202-3.View ArticlePubMedGoogle Scholar
- Gygi SP, Rochon Y, Franza BR, Aebersold R: Correlation between protein and mRNA abundance in yeast. Mol Cell Biol. 1999, 19: 1720-1730.PubMed CentralView ArticlePubMedGoogle Scholar
- Winston PW, Bates PS: Saturated salt solutions for the control of humidity in biological research. Ecology. 1960, 41: 232-237. 10.2307/1931961.View ArticleGoogle Scholar
- Ewing B, Green P: Basecalling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8: 186-194.View ArticlePubMedGoogle Scholar
- Ewing B, Hillier L, Wendl M, Green P: Basecalling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8: 175-185.View ArticlePubMedGoogle Scholar
- Huang X, Madan A: CAP3: a DNA sequence assembly program. Genome Res. 1999, 9: 868-877. 10.1101/gr.9.9.868.PubMed CentralView ArticlePubMedGoogle Scholar
- Gene Ontology Consortium: Gene Ontology: tool for the unification of biology. Nature Genet. 2000, 25: 25-29. 10.1038/75556.View ArticleGoogle Scholar
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.View ArticlePubMedGoogle Scholar
- Livak KJ, Smittgen TD: Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2-ΔΔCt Method. Methods-Elseveir Science (USA). 2001, 25: 402-408.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.