Correction
14 Sep 2023: The PLOS Genetics Staff (2023) Correction: A novel nematode species from the Siberian permafrost shares adaptive mechanisms for cryptobiotic survival with C. elegans dauer larva. PLOS Genetics 19(9): e1010943. https://doi.org/10.1371/journal.pgen.1010943 View correction
Figures
Abstract
Some organisms in nature have developed the ability to enter a state of suspended metabolism called cryptobiosis when environmental conditions are unfavorable. This state-transition requires execution of a combination of genetic and biochemical pathways that enable the organism to survive for prolonged periods. Recently, nematode individuals have been reanimated from Siberian permafrost after remaining in cryptobiosis. Preliminary analysis indicates that these nematodes belong to the genera Panagrolaimus and Plectus. Here, we present precise radiocarbon dating indicating that the Panagrolaimus individuals have remained in cryptobiosis since the late Pleistocene (~46,000 years). Phylogenetic inference based on our genome assembly and a detailed morphological analysis demonstrate that they belong to an undescribed species, which we named Panagrolaimus kolymaensis. Comparative genome analysis revealed that the molecular toolkit for cryptobiosis in P. kolymaensis and in C. elegans is partly orthologous. We show that biochemical mechanisms employed by these two species to survive desiccation and freezing under laboratory conditions are similar. Our experimental evidence also reveals that C. elegans dauer larvae can remain viable for longer periods in suspended animation than previously reported. Altogether, our findings demonstrate that nematodes evolved mechanisms potentially allowing them to suspend life over geological time scales.
Author summary
Survival in extreme environments for prolonged periods is a challenge that only a few organisms, are capable of. It is not well understood, which molecular and biochemical pathways are utilized by such cryptobiotic organisms, and how long they might suspend life. Here, we show that a soil nematode Panagrolaimus kolymaensis, suspended life for 46,000 years in the Siberian permafrost. Through comparative analysis, we find that P. kolymaensis and model organism C. elegans utilize similar adaptive mechanisms to survive harsh environmental conditions for prolonged periods. Our findings here are important for the understanding of evolutionary processes because generation times could be stretched from days to millennia, and long-term survival of individuals of species can lead to the refoundation of otherwise extinct lineages.
Citation: Shatilovich A, Gade VR, Pippel M, Hoffmeyer TT, Tchesunov AV, Stevens L, et al. (2023) A novel nematode species from the Siberian permafrost shares adaptive mechanisms for cryptobiotic survival with C. elegans dauer larva. PLoS Genet 19(7): e1010798. https://doi.org/10.1371/journal.pgen.1010798
Editor: Gregory P. Copenhaver, The University of North Carolina at Chapel Hill, UNITED STATES
Received: February 17, 2023; Accepted: May 24, 2023; Published: July 27, 2023
Copyright: © 2023 Shatilovich et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All the data files are available at https://doi.org/10.5281/zenodo.6590382.
Funding: This work was supported by the Russian Foundation for Basic Research (19-29-05003-mk) to AS and ER. VRG and TVK acknowledge the financial support from the Volkswagen Foundation (Life research grant 92847). PHS and TTH are supported by a DFG ENP grant to PHS (DFG project 434028868). GMH is funded by a UCD Ad Astra Fellowship. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Organisms from diverse taxonomic groups can survive extreme environmental conditions, such as the complete absence of water or oxygen, high temperature, freezing, or extreme salinity. The survival strategies of such organisms include a state known as suspended animation or cryptobiosis, in which they reduce metabolism to an undetectable level [1]. Spectacular examples of long-term cryptobiosis include a Bacillus spore that was preserved in the abdomen of bees buried in amber for 25 to 40 million years [2], and a 1000 to 1500 years-old Lotus seed, found in an ancient lake, that was subsequently able to germinate [3]. Metazoans such as tardigrades, rotifers, and nematodes are also known for remaining in cryptobiosis for prolonged periods [4,5]. The longest records of cryptobiosis in nematodes are reported for the Antarctic species Plectus murrayi [6] (25.5 years in moss frozen at -20°C), and Tylenchus polyhypnus [7] (39 years desiccated in an herbarium specimen).
Intensive research during the last decade has demonstrated that permafrost (perennially frozen sediments) are unique ecosystems preserving life forms at sub-zero temperatures over thousands of years [8,9,10,11]. Permafrost remains are an exceptional source for discovering a wide variety of unicellular and multicellular living organisms surviving in cryptobiosis for prolonged periods [1,12,13]. The Siberian permafrost is a unique repository for preserving organisms in sub-zero temperatures for millions of years. Expeditions in the past decade have resulted in the revival of several organisms across various taxa from the Siberian permafrost [14,15,16,17]. The possibility to exploit permafrost as a source for reanimating multicellular animals was recognized as early as 1936. A viable Cladocera crustacean, Chydorus sphaericus, preserved in the Transbaikalian permafrost for several thousand years [18,19], was discovered by P. N. Kapterev, who worked at the scientific station Skovorodino as a GULAG prisoner. Unfortunately, this observation remained unnoticed for many decades. We recently reanimated soil nematodes that were preserved in Siberian permafrost for potentially thousands of years, and initial morphological observations provisionally described them as belonging to the genera Panagrolaimus and Plectus. Previous studies demonstrated several species of Panagrolaimus can undergo cryptobiosis in the form of anhydrobiosis (through desiccation) and cryobiosis (through freezing) [20,21,22,23,24]. In various nematodes, entry into anhydrobiosis is often accompanied by a preparatory phase of exposure to mild desiccation, known as preconditioning [22,25]. This induces a specific re-modelling of the transcriptome, the proteome, and metabolic pathways that enhances survival ability [26,27,28]. Some panagrolaimids possess adaptive mechanisms for rapid desiccation where most of the cellular water is lost, while others possess freezing tolerance without loss of water at sub-zero temperatures by inhibiting the growth and recrystallisation of ice crystals [22].
Here, we present a high-quality genome assembly, detailed morphological phylogenetic analysis, and define a novel species, Panagrolaimus kolymaensis. Precise radiocarbon dating indicates that P. kolymaensis remained in cryptobiosis for about 46,000 years, since the late Pleistocene. Making use of the model organism C. elegans, we demonstrate that C. elegans dauer larvae and Panagrolaimus kolymaensis utilize comparable molecular mechanisms to survive extreme desiccation and freezing, i.e. upregulation of trehalose biosynthesis and gluconeogenesis.
Results
Taxonomy
Panagrolaimus kolymaensis Shatilovich, Gade, Pippel, Hoffmeyer, Tchesunov, Stevens, Winkler, Hughes, Traikov, Hiller, Rivkina, Schiffer, Myers & Kurzchalia sp. nov.
urn:lsid:zoobank.org:act:57A9E39B-5603-46B6-A035-4B9BDBC1C441
Discovery site and radiocarbon dating
Previously, we had shown that nematodes from the Siberian permafrost with morphologies consistent with the genera Panagrolaimus and Plectus could be reanimated thousands of years after they had been frozen. Several viable nematode individuals were found in two of the more than 300 studied samples of permafrost deposits spanning different ages and genesis. Samples were collected by researchers of the Soil Cryology Lab, Pushchino, Russia, during perennial paleo-ecological expeditions carried out in the coastal sector of the northeastern Arctic [12]. The detailed description of the study site (outcrop Duvanny Yar, Kolyma River, Fig 1A), sampling and revitalizing procedures are provided in S1 Text. Like other late Pleistocene permafrost formations in the northeastern Arctic, Duvanny Yar is composed of permanently frozen ice-rich silt deposits riddled with large polygonal ice wedges that divide them into mineral blocks [29,30] (Fig 1B). Sediments include sandy alluvial layers, peat lenses, buried paleosols and Pleistocene rodent burrows (Fig 1C). The burrow (P-1320), in which Panagrolaimus nematodes were found (Fig 1D), has been taken from the frozen outcrop wall at a depth of about 40 m below the surface and about 11 m above river water level in undisturbed and never thawed late Pleistocene permafrost deposits. The fossil burrow left by arctic gophers of the genus Citellus consists of an entrance tunnel and large nesting chamber up to 25 cm in diameter [29].
(A) Location of the Duvanny Yar outcrop on the Kolyma River, northeastern Siberia, Russia. https://climate.copernicus.eu/sites/default/files/inline-images/C3S_indicator_sea_ice_sidebar_figS2_branded.png) (https://climate.copernicus.eu/data-protection-and-privacy-statement) (B) View of the upper part of outcrop composed of ice wedges and permafrost silty deposits. (C) Lithostratigraphic scheme of deposits, showing location of studied rodent borrow (red circle). (D) fossil rodent burrow with herbaceous litter and seeds buried in permafrost deposits; m a.r.l. = meters above river level.
https://doi.org/10.1371/journal.pgen.1010798.g001
The sterility of permafrost sampling and age of cultivated biota have been discussed in detail in several reviews [9,31,32]. Based on previous reports, the age of the organisms found in a burrow is equal to the freezing time and corresponds to the age of organic matter conserved in the syncryogenic sediments. This makes it possible to use radiocarbon dating of organic matter to establish the age of organisms. We performed Accelerator Mass Spectrometry (AMS) radiocarbon analysis of plant material obtained from studied borrow P-1320 and determined a direct 14C age of 44,315±405 BP (Institute of Geography, RAS; sample IGANAMS 9137). Calibrated age range is 45,839–47,769 cal BP (95.4% probability) (S1 Fig).
Like other parthenogenetic Panagrolaimus, the newly discovered species is triploid
The revived animal was cultivated in the laboratory for over 100 generations and initially described as Panagrolaimus aff. detritophagus [33] based on morphology. We conducted a detailed morphological analysis of the revived animal (Figs 2 and S2, and Table A in S1 Text and Box 1), which confirmed unambiguously that the animal belongs to the genus of Panagrolaimus, in agreement with a previous phylogenetic analysis of the 18S ribosomal RNA sequence [33]. However, due to the morphological uniformity of Panagrolaimus, unusual even for nematodes, morphology and molecular analysis of a single ribosomal RNA sequence is insufficient to describe a species. We found the species to be parthenogenetic, which further complicates description under most species concepts. Due to these limitations, we decided to refer to the phylogenetic species concept, using phylogenetic trees based on multiple genes as markers.
Scanning electron pictures (A, C), light microscopy photographs (E, F) and graphic presentations (B, D, G) of holotype: A, B) entire body, C, D) anterior ends, E) anterior body, F) perivulvar body region, G) tail. Abbreviations: l.f.–lateral field, ov–ovary, pro–procorpus of the pharynx, t.b.–terminal bulb of the pharynx, u–uterus with eggs, v–vulva, v.p.–ventral pore. Scale bars: A, D, E, F, G– 20 μm, B– 100 μm, C– 2 μm.
https://doi.org/10.1371/journal.pgen.1010798.g002
To obtain comprehensive molecular data for phylogenomic species determination, we generated a genome assembly using PacBio HiFi sequencing with long reads (84X coverage, mean length 14,425 bp). Our analysis of repeat and gene content is described in Table C in S1 Text. K-mer analysis of the reads clearly indicated that this animal has a triploid genome (Fig 3A), like other parthenogenetic Panagrolaimus species [34]. Despite the challenges that a triploid genome poses for assembly, we obtained a highly contiguous contig assembly of the three pseudohaplotypes that comprise almost 266 Mb and thus have a similar genome size as other parthenogenetic Panagrolaimus species [34]. The contig N50 value of all three pseudohaplotypes is 3.8 Mb. Since these pseudohaplotypes exhibited a noticeable degree of divergence, we further investigated their relationship by using the apparent homeologs in our gene predictions to align the longest continuous contigs based on micro synteny (Fig 3B). Links between the contigs clearly show the triploid state of the genome.
A) Kmer spectra of the P. kolymaensis PacBio HiFi data. Kmers of length 19 were counted using Jellyfish. B) Circos plot showing the triploid structure of the P. kolymaensis genome. Lines represent the position of 6,715 homeologs in eight contigs that comprise 39.9 Mb (15%) of the assembly. Homeologs were identified by clustering protein-coding genes into orthogroups using OrthoFinder and selecting groups containing three sequences. Contig IDs and scale is shown. C) Inferred species tree for all taxa. The maximum likelihood tree inferred using a concatenated supermatrix (18S and 28S genes) with bootstrap support values is displayed. All genera are represented as monophyletic clades. P. kolymaensis is highlighted in red and basal to all other Panagrolaimus taxa. Internal nodes, where all subsequent branches represent identical sequences, are displayed with a black star. D) P. kolymaensis possesses C. elegans gene orthologs to enzymes required for TCA cycle, glyoxylate shunt, glycolysis, gluconeogenesis, trehalose synthesis, and polyamine synthesis. Black filled circles: Ortholog presence suggested by orthogroup clustering, phylogenetic analysis, and domain architecture. White filled circles: No ortholog found via current analysis. Coloured filled circles: Presence of P. kolymaensis gene(s), related to several C. elegans genes (all genes of same colour) that are all co-orthologous to that gene (those genes). Label: C. elegans enzyme names and orthogroup that contains that gene according to our orthogroup clustering.
https://doi.org/10.1371/journal.pgen.1010798.g003
Applying a phylogenomic concept to define P. kolymaensis
To place the species in the genus Panagrolaimus, we conducted a broader multi-gene phylogenomic analysis using Maximum likelihood methods. Our analysis of a concatenated, partitioned alignment of 60 genes, and a coalescence-based approach using a broader set of 12,295 gene trees, retrieves the revived animal as sister to all other sequenced Panagrolaimus species, but as an ingroup to Propanagrolaimus [35] (Figs 3C, S3B, and S3C). Thus, the phylogenetic placement provides strong evidence that this animal represents a novel species. Furthermore, there is substantial sequence divergence between this novel species, and Panagrolaimus sp. PS1159 and Panagrolaimus sp. ES5, estimated to be on average 2.06 and 2.11 amino acid substitutions per site in our concatenated alignment, respectively. The substantial divergence is in line with previous data on ages of Panagrolaimus nematodes [34], and more broadly seen in nematodes, which can be hyper-diverse [36,37]. Our data also contradicts the assumption that all parthenogenetic strains in the genus Panagrolaimus have a monophyletic origin [34] (Fig 3C). Using the gene-tree reconciliation approach implemented in GRAMPA (Gene-tree Reconciliation Algorithm with MUL (Multi labelled)-trees for Polyploid Analysis), we explored whether the additional set of proteins we found in the parthenogenetic species is a result of auto- or allopolyploidy. Finding the additional proteins basally branching, outside of the lineage containing both parthenogenetic and sexual species, suggests that an allopolyploid origin of these extra gene-copies (S6 Fig).
Based on the Kolyma River location where the animal was unearthed, we propose the following taxonomic classification and species name:
- Phylum Nematoda Potts, 1932
- Class Chromadorea Inglis 1983
- Suborder Tylenchina Thorne, 1949
- Family Panagrolaimidae Thorne, 1937
- Panagrolaimus kolymaensis
C. elegans dauer larvae and P. kolymaensis might utilize partially similar mechanisms to enter and remain in cryptobiotic state for prolonged periods of time
In the absence of established genetic methods in P. kolymaensis, we referred to the model C. elegans as a comparator system to gain insights into possible pathways for long term survival [25,27,28,33]. The high-quality genome of P. kolymaensis allowed us to compare its molecular toolkit for cryptobiosis with that of C. elegans. We used orthology clustering and phylogenetics to investigate whether the genome of P. kolymaensis contains genes previously implicated in cryptobiosis in the C. elegans dauer larva. Our analysis showed that, like other Panagrolaimus species [38,39], P. kolymaensis also encodes orthologs to a C. elegans trehalose phosphate synthase gene (tps-2) and to a trehalose phosphatase gene (gob-1) (Fig 3D and S1 Orthology analysis). Furthermore, we found orthologs to all C. elegans enzymes required for polyamine biosynthesis, the TCA cycle, glycolysis, gluconeogenesis, and glyoxylate shunt (Fig 3D and S1 Orthology analysis) suggesting that P. kolymaensis might partially utilize similar molecular mechanisms as C. elegans to facilitate survival of unfavorable conditions.
Our earlier findings established that amongst several developmental stages of C. elegans, only the dauer larva, formed during unfavorable conditions (such as low nutrients and high population density), could survive anhydrobiosis and exposure to freezing [25,33]. The dauer larva is in a hypometabolic state with distinct metabolic properties such as reduced oxygen consumption and heat dissipation in comparison to other larval stages of C. elegans. To survive extreme desiccation, C. elegans dauer larvae (in its hypometabolic state) need to be first preconditioned at high relative humidity (98% RH) for 4 days [25]. During preconditioning, dauer larvae upregulate trehalose biosynthesis that ensures their survival to harsh desiccation [25,28]. We tested whether the survival of P. kolymaensis is also facilitated by preconditioning. As there is no dauer stage in the Panagrolaimus life cycle, we performed our experiments with a population of all the larval stages and adults. Although a small proportion of P. kolymaensis individuals survive harsh desiccation and freezing without preconditioning (Fig 4A), the mixture of all the larval stages and adults of P. kolymaensis survive significantly higher (p value< 0.0001) in proportion to harsh desiccation upon preconditioning (Fig 4A). Similarly, preconditioning and desiccation further enhanced survival rate of P. kolymaensis to freezing (-80°C) (Fig 4A). Like C. elegans dauer larva, P. kolymaensis upregulates trehalose levels up to 20-fold upon preconditioning (Fig 4B). We previously reported that to upregulate trehalose levels upon preconditioning, C. elegans dauer larva dissipate their fat reserves (Triacylglycerols) by activating the glyoxylate shunt and gluconeogenic pathway [28] Upon preconditioning, we found that triacylglyceride (TAG) levels are significantly decreased in P. kolymaensis (S5A and S5B Fig). To further investigate whether acetyl-CoA derived from the degradation of TAGs culminates in trehalose, we applied the previously developed method of metabolic labelling with 14C-acetate in combination with 2D-TLC [1,33]. The 14C-acetate metabolized by the worms is incorporated into TAGs. Upon degradation of TAGs, 14C-acetyl CoA is released which acts as a precursor for trehalose biosynthesis. As shown in Fig 4C, preconditioning led to a huge increase of radioactivity in trehalose and to a small increase in some amino acids (glycine/serine, phenylalanine; panels C and D). Interestingly, P. kolymaensis displayed an additional spot (Fig 4D, enumerated as 7), that was not found in C. elegans. We identified this spot as trehalose-6-phosphate (S5C–S5H Fig), a precursor of trehalose, based on the fragmentation pattern of the molecule, using mass spectrometry. Thus, to resist harsh desiccation, like C. elegans dauer larvae, P. kolymaensis might utilize the glyoxylate shunt and consequently acetate derived from TAGs to synthesize trehalose. Detection of the immediate precursor (trehalose-6-phosphate) suggests that the flux of metabolites is intense in P. kolymaensis. Finally, we investigated whether C. elegans dauer larvae can also survive in prolonged cryptobiotic state. Despite preconditioning, the survival ability of desiccated dauer larvae at room temperature declines very rapidly, with most larvae dead after almost 10 days (Fig 4E) (S4A and S4B Fig). Direct freezing without any cryoprotectants at -80 °C leads to instant death of the animals. To test whether combining these conditions could extend the viability of dauer larvae (S4A and S4B Fig), we transferred the desiccated larvae to -80 °C. Remarkably, under these conditions, there was no significant decline in viability even after 480 days (Fig 4E). Moreover, after thawing the animals resumed reproductive growth and produced progeny in numbers like those of animals kept under control conditions (Fig 4F). Since we did not observe any reduction in the survival at any time points, this suggests that the combination of anhydrobiosis and freezing can prolong the survival ability of dauer larvae. Thus, C. elegans dauer larvae, when exposed to combination of cryptobiotic states can survive for extremely long periods of time.
A) Survival rate of P. kolymaensis nematodes to desiccation and freezing (-80°C). Error bars indicate standard error of mean of two independent experiments with two technical replicates performed on two different days. Statistical comparison was performed using unpaired t-test with Welch correction. n.s. p > 0.05, ****p < 0.0001. For desiccation (non-preconditioned) n = 289, freezing (non-preconditioned) n = 675, desiccation (preconditioned at 98%RH) n = 953 and freezing (preconditioned at 98%RH) n = 1295. B) P. kolymaensis nematodes and daf-2(e1370) dauer larvae upregulate trehalose levels upon preconditioning at 98%RH. Error bars indicate standard error of mean of two independent experiments with three replicates performed on two different days. Statistical comparison was performed using two-way ANOVA with Holm-Sidak’s multiple comparison test, ****p < 0.0001. C-D) 2D-thin layer chromatography of 14C-acetate labelled metabolites from P. kolymaensis that were non-preconditioned and preconditioned at 98%RH. Enumerated spots indicate trehalose (1), glucose (2), glutamate (3), glutamine (4), serine/glycine (5) and phenylalanine (6). Representative images from at least two independent experiments performed on two different days. E) Desiccated daf-2 (e1370) dauer larvae survive to freezing (-80°C) for an extremely long period. Error bars indicate standard error of mean of two independent experiments with two technical replicates performed on two different days. F) Brood size of desiccated dauer larvae exposed to freezing remain like that non- desiccated dauer larvae. Average brood size is the mean of seven dauer larvae per each condition. Statistical comparison was performed by using non-parametric Kolmogorov-smirnov test n.s p>0.05.
https://doi.org/10.1371/journal.pgen.1010798.g004
Discussion
The new nematode species from permafrost can now be placed into the genus Panagrolaimus [40], which contains several described parthenogenetic and gonochoristic species [34,41]. Many Panagrolaimus display adaptation to survival in harsh environments [22] and the genus includes the Antarctic species P. davidi [23]. The genus Panagrolaimus is exceptional in its morphological uniformity even among nematode species that are hard to classify based on morphology in general. Thus, species designation via microscopic (including SEM) analysis is unreliable, which is further complicated by the absence of males in parthenogenetic species. Males have an important diagnostic feature such as spicules and pericloacal papillae, females differ from one species to another mainly by morphometrics, where interspecies differences (absolute measures and ratios) might be subtle. Our specimens are similar based on absolute sizes and ratios to females of the bisexual species Panagrolaimus detritophagus [42]. The only non-overlapping morphometric character is index “b” (body length: pharynx length): 5.6–6.8 in P. kolymaensis versus 4.4–5.1 in P. detritophagus.
Consequently, we turned to phylogenomic methods under the phylogenetic species concept to place the species on the tree. This showed that this species is an outgroup to other known Panagrolaimus species, raising the possibility of a second independent evolution of parthenogenesis in the genus, in contrast to previous findings [34,35,41]. Alternatively, the hybrid origin of parthenogenetic Panagrolaimus could influence the phylogenetic positioning of strains, raising the possibility that the new species is a true sister to the other parthenogenetic strains. To fully resolve the phylogenetic positioning further, extensive sampling, and genome sequencing of Panagrolaimus species is needed. We found P. kolymaensis to be triploid and thus a hybrid origin is possible, as seen in other parthenogenetic Pangrolaimus [34].
The highly contiguous genome of P. kolymaensis will allow for analyses of this trait in comparison to other Panagrolaimus species currently being genome sequenced. Our results provide a deeper insight into the homology of molecular and biochemical mechanisms between C. elegans and P. kolymaensis, which are not only taxonomically but also ecologically distinct. C. elegans can mostly be found in rotting fruits and plants in temperate regions [43,44], while Panagrolaimus species are globally distributed and prevalent in leaf litter and soil [41], including in harsh environments [22].
We show through orthology analysis that the well-studied molecular pathways used by C. elegans larvae to enter the dauer state, such as insulin [45,46] (DAF-11, DAF-2 & DAF-16), TGF-β [47] (DAF-7), steroid [48] (DAF-9, DAF-12) are present in the genome of the P. kolymaensis (S4C Fig). The presence of homologous genes in two species does not necessarily demonstrate their functionality in both. Therefore, further functional analyses are needed to study molecular pathways in detail. Trehalose accumulation (Fig 4B) and depletion of triacylglycerols (S5A and S5B Fig) ensure the functionality of the trehalose biosynthesis pathway and utilization of glyoxylate shunt during desiccation in P. kolymaensis. Without the activity of the enzyme TPS-2 and glyoxylate shunt, it is unfeasible to synthesize trehalose in nematodes. We do not eliminate the possibility of other biochemical features that might contribute to the desiccation survival ability of P. kolymaensis, but with regards to trehalose biosynthesis and the glyoxylate shunt, our data suggest that the molecular tool kit is partially orthologous. In our future studies, we intend to perform RNAi-based experiments to infer the concrete mechanisms. Our results hint at convergence or parallelism in the molecular mechanisms organizing dauer formation and cryptobiosis.
As mentioned above, preconditioning enhances the survival of P. kolymaensis by rendering them desiccation tolerant. We previously reported that preconditioning elevates trehalose biosynthesis in C. elegans dauer larvae and the elevated trehalose renders desiccation tolerance by protecting the cellular membranes [25]. It is not surprising that P. kolymaensis upregulates trehalose, however the magnitude of trehalose elevation is higher than C. elegans dauer larvae. This indicates that central regulators (DAF-16, DAF-12) of trehalose upregulation may differentially regulate tps-2 in P. kolymaensis [38,49,50]. Although P. kolymaensis utilizes the glyoxylate shunt and gluconeogenesis to upregulate trehalose levels, it is intriguing to observe that they accumulate substantial levels of trehalose-6-phosphate. Further investigation of this observation using RNAi or inhibitor-based experiments will provide insights into molecular mechanisms of metabolic regulation in P. kolymaensis upon preconditioning. Our findings for the first time demonstrate that C. elegans dauer larvae possess an inherent ability to survive freezing for prolonged periods if they undergo anhydrobiosis. It is tempting to speculate that undergoing anhydrobiosis might be a survival strategy of C. elegans to survive the seasonal changes in nature.
In summary, our findings indicate that by adapting to survive cryptobiotic state for short time frames in environments like permafrost, some nematode species gained the potential for individual worms to remain in the state for geological timeframes. This raises the question of whether there is an upper limit to the length of time an individual can remain in the cryptobiotic state. Long timespans may be limited only by drastic changes to the environment such as strong fluctuations in ambient temperature, natural radioactivity, or other abiotic factors. These findings have implications for our understanding of evolutionary processes, as generation times may be stretched from days to millennia, and long-term survival of individuals of species can lead to the refoundation of otherwise extinct lineages. This is particularly interesting in the case of parthenogenetic species, as each individual can find a new population without the need for mate finding, i.e. evading the cost of sex. Finally, understanding the precise mechanisms of long-term cryptobiosis and cues that lead to successful revivals can inform new methods for long term storage of cells and tissues.
Methods
Nomenclatural Acts
The electronic edition of this article conforms to the requirements of the amended International Code of Zoological Nomenclature, and hence the new names contained herein are available under that Code from the electronic edition of this article. This published work and the nomenclatural acts it contains have been registered in ZooBank, the online registration system for the ICZN. The ZooBank LSIDs (Life Science Identifiers) can be resolved and the associated information viewed through any standard web browser by appending the LSID to the prefix “http://zoobank.org/”. The LSID for this publication is: urn:lsid:zoobank.org:pub:150A5436-941D-4244-8884-89E39070BFE5. The electronic edition of this work was published in a journal with an ISSN, and has been archived and is available from the following digital repositories: LOCKSS.
Materials and C. elegans strains
[1-14C] -acetate (sodium salt) was purchased from Hartmann Analytic (Braunschweig, Germany). All other chemicals were purchased from Sigma-Aldrich (Taufkirchen, Germany). The Caenorhabditis Genetic Centre (CGC) which is funded by NIH Office of Research Infrastructure Programs (P40 OD010440) provided daf-2(e1370) and E. coli NA22 strains.
Genomic DNA isolation from P. kolymaensis nematodes
After isolation (S1 Text), to ensure our strain P. kolymaensis (Pn2-1) can adapt to different laboratories, we grew them for multiple generations. The strain was in culture for several generations during genomic DNA isolation and was frozen after genomic DNA isolation was performed. P. kolymaensis nematodes (isofemale strain Pn2-1) were grown on several plates of NGM agar plated with E. coli NA22 bacteria at 20°C. Worms were collected from the plates, washed with water at least three to five times by centrifugation at 1000 g to remove any residual bacteria and any debris. The worm pellet was dissolved in 5 volumes of worm lysis buffer (0.1M Tris-HCl pH = 8.5, 0.1M NaCl, 50mM EDTA pH = 8.0) and distributed in 1.5 ml of microcentrifuge tubes. These tubes are incubated at -80°C for 20 minutes. 100 μl of Proteinase ‘K’ (20 mg/ml) was added to each tube and they are incubated at 60°C overnight. 625 μl of cold GTC buffer (4M Guanidinium Thiocynate, 25mM Sodium citrate, 0.5% (v\v) N-lauroylsarcosine, 7%(v/v) Beta Mercaptoethanol) was added to the tube, incubated on ice 30 min, and mixed by inverting every 10 min. 1 volume of phenol–chloroform-isoamyl alcohol (pH = 8) was added to the lysate and mixed by inverting the tube 10–15 times. Tubes were centrifuged for 5 min at 10,000 g at 4°C to separate the phases. The upper aqueous phase was carefully collected into a fresh tube. One volume of fresh chloroform was added and mixed by inverting the tubes for 10–15 times and centrifuged for 5 min at 10,000 g at 4°C to separate the phases. One volume of cold 5 M NaCl was added, mixed by inverting the tubes and incubated on ice for 15 min. After incubation these tubes were centrifuged for 15 min at 12,000–16,000 g at 4°C. The supernatant containing the nucleic acids were slowly transferred into a fresh tube. One volume of isopropanol was added to the tube, inverted few times, and incubated on ice for 30 minutes. After incubation, the tubes were centrifuged at 3000 g for 30–45 min at 25°C and the supernatant was discarded without disturbing the pellet. The pellet was washed twice with 1 ml of 70% ethanol, tubes were centrifuged at 3000 g for 5 min and supernatant was discarded and incubated at 37°C for 10–15 min to dry the pellet. The pellet was resuspended carefully in TE buffer. The quality of the genomic DNA was analyzed with pulse field gel electrophoresis.
Genome sequencing and assembly
The long insert library was prepared as recommended by Pacific Biosciences according to the ‘Procedure & Checklist-Preparing gDNA Libraries Using the SMRTbellExpress Template Preparation Kit 2.0’ protocol. In summary, RNAse treated HMW gDNA was sheared to 20 kb fragments on the MegaRuptordevice (Diagenode) and 10 μg sheared gDNA was used for library preparation. The PacBio SMRTbell library was size selected in two fractions (9-13kb, > 13kb) using the BluePippin device with cassette definition of 0.75% DF MarkerS1 3–10 kb Improved Recovery. The second fraction of the size-selected library was loaded with 95 pM on plate on a Sequel SMRT cell (8M). Sequel polymerase 2.0 was used in combination with the v2 PacBio sequencing primer and the Sequel sequencing kit 2.0EA, with a runtime of 30 hours. We created PacBio CCS reads from the subreads.bam file using PacBio’s ccs command linetool (version4.2.0), outputting 8.5Gb of high-quality CCS reads (HiFi reads N50 of 14.4 kb). HiCanu (version 2.2) [51] was used to create the contig assembly. Blobtools [52] (version 1.1.1) was used to identify and remove bacterial contigs. The final triploid contig assembly consists of 856 contigs has a N50 of 3.82 Mb and a size of 266Mb. The mitochondrial genome was created with the mitoHifi pipeline (version 2, https://github.com/marcelauliano/MitoHiFi based on the assembled contigs and the closely related reference mitochondrial genome of Panagrellus redivivus (strain: PS2298/MT8872, ENAaccession: AP017464). The mitoHifi pipeline identified 49 mitochondrial contigs ranging from 13-32Kb. The final annotated circular mitochondrial genome has a length of 17467 bp.
To identify pseudohaplotypes in the P. kolymaensis genome assembly, we selected the longest isoform of each predicted protein-coding gene in our assembly and in the C. elegans genome (downloaded from WormBase Parasite, release WBPS15) using AGAT (version 0.4.0) and clustered them into orthologous groups (OGs) using OrthoFinder (version 2.5.2). We identified OGs that contained three Panagrolaimus sequences (i.e. groups that were present as single-copy in all three pseudohaplotypes) and used these to identify trios of multi-megabase size contigs derived from the three pseudohaplotypes. Synteny between the three pseudohaplotypes was visualized using Circos to plot the positions of each homeolog (version 0.69–8).
Genome annotation
RepeatModeler 1.0.8 (http://www.repeatmasker.org/) was used with parameter ‘-engine ncbi’ to create a library of repeat families which was used with RepeatMasker 4.0.9 to soft-mask the Panagrolaimus genome. To annotate genes, we cross mapped protein models from an existing Panagrolaimus as external evidence in the Augustus based pipeline. The completeness of our predictions was evaluated using BUSCO on the gVolante web interface.
Orthology analysis
We conducted a gene orthology analysis using genomic data from P. kolymaensis, the Plectid nematode species from the permafrost, as well as genomic data from WormBase Parasite (https://parasite.wormbase.org; accessed 17/12/2020): Caenorhabditis elegans, Diploscapter coronatus, Diploscapter pachys, Halicephalobus mephisto, Panagrellus redivivus, Panagrolaimus davidi, Panagrolaimus sp. ES5, Panagrolaimus sp. PS1159, Panagrolaimus superbus, Plectus sambesii, and Propanagrolaimus sp. JU765. For plectids, genomic resources are scarce. We therefore added transcriptome data of Plectus murrayi, Anaplectus granulosus, Neocamacolaimus parasiticus, and Stephanolaimus elegans, with the latter three transcriptomes kindly provided by Dr. Oleksandr Holovachov (Swedish museum of natural history). Transcriptomes for Anaplectus granulosus, and Neocamacolaimus parasiticus have been published and are readily available [53,54]. All three transcriptomes were assembled de novo with Trinity [55]. The exact procedures are described in the respective publications [53,54]. The Stephanolaimus elegans transcriptome was assembled using the same methodologies as Neocamacolaimus parasiticus.
The Plectus murrayi transcriptome was built from raw reads deposited at NCBI (https://sra-downloadb.be-md.ncbi.nlm.nih.gov/sos2/sra-pub-run-13/SRR6827978/SRR6827978.1; accessed 22.12.2020) and assembled using Galaxy Trinity version 2.9.1 [55,56]. All default options were used including in silico normalization of reads before assembly. Transdecoder (conda version 5.5.0) [57] was used to translate to amino acid sequence. Identical reads were removed with cd-hit version 4.8.1 [58,59], with shorter isoforms removed using the Trinity get_longest_isoform_seq_per_trinity_gene.pl command [57] (Trinity conda version 2.8.5; Anaconda Software Distribution, Conda, Version 4.9.2, Anaconda, Nov. 2020). Amino acid translations of the longest isoforms were extracted with AGAT (Dainat, https://www.doi.org/10.5281/zenodo.3552717) from genome assembly FASTA files and genome annotation GFF3 files using the ‘agat_convert_sp_gxf2gxf.pl’, ‘agat_sp_keep_longest_isoform.pl’ and ‘agat_sp_extract_sequences.pl’ scripts, respectively. All FASTA headers were modified to allow for simple species assignment of each sequence in subsequent analysis. Orthology analysis was conducted with OrthoFinder v. 2.5.1 [60,61] using default settings. For genes of interest, we constructed alignments with MAFFT v. 7.475 [62] using the localpair and maxiterate (1000) functions. Spurious sequences and areas that were not well aligned were removed with Trimal v. 1.4.rev22 [63] (procedure stated in S1 Orthology analysis below each phylogeny). We then ran phylogenetic analysis with Iqtree2 v. 2.0.6 [64], with -bb 1000 option, testing the model for each analysis (models eventually used stated in S1 Orthology analysis). PFAM domains were explored using Interproscan v. 5.50–84.0 [65]. The phylogenies were visualized with Dendroscope 3.7.6 [66] and figures were created with Inkscape (https://inkscape.org). Most of our analysis was performed on the HPC RRZK CHEOPS of the Regional Computing Centre (RRZK) of the University of Cologne.
Phylogenomics
Sequences of 18S and 28S genes from 44 taxa across the Propanagrolaimus, Panagrolaimus, Panagrellus and Halicephalobus genera (all listed in S1 Text were aligned (MAFFT L-INS-I v7.475) [62], concatenated [67]and used to infer a species tree using maximum likelihood via (IQTREE) [68] and partitioned by best-fit models of sequence evolution for both [69]. Nodal support was determined using 1000 bootstrap pseudoreplicates. A further 60 genes from 101 taxa (all listed in S1 Text) were used to confirm the taxonomic position using the supermatrix concatenation methods outlined above. Given the limitations of differential gene sampling, we expanded our phylogenomic analyses to include a coalescence approach using 12,295 ML gene trees inferred for orthogroups containing the target animal. Instances of multiple genes per species per group were treated as paralogs/orthologs and analysed using ASTRAL-Pro [70]. Given the number of copies of genes per orthogroup, we explored whether auto or alloploidy was the source of extra genes observed using the gene-tree reconciliation approach implemented in GRAMPA (Gene-tree Reconciliation Algorithm with MUL (Multi labelled)-trees for Polyploid Analysis) [71]. All gene trees rooted at the midpoint and the final ASTRAL-pro species tree were used as inputs, with the most parsimonious result analyzed further.
Desiccation survival assay
C. elegans dauer larvae desiccation assays were performed as described in [25]. P. kolymaensis desiccation assays were performed similarly as described in [25] with mixed population (Mixture of all larval stages and adults) of the nematodes.
Exposure of nematodes to extreme environments
C. elegans dauer larvae or mixed population (Mixture of all larval stages and adults) of P. kolymaensis nematodes were preconditioned and desiccated as described in [28], then transferred to elevated temperature of 34°C, freezing (-80°C) and anoxia. Anoxic environment was generated in a desiccation chamber at 60%RH by flushing the Nitrogen gas into the chamber. The concentration of oxygen inside the chamber was monitored. After each timepoint they were rehydrated with 500 μl of water for 2–3 hours. Rehydrated worms were transferred to NGM agar plates with E. Coli NA22 as food. Survivors were counted after overnight incubation at 15°C. Each experiment was performed on two different days with at least two technical replicates.
Trehalose quantification from nematode lysates
Trehalose measurements were performed as described in previous reports [27].
Radiolabeling, metabolite extraction and 2D-TLC
The above-mentioned procedures were performed according to previous reports [27,28].
Identification of trehalose-6-phosphate from TLC plates
Normalized aqueous fractions from the non-preconditioned and preconditioned samples were separated by high performance thin layer chromatography (HPTLC), using 1-propanol-methanol-ammonia (32%)-water (28:8:7:7 v/v/v/v) as first, dried for 15 min and 1-butanol-acetone-glacial acetic acid–water (35:35:7:23 v/v/v/v) second dimension respectively. Using the trehalose as a standard on both dimensions of the TLC, the regions of interest were scrapped out from the TLCs. The scraped-out silica was extracted with 10 ml of 50% methanol twice. The fractions were combined, dried under vacuum and dissolved in 100 μl of MS mix solution containing 4:2:1 (Isopropanol:Methanol:Chloroform) with 7.5 mM ammonium formate. Mass spectrometric analysis was performed on a Q Exactive instrument (Thermo Fischer Scientific, Bremen, DE) equipped with a robotic nanoflow ion source TriVersa NanoMate (Advion BioSciences, Ithaca, USA) using nanoelectrospray chips with a diameter of 4.1 μm. The ion source was controlled by the Chipsoſt 8.3.1 soſtware (Advion BioSciences). Ionization voltage was + 0.96 kV in negative mode; backpressure was set at 1.25 psi. The temperature of the ion transfer capillary was 200°C; S-lens RF level was set to 50%. FT MS spectra were acquired within the range of m/z 50–750 at the mass resolution of R m/z 200 = 140000; automated gain control (AGC) of 3×106 and with the maximal injection time of 3000 ms. FT MS/MS spectra were acquired within the range of m/z 50–750 at the mass resolution of R m/z 200 = 140000; automated gain control (AGC) of 3×104 and with a maximal injection time of 30 s.
Triacylglycerols measurement from P. kolymaensis lysates
Non-preconditioned and preconditioned pellets were lysed in 200 μl of isopropanol with 0.5 mm Zircornium beads twice for 15 min. The lysates were centrifuged at 1300 g for 5 min. The supernatant was carefully collected without any debris, 20 μl of the lysate was used for protein estimation. Normalization was performed according to soluble protein levels, supernatant volumes corresponding to 50–100 μg of proteins were dried in the desiccator. 700 μl of IS ((10:3 (Methyl tert-butyl ether: ethanol)) mix (warmed to room temperature) was added to dried samples and left on the shaker for 1 hour. The samples were centrifuged at 1400 rpm and 4°C. 140 μl of water was added and left on the shaker for 15 min. These samples were centrifuged at 13400 rpm for 15 min. The upper organic fraction was collected and transferred to 1.5 ml glass vial and left for drying in the desiccator. The dried samples were reconstituted in a volume of 300 μl of 4:2:1 (Isopropanol:Methanol:Chloroform). Volume corresponding to 1 μg was used for injection.
LC-MS/MS analysis was performed on a high-performance liquid chromatography system (Agilent 1200 HPLC) coupled to a Xevo G2-S QTof (Waters). The samples were resolved on a reverse phase C18 column (Cortecs C18 2.7um from Waters) with 50:50:0.1:1% (Water:Methanol:Formicacid:1MAmmoniumformate) and 25:85:0.1:1% (Acetonitrile:Isopropanol:Formic acid:1M Ammonium formate) as mobile phase. The following gradient program was used: Eluent B from 0% to 100% within 12 min; 100% from 12 min to 17min; 0% from 17 min to 25 min. The flow rate was set at 0.3 ml/min. The samples were normalised according to the total protein concentration and the worm numbers. TAG 50:00:00 was used as internal standard. Skylinesoftware (https://skyline.ms/project/home/software/Skyline/begin.view) was used to analyse the raw data. TAGs were extracted from Lipidmaps (https://www.lipidmaps.org/) database.
Supporting information
S1 Fig. Calibration of a radiocarbon (14C) date.
Radiocarbon date (44,315±405 BP) and calibrated age (45,839–47,769 cal BP) of plant material collected from buried borrow P-1320. Radiocarbon ages were converted to calendar age equivalents with the OxCal V.4.4 program using the IntCal20 calibration curve. Pink-shaded area—radiocarbon date with standard deviation; grey-shaded area—radiocarbon date projection on the calibration curve with 95.4% probability.
https://doi.org/10.1371/journal.pgen.1010798.s001
(PDF)
S2 Fig. Morphology of P. kolymaensis female.
Graphic presentations of holotype (A, B) and SEM pictures (C-I): a) anterior body, B) female reproductive branch, C–E) anterior end of three different female specimens, F) anterior part of the lateral ridge, G) vulva, H) ventral excretory/secretory pore, I) posterior body with anus and lateral ridge. Scale bars: a,b—50 μm, c,i—3 μm, d—2 μm, e –1 μm, f,h—5 μm, g—10 μm,.
https://doi.org/10.1371/journal.pgen.1010798.s002
(PDF)
S3 Fig. Phylogenies inferred for genes sets using both concatenation and coalescence approaches.
A) The ploidy level was analysed with Smudgeplot v0.2.173. KMC version 3.1.074 was used to count the 21-mers in the PacBio CCS reads. Then we ran smudgeplot.py to determine the lower and upper coverage cut-offs. These were determined to be 14 and 380. 21-mers with a coverage between 14 and 380 were filtered with kmc_tools. Then, we computed the k-mer pairs from the filtered 21-mers by running smudgeplot.py hetkmers. Finally, the produced a smudgeplot shows an estimated ploidy of 3. B) Species tree inferred for 12,295 gene trees using coalescence approach. The species tree implemented using the coalescence approach with orthogroup gene tree is displayed. Novel species in this study are displayed in red. All nodes have a posterior probability of 1. C) Species tree inferred for 102 taxa. The maximum likelihood tree inferred using a concatenated supermatrix of 60 genes is displayed. Bootstrap values are only displayed for nodes with less than 100% support. The platyhelminth Macrostum lignano serves as an outgroup for rooting. The ancestral Panagrolaimus, sister to all others within the Panagrolaimus genus, is highlighted in red.
https://doi.org/10.1371/journal.pgen.1010798.s003
(PDF)
S4 Fig. Combination of cryptobiotic states enhances survival of C. elegans dauer larvae.
A) Desiccated dauer larvae manifest enhanced survival rate to heat stress (34°C). Error bars indicate standard deviation of two independent experiments with two technical replicates. B) Desiccated dauer larvae display enhances survival rate to anoxia. Error bars indicate standard deviation of two independent experiments with two technical replicates. Statistical comparison was performed by paired two tailed t-test. *p<0.05. C) P. kolymaensis possesses gene orthologs to most genes implicated in dauer formation and metabolism in C. elegans. Black filled circles: Ortholog presence suggested by orthogroup clustering, phylogenetic analysis, and domain architecture. White filled circles: No ortholog found via current analysis (in all cases these C. elegans genes did not cluster with any Panagrolaimus genes in the orthogroup clustering). Label: C. elegans enzyme names and orthogroup that contains that gene according to our orthogroup clustering.
https://doi.org/10.1371/journal.pgen.1010798.s004
(PDF)
S5 Fig. P. kolymaensis reduces triacylglycerols (TAGs) levels and accumulates trehalose-6-phosphate upon preconditioning at 98%RH.
A) 1D-Thin layer chromatography of acetate labelled organic fractions of non-preconditioned (1) and preconditioned (2) P. kolymaensis. B) Mass spectrometric quantification of TAG levels of non-preconditioned (1) and preconditioned (2) P. kolymaensis. Error bars indicate standard deviation of two biological replicates with two technical replicates. Statistical analysis was performed using unpaired t-test with Welch correction **p< 0.001. C-D) non-preconditioned and preconditioned mass spectrum of an empty region, e-f) spot 1 (trehalose), G-H) spot 7 (trehalose-6-phosphate) scraped out and extracted from the 2D-TLC.
https://doi.org/10.1371/journal.pgen.1010798.s005
(PDF)
S6 Fig. Exploring a possible allopolyploid origin for extra proteins.
Gene-tree reconciliation was used to determine whether extra sets of proteins across orthogroups originate through auto- or allopolyploidy. Different copies of proteins (designated by ‘+’ and ‘*’) suggest an allopolyploid origin. Parthenogenetic species are highlighted in bold.
https://doi.org/10.1371/journal.pgen.1010798.s006
(PDF)
S1 Text. Supplementary Information and methods.
https://doi.org/10.1371/journal.pgen.1010798.s007
(PDF)
S1 Data. Desiccated dauer larvae manifest enhanced survival rate to heat stress (34°C).
https://doi.org/10.1371/journal.pgen.1010798.s008
(XLSX)
S2 Data. Desiccated dauer larvae display enhances survival rate to anoxia.
https://doi.org/10.1371/journal.pgen.1010798.s009
(XLSX)
S3 Data. Mass spectrometric quantification of TAG levels of non- preconditioned (1) and preconditioned (2) P. kolymaensis.
https://doi.org/10.1371/journal.pgen.1010798.s010
(XLSX)
S4 Data. Survival rate of P. kolymaensis nematodes to desiccation and freezing (-80°C).
https://doi.org/10.1371/journal.pgen.1010798.s011
(XLSX)
S5 Data. P. kolymaensis nematodes and daf-2(e1370) dauer larvae upregulate trehalose levels upon preconditioning at 98%RH.
https://doi.org/10.1371/journal.pgen.1010798.s012
(XLSX)
S6 Data. Desiccated daf-2 (e1370) dauer larvae survive to freezing (-80°C) for an extremely long period.
https://doi.org/10.1371/journal.pgen.1010798.s013
(XLSX)
S7 Data. Brood size of desiccated dauer larvae exposed to freezing remain like that non- desiccated dauer larvae.
https://doi.org/10.1371/journal.pgen.1010798.s014
(XLSX)
S1 Orthology analysis. Supplementary file orthology analysis.
https://doi.org/10.1371/journal.pgen.1010798.s015
(PDF)
Acknowledgments
VG is thankful to Andrej Shevchenko, members of the Volkswagen grant and Kurzchalia lab for helpful discussions and the core facilities of MPI-CBG for assistance. We are grateful to Dr. S. Gubin for field study and sampling, our colleagues in Soil Cryology Lab, Pushchino and North-East Scientific Station in Chersky, Republic of Sakha (Yakutia) for their help and cooperation. The authors thank Long Read Team of the DRESDEN-concept Genome Center, DFG NGS Competence Center, part of the Center for Molecular and Cellular Bioengineering (CMCB), Technische Universität Dresden and MPI-CBG. We thank https://www.copernicus.eu/en for the base map in Fig 1A. The authors are thankful to Richard Roy and Jens Bast for critically reading the manuscript. We thank Iain Pattern for suggestions on writing the manuscript.
References
- 1. Keilin D. The Leeuwenhoek Lecture—The problem of anabiosis or latent life: history and current concept. Proc R Soc Lond Ser B—Biol Sci. 1959;150: 149–191.
- 2. Cano R, Borucki M. Revival and identification of bacterial spores in 25- to 40-million-year-old Dominican amber. Science. 1995;268: 1060–1064. pmid:7538699
- 3. Shen-Miller J. Sacred lotus, the long-living fruits of China Antique. Seed Sci Res. 2002;12: 131–143.
- 4. Guidetti R, Jönsson KI. Long-term anhydrobiotic survival in semi-terrestrial micrometazoans. J Zool. 2002;257: 181–187.
- 5. Tsujimoto M, Imura S, Kanda H. Recovery and reproduction of an Antarctic tardigrade retrieved from a moss sample frozen for over 30 years. Cryobiology. 2016;72: 78–81. pmid:26724522
- 6. Kagoshima H, Kito K, Aizu T, Shin-i T, Kanda H, Kobayashi S, et al. Multi-decadal survival of an Antarctic nematode, Plectus murrayi, in a -20°C stored moss sample. CryoLetters. 2012;33: 280–288.
- 7. Steiner G, Florence EA. Resuscitation of the nematode Tylenchus polyhypnus, n. sp., after almost 39 years’ dormancy. J Wash Acad Sci. 36: 97–99.
- 8. Gilichinsky DA, Wilson GS, Friedmann EI, Mckay CP, Sletten RS, Rivkina EM, et al. Microbial Populations in Antarctic Permafrost: Biodiversity, State, Age, and Implication for Astrobiology. Astrobiology. 2007;7: 275–311. pmid:17480161
- 9. Rivkina E, Abramov A, Spirina E, Petrovskaya L, Shatilovich A, Shmakova L, et al. Earth’s perennially frozen environments as a model of cryogenic planet ecosystems. Permafr Periglac Process. 2018;29: 246–256.
- 10. Sipes K, Almatari A, Eddie A, Williams D, Spirina E, Rivkina E, et al. Eight Metagenome-Assembled Genomes Provide Evidence for Microbial Adaptation in 20,000- to 1,000,000-Year-Old Siberian Permafrost. Kelly RM, editor. Appl Environ Microbiol. 2021;87. pmid:34288700
- 11. Liang R, Li Z, Lau Vetter MCY, Vishnivetskaya TA, Zanina OG, Lloyd KG, et al. Genomic reconstruction of fossil and living microorganisms in ancient Siberian permafrost. Microbiome. 2021;9: 110. pmid:34001281
- 12. Gilichinsky DA, Rivkina EM. Permafrost Microbiology. Encyclopedia of Geobiology. 2011:726–732.
- 13. Yashina S, Gubin S, Maksimovich S, Yashina A, Gakhova E, Gilichinsky D. Regeneration of whole fertile plants from 30,000-y-old fruit tissue buried in Siberian permafrost. Proc Natl Acad Sci. 2012;109: 4008–4013. pmid:22355102
- 14. Shmakova L, Malavin S, Iakovenko N, Vishnivetskaya T, Shain D, Plewka M, et al. A living bdelloid rotifer from 24,000-year-old Arctic permafrost. Curr Biol. 2021;31: R712–R713. pmid:34102116
- 15. Malavin S, Shmakova L, Claverie J-M, Rivkina E. Frozen Zoo: a collection of permafrost samples containing viable protists and their viruses. Biodivers Data J. 2020;8: e51586. pmid:32733138
- 16. Krivushin K, Kondrashov F, Shmakova L, Tutukina M, Petrovskaya L, Rivkina E. Two Metagenomes from Late Pleistocene Northeast Siberian Permafrost. Genome Announc. 2015;3: e01380–14. pmid:25555741
- 17. Zhou J, Davey ME, Figueras JB, Rivkina E, Gilichinsky D, Tiedje JM. Phylogenetic diversity of a bacterial community determined from Siberian tundra soil DNA. Microbiology. 1997;143: 3913–3919. pmid:9421915
- 18. Kapterev P. New data on revitalization of organisms from perpetually frozen grounds. CR Acad Sci URSS. 1938;20: 315–317.
- 19. Kapterev PN. Anabiosis in the conditions of permanent congelation. Bulletin de L’Academie des Sciences de L’URSS. 1936;6: 1073–1088.
- 20. Adhikari BN, Wall DH, Adams BJ. Desiccation survival in an Antarctic nematode: molecular analysis using expressed sequenced tags. BMC Genomics. 2009;10: 69. pmid:19203352
- 21. Tyson T, O’Mahony Zamora G, Wong S, Skelton M, Daly B, Jones JT, et al. A molecular analysis of desiccation tolerance mechanisms in the anhydrobiotic nematode Panagrolaimus superbus using expressed sequenced tags. BMC Res Notes. 2012;5: 68. pmid:22281184
- 22. McGill LM, Shannon AJ, Pisani D, Félix M-A, Ramløv H, Dix I, et al. Anhydrobiosis and Freezing-Tolerance: Adaptations That Facilitate the Establishment of Panagrolaimus Nematodes in Polar Habitats. Liang W, editor. PLOS ONE. 2015;10: e0116084. pmid:25747673
- 23. Thorne MAS, Kagoshima H, Clark MS, Marshall CJ, Wharton DA. Molecular Analysis of the Cold Tolerant Antarctic Nematode, Panagrolaimus davidi. Winter CE, editor. PLoS ONE. 2014;9: e104526. pmid:25098249
- 24. Sandhove J, Spann N, Ristau K. The Anhydrobiotic Potential of the Terrestrial Nematodes Plectus parietinus and Plectus velox: Anhydrobiosis of Plectus Parietinus and P. Velox. J Exp Zool Part Ecol Genet Physiol. 2016;325: 434–440.
- 25. Erkut C, Penkov S, Khesbak H, Vorkel D, Verbavatz J-M, Fahmy K, et al. Trehalose Renders the Dauer Larva of Caenorhabditis elegans Resistant to Extreme Desiccation. Curr Biol. 2011;21: 1331–1336. pmid:21782434
- 26. Erkut C, Vasilj A, Boland S, Habermann B, Shevchenko A, Kurzchalia TV. Molecular Strategies of the Caenorhabditis elegans Dauer Larva to Survive Extreme Desiccation. Braendle C, editor. PLoS ONE. 2013;8: e82473. pmid:24324795
- 27. Gade VR, Traikov S, Oertel J, Fahmy K, Kurzchalia TV. C. elegans possess a general program to enter cryptobiosis that allows dauer larvae to survive different kinds of abiotic stress. Sci Rep. 2020;10: 13466.
- 28. Erkut C, Gade VR, Laxman S, Kurzchalia TV. The glyoxylate shunt is essential for desiccation tolerance in C. elegans and budding yeast. eLife. 2016;5: e13614. pmid:27090086
- 29. Zanina OG, Gubin SV, Kuzmina SA, Maximovich SV, Lopatina DA. Late-Pleistocene (MIS 3–2) palaeoenvironments as recorded by sediments, palaeosols, and ground-squirrel nests at Duvanny Yar, Kolyma lowland, northeast Siberia. Quat Sci Rev. 2011;30: 2107–2123.
- 30. Murton JB, Goslar T, Edwards ME, Bateman MD, Danilov PP, Savvinov GN, et al. Palaeoenvironmental Interpretation of Yedoma Silt (Ice Complex) Deposition as Cold-Climate Loess, Duvanny Yar, Northeast Siberia: Palaeoenvironmental Interpretation of Yedoma Silt, Duvanny Yar. Permafr Periglac Process. 2015;26: 208–288.
- 31. Abramov A, Vishnivetskaya T, Rivkina E. Are permafrost microorganisms as old as permafrost? FEMS Microbiol Ecol. 2021;97: fiaa260. pmid:33601419
- 32. Gilichinsky DA, Wagener S, Vishnevetskaya TA. Permafrost microbiology. Permafr Periglac Process. 1995;6: 281–291.
- 33. Shatilovich AV, Tchesunov AV, Neretina TV, Grabarnik IP, Gubin SV, Vishnivetskaya TA, et al. Viable Nematodes from Late Pleistocene Permafrost of the Kolyma River Lowland. Dokl Biol Sci. 2018;480: 100–102. pmid:30009350
- 34. Schiffer PH, Danchin EGJ, Burnell AM, Creevey CJ, Wong S, Dix I, et al. Signatures of the Evolution of Parthenogenesis and Cryptobiosis in the Genomes of Panagrolaimid Nematodes. iScience. 2019;21: 587–602. pmid:31759330
- 35. Schiffer PH, Nsah NA, Grotehusmann H, Kroiher M, Loer C, Schierenberg E. Developmental variations among Panagrolaimid nematodes indicate developmental system drift within a small taxonomic unit. Dev Genes Evol. 2014;224: 183–188. pmid:24849338
- 36. Dey A, Chan CKW, Thomas CG, Cutter AD. Molecular hyperdiversity defines populations of the nematode Caenorhabditis brenneri. Proc Natl Acad Sci. 2013;110: 11056–11060. pmid:23776215
- 37. Cutter AD. Divergence Times in Caenorhabditis and Drosophila Inferred from Direct Estimates of the Neutral Mutation Rate. Mol Biol Evol. 2008;25: 778–786. pmid:18234705
- 38. Seybold AC, Wharton DA, Thorne MAS, Marshall CJ. Investigating trehalose synthesis genes after cold acclimation in the Antarctic nematode Panagrolaimus sp. DAW1. Biol Open. 2017; bio.023341.
- 39. Evangelista CCS, Guidelli GV, Borges G, Araujo TF, de Souza TAJ, da C Neves UP, et al. Multiple genes contribute to anhydrobiosis (tolerance to extreme desiccation) in the nematode Panagrolaimus superbus. Genet Mol Biol. 2017;40: 790–802. pmid:29111563
- 40.
Andrássy I. Klasse Nematoda. Gustav Fischer Verlag; 1984.
- 41. Lewis SC, Dyal LA, Hilburn CF, Weitz S, Liau W-S, LaMunyon CW, et al. Molecular evolution in Panagrolaimus nematodes: origins of parthenogenesis, hermaphroditism and the Antarctic species P. davidi. BMC Evol Biol. 2009;9: 15. pmid:19149894
- 42.
Andrássy I. Free-living Nematodes of Hungary. Hungarian Natural History Museum; 2005.
- 43. Schulenburg H, Félix M-A. The Natural Biotic Environment of Caenorhabditis elegans. Genetics. 2017;206: 55–86. pmid:28476862
- 44. Frézal L, Félix M-A. C. elegans outside the Petri dish. eLife. 2015;4: e05849. pmid:25822066
- 45. Gems D, Sutton AJ, Sundermeyer ML, Albert PS, King KV, Edgley ML, et al. Two pleiotropic classes of daf-2 mutation affect larval arrest, adult behavior, reproduction and longevity in Caenorhabditis elegans. Genetics. 1998;150: 129–155.
- 46. Vowels JJ, Thomas JH. Multiple chemosensory defects in daf-11 and daf-21 mutants of Caenorhabditis elegans. Genetics. 1994;138: 303–316. pmid:7828815
- 47. Ren P, Lim C-S, Johnsen R, Albert PS, Pilgrim D, Riddle DL. Control of C. elegans Larval Development by Neuronal Expression of a TGF-β Homolog. Science. 1996;274: 1389–1391. pmid:8910282
- 48. Gerisch B, Weitzel C, Kober-Eisermann C, Rottiers V, Antebi A. A Hormonal Signaling Pathway Influencing C. elegans Metabolism, Reproductive Development, and Life Span. Dev Cell. 2001;1: 841–851. pmid:11740945
- 49. Penkov S, Kaptan D, Erkut C, Sarov M, Mende F, Kurzchalia TV. Integration of carbohydrate metabolism and redox state controls dauer larva formation in Caenorhabditis elegans. Nat Commun. 2015;6: 8060. pmid:26290173
- 50. Penkov S, Raghuraman BK, Erkut C, Oertel J, Galli R, Ackerman EJM, et al. A metabolic switch regulates the transition between growth and diapause in C. elegans. BMC Biol. 2020;18: 31. pmid:32188449
- 51. Nurk S, Walenz BP, Rhie A, Vollger MR, Logsdon GA, Grothe R, et al. HiCanu: accurate assembly of segmental duplications, satellites, and allelic variants from high-fidelity long reads. Genome Res. 2020;30: 1291–1305. pmid:32801147
- 52. Laetsch DR, Blaxter ML. BlobTools: Interrogation of genome assemblies. F1000Research. 2017;6: 1287.
- 53. Smythe AB, Holovachov O, Kocot KM. Improved phylogenomic sampling of free-living nematodes enhances resolution of higher-level nematode phylogeny. BMC Evol Biol. 2019;19: 121. pmid:31195978
- 54. Ahmed M, Adedidran F, Holovachov O. A draft transcriptome of a parasite Neocamacolaimus parasiticus (Camacolaimidae, Plectida). J Nematol. 2021;53: 1–4. pmid:33860270
- 55. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29: 644–652. pmid:21572440
- 56. Goecks J, Nekrutenko A, Taylor J, Galaxy Team T. Galaxy: a comprehensive approach for supporting accessible, reproducible, and transparent computational research in the life sciences. Genome Biol. 2010;11: R86. pmid:20738864
- 57. Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 2013;8: 1494–1512. pmid:23845962
- 58. Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22: 1658–1659. pmid:16731699
- 59. Fu L, Niu B, Zhu Z, Wu S, Li W. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics. 2012;28: 3150–3152. pmid:23060610
- 60. Emms DM, Kelly S. OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. Genome Biol. 2015;16: 157. pmid:26243257
- 61. Emms DM, Kelly S. OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 2019;20: 238. pmid:31727128
- 62. Katoh K, Standley DM. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability. Mol Biol Evol. 2013;30: 772–780. pmid:23329690
- 63. Capella-Gutierrez S, Silla-Martinez JM, Gabaldon T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009;25: 1972–1973. pmid:19505945
- 64. Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, et al. IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. Teeling E, editor. Mol Biol Evol. 2020;37: 1530–1534. pmid:32011700
- 65. Jones P, Binns D, Chang H-Y, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30: 1236–1240. pmid:24451626
- 66. Huson DH, Scornavacca C. Dendroscope 3: An Interactive Tool for Rooted Phylogenetic Trees and Networks. Syst Biol. 2012;61: 1061–1067. pmid:22780991
- 67. Kück P, Meusemann K. FASconCAT: Convenient handling of data matrices. Mol Phylogenet Evol. 2010;56: 1115–1118. pmid:20416383
- 68. Nguyen L-T, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies. Mol Biol Evol. 2015;32: 268–274. pmid:25371430
- 69. Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14: 587–589. pmid:28481363
- 70. Zhang C, Scornavacca C, Molloy EK, Mirarab S. ASTRAL-Pro: Quartet-Based Species-Tree Inference despite Paralogy. Thorne J, editor. Mol Biol Evol. 2020;37: 3292–3307. pmid:32886770
- 71. Thomas GWC, Ather SH, Hahn MW. Gene-Tree Reconciliation with MUL-Trees to Resolve Polyploidy Events. Renner S, editor. Syst Biol. 2017;66: 1007–1018. pmid:28419377