Development of two new sets of PCR primers for eDNA metabarcoding of brittle stars (Echinodermata, Ophiuroidea)

Abstract

Brittle stars (class Ophiuroidea) are marine invertebrates comprising approximately 2, 100 extant species, and are considered to constitute the most diverse taxon of the phylum Echinodermata. As a non-invasive method for monitoring biodiversity, we developed two new sets of PCR primers for metabarcoding environmental DNA (eDNA) from brittle stars. The new primer sets were designed to amplify 2 short regions of the mitochondrial 16S rRNA gene, comprising a conserved region (111–115 bp, 112 bp on average; named “16SOph1”) and a hyper-variable region (180–195 bp, 185 bp on average; named “16SOph2”) displaying interspecific variation. The performance of the primers was tested using eDNA obtained from two sources: a) rearing water of an 2.5 or 170 L aquarium tanks containing 15 brittle star species and b) from natural seawater collected around Misaki, the Pacific coast of central Japan, at depths ranging from shallow (2 m) to deep (> 200 m) sea. To build a reference library, we obtained 16S rRNA sequences of brittle star specimens collected from around Misaki and from similar depths in Japan, and sequences registered in International Nucleotide Sequence Database Collaboration. As a result of comparison of the obtained eDNA sequences with the reference library 37 (including cryptic species) and 26 brittle star species were detected with certain identities by 16SOph1 and 16SOph2 analyses, respectively. In shallow water, the number of species and reads other than the brittle stars detected with 16SOph1 was less than 10% of the total number. On the other hand, the number of brittle star species and reads detected with 16SOph2 was less than half of the total number, and the number of detected non-brittle star metazoan species ranged from 20 to 46 species across 6 to 8 phyla (only the reads at the “Tank” were less than 0.001%). The number of non-brittle star species and reads at 80 m was less than 10% with both of the primer sets. These findings suggest that 16SOph1 is specific to the brittle star and 16SOph2 is suitable for a variety of marine metazoans. It appears, however, that further optimization of primer sequences would still be necessary to avoid possible PCR dropouts from eDNA extracts. Moreover, a detailed elucidation of the brittle star fauna in the examined area, and the accurate identification of brittle star species in the current DNA databank is required.