Expression profiling and cross-species RNA interference (RNAi) of desiccation-induced transcripts in the anhydrobiotic nematode Aphelenchus avenae
- Wesley Reardon1,
- Sohini Chakrabortee†2,
- Tiago Campos Pereira†2, 3,
- Trevor Tyson†1,
- Matthew C Banton2,
- Katharine M Dolan1, 4,
- Bridget A Culleton1,
- Michael J Wise5,
- Ann M Burnell1Email author and
- Alan Tunnacliffe2Email author
© Reardon et al; licensee BioMed Central Ltd. 2010
Received: 5 August 2009
Accepted: 19 January 2010
Published: 19 January 2010
Some organisms can survive extreme desiccation by entering a state of suspended animation known as anhydrobiosis. The free-living mycophagous nematode Aphelenchus avenae can be induced to enter anhydrobiosis by pre-exposure to moderate reductions in relative humidity (RH) prior to extreme desiccation. This preconditioning phase is thought to allow modification of the transcriptome by activation of genes required for desiccation tolerance.
To identify such genes, a panel of expressed sequence tags (ESTs) enriched for sequences upregulated in A. avenae during preconditioning was created. A subset of 30 genes with significant matches in databases, together with a number of apparently novel sequences, were chosen for further study. Several of the recognisable genes are associated with water stress, encoding, for example, two new hydrophilic proteins related to the late embryogenesis abundant (LEA) protein family. Expression studies confirmed EST panel members to be upregulated by evaporative water loss, and the majority of genes was also induced by osmotic stress and cold, but rather fewer by heat. We attempted to use RNA interference (RNAi) to demonstrate the importance of this gene set for anhydrobiosis, but found A. avenae to be recalcitrant with the techniques used. Instead, therefore, we developed a cross-species RNAi procedure using A. avenae sequences in another anhydrobiotic nematode, Panagrolaimus superbus, which is amenable to gene silencing. Of 20 A. avenae ESTs screened, a significant reduction in survival of desiccation in treated P. superbus populations was observed with two sequences, one of which was novel, while the other encoded a glutathione peroxidase. To confirm a role for glutathione peroxidases in anhydrobiosis, RNAi with cognate sequences from P. superbus was performed and was also shown to reduce desiccation tolerance in this species.
This study has identified and characterised the expression profiles of members of the anhydrobiotic gene set in A. avenae. It also demonstrates the potential of RNAi for the analysis of anhydrobiosis and provides the first genetic data to underline the importance of effective antioxidant systems in metazoan desiccation tolerance.
Water is essential for life, but some organisms have the ability to survive extreme desiccation by entering a state of suspended animation in a process known as anhydrobiosis ("life without water") . In the dry state, metabolism and life processes come to a halt. Anhydrobiotic organisms are capable of surviving in this state for an indefinite period of time after which, when environmental conditions become favourable, the organism can resume normal metabolic activity [2, 3]. These organisms are widespread throughout nature and include yeasts, bdelloid rotifers, tardigrades, plants and nematodes. For example, the resurrection plant Craterostigma plantagineum can recover from loss of most of its cellular fluid within 24 h of contact with water [4, 5]. Many nematode species are also able to enter anhydrobiosis, for example, Aphelenchus avenae and Panagrolaimus superbus[6, 7]. Some steinernemitid species also exhibit partial desiccation tolerance [8, 9].
Anhydrobiosis is often associated with the accumulation of disaccharides and other oligosaccharides, chiefly trehalose in animals and yeasts, and sucrose in plants [1, 4, 10]. These sugars have been suggested to play a role in protecting membranes and proteins by replacing structural water or by forming a stabilising intracellular glass that would inhibit membrane fusion and protein denaturation; they might also be involved in maintaining membrane fluidity [11, 12]. However, it has been shown in A. avenae, one of the best-characterised anhydrobiotic nematodes, that while it accumulates trehalose prior to anhydrobiosis, probably at least partly due to upregulation of trehalose synthase genes , this alone is not sufficient . In other desiccation-tolerant invertebrates, such as bdelloid rotifers [15–17] and some tardigrades , disaccharide accumulation is not detectable. In the yeast Saccharomyces cerevisiae, trehalose biosynthesis can be abolished by mutation while retaining 50-100% capacity to survive desiccation, depending on the strain . This suggests that while trehalose might play a role in anhydrobiosis, at least in some species, there must be other important biological processes involved . These processes might include changes in primary metabolism; alterations to cell membranes; the accumulation of compatible solutes and hydrophilic proteins; and activation of antioxidant and molecular chaperone systems.
One group of proteins receiving increasing attention is the LEA (late embryogenesis abundant) proteins. These were first associated with desiccation tolerance in plant seeds and resurrection plants (reviewed in ), but recently homologues of Group 3 LEA proteins have been linked to anhydrobiosis in certain nematodes [22–25] and other invertebrates [17, 26–29]. Group 3 LEA proteins are characterized by extreme hydrophilicity and 11-amino acid periodicity; they have relatively little secondary structure in the hydrated state, but can become more structured on drying [21, 30]. A nematode Group 3 LEA protein from A. avenae, AavLEA1, has been shown to function as a protein anti-aggregant both in vitro and in vivo[31, 32]. It has therefore been postulated that many LEA proteins possess a molecular shield activity whereby they reduce inappropriate interactions between other proteins by an electrosteric mechanism [31, 33]. LEA proteins also exhibit other activities consistent with a role in protection of cells against abiotic stress, including membrane protection, ion binding, antioxidant functions, hydration buffering and nucleic acid binding (references in  and ).
In this study we have constructed an EST (expressed sequence tag) library of genes upregulated on desiccation in the free-living mycophagous nematode A. avenae and have investigated their expression profiles under various stress conditions. We have identified a new LEA protein gene, as well as other stress response genes, and have used a cross-species gene silencing approach to assess the importance of a subset of the EST panel for anhydrobiosis. Our findings indicate that glutathione peroxidases represent a key class of enzymes required for desiccation tolerance in the nematodes tested and emphasise the need for redox balancing in anhydrobiosis.
An EST dataset from A. avenae enriched for dehydration-induced sequences
ESTs, with similarity to known sequences, identified by enrichment for dehydration-induced gene expression in A. avenae.
Closest BLAST hit (organism)a
Homologue accession no.
LEA protein (Caenorhabditis elegans)
*LEA-like hypothetical protein (C. elegans)
C-type lectin (clec-49) (C. elegans)
*C-type lectin (clec-202) (C. elegans)
polygalacturonase (Aspergillus clavatus)
polygalacturonase (Penicillium griseorseum)
AHA-1, activator of Hsp90 (Bombyx mori) myo-inositol-1-phosphate synthase (C. elegans)
GABA/glycine receptor (ggr-3) (C. elegans)
hypothetical protein Y53G8B.4 (C. elegans) C2H2-type zinc finger protein (Brugia malayi)
*MSFI-like/14-3-3 protein (C. elegans)
vesicle transport (unc-16) (C. elegans)
DJ-1 like (Tribolium castaneum)
*saposin (spp-3) (C. elegans)
Sm-like protein (lsm-5) (C. elegans)
dual specificity phosphatase F28C6.8 (C. elegans)
Appr-1-p processing enzyme (B. malayi)
hypothetical protein F42A10.7 (C. elegans)
aspartate protease inhibitor (Onchocercavolvulus)
GATA-type zinc finger protein (B. malayi)
tropomyosin (Heterodera glycines)
transthyretin-related (Radopholus similis)
collagen family (col-14) (C. elegans)
glutathione peroxidase (Globodera rostochiensis)
high mobility group protein (hmg1.1) (C. elegans)
zinc/RING finger protein (mnat-1) (C. elegans)
Predicted ion channel (ced-11) (C. elegans)
rhoGAP domain (B. malayi)
hypothetical protein ZK228.3 (C. elegans)
A second example from the EST dataset, EF026240, also encodes a hydrophilic protein but it is less clearly identifiable as an LEA protein. The cDNA is again full length, with spliced leader (SL1a) and polyA tail, whose 464 bp insert encodes a protein of 102 amino acids in its longest open reading frame, with Mr 11656 and pI 6.53. The protein has a GRAVY score of -1.38, and a hydropathy plot shows it to be hydrophilic along almost its entire length. Disorder prediction algorithms are consistent with it being unfolded - VL-XT, Foldindex and IUPred predict 78%, 100% and 100% disorder respectively - and the Uversky plot positions the protein within natively unfolded space (Figure 1). POPP analysis shows clustering with Group 2 and Group 3 LEA proteins, but only at the level of band 3, indicating that this association is rather weak. Furthermore, there is no K11, and limited A11, periodicity and the predicted protein sequence includes one cysteine and one tryptophan, which it is unusual to find in Group 3 LEA proteins. A BLASTX run against the C. elegans genome returns lea-1 as the top hit, although the score is weak (0.002); a BLAST enquiry against the whole sequence collection returns as best hit (0.029) an LEA-like sequence from the Antarctic bacterium Polaromonas sp. JS666. Intriguingly, though, using BLAST to compare the EF026240 protein against AavLEA2, described above, gives a score of 5E-14, indicating a match of respectable significance. When aligned, the two predicted protein sequences are given as 34.3% identical with the ALIGN algorithm. In summary, the bioinformatics are consistent with this example being a hydrophilic, intrinsically disordered protein which probably has a distant evolutionary relationship to the Group 3 LEA proteins. For the moment, therefore, we consider it to be an LEA-like hypothetical protein.
While the presence of LEA and LEA-like protein gene sequences among the dehydration-induced EST dataset is in line with a presumed role for these polypeptides in the response to water stress, several other genes associated with stress responses were also found. For example, genes encoding: myo-inositol-1-phosphate synthase which controls synthesis of an osmolyte, myo-inositol, implicated in desiccation tolerance [37, 38]; AHA-1, a stress-regulated activator of the ATPase activity of molecular chaperone Hsp90 ; and glutathione peroxidase, which catalyses the reduction of hydroxyperoxides by glutathione and is important in the response to oxidative stress (reviewed in ).
The control of the response to desiccation is very poorly understood and therefore it is of interest that a number of genes involved in signal transduction and the regulation of gene activity were recovered from the screen. For example, three zinc finger protein genes were identified (Table 1), encoding a C2H2 zinc finger family member (GR463901), a GATA-type zinc finger (GR463910), and a RING finger domain protein (DR121009). The latter is most likely involved in the regulation of protein activity, rather than gene regulation directly; it gives a best BLAST hit with mnat-1, a C. elegans gene whose protein is probably involved in the ubiquitination pathway as an E3-type ubiquitin ligase. It is also implicated in stabilising the CDK-activating kinase intrinsic to cell cycle progression through its MAT1 domain and is functionally associated with TFII, the transcription/DNA repair factor (http://www.wormbase.org; ). The zinc finger proteins represented by GR463901 and GR463910 would merit further investigation as candidate anhydrobiotic gene set regulators: several transcription factors were found to be upregulated by water loss in Steinernema feltiae, S. carpocapsae and Plectus murrayi, and a GATA-type transcription factor (ELT-2) has been found to regulate C. elegans transcription during the innate immune response and osmotic stress ; in plants, many gene regulators have been characterised in the response to drought (e.g. reviewed in ). Another gene, represented by GR463906, encodes a dual specificity phosphatase able to hydrolyse phosphate from Tyr and Ser/Thr residues; such phosphatases are key players in signal transduction, for example, in MAPK pathways, and in cell cycle regulation . DR121003 encodes a homologue of C. elegans UNC-16 which is primarily involved in vesicle transport, but whose protein is also known to interact with the MAPK, JNK, and also with JNK kinases; it thus functionally intersects with signal transduction pathways. Further investigation of these genes could therefore lead to improved understanding of the dehydration stress response in nematodes.
Some predicted physico-chemical parameters of the novel ESTs identified by enrichment for dehydration-induced gene expression in A. avenae and their relative expression following exposure of the nematodes to four environmental stresses.
No. amino acids
Predicted GRAVY score
Evidence for horizontal gene transfer in A. avenae
Homologues of two polygalacturonase genes in the fungi Aspergillus awamori and Penicillium griseoroseum[46, 47] were found in the A. avenae EST set. These genes seem to be derived from the A. avenae genome and are not a result of fungal contamination since preliminary Southern hybridisation experiments show they are not present in Rhizoctonia solani, the food source used (data not shown). Polygalacturonases (EC:126.96.36.199) hydrolyse 1,4-alpha-D-galactosiduronic linkages in pectin, a major component of the cell walls of plants . Plant pathogenic bacteria and fungi secrete hydrolytic enzymes such as polygalacturonases and cellulases to degrade plant cell walls. Plant parasitic nematodes also secrete a great diversity of plant cell-wall degrading enzymes including cellulases and galacturonases. Phylogenetic analyses support the interpretation that the cellulases and polygalacturonases present in cyst and root-knot nematodes have been acquired by horizontal gene transfer from bacteria [49–53]. In contrast, the polygalacturonase genes of A. avenae appear to be derived by horizontal transfer from fungi (Additional File 1; Table 1). Interestingly a polygalacturonase gene found in the genome of the rice weevil Sitophilus oryzae also appears to have been incorporated into the genome of this insect by horizontal gene transfer from a fungal source . While the cellulase and polygalacturonase genes of the root-knot nematode M. incognita are most closely related to bacterial sequences, genome-wide analysis indicates that three other plant-cell wall degrading/modifying enzymes in this nematode (pectate lyase, arabinase and expansin) may be of fungal origin , and it has also been shown that the cellulases of the pine wood nematode are most similar to cellulases from fungi . A. avenae is routinely cultured in the laboratory on fungi such as Rhizoctonia solani or Botrytus cinerea, but A. avenae is also capable of growing and reproducing in vitro on tissue cultures prepared from several plant species  and it has also been shown that A. avenae will feed and multiply on plant roots, and suppress growth . The occurrence of a plant cell wall degrading enzyme in A. avenae provides further evidence that in addition to feeding on fungi A. avenae might also utilise plant tissues as a food source in nature.
Multiple variant spliced leaders in A. avenae mRNAs
Twelve different spliced leader sequences in A. avenae mRNAs.
EST accession no.
EF026240, GR463910, GR463914, GR463901
SL1 [C. elegans]
SL2 [C. elegans]
Regulation of gene expression after exposure to environmental stress
The expression of 42 A. avenae genes in response to evaporative water loss (90% RH, 24 h), increased osmotic potential (1 M sucrose, 24 h), cold (4°C, 24 h) and heat (32°C, 24 h) was evaluated by quantitative PCR (Table 1; Table 2). Two thirds of the recognisable genes were found to be upregulated by drying, thus validating the selection procedure and implicating the EST set in anhydrobiosis. Interestingly, the large majority of genes analysed were also responsive to osmotic stress and cold (83% and 77%, respectively), probably reflecting similar effects on water activity and overlapping stress responses to all three environmental conditions. In contrast, heat stress caused relatively few genes to be induced, with 30% responding to elevated temperature, implicating a rather different response mechanism.
The degree of expression varied widely between genes although ANOVA analysis showed that, in most cases, a 1.5-fold increase in expression levels was significant (Table 1). Nevertheless, some genes were highly upregulated: the LEA protein and LEA-like genes (EF026241 and EF026240) showed a 43-fold and 36-fold increase in expression level, respectively, in response to desiccation, and a response to osmotic stress in excess of 20-fold. These genes were also substantially upregulated on exposure to low temperature (8-fold and 47-fold respectively). Other genes were not induced to the same levels in most cases, but notable results include that of a glutathione peroxidase (GR463914) whose mRNA concentration increased 32- and 14-fold in response to desiccation and osmotic stress, respectively; a DJ-1-like gene (GR463903) showed a 10-fold induction on desiccation; and an AHA1 activator gene (GR463897) also showed high relative expression levels following exposure to desiccation and osmotic stress. The above results are consistent with a requirement for increased antioxidant and molecular shield/chaperone activity during water stress, as has been noted previously (e.g. ). Of the novel genes (Table 2), four showed more than 10-fold upregulation on osmotic stress and two displayed over five-fold induction in response to desiccation.
A. avenae is recalcitrant to RNA interference
Although the selective cloning procedure and the expression data of Tables 1 and 2 implicate the large majority of the gene set in anhydrobiosis, it would be desirable to obtain more direct evidence that they are required for survival of the dry state. To this end, we sought to establish RNA interference (RNAi) in A. avenae as a tool to selectively eliminate or reduce gene activity. Populations of A. avenae in which expression of particular genes had been reduced by gene silencing could then be assessed for desiccation tolerance; where survival was compromised, the respective genes would be strongly implicated in anhydrobiosis. RNAi occurs in C. elegans after injection of double-stranded RNA (dsRNA) corresponding to the target gene into worms , after soaking nematodes in a dsRNA solution  or after feeding nematodes on bacteria expressing dsRNA . The last of these methods is not available to us since A. avenae do not feed on bacteria. Therefore, attempts were made at injecting dsRNA into A. avenae but these were unsuccessful due to technical difficulties relating to the small size of the worms. The injection method is also laborious and is not amenable to studies at the population level, which is usually required for assessment of desiccation tolerance. A soaking method was therefore attempted using dsRNA corresponding to the RNA polymerase II large (α) subunit gene, ama-1, which is known to give a lethal phenotype when silenced in C. elegans and other nematodes such as the human parasitic species Brugia malayi. Two cDNA sequences corresponding to the A. avenae ama-1 transcript were cloned using primers directed against conserved peptide sequences identified in an alignment of nematode and other RNA polymerase II large subunit genes . Soaking of A. avenae for 24 h in 1 mg/ml ama-1 dsRNA did not, however, show a phenotype compared to negative control experiments (data not shown). Positive technical controls were performed using C. elegans and cognate ama-1 sequences, resulting in sick animals within 24 h. Further controls were also carried out using unc-22 and GFP dsRNA on wild-type and GFP-expressing C. elegans (PD4251) to confirm the effectiveness of RNAi.
Cross-species RNAi of anhydrobiosis genes
One strategy which might partially circumvent difficulties with RNAi in A. avenae is cross-species gene silencing using A. avenae sequences to reduce expression of homologous transcripts in another anhydrobiotic nematode. Cross-species RNAi has been shown using sequences from the animal parasitic nematode, Ascaris suum, in C. elegans, for example . Recent work by one of our groups (AMB) has demonstrated efficient systemic RNAi in Panagrolaimus superbus, a nematode in the same phylogenetic group as A. avenae with a marked resistance to desiccation. Indeed, since P. superbus is a bacterivore, it is susceptible to large scale RNAi in nematode populations using a feeding method very similar to that used routinely in C. elegans.
From the A. avenae EST set, 20 candidates (indicated in bold in Tables 1 and 2) were selected for RNAi in P superbus. The majority of these were chosen because they showed a high degree of similarity to sequences in C. elegans and/or Caenorhabditis briggsae, the rationale being that they were then more likely to have similar counterparts in P. superbus: both P. superbus and A. avenae are tylenchid nematodes, and thus more closely related in evolutionary terms than either is to the genus Caenorhabditis. For comparison, three of the 20 ESTs chosen were novel sequences (GR463919, GR463921, GR463927).
P. superbus were fed bacteria expressing dsRNA corresponding to individual A. avenae sequences, after which worms were first pre-conditioned at 90% RH for 24 h followed by 24 h desiccation at 10% RH. The worms were then rehydrated overnight in buffer before scoring for survival. Worms fed dsRNA corresponding to GFP and treated similarly were taken as negative controls. Prior to desiccation, nematodes were also subjected to phenotypic analysis, but none of the experiments performed revealed any extraordinary change in morphology or behaviour, either of adults, larvae or embryos (Additional File 2).
Glutathione peroxidase gene activity contributes to anhydrobiosis in P. superbus
If the cross-species strategy used above is valid, then genes that are upregulated in response to desiccation in A. avenae, and that have well-conserved sequence counterparts in P. superbus, should be part of an anhydrobiotic gene set in both nematodes. Therefore, we attempted to identify genes in P. superbus which were homologues or paralogues of the two genes indicated by the cross-species RNAi experiments, in order to carry out similar experiments with cognate P. superbus sequences. Consequently, an EST library of ~9,000 P. superbus cDNAs (Tyson et al., in preparation) was searched for sequences corresponding to A. avenaeGR463914 and GR463921. Although no matches were found for the novel sequence represented by GR463921, two different sequences showing similarity to GR463914 were obtained, Ps92 (accession number GR881192) and Ps114 (accession number GR881191), both of which are similar to glutathione peroxidase sequences in GenBank and other databases. For example, both Ps92 and Ps114 gave highly significant matches with the C. elegans glutathione peroxidase gene C11E4.1, showing BLASTX scores of 2E-97 and 1E-86, respectively. The ESTs probably represent different genes in P. superbus since they are only 66% identical using the ALIGN tool.
Anhydrobiosis is one of the most fascinating biological phenomena, yet the molecular mechanisms involved are still relatively poorly understood. The application of recently-developed molecular techniques to the study of anhydrobiosis has been hampered by the lack of a suitable desiccation-tolerant model organism, certainly among the metazoans, for which genome sequence and genetic tools are not yet available. Perhaps the best-characterised metazoan example is the nematode, A. avenae, in which the physiology and biochemistry of anhydrobiosis has been studied since the 1970s (e.g. ). To gain insight into the set of genes employed by A. avenae to achieve the anhydrobiotic state, which can be termed the anhydrobiotic gene set, or the desiccome , we have used a selective cloning technique to identify genes that are upregulated by evaporative water loss. The majority of genes identified were confirmed as desiccation-inducible using real-time qPCR, thus validating the cloning procedure; many examples of this gene set were also found to be responsive to osmotic stress and reduced temperature, possibly due to a common effect on water activity in all three stress conditions. Relatively few members of the EST set were responsive to elevated temperature, however. This fits with general observations in other invertebrates and plants where a large overlap in the stress responses to dehydration, osmotic stress and cold is observed, but where heat stress activates a different response pathway . The complexity of the response to water loss is illustrated by the involvement of molecular chaperone/shield proteins, redox balancing systems, osmolyte biosynthetic enzymes, C-type lectins, protein and polysaccharide metabolic enzymes, as well as genes controlling gene transcription and signal transduction pathways. In particular, new LEA and LEA-like protein genes were identified suggesting that several proteins of this type are involved in the response to water stress in invertebrates, as in plants. It seems unlikely, however, that the diversity of the LEA proteome in invertebrates will rival that of plants, where >50 LEA protein genes can be found in a single organism [73, 74]; in C. elegans, for example, only three such genes have been noted .
While such information on transcriptional stress responses can be extremely useful in assessing which genes and regulatory systems are involved in combating the stress imposed, it is becoming clear, at least in some organisms such as yeast, that not all stress-induced genes are essential for stress tolerance. Thus, in Saccharomyces cerevisiae, for example, very few genes required for DNA repair were induced by DNA-damaging agents . Other researchers report similar conclusions for yeast subjected to osmotic stress , oxidative stress , and anaerobic growth conditions . It is therefore important, if possible, to go beyond a purely correlative study, where stress-induced genes are catalogued, to experiments where genes of interest can be mutated or silenced and the effect of this on stress tolerance determined.
Silencing of specific gene targets by RNAi has been invaluable in nematodes, in particular C. elegans, as an alternative to mutagenesis that allows gene function to be probed on a large scale [79, 80]. Such high throughput functional genomics would be ideal for identifying the anhydrobiotic gene set in a nematode such as A. avenae. Unfortunately, attempts to perform RNAi on A. avenae by microinjection, or by soaking in dsRNA corresponding to the gene target, were not successful. This might be due to resistance or partial resistance to RNAi in this species, as has been found in some animal parasitic nematodes . However, with the soaking method, recalcitrance might simply reflect poor ingestion of surrounding medium by A. avenae, since fluorescently labelled dsRNA was not visible in its intestine. Alternative RNAi strategies might be to attempt a biolistic delivery method, as recently demonstrated in dehydration-stressed barley leaves , or to introduce dsRNA by electroporation, as described for the gastrointestinal parasite Trichostrongylus colubriformis. However, such RNAi techniques would be difficult to adapt for a high throughput screening experiment.
The cross-species method developed here, which exploits the relative ease with which RNAi can be performed in P. superbus, was extremely useful in allowing a limited screen of A. avenae sequences. Only two out of 20 A. avenae ESTs tested gave a marked reduction in P. superbus desiccation tolerance, but this probably reflects a requirement by the RNAi machinery for precise sequence matching between the short interfering RNAs generated from dsRNA and the target mRNA ; sequence divergence of homologous RNA targets will reduce the efficacy of RNAi across species. It is also formally possible that, as in yeast, many genes induced by water stress in A. avenae are not required for stress tolerance in P. superbus, or indeed in A. avenae. Nevertheless, the observation that silencing of two particular transcripts reduces survival of drying in P. superbus is the first such demonstration in any desiccation-tolerant organism and demonstrates the potential for RNAi in dissecting anhydrobiosis and determining the desiccome. This result also demonstrates the advantage of P. superbus as a model anhydrobiote. Accordingly, a large cognate EST dataset has already been assembled for high throughput screening by RNAi for genes required in anhydrobiosis in this nematode; a number of interesting candidates have been identified and are currently being characterised (Tyson et al., unpublished). The availability of a panel of other Panagrolaimus species of varying degrees of desiccation tolerance will also be useful for comparative studies of the phenomenon .
Of the two A. avenae ESTs that gave a phenotype in the cross-species RNAi desiccation experiments, one was a novel sequence not represented in sequence collections. This EST, GR463921, is only 350 bp in length and is not full length, since it lacks clearly defined 5' and 3' ends, as would be signified by a spliced leader and polyA tail, respectively (although a proportion of nematode mRNAs lack spliced leaders; ). Nevertheless, it shows modest upregulation during desiccation and osmotic stress, and merits further investigation by isolation of full length cDNA clones from A. avenae and a search for cognate sequences in P. superbus; the latter must exist for the cross-species RNAi to have been successful. However, an alternative role for GR463921 might be as a non-coding RNA involved in regulating anhydrobiosis in nematodes. A precedent for this exists in the desiccation-tolerant resurrection plant, C. plantagineum, where an abundant, small non-coding RNA (CDT-1) is able to induce tolerance in callus tissue . CDT-1 is a member of a family of retroelements which govern synthesis of siRNA involved in control of desiccation tolerance pathways . Formation of dsRNA is probably important for CDT-1 function and, intriguingly, the GR463921 sequence is predicted to form a high degree of secondary structure (Additional File 3) by the RNAfold algorithm (; http://rna.tbi.univie.ac.at). It will therefore be of interest to investigate this transcript further with a view to understanding its precise role in nematode anhydrobiosis.
The second EST giving an RNAi phenotype on desiccation, which matches glutathione peroxidase sequences in GenBank and other sequence collections, was a full length cDNA of 649 bp including its spliced leader (SL1a) but excluding its polyA tail. It appears to contain an open reading frame of 151 amino acids (Mr 17480) bounded by a UGA stop codon. However, some glutathione peroxidases are known to contain selenocysteine [88, 89] which is not directly represented in the genetic code, but is incorporated into nascent protein when the ribosome reads UGA in the presence of a selenocysteine insertion element in the mRNA . Therefore, the predicted open reading frame might be shorter than the actual protein length; indeed, there are two in-frame UGA codons before a third, TAA, stop codon is encountered, consistent with up to two selenocysteines being incorporated into the A. avenae glutathione peroxidase. The two P. superbus glutathione peroxidase ESTs, GR881191 (Ps114) and GR881192 (Ps92), are somewhat longer than the A. avenae example, at 735 bp and 808 bp, respectively. Ps114 does not quite extend to the 3' end, since its open reading frame of 230 amino acids lacks a stop codon, while the reading frame of Ps92 is likely complete at 228 residues. Interestingly, both predicted proteins from P. superbus have a hydrophobic sequence of about 20 amino acids at their N-terminal ends consistent with an endoplasmic reticulum translocation signal. Both proteins are therefore probably secreted; indeed, the best BLAST matches are with cuticular glutathione peroxidases from other nematodes. For example, the Ps114 protein gives a score of 7E-70 with the cuticular glutathione peroxidase from Brugia pahangi, also known as its major surface antigen gp29. In contrast, the A. avenae protein does not contain a signal peptide and therefore is probably intracellular. As such, the A. avenae and P. superbus glutathione peroxidases described are likely to have different locations within each nematode, and to have different substrate preferences, with the secreted P. superbus examples not metabolising hydrogen peroxide but instead larger hydroperoxide substrates such as lipid peroxides. The fact that both types of enzyme are implicated in anhydrobiosis is consistent with the need to protect both intracellular and membrane targets from oxidative stress during desiccation.
To further investigate the role of antioxidants in nematode anhydrobiosis, we also identified a third P. superbus EST, Ps41 (accession number GR881190), which is likely to encode a peroxiredoxin (also called thioredoxin peroxidase) since it gave a highly significant BLASTX match (6E-47) with the C. elegans prdx-6 gene. Like glutathione peroxidases, peroxiredoxins are also thiol peroxidases involved in the control of cellular peroxide levels, and both classes of enzyme are evolutionarily related . Preliminary, unpublished results suggest that RNAi treatment with peroxiredoxin dsRNA also reduces survival of desiccation in P. superbus, consistent with a role for a second type of thiol peroxidase in nematode anhydrobiosis.
Oxidative stress is experienced by organisms undergoing a wide range of abiotic stress conditions, probably due to malfunction of mitochondria (and chloroplasts in plants) resulting in generation of reactive oxygen species . Accordingly, an oxidative stress response forms part of the environmental stress response defined by Gasch et al.  in yeast and is part of the minimal stress response described by Kültz . Desiccation stress is no exception and it has long been recognised that antioxidants are likely to play a role in desiccation tolerance (reviewed in [95–97]). Recent global studies in plants have reinforced this perception (reviewed in ) and, for example, peroxiredoxin and glutathione peroxidase, amongst other antioxidants, are induced by dehydration in desiccation-tolerant resurrection plants [99, 100].
The demonstration of a role for glutathione peroxidase in nematode anhydrobiosis is consistent with, and validates, the many correlative studies which have previously implicated antioxidants. In nematodes, for instance, three studies have identified dehydration-induced genes encoding antioxidant enzymes: Gal et al.  described a number of ESTs from S. feltiae, a desiccation-resistant entomopathogenic species, including a glutathione peroxidase; in A. avenae, a glutaredoxin gene was induced by drying, together with an LEA protein gene and a gene encoding another hydrophilic protein, anhydrin ; and, in a recent large-scale study of the Antarctic nematode P. murrayi by Adhikari et al. (; and see review by Wharton and Marshall ), genes encoding glutathione-S-transferase, superoxide dismutase, peroxiredoxin and glutathione peroxidase were included in a dehydration-induced EST panel. The present study takes our understanding of anhydrobiosis a step further by showing a reduced level of desiccation tolerance after silencing glutathione peroxidase transcripts, thus clearly demonstrating a stress-combative role for antioxidants. Further work will explore gene silencing more extensively in P. superbus to determine the desiccome in this species in more detail.
Using a combination of differential cloning, expression profiling and cross-species RNAi, we have identified a panel of genes implicated in the response to desiccation in the anhydrobiotic nematode, A. avenae. The gene set includes novel hydrophilic proteins, one of which is a previously-undescribed LEA protein, and other stress response genes controlling synthesis of compatible solutes, chaperone systems and antioxidants. The large majority of these genes were upregulated by evaporative water loss, and many also by osmotic upshift and low temperature, conditions which also involve water stress. We also presented evidence for RNAi recalcitrance in A. avenae, but were able to show cross-species RNAi of two A. avenae sequences in another anhydrobiotic nematode, P. superbus. One of these was a novel sequence, which is potentially a non-coding regulatory RNA, and the second encoded an intracellular glutathione peroxidase: both reduced survival of desiccation of P. superbus after RNAi treatment. Two cognate glutathione peroxidase sequences from P. superbus were also identified and, when used for RNAi, were shown to have a similar effect on desiccation tolerance. These results therefore demonstrate the complex nature of desiccation tolerance and emphasise the need for effective antioxidant systems in anhydrobiosis.
Culture, desiccation and stress treatment of nematodes
Aphelenchus avenae was grown in the dark at 20°C on the fungus Rhizoctonia solani, itself grown on a substrate of wheat . After a growth period of 2-3 weeks (before swarming of the nematodes in the jars), nematodes were rinsed with tap water and brought to a concentration of approximately 3,300 worms per ml. Nematodes were cleaned in 30% (w/v) sucrose solution, followed by three washes in tap water and the nematodes were transferred to 25 mm Supor® Membrane Disc Filters (0.45 μm, Pall Corporation, Michigan, USA; 30,000 nematodes per filter) by vacuum filtration. These were placed in a desiccation chamber at 20°C for 24 h over a saturated solution of MgSO4.7H2O (90% RH), or over a saturated solution of KCl (85% RH) .
For the qPCR experiments the nematodes were grown and harvested as described above but they were cleaned by placing the nematode pellet in a sieve containing a layer of Kleenex® allowing them to crawl through the paper (for 2-3 h) into the water filled dish in which the sieve had been placed. A. avenae populations were exposed to individual stresses as follows. Heat stress: one replicate consisted of 1 ml of the above worm suspension placed in a 30 mm Petri dish incubated at 32°C for 24 h. Cold stress: one replicate consisted of 1 ml of worm suspension in a 30 mm Petri dish incubated at 4°C for 24 h. Desiccation: one replicate consisted of 1 ml of worm suspension vacuum filtered onto a 25 mm Supor® membrane disc and placed in a desiccation chamber set at 90% RH (using a saturated solution of MgSO4.7H2O) for 24 h. Osmotic stress: one replicate consisted of 1 ml of worm suspension added to 1 ml of 1 M sucrose solution in a 30 mm Petri dish and placed at 20°C for 24 h. Control nematodes consisted of 1 ml of worm suspension in 30 mm Petri dishes incubated at 20°C for 24 h. Three biological replicates were set up for each stress and for control worms. These stress conditions yielded an 80% survival rate for each stress compared to a 90% survival rate for the control worms.
P. superbus and C. elegans were grown on NGM plates and fed Escherichia coli strain OP50, or strain HT115 carrying an empty L4440 plasmid, in the dark at 25°C.
Construction of a desiccation induced EST library
Total RNA was extracted from desiccated and control nematodes using TRI® Reagent (Sigma-Aldrich GmbH, Steinheim, Germany). Total RNA from desiccated nematodes was converted to 5'-oligo-capped cDNA using the GeneRacer™ Kit (Invitrogen, Carlsbad, California, USA), following the manufacturer's instructions. The GeneRacer™ Oligo dT primer was used in the reverse transcription reaction. PCR was carried out on the 5'-oligo-capped cDNA with the GeneRacer™ 5'-primer (CGACTGGAGCACGAGGACACTGA) and the GeneRacer™ 3'-primer (GCTGTCAACGATACGCTACGTAACG) using Platinum®Taq DNA Polymerase (Invitrogen). The product of this reaction was cloned into the pCR2.1 TOPO vector using a TOPO-TA Cloning kit (Invitrogen).
Reverse Northern hybridization
Recombinant clones were PCR amplified in a 96 well plate using standard M13 forward and reverse primers and the presence of an insert was confirmed by electrophoresis on a 1% agarose gel. The PCR products were denatured by adding NaOH to a total concentration of 0.4 M and EDTA to a total concentration of 0.04 M and heating at 99°C for 10 min. The denatured DNA was blotted onto Hybond™-N+ transfer membrane (Amersham Biosciences, Buckinghamshire, UK) in duplicate using a SRC 60-D Dot Blot Minifold (Schleicher & Schuell, Dassel, Germany). The membranes were air dried at 55°C for 30 min and UV-crosslinked using a UV Stratalinker® 1800 (1200 mJ; Stratagene, California, USA).
Screening the EST library
Total RNA was extracted from desiccated (90% RH for 24 h) and control nematodes as before . The RNA was resuspended in 30 μl DEPC-treated H2O, concentration determined by measuring absorbance at 260 nm and adjusted to 1 μg/μl. 5 μg of RNA was converted to cDNA using the SuperScript™ First-Strand Synthesis System (Invitrogen). Experimental (desiccated) and control probes were constructed by radiolabelling this cDNA with [32P]-dCTP using the Prime-a-Gene® Labelling System (Promega, Southampton, UK) in accordance with the manufacturer's instructions. One replicate of each dot blot membrane was hybridized with the experimental probe and the corresponding replicate was hybridized with the control probe. Hybridisation was carried out in hybridisation buffer (5% dextran sulfate, 1 M NaCl, 1% SDS in water) at 60°C for 18 h in glass tubes with gentle rotation. The membranes were washed twice in wash buffer (10× SSC, 0.5% SDS in water) for 4 min at room temperature followed by two washes for 15 min at 60°C. The membranes were then exposed on X-OMAT AR photography paper (Kodak, New York, USA) for 24 h. The control and experimental dot blots were digitally scanned and analyzed using the ImageQuant TL Array Analysis software (Amersham Biosciences). Clones considered to be upregulated in response to desiccation were selected for sequencing and sent to Agowa GmbH (Berlin, Germany). The resulting sequences were assembled into contigs using the online sequence assembly program CAP3 http://pbil.univ-lyon1.fr/cap3.php, yielding a dataset of 88 unique ESTs, enriched in those upregulated in response to desiccation.
The POPPs  is a suite of applications, one of which clusters protein sequences based on similarities in their peptide compositions, in this case to see whether the input sequence clusters with any of the LEA protein groups and, if so, which. The GRAVY calculation was undertaken by a small Python application which uses the Eisenberg Consensus scale  in preference to the Kyte-Doolittle scale , which has well known anomalies. Three predictors of unfolded proteins were used in this study: Foldindex (; http://bip.weizmann.ac.il/fldbin/findex/, VL-XT (; http://www.pondr.com) and IUPred (; http://iupred.enzim.hu). Uversky plots  of charge vs. hydropathy were performed at pondr.com and raw data downloaded to construct graphical output in GraphPad Prism version 4.0b (GraphPad Software, California, USA). An application called Perdix was used to search the input sequences for amino acid periodicity, i.e. amino acids that recur every N amino acids. So, for example, in the peptide GAGPG the G recurs with period 2. The application Perdix, and the earlier application for computing hydrophobicity and other peptide statistics, both currently unpublished, can be obtained from Wise. The pairwise sequence alignment application ALIGN  is part of the San Diego Supercomputer Workbench http://workbench.sdsc.edu.
Quantitative real time PCR
The replicate worm samples from each individual stress treatment were pooled and RNA was extracted using Trizol reagent for each stressed and control worms followed by treatment with DNAse I (RNase free, Invitrogen). The concentration and quality of RNA was determined using a Thermo Scientific Nanodrop Spectrophotometer. Total RNA (1 μg per reaction) was converted to cDNA using the Roche Transcriptor First Strand cDNA Synthesis Kit. Primers were designed for each gene to amplify a fragment of approximately 125 bp (see Additional File 4). Real time qPCR reactions were carried out with a Roche LightCycler 480 thermocycler using Roche SYBR I Master Mix with 0.0025 pmol primers (each), 500 ng of sample cDNA template and 0.5 U heat-labile uracil-N-glycosolase per reaction. Relative expressions were calculated using the second derivative maximum method with A. avenae 5.8S rRNA and ama-1 genes as reference genes. Statistically significant differences in expression levels were confirmed using ANOVA with a confidence level of 99.9% and Dunnett's comparison test with α = 1.
Assessing A. avenae for RNAi capability
Two ama-1 cDNA fragments were obtained from A. avenae by RT-PCR. The short fragment (322 bp), was obtained with the primers AaAMA1s (5'-TTTTTCATGCAATGGGCGGTC-3'; forward) and AaAMA1as (5'-ACATAAACCTCTCTCACGACG-3'; reverse). Reverse transcription was performed using SuperScript II (Invitrogen) according to manufacturer's instructions, and the PCR was done under the following conditions: 95°C 5 min; then 35 cycles of 95 C 30 s, 57°C 30 s and 72°C 1 min; with a final extension at 72°C for 10 min. A longer version of ama-1 cDNA (800 bp) was cloned using the same RT-PCR conditions with primers Aa.AMA1.s2 (5'-ACGTCGACTGATTAAGGCCAT-3'; forward) and Aa.AMA1.as3 (5'-TGATGATCTCCTTCAGACGAG-3'; reverse). PCR products underwent a second amplification round using T7-tagged primers. These new amplicons were then used for in vitro transcription using Ampliscribe T7 Flash Transcription Kit (Epicentre) according to manufacturer's instructions. For RNAi attempts, worms were soaked in 20 μl soaking buffer (5 mM Tris-HCl pH7.0) containing dsRNA (1 μg/μl) for 24 h and then plated on LB agar plates (50 μg/ml carbenicillin). Worms were placed over a small spot of fungus (Botrytis cinerea) previously grown on the LB agar plate for 24 h. Survival analysis was performed by observation under a microscope. C. elegans was treated in a similar fashion except worms were incubated with GFP and unc-22 dsRNA for 24 h at ambient temperature and recovered on E. coli-seeded NGM plates until phenotypes were observed. For dsRNA uptake analysis, dsRNA molecules were first synthesised by in vitro transcription from undigested L4440-GFP plasmids. Then 5 μg dsRNA was labelled using Label IT CX-rhodamine Nucleic Acid Labeling Kit according to manufacturer's instructions (Mirus/Cambridge Bioscience). Worms were soaked in 20 μl soaking buffer for 22 h at 15°C in the dark. Uptake efficiency was assessed using a confocal microscope (Zeiss LSM510 META).
Cross-species RNA interference and desiccation tolerance
Plasmid DNA of selected A. avenae EST clones (in TOPO vector; Invitrogen) were digested using Sac I (NEB) and Xho I (NEB), inserts were gel-purified and ligated into Sac I/Xho I -digested L4440 feeding vector . Competent HT115 cells were transformed with ligation products and selected on LB ampicillin (100 μg/mL). Subcloning was confirmed by digesting L4440 plasmids with Sac I/Xho I. Transformed HT115 clones containing individual plasmids for RNAi were each grown in 10 ml LB ampicillin (50 μg/mL) overnight, recovered by centrifugation and resuspended in 300 μl plain LB medium. Bacteria were transferred to Petri dish plates (9 cm diameter) containing NGM agar plus IPTG (240 μg/ml) and ampicillin (50 μg/ml). Plates were left for 24 h at ambient temperature to allow dsRNA induction. P. superbus was prepared for RNAi as follows: a large (140 mm) NGM plate containing mixed stage animals with many gravid adults was harvested, worms washed 3× in M9 buffer and synchronised by bleaching (3 min in 1% NaOCl, 0.5 M NaOH). After bleaching, worms were washed 3× in M9 and left to hatch in M9 on a shaker for 48 h. Approximately 1500 worms were placed onto each feeding plate (three per test gene) where they remained for approximately 12 days in the dark at 25°C before desiccation. Worms were then collected, washed 3× in M9 buffer, counted and placed onto three 25 mm Supor® Membrane Disc Filters by vacuum filtration (approximately 300 RNAi-treated worms per filter) and processed as follows: 24 h at 90% RH (saturated BaCl2); 24 h at 10% RH (dry silica gel) and overnight in M9 buffer. Non-desiccated control worms were processed as follows: 48 h at 100% RH (UHP H2O) and 24 h in M9 buffer. All treatments were in 350 ml sealed plastic boxes. After rehydration, viability was assessed under the microscope, where 80 worms were randomly chosen for analysis from each filter (spontaneously moving worms were considered alive). C. elegans fed with a clone expressing unc-22 dsRNA was used as a positive control of RNAi treatment (twitching phenotype). Negative controls consisted of P. superbus and C. elegans worms fed with a clone expressing GFP (green fluorescent protein) dsRNA. Statistical relevance was determined by one-way ANOVA and a Tukey post hoc test using InStat3 (GraphPad Software, San Diego, CA).
We would like to thank Chiara Boschetti for her help with confocal microscopy. The support of the Leverhulme Trust (grants F/09 717/A and F/09 717/B) and Science Foundation Ireland (02/IN/B032) is gratefully acknowledged. MCB was funded by a BBSRC CASE studentship (BBS/S/L/2004/11164). BAC was funded by an Irish Research Council for Science Technology and Engineering EMBARK post graduate fellowship. TCP was in receipt of a fellowship from the CAPES Funding Agency of the Ministry of Education and Culture, Brazil.
- Crowe JH, Hoekstra FA, Crowe LM: Anhydrobiosis. Ann Rev Physiol. 1992, 54: 579-599. 10.1146/annurev.ph.54.030192.003051..Google Scholar
- Billi D, Potts M: Life and death of dried prokaryotes. Res Microbiol. 2002, 153: 7-12. 10.1016/S0923-2508(01)01279-7PubMedGoogle Scholar
- Clegg JS: Cryptobiosis - a peculiar state of biological organization. Comp Bioch Physiol B. 2001, 128: 613-624. 10.1016/S1096-4959(01)00300-1..Google Scholar
- Ingram J, Bartels D: The molecular basis of dehydration tolerance in plants. Ann Rev Plant Physiol. 1996, 47: 377-403. 10.1146/annurev.arplant.47.1.377..Google Scholar
- Bartels D, Schneider K, Terstappen G, Piatkowski D, Salamini F: Molecular cloning of abscisic acid modulated genes which are induced during desiccation of the resurrection plant Craterostigma plantagineum. Planta. 1990, 181: 27-34. 10.1007/BF00202321..PubMedGoogle 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..Google Scholar
- Shannon AJ, Browne JA, Boyd J, Fitzpatrick DA, Burnell AM: The anhydrobiotic potential and molecular phylogenetics of species and strains of Panagrolaimus (Nematoda, Panagrolaimidae). J Exp Biol. 2005, 208: 2433-2445. 10.1242/jeb.01629PubMedGoogle 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-3105PubMedGoogle Scholar
- Tyson T, Reardon W, Browne JA, Burnell AM: Gene induction by desiccation stress in the entomopathogenic nematode Steinernema carpocapsae reveals parallels with drought tolerance mechanisms in plants. Int J Parasitol. 2007, 37: 763-776. 10.1016/j.ijpara.2006.12.015PubMedGoogle Scholar
- Elbein AD, Pan YT, Pastuszak I, Carroll D: New insights on trehalose: a multifunctional molecule. Glycobiology. 2003, 13: 17R-27R. 10.1093/glycob/cwg047PubMedGoogle Scholar
- Carpenter JF, Crowe LM, Crowe JH: Stabilization of phosphofructokinase with sugars during freeze-drying: characterization of enhanced protection in the presence of divalent cations. Biochim Biophys Acta. 1987, 923: 109-115.PubMedGoogle Scholar
- Crowe JH, Carpenter JF, Crowe LM: The role of vitrification in anhydrobiosis. Ann Rev Physiol. 1998, 60: 73-103. 10.1146/annurev.physiol.60.1.73..Google Scholar
- Goyal K, Browne JA, Burnell AM, Tunnacliffe A: Dehydration-induced tps gene transcripts from an anhydrobiotic nematode contain novel spliced leaders and encode atypical GT-20 family proteins. Biochimie. 2005, 87: 565-574. 10.1016/j.biochi.2005.01.010PubMedGoogle Scholar
- Higa LM, Womersley CZ: New insights into the anhydrobiotic phenomenon: the effects of trehalose content and differential rates of evaporative water-loss on the survival of Aphelenchus avenae. J Exp Zool. 2003, 267: 120-129. 10.1002/jez.1402670205..Google Scholar
- Lapinski J, Tunnacliffe A: Anhydrobiosis without trehalose in bdelloid rotifers. FEBS Lett. 2003, 553: 387-390. 10.1016/S0014-5793(03)01062-7PubMedGoogle Scholar
- Caprioli M, Katholm AK, Melone G, Ramløv H, Ricci C, Santo N: Trehalose in desiccated rotifers: a comparison between a bdelloid and a monogonont species. Comp Biochem Physiol A. 2004, 139: 527-532. 10.1016/j.cbpb.2004.10.019..Google Scholar
- McGee BM: Hydrophilic proteins in the anhydrobiosis of bdelloid rotifers. PhD thesis. 2006, University of Cambridge, Institute of BiotechnologyGoogle Scholar
- Hengherr S, Heyer AG, Koehler HR, Schill RO: Trehalose and anhydrobiosis in tardigrades - evidence for divergence in responses to dehydration. FEBS J. 2008, 275: 281-288.PubMedGoogle Scholar
- Ratnakumar S, Tunnacliffe A: Intracellular trehalose is neither necessary nor sufficient for desiccation tolerance in yeast. FEMS Yeast Res. 2006, 6: 902-913. 10.1111/j.1567-1364.2006.00066.xPubMedGoogle Scholar
- Tunnacliffe A, Lapinski J: Resurrecting Van Leeuwenhoek's rotifers: a reappraisal of the role of disaccharides in anhydrobiosis. Phil Trans R Soc B. 2003, 358: 1755-1771. 10.1098/rstb.2002.1214PubMedPubMed CentralGoogle Scholar
- Tunnacliffe A, Wise MJ: The continuing conundrum of the LEA proteins. Naturwissenschaften. 2007, 94: 791-812. 10.1007/s00114-007-0254-yPubMedGoogle Scholar
- Browne JA, Tunnacliffe A, Burnell A: Plant desiccation gene found in a nematode. Nature. 2002, 416: 38- 10.1038/416038aPubMedGoogle Scholar
- Gal TZ, Glazer I, Koltai H: An LEA group 3 family member is involved in survival of C. elegans during exposure to stress. FEBS Lett. 2004, 577: 21-26. 10.1016/j.febslet.2004.09.049PubMedGoogle 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. Eucaryot Cell. 2004, 3: 966-75. 10.1128/EC.3.4.966-975.2004..Google Scholar
- Adhikari BN, Wall DH, Adams BJ: Desiccation survival in an Antarctic nematode: molecular analysis using expressed sequenced tags. BMC Genomics. 2009, 10: 69- 10.1186/1471-2164-10-69PubMedPubMed CentralGoogle Scholar
- Kikawada T, Nakahara Y, Kanamori Y, Iwata K-I, Watanabe M, McGee B, Tunnacliffe A, Okuda T: Dehydration-induced expression of LEA proteins in an anhydrobiotic chironomid. Biochem Biophys Res Comm. 2006, 348: 56-61. 10.1016/j.bbrc.2006.07.003PubMedGoogle Scholar
- Hand SC, Jones D, Menze MW, Witt TL: Life without water: expression of plant LEA genes by an anhydrobiotic arthropod. J Exp Zool. 2006, 305A: 1-5. 10.1002/jez.a.261..Google Scholar
- Pouchkina-Stantcheva NN, McGee BM, Boschetti C, Tolleter D, Chakrabortee S, Popova AV, Meersman F, Macherel D, Hincha DK, Tunnacliffe A: Functional divergence of former alleles in an ancient asexual invertebrate. Science. 2007, 318: 268-271. 10.1126/science.1144363PubMedGoogle Scholar
- Bahrndorff S, Tunnacliffe A, Wise MJ, McGee B, Holmstrup M, Loeschcke V: Bioinformatics and protein expression analyses implicate LEA proteins in the drought response of Collembola. J Insect Physiol. 2009, 55: 210-217. 10.1016/j.jinsphys.2008.11.010PubMedGoogle Scholar
- Tolleter D, Jaquinod M, Mangavel C, Passirani C, Payet N, Avelange-Macherel M-H, Macherel D: Drying reveals structure and function of a plant mitochondrial protein. Plant Cell. 2007, 19: 1580-1589. 10.1105/tpc.107.050104PubMedPubMed CentralGoogle Scholar
- Goyal K, Walton LJ, Tunnacliffe A: LEA proteins prevent protein aggregation due to water stress. J Biochem. 2005, 388: 151-157. 10.1042/BJ20041931..Google Scholar
- Chakrabortee S, Boschetti C, Walton LJ, Sarkar S, Rubinsztein DC, Tunnacliffe A: Hydrophilic protein associated with desiccation tolerance exhibits broad protein stabilization function. Proc Natl Acad Sci USA. 2007, 104: 18073-18078. 10.1073/pnas.0706964104PubMedPubMed CentralGoogle Scholar
- Wise MJ, Tunnacliffe A: POPP the question: what do LEA proteins do?. Trends Plant Sci. 2004, 9: 13-17. 10.1016/j.tplants.2003.10.012PubMedGoogle Scholar
- Tunnacliffe A, Hincha DK, Leprince O, Macherel D: LEA proteins: versatility of form and function. Dormancy and Resistance in Harsh Environments. Topics in Current Genetics. Edited by: Lubzens E, Cerda J, Clark MS. 2009, Berlin: SpringerGoogle Scholar
- Wise MJ: The POPPs: clustering and searching using peptide probability profiles. Bioinformatics. 2002, 18: S38-S45.PubMedGoogle Scholar
- Wise MJ: LEAping to conclusions: a computational reanalysis of late embryogenesis abundant proteins and their possible roles. BMC Bioinformat. 2003, 4: 52-10.1186/1471-2105-4-52. 10.1186/1471-2105-4-52Google Scholar
- Womersley C, Smith L: Anhydrobiosis in nematodes. 1. The role of glycerol myoinositol and trehalose during desiccation. Comp Biochem Physiol B. 1981, 70: 579-586. 10.1016/0305-0491(81)90299-6..Google Scholar
- Lehner A, Chopera DR, Peters SW, Keller F, Mundree SG, Thomson JA, Farrant JM: Protection mechanisms in the resurrection plant Xerophyta viscosa: cloning, expression, characterisation and role of XvINO1, a gene coding for a myo-inositol 1-phosphate synthase. Funct Plant Biol. 2008, 35: 26-39. 10.1071/FP07142..Google Scholar
- Meyer P, Prodromou C, Liao C, Hu B, Roe SM, Vaughan CK, Vlasic I, Panaretou B, Piper PW, Pearl LH: Structural basis for recruitment of the ATPase activator Aha1 to the Hsp90 chaperone machinery. EMBO J. 2004, 23: 1402-10. 10.1038/sj.emboj.7600141PubMedGoogle Scholar
- Dayer R, Fischert BB, Eggen RIL, Lemaire SD: The peroxiredoxin and glutathione peroxidase families in Chlamydomonas reinhardtii. Genetics. 2008, 179: 41-57. 10.1534/genetics.107.086041PubMedPubMed CentralGoogle Scholar
- Gervais V, Busso D, Wasielewski E, Poterszman A, Egly JM, Thierry JC, Kieffer B: Solution structure of the N-terminal domain of the human TFIIH MAT1 subunit - new insights into the RING finger family. J Biol Chem. 2001, 276: 7457-7464. 10.1074/jbc.M007963200PubMedGoogle Scholar
- Kerry S, TeKippe M, Gaddis NC, Aballay A: GATA transcription factor required for immunity to bacterial and fungal pathogens. PLoS One. 2006, 1: e77- 10.1371/journal.pone.0000077PubMedPubMed CentralGoogle Scholar
- Umezawa T, Fujita M, Fujita Y, Yamaguchi-Shinozaki K, Shinozaki K: Engineering drought tolerance in plants: discovering and tailoring genes to unlock the future. Curr Op Biotechnol. 2006, 17: 113-122.Google Scholar
- Patterson KI, Brummer T, O'Brien PM, Daly RJ: Dual-specificity phosphatases: critical regulators with diverse cellular targets. Biochem J. 2009, 418: 475-489.PubMedGoogle Scholar
- Nehrke K: A reduction in intestinal cell pHi due to loss of the Caenorhabditis elegans Na+ /H+ exchanger NHX-2 increases life span. J Biol Chem. 2003, 278: 44657-44666. 10.1074/jbc.M307351200PubMedGoogle Scholar
- Nagai M, Ozawa A, Katsuragi T, Kawasaki H, Sakai T: Cloning and heterologous expression of gene encoding a polygalacturonase from Aspergillus awamori. Biosci Biotechnol Biochem. 2000, 64: 1580-1587. 10.1271/bbb.64.1580PubMedGoogle Scholar
- Ribon AOB, Queiroz MV, Coelho JLC, de Araujo EF: Differential expression of polygalacturonase-encoding genes from Penicillium griseoroseum in different carbon sources. J Ind Microbiol Biotech. 2002, 29: 145-148. 10.1038/sj.jim.7000296..Google Scholar
- Ridley BL, O'Neill MA, Mohnen DA: Pectins: structure, biosynthesis, and oligogalacturonide-related signaling. Phytochemistry. 2001, 57: 929-967. 10.1016/S0031-9422(01)00113-3PubMedGoogle Scholar
- Yan YT, Smant G, Stokkermans J, Qin L, Helder J, Baum T, Schots A, Davis E: Genomic organization of four beta-1, 4-endoglucanase genes in plant-parasitic cyst nematodes and its evolutionary implications. Gene. 1998, 220: 61-70. 10.1016/S0378-1119(98)00413-2PubMedGoogle Scholar
- Jaubert S, Laffaire J-B, Abad P, Rosso M-N: A polygalacturonase of animal origin isolated from the root-knot nematode Meloidogyne incognita. FEBS Lett. 2002, 522: 109-112. 10.1016/S0014-5793(02)02906-XPubMedGoogle Scholar
- Jones JT, Furlanetto C, Kikuchi T: Horizontal gene transfer from bacteria and fungi as a driving force in the evolution of plant parasitism in nematodes. Nematology. 2005, 7: 641-646. 10.1163/156854105775142919..Google Scholar
- Abad P, Gouzy J, Aury JM, Castagnone-Sereno P, Danchin EGJ, Deleury E, Perfus-Barbeoch L, Anthouard V, Artiguenave F, Blok VC, Caillaud MC, Coutinho PM, Dasilva C, De Luca F, Deau F, Esquibet M, Flutre T, Goldstone JV, Hamamouch N, Hewezi T, Jaillon O, Jubin C, Leonetti P, Magliano M, Maier TR, Markov GV, McVeigh P, Pesole G, Poulain J, Robinson-Rechavi M, Sallet E, Segurens B, Steinbach D, Tytgat T, Ugarte E, van Ghelder C, Veronico P, Baum TJ, Blaxter M, Bleve-Zacheo T, Davis EL, Ewbank JJ, Favery B, Grenier E, Henrissat B, Jones JT, Laudet V, Maule AG, Quesneville H, Rosso MN, Schiex T, Smant G, Weissenbach J, Wincker P: Genome sequence of the metazoan plant-parasitic nematode Meloidogyne incognita. Nat Biotech. 2008, 26: 909-915. 10.1038/nbt.1482..Google Scholar
- Kyndt T, Haegeman A, Gheysen G: Evolution of GHF5 endoglucanase gene structure in plant-parasitic nematodes: no evidence for an early domain shuffling event. BMC Evol Biol. 2008, 8: 305- 10.1186/1471-2148-8-305PubMedPubMed CentralGoogle Scholar
- Shen Z, Denton M, Mutti N, Pappan K, Kanost MR, Reese JC, Reeck GR: Polygalacturonase from Sitophilus oryzae: Possible gene transfer of a pectinase gene from fungi to weevils. J Insect Sci. 2003, 3: 1-9. 10.1672/1536-2442(2003)003[0001:TMOAAE]2.0.CO;2Google Scholar
- Kikuchi T, Jones JT, Aikawa T, Kosaka H, Ogura N: A family of glycosyl hydrolase family 45 cellulases from the pine wood nematode Bursaphelenchus xylophilus. FEBS Lett. 2004, 572: 201-205. 10.1016/j.febslet.2004.07.039PubMedGoogle Scholar
- Barker KR, Darling HM: Reproduction of Aphelenchus avenae on plant tissiues in culture. Nematologica. 1965, 11: 162-166.Google Scholar
- Hodda M: Aphelenchida. Australian Faunal Directory. Australian Biological Resources Study, Canberra. 2003, Viewed 20 July 2009., http://www.environment.gov.au/biodiversity/abrs/online-resources/fauna/afd/taxa/APHELENCHIDAEGoogle Scholar
- Hastings KEM: SL trans-splicing: easy come or easy go?. Trends Genet. 2005, 21: 240-247. 10.1016/j.tig.2005.02.005PubMedGoogle Scholar
- Giuliano DB, Blaxter ML: Operon conservation and the evolution of trans-splicing in the phylum Nematoda. PLoS Genet. 2006, 2: e198- 10.1371/journal.pgen.0020198Google Scholar
- Pettitt J, Muller B, Stansfield I, Connolly B: Spliced leader trans-splicing in the nematode Trichinella spiralis uses highly polymorphic, noncanonical spliced leaders. RNA. 2008, 14: 760-770. 10.1261/rna.948008PubMedPubMed CentralGoogle Scholar
- Fire A, Xu S, Montgomery MK, Kostas SA, Driver SE, Mello CC: Potent and specific genetic interference by double-stranded RNA in Caenorhabditis elegans. Nature. 1998, 391: 806-811. 10.1038/35888PubMedGoogle Scholar
- Tabara H, Grishok A, Mello CC: RNAi in C. elegans: soaking in the genome sequence. Science. 1998, 282: 430-431. 10.1126/science.282.5388.430PubMedGoogle Scholar
- Timmons L, Fire A: Specific interference by ingested dsRNA. Nature. 1998, 395: 854- 10.1038/27579PubMedGoogle Scholar
- Aboobaker AA, Blaxter ML: Use of RNA interference to investigate gene function in the human filarial nematode parasite Brugia malayi. Mol Biochem Parasitol. 2003, 129: 41-51. 10.1016/S0166-6851(03)00092-6PubMedGoogle Scholar
- Sidow A, Thomas WK: A molecular evolutionary framework for eukaryotic model organisms. Curr Biol. 1994, 4: 596-603. 10.1016/S0960-9822(00)00131-7PubMedGoogle Scholar
- Urwin PE, Lilley CJ, Atkinson HJ: Ingestion of double-stranded RNA by preparasitic juvenile cyst nematodes leads to RNA interference. Mol Plant Microbe Interact. 2002, 15: 747-52. 10.1094/MPMI.2002.15.8.747PubMedGoogle Scholar
- Sukno SA, McCuiston J, Wong MY, Wang XH, Thon MR, Hussey R, Baum T, Davis E: Quantitative detection of double-stranded RNA-mediated gene silencing of parasitism genes in Heterodera glycines. J Nematol. 2007, 39: 145-152.PubMedPubMed CentralGoogle Scholar
- Ragsdale EJ, Crum J, Ellisman MH, Baldwin JG: Three-dimensional reconstruction of the stomatostylet and anterior epidermis in the nematode Aphelenchus avenae (Nematoda: Aphelenchidae) with implications for the evolution of plant parasitism. J Morphol. 2008, 269: 1181-1196. 10.1002/jmor.10651PubMedPubMed CentralGoogle Scholar
- Gao G, Raikar S, Davenport B, Mutapcic L, Montgomery R, Kuzmin E, Bennett KL: Cross-species RNAi: Selected Ascaris suum dsRNAs can sterilize Caenorhabditis elegans. Mol Biochem Parasitol. 2006, 146: 124-128. 10.1016/j.molbiopara.2005.11.003PubMedGoogle Scholar
- Shannon AJ, Tyson T, Dix I, Boyd J, Burnell AM: Systemic RNAi mediated gene silencing in the anhydrobiotic nematode Panagrolaimus superbus. BMC Mol Biol. 2008, 9: 58- 10.1186/1471-2199-9-58PubMedPubMed CentralGoogle Scholar
- Potts M, Slaughter SM, Hunneke F-U, Garst JF, Helm RF: Desiccation tolerance of prokaryotes: application of principles to human cells. Integr Comp Biol. 2005, 45: 800-809. 10.1093/icb/45.5.800..PubMedGoogle Scholar
- Roelofs D, Aarts MGM, Schat H, van Straalen NM: Functional ecological genomics to demonstrate general and specific responses to abiotic stress. Funct Ecol. 2008, 22: 8-18.Google Scholar
- Hundertmark M, Hincha DK: LEA (late embryogenesis abundant) proteins and their encoding genes in Arabidopsis thaliana. BMC Genomics. 2008, 9: 118- 10.1186/1471-2164-9-118PubMedPubMed CentralGoogle Scholar
- Bies-Ethève N, Gaubier-Comella P, Debures A, Lasserre E, Jobet E, Raynal M, Cooke R, Delseny M: Inventory, evolution and expression profiling diversity of the LEA (late embryogenesis abundant) protein gene family in Arabidopsis thaliana. Plant Mol Biol. 2008, 67: 107-124. 10.1007/s11103-008-9304-xPubMedGoogle Scholar
- Birrell GW, Brown JA, Wu HI, Giaever G, Chu AM, Davis RW, Brown JM: Transcriptional response of Saccharomyces cerevisiae to DNA-damaging agents does not identify the genes that protect against these agents. Proc Natl Acad Sci USA. 2002, 99: 8778-8783. 10.1073/pnas.132275199PubMedPubMed CentralGoogle Scholar
- Warringer J, Ericson E, Fernandez L, Nerman O, Blomberg A: High-resolution yeast phenomics resolves different physiological features in the saline response. Proc Natl Acad Sci USA. 2003, 100: 15724-15729. 10.1073/pnas.2435976100PubMedPubMed CentralGoogle Scholar
- Thorpe GW, Fong CS, Alic N, Higgins VJ, Dawes IW: Cells have distinct mechanisms to maintain protection against different reactive oxygen species: oxidative-stress-response genes. Proc Natl Acad Sci USA. 2004, 101: 6564-6569. 10.1073/pnas.0305888101PubMedPubMed CentralGoogle Scholar
- Tai SL, Snoek I, Luttik MAH, Almering MJH, Walsh MC, Pronk JT, Daran J-M: Correlation between transcript profiles and fitness of deletion mutants in anaerobic chemostat cultures of Saccharomyces cerevisiae. Microbiology. 2007, 153: 877-886. 10.1099/mic.0.2006/002873-0PubMedPubMed CentralGoogle Scholar
- Maeda I, Kohara Y, Masayuki M, Sugimoto A: Large-scale analysis of gene function in Caenorhabditis elegans by high-throughput RNAi. Curr Biol. 2001, 11: 171-176. 10.1016/S0960-9822(01)00052-5PubMedGoogle Scholar
- Kamath RS, Fraser AG, Dong Y, Poulin G, Durbin R, Gotta M, Kanapin A, Le Bot N, Moreno S, Sohrmann M, Welchman DP, Zipperlen P, Ahringer J: Systematic functional analysis of the Caenorhabditis elegans genome using RNAi. Nature. 2003, 421: 231-237. 10.1038/nature01278PubMedGoogle Scholar
- Knox DP, Geldhof P, Visser A, Britton C: RNA interference in parasitic nematodes of animals: a reality check?. Trends Parasitol. 2007, 23: 105-107. 10.1016/j.pt.2007.01.007PubMedGoogle Scholar
- Marzin S, Mihaly R, Pauk J, Schweizer P: A transient assay system for the assessment of cell-autonomous gene function in dehydration-stressed barley. J Exp Bot. 2008, 59: 3359-3369. 10.1093/jxb/ern186PubMedPubMed CentralGoogle Scholar
- Issa Z, Grant WN, Stasiuk S, Shoemaker CB: Development of methods for RNA interference in the sheep gastrointestinal parasite, Trichostrongylus colubriformis. Int J Parasitol. 2005, 35: 935-940. 10.1016/j.ijpara.2005.06.001PubMedGoogle Scholar
- Elbashir SM, Martinez J, Patkaniowska A, Lendeckel W, Tuschl T: Functional anatomy of siRNAs for mediating efficient RNAi in Drosophila melanogaster embryo lysate. EMBO J. 2001, 20: 6877-6888. 10.1093/emboj/20.23.6877PubMedPubMed CentralGoogle Scholar
- Furini A, Koncz C, Salamini F, Bartels D: High level transcription of a member of a repeated gene family confers dehydration tolerance to callus tissue of Craterostigma plantagineum. EMBO J. 1997, 16: 3599-3608. 10.1093/emboj/16.12.3599PubMedPubMed CentralGoogle Scholar
- Hilbricht T, Varotto S, Sgaramella V, Bartels D, Salamini F, Furini A: Retrotransposons and siRNA have a role in the evolution of desiccation tolerance leading to resurrection of the plant Craterostigma plantagineum. New Phytol. 2008, 179: 877-887. 10.1111/j.1469-8137.2008.02480.xPubMedGoogle Scholar
- Mathews DH, Sabina J, Zucker M, Turner H: Expanded sequence dependence of thermodynamic parameters provides robust prediction of RNA secondary structure. J Mol Biol. 1999, 288: 911-940. 10.1006/jmbi.1999.2700PubMedGoogle Scholar
- Rotruck JT, Pope AL, Ganther HE, Swanson AB, Hafeman DG, Hoekstra WG: Selenium: biochemical role as a component of glutathione peroxidase. Science. 1973, 179: 588-590. 10.1126/science.179.4073.588PubMedGoogle Scholar
- Kryukov GV, Castellano S, Novoselov SV, Lobanov AV, Zehtab O, Guigó R, Gladyshev VN: Characterization of mammalian selenoproteomes. Science. 2003, 300: 1439-1443. 10.1126/science.1083516PubMedGoogle Scholar
- Walczak R, Westhof E, Carbon P, Krol A: A novel RNA structural motif in the selenocysteine insertion element of eukaryotic selenoprotein mRNAs. RNA. 1996, 2: 367-379.PubMedPubMed CentralGoogle Scholar
- Passardi F, Theiler G, Zamocky M, Cosio C, Rouhier N, Teixera F, Margis-Pinheiro M, Ioannidis V, Penel C, Falquet L, Dunand C: PeroxiBase: the peroxidase database. Phytochemistry. 2007, 68: 1605-1611. 10.1016/j.phytochem.2007.04.005PubMedGoogle Scholar
- Leprince O, Atherton NM, Deltour R, Hendry G: The involvement of respiration in free radical processes during loss of desiccation tolerance in germinating Zea mays (an electron paramagnetic resonance study). Plant Physiol. 1994, 104: 1333-1339.PubMedPubMed CentralGoogle Scholar
- Gasch AP, Spellman PT, Kao CM, Carmel-Harel O, Eisen MB, Storz G, Botstein D, Brown PO: Genomic expression programs in the response of yeast cells to environmental changes. Mol Biol Cell. 2000, 11: 4241-4257.PubMedPubMed CentralGoogle Scholar
- Kültz D: Molecular and evolutionary basis of the cellular stress response. Ann Rev Physiol. 2005, 67: 225-257. 10.1146/annurev.physiol.67.040403.103635..Google Scholar
- Leprince O, Hendry GAF, McKersie BD: The mechanisms of desiccation tolerance in developing seeds. Seed Sci Res. 1993, 3: 231-246. 10.1017/S0960258500001859..Google Scholar
- Kranner I, Beckett RP, Wornik S, Zorn M, Pfeifhofer HW: Revival of a resurrection plant correlates with its antioxidant status. Plant J. 2002, 31: 13-24. 10.1046/j.1365-313X.2002.01329.xPubMedGoogle Scholar
- França MB, Panek AD, Eleutherio ECA: Oxidative stress and its effects during dehydration. Comp Biochem Physiol A. 2007, 146: 621-631. 10.1016/j.cbpa.2006.02.030..Google Scholar
- Moore JP, Tuan Le N, Brandt WF, Driouich A, Farrant JM: Towards a systems-based understanding of plant desiccation tolerance. Trends Plant Sci. 2008, 14: 110-117. 10.1016/j.tplants.2008.11.007..Google Scholar
- Ingle RA, Schmidt UG, Farrant JM, Thomson JA, Mundree SG: Proteomic analysis of leaf proteins during dehydration of the resurrection plant Xerophyta viscosa. Plant Cell Environ. 2007, 30: 435-446. 10.1111/j.1365-3040.2006.01631.xPubMedGoogle Scholar
- Jiang G, Wang Z, Shang H, Yang W, Hu Z, Phillips J, Deng X: Proteome analysis of leaves from the resurrection plant Boea hygrometrica in response to dehydration and rehydration. Planta. 2007, 225: 1405-1420. 10.1007/s00425-006-0449-zPubMedGoogle Scholar
- Wharton DA, Marshall CJ: How do terrestrial Antarctic organisms survive in their harsh environment?. J Biol. 2009, 8: 39-PubMedPubMed CentralGoogle Scholar
- Fisher JM, Evans AAF: Penetration and feeding by Aphelenchus avenae. Nematologica. 1967, 13: 425-428.Google Scholar
- Winston PW, Bates DH: Saturated solutions for the control of humidity in biological research. Ecology. 1960, 41: 232-237. 10.2307/1931961..Google Scholar
- Eisenberg D: Three-dimensional structure of membrane and surface proteins. Ann Rev Biochem. 1984, 53: 595-623. 10.1146/annurev.bi.53.070184.003115PubMedGoogle Scholar
- Kyte J, Doolittle RF: A simple method for displaying the hydropathic character of a protein. J Mol Biol. 1982, 157: 105-132. 10.1016/0022-2836(82)90515-0PubMedGoogle Scholar
- Prilusky J, Felder CE, Zeev-Ben-Mordehai T, Rydberg E, Man O, Beckmann JS, Silman I, Sussman JL: FoldIndex©: a simple tool to predict whether a given protein sequence is intrinsically unfolded. Bioinformatics. 2005, 21: 3435-3438. 10.1093/bioinformatics/bti537PubMedGoogle Scholar
- Li X, Romero P, Rani M, Dunker AK, Obradovic Z: Predicting protein disorder for N-, C-, and internal regions. Genome Informat. 1999, 10: 30-40.Google Scholar
- Dosztányi Z, Csizmók V, Tompa P, Simon I: The pairwise energy content estimated from amino acid composition discriminates between folded and intrinsically unstructured proteins. J Mol Biol. 2005, 347: 827-839. 10.1016/j.jmb.2005.01.071PubMedGoogle Scholar
- Uversky VN, Gillespie JR, Fink AL: Why are "natively unfolded" proteins unstructured under physiologic conditions?. Proteins Struct Funct Genet. 2000, 41: 415-427. 10.1002/1097-0134(20001115)41:3<415::AID-PROT130>3.0.CO;2-7PubMedGoogle Scholar
- Pearson WR: Rapid and sensitive sequence comparison with FASTP and FASTA. Methods Enzymol. 1990, 183: 63-98. full_textPubMedGoogle 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.