Revealing the Hyperdiverse Mite Fauna of Subarctic Canada through DNA Barcoding

Abstract

Although mites are one of the most abundant and diverse groups of arthropods, they are rarely targeted for detailed biodiversity surveys due to taxonomic constraints. We address this gap through DNA barcoding, evaluating acarine diversity at Churchill, Manitoba, a site on the tundra-taiga transition. Barcode analysis of 6279 specimens revealed nearly 900 presumptive species of mites with high species turnover between substrates and between forested and non-forested sites. Accumulation curves have not reached an asymptote for any of the three mite orders investigated, and estimates suggest that more than 1200 species of Acari occur at this locality. The coupling of DNA barcode results with taxonomic assignments revealed that Trombidiformes compose 49% of the fauna, a larger fraction than expected based on prior studies. This investigation demonstrates the efficacy of DNA barcoding in facilitating biodiversity assessments of hyperdiverse taxa.

About the Authors:

Monica R. Young

* E-mail: [email protected]

Affiliation: Biodiversity Institute of Ontario and Department of Integrative Biology, University of Guelph, Guelph, Ontario, Canada

Valerie M. Behan-Pelletier

Affiliation: Agriculture and Agri-Food Canada, Ottawa, Ontario, Canada

Paul D. N. Hebert

Affiliation: Biodiversity Institute of Ontario and Department of Integrative Biology, University of Guelph, Guelph, Ontario, Canada

Introduction

Species identification and discovery has been greatly accelerated by DNA barcoding, the analysis of sequence variation in a 648 base pair segment of the mitochondrial CO1 gene [1]. DNA barcoding has been successful in many animal groups [1]–[4], reflecting the fact that intraspecific sequence variation is consistently low, typically a fraction of a percent, while interspecific divergence usually exceeds 2%. When deep intraspecific variation is detected, cryptic species are often subsequently revealed through ecological or morphological study [5], [6].

The congruence in patterns of sequence variation across different taxonomic lineages allows the use of DNA barcodes to explore biodiversity in groups which lack a well-developed taxonomic framework. It facilitates rapid diversity assessment in such cases by enabling the delineation of MOTUs, molecular operational taxonomic units [7]. Because the quantification of biodiversity is transparent and reproducible, DNA barcoding is becoming a standard practice for assessing diversity patterns in poorly known taxa [4], [8], [9].

Although only 45,000 species have been described, Acari (mites) are believed to be one of the most diverse groups of arthropods, perhaps including more than 1 million species [10]. They are certainly one of the most abundant groups of arthropods as mite densities reach nearly 2 M individuals/m2 in temperate deciduous forest sites [11], nearly 0.5 M/m2 in dry tropical forests [11], and more than 0.1 M/m2 in northern sites [12]. Although they are often treated as members of the soil fauna, mites are associated with varied substrates [13] forming distinct assemblages on tree trunks, in soils, in surface litter, on fungi, and in aquatic habitats [10], [13]–[15].

The three dominant orders of soil mites have varied feeding modes, genetic systems and dispersal mechanisms [16]. The Sarcoptiformes are generally mycophagous or saprophagous feeders [16] with long adult lifespans [17] and a thelytokous parthenogenetic genetic system [18]. By contrast, the Mesostigmata tend to be free-living predators or parasites [16] with short adult lifespans and haplodiploid genetic systems [18]. Members of a third order, the Trombidiformes, show the greatest diversity in feeding mode (animal and plant parasites, free living predators, free living detritivores), and in genetic systems [16], [18].

Despite their diversity and abundance, mites are rarely included in biodiversity assessments because of serious taxonomic barriers. The status of many species is uncertain due to synonymies [19], morphotypes which are distinct species [20], and sexual dimorphisms [21]. Immature life stages are also excluded from surveys as they lack diagnostic morphological characters. Aside from these challenges, there is a scarcity of taxonomic experts. Consequently, surveys are often limited to higher level taxonomic assignments [22]–[25], or to assessments of a particular group [13], [26]–[28]. These factors preclude detailed assessments of the fauna, such as the examination of species turnover in space or time. DNA barcoding has the potential to radically advance our understanding of both the extent and patterns of species diversity in mites by providing a transparent, consistent method for delineating species which allows the inclusion of all life stages and both sexes. DNA barcoding has been successfully used in delimiting mite species [29], but prior work has focused on phylogenetic studies of a few species or genera [30]–[32].

Our work assesses the diversity of the mite fauna at one site in the Canadian subarctic and determines the extent of faunal divergence between major habitats. As such, it represents the first comprehensive assessment of mite diversity using molecular methods in any geographic setting.

Materials and Methods

Ethics Statement

No specific permits were required for the described field studies. Additionally, no specific permissions were required for these collection locations/activities as they are not privately owned or protected, and did not involve the collection of endangered or protected species.

Study Site/Sampling Design

Specimens were collected in the vicinity of Churchill during the snow-free season from June through August 2008 to 2011 by sweep netting, pitfall traps, Berlese funnel extractions and aspirators (Figure 1). A more regimented survey in 2010 included systematic sampling from 7 substrates at 10 locations over a 6 week period in boreal forest, bog, fen, tundra, marine beach, and rock bluff habitats. Seven substrates were sampled at each locale including moss, soil, litter, woody debris and lichens (Cladina spp., Peltigera leucophlebia, Parmelia/Hypogymnia). Approximately 500 mL of material from each substrate was collected, and extracted using modified Berlese funnels into 95% ethanol (EtOH). As well, two transects of five pitfall traps (with 95% EtOH) were deployed at each site and visited every 3 days for a total period of 9 days. Specimens were removed at each visit and placed into fresh 95% EtOH. Each of the pitfall transects and each substrate sample was treated as a separate analytical unit.

[Figure omitted. See PDF.]

Figure 1. Sampling locations at Churchill.

Maps depict a) the location of Churchill in Canada, and b) all sample locations along with specific sampling sites in the Churchill region.

https://doi.org/10.1371/journal.pone.0048755.g001

Sorting/Identifications

The specimens in each analytical unit were sorted into morphospecies, and 3–5 specimens of each were selected for sequence analysis. In total 8240 specimens (approximately 14% of the total catch) were selected for analysis. All specimens were identified to a family level using keys in Krantz and Walter [16], and sarcoptiform mites were identified to genus.

Barcoding Methodology

Each specimen was photographed and subsequently placed in a well containing 50 µl of 95% EtOH in a 96 well microplate. Collection details for each specimen together with its taxonomic assignment and its photograph are provided in a single data set on BOLD. The records can be retrieved using a DOI (http://dx.doi.org/10.5883/DATASET-MTBAR12N), a novel feature on BOLD enabling easy access and citability of barcode data [33]. Specimens were sequenced for the barcode region of the COI gene using standard protocols at the Canadian Centre for DNA Barcoding (CCDB) [34], using a cocktail of LepF1/LepRI [5] and LCO1490/HCO2198 [35] primers. Failed amplification reactions were further processed using the MLepF1 (Hebert unpublished) and MLepR2 (Prosser unpublished) primers. Glass fibre extraction was employed followed by voucher recovery [36]. DNA extracts were placed in archival storage at −80°C at the Biodiversity Institute of Ontario (BIO). Vouchered specimens were stored in 95% EtOH or slide mounted in Canada balsam, and deposited at both BIO and the Canadian National Collection of Insects, Arachnids and Nematodes.

Contigs were assembled and edited using CodonCode Aligner v. 3.0.1, and aligned by eye in MEGA 5.03 [37]. Each sequence with a length greater than 500 base pairs (bp) and with less than 1% ambiguous sites (Ns) was assigned a Barcode Index Number (BIN) by BOLD (Ratnasingham and Hebert 2012, in prep). All sequences with a length of 300 bp or longer, and with less than 1% Ns were also assigned to MOTUs using jMOTU [38] with a threshold of 15 nucleotide changes (2.3%), which is generally consistent with BIN assignments by BOLD. These MOTUs were used for further analysis, and are referred to interchangeably as BINs.

All sequence records together with trace files and images are available on BOLD as a single citable dataset (http://dx.doi.org/10.5883/DATASET-MTBAR12N). The sequences are also available on GenBank (Accessions GC680425–GU680432, GU680434–GU680497, GU702808–GU702809, HM405807–HM405810, HM405830–HM405857, HM431992–HM431993, HM431995–HM431998, HM904908, HM907069–HM907086, HM907088–HM907127, HM907130–HM907134, HM907138–HM907180, HM907182–HM907327, HM907329–HM907484, HQ558324–HQ558388, HQ558390–HQ558476, HQ558478–HQ558511, HQ558513–HQ558537, HQ558539–HQ558542, HQ558544–HQ558611, HQ558613–HQ558625, HQ558627, HQ558629–HQ558669, HQ558671–HQ558720, HQ558722–HQ558791, HQ941470–HQ941573, HQ941576–HQ941579, HQ966220–HQ966228, HQ966230–HQ966236, HQ966238–HQ966247, JX833624–JX838789).

Assessing Richness

We constructed specimen-based accumulation curves using random sampling and 1000 iterations for overall diversity, for each order, and for each family with more than 100 specimens, or with more than 10 BINs. This was done to assess diversity, and to determine which groups were undersampled. The slope of the accumulation curve for the last 10 specimens on the curve was calculated for each order, family, and others (families with fewer than 100 specimens or 10 BINs were pooled) to assess the completeness of sampling [39]. Clades with a slope >0.1 were viewed as very undersampled, while those with a slope >0.01 indicated lineages with modest undersampling. Predictions of total mite richness and richness of each order were also calculated using specimen-based Chao’s species richness estimator [40] using the vegan package in R [41], [42].

Faunal Similarity

Faunal similarity was assessed for samples collected systematically between the previously outlined sites and substrates. The similarity in community composition was visualized using cluster dendrograms computed from complete linkage hierarchial clustering method on Hellinger transformed abundances [43] and Bray Curtis dissimilarities using the vegan package in R [41], [42]. To test the significance of the clustering pattern between site type (forested or non-forested), we conducted an analysis of similarity (ANOSIM) with 999 permutations using the vegan package in R [41], [42].

Results

Assessing Richness

Barcode sequences were recovered from 6365 of the 8240 specimens, a success rate of 77.2% (Table 1). However, there was significant variation (x22 = 60.7, p<0.001) in recovery success among the three orders with a high of 80.4% for Sarcoptiformes and a low of 68.2% for Trombidiformes (Table 1). Most sequences (98.6%) were greater than 300 bp in length with less than 1% Ns. There were representatives of 899 BINs, with an average of 6.9 specimens per BIN. Mesostigmata was the least diverse order, accounting for 15% of the total diversity, while Sarcoptiformes composed 36% of the total diversity (Table 1). Trombidiformes were the most diverse order with 437 BINs (Table 1), accounting for 49% of the total diversity.

[Figure omitted. See PDF.]

Table 1. Observed and expected BIN richness for each order, as calculated by Chao’s estimator in R including error (±SE) estimates.

https://doi.org/10.1371/journal.pone.0048755.t001

Although the overall BIN accumulation curve did not reach an asymptote (Figure 2), there was sufficient data to estimate total mite diversity at Churchill as 1229 (±49) BINs (Table 1). The Sarcoptiformes was the best sampled order with a final accumulation curve slope of 0.022, and an estimated BIN richness of 423 (±27), while Mesostigmata were moderately well sampled with a final accumulation curve slope of 0.049, and an estimated richness of 173 (±16) BINs. The Trombidiformes was the least well-sampled order with a final accumulation curve slope of 0.093 and an estimated BIN richness of 633 (±38) (Table 1).

[Figure omitted. See PDF.]

Figure 2. Bin accumulation curves.

Curves represent a) the overall dataset, with 899 BINs from 6279 samples, and b) curves for each order - Trombidiformes (437 BINs, n = 1933), Mesostigmata (135 BINs, n = 849), and Sarcoptiformes (327 BINs, n = 3497).

https://doi.org/10.1371/journal.pone.0048755.g002

Among the five Mesostigmata families which met the requirements for detailed analysis, only the Zerconidae had a terminal slope of <0.01 (Table 2). Two families, Phytoseiidae and Melicharidae, had a slope >0.1 indicating that they include a very high number of uncollected species. Eleven of the 30 families of Sarcoptiformes met the requirements for analysis and 7 had a terminal slope >0.01 (Table 2). Four families (Ceratozetidae, Haplozetidae, Mycobatidae, Trhypochthoniidae) were well sampled, while two (Suctobelbidae, Brachychthoniidae) had terminal slopes ≥0.1 (Table 2). Ten of the 21 families of Trombidiformes met the requirements for analysis and all were under-sampled (>0.01 slope) with six families exceeding the 0.1 slope (Table 2). The Scutacaridae and Siteroptidae exhibited very unsaturated accumulation curves, with slopes of 0.50 and 0.37 respectively (Table 2).

[Figure omitted. See PDF.]

Table 2. Richness and terminal slope for the accumulation curve of selected families in three orders of Acari.

https://doi.org/10.1371/journal.pone.0048755.t002

The inclusion of qualitative samples from 2011 increased overall BIN richness by 30% (206 BINs), and increased the overall estimate of richness by 26% (253 BINs) (Table 3). Mite richness increased similarly among the orders, ranging from 28–36% (Table 3).

[Figure omitted. See PDF.]

Table 3. Effect of sampling patchy, short-lived habitats.

https://doi.org/10.1371/journal.pone.0048755.t003

Faunal Similarity

Mites showed high turnover between samples from different sites and substrates, with mean Bray-Curtis Dissimilarities ranging from 0.73 to 0.93 among the three orders. All three orders showed similar clustering patterns among sites, with distinct separation between forested and non-forested sites (ANOSIM Mesostigmata R = 0.396, p = 0.01; Sarcoptiformes R = 0.616, p = 0.004; Trombidiformes R = 0.620, p = 0.004) (Figure 3). When looking at faunal similarity between substrates, slightly different patterns were revealed. In all three orders the fauna from the arboreal lichens (Parmelia, Hypogymnia) was very distinct, as well as the fauna from woody debris (Figure 4). However, the trombidiform and mesostigmatan faunas from pitfall traps were also highly dissimilar to those from other substrates, while the Sarcoptiformes from pitfall traps was similar to the fauna found on forest floor lichens such as Peltigera and Cladonia (Figure 4).

[Figure omitted. See PDF.]

Figure 3. Cluster dendrograms showing the similarity of species assemblages for three mite orders among 10 sites in (F) and non-forested (N) settings.

https://doi.org/10.1371/journal.pone.0048755.g003

[Figure omitted. See PDF.]

Figure 4. Cluster dendrograms of showing the similarity of species assemblages for three mite orders among 8 substrates.

https://doi.org/10.1371/journal.pone.0048755.g004

Discussion

Sequencing Success

Varied primer binding and size differences between the major groups of mites may be responsible for variation in sequence recovery. Most species of Trombidiformes, the order with the lowest success, are very small so DNA concentrations may have been too low for successful PCR amplification. Some Mesostigmata are more heavily sclerotized than the other lineages, perhaps also reducing DNA recovery. Okassa et al. [44] reported problems in recovery of Cyt B sequences in their work on phytoseiid mites when DNA concentrations were less than 2.33 ng/ul. A shift to smaller elution volumes might improve success by producing higher DNA concentrations. Homogenizing specimens should aid DNA recovery [45], but it would lead to the destruction of specimens preventing their subsequent taxonomic study. Designing and utilizing taxa specific primers might also increase amplification success [46], but primer design requires prior taxonomic knowledge and affiliated reference sequences, both of which are usually unavailable for mites.

Assessing Richness

Our work has revealed the extreme diversity of the mite fauna at Churchill, despite our failure to collect and sequence all taxa. Sarcoptiform mites were better sampled than the other two orders as only 23% of the expected fauna remains undocumented. Lower sequencing success may at least partially account for lower completeness of species coverage for the other two orders. Sarcoptiform mites were the most abundant group in our samples, suggesting that our collection methods and scale of analysis were adequate to encounter most taxa. Alternatively, the Sarcoptiformes may include fewer rare species or show less local structure [47]. Mesostigmatan mites were moderately sampled, but 22% of the predicted fauna awaits collection. Their lesser coverage may be an artefact of their low abundance, a more patchy distribution, higher proportion of rare species [47], or lower success in sequencing. The accumulation curve of Trombidiformes is the least saturated, indicating 31% of trombidiform species await detection. With hyperdiverse groups increasing sample size will eliminate some of the singletons in the data but will invariably add new ones [48].

Similar trends in BIN richness and sampling saturation were also evident in the family accumulation curves (Figures 5, 6, 7) with most mesostigmatan families except the Zerconidae showing evidence of undersampling (Figure 5). Additionally, accumulation curves for all trombidiform families are unsaturated (Figure 7), implying general undersampling of this order. By contrast, all families of Sarcoptiformes except the Suctobelbidae and Brachychthoniidae showed a close approach to asymptotic diversity (Figure 6). However, differences in sampling saturation could be a result of differing sequencing success rates.

[Figure omitted. See PDF.]

Figure 5. BIN accumulation curves for 5 families of Mesostigmata.

https://doi.org/10.1371/journal.pone.0048755.g005

[Figure omitted. See PDF.]

Figure 6. BIN accumulation curves for 11 families of Sarcoptiformes.

https://doi.org/10.1371/journal.pone.0048755.g006

[Figure omitted. See PDF.]

Figure 7. BIN accumulation curves for 10 families of Trombidiformes.

https://doi.org/10.1371/journal.pone.0048755.g007

Despite incomplete sampling at Churchill, we found much greater acarine richness than recorded in past studies. Danks [49] reported 342 species of mites from the North American arctic, with 76 species of Mesostigmata, 144 Sarcoptiformes, and 122 Trombidiformes. We found nearly three times as many taxa at a single site, Churchill, with numbers increased by 78%, 127% and 258%, respectively from those reported by Danks [49]. Our sampling methods were likely to encounter a larger fraction of the species than the studies reviewed by Danks [49] which mainly assessed acarine diversity in soil cores. A few prior acarine surveys have generated species lists for entire countries [50], [51], but most are restricted to an order or family [52], [53]. Interestingly, the number of mite species which we detected using molecular methods at Churchill is close to the counts for New Zealand (1200 species [51]), and the UK (1700 species [50]).

No members of two other mite orders (Opiliocariformes, Holythrida) have been reported from the arctic [49] and they were absent from our samples. Their absence is unsurprising as both are small orders found in tropical/warm temperate climates. The Trombidiformes accounted for 49% of the species in our samples, a considerably higher level than 35% reported for the Canadian arctic [49], 23% for the Canadian subarctic [54], and 39% for the UK [50]. However, our results do conform with estimates from the high arctic (40% to 63% [54]), and with global acarine species descriptions as 48% are Trombidiformes [50]. Trombidiformes also compose 56% of the descriptions in North America, and 49% of the Australian fauna [50]. Trombidiformes have generally been thought to be a less important component of the acarine fauna in the subarctic, but our results challenge this conclusion. It is possible that insufficient efforts have been made to characterise the Trombidiformes of these regions, or that the resident species are more morphologically cryptic.

Sampling techniques can have an important impact on faunal discovery. Regimented sampling methods such as transects often overlook taxa that are rare or patchily distributed among sites. For example, a strict sampling technique revealed only 63% of the ant fauna known from La Selva [48]. More importantly, the species accumulation curve prematurely reached an asymptote, underestimating the true species richness [48]. Systematic sampling tends to capture dominant species, but often overlooks rare or transient species [48], [55]. By sampling patchy and temporary habitats, we dramatically increased the discovery of mite species at Churchill, most markedly in the Mesostigmata, reinforcing the notion that they are patchily distributed. The importance of microhabitats as potential refugia for rare species has been demonstrated for mesostigmatan mites, where most microhabitats contained only 2–3% of the collected species [56]. Behan-Pelletier [57] captured only one fifth of the fauna using systematic biodiversity sampling, while the rest of the fauna were uncovered by qualitative sampling of patchy habitats. This demonstrates that strict biodiversity surveys do not capture the complete fauna of a region, and emphasises the importance of qualitative sampling of ephemeral habitats to capture rare species.

Faunal Similarity

We found marked divergence in the mite faunas from forested and tundra settings. Prior studies have established that the species composition of mite communities can be influenced by vegetation type [58]. However, it is generally thought that the composition of soil mite communities is more strongly correlated with soil moisture, although vegetation type typically covaries with moisture [28], [58]. Rouse [59] showed that amount and seasonal patterns of soil moisture in Churchill are significantly different between forest and tundra sites, variation which may explain the distinctness of their mite communities. One exception to this pattern was the fauna of fens which are wet, but treeless habitats. The sarcoptiform fauna from the fen grouped more closely with the forested sites, whereas the trombidiform fauna was more similar to the tundra sites. The Trombidiformes may have been less impacted by the high soil moisture of the fen as they tend to be active surface predators [16], while the Sarcoptiformes are less mobile and are influenced by heterogeneity in soil [60].

Mite faunal similarity patterns between substrates were less obvious than those between sites. Our single arboreal substrate samples did not allow for a clear comparison between arboreal and forest floor substrates such as described by Lindo and Winchester [61]. However, arboreal lichens (Parmelia/Hypogymnia), woody debris, and pitfall samples formed the most distinct communities. The pitfall traps likely catch a functionally different fauna [55], [62], such as the fast active predators moving across the soil surface, as well as those that may be phoretic on other insects caught in the trap. Woody debris is slightly more ephemeral in nature, perhaps providing a unique microhabitat for transient species or more specialized species. On the other hand, the arboreal lichens potentially represent a distinctly different faunal community living in arboreal substrates [61], and represent a largely undocumented source of acarine diversity in Churchill.

Conclusion

Our study has revealed the power of DNA barcoding to provide insights into the diversity and distributional patterns of mites that could not have been gained through morphological approaches. Because of its use, we were able to analyze all life stages and both sexes, revealing a mite fauna with a species richness rivalling the most diverse of temperate habitats. Our work has also indicated the value of supplementing systematic sampling designs with qualitative sampling to ensure the examination of novel habitat types. The vouchered specimens generated through this study represent a valuable resource for future taxonomic research, particularly since specimens are partitioned into genetically cohesive assemblages. The utility of DNA barcoding for local biodiversity assessments is clear, but it will bring particular power to analyses which seek a deeper understanding of beta diversity patterns.

Acknowledgments

Alex Smith and Karl Cottenie provided advice on data analysis, while Christina Carr and Sarah Adamowicz aided with the molecular techniques. David Porco activated early work on mite barcoding at Churchill, while staff at the CCDB contributed importantly to specimen processing and sequence acquisition. Hans Klompen and other instructors at the Summer Acarology Program at The Ohio State University, Columbus, Ohio, aided taxonomic assignments.

Author Contributions

Conceived and designed the experiments: MRY PDNH. Performed the experiments: MRY. Analyzed the data: MRY. Contributed reagents/materials/analysis tools: PDNH. Wrote the paper: MRY PDNH VMBP. Contributed to specimen identification: VMBP MRY.

Citation: Young MR, Behan-Pelletier VM, Hebert PDN (2012) Revealing the Hyperdiverse Mite Fauna of Subarctic Canada through DNA Barcoding. PLoS ONE 7(11): e48755. https://doi.org/10.1371/journal.pone.0048755

References

1. Hebert PDN, Cywinska A, Ball SL, DeWaard JR (2003) Biological identifications through DNA barcodes. Proceedings of the Royal Society B: Biological Sciences 270: 313–322.

2. Kerr KCR, Stoeckle MY, Dove CJ, Weigt LA, Francis CM, et al. (2007) Comprehensive DNA barcode coverage of North American birds. Molecular Ecology Notes 7: 535–543.

3. Ward RD, Zemlak TS, Innes BH, Last PR, Hebert PDN (2005) DNA barcoding Australia’s fish species. Philosophical Transactions of the Royal Society B: Biological Sciences 360: 1847–1857.

4. Zhou X, Adamowicz SJ, Jacobus LM, DeWalt RE, Hebert PDN (2009) Towards a comprehensive barcode library for arctic life – Ephemeroptera, Plecoptera, and Trichoptera of Churchill, Manitoba, Canada. Frontiers in Zoology 6: 30 .

5. Hebert PDN, Penton EH, Burns JM, Janzen DH, Hallwachs W (2004) Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. Proceedings of the National Academy of Sciences of the United States of America 101: 12–17.

6. Burns JM, Janzen DH, Hajibabaei M, Hallwachs W, Hebert PDN (2008) DNA barcodes and cryptic species of skipper butterflies in the genus Perichares in Area de Conservación Guanacaste, Costa Rica. Proceedings of the National Academy of Sciences 105: 6350–6355.

7. Smith MA, Fisher BL, Hebert PDN (2005) DNA barcoding for effective biodiversity assessment of a hyperdiverse arthropod group: the ants of Madagascar. Proceedings of the Royal Society B: Biological Sciences 360: 1825–1834.

8. Zaldívar-Riverón A, Martínez JJ, Ceccarelli FS, De Jesús-Bonilla VS, Rodríguez-Pérez A, et al. (2010) DNA barcoding a highly diverse group of parasitoid wasps (Barconidae: Doryctinae) from a Mexican nature reserve. Mitochondrial DNA 21: 18–23.

9. Smith MA, Fernandez-Triana J, Roughley R, Hebert PDN (2009) DNA barcode accumulation curves for understudied taxa and areas. Molecular Ecology Resources 9: 208–216.

10. Walter DE, Proctor HC (1999) Mites: Ecology Evolution and Behaviour. New York: CABI Publishing. 352 p.

11. Peterson H, Luxton M (1982) A comparative analysis of soil fauna populations and their role in decomposition processes. Oikos 39: 288–388.

12. Hammer M (1952) Investigations on the microfauna of Northern Canada. Part I: Oribatida. Acta Arctica 4: 1–108.

13. Walter DE, Proctor HC (1998) Predatory mites in tropical Australia: Local species richness and complementarity. Biotropica 30: 72–81.

14. Proctor HC, Montgomery KM, Rosen KE, Kitching RL (2002) Are tree trunks habitats or highways? A comparison of oribatid mite assemblages from hoop-pine bark and litter. Australian Journal of Entomology 41: 294–299.

15. Wallwork JA (1972) Distribution patterns and population dynamics of the micro-arthropods of a desert soil in southern California. Journal of Animal Ecology 41: 291–310.

16. Krantz GW, Walter DE, editors (2009) A Manual of Acarology. Lubbock: Texas Tech University Press. 807 p.

17. Norton RA (1994) Chapter 5: Evolutionary aspects of Oribatid mite life histories and consequences for the origin of the Astigmata. In: MA Houck. Mites: Ecological and Evolutionary Analyses of Life-History Patterns. New York, Chapman and Hall. 99–135.

18. Norton RA, Kethley JB, Johnston DE, OConnor BM (1993) Phylogenetic perspectives on genetic systems and reproductive modes of mites. In DL Wrench, MA Ebbert. Evolution and Diversity of Sex Ratio in Insects and Mites. Chapman and Hall, New York.

19. Pfingstl T, Schäffer S, Krisper G (2010) Re-evaluation of the synonymy of Latovertex Mahunka, 1987 and Exochocepheus Woolley and Higgins, 1968 (Acari: Oribatida: Scutoverticidae). International Journal of Acarology 36: 327–342.

20. Laumann M, Norton RA, Weigmann G, Scheu S, Maraun M, et al. (2007) Speciation in the parthenogenetic oribatid mite genus Tectocepheus (Acari, Oribatida) as indicated by molecular phylogeny. Pedobiologia 51: 111–122.

21. Colloff MJ, Cameron SL (2009) Revision of the oribatid mite genus Austronothrus Hammer (Acari: Oribatida): sexual dimorphism and a re-evaluation of the phylogenetic relationships of the family Crotoniidae. Invertebrate Systematics 23: 87–110.

22. Filzek PDB, Spurgeon DJ, Broll G, Svendsen C, Hankard PK, et al. (2004) Metal effects on soil invertebrate feeding: Measurements using the bait lamina method. Ecotoxicology 13: 807–816.

23. Franklin E, Magnusson WE, Luizão FJ (2005) Relative effects of biotic and abiotic factors on the composition of soil invertebrate communities in an Amazonian savannah. Applied Soil Ecology 29: 259–273.

24. Sjursen H, Michelsen A, Holmstrup M (2005) Effects of freeze-thaw cycles on microarthropods and nutrient availability in a sub-arctic soil. Applied Soil Ecology 28: 79–93.

25. Sjursen Konestabo H, Michelsen A, Holmstrup M (2007) Responses of springtail and mite populations to prolonged periods of soil freeze-thaw cycles in a sub-arctic ecosystem. Applied Soil Ecology 36: 136–146.

26. Seyd EL (1981) Studies on the moss mites of Snowdonia (Acari: Oribatei). 2. The Cnicht. Biological Journal of the Linnean Society 15: 287–298.

27. Stamou GP, Sgardelis SP (1989) Seasonal distribution patterns of oribatid mites (Acari: Cryptostigmata) in a forest ecosystem. Journal of Animal Ecology 58: 893–904.

28. Minor M, Cianciolo J (2007) Diversity of soil mites (Acari: Oribatida, Mesostigmata) along a gradient of land use types in New York. Applied Soil Ecology 35: 140–153.

29. Hinomoto N, Tran DP, Pham AT, Ngockk Le TB, Tajima R, et al. (2007) Identification of spider mites (Acari: Tetranychidae) by DNA sequences: A case study in northern Vietnam. International Journal of Acarology 33: 53–60.

30. Dabert J, Ehrnsberger R, Dabert M (2008) Glaucalges tytonis sp. N. (Analgoidea, Xolalgidae) from the barn owl Tyto alba (Strigiformes, Tytonidae): compiling morphology with DNA barcode data for taxon description in mites (Acari). Zootaxa 1719: 41–52.

31. Ros VID, Breeuwer JAJ (2007) Spider mite (Acari: Tetranychidae) mitochondrial COI phylogeny reviewed: host plant relationships, phylogeography, reproductive parasites and barcoding. Experimental and Applied Acarology 42: 239–262.

32. Heethoff M, Laumann M, Weigmann G, Raspotnig G (2011) Integrative taxonomy: Combining morphological, molecular and chemical data for species delineation in the parthenogenetic Trhypochthonius tectorum complex (Acari, Oribatida, Trhypochthoniidae). Frontiers in Zoology 8: doi:10.1186/1742-9994-8-2.

33. Ratnasingham S, Hebert PDN (2007) BOLD: The Barcode of Life Data System (www.barcodinglife.org). Molecular Ecology Notes. 7: 355–364.

34. CCDB. The Canadian Centre for DNA Barcoding website. Available: www.ccdb.ca. Accessed 6 September 2012.

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

36. Porco D, Rougerie R, Deharveng L, Hebert P (2010) Coupling non-destructive DNA extraction and voucher retrieval for small soft-bodied Arthropods in a high-throughput context: the example of Collembola. Molecular Ecology Resources 10: 942–945.

37. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, et al. (2011) MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Molecular Biology and Evolution 28: 2731–2739.

38. Jones M, Ghoorah A, Blaxter M (2011) jMOTU and Taxonerator: Turning DNA barcode sequences into annotated operational taxonomic units. PLoS one 6: e19259 .

39. Hortal J, Lobo JM (2005) An ED-based protocol for optimal sampling of biodiversity. Biodiversity and Conservation 14: 2913–2947.

40. Chao A (1987) Estimating the population size for capture-recapture data with unequal catchability. Biometrics 43: 783–791.

41. Oksnen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, et al.. (2011) vegan: Community Ecology Package. R Package V. 2.0. Available: http://CRAN.R-project.org/package=vegan. Accessed 6 September 2012.

42. R Development Core Team (2011) R: A language and Environment for Statistical Computing. R Foundation for Statistical Computing. Vienna, Austria.

43. Rao R (1995) A review of canonical coordinates and an alternative to correspondence analysis using hellinger distance. QÜESTIIÓ 19: 23–63.

44. Okassa M, Kreiter S, Tixier M-S (2012) Obtaining molecular data for all life stages of Typhlodromus (Typhlodromus) exhilarates (Mesostigmata: Phytoseiidae): consequences for species identification. Experimental and Applied Acarology 57: 105–116.

45. Cruickshank RH (2002) Molecular markers for the phylogenetics of mites and ticks. Systematic and Applied Acarology 7: 3–4.

46. Dabert M, Witalinski W, Kasmierski A, Olszanowski Z, Dabert J (2010) Molecular phylogeny of acariform mites (Acari, Arachnida): strong conflict between phylogenetic signal and long-branch attraction artifacts. Molecular Phylogenetics and Evolution 56: 221–241.

47. Thompson GG, Withers PC (2003) Effect of species richness and relative abundance on the shape of the species accumulation curve. Austral Ecology 28: 355–360.

48. Longino JT, Colwell RK (1997) Biodiversity assessment using structured inventory: Capturing the ant fauna of a tropical rainforest. Ecological Applications 7: 1263–1277.

49. Danks HV (1981) Arctic Arthropods: A review of systematics and ecology with particular reference to the North American fauna. Ottawa: Entomological Society of Canada. 608 p.

50. Halliday RB, OConnor BM, Baker AS (2000) Global diversity of mites. In Raven PH, Williams T. Nature and Human Society: The Quest for a Sustainable World. Proceedings of the 1997 Forum on Biodiversity. Washington DC: National Academy Press. 192–212.

51. Zhang Z-Q, Rhode BE (2003) A faunistic summary of acarine diversity in New Zealand. Systematic and Applied Acarology 8: 75–84.

52. Maraun M, Schatz H, Scheu S (2007) Awesome or ordinary? Global diversity patterns of oribatid mites. Ecography 30: 209–216.

53. de Castro TMMG, de Moraes GJ (2010) Diversity of phytoseiid mites (Acari: Mesostigmata: Phytoseiidae) in the Atlantic Forest of São Paulo. Systematics and Biodiversity 8: 301–307.

54. Behan VM, Hill SB (1980) Distribution and diversity of North American Arctic soil acari. In D. Dindal. Soil Biology as Related to Land Use Practices. Proceedings of the VII International Soil Zoology Colloquium, July 29-Aug. 3, 1979. Syracuse: Office of Pesticide and Toxic Substances, Environmental Protection Agency. 717–749.

55. Longino JT, Coddington J, Colwell RK (2002) The ant fauna of a tropical rain forest: Estimating species richness three different ways. Ecology 83: 689–702.

56. Madej G, Barczyk G, Gawenda I (2011) Importance of microhabitats for preservation of species diversity, on the basis of mesostigmatid mites (Mesostigmata, Arachnida, Acari). Polish Journal of Environmental Studies 20: 961–968.

57. Behan-Pelletier VM (1998) Ceratozetoidea (Acari: Oribatida) of lowland tropical rainforest, La Selva, Costa Rica. Acarologia 39: 349–381.

58. Maclean SF, Douce GK, Morgan EA, Skeel MA (1977) Community organization in the soil invertebrates of Alaskan arctic tundra. Ecological Bulletins 25: 90–101.

59. Rouse WR (1984) Microclimate of arctic tree line 2. Soil microclimate of tundra and forest. Water Resources Research 40: 67–73.

60. Nielsen UN, Osler GHR, Campbell CD, Neilson R, Burslem DRFP, et al. (2010) The enigma of soil animal species diversity revisited: The role of small-scale heterogeneity. PLoS ONE 5: e11567 .

61. Lindo Z, Winchester NN (2009) Spatial and environmental factors contributing to patterns in arboreal and terrestrial oribatid mite diversity across spatial scales. Oecologia 160: 817–825.

62. Fisher BL (1999) Improving inventory efficiency: A case study of leaf-litter ant diversity in Madagascar. Ecological Applications 9: 714–731.