- Methodology article
- Open Access
Environmental DNA sequencing primers for eutardigrades and bdelloid rotifers
BMC Ecologyvolume 9, Article number: 25 (2009)
The time it takes to isolate individuals from environmental samples and then extract DNA from each individual is one of the problems with generating molecular data from meiofauna such as eutardigrades and bdelloid rotifers. The lack of consistent morphological information and the extreme abundance of these classes makes morphological identification of rare, or even common cryptic taxa a large and unwieldy task. This limits the ability to perform large-scale surveys of the diversity of these organisms.
Here we demonstrate a culture-independent molecular survey approach that enables the generation of large amounts of eutardigrade and bdelloid rotifer sequence data directly from soil. Our PCR primers, specific to the 18s small-subunit rRNA gene, were developed for both eutardigrades and bdelloid rotifers.
The developed primers successfully amplified DNA of their target organism from various soil DNA extracts. This was confirmed by both the BLAST similarity searches and phylogenetic analyses. Tardigrades showed much better phylogenetic resolution than bdelloids. Both groups of organisms exhibited varying levels of endemism.
The development of clade-specific primers for characterizing eutardigrades and bdelloid rotifers from environmental samples should greatly increase our ability to characterize the composition of these taxa in environmental samples. Environmental sequencing as shown here differs from other molecular survey methods in that there is no need to pre-isolate the organisms of interest from soil in order to amplify their DNA. The DNA sequences obtained from methods that do not require culturing can be identified post-hoc and placed phylogenetically as additional closely related sequences are obtained from morphologically identified conspecifics. Our non-cultured environmental sequence based approach will be able to provide a rapid and large-scale screening of the presence, absence and diversity of Bdelloidea and Eutardigrada in a variety of soils.
Micro-invertebrates, though very important to the soil biocenose (self-regulating ecological communities) and energy flux of a system, are still poorly understood in terms of their taxonomy and geographical distributions [1–4]. Like many microfaunal organisms, Rotifera and Tardigrada pose problems for taxonomists and evolutionary biologists due to the difficulties associated with isolation, identification and enumeration of organisms that do not preserve any discernable morphological characters. Even when it is possible to successfully culture these organisms, limited phenotypic differentiation among taxa and cyclomorphosis (seasonal change in body shape; ) confound accurate taxonomy. This lack of consistent morphological information and the extreme abundance of meiofaunal organisms makes identification of rare, or even common, cryptic taxa a large and unwieldy task [6, 4] as only painstaking microscopy can be used to identify synapomorphies.
Environmental sequencing is valuable for performing large-scale surveys of the diversity of organisms that cannot be cultured or grown in the laboratory or when species are difficult to distinguish using phenotypic characters. These issues argue for culture independent molecular surveys of meiofaunal diversity in natural ecosystems. Microbiologists have faced many of the same problems and solved them by turning to conserved DNA sequences as a means of describing communities [7, 8]. Instead of isolating and culturing individuals, communities are characterized by extracting all of the DNA in a particular sample (soil, water, air), amplifying a specific gene using PCR, cloning individual PCR products, and then sequencing individual clones. This environmental DNA approach has revolutionized microbiology. For example, these techniques have been successfully used to provide new insights into fungi [9, 10], novel Chloroflexi , abundance and distribution of Psychrobacter and Exiguobacterium  and have been used to provide information about the structure and function of alpine and arctic soil microbial communities .
Our survey focuses on the 18S rRNA gene, commonly used for phylogenetic inference of eukaryotes due to its highly conserved sequence and ability to resolve relatively deep nodes. This is the first description of the general utility of environmental DNA sequencing approaches for characterizing difficult to study ecological communities of eutardigrades and bdelloid rotifers.
Environmental sequencing as described here differs from other molecular survey methods [6, 14] in that there is no need to pre-isolate the bdelloid rotifers or eutardigrades of interest from soil (or other mediums) before amplifying their DNA. The successful development of clade-specific 18s SSU primers has shown to be effective when surveying the diversity of targeted groups of organisms. For example, clade specific 18s SSU primers have been used to describe soil metazoans[15, 16] and reveal the hidden diversity and biogeographic endemism of kinetoplastids (flagellate protozoa) .
The use of 18S rDNA allows for sequences to be combined into already existing 18S and 16S rDNA databases, including those being developed by microbial ecologists through their large scale molecular surveys as referred to above. Here we describe the utility of screening for bdelloid rotifer and eutardigrade diversity in two very distinct sample sites with targeted 18S primers: the high-elevation sites located within the Niwot Ridge Long Term Ecological Research (LTER) site in the Colorado Rockies, and the low-elevation sites located within the Calhoun Experimental Forest in South Carolina.
We developed two forward primers for taxon specific amplification of eutardigrades and bdelloid rotifers. These primers were used in combination with a universal reverse 18S rDNA primer to specifically characterize the diversity of these two groups from several environments. PCR, BLAST and phylogenetic analysis confirmed that each set of primers amplifies the targeted groups with fidelity and specificity (Figures 1 &2). We have observed many invertebrates within the soils prior to DNA extraction and amplification, including mites, nematodes, and insects; none of these were observed within the sequencing data produced using the specific primers in this study. Thus, our primers are shown to be specific to the targeted groups of organisms. The closest known sequences or clades to the environmental sequences are noted below. Note that we do not infer that the environmental sequences are of the same species or genera to those closest to them.
Out of 1,814 nucleotide positions there were 900 variable sites, of which 677 were phylogenetically informative, comprising 68 unique phylotypes. Phylogenetic analysis clearly separates the two main groups of tardigrades: the Heterotardigrada and the Eutardigrada (Figure 3). Many of the environmental sequences from the high-elevation talus sites clustered into distinct clades, suggesting each clade may comprise a separate species. Eutardigrade sequences from soils near the Arikiree Glacier (AGL) grouped within the Macrobiotoidea and Hypsibiodea groups. Those within the Macrobiotidea are most closely related to Richtersius coronifer, a cosmopolitan species sampled from high elevation and arctic habitats . The AGL sequences that grouped within the Hypsibiodea are related to those of the englacial dominating Hypsibius genus. These Hypsibius sequences from the AGL site are nearby and grouped with the two talus sites (T1T2 & T3T6).
The Calhoun Hardwood site sequences cluster closest to Isohypsibius papillifer typically found in Europe, Asia, Australia, & South America . The genus Isohypsibius is composed of species that are widespread and has been documented circumglobally as well [19, 20], (GBIF Sweden, 17 records; National Museum of Natural History, 10 records; Australian Antarctic Data Centre, 3 records).
The Calhoun Grassland sequences cluster basally with the Arikiree and Talus sites within the Hypsibius group, noted as "Acutuncus/Hypsibius" in contrast to another group labeled "Acutuncus/Calohypsibius" in Figure 3, (see [21, 22] for clarification about taxonomic identification issues with Hypsibius and Acutuncus).
Out of 1638 sites 896 were variable and 718 where phylogenetically informative. The environmentally obtained sequences totaled 54 unique phylotypes (49 from this study). Phylogenetic analysis clearly separates all of the main clades of rotifers: Seisonidea, Monogononta and Bdelloidea (Figure 4). All of the environmental sequences we sampled grouped within the Bdelloidea. We also discovered three relatively diverse clades. The first is dominated by Niwot Ridge sequences (Clade A). One of the clades within Clade A (Sub A) is mainly dominated by sequence types from the T1T2 site. The second clade (Clade B) is dominated by those sequences from the Calhoun sites. What is interesting here is that the most derived cluster within Clade B contains uncultured sequences from Japan (Ibaraki upland soils) along with sequences from a high elevation site in Socompa, South America (Sub B) .
The final main group of sequences, Clade C, contains sequences from several locales, but mostly those from the AGL site. Again, like in Sub B, we observe uncultured sequence data from Japan (Fukushima and Ibaraki) clustering with a sequence from Socompa.
The lack of 18S rDNA sequence information in online data bases (14 bdelloid sequences in GenBank  as of this writing), makes the identification of environmentally obtained sequences even more difficult.
The development of clade-specific primers that allows characterization of eutardigrade and bdelloid rotifer communities from environmental samples should greatly increase our ability to discern the community diversity of these taxa in environmental samples. Moreover, the rDNA sequence data can be directly stored (within software packages like ARB ) and compared with other surveys that attempt to characterize invertebrate community composition [16, 26].
We anticipate that environmental DNA surveys using clade-specific primers, like those we have developed, will be used to complement more directed studies that cultivate individual micro-eukaryotes as a means of more fully describing the diversity of ecological communities. We have yet to assess whether isolation of individuals and environmental DNA surveys yield different estimates of community composition, as is the case for surveys of bacteria (but see [16, 26]) and bdelloid rotifers .
Environmental sequencing as shown here differs from other molecular survey methods [6, 14, 16, 26] in that there is no need to pre-isolate the organisms of interest from soil (or other media), in order to amplify their DNA. Here, we simply extract total cellular DNA from all organisms in the soil and use targeted primers for the group of interest. This allows for a single DNA extraction prep instead of one DNA extraction prep for each targeted organism of interest.
Although there are too few data to make robust biological inferences, several results are noteworthy. We found sequences from the highest elevation site in Colorado (near the Arikiree glacier) that grouped together with R. coronifer, a cosmopolitan morpho-species known to exist in high mountain and arctic habitats which is also known to survive extreme desiccation and temperatures down to -196°C . Additionally, several sequences from the Calhoun hardwood forest were very similar to Isohypsibius papillifer, a widespread European species. Moreover, the genus Isohypsibius is ubiquitous, distributed from North America, Northern Europe, and Asia, all the way to Antarctica.
Interestingly, the sequences from the AGL site seem to have the most distant set of sequences compared to the other sites. One set of sequences is from within the Macrobiotidae, Richtersius group and the other from the Hypsibiodea, Hypsibius group. This is probably due to the longer duration of moist and wet soils that allows for a greater diversity of eutardigrade groups.
It is not too suprising that the majority of the Eutardigrade sequences amplified from the Talus and glacier sites are dominated by Hypsibius-related sequences. The Hypsibidae are known to dominate englacial habitats and are the dominate family of polar and cryoconite tardigrades. Hypsibius species are hydrophilic and are composed of bacteriophagous and/or algivorous feeding types. These biological factors aid in the colonization of nunatuks and glacial habitats (as reviewed in ).
However, several sequences from the Macrobiotidae were also found within the glacial habitat of the AGL site. Macrobiotidae are traditionally considered cosmopolitan occurring in many habitats, including those that are periodically frozen . The AGL sequences cluster closest to the known sequences of Richtersius sp. (Figure 3). Richtersius have been the focus of many anhydrobiosis studies and have shown significant improvements in desiccation survival when many individuals aggregate together during anhydrobiosis . This could lead to positive density dependence and even allow these animals to achieve greater monopolization [as reviewed in ] to local habitats that encounter extreme desiccation events like the high elevation AGL and talus sites. However, aggregation can create problems with environmental sequencing strategies like the one proposed here. If aggregation in the wild occurs within other eutardigrade groups then environmental sequencing may lead to amplification of only those extremely high-abundant clusters of animals.
In contrast to the tardigrades, there was less agreement of support between the two different phylogenetic reconstruction methods of Bayesian and parsimony analysis for bdelloid rotifers. It was not possible to identify what bdelloids the environmental sequences were related to due to lack of abundant reference sequences. However, while it was possible to make some general statements about the bdelloid communities at the listed sample sites, the lack of resolution of 18S rDNA compared to 28S rDNA  makes it difficult to delineate the more recent clades of Bdelloidea (Figure 4). In fact, a similar level of poor resolution of bdelloids is also seen from phylogenies produced via cytochrome oxidase subunit 1 sequence data, wherein the early nodes are mostly saturated with polytomies (Robeson & Birky unpublished). Better resolution of this group at the tips of the phylogeny is often seen regardless of the phylogenetic reconstruction method chosen.
It is interesting that sequences from Socompa  cluster with the Calhoun sequences as opposed to other high elevation sites like the dry Talus, in Niwot Ridge. Although Socompa is a very high elevation site (5824 m above sea level), it is most likely similar in its microhabitat to the Calhoun sites, where there is greater moisture compared to the dry Talus. The Socompa site is characterized as a fumerole environment . Typically fumaroles are areas where steam and volcanic gases vent out of the earth's crust due to the degassing of magma and/or geothermal heating of shallow ground water. This particular fumerole site is weakly active, creating an environment in which communities of mosses and liverworts are sustained by warm water vapor. The potentially similar microhabitats may be the reason for finding such similar sequence types in very different locales.
Bdelloid rotifers in particular show evidence for geographic structure among clades. Whether this apparent pattern reflects environmental filtering, priority effects (differences in arrival time that can have a lasting effect on differences in species dominance), or some other process remains to be seen. Nonetheless, the data presented here support the contention of , in which instances of endemism are seen (Clade A & B), with a few phylogenetic clusters of widespread bdelloids sampled from very different locales (Clade C and Sub B). It may be that harsher conditions in which there are very ephemeral moments of soil moisture creates higher levels of endemism of bdelloids, whereas environments in which soil moisture is sustained for longer periods of time allow for increased chances of long distance dispersal to suitable habitats and persistence. The location of the Socompa fumerole sites in the phylogeny (Figure 4) and its high similarity to sequences from Japan and within the Calhoun sites (Clade B & C) may be an indication of the latter point. One caveat here is that the 18s rDNA sequences are more conserved than their cytochrome oxidase subunit 1 counterparts [4, 33] preserving more ancient than contemporary relatedness.
Large-scale surveys of rotifer and tardigrade diversity using traditional approaches makes for a large and unwieldy set of tasks (i.e. difficulties associated with isolation, identification and enumeration of organisms that do not preserve any discernable morphological characters).
Environmental sequencing is valuable for performing large-scale surveys of the diversity of organisms that cannot be cultured or grown in the laboratory or in which species are difficult to distinguish using phenotypic characters. The DNA sequences obtained from non-cultured based methods can be identified post-hoc (placed phylogentically) as closely related sequences are obtained from morphologically identified conspecifics. Our environmental sequence based approach, which does not require culturing or isolation of animals from soils, provides a rapid and large-scale screening for the presence, absence and diversity of Bdelloidea and Eutardigrada in a variety of soils.
We have shown that targeted amplification of eutardigrades and bdelloid rotifers are possible from a range of soil types. This sequence data can be used to quickly assess the peculiar biogeography [31, 34] and genetic diversity of soil samples, more often informing us of dominate groups within each sample.
It should also be emphasized that environmental sequencing strategies like this are not intended to replace, but instead complement ongoing morphological work, explore the possible effects of heterogeneity within individuals, and the effect of this variation on phylogenetic analysis . This highlights the need for morphological taxonomists and molecular ecologists to work together in order to make environmental sequencing methods, like the one proposed here, more robust. In particular, studies such as these are most empowered by the cataloging of sequence data from vouchered specimens.
Soil DNA extraction
Soil samples (~5 g) were taken from all sites. Three sites from within the Niwot Ridge Long Term Ecological Research (LTER) area in the Front Range of the Colorado Rocky Mountains, United States of America (40° 03' N, 105° 35' W). These sites are: the Arikiree Glacier (AGL), Talus site 1 (T1T2), and Talus site 2 (T3T6) as described previously by . Other soil samples were also obtained from the Calhoun Experimental Forest (managed by the US Department of Agriculture located in northwestern South Carolina in the Piedmont region, 34.5°N, 82°W), these sites are: Hardwood (H), Grassland (G), and Cultivated (C). Total cellular DNA was extracted from soil using the PowerSoil DNA Isolation Kit #12888 (Mo Bio Laboratories, Inc, Carlsbad, CA).
Only forward 18S SSU primers were developed to target specific groups (bdelloids and eutardigrades). Primer development entailed downloading all available target sequences of interest along with their closest set of outgroup taxa from GenBank  and aligned using Muscle  and edited in ARB  to align conserved regions only. A region of bases unique to the target group that excluded as many matches as possible to the outgroup taxa were chosen for primer development. Bdel_2: 5'-CGG CTC ATT ACA TCA GCT ATA ACT T-3' was used for bdelloid rotifers, and Tard_1: 5'-TCT CAG TAC TTG CTT TAA CAA GGC-3' was used for eutardigrades. Amplicon products produced were ~1700 base pairs in length. All eutardigrade and bdelloid rotifer environmental sequences had a sequence identity to those in GenBank ranging from 91 to 98% with a query coverage of 99 to 100% and 95-99% with a query coverage of 97-100% respectively.
Other 'universal' primers used in this study were taken or derived from [38–40] and are listed here as follows: 18S2a: 5'-GAT CCT TCC GCA GGT TCA CC-3'; 18S3: 5'-GAC TCA ACA CGG GAA ACC TCA CC-3'; 18S10: 5'-CTA AGG GCA TCA CAG ACC-3'
The reverse primer 18S2a was used in conjunction with either the Tard_1 or Bdel_2 primer in order to amplify the DNA of either eutardigrades or bdelloid rotifers directly from soil. The PCR cycling conditions were as follows: initial denaturation at 94°C for 2 min, followed by 40 cycles of: 94°C for 30", 60°C for 30", 72°C for 2', with a final extension at 72°C for 10'. PCR reaction contained (all reagents from Invitrogen, Carlsbad, CA, USA) 1× PCR Buffer, 1.5 mM MgCl2, 0.2 μM dNTPs, 0.4 μM of each primer, Taq polymerase (0.5 units), template DNA: 2 μL.
Cloning & Sequencing
The final PCR product was purified using the Wizard SV Gel and PCR Clean-up System (Promega, Madison, WI) or the QIAquick Gel Extraction Kit 28704 (QIAGEN, Valencia, CA). Purified PCR product was then cloned using the Invitrogen TOPO TA Kit (with pCR2.1-TOPO vector) with One Shot TOP10 Chemically Competent E. coli (K4500-01). Pelleted cells were sent to Functional Biosciences, Inc (Madison, WI) for sequencing. The 18S3 and 18S10 primers were only used at this step for internal sequencing along with M13 primers to generate robust sequence data for contig assembly.
Sequence data was assembled, vector and primer sequence removed, then edited by hand using Sequencher 4.7 (Gene Codes Cooporation, Ann Arbor MI). Sequences where chimera-checked using the Bellerophon server  and determined that no chimeras by sample site amplicons were detected. Usable data were then exported for BLAST  searches. All sequences produced and/or used in this study are listed by accession in Table 1.
Pre-aligned guide and outgroup sequences were downloaded from the SILVA database . The SILVA aligner was used to align the environmental 18s rDNA SSU sequence data according to secondary structure . The data was further edited by eye and exported from ARB  using an 'in-house' filter to remove highly ambiguous regions of the alignment. All terminal gaps in the alignment were converted to missing (i.e. as '?' characters) and gaps '-' counted as a 5th character state. TNT  and a multi-core version of MrBayes  were used to confirm the phylogenetic placement of environmentally obtained sequences. Parsimony analysis was performed by generating 1000 bootstrap replicates. Before re-sampling, the trees were collapsed using TBR. Each bootstrap replicate was composed of twenty iterations of 'Wagner addition trees' (trees formed by sequentially adding the taxa at the best available position, using Fitch parsimony) followed by swapping with TBR, the single best tree was then used for random sector searches and trees saved. MrBayes was used to perform 5 and 8 million generations using the GTR + G + I model of evolution as specified by MultiPhyl Online on the bdelloid and eutardigrade data sets respectively .
Kutikova LA: Bdelloid rotifers (Rotifera, Bdelloidea) as a component of soil and land biocenoses. Biol Bull. 2003, 30 (3): 271-274. 10.1023/A:1023811929889. Translated from Izvestiya Akademii Nauk, Seriya Biologicheskaya 2003, 3:332-336
Stead TK, Schmid-Araya JM, Schmid PE, Hildrew AG: The distribution of body size in a stream community: one system, many patterns. J Anim Eco. 2005, 74: 475-487. 10.1111/j.1365-2656.2005.00943.x.
Woodward G, Ebenman B, Ernmerson M, Montoya JM, Olesen JM, Valido A, Warren PH: Body size in ecological networks. Trends Ecol Evol. 2005, 20: 402-409. 10.1016/j.tree.2005.04.005.
Fontaneto D, Herniou EA, Barraclough TG, Ricci C: On the global distribution of microscopic animals: new worldwide data on bdelloid rotifers. Zool Stud. 2007, 46 (3): 336-346.
Gomez A: Molecular Ecology of rotifers: from population differentiation to speciation. Hydrobiologia. 2005, 546: 83-99. 10.1007/s10750-005-4104-7.
Blaxter M, Mann J, Chapman T, Thomas F, Whitton C, Floyd R, Eyualem-Abebe : Defining operational taxonomic units using DNA barcode data. Phil Trans R Soc. 2005, 360 (1462): 1935-1943. 10.1098/rstb.2005.1725.
Pace N: A molecular view of microbial diversity and the biosphere. Science. 276: 734-740. 10.1126/science.276.5313.734.
Reysenbach A, Giver LJ, Wickham GS, Pace NR: Differential amplification of rRNA genes by polymerase chain reaction. Appl Environ Microbiol. 1992, 58: 3417-3418.
Schadt C, Martin AP, Lipson DA, Schmidt SK: Seasonal dynamics of previously unknown fungal lineages in tundra soils. Science. 2003, 301: 1359-1361. 10.1126/science.1086940.
Schmidt SK, Wilson KL, Meyer AF, Porter TM, Schadt CW, Moncalvo JM: The missing fungi - New insights from culture-independent molecular studies of soil. Accessing Uncultivated Microorganisms. Edited by: Karsten Zengler. 2007, Washington DC: ASM Press, 55-66.
Costello EK, Schmidt SK: Microbial Diversity in Alpine Tundra Wet Meadow Soil: Novel Chloroflexi from a Cold, Water-Saturated Environment. Environ Microbiol. 2006, 8: 1471-1486. 10.1111/j.1462-2920.2006.01041.x.
Rodrigues DF, Jesus E, Ayala-del-Río H, Pellizari VH, Gilichinksy D, Sepulveda-Torres L, Tiedje JN: Biogeography of two cold-adapted genera: Psychrobacter and Exiguobacterium. ISMEJ. 2009, 3: 658-665. 10.1038/ismej.2009.25.
Nemergut DR, Costello E, Meyer AF, Pescador MY, Weintraub MN, Schmidt SK: Structure and function of alpine and arctic soil microbial communities. Res Microbiol. 2005, 156: 775-784. 10.1016/j.resmic.2005.03.004.
Sands CJ, Convey P, Linse K, McInnes SJ: Assessing the meiofaunal variation among individuals utilizing morphological and molecular approaches: an example using the Tardigrada. BMC Ecol. 2008, 8: 7-10.1186/1472-6785-8-7.
Waite IS, O'Donnell AG, Harrison A, Davies JT, Colvan SR, Ekschmitt K, Dogan H, Wolters V, Bongers T, Bongers M, Bakonyi G, Nagy P, Papatheodorou EM, Stamou GP, Bostrom S: Design and evaluation of nematode 18S rDNA primers for PCR and denaturing gradient gel electrophoresis (DGGE) of soil community DNA. Soil Biol Biochem. 2003, 35 (9): 1165-1173. 10.1016/S0038-0717(03)00177-9.
Wu TH, Ayres E, Li G, Bardgett RD, Wall DH, Garey JR: Molecular profiling of soil animal diversity in natural ecosystems: Incongruence of molecular and morphological results. Soil Biol & Biochem. 2009, 41: 849-857. 10.1016/j.soilbio.2009.02.003.
Heyden von der S, Cavalier-Smith T: Culturing and environmental DNA sequencing uncover hidden kinetoplastid biodiversity and a major marine clade within ancestrally freshwater Neobodo designis. Int J Syst Evol Micr. 2005, 55: 2605-2621. 10.1099/ijs.0.63606-0.
Ramlov H, Westh P: Cryptobiosis in the eutardigrade Adorybiotus (Richtersius) coronifer: tolerance to alcohols, temperature and de novo protein synthesis. Zool Anz. 2001, 240: 517-523. 10.1078/0044-5231-00062.
McInnes SJ: Zoogeographic distribution of terrestrial/freshwater tardigrades from current literature. J Nat Hist. 1994, 28: 257-352. 10.1080/00222939400770131.
GBIF Data Portal. 09-24-2009, [http://www.gbif.net]
Dastych H: Redescription of Hypsibius antarcticus (Richters, 1904), with some notes on Hypsibius arcticus (Murray, 1907) (Tardigrada). Mitt Ham Zool Mus Inst. 1991, 88: 141-159.
Sands CJ, McInnes SJ, Marley NJ, Goodall-Copestake WP, Convey P, Linse K: Phylum Tardigrada: an "individual" approach. Cladistics. 2008, 24: 861-871. 10.1111/j.1096-0031.2008.00219.x.
Costello EK, Halloy SRP, Reed SC, Sowell P, Schmidt SK: Fumarole-supported islands of biodiversity within a hyperarid, high-elevation landscape on Socompa volcano, Puna de Atacama, Andes. App Env Micro. 2009, 75 (3): 735-747. 10.1128/AEM.01469-08.
Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL: GenBank. Nucleic Acids Res. 2005, 33: D34-D38. 10.1093/nar/gki063.
Ludwig W, Strunk O, Westram R, Richter L, Meier H, Yadhukumar , Buchner A, Lai T, Steppi S, Jobb G, Förster W, Brettske I, Gerber S, Ginhart AW, Gross O, Grumann S, Hermann S, Jost R, König A, Liss T, Lüßmann R, May M, Nonhoff B, Reichel B, Strehlow R, Stamatakis A, Stuckmann N, Vilbig A, Lenke M, Ludwig T, Bode A, Schleifer KH: ARB: a software environment for sequence data. Nucleic Acids Res. 2004, 32 (4): 1363-1371. 10.1093/nar/gkh293.
Porazonska DL, Giblin-Davis RM, Faller L, Farmerie W, Kanzaki N, Morris K, Powers TO, Tucker AE, Sung W, Thomas WK: Evaluating high-throughput sequencing as a method for metagenomic analysis of nematode diversity. Molec Ecol Res. 2009, 9 (6): 1439-1450. 10.1111/j.1755-0998.2009.02611.x.
Kaya M, Herniou EA, Barraclough TG, Fontaneto D: Inconsistent estimates of diversity between traditional and DNA taxonomy in bdelloid rotifers. Org, Div & Evol. 2009, 9: 3-12. 10.1016/j.ode.2008.10.002.
Ramlov H, Westh P: Survival of the cryptobiotic eutardigrade Adorybiotus coronifer during cooling to -196°C: effect of cooling rate, trehalose level and short term acclimation. Cryobiology. 1992, 29: 125-130. 10.1016/0011-2240(92)90012-Q.
Pugh PJA, McInnes SJ: The origin of Arctic terrestrial and freshwater tardigrades. Polar Biol. 1998, 19: 177-182. 10.1007/s003000050232.
Ivarsson H, Jönsson KI: Aggregation effects on anhydrobiotic survival in the tardigrade Richtersius coronifer. J Exp Zool. 2004, 301A: 195-199. 10.1002/jez.a.20018.
Pilato G, Binda MG: Biogeography of limno-terrestrial tardigrades: Are they truly incompatible binomials?. Zool Anz. 2001, 240: 511-516. 10.1078/0044-5231-00061.
Fontaneto D, Herniou HA, Boschetti C, Capriolo M, Melone G, Ricci C, Barraclough TG: Independently evolving species in asexual bdelloid rotifers. PLoS Biology. 2007, 5 (4): e87-10.1371/journal.pbio.0050087.
Fontaneto D, Barraclough TG, Chen K, Ricci C, Herniou EA: Molecular ecidence for broad-scale distributions in bdelloid rotifers: everything is not everywhere but most things are very widespread. Mol Ecol. 2008, 17: 3136-3146. 10.1111/j.1365-294X.2008.03806.x.
Nelson DR: Current Status of the Tardigrada: evolution and ecology. Integ Comp Biol. 2002, 42: 652-659. 10.1093/icb/42.3.652.
Carranza S, Giribet G, Ribera C, Baguii J, Riutort M: Evidence that Two Types of 18s rDNA Coexist in the Genome of Dugesia (Schmidtea) mediterranea (Platyhelminthes, Turbellaria, Tricladida). Mol Biol Evol. 1996, 13 (6): 824-832.
Freeman KR, Pescador MY, Reed SC, Costello EK, Robeson MS, Schmidt SK: Soil CO2 Flux and photoautotrophic community composition in high-elevation, "barren" soil. Environ Microbiol. 2009, 11 (3): 674-686. 10.1111/j.1462-2920.2008.01844.x.
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-97. 10.1093/nar/gkh340.
Winnepenninckx B, Bakeljau T, De Wachter R: Small ribosomal subunit RNA and the phylogeny of Mollusca. Nautilus. 1994, 98-100. Suppl 2
Holland ND, Campbell TG, Garey JR, Holland LZ, Wilson NG: The Florida amphioxus (Cephalochordata) hosts larvae of the tapeworm Acanthobothrium brevissime: natural history, anatomy and taxonomic identification of the parasite. Acta Zool-Stockholm. 2008, 90 (1): 75-86. 10.1111/j.1463-6395.2008.00343.x.
Nishida H, Sugiyama J: Phylogenetic relationships among Taphrina, Saitoella, and other higher fungi. Mol Biol Evol. 1993, 10 (2): 431-6.
Huber T, Faulkner G, Hugenholtz P: Bellerophon; a program to detect chimeric sequences in multiple sequence alignments. Bioinformatics. 2004, 20: 2317-2319. 10.1093/bioinformatics/bth226.
Zhang Z, Schwartz S, Wagner L, Miller W: A greedy algorithm for aligning DNA sequences. J Comput Biol. 2000, 7 (1-2): 203-14. 10.1089/10665270050081478.
Pruesse E, Quast C, Knittel K, Fuchs B, Ludwig W, Peplies J, Glöckner FO: SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nuc Acids Res. 2007, 35 (21): 7188-7196. 10.1093/nar/gkm864.
Goloboff P, Farris J, Nixon K: TNT, a free program for phylogenetic analysis. Cladistics. 2008, 24: 1-13. 10.1111/j.1096-0031.2007.00173.x. [http://www.ingentaconnect.com/content/bpl/cla/2008/00000024/00000005/art00009]
Altekar G, Dwarkadas S, Huelsenbeck JP, Ronquist F: Parallel metropolis coupled Markov chain Monte Carlo for Bayesian phylogenetic inference. Bioinformatics. 2004, 20 (3): 407-415. 10.1093/bioinformatics/btg427.
Keane TM, Naughton TJ, McInerney JO: MultiPhyl: A high-throughput phylogenomics webserver using distributed computing. Nucleic Acids Res. 2007, 35: W33-W37. 10.1093/nar/gkm359.
We thank Andrew King, Noah Fierer, David Mark Welch, and anonymous reviewers for providing helpful comments on the manuscript. Noah Fierer provided several soil samples and David Mark Welch supplied bdelloid rotifer genomic controls. The Ecology and Evolutionary Biology Department at the University of Colorado at Boulder for provided student funding for this project. This project was supported by starter funds from the NSF Microbial Observatories Program (MCB-0455606).
MSR conceived of and directed the project as well as developed the clade-specific primers. MSR, KRF, JW, & BA sampled, extracted and/or amplified and sequenced DNA from several sites or individual organisms. MSR, EKC, APM & SKS participated in the design and coordination of the study. APM, SKS, & BA guided and provided suggestions throughout the project and aided in the interpretation of the data. All authors helped to draft the manuscript. All authors read and approved the final manuscript