Skip to main content

Assessing meiofaunal variation among individuals utilising morphological and molecular approaches: an example using the Tardigrada



Meiofauna – multicellular animals captured between sieve size 45 μm and 1000 μm – are a fundamental component of terrestrial, and marine benthic ecosystems, forming an integral element of food webs, and playing a critical roll in nutrient recycling. Most phyla have meiofaunal representatives and studies of these taxa impact on a wide variety of sub-disciplines as well as having social and economic implications. However, studies of variation in meiofauna are presented with several important challenges. Isolating individuals from a sample substrate is a time consuming process, and identification requires increasingly scarce taxonomic expertise. Finding suitable morphological characters in many of these organisms is often difficult even for experts. Molecular markers are extremely useful for identifying variation in morphologically conserved organisms. However, for many species markers need to be developed de novo, while DNA can often only be extracted from pooled samples in order to obtain sufficient quantity and quality. Importantly, multiple independent markers are required to reconcile gene evolution with species evolution. In this primarily methodological paper we provide a proof of principle of a novel and effective protocol for the isolation of meiofauna from an environmental sample. We also go on to illustrate examples of the implications arising from subsequent screening for genetic variation at the level of the individual using ribosomal, mitochondrial and single copy nuclear markers.


To isolate individual tardigrades from their habitat substrate we used a non-toxic density gradient media that did not interfere with downstream biochemical processes. Using a simple DNA release technique and nested polymerase chain reaction with universal primers we were able amplify multi-copy and, to some extent, single copy genes from individual tardigrades. Maximum likelihood trees from ribosomal 18S, mitochondrial cytochrome oxidase subunit 1, and the single copy nuclear gene Wingless support a recent study indicating that the family Hypsibiidae is a non-monophyletic group. From these sequences we were able to detect variation between individuals at each locus that allowed us to identify the presence of cryptic taxa that would otherwise have been overlooked.


Molecular results obtained from individuals, rather than pooled samples, are a prerequisite to enable levels of variation to be placed into context. In this study we have provided a proof of principle of this approach for meiofaunal tardigrades, an important group of soil biota previously not considered amenable to such studies, thereby paving the way for more comprehensive phylogenetic studies using multiple nuclear markers, and population genetic studies.


Assessing variation within and between species provides information relating to taxonomic relationships [1, 2], as well as population structure, demographics and biogeographic patterns [3, 4]. Microscopic animals that are collectively known as the meiofauna present a series of challenges to studies of variation. Despite their small size (generally captured in sieve mesh sizes between 45 μm and 1000 μm), meiofauna form an integral and vital component of the soil food web, playing a fundamental role in soil ecosystem processes, particularly in nutrient recycling and decomposition processes. Twenty of the 34 recognised animal phyla have meiofaunal representatives, 5 of which are exclusively meiofaunal [5], the implication being that in a handful of soil or sediment there is likely to be a high level of biodiversity. In more extreme environments, such as those of the Antarctic, meiofauna may constitute the majority of, and in some cases the only, metoazoan element present in the soil ecosystem [68]. Meiofauna are a focus of research from many sub-disciplines (including epidemiology, ecology, soil science, agriculture, aquaculture, and pollution monitoring), as well as providing model organisms for studies of evolutionary development (the nematode Caenorhabditis elegans) and the evolution of sex (bdelloid rotifers). These tiny animals are difficult to see, often smaller than the particulate matter that comprises the substrata they inhabit, and thus difficult to separate from sediment, detritus and non-target species. Within taxonomic groupings there may be few visible distinguishing characters, and those that are present may be subtle, requiring specialist and increasingly scarce taxonomic expertise [9]. Nevertheless, accurate species identification is an essential first step to any scientific study.

Molecular phylogenetic techniques take advantage of developments in our ability to detect variation in DNA, effectively increasing the resolution available in comparison with morphological or phenotypic variation. This can be particularly useful when dealing with morphologically conserved groups. Each base change in a DNA sequence is analogous to a distinct morphological character in a phylogenetic analysis, which means many hundreds of characters can potentially be included in a single gene analysis. Models have been developed to take rates of evolution and homoplasy into account [10, 11] providing molecular phylogenetic inference with a robust theoretical basis [12]. The concept that variation in a suitable and defined length of DNA sequence may be used as a bar code for species identification is becoming increasingly appreciated [13, 14].

An important advantage of molecular data is that independent replication is possible. A phylogeny based on a single gene provides information regarding the history of and relationships between the taxa sampled. However, the information is limited to the occurrence of mutations and results in what is known as the genetree/species tree problem [15]. A significant event in the history of a lineage will only be recorded in a phylogeny if a mutation occurs at or just after the event. Conversely a high mutation rate may lead to a confused signal in cases where, for instance, an informative mutation reverts back to its ancestral state, or further changes to a state present in another lineage (homoplasy). Furthermore, there are processes, such as selection, that may affect the evolution of a gene, thereby confusing the signal of evolutionary history of the organism. These caveats may be at least partly countered by conducting independent analyses on several unlinked genes in order to generate a general consensus of phylogenies that more closely reflects the "true" evolutionary history, or by concatenating sequences to produce a super phylogeny [16].

It is often not appreciated that transferring molecular techniques from the few "model" organisms that are the focus of intense laboratory research to the more numerous but less studied groups is not straightforward. Conversely, molecular ecologists who work on meiofauna and other invertebrates may be surprised when they find their difficulties are not appreciated by the wider scientific community. DNA quantity is a major issue with tiny organisms. To increase DNA yield, the whole organism or, more often several pooled individuals, are used in extractions and, thus, problems with contaminants from gut contents or commensal organisms arise [17, 18]. The difficulties working with understudied species are often only evidenced by the lack of published literature available. For example, as far as we are aware, there have only been three assessments made of intraspecific variation in any species of the meiofaunal phylum Tardigrada [1921]. This may be partially due to tardigrades and other meiofaunal groups falling outside economic and charismatic categories that influence the direction of science programs, but it is also likely to be due to the difficulties involved in obtaining genetic data from these members of the meiofaunal community. Recently a variety of techniques have been developed or demonstrated as suitable for DNA extraction from a variety of meiofauna [2, 2225] providing a basis for basic evolutionary studies and the development of molecular tools for identification or classification [20, 22, 23, 26].

Tardigrades are found in most terrestrial, freshwater and marine habitats, including some remote Antarctic nunataks where they are found in the absence of the otherwise ubiquitous nematodes [6]. They are one of the few phyla with representatives found from the highest and coldest terrestrial environments to the deepest oceans. They have featured in high impact publications mostly due to their uncertain phylogenetic position in the tree of life [2730], although their ubiquitous distribution makes them ideal candidates for historical biogeographic reconstruction [3133]. Molecular phylogenetic work to date is strongly concordant with morphological based systematics [29], supporting both tardigrade monophyly and monophyly in the constituent classes Eutardigrada and Heterotardigrada [17, 24, 34].

In this paper we describe a method that enables the quantitative and qualitative assessment of morphology and genetic variation among individual tardigrades. The method includes a novel technique for separating all organisms from their substrate as compared with "cherry-picking" visible organisms under a dissecting microscope, a simple DNA release technique applicable at the individual level, and a general protocol for amplifying genomic DNA from multicopy and single copy genes. Previous studies have amplified multi-copy 18S and high-copy number (CO1) genes from tardigrades [2, 20, 25, 26], but this is the first time single copy nuclear genes have been amplified from genomic DNA extracted from a single tardigrade. We demonstrate that this method allows for informed re-assessment of morphological variation, allows independent replication of phylogenetic analyses and is suitable for assessing within population variation enabling population genetic studies.

Results and Discussion

Sample collection and preparation

Separating meiofauna from the substrate can be done in two ways: "cherry-picking" individuals using a pipette or Irwin loop under a microscope, or mechanically separation, usually using a density gradient. "Cherry-picking" is a fast technique if only a few specifically targeted individuals are required and the organisms are easily identifiable from the substrate. However, if quantitative biomass and diversity data are required this technique will be biased towards larger, vagile and more or less pigmented organisms (depending upon which contrasts with the substrate). Where quantitative and qualitative results are required to extract all individuals from a substrate the cherry-picking approach is inappropriate and such studies have adopted mechanical separation using density gradients [see [35]]. However, the media used, such as Ludox AM (Dupont, France) [3537], Percoll (Pharmacia, Uppsala, Sweden) [38] or 50% sucrose [39] may have detrimental downstream effects, particularly for PCR. Using OptiPrep™ – Density Gradient Media (Axis-Shield, UK) we were able to isolate and identify all individual tardigrades and eggs from habitat substrata (see methods). Substrata used included fresh, frozen and dried herbarium moss specimens. Contrasting with previous protocols OptiPrep is non-ionic, non-toxic and does not require washing to remove the media. It has been used to fractionate sub-cellular organelles and does not interfere with marker enzyme activities, allowing fractions of cell organelles to be analyzed without removal of the gradient medium [40]. In the current study it proved to be clean and efficient at extracting meiofauna from substrata, allowing biodiversity and biomass from each sample to be determined, and detailed morphological examination without affecting subsequent DNA extraction or downstream biochemical processes. Furthermore, tardigrades extracted from fresh and frozen samples remained alive after the extraction process.

DNA extraction

Taking into account the limited amount of DNA available in a single tardigrade we opted for a DNA release method rather than an extraction method to maximise the amount of DNA available. We trialled several different release methods (NaOH digestion, [22]; TE/ddH2O boil, [41]; Proteinase K freeze thaw cycles, [23]). All these techniques were sufficient for amplification of 18S rDNA, as many more expensive commercial kits have proven to be [24]. However, we found a 20 min boil in 40 μL of 5% chelex the most reliable technique across all gene regions, possibly because the chelex beads inactivate inhibitors that would otherwise prevent reliable amplification. Obviously the risk of contamination and competition in PCR is significant and all precautions are advised to prevent inadvertent introduction of foreign DNA. Fortuitously this technique is cheap and rapid and the only reagent necessary (Chelex 100 – BIORAD) is inexpensive and readily available in many molecular laboratories.

DNA quantification

In order to quantify DNA concentration we used a PicoGreen™ assay (Molecular Probes), a sensitive technique for determining concentrations of DNA in a solution. After extensive optimization to reduce the volumes of template required and level of background fluorescence, we were unable to detect DNA in any of our extractions. However, we were able to successfully amplify ribosomal 18S and mitochondrial cytochrome c oxidase subunit 1 (CO1) products from most of our tardigrade extractions – even from animals extracted from rehydrated herbarium samples that had been stored dry for 15 years. Recently Kiel et al [42] extracted 8.4 μg of genomic DNA from ~2000 pooled Hypsibius klebelsbergi, indicating that each tardigrade contributed ~4.2 ng of DNA to the extract. For amplification using polymerase chain reaction, 20 – 50 ng of DNA is usually recommended in a 1 μL volume. Clearly this will never be possible at individual level for many meiofaunal groups – particularly tardigrades – making studies of individuals rather than pooled samples challenging. Even with the advent of whole genome amplification, to achieve representative amplification of the genome a minimum starting quantity of 1 μl of 10 ng/μL is recommended (GE HealthCare: GenomiPhi™ instruction manual).

Polymerase Chain Reaction (PCR)

Our strategy for amplifying genomic DNA was to use nested or hemi-nested PCR (see Table 1 for primer sequences). Rarely was product visible on a gel after a single round of PCR. However, re-amplifying from the first reaction using at least one internal primer resulted in reasonable amplification of the desired product. Large gene regions (e.g. 18S) were amplified in overlapping fragments; first the whole region was amplified, then using this reaction mix as template three overlapping fragments were amplified. Products amplified in this way generally returned clean sequences. Two of the three 18S fragments amplified reliably on all templates, fragment one was difficult to amplify in heterotardigrades, while CO1 was less reliable and Wingless amplified only template from fresh or recently frozen samples. We were also able to amplify Alpha Spectrin (aspec) and Elongation Factor 1 alpha (EF1 a) from some individuals. We were able to verify that the aspec sequences we obtained were homologues but we were unable to verify these as originating from Tardigrada. We were able to obtain EF1a sequences from most individuals, however, up to 8 different paralogues were identified, each with remarkably conserved coding sequence making orthologue-specific primer redesign difficult (data not presented). GenBank accession numbers and sample details can be found in the additional information file associated with this paper.

Table 1 Primers used in each reaction combination

Phylogenetic reconstruction

Using maximum likelihood methods we have constructed phylogenies from nuclear ribosomal 18S, mitochondrial CO1 and nuclear Wingless regions (Figures 1 and 2). As we were unable to amplify heterotardigrades using Wingless primers we used Milnesium sp. as an out-group in this case as it is believed to represent a basal eutardigrade [17, 23, 28]. For the 18S and CO1 phylogenies we used heterotardigrades as out-groups to root the trees and explore the relationships within the eutardigrades. Although all studies to date have relied on a small number of taxa (11 tardigrade taxa in [17]), results of both molecular and morphological work indicate monophyly of heterotardigrades and eutardigrades, and Order Apochaela holding a basal position within the Eutardigrada [17, 23, 28]. Our results are consistent with these findings. The most striking feature contained in these phylogenies is the non-monophyletic relationship between the eutardigrade families Hypsibiidae and Macrobiotidae and Murrayidae, with Macrobiotidae and Murrayidae a clade nested within Hypsibiidae. Although not discussed, this is also evident in the phylogeny of Jørgensen and Kristensen [17] and has recently been independently identified by Kiehl et al 07 [34]. By utilizing our current method, further work is planned to fully explore phylogenetic relationships among and within tardigrade families.

Figure 1
figure 1

Tardigrade phylogeny constructed from 18S rRNA. Phylogeny produced using maximum likelihood analyses under a GTR+I+Γ model using 884 bases of 18S rRNA. Numbers at nodes are support values generated from 1000 bootstrap pseudoreplicates. Terminal labels associated with unique numbers were identified as: Macro 001–008 – Macrobiotus furciger, Macro 048-52,68-077 – Macrobiotus sp. (hufelandi type), Dacty – Dactylobiotus sp., Acut – Acutuncus antarcticus, Diph – Diphascon sp., RamO – Ramazzottius oberhaeuseri, Isohyp – Isohypsibius asper, MilnT – Milnesium tardigradum, MilnC – Milnesium sp. "charcot", Bryod – Bryodelphax sp., EchinC – Echiniscus sp., EchinT – Echiniscus testudo.

Figure 2
figure 2

Tardigrade phylogeny constructed from CO1 and Wingless. Phylogenies produced using maximum likelihood analyses under GTR+I+Γ models using: A. 534 bases of mitochondrial CO1 and B. 321 bases of Wingless. Terminal labels correspond to those described in Figure 1. Support at nodes was generated using 1000 bootstrap pseudoreplicates.

Identifying molecular variation

An important advantage of analysing sequences from individuals rather than a pooled sample is that variation that is difficult to distinguish using morphology is likely to be identified. For example, Macrobiotus sp. is clustered into two groups, one containing individuals from near Ville de Jumelles, St Maur Créteil, France, the other containing individuals from Charcot Island, Antarctica (Figure 1, Additional File 1). The individual Macro-048 was obtained from near Ville de Jumelles, St Maur Créteil, and yet is clearly different to other individuals identified as conspecifics from the same sample of moss. Similarly Acutuncus sp. was collected from Jubany Station, King George Island, Antarctica. There were clearly two different taxa represented at this one site – one taxon likely to be Hypsibius.

Misidentification in these kinds of studies is not surprising as identification can only be carried out under low power in order not to compromise material for subsequent molecular study (we used 400× inverted microscope). Rigorous taxonomic identification requires permanent mounting and viewing under 1000× magnification (see Methods regarding voucher specimens). Thus the limitation of potential misidentification is common to all current molecular studies of tardigrades. Our protocol allows for individual molecular variation to be assessed. Then, where unexpected variation occurs, an informed re-assessment of morphology can be made. To demonstrate this, we took a second sample of substratum from Jubany Station, re-extracted the tardigrades present and mounted multiple individuals to be observed under high power magnification. In this sample Acutuncus antarcticus and two undescribed species of Hypsibius were identified (data not presented). Furthermore an egg of Acutuncus antarcticus, (for which identification is considerably easier due to ornate sculpturing on the case) was taken from the sample and DNA extracted, amplified and sequenced, and the sequence then used to verify which of our samples were Acutuncus and which were the unknown (labelled as Hyp/Acut in Figures 1 and 2 to indicate the uncertain identification).

One other study has used individual tardigrades to explore diversity in environmental samples [20]. In this study remarkable unexpected diversity was present, even between morphologically similar groups, which would have been confused if samples had have been pooled. Pooling mixed taxon samples would at best mask the diversity present or, at worst, generate false sequence due to PCR recombination [43]. Indeed, PCR recombination was detected in our study. A blast search of the full sequences of the two aberrant Echiniscus testudo 18S specimens (EchinT 058 and EchinT 060: Fig 1) closely matched Echiniscus sequence previously deposited on GENBANK. However, blasting a small (~100 base) anomalous region of the sequence perfectly matched the yeast Candida sp. indicating that the variation in these PCR products was a chimeric artefact generated from two competing templates.

Utility for population comparisons

Traditionally in population studies a suite of variable markers are developed and used to screen many individuals to explore intraspecific population structure using allele frequency (genic) and genotype frequency (genotypic) approaches. More recently, studies have taken advantage of the genealogical information available in sequence data to investigate contemporary and historical population processes [3, 4447]. Our method potentially enables multiple gene sequences to be obtained from individuals allowing genic, genotypic and genealogic analyses [48]. There has only been a single molecular population genetic study of tardigrades [21]. This study provides an important contribution to the field as it indicates that Echiniscus testudo has a high migration potential. Importantly though, the data were collected from pooled individuals and are based on a single gene (CO1) which may not accurately reflect the true population structure [49, 50]. To thoroughly explore population structure, multiple independent markers are required and intrapopulation variation needs to be evaluated [51]. From our example we can demonstrate intraspecific variation in Dactylobiotus sp. samples in all three gene regions examined (the 8 Dactylobiotus specimens came from a single portion of moss and are presumed to be conspecifics) and 5 sequence variable alleles at the Wingless locus for Ramazzotius oberhaeuseri from a single clump of moss.


We have described a method that allows for extraction of individual meiofaunal specimens from their substratum. Our protocol for assessing between-individual variation is guided by but independent of morphology. We have demonstrated that phylogenetic analyses using data gained from individuals rather than pooled samples allow the detection of hidden diversity and assists in distinguishing true diversity from errors. Sequence data obtained from multiple genes will also have application in improving resolution in population level studies. We are currently employing the techniques described here to conduct large-scale phylogenetic and biogeographic studies of tardigrades.


Sample collection and specimen isolation

Fresh substrata samples of mosses and lichens were collected from walls in St. Maur des Fossés, Paris, France, a type locality for several tardigrade species. Herbarium and frozen samples from sub- and maritime Antarctic islands and from mainland Antarctic nunataks were also examined (See Additional File 1 for details concerning species, locations and genbank accession numbers). Samples were soaked in double-distilled water and transferred to a mini cup blender (Waring). Individual samples of approximately 2 cm3 were homogenised briefly (two short 2 second bursts at the lower speed) then filtered and washed with double-distilled water through a coarse (~1 mm) mesh sieve. The resulting sediment and water were transferred to 1.5 mL micro-centrifuge tubes prepared with OptiPrep™ – Density Gradient Media (Axis-Shield, UK). A layer of 50% OptiPrep™ solution was floated over a 100% solution and the sample added above these layers. The tubes were spun in a bench-top centrifuge (Griffin) at 1000 rpm for 1 minute. The material from the interface between sample water and the 50% OptiPrep™ solution was removed and filtered through a 45 μm mesh sieve, using double-distilled water. The sieve was washed out into a 50 mL Petri dish and viewed under a Wild M5 dissecting microscope. Tardigrades, and eggs where present, were lifted from the Petri dish via Irwin loop (a non-corrosive nickel-chromium wire loop, approx. 200 μm × 500 μm) into a cavity slide filled with double-distilled water. The slide was viewed on a Nikon Diaphot inverted microscope with a 40 to 400× magnification range to identify the tardigrades to genus, species, or type. The individual tardigrades, or eggs, were lifted from the slide via Irwin loop and transferred to 5 μL of double-distilled water in a labelled 0.5 mL microcentrifulge tube and frozen at -80°C.

Voucher specimens

In any given moss several different taxa were often present. Representatives of each were mounted directly onto slides with de Faure's mounting media [52], and individuals, where possible, were identified to species under high magnification (Olympus BX50 – max. 1000× magnification). These voucher specimens were deposited at the British Antarctic Survey Data Resource Centre.

DNA isolation and PCR conditions

Individual specimens were subjected to two rounds of freeze-thaw cycling (-80°C to 55°C) in sealed 0.5 mL microcentrifuge tubes to assist in rupturing the cuticle. To each microcentrifuge tube 40 μL of vortexed 5% Chelex 100 (Biorad) was added and the specimens incubated at 99°C for 20 min. Samples were centrifuged for 1 min at maximum speed in a microcentrifuge (Eppendorf) and stored at -20°C. A blank (5 μL H2O, 40 μL 5% chelex) was treated in the same way and used as a negative control in PCR.

We attempted to quantify the amount of DNA present in the extract using a PicoGreen™ (Molecular Probes) assay. We optimised the assay for 10 μL total volume (5 μL of sample) to minimise background fluorescence with a standard ranging from 1 ng μL-1 to 0.00001 ng μL-1 to be detected on a Q-PCR thermo-cycler (Stratagene). The standard curve detected DNA fluorescence to 0.0005 ng μL-1 after which the relationship between fluorescence and DNA concentration was no longer linear. Our samples did not differ from background fluorescence (<0.0001 ng μL-1).

To increase the yield of PCR products we labelled each primer with an M13 tail, the forward primers with M13 REV and the reverse primers with M13 -21 [53]. Primer details are given in Table 1.

We amplified the nuclear "multi-copy" gene 18S rDNA (18S) and the mitochondrial (thus high copy number) cytochrome c oxidase subunit 1 (CO1) gene for which there were data available on GENBANK to verify that the products we were amplifying were from the target organisms as opposed to contamination from airborne substrata, gut contents or commensal organisms associated with the tardigrade cuticle [17]. For 18S we used the primers SSU 1F and modified SSU 82R [54, 55] to amplify approximately 1800 bases in 10 μL volumes containing 16 μM ammonium sulphate, 68 mM Tris-HCl (pH 8), 10 mM β-mercaptoethanol, 5% bovine serum albumin 10 mg/mL (Sigma), 3 mM magnesium chloride, 200 μM each dNTP, 0.5 μM each primer, 0.5 units of Taq DNA Polymerase (Bioline), and 1 μL template DNA. Cycling conditions were 94°C for 2 min followed by 35 cycles of 94°C for 1 min, 60°C for 30 sec and 72°C for 1 min. This was followed by a 4 min extension at 72°C.

Using 2 μL of this reaction mixture we conducted three separate reactions to amplify three overlapping fragments of the gene using the primers SSU 1F/SSU 26R, SSU 22F/SSU 13R and SSU 26F/SSU 82R [55] (see Table 1). All three fragments were amplified in 40 μL volumes using the above conditions with the exception of magnesium chloride concentration, which was lowered to 2 mM, and annealing temperature, which was raised to 65°C, to increase specificity.

The mitochondrial CO1 gene was amplified by hemi-nested PCR. The first amplification was with the primers LCO_1490 [56] and mtD9 [57] in 10 μL volumes containing 16 μM ammonium sulphate, 68 mM Tris-HCl (pH 8), 10 mM β-mercaptoethanol, 5% bovine serum albumin 10 mg/mL (Sigma), 2 mM magnesium chloride, 200 μM each dNTP, 0.5 μM each primer, 0.5 units of Taq DNA Polymerase (Bioline), and 1 μL template DNA. Cycling conditions were 94°C for 2 min followed by 35 cycles of 94°C for 30 sec, 45°C for 30 sec and 72°C for 45 sec. This was followed by a 4 min extension at 72°C.

The second amplification was with the primers LCO_1490 and HCO_2198 [56]. Reaction conditions were in 40 μL volumes using the above concentrations. Cycling conditions were as above but with annealing temperature raised to 50°C.

A similar nested amplification strategy was applied to the genes Wingless, Alpha Spectrin and Elongation Factor 1 alpha, using the initial PCR cycling conditions suggested by Regier [58]. Amplifications were conducted in 10 μL (first amplification) and 40 μL (second amplification) reactions using reaction mix concentrations as above. The cycling conditions for the first amplification were 94°C for 2 min followed by 24 cycles starting 94°C 30 sec, 56°C 30 sec, 72°C 1 min but decreasing annealing temperature from 56°C to 45°C by 0.4°C each cycle and increasing extension time by 2 seconds each cycle. This was followed by a further 12 cycles of 94°C 30 sec, 45°C 30 sec 72°C 2 min increasing the extension time by 3 seconds each cycle. The reaction was terminated with a 4 min extension at 72°C. The second amplification used 2 μL of the first reaction as template but used more conventional cycling conditions of 94°C for 2 min followed by 35 cycles of 94°C for 30 sec, 45°C for 30 sec and 72°C for 45 sec. This was followed by a 4 min extension at 72°C. Amplified products were sequenced by Macrogen inc. Seoul, Korea. All sample information including GenBank accession numbers are available as an additional file.


Sequence trace files were base-called and aligned using Codoncode Aligner V1.6.3 (CodonCode Corp). Coding genes were checked for open reading frame and blast searched (tblastx) to assess gene homology. Where highly divergent nucleotide sequences proved difficult to align, sequences were converted to amino acids and aligned by eye and back translated to nucleotide in SE-AL [59]. Ambiguous alignment regions were excised from the analyses.

To determine an appropriate model of sequence evolution we used the iterative optimisation approach suggested by Swofford [60] and validated by Sullivan [61]. The data was used to construct a neighbour joining (NJ) tree in PAUP*4 10b [62] from which the likelihood parameters were estimated and used in a heuristic likelihood tree search (NJ starting tree, TBR). After a minute the search was stopped and likelihood parameters were estimated from the new tree. These parameters were used to conduct another likelihood search that was again stopped after a minute. This procedure was repeated until the likelihood scores stabilised. Once parameters were optimised for the full model, we systematically reduced the number of parameters to see if a simpler model could be used to describe the patterns of evolution without adversely affecting the overall likelihood scores. Using the appropriate model (GTR+I+Γ) we conducted heuristic searches with 100 random starting addition sequences using maximum likelihood option in PAUP*. We conducted 1000 bootstrap replicates (1 random addition sequence) to assess the strength of each phylogenetic inference.


  1. Guidetti R, Rebecchi L, Bertolani R: Cuticle structure in the systematics of Macrobiotidae (Tardigrada, Eutardigrada). Acta Zool. 2000, 81: 27-36. 10.1046/j.1463-6395.2000.00034.x.

    Article  Google Scholar 

  2. Guidetti R, Gandolfi A, Rossi V, Bertolani R: Phylogenetic analysis of Macrobiotidae (Eutardigrada, Parachela): a combined morphological and molecular approach. Zoologica Scripta. 2005, 34: 235-244. 10.1111/j.1463-6409.2005.00193.x.

    Article  Google Scholar 

  3. Garrick RC, Sands CJ, Rowell DM, Tait NN, Greenslade P, Sunnucks P: Phylogeography recapitulates topography: very fine-scale local endemism of a saproxylic 'giant' springtail at Tallaganda in the Great Dividing Range of south-east Australia. Molecular Ecology. 2004, 13: 3315-3330. 10.1111/j.1365-294X.2004.02340.x.

    Article  Google Scholar 

  4. Sunnucks P, Blacket MJ, Taylor JM, Sands CJ, Ciavaglia SA, Garrick RC, Tait NN, Rowell DM, Pavlova A: A tale of two flatties: different responses of two terrestrial flatworms to past environmental climatic fluctuations at Tallaganda in montane south-eastern Australia. Molecular Ecology. 2006, 15: 4513-4531. 10.1111/j.1365-294X.2006.03107.x.

    CAS  Article  PubMed  Google Scholar 

  5. Higgins RP, Thiel H: Introduction to the study of meiofauna. 1988, Smithsonian Institution Press. Washington DC, USA, 488-

    Google Scholar 

  6. Convey P, McInnes SJ: Exceptional tardigrade-dominated ecosystems in Ellsworth Land, Antarctica. Ecology. 2005, 86: 519-527. 10.1890/04-0684.

    Article  Google Scholar 

  7. Maslen NR, Convey P: Nematode diversity and distribution in the southern maritime Antarctic – clues to history?. Soil Biology and Biochemistry. 2006, 38: 3141-3151. 10.1016/j.soilbio.2005.12.007.

    CAS  Article  Google Scholar 

  8. Convey P: Antarctic Ecosystems. Encyclopedia of biodiversity. Edited by: Levin S. Elsevier,

  9. Godfray HCJ: Challenges for taxonomy. Nature. 2002, 417: 17-19. 10.1038/417017a.

    CAS  Article  PubMed  Google Scholar 

  10. Yang Z: Maximum-likelihood models for combined analysis of multiple sequence data. Journal of Molecular Evolution. 1996, 42: 587-596. 10.1007/BF02352289.

    CAS  Article  PubMed  Google Scholar 

  11. Yang Z: Among site rate variation and its impact on phylogenetic analyses. Trends in Ecology and Evolution. 1996, 11: 367-372. 10.1016/0169-5347(96)10041-0.

    CAS  Article  PubMed  Google Scholar 

  12. Sullivan J, Joyce P: Model selection in phylogenetics. Annual Review of Ecology Evolution and Systematics. 2005, 36: 445-466. 10.1146/annurev.ecolsys.36.102003.152633.

    Article  Google Scholar 

  13. Blaxter M, Mann J, Chapman T, Thomas F, Whitton C, Floyd R, Abebe E: Defining operational taxonomic units using DNA barcode data. Philosophical Transactions of the Royal Society of London Series B. 2005, 360: 1935-1943. 10.1098/rstb.2005.1725.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  14. Hebert PDN, Cywinska A, Ball SL, DeWaard JR: Biological identifications through DNA barcodes. Proceedings of the Royal Society of London Series B Biolological Sciences. 2003, 270: 313-321. 10.1098/rspb.2002.2218.

    CAS  Article  Google Scholar 

  15. Page RDM, Charleston MA: From gene to organismal phylogeny: Reconciled trees and the gene tree/species tree problem. Molecular Phylogenetics and Evolution. 1997, 7: 231-240. 10.1006/mpev.1996.0390.

    CAS  Article  PubMed  Google Scholar 

  16. Gadagkar DR, Rosenberg MS, Kumar S: Inferring species phylogenies from multiple genes: Concatenated sequence tree versus consensus gene tree. Journal of Experimental Zoology. 2004, 304: 64-74.

    Google Scholar 

  17. Jørgensen A, Kristensen RM: Molecular phylogeny of Tardigrada – investigation of the monophyly of Heterotardigrada. Molecular Phylogenetics and Evolution. 2004, 32: 666-670. 10.1016/j.ympev.2004.04.017.

    Article  PubMed  Google Scholar 

  18. Bridge PD, Roberts PJ, Spooner BM, Pancahl G: On the unreliability of published DNA sequences. New Phytologist. 2003, 160: 43-48. 10.1046/j.1469-8137.2003.00861.x.

    CAS  Article  Google Scholar 

  19. Rebecchi L, Rossi V, Altiero T, Bertolani R, Menozzi P: Reproductive modes and genetic polymorphism in the tardigrade Richtersius coronifer (Eutardigrada, Macrobiotidae). Invertebrate Biology. 2003, 122: 19-27.

    Article  Google Scholar 

  20. Blaxter M, Elsworth B, Daub J: DNA taxonomy of a neglected animal phylum: an unexpected diversity of tardigrades. Proceedings of the Royal Society of London B. 2004, 271: 189-192. 10.1098/rsbl.2003.0130.

    Article  Google Scholar 

  21. Jørgensen A, Møbjerg N, Kristensen RM: A molecular study of the tardigrade Echiniscus testudo (Echiniscidae) reveals low DNA sequence diversity over a large geographic area. Proceedings of the Tenth International Symposium on Tardigrada, Journal of Limnology. 2007, 66: 77-83.

    Google Scholar 

  22. Floyd R, Eyualem A, Papert A, Blaxter M: Molecular barcodes for soil nematode identification. Molecular Ecology. 2002, 11: 839-850. 10.1046/j.1365-294X.2002.01485.x.

    CAS  Article  PubMed  Google Scholar 

  23. Leutbecher C: A routine method of DNA-extraction from extremely small metazoans, e.g. single rotifer specimens for RAPD-PCR analysis. Hydrobiologia. 2004, 437: 133-137. 10.1023/A:1026546705582.

    Article  Google Scholar 

  24. Nichols PB, Nelson DR, Garey JR: A family level analysis of tardigrade phylogeny. Hydrobiologia. 2006, 558: 53-60. 10.1007/s10750-005-1414-8.

    Article  Google Scholar 

  25. Schill RO: Comparison of different protocols for DNA preparation and PCR amplification of mitochondrial genes of tardigrades. Proceedings of the Tenth International Symposium on Tardigrada, Journal of Limnology. 2007, 66: 164-170.

    Google Scholar 

  26. Schill RO, Steinbrück G: Identification and differentiation of Heterotardigrada and Eutardigrada species by riboprinting. Journal of Zoological Systematic Evolutionary Research. 2007, 10.1111/j.1439-0469.2007.00409.x.

    Google Scholar 

  27. Philippe H, Lartilot N, Brinkmann H: Multigene analyses of bilaterian animals corroborate the monophyly of Ecdysozoa, Lophotrochozoa and Protostomia. Molecular Biology and Evolution. 2005, 22: 1246-1253. 10.1093/molbev/msi111.

    CAS  Article  PubMed  Google Scholar 

  28. Garey JR, Krotec M, Nelson DR, Brooks J: Molecular analysis supports a tardigrade arthropod association. Invertebrate Biology. 1996, 115: 79-88. 10.2307/3226943.

    Article  Google Scholar 

  29. Garey JR, Nelson DR, Mackey LY, Li J: Tardigrade phylogeny: congruency of morphological and molecular evidence. Zoologischer Anzeiger. 1999, 238: 205-210.

    Google Scholar 

  30. Giribet G, Carranza S, Baguna J, Riutort M, Ribera C: First molecular evidence for the existence of a Tardigrada+Arthropoda clade. Molecular Biology and Evolution. 1996, 13: 76-84.

    CAS  Article  PubMed  Google Scholar 

  31. McInnes SJ, Pugh PJA: Biogeography of limno-terrestrial Tardigrada, with particular reference to the Antarctic fauna. Journal of Biogeography. 1998, 25: 31-36. 10.1046/j.1365-2699.1998.251176.x.

    Article  Google Scholar 

  32. McInnes SJ, Pugh PJA: An attempt to revisit the global biogeography of limno-terrestrial Tardigrada. Proceedings of the Tenth International Symposium on Tardigrada, Journal of Limnology. 2007, 66: 90-96.

    Google Scholar 

  33. Pilato G, Binda MG: Biogeography and limno-terrestrial tardigrades: are they truly incompatible binomials?. Zoologischer Anzeiger. 2001, 240: 511-516. 10.1078/0044-5231-00061.

    Article  Google Scholar 

  34. Kiehl E, Dastych H, D'Haese J, Greven H: The 18S rDNA sequences support polyphyly of the Hypsibiidae (Eutardigrada). Proceedings of the Tenth International Symposium on Tardigrada, Journal of Limnology. 2007, 66: 21-25.

    Google Scholar 

  35. Bartels PJ, Nelson DR: A large-scale, multihabitat inventory of the phylum Tardigrada in the Great Smoky Mountains National Park, USA: A preliminary report. Hydrobiologia. 2006, 558: 111-118. 10.1007/s10750-005-1405-9.

    Article  Google Scholar 

  36. Nelson DR, McInnes SJ: Tardigrades. Freshwater Meiofauna: Biology and Ecology. Edited by: Rundle SD, Robertson AL, Schmid-Araya JM. 2002, Buckhuys, Leiden, 177-215.

    Google Scholar 

  37. Burgess R: An improved protocol for separating meiofauna from sediments using colloidal silica sols. Marine Ecology Progress Series. 2001, 214: 161-165. 10.3354/meps214161.

    Article  Google Scholar 

  38. Bernhard JM, Buck KR, Farmer MA, Bowser SS: The Santa Barbara Basin is a symbiosis oasis. Nature. 2000, 403: 77-80. 10.1038/47476.

    CAS  Article  PubMed  Google Scholar 

  39. Hallas TE: A mechanical method for the extraction of Tardigrada. Memorie dell'Istituto Italiano di Idrobiologia. 1975, 32 (Suppl): 153-157.

    Google Scholar 

  40. Graham J, Ford T, Rickwood D: The preparation of subcellular organelles from mouse liver in self-generated gradients of iodixanol. Analytical Biochemistry. 1994, 220: 367-73. 10.1006/abio.1994.1351.

    CAS  Article  PubMed  Google Scholar 

  41. Hosek J, Svastova P, Moravkova M, Pavlik I, Bartos M: Methods of mycobacterial DNA isolation from different biological material: a review. Veterinarni Meicina. 2006, 51: 180-192.

    CAS  Google Scholar 

  42. Kiehl E, Dastych H, D'Haese J, Greven H: A cDNA library of the eutardigrade Hyspsibius klebelsbergi Mihelčič, 1959 and analysis of the actin gene. Proceedings of the Tenth International Symposium on Tardigrada, Journal of Limnology. 2007, 66: 152-157.

    Google Scholar 

  43. Cronn R, Cedroni M, Haselkorn T, Grover C, Wendel JF: PCR-mediated recombination in amplification products derived from polyploidy cotton. Theoretical and Applied Genetics. 2002, 104: 482-489. 10.1007/s001220100741.

    CAS  Article  PubMed  Google Scholar 

  44. Templeton AR, Routman E, Phillips CA: Separating population substructure from population history: a cladistic analysis of the geographical distribution of mitochondrial DNA haplotypes in the tiger salamander, Ambystoma tigrinum. Genetics. 1996, 140: 767-782.

    Google Scholar 

  45. Castens BC, Stevenson AL, Degenhardt JD, Sullivan J: Testing nested phylogenetic and phylogeographic hypotheses in the Plethodon vandykei species group. Systematic Biology. 2004, 53: 781-792. 10.1080/10635150490522296.

    Article  Google Scholar 

  46. Castens BC, Brunsfeld SJ, Demboski JR, Good JM, Sullivan J: Investigating the evolutionary history of the Pacific northwest mesic forest ecosystem: hypothesis testing within a comparative phylogeographic framework. Evolution. 2005, 59: 1639-1652.

    Google Scholar 

  47. Garrick RC, Sands CJ, Rowell DM, Hillis DM, Sunnucks P: Catchments catch all: long-term population history of a giant hexapod from the southeast Australian highlands – a multi-gene approach. Molecular Ecology. 2007, 16: 1865-1882. 10.1111/j.1365-294X.2006.03165.x.

    CAS  Article  PubMed  Google Scholar 

  48. Garrick RC, Sunnucks P: Development and application of three-tiered nuclear genetic markers for basal Hexapods using single-stranded conformation polymorphism coupled with targeted DNA sequencing. BMC Genetics. 2006, 7: 11-10.1186/1471-2156-7-11.

    PubMed Central  Article  PubMed  Google Scholar 

  49. Shaw K: Conflict between nuclear and mitochondrial DNA phylogenies of a recent species radiation: what mtDNA reveals and conceals about modes of speciation in Hawaiian crickets. Proceedings of the National Acadamy of Science of the United States of America. 2002, 99: 16122-16127. 10.1073/pnas.242585899.

    CAS  Article  Google Scholar 

  50. McCracken KG, Sorenson MD: Is homoplasy or lineage sorting the source of incongruent mtDNA and nuclear gene trees in the stiff-tailed ducks (Nomonyx-Oxyura)?. Systematic Biology. 2005, 54: 35-55. 10.1080/10635150590910249.

    Article  PubMed  Google Scholar 

  51. Sunnucks P: Efficient genetic markers for population biology. Trends in Ecology and Evolution. 2000, 15: 199-203. 10.1016/S0169-5347(00)01825-5.

    Article  PubMed  Google Scholar 

  52. Gray P: The Microtomist's Formulary and Guide. 1954, New York, Blakiston Co. Inc

    Google Scholar 

  53. Regier JC, Shi D: Increased yield of PCR product from degenerate primers with nondegenerate, nonhomologous 5' tails. BioTechniques. 2005, 38: 34-38.

    CAS  Article  PubMed  Google Scholar 

  54. Medlin L, Elwood HJ, Stickel S, Sogin ML: The characterization of enzymatically amplified eukaryotic 16S-like rRNA-coding regions. Gene. 1988, 71: 491-499. 10.1016/0378-1119(88)90066-2.

    CAS  Article  PubMed  Google Scholar 

  55. Nematode Phylogenetics Primers. []

  56. Folmer O, Black M, Hueh W, Lutz R, Vrijenhoek R: DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Molecular Marine Biology and Biotechnology. 1994, 3: 294-299.

    CAS  PubMed  Google Scholar 

  57. Simon C, Frati F, Beckenbach A, Crespi B, Liu H, Flook P: Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Annals of the Entomological Society of America. 1994, 87: 651-701.

    CAS  Article  Google Scholar 

  58. Regier JC: Protocols, Concepts, and Reagents for preparing DNA sequencing templates. []

  59. Rambaut A: Se-Al: Sequence Alignment Editor. []

  60. Swofford DL, Olsen GJ, Waddell PJ, Hillis DM: Phylogenetic inference. Molecular Systematics. Edited by: Hillis DM, Moritz C, Mable B. 1996, Massachusetts. Sinauer Associates, Sunderland, 407-514. 2

    Google Scholar 

  61. Sullivan J, Abdo Z, Joyce P, Swofford D: Evalutating the performance of a successive-approximations approach to parameter optimization in maximum-likelihood phylogeny estimation. Molecular Biology and Evolution. 2005, 22: 1386-1392. 10.1093/molbev/msi129.

    CAS  Article  PubMed  Google Scholar 

  62. Swofford DL: PAUP*: Phylogenetic Analysis Using Parsimony (* and Other Methods). Version 4. 2002, Sinauer Associates, Sunderland, Massachusetts

    Google Scholar 

Download references


We thank Will Goodall-Copesake and Nigel Marley for discussions that led to improvements to the text. We thank the authors of [34] for sharing their information and manuscript prior to its publication. We also thank anonymous reviewers for their helpful and constructive comments on the manuscript. This work comes under the British Antarctic Survey "BIOdiversity dynamics: Phylogeography, Evolution And Radiation of Life" (BIOPEARL) project.

Author information



Corresponding author

Correspondence to Chester J Sands.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

CJS performed all molecular procedures and analyses and prepared the initial manuscript. PC collected samples from Antarctic environs, contributed to the project design and the manuscript. KL initiated and led the project, contributed the project design, and contributed to the text. SJM collected samples from France, conducted all isolation and taxonomic identification and contributed to the manuscript.

Electronic supplementary material


Additional File 1: Sample information. Sample information provided includes taxonomic identification, location, sample type and GenBank accession numbers. (XLS 40 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Sands, C.J., Convey, P., Linse, K. et al. Assessing meiofaunal variation among individuals utilising morphological and molecular approaches: an example using the Tardigrada. BMC Ecol 8, 7 (2008).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Single Copy Nuclear Gene
  • Density Gradient Medium
  • Meiofaunal Group
  • Alpha Spectrin
  • Polymerase Chain Reaction Recombination