Biodiversity Data Journal 7: e30970 'e @) doi: 10.3897/BDJ.7.e30970 open access Research Article Reference DNA barcodes and other mitochondrial markers for identifying Caribbean Octocorals Jaime G. Morint, Dagoberto E. Venera-Pontén$, Amy C. Driskelll, Juan A. Sanchez!, Howard R. Lasker*, Rachel Collin® + Laboratorio de Sistematica Molecular y Filogeografia, Facultad de Ciencias Bioldgicas, Universidad Nacional Mayor de San Marcos, Lima, Peru § Smithsonian Tropical Research Institute, Panama City, Panama | Laboratories of Analytical Biology, National Museum of Natural History, Smithsonian Institution, Washington, D.C., United States of America q Laboratorio de Biologia Molecular Marina - BIOMMAR, Bogota, Colombia # Department of Geology, University at Buffalo, Buffalo, United States of America o& Smithsonian Tropical Research Institute, Balboa, Panama Corresponding author: Rachel Collin (collinr@si.edu) Academic editor: Danwei Huang Received: 31 Oct 2018 | Accepted: 04 Feb 2019 | Published: 20 Feb 2019 Citation: Morin J, Venera-Ponton D, Driskell A, Sanchez J, Lasker H, Collin R (2019) Reference DNA barcodes and other mitochondrial markers for identifying Caribbean Octocorals. Biodiversity Data Journal 7: e30970. https://doi.org/10.3897/BDJ.7.e30970 Abstract DNA barcoding is a useful tool for documenting the diversity of metazoans. The most commonly used barcode markers, 16S and COl, are not considered suitable for species identification within some "basal" phyla of metazoans. Nevertheless metabarcoding studies of bulk mixed samples commonly use these markers and may obtain sequences for "basal" phyla. We sequenced mitochondrial DNA fragments of cytochrome oxidase c subunit | (COI), 16S ribosomal RNA (16S), NADH dehydrogenase subunits 2 (16S-ND2), 6 (ND6- ND3) and 4L (ND4L-MSH) for 27 species of Caribbean octocorals to create a reference barcode dataset and to compare the utility of COl and 16S to other markers more typically used for octocorals. The most common genera (Erythropodium, Ellisella, Briareum, Plexaurella, Muriceopsis and Pterogorgia) were effectively distinguished by small differences (5 or more substitutions or indels) in COl and 16S sequences. Gorgonia and Antillogorgia were effectively distinguished from each other by unique haplotypes, but the small genetic differences make distance approaches ineffective for these taxa. Plexaura, This is an open access article distributed under the terms of the CCO Public Domain Dedication. 2 Morin J et al Pseudoplexaura and Eunicea were indistinguishable from each other but were generally effectively distinguished from other genera, further supporting the idea that these genera have undergone a rapid endemic radiation in the Caribbean. Keywords cytochrome oxidase |, 16S, Panama, Bocas del Toro, gorgonian, DNA barcoding Introduction DNA barcoding is a useful tool for documenting the diversity of most metazoan groups (Bucklin et al. 2011). Short, easily amplifiable fragments that vary amongst closely related species are sequenced from specimens identified by experts and then used as a tool to identify sequences from unidentified samples (Hebert and Gregory 2005, Bucklin et al. 2011, Hajibabaei et al. 2007). This approach is useful in a variety of contexts (Bucklin et al. 2011, Bucklin et al. 2016, Fonseca et al. 2010, Hajibabaei et al. 2007); however the most commonly used barcode marker, a fragment of cytochrome c oxidase subunit | (COl hereafter), is not suitable for species identification within some metazoan groups such as Porifera and Anthozoa (Shearer et al. 2002, Huang et al. 2008). Consequently, less effort has been invested in building reference databases of COI and 16S barcodes within these taxa, as opposed to those groups where these markers are perceived as useful for species identification. DNA barcoding has been used extensively to detect cryptic diversity within clades, with specific studies generally focusing on diversity within a single family, class or order. When COI is problematic, alternate barcodes are used, for example, 16S is commonly used as the molecular barcode in hydrozoans (Moura et al. 2011Zheng et al. 2014). For some clades, significant effort has been spent on identifying reliable, divergent markers and building extensive databases of these sequences for identification purposes (e.g. van der Ham et al. 2009, McFadden et al. 2014). Recent assessments of biodiversity are increasingly focusing on metabarcoding analyses of bulk mixed samples, such as gut contents (Leray et al. 2013, Oliverio et al. 2009), settlement plates (Zaiko et al. 2016), plankton samples (Bucklin et al. 2016), sediment samples (Fonseca et al. 2010) or environmental DNA from water samples (Stat et al. 2017, Thomsen et al. 2012). These samples are typically analysed using markers that amplify across the most diverse set of taxa; most commonly COI and sometimes 16S for studies of species diversity or the nuclear SSU 18S gene in studies aimed at higher taxonomic levels. Due to the non-selective nature of metabarcoding, this approach generates sequences for taxa where the marker is not useful for species-level identification as well as for those where it is useful. Therefore, a reference set of COI and 16S sequences for these taxa could at least provide the lowest possible identification when these problematic taxa occur in mixed samples. Reference DNA barcodes and other mitochondrial markers for identifying ... 3 Most Caribbean gorgonian octocorals are endemic and closely related (Sanchez et al. 2003, Sanchez 2016, Wirshing et al. 2005). This results in many taxa with morphologically similar sister species making it extremely difficult to differentiate a number of the species with morphological characters (Sanchez and Wishing 2005). In addition, octocorals are notoriously difficult to distinguish with molecular data (France and Hoover 2002, McFadden et al. 2010, McFadden et al. 2014, McFadden et al. 2004, van der Ham et al. 2009) since their mitochondria have a unique DNA repair mechanism that is thought to be responsible for the slow evolution of mitochondrial genes (Bilewitch and Degnan 2011). Nevertheless, gorgonians are an abundant and diverse part of the macro-fauna of Caribbean reefs, with as many as 40 species occurring on a single Caribbean reef (Lasker and Coffroth 1983, Sanchez et al. 1998, Sanchez et al. 1997, Etnoyer et al. 2010, Tsounis et al. 2018). In some places, they contribute significantly to three dimensional reef structures and provide microhabitats for macro-invertebrates and fishes (Sanchez 2016). Furthermore, octocoral abundances have increased at some locations in contrast to the overall decline in scleractinians (Edmunds and Lasker 2016, Ruzicka et al. 2013). Therefore, tools that facilitate their inclusion in biodiversity assays would be a useful addition to researchers’ toolkits. Here we had two main objectives: (1) to generate a reference set of DNA barcode sequences for common octocorals from the Southern Caribbean and (2) determine how the 16S and COI barcode fragments compare with other mitochondrial markers in their ability to distinguish genera and, in some cases, species in the fauna of Bocas del Toro, Panama. Materials and Methods A total of 180 octocoral tissue samples were collected by SCUBA-diving during the summer of 2007 from the shallow waters around the Bocas del Toro Archipelago on the Caribbean coast of Panama. Tissue samples consisted of 10-20 cm sections clipped from a distal branch for each colony. Samples were identified to species using a combination of visual identification in the field (colony morphology) and microscopic examination of spicule preparations. In some cases, individuals were intermediate in morphology or could not be identified to species and were only identified to genus. Dry tissue vouchers were deposited at the USNM and the University of Panama. Details about individual samples are provided in the BoLD project workbench (Ratnasingham and Herbert 2007). Genomic DNA was extracted from small samples of tissue (0.5 cm’) from each specimen using the mouse tail kit on a Biosprint 96 (Qiagen). The resuspension volume was 200 ul. We used PCR to amplify fragments of the mitochondrial genes COI, 16S, ND6-ND3, 16S- ND2 and ND4L-MSH with the primers listed in Table 1. When the barcode (“Folmer”) fragment of COI could not be amplified with the primers derived from Folmer et al. 1994, we amplified a larger fragment that includes part of the COIl gene. We then extracted the barcode sequence from this longer sequence. When neither of these amplifications was successful, we amplified a shorter COI fragment with gorgonian-specific primers. The 10 ul PCR cocktail included 0.5 ul 50 mM MgClo, 0.1 ul 20 ug/ul BSA, 0.5 ul dNTPs (2.5 mM 4 Morin J et al each), 0.3 yl each 10 mM primer and Biolase taq polymerase (Bioline). Annealing temperatures were 46°C for the barcode fragment of COI and ND4L-MSH, 48°C for 16S and the gorgonian-specific COI fragment and 52°C for ND6-ND3, 16S-ND2 and the long fragment of COI. Amplified PCR products were cleaned using the ExoSAP-IT protocol (Affymetrix) and sequenced using the Big Dye protocol (Thermofisher). Table 1. Primers used for amplification and sequencing of loci in this study and size of amplified fragments. *Lengths of fragments that include non-coding regions as well as ribosomal sequences can vary. | McFadden et al. 2004; @France and Hoover 2001; Meyer 2003; “this study; °Palumbi et al. 1991; © Brugler and France 2008; ’Sanchez et al. 2003. Locus Primers Expected Size* COl long COII-8068F'’ CCATAACAGGACTAGCAGCATC 940 COIOCTR? ATCATAGCATAGACCATACC COI Folmer dglLCO1490? GGTCAACAAATCATAAAGAYATYGG 655 dgHCO21987 TAAACTTCAGGGTGACCAAARAAYCA COl short GorgM13_F* CACGACGTTGTAAAACGACGTATGT TAGGAGATGATCATCTATAT 468 GorgM13_R* GGATAACAATTTCACACAGGGAATGTTGTATTAAAATTYCTRTCTGT 16S 16Sar° C@CCTGTTTATCAAAAACAT 668 16Sbr° CCGGTCTGAACTCAGATCACGT ND6-ND3_—- NDG-1487F' TTTGGTTAGTTATTGCCTTT ae ND3-2126R' CACATTCATAGACCGACACTT 16S-ND2 = 16S-647F' TTTGGTTAGTTATTGCCTTT 758 ND2-1418R' ACATCGGGAGCCCACATA ND4L-MSH ND4L-2475F® TAGTTTTACTGGCCTCTAC 870 MUT-3458R’ TSGAGCAAAAGCCACTCC Sequences were screened for quality and contigs of forward and reverse sequences were produced using Sequencher 5.4.6 (Gene Codes). Only sequences with a Phred quality score of at least 30 for more than 85% of the bases were combined into contigs and used for analysis. Sequences were compared internally across our dataset and BLASTned against GenBank sequences to check for contamination or mislabelling. After this step, through an abundance of caution, all 27 ND4L-MSH sequences from a single plate were excluded from the analyses because a subset of them showed unexplained divergences from previously published GenBank sequences. No other plate of sequences in our analysis showed similar problems. Sequences of each marker were aligned with ClustalX (gap-opening penalty: 15, gap-extension penalty: 6.66, transition weight: 0.5, delay divergent cutoff: 30%) and used to generate a matrix of pairwise differences [including the number of substitutions and indels (insertions/deletions)] between all the sequenced specimens. This matrix was then used to build heatmaps displaying the average pairwise differences between species. The patterns observed in heatmaps were contrasted with those observed in maximum parsimony and neighbour-joining trees (BIONJ, Gascuel 1997) built with the same matrix of pairwise differences. Reference DNA barcodes and other mitochondrial markers for identifying ... 5 Most markers showed two groups within which the pairwise differences amongst species (and genera) were consistently small. For these groups, we constructed haplotype networks for all the markers using Haplotype Viewer (Center for Integrative Bioinformatics Vienna; hito:/www.cibiv.at/) to determine if haplotypes were unique or related to particular species or genera. Data Resources The data underpinning the analysis reported in this paper are deposited in the Barcode of Life Database (dx.doi.org/10.5883/DS-OCTOCORA) (Ratnasingham and Herbert 2007) and the barcode markers have the following GenBank numbers COI:MK153303-MK153482 and 16S:MK153483-MK153626. Results and Discussion We were able to collect and successfully sequence 28 species represented by 180 individuals of octocorals from Bocas del Toro (Table 2). This represents 45% of the 61 species and includes 12 of the sixteen genera and all seven families reported for the region (Sanchez and Wishing 2005). Most markers were sequenced successfully for most species but 16S-ND2 could not be amplified for the genus Erythropodium, whereas the ND4L-MSH fragment failed to amplify in 5 species. Likewise, samples of Plerogorgia anceps almost entirely failed to amplify. Since all the specimens were processed with the same methods, differences in sequencing success may be due to taxon-specific secondary metabolites interfering with PCR. Table 2. Number of individuals collected and successfully sequenced for each genetic marker for samples identified to species. The colour of the cells indicates the availability of sequences for the same marker and species in GenBank. Dark grey: A new contribution; no available sequences overlap our fragment by 2100 bp. Light grey: Partial sequence is already available; our sequences overlap existing sequences by >350 bp; our data extends this by >50 bp. White: Sequence already available; sequences available in GenBank overlap ours completely with <50bp additional data. X mark indicates species that failed to amplify or resulted in unusable or suspicious sequences. “Sequences were also generated for an additional 4 Antillogrogia, 3 Erythopodium, 44 Eunicea, 3 Muricea, 1 Plexaura, 5 Plexaurella, 16 Pseudoplexaura and 9 Pterogrogia that were not identified to species. Species Individuals collected* COl 16S 16S-ND2 ND6-ND3 ND4L-MSH Antillogorgia acerosa 4 3 3 3 3 3 Antillogorgia americana 4 4 4 4 4 3 Antillogorgia bipinnata 8 5 5 5 5 4 Antillogorgia rigida 2 1 1 1 1 X 6 Morin J et al Species Individuals collected* COl 16S 16S-ND2 ND6-ND3 ND4L-MSH Briareum asbestinum S) 8 7 7 3 7 Ellisella schmitti 2 2 2 1 2 2 Erythropodium caribaeorum =. 3 3 1 x 3 x Eunicea calyculata 3 3 2 3 3 3 Eunicea clavigera 1 1 1 1 j 1 Eunicea flexuosa 9 8 8 8 8 8 Eunicea fusca 5 4 4 4 4 4 Eunicea laxispica 1 1 1 1 1 1 Eunicea succinea 1 1 1 1 1 X Eunicea tourneforti 6 6 4 6 6 3 Gorgonia mariae 8 8 8 8 8 X Gorgonia ventalina 4 4 4 4 4 Xx Muricea laxa 7 7 3 6 7 6 Muricea muricata 3 3 3 3 3 3 Muricea pinnata 5 5 3 5 S 5 Muriceopsis bayeriana 6 6 6 6 6 X Muriceopsis flavida 4 3 4 3 3 Xx Plexaura homomalla 14 8 6 8 8 ¥ Plexaura kuna 4 4 4 4 4 2 Plexaurella dichotoma 3 3 3 3 3 3 Plexaurella nutans 4 2 2 2 2 2 Pseudoplexaura porosa 2 2 1 2 2 2 Pseudoplexaura wagenaari 1 1 Xx 1 xX 1 Pterogorgia anceps 4 x X 1 x X Pterogorgia citrina 2 2 2 2 2 x New additions to GenBank from this study include COI for 19 species, 16S for 20 species, 16S-ND2 and ND4L-MSH for 3 species and 16S-ND2 and ND6-ND3 for 13 species (Table 2). Of those markers that were already in GenBank for a particular species, >40 of our sequences extended significantly beyond the previously reported fragment. Our samples also increased the geographic coverage of these species. For example, the majority of 16S-ND2 and ND4L-MSH sequences already in GenBank were collected from the Bahamas (Sanchez et al. 2003), so the addition of sequences from Panama provides new references for other areas of the species ranges. As expected for octocorals, none of our 5 mitochondrial markers was suitable for distinguishing amongst congeneric species across the entire dataset, but they were generally useful to distinguish amongst most genera except those in the Plexaura- Reference DNA barcodes and other mitochondrial markers for identifying ... 7 Pseudoplexaura-Eunicea group and the Antillogorgia-Gorgonia group (Figs 1, 2, 3, 4Fig. 5). The studied markers were able to distinguish species, but their patterns of pairwise differences varied amongst the different markers. COI, 16S-ND2 and ND4L-MSH displayed higher divergences overall than did the other two markers. The genera Briareum, Plexaurella, Erythropodium, Muriceopsis and Pterogorgia were clearly distinguishable from each other. Other taxa were partially distinguishable. For example, species of Plexaura, Pseudoplexaura and Eunicea were often indistinguishable from each other but could generally be distinguished from species in other genera across the dataset. The same situation occurred with Antillogorgia and Gorgonia which had very similar haplotypes. Ellisella and Erythropodium were clearly distinguishable with COI, 16S and 16S-ND2 but surprisingly appeared to fall within the Plexaura-Pseudoplexaura-Eunicea clade with ND4L- MSH (Ellisella) and ND6-ND3 (both). We could find no problems with our sequences or data flow which would account for this, but this result should be independently confirmed for ND4L-MSH and ND6-ND3 in these two genera. > rrrBPerrrtrrrrermerrreoanaz2zazza2aztsasaTsT7 SF z i Fs 85 2 g COI = ZS se & = 2 5 = 2 » 5 Fs 822 £ = = 2 = = 5 2 8 5 = Figure 1. EES Heatmap of the mean pairwise differences between species of Caribbean octocorals, based on nucleotide substitutions in their COl sequences. Colours on the diagonal indicate the mean intra-species differences. Species with a dash in the diagonal were represented by only one individual and thus their intra-species differences could not be calculated. Haplotype network analyses of the Plexaura-Pseudoplexaura-Eunicea group and the Antillogorgia and Gorgonia group, were able to distinguish some genera and species which were not easily distinguished by distance analysis as visualised in the heatmaps. All of the markers except for 16S show that Gorgonia and Antillogorgia can be separated, as each is comprised of unique haplotypes Fig. 6, although the haplotypes generally differ by only one or a few base pairs. In addition, COI appeared to distinguish the different species within each genus. For example, the COI network clearly separated G. mariae from G. ventalina but did not clearly separate the species of Antillogorgia from each other. However, our small sample sizes are not a strong test of the efficacy of each marker to uniquely identify 8 Morin J et al each species. Larger samples from the across the Caribbean should be sequenced before these haplotypes could be considered diagnostic. rr rPrReanananananmnnnmnroaoeg zs Zs =P? BBE ER BSG ft eee ee ece s g BE se ea ZBEiEi2 = oe 8 2 8.8 oe oo e228 EBEEE s See ere a Ae Be oe. = 8 & 6 Be Ee E 8388 88 388 § Fes 3 sss tat le OR fete a Pes Sli. FB - 5 Fe FE BE es s & a 8 ¢ eEsgsea#s [eg ear Ses s B & 4823 2 FA a's & chal PRIN IC ees ga = FO — Aa Se a Bo. 3 Ss Zz =e-es & 16S zee s 282° 237 8 FB @¢ eee & 88s 8 2 ER Se Fs 5 sake gR5 EE = = a = » $* 5 8 ae 5 2 5 & @ 6 3 = 8 | ¢ Ee = s Zz 5 & = 3 & = é c 5 $ FI = = = Col Bp 1s a a 0 15 20 30 Erythropodium_caribaeorum 40 Eun eulata 50 Plexaurella_nutans Pseudoplexaura_porosa Pterogorgia_citrina Figure 2. EES Heatmap of the mean pairwise differences between species of Caribbean octocorals, based on nucleotide (substitutions) and indel (insertions/deletions/gaps) differences of their 16S sequences. Colours on the diagonal indicate the mean intra-species differences. Species with a dash in the diagonal were represented by only one individual and thus their intra-species differences could not be calculated. i] wuy qnuy way Wwuy id Ss 53538 rre 2: 15. 3. xe] EEE E & BE ss s & &® 8 5 Vid eases, | iE 2S) 5s = Ss — » 2 e's ‘> |x 3% 2 FPst* saga 8 F - a ar = g 8 | s 8h 8 e585 Fe fad s @ £ $5 6 #3 » $322 & = S 2828 g 5 » 2 3% Be s gs g * 3 * 2 $ Fs 5 = = Figure 3. EE Heatmap of the mean pairwise differences between species of Caribbean octocorals, based on nucleotide (substitutions) and indel (insertions/deletions/gaps) differences of their 16S-ND2 sequences. Colours on the diagonal indicate the mean intra-species differences. Species with a dash in the diagonal were represented by only one individual and thus their intra-species differences could not be calculated. Reference DNA barcodes and other mitochondrial markers for identifying ... PREREERES PERS ERY SRS SSSR ier riz 2 &. ¢s q B2BBEBE FHS SEB SE EBS SS SEBS SEESEEEBERCKEERE| SB EEEE Sr FFF FF FFE BEARER EEE EE S s¢$s 8 § =32e8e88:88e:8882882888 8888 & 8B E EF EE SR a a a a ee ee Ss ae ls — Se -* Ss 5 6 5 EL = Ss eeéaes = = = J s £ 5 Sess 8 $$ R22 & 5? » 3 = a 2 E 6 * & 8 2yrs *® & a. i. S 5 = . s = 2 8 8 ma 3 5 Fe $ 5 ¢ rH a 4 z 8 g 3 5 gS Ss & 5 3 = Antillogorgia_acerosa Antillogorgia_americana Antillogorgia_bipinnata Antillogorgia_rigida Briareum_asbestinum Ellisella_schmitti Erythropodium_caribaeorum Eunicea_calyculata Eunicea_calyculata_coronata Eunicea_clavigera Eunicea_fexuosa Eunicea_fusea Eunicea_laxispica Eunicea_succinea Eunicea_tourneforti Gorgonia_mariae Gorgonia_ventalina Muricea_laxa Muricea_muricata Muricea_pinnata Muriceopsis_bayeriana Muriceopsis_flavida Plexaura_homomalla Plexaura_kuna Plexaurella_dichotoma Plexaurella_nutans Pseudoplexaura_porosa Pterogorgia_citrina Figure 4. ERY Heatmap of the mean pairwise differences between species of Caribbean octocorals, based on nucleotide (substitutions) and indel (insertions/deletions/gaps) differences of their ND6- ND3 sequences. Colours on the diagonal indicate the mean intra-species differences. Species with a dash in the diagonal were represented by only one individual and thus their intra- species differences could not be calculated. rPrrraearaenmnanrnnrnnrnnmpezazz Zzawttztwet7s7s Fy, Boce. (Bo Sete 8. 8 eS ce ce ee ce ee BEEeaBPERERRBRS RRB RBER SER E EB EE Ee & E & 5 = oh oo of *ee8f Ee 8 88 8&8 8 & & Sk FS & EEE a & Seeger eee eee se esse gee = = 2 ‘ee 2 eos oot le aa | = =e & moe, & & 8 F FES = 8 by. > pr F BP KR wz 6 & 2 Bs g 4 &©€ ® ES Ss Erte & s Ss & 52327286 5&5 5B @ 8 & & ‘a ee Ss 8B e668 & » 6 3 228 Fes § Be Ee FEES 8 = © @ 8 ¢ os 8 38 6 BSaS FF 3 Bw gS & Ss 8 = 8 I - S y= ee — Ss S& 8 > = s gs = 5 © = 5 ao 5 s s = ££ 4 i Tl ea i & = s SE =» $s 8 2 = 3 6 8 5s ££ $ * = 8 ) é = 5 Es x a. = 7 Antillogorgia_acerosa Antillogorgia_americana Antillogorgia_bipinnata Briareum_asbestinum Ellisella_schmitti Eunicea_calyculata Eunicea_calyculata_coronata Si- Eunicea_clavigera 5.0 5.0|_3.4 S$ 40 25 25 3.5| 0.7) Eunicea_flexuosa Eunicea_fusca 3.5 Eunicea_laxispica 20 20 30 25 3.5 0.0 0.0 5.0 25 2: Eunicea_tourneforti Muricea_laxa Muricea_muricata Muricea_pinnata Plexaura_homomalla Plexaura_kuna Plexaurella_dichotoma Plexaurella_nutans Pseudoplexaura_porosa Pseudoplexaura_wagenaari Figure 5. EE Heatmap of the mean pairwise differences between species of Caribbean octocorals, based on nucleotide (substitutions) and indel (insertions/deletions/gaps) differences of their ND4L- MSH sequences. Colours on the diagonal indicate the mean intra-species differences. Species with a dash in the diagonal were represented by only one individual and thus their intra-species differences could not be calculated. 10 Morin J et al 16S-ND2 ND6-ND3 o @ r=) a e@e@3—_— Antillogorgia americana Antillogorgia acerosa Gorgonia mariae _ Antillogorgia bipinnata Gorgonia ventalina Antillogorgia rigida Figure 6. EE Haplotype networks involving Gorgonia and Antillogorgia DNA sequences. Each haplotype is represented by a circle. The number of individuals within each haplotype is indicated within the corresponding circle. The circle's sizes are proportional to the number of individuals. The white areas within some circles represent individuals identified only to the genus level (species unknown), whereas other colours are species-specific. The length between haplotypes represents the number of nucleotide differences between them: 1 length unit = 1 substitution or indel (insertion/deletion). In contrast, network analysis of the Plexaura-Pseudoplexaura-Eunicea group further reinforces the results of the heatmaps, indicating that these genera cannot be distinguished using any of the markers Fig. 7. For all of the makers, there were single haplotypes that were shared by Eunicea, Plexaura and Pseudoplexaura. Muricea was not clearly distinguished from the other genera using COI and 16S with the pairwise differences, but species in this genus were consistently distinguished with the network approach. These results further support the idea that the Plexaura-Pseudoplexaura-Eunicea group has undergone a rapid endemic radiation in the Caribbean making them particularly difficult to distinguish (Sanchez et al. 2003, Sanchez and Wishing 2005, Sanchez 2016). Likewise, both Antillogorgia and Eunicea include exemplary cases of depth-induced ecological speciation and recent divergence (Prada and Hellberg 2013, Calixto-Botia and Sanchez 2017). As has been previously demonstrated for octocorals from other regions (Baco and Cairns 2012,France and Hoover 2002, France and Hoover 2001, McFadden et al. 2014, McFadden et al. 2004, McFadden et al. 2010), we found limited ability of the 5 markers to distinguish species based on distance methods typically employed by barcode studies Figs 1, 2, 3, 4Fig. 5. Nevertheless, metabarcoding studies are liable to obtain COI and 16S sequences from octocorals and researchers are likely to search for species identifications in BoLD or GenBank. Our sequences are the first COl and 16S sequences for many Caribbean octocoral species in GenBank and our results demonstrate that the COl barcode sequence can distinguish Caribbean species as effectively as other mitochondrial Reference DNA barcodes and other mitochondrial markers for identifying ... 11 sequences more commonly used in octocoral-focused studies. Despite this, the BoLD database assigned our 28 species to only 9 Barcode Index Numbers (BINs; Ratnasingham and Hebert 2013) highlighting that alternate approaches, such as haplotype network analysis or character-based DNA barcoding (Rach et al. 2008), designed for small divergences, are more appropriate for attempts to use DNA sequences to automatically identify octocorals (McFadden et al. 2010Baco and Cairns 2012). COI 16S 2 1 r 2 IGN “4 . 20 a ae wy 16S-ND2 ND6-ND3 4 » @ * | hm iy 1 Eunicea Plexaura Pseudoplexaura Muricea Ellisella Erythropodium Figure 7. EE Haplotype networks of DNA sequences belonging to genera that could not be clearly separated by their nucleotide differences (heatmaps). For each gene, only the genera that were not clearly different in the heatmaps were included. Each haplotype is represented by a circle and each genus is represented by a unique colour. The number of individuals belonging to each haplotype are indicated within the corresponding circle. The circle's sizes are proportional to the number of individuals. The length between haplotypes represents the number of nucleotide differences (Substitutions or indels) between them. Haplotypes present in multiple species, all belonging to the same genus, are indicated with a black star and haplotypes with all individuals belonging to the same species are indicated with underlined numbers. To be conservative, where samples only identified to genus were included in a circle, we did not consider them conspecific with samples identified to species, although they could have been. Acknowledgements We thank the Autoridad de los Recursos Acuaticos de Panama for issuing permits for this work, the participants in the 2007 workshop on Gorgonian Biology and Systematics in Bocas del Toro and C. McFadden whose suggestions significantly improved this manuscript. This research was supported by a Grand Challenges Level 1 grant from the Smithsonian Institution. All or portions of the laboratory and computer work were conducted in and with the support of the L.A.B. facilities of the National Museum of Natural History. 12 Morin J et al Author contributions JM- contributed to the analysis and manuscript preparation; RC- contributed to the analysis and manuscript preparation; DVP-— contributed to the analysis and manuscript preparation; ACD-— generated the sequences and contributed to the analysis; JAS- Collected and identified the samples; HRL- Collected and identified the samples References ° Baco AR, Cairns SD (2012) Comparing molecular variation to morphological species designations in the deep-sea coral Narella reveals new insights into seamount coral ranges. PLoS One 7 (9): 1-15. https://doi.org/doi.org/10.1371/journal.pone.0045555 ° Bilewitch JP, Degnan SM (2011) A unique horizontal gene transfer event has provided the octocoral mitochondrial genome with an active mismatch repair gene that has potential for an unusual self-contained function. BMC Evolutionary Biology 11 (1): 228. https://doi.org/10.1186/1471-2148-1 1-228 ° Brugler MR, France SC (2008) The mitochondrial genome of a deep-sea bamboo coral (Cnidaria, Anthozoa, Octocorallia, Isididae): genome structure and putative origins of replication are not conserved among octocorals. Journal of Molecular Evolution 67 (2): 125-136. https://doi.org/10.1007/s00239-008-9116-2 ° Bucklin A, Steinke D, Blanco-Bercial L (2011) DNA barcoding of marine Metazoa . Annual Review of Marine Science 3 (1): 471-508. https://doi.org/10.1146/annurev- marine-120308-080950 ° Bucklin A, Lindeque PK, Rodriguez-Ezpeleta N, Albaina A, Lehtiniemi M (2016) Metabarcoding of marine zooplankton: prospects, progress and pitfalls. Journal of Plankton Research 38 (3): 393-400. https://doi.org/10.1093/plankt/fow023 ° Calixto-Botia |, Sanchez JA (2017) A case of modular phenotypic plasticity in the depth gradient for the gorgonian coral Antillogorgia bipinnata (Cnidaria: Octocorallia). BMC Evolutionary Biology 17 (1): 55. https://doi.org/10.1186/s12862-017-0900-8 ° Edmunds PJ, Lasker HR (2016) Cryptic regime shift in benthic community structure on shallow reefs in St. John, US Virgin Islands. Marine Ecology Progress Series 559: 1-12. https://doi.org/10.3354/meps11900 ° Etnoyer P, Wirshing H, Sanchez J (2010) Rapid assessment of octocoral diversity and abitat on Saba Bank, Netherlands Antilles. PLoS One 5 (5): e10668. https:// doi.org/10.1371/journal.pone.0010668 ° Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R (1994) DNA primers for amplification of mitochondrial cytochrome c oxidase subunit | from diverse metazoan invertebrates. Molecular marine biology and biotechnology 3 (5): 294-299. ° Fonseca V, Carvalho G, Sung W, Johnson H, Power D, Neill S, Packer M, Blaxter M, D. Lambshead PJ, Thomas WK, Creer S (2010) Second-generation environmental sequencing unmasks marine metazoan biodiversity. Nature Communications 1 (7): 98. https://doi.org/10.1038/ncomms1095 ° France S, Hoover L (2002) DNA sequences of the mitochondrial COI gene have low levels of divergence among deep-sea octocorals (Cnidaria: Anthozoa). Hydrobiologia 471: 149-155. https://doi.org/10.1023/A:1016517724749 Reference DNA barcodes and other mitochondrial markers for identifying ... 13 France SC, Hoover LL (2001) Analysis of variation in mitochondrial DNA sequences (ND3, ND4L, MSH) among Octocorallia (=Alcyonaria) (Cnidaria: Anthozoa). Bulletin of the Biological Society of Washington 10: 110-118. Gascuel O (1997) BIONJ: an improved version of the NJ algorithm based on a simple model of sequence data. Molecular Biology and Evolution 14 (7): 685-695. https:// doi.org/10.1093/oxfordjournals.molbev.a025808 Hajibabaei M, Singer GA, Hebert PD, Hickey DA (2007) DNA barcoding: how it complements taxonomy, molecular phylogenetics and population genetics. Trends in Genetics 23 (4): 167-172. https://doi.org/10.1016/j.tig.2007.02.001 Hebert PN, Gregory TR (2005) The promise of DNA barcoding for taxonomy. Systematic Biology 54 (5): 852-859. https://doi.org/10.1080/10635150500354886 Huang D, Meier R, Todd PA, Chou LM (2008) Slow mitochondrial COI sequence evolution at the base of the metazoan tree and its implications for DNA barcoding. Journal of Molecular Evolution 66 (2): 167-174. htips://doi.org/10.1007/ s00239-008-9069-5 Lasker HR, Coffroth MA (1983) Octocoral distributions at Carrie Bow Cay, Belize. Marine Ecology Progress Series 13: 21-28. https://doi.org/10.3354/meps013021 Leray M, Yang J, Meyer C, Mills S, Agudelo N, Vincent R, Boehm J, Machida R (2013) Data from: A new versatile primer set targeting a short fragment of the mitochondrial COI region for metabarcoding metazoan diversity: application for characterizing coral reef fish gut contents. Frontiers in Zoology 10: 34. httos://doi.org/10.5061/ DRYAD.6GD51 McFadden CS, Tullis ID, Hutchinson MB, Winner K, Sohm JA (2004) Variation in coding (NADH dehydrogenase subunits 2, 3, and 6) and noncoding intergenic spacer regions of the mitochondrial genome in Octocorallia (Cnidaria: Anthozoa). Marine Biotechnology 6 (6): 516-526. https://doi.org/10.1007/s10126-002-0102-1 McFadden CS, Benayahu Y, Pante E, Thoma JN, Nevarez PA, France SC (2010) Limitations of mitochondrial gene barcoding in Octocorallia . Molecular Ecology Resources 11 (1): 19-31. https://doi.org/10.1111/j.1755-0998.2010.02875.x McFadden CS, Brown AS, Brayton C, Hunt CB, van Ofwegen LP (2014) Application of DNA barcoding in biodiversity studies of shallow-water octocorals: molecular proxies agree with morphological estimates of species richness in Palau. Coral Reefs 33 (2): 275-286. https://doi.org/10.1007/s00338-013-1123-0 Meyer CP (2003) Molecular systematics of cowries (Gastropoda: Cypraeidae) and diversification patterns in the tropics. Biological Journal of the Linnean Society 79 (3): 401-459. https://doi.org/10.1046/j.1095-8312.2003.00197.x Moura CJ, Cunha MR, Porteiro FM, Rogers AD (2011) The use of the DNA barcode gene 16S mRNA for the clarification of taxonomic problems within the family Sertulariidae (Cnidaria, Hydrozoa). Zoologica Scripta 40 (5): 520-537. htips:// doi.org/10.1111/j.1463-6409.2011.00489.x Oliverio M, Barco A, Modica MV, Richter A, Mariottini P (2009) Ecological barcoding of corallivory by second internal transcribed spacer sequences: hosts of coralliophiline gastropods detected by the cnidarian DNA in their stomach. Molecular Ecology Resources 9 (1): 94-103. https://doi.org/10.1111/j.1755-0998.2008.02388.x Palumbi SR, Martin A, Romano S, McMillan WO, Stice L, Grabowski G (1991) The Simple Fool’s Guide to PCR, Version 2. University of Hawaii Press, Honolulu 14 Morin J et al Prada C, Hellberg ME (2013) Long prereproductive selection and divergence by depth in a Caribbean candelabrum coral. Proceedings of the National Academy of Sciences 110 (10): 3961-3966. https://doi.org/10.1073/pnas.1208931110 Rach J, DeSalle R, Sarkar IN, Schierwater B, Hadrys H (2008) Character-based DNA barcoding allows discrimination of genera, species and populations in Odonata . Proceedings of the Royal Society B: Biological Sciences 275 (1632): 237-247. httops:// doi.org/10.1098/rspb.2007.1290 Ratnasingham S, Herbert PDN (2007) BOLD: The Barcode of Life Data system (http:// www.barcodinglife.org). Molecular Ecology Notes 7 (3): 355-364. https:// doi.org/10.1111/j.1471-8286.2007.01678.x Ratnasingham S, Hebert PN (2013) A DNA-Based Registry for All Animal Species: The Barcode Index Number (BIN) System. PLoS ONE 8 (7): e66213. htips://doi.org/10.1371/ journal.pone.0066213 Ruzicka RR, Colella MA, Porter JW, Morrison JM, Kidney JA, Brinkhuis V, Lunz KS, Macaulay KA, Bartlett LA, Meyers MK, Colee J (2013) Temporal changes in benthic assemblages on Florida Keys reefs 11 years after the 1997/1998 El Niho. Marine Ecology Progress Series 489: 125-141. https://doi.org/10.3354/meps10427 Sanchez JA, Zea S, Diaz JM (1997) Gorgonian communities of two contrasting environments from oceanic Caribbean atolls. Bulletin of Marine Science 61: 61-72. Sanchez JA, Diaz JM, Zea S (1998) Octocoral and black coral distribution patterns on the barrier reef-complex of Providencia Island, southwestern Caribbean. Caribbean Journal of Science 34: 250-264. Sanchez JA, McFadden CS, France SC, Lasker HR (2003) Molecular phylogenetic analyses of shallow-water Caribbean octocorals. Marine Biology 142 (5): 975-987. https://doi.org/10.1007/s00227-003-1018-7 Sanchez JA, Wishing HH (2005) A Field key to the identification of tropical western Atlantic zooxanthellate octocorals (Octocorallia: Cnidaria). Caribbean Journal of Science 41 (3): 508-522. Sanchez JA (2016) Diversity and evolution of octocoral animal forests at both sides of tropical America. Marine Animal Forests. htips:// doi.org/10.1007/978-3-319-17001-5 39-1 Shearer TL, van Oppen MJH, Romano SL, Woérheide G (2002) Slow mitochondrial DNA sequence evolution in the Anthozoa (Cnidaria). Molecular Ecology 11 (12): 2475-2487. https://doi.org/10.1046/j.1365-294x.2002.01652.x Stat M, Huggett MJ, Bernasconi R, DiBattista JD, Berry TE, Newman Su, Harvey ES, Bunce M (2017) Ecosystem biomonitoring with eDNA: metabarcoding across the tree of life in a tropical marine environment. Scientific Reports 7 (1): 12240. htips:// doi.org/10.1038/s41598-017-12501-5 Thomsen PF, Kielgast J, lversen LL, Maller PR, Rasmussen M, Willerslev E (2012) Detection of a diverse marine fish fauna using environmental DNA from seawater samples. PLoS One 7 (8): e41732. https://doi.org/10.1371/journal.pone.0041732 Tsounis G, Edmunds PJ, Bramanti L, Gambrel B, Lasker HR (2018) Variability of size structure and species composition in Caribbean octocoral communities under contrasting environmental conditions. Marine Biology 165 (2): 1-14. https:// doi.org/10.1007/s00227-018-3286-2 Reference DNA barcodes and other mitochondrial markers for identifying ... 15 van der Ham JL, Brugler MR, France SC (2009) Exploring the utility of an indel-rich, mitochondrial intergenic region as a molecular barcode for bamboo corals (Octocorallia: Isididae). Marine Genomics 2: 183-192. https://doi.org/10.1016/j.margen.2009.10.002 Wirshing HH, Messing CG, Douady CJ, Reed J, Stanhope My, Shivji MS (2005) Molecular evidence for multiple lineages in the gorgonian family Plexauridae (Anthozoa: Octocorallia). Marine Biology 147 (2): 497-508. https://doi.org/10.1007/ $00227-005-1592-y Zaiko A, Schimanski K, Pochon X, Hopkins GA, Goldstien S, Floerl O, Wood SA (2016) Metabarcoding improves detection of eukaryotes from early biofouling communities: implications for pest monitoring and pathway management. Biofouling 32 (6): 671-684. https://doi.org/10.1080/08927014.2016.1186165 Zheng L, He J, Lin Y, Cao W, Zhang W (2014) 16S rRNA is a better choice than COI for DNA barcoding hydrozoans in the coastal waters of China. Acta Oceanologica Sinica 33 (4): 55-76. https://doi.org/10.1007/s13131-014-0415-8