A needle in a haystack: A new metabarcoding approach to survey diversity at the species level of Arcellinida (Amoebozoa: Tubulinea).

Environmental DNA-based diversity studies have increased in popularity with the development of high throughput sequencing technologies. This permits the potential simultaneous retrieval of vast amounts of molecular data from many different organisms and species, thus contributing to a wide range of biological disciplines. Environmental DNA protocols designed for protists often focused on the highly conserved small subunit of the ribosome gene, that does not permit species-level assignments. On the other hand, eDNA protocols aiming at species-level assignments allow a fine level ecological resolution and reproducible results. These protocols are currently applied to organisms living in marine and shallow lotic freshwater ecosystems, often in a bioindication purpose. Therefore, in this study, we present a species-level eDNA protocol designed to explore diversity of Arcellinida (Amoebozoa: Tubulinea) testate amoebae taxa that is based on mitochondrial cytochrome oxidase subunit I (COI). These organisms are widespread in lentic water bodies and soil ecosystems. We applied this protocol to 42 samples from peatlands, estuaries and soil environments, recovering all the infraorders in Glutinoconcha (with COI data), except for Hyalospheniformes. Our results revealed an unsuspected diversity in morphologically homogeneous groups such as Cylindrothecina, Excentrostoma or Sphaerothecina. With this protocol we expect to revolutionize the design of modern distributional Arcellinida surveys. Our approach involves a rapid and cost-effective analysis of testate amoeba diversity living in contrasted ecosystems. Therefore, the order Arcellinida has the potential to be established as a model group for a wide range of theoretical and applied studies.