# Living on the edge: Assessing the diversity of South African Pocillopora on the margins of the Southwestern Indian Ocean

Brent Chiazzari, Hélène Magalon, Pauline Gélin, Angus Macdonald

► To cite this version:

Brent Chiazzari, Hélène Magalon, Pauline Gélin, Angus Macdonald. Living on the edge: Assessing the diversity of South African Pocillopora on the margins of the Southwestern Indian Ocean. PLoS ONE, 2019, 14 (8), pp.e0220477. 10.1371/journal.pone.0220477. hal-03330875

HAL Id: hal-03330875

https://hal.univ-reunion.fr/hal-03330875

Submitted on 1 Sep 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Living on the edge: Assessing the diversity of South African *Pocillopora* on the margins of the Southwestern Indian Ocean

Brent Chiazzari<sup>1*</sup>, Hélène Magalon<sup>2</sup>, Pauline Gélin<sup>2</sup>, Angus Macdonald<sup>1</sup>

1 School of Life Sciences, University of KwaZulu-Natal, Westville, KwaZulu-Natal, South Africa, 2 UMR ENTROPIE (Université de La Réunion, IRD, CNRS), Laboratoire d'excellence CORAIL, Faculté des Sciences et Technologies, La Réunion, France

* brent.chiazzari@gmail.com

<figure><img src="image_2.png" /></figure>

## Abstract

Scleractinia of the Maputaland reef complex (MRC) in South Africa exist at the margins of the Western Indian Ocean (WIO) coral distribution and are the only substantial hermatypic coral communities in South Africa. *Pocillopora* species occupy a conspicuous component of the MRC, and previous investigations identified three species of *Pocillopora* utilizing conventional taxonomy. Thus, our aims were four-fold: to elucidate *Pocillopora* species diversity using genetic techniques, primarily using species delimitation methods based on the ORF gene; to test for the presence of hybridisation within the *Pocillopora* community on the South-West margin of distribution in the Indian Ocean using two nuclear and two mitochondrial markers; to test the presence of cryptic species, using 13 microsatellite markers, finally, to elucidate the degree of genetic diversity within each *Pocillopora* species found and compare this to communities in lower latitudes. We illustrate taxonomic inconsistencies between these inventories and our phylogenetic data. The MRC harbours unique populations of *Pocillopora*, consisting of three species hypothetically co-occurring throughout the south WIO, namely: *P. meandrina*/*P. eydouxi*, commonly misidentified as *P. verrucosa*, *P. verrucosa*, sometimes correctly identified, but also commonly misidentified as *P. damicornis* *sensu lato*, and *P. villosa*, almost always misidentified as *P. eydouxi*. The hypothesis that hybrid swarms of *Pocillopora* occur in marginal environments such as the MRC was not supported, with low levels of introgressive hybridization reported instead. Analyses illustrate low genetic diversity at the species and population resolutions, suggesting a small founder population for each species. Nevertheless, these populations are demographically unique, exhibiting high levels of ITS2 haplotype endemism compared to higher latitude populations and the rest of the WIO. *Pocillopora* diversity on the MRC represents a unique assemblage and warrants further protection.

<figure><img src="image_3.png" /></figure>

### OPEN ACCESS

**Citation:** Chiazzari B, Magalon H, Gélin P, Macdonald A (2019) Living on the edge: Assessing the diversity of South African *Pocillopora* on the margins of the Southwestern Indian Ocean. PLoS ONE 14(8): e0220477. https://doi.org/10.1371/journal.pone.0220477

**Editor:** Chaolun Allen Chen, Academia Sinica, TAIWAN

**Received:** March 28, 2019

**Accepted:** July 16, 2019

**Published:** August 2, 2019

**Copyright:** © 2019 Chiazzari et al. 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 author and source are credited.

**Data Availability Statement:** The data underlying the results presented in the study are available from NCBI GenBank (accession numbers: ORF: MK359369 - MK359422, ITS2: MK385611 - MK385625 and MK359423 - MK359433, CR: MK429833 - MK429834, HSP70A: MK429806 - MK429832).

**Funding:** This work was supported by a personal grant to BC from the National Research Foundation of South Africa, and a portion of the laboratory work was funded by SANBI as part of the SeaKeys
project (FBI Marine grant). HM funded the microsatellite laboratory work through personal university stipend. AM funded the bulk of laboratory work and field work. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

## Introduction

Modern developments in phylogenetics and the development of genetics-based methods for delineating Scleractinians has improved our understanding of metazoan phylogeny [1]. *Pocillopora* species are characteristically plastic in colony morphology [2] and have historically challenged taxonomists [3]. Indeed, using morphological characters alone, notably the use of corallum macromorphology has proved troublesome [4] and as such, the use of genetics-based methods is useful to align species identification. Correlating genetic signatures to morphology in *Pocillopora* was proved difficult due to incongruencies between whole colony morphology and genetic lineages [5]. Some doubts have been expressed about the sole use of either of these techniques [6]. Nevertheless, research has aimed to reconcile phenotypic and genotypic traits with success using micromorphometrics (columella proportions) and genetic data in the *Pocillopora damicornis* complex [7,8]. These genetic studies have identified high species diversity in the *Pocillopora* genus; however, they have also highlighted that there may be no practical identification of these plastic species with grading morphs *in situ*.

Some studies [4,5,9] avoided primary morphometric species delineation and focussed on applying contemporary genetic algorithms without *a priori* species designations. These approaches have identified new cryptic lineages from different global biogeographic ocean provinces and linked conspecifics, primarily identified and grouped by morphology [3,7,9-13], though most of these studies were undertaken in isolation. In addition, the fast-paced identification and proposal of new species, the revival of junior synonyms and possible splitting of regional conspecifics, have all contributed to the knowledge about *Pocillopora* lineages and diversification. It has also demonstrated possible inconsistencies in nomenclature and allowed for taxonomic congruency of newly recognised genetic strains between studies. Using methodologies previously outlined [14], original and previously published mitochondrial and nuclear data were synthesized into haplotype collections without a primary focus on species membership [4]. The use of multiple markers and methods to assign Primary Species Hypotheses (PSH), and the addition of further analyses (microsatellite markers) to assign Secondary Species Hypotheses (SSH), form a robust basis for delimitation of true evolutionary units. These methods and the naming and cataloguing of the genetic diversity of the *Pocillopora*, provide a collated reference dataset of genetic marker information on the phylogenetics of *Pocillopora* including microsatellite data. This accessioned dataset provides a powerful referencing tool for aligning localised genetic diversity of *Pocillopora*, even if species designations are not yet known, so as to avoid pre-emptive designations based on morphology alone.

Scleractinia of the Maputaland reef complex (MRC) in South Africa exist at the margins of the Western Indian Ocean (WIO) reef coral range and are the most diverse hermatypic coral communities in South Africa. These coral communities are linked to the WIO centre of diversity by Mozambican Channel eddies that move warm water southward from the Equator transporting coral larvae southward [15]. This and the favourable stable East African marine environment may account for the comparatively diverse coral assemblages found on this high latitude reef environment. *Pocillopora* species occupy a conspicuous component of the MRC, covering about 1% of total reef cover, with hard coral abundance having increased in the last 15 years [16]. The diversity of South African hermatypic corals has been previously investigated [17,18], and they identified three species of *Pocillopora* (*P. damicornis*, *eydouxi*, and *verrucosa*) using conventional taxonomy (i.e. based on morphological criteria). However, these communities have yet to be investigated using molecular methods, which would improve the understanding of these corals. We predicted that evolutionary novelties could occur within this marginal reef assemblage where, for example, hybridisation may favour adaptation in this marginal environment [19-22]. It has been reported that genetic structuring of two *Acropora*
species occurring in the MRC indicates a bidirectional flow of genetic diversity [23,24] as is the case of other species' passive dispersal against the Agulhas Current possibly due to strong near-shore counter currents and other oceanographic phenomena [25]. These phenomena could have an effect on the perceived southward attenuation in genetic diversity in *Pocillopora*.

Species delineation of *Pocillopora* using genetics has improved in the past few years, with novel and complete datasets [3-5,7-9,11,12,26-30], and methods available for scientists to use [3-5,7,9-12,29,31-34]. Using these methods, the objective was to produce DNA sequences for the mitochondrial ORF and D-loop, the nuclear HSP70A and ITS2 and microsatellite data, a *priori* of morphological designation. Our aims were four-fold: First; to elucidate *Pocillopora* species diversity, primarily based on the ORF gene defining Primary Species Hypotheses (PSH), second; to test for the presence of hybridisation within the *Pocillopora* community on the South-West margin of distribution in the Indian Ocean using two nuclear and two mitochondrial markers, third; to define Secondary Species Hypotheses (SSH) and the presence of cryptic species, using 13 microsatellite markers, and finally forth: elucidate the degree of genetic diversity within each *Pocillopora* species found using all markers. These data are compared to already collated datasets, thus allowing investigation of the unique position and structuring of the MRC *Pocillopora*. We hypothesise that the MRC may hold similar levels of genetic diversity for the *Pocillopora* as lower latitude reef systems of the WIO region and harbour unique lineages new to science. Here we present an investigation into the diversity of *Pocillopora* using contemporary markers and techniques.

# Materials and methods

## Identification, collection, and preservation

Fifty-seven *Pocillopora* colonies were photographed and sampled from the Northern, Central, and Southern MRC within the Isimangaliso Wetland Park (Fig 1) covering a distance of 100 km, between 2014 and 2017 using the collection permits RES2012/66 and RES2014/06 from the South African Department of Forestry and Fisheries. Emphasis was placed on documenting as many colony forms as possible. To reduce the chance of collecting clonemates, only colonies that were more than 5 m apart were collected. Apical branches of *Pocillopora* colonies (~10 cm length) were collected using SCUBA and free diving from depths ranging from 0 to 35 m and encompassing all reef habitats. Finally, colony fragments were stored at room temperature in 100% ethanol for processing.

Each colony was assigned to a morphotype (named numerically) based on corallum macro-morphology using *in situ* photographs, according to widely accepted descriptions [6]. Using colony form as a guideline, and through personal communication and inspection with a previous checklist author (Dr Louis Celliers) and coral field-expert (Dr David Obura), colonies were compared and aligned with the morphotype and species description of previous checklists focused on the MRC [17,18].

## DNA extraction, PCR amplification, and sequencing

Whole genomic DNA was extracted from each colony using a Zymogen Quick-gDNA Mini-Prep extraction kit following overnight proteinase-K digestion, according to the manufacturer's protocol (Zymo Research Corporation, California, USA). PCR was conducted to amplify the Open Reading frame (ORF) amplified using the FATP6.1 and RORF primers [31], the putative control region (CR) in the mtDNA D-loop structure using the FNAD5.2deg and RCOI primers [31], Heat Shock Protein (HSP70A) using the HSP70A-F and HSP70A-R primers [10], and the Internal Transcribed Spacer region 2 (ITS2) using the ITSc2-5 and R28S1 primers [31].
<figure><img src="image_4.png" /><figcaption>Fig 1. Map of the Maputaland Reef Complex within South Africa's Isimangaliso Wetland Park from which all samples were collected.</figcaption></figure>

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

PCRs were performed in 25 µL of 1X OneTaq Quick-Load (New England Biolabs) 1X Master Mix, 0.3 µM of both forward and reverse primers, 1 µM of bovine serum albumin (BSA), 2–20 ng/µL of extracted DNA template. The PCR thermal cycle was: [95°C for 5 min], 30 × [(94°C for 30 s) (55°C for 45 s for ORF, 58°C for 45 s for D-loop, and 58°C for 45 s for ITS2) (72°C for 45 s)] and [72°C for 10 min], [4°C for ∞]. PCR product for each sample was checked for quality and quantity using a Nanodrop 3000 spectrophotometer and subsequently run on a 1% agarose gel, using 1 µL of loading dye and 5 µL of PCR product per sample. Samples were then sequenced using an ABI 3730 capillary sequencer using Big End Dye technology at Inqaba Biotechnical Industries (Pty) Ltd, Hatfield, South Africa.

Colonies were further genotyped using 13 microsatellite loci: Pd2-001, Pd3-004, Pd3-005, Pd2-006, Pd3-008, Pd3-009 [35], PV2, PV7 [36], Poc40 [9], Pd4, Pd11, Pd13 [37] and Pd3-EF65 [38]. Forward primers were indirectly fluorochrome labelled (6-FAM, VIC, NED, PET) by adding a universal M13-tail at the 5'-end (5' -ACGACGTTGTAAAACGAC-3') [39] and were multiplexed post-PCR. Each amplification reaction was performed in a total volume of 10 µL: 1X of MasterMix Applied (Applied Biosystems), 0.025 µM of forward primer tagged with the M13-tail, 0.25 µM of reverse primer, 0.25 µM of fluorescent dyed M13-tail and 2 ng/
µL of genomic DNA. The thermocycling program was: [94°C for 5 min], 7 × [94°C for 30 s, 62°C (-1°C at each cycle) for 30 s, 72°C for 30 s], 30 × [94°C for 30 s, 55°C for 30 s, 72°C for 30 s], 8 × [94°C for 30 s, 56°C for 30 s, 72°C for 30 s] for the final hybridization of fluorescent dyed M13-tail + 72°C for 5 min. PCR products were genotyped using an ABI 3730XL sequencer at the Plateforme Gentyane (INRA, Clermont- Ferrand, France). Allelic sizes were determined with GENEMAPPER 4.0 (Applied Biosystems) using an internal size standard (Genescan LIZ-500, Applied Biosystems).

## Sequence editing and phylogenetic analyses

Sequences were checked, edited and aligned using Bioedit 7.2.6 [40]. A collated list of haplotypes of the ORF (55 haplotypes), D-loop (32 haplotypes), and ITS2 (128 haplotypes) genes [4], were used as a reference for comparison to local South African congeners. ORF has consistently demonstrated suitability in defining species boundaries [3–5,11,34], thus we present our PSH findings primarily based on the ORF gene. The remaining genes were used in conjunction with the ORF gene, mostly for corroboration to describe new genetic lineages of the genes within the corals studied.

The ITS2 sequences were considered autologous and double peaks were resolved with reverse sequences and then by using PHASE 2.1.1 [41]. ITS2 haplotype sets were reconstructed for all sequences using Phase 2.1.1 with 100 burn-in iterations, followed by 100 iterations, all with different random starting seeds and a threshold probability of 0.9. Haplotype sets were constructed in DnaSP v6 [42], and final ITS2 haplotype sets all had single nucleotide polymorphisms (SNPs) of probability 1.0. Gaps in all sequence sets were counted and all variable sites were used.

MrModeltest [43] was used to determine the appropriate substitution model of evolution for each gene (ORF: K80, CR: K80, ITS2: K80, HSP70A: HKY) for Bayesian phylogenies using MrBayes [44]. All haplotypes were assigned the same weight. The Bayesian trees constructed with haplotypes were created using four Markov Chain Monte Carlo (MCMC) iterations of 20×10⁶ generations each, sampled every 10 generations. The first 20×10³ trees were discarded as burn-in, with the rest of the genealogies used to construct a 50% majority-rule consensus tree. Ultrametric phylogenies using BEAST 2.0 were completed using the GTR+I+G model for each marker [45]. Four MCMC chains were run for 20×10⁶ generations and sampled every 20×10³ generations. Parameter convergence (ESS values were at least 800 for all trees) was verified using Tracer [46], and tree convergence topology was checked using AWTY [47] for all Bayesian analyses. The consensus trees (maximum clade credibility tree; 10% burn-in) were constructed with TREEANNOTATOR 1.7 [48].

Maximum Likelihood (ML) methods were conducted using phyML [49], using the same substitution models calculated in MrModeltest [43], and results were used for the Automatic Barcode Gap Discovery (ABGD) species delimitation method by analysing the distribution of pairwise differences within a dataset. Trees were rooted using two sister genera within the Pocilloporidae; *Stylophora* [50] and *Seriatopora* [51], accession numbers: KX618661–KX618678 (for ORF, CR, and ITS2) and JX869141–JX869151 for HSP70A, after a BLAST search [52]. Bootstrap values greater than 75% and posterior probabilities of 0.95 were accepted as robust and included in trees [53]. Outgroups were pruned for species delimitation analyses and haplotype networks. Haplotype networks were created using PopART [54] using the minimum spanning network construction [55].

## Species delimitation

The methods used to delimit putative species boundaries for the ORF, D-loop, and HSP70A genes, were the Automatic Barcode Gap Discovery (ABGD [56]), Poisson Tree Processes
model (PTP [57]), Bayesian PTP (bPTP), and the generalized mixed Yule coalescent model (GMYC; [58]). Each analysis employed equal weighting of one individual per haplotype, as this has been shown to give the most parsimonious and congruent results among delimitation methods [58]. The four species delimitation methods were used in combination to identify species boundaries using the PSH method, considering how often each method partitioned haplotypes into species. Fields for recombination (FFR) were used to distinguish species breaks using the ITS2 gene, following published criteria [59].

ABGD was used to develop genetic distance-based methods to compare non-overlapping intra and interspecific genetic distance, using the haplotype sets uploaded to the online server: http://wwwabi.snv.jussieu.fr/public/abgd/abgdweb.html. The priors set for each marker were the same: Pmin = 0.001, Pmax = 0.08 and X = 1.5, using the Kimura (K80-2P) model of evolution [60]. PTP methods were used to test the number of substitutions among leaves, where the probability of a substitution tested gives rise to a species, assuming a Poisson distribution. PTP and Bayesian implementation of PTP (bPTP) methods were carried out using the ML tree as per suggestions on the PTP online server: http://species.hits.org/ptp/ [57]. The bPTP analysis was run for 2×10⁶ MCMC generations and a burn-in of 5×10⁵, thinning every 200 generations, and convergence was checked visually as recommended [57]. GMYC methods were carried out using the GMYC online server: http://species.h-its.org/gmyc/ using previously described methods [61]. Ultrametric phylogenies were first calculated using BEAST 2 [45]. Priors were set in BEAUTi [45] following a strict Yule process, i.e. no extinction using a relaxed log-normal clock and assuming a neutral coalescence tree [62]. Three MCMC iterations were run for 20×10⁶ generations and sampled every 2×10³ generations. Chain convergence was assessed using TRACER 1.6 [46] and tree convergence topology was checked using AWTY [47]. The bGMYC results were run using the cut-off probability of 80% [bGMYC (80)]. The consensus trees (maximum clade credibility tree; 10% burn-in; tree not presented) were constructed with TREEANNOTATOR 1.7 [48].

## SSH delimitation using microsatellite data

Genotyped individuals were assigned to an SSH with reference to a comprehensive published dataset [4]. The sixteen PSHs identified were used as a point of reference for species clusters in our data set. Population structure was assessed using Bayesian clustering analysis in the program Structure 2.3.4 [63]. The analysis was run for the whole dataset using five MCMC simulations per run, with 10⁶ iterations and 5×10⁵ steps as burn-in, to test possible clusters (K) between one and ten. Structure harvester was used to determine the most likely combination of homogenous clusters [64].

## Results

### Sequence variability and haplotype assignment

The fifty-seven colonies collected from the MRC yielded three ORF haplotypes (24 ORF27, 19 ORF46, 11 ORF39), of which three individuals did not amplify. These haplotypes are previously published [4], with 15 variable sites over 842 base pairs (bp) (π = 0.008) and a G+C content of 40.0%. The trimmed and aligned CR made up 464 bp of the D-loop region, beginning from position 441 bp of D-loop. The CR sequences yielded two unique haplotypes, dissimilar from haplotypes of previously published lists, due to a single adenine insertion at position 322. Disregarding this insertion, all but four individuals would belong to previously recorded D-loop04 (n = 19) and D-loop12 (n = 28) sensu [4]. The CR haplotypes varied over 9 sites of the 464 bp (π = 0.005) and a G+C content of 40.2%. ITS2 haplotypes consisted of 31 haplotypes with four previously published haplotypes and 27 new ones, and varied over 11 sites of the 710
bp (π = 0.003) with a G+C content of 57.6%. Nuclear HSP70A haplotypes varied over 9 sites of the 741 bp (π = 0.003) and had a G+C content of 48.7%, these sequences could not be compared to established lists due to their short length, compared to the whole HSP70 region (HSP70A + HSP70B). Forty-nine colonies had the CR and ORF, and HSP70A and ITS2 regions concatenated for analysis. The mtDNA (CR+ORF) sequences had 24 variable sites over 1724 bp (π = 0.0068), with a G+C content of 39.4%. The phased nDNA (HSP70A+ITS2) sequence set consisted of 52 reconstructed haplotypes, and 36 variable sites over 1230 bp (π = 0.0061), with a G+C content of 52.8%. Unique sequences obtained using the four genes were deposited into GenBank [65] with the accession numbers: ORF: MK359369—MK359422, ITS2: MK385611—MK385625 and MK359423—MK359433, CR: MK429833—MK429834, HSP70A: MK429806—MK429832.

**Mitochondrial DNA phylogenetic analysis and species delineation**

Bayesian, ultrametric and ML phylogenies for the ORF were congruent using single and multi-
ple individuals per haplotype for the ABGD analysis when using GTR and HKY model param-
eters. The time-calibrated ORF phylogram formed three significantly distinct and robust
clades, for each of the three haplotypes, and correlated with the *a priori* morphotypes (Fig 2).
Calculated divergence time for ORF27 from ORF39 and ORF46 was approximately 2 million
years ago (MYA) with a 13 bp difference, while ORF39 and ORF46 shared a divergence time

<figure><img src="image_5.png" /><figcaption>**Fig 2.** *In situ* photographs of the three colony morphotypes collected at the Maputaland reef complex, South Africa. Coloured borders correspond to the assigned Secondary Species Hypotheses (SSH) (green: SSH09a, orange: SSH13a, blue: SSH13b).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0220477.g002
<figure><img src="image_6.png" /><figcaption>**Fig 3.** Phylogenetic analysis of *Pocillopora* sp. of South Africa in the Southwestern Indian Ocean (SWIO), using the mtDNA region: Open Reading Frame (ORF). (A) Ultrametric phylogram of the three ORF haplotypes. Node values indicate Ultrametric probability/Bayesian probability/ Maximum Likelihood values. (B) Species delimitation methods using ABGD, ML PTP, B PTP, and GMYC methods using Primary Species Hypothesis (PSH) criteria, and Secondary Species Hypothesis (SSH) using thirteen microsatellite markers. (C) Haplotype network coloured according to PSH delineation and sized according to the number of colonies. Scale bar represents millions of years before present.</figcaption></figure>

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

of approximately 1.8 MYA (Fig 3) with a 3 bp difference. The control region (CR) data resulted in a similarly well-supported phylogram of three haplotypes, although it resolved weaker relationships between our *a priori* morphotype allocations and PSH than the ORF (S1 Fig). Divergence times using the CR were identical to the ORF. CR could not be compared to previous studies due to shorter sequence lengths in this dataset. Nevertheless, the 1 bp insert at position 322 is enough to determine that at least one unique D-loop haplotype exists in the studied colonies.

All delimitation methods resolved the same two putative species for the ORF gene: the first species corresponding to ORF46 and ORF39 colonies and the second one to ORF27 colonies. ABGD resolved these two clades, albeit with low statistical significance (p = 0.68). Similar results were obtained by employing GMYC methods. These PSHs were named in accordance with previous species delimitations [4]: PSH09 and PSH13. These PSHs show barcoding gaps with a comprehensive range of haplotypes not present in South Africa, and have further divisions corresponding to SSH using assignment tests with described microsatellite data [4]. Thus, the haplotypes in this study are represented with reference to these SSHs, namely SSH09a (ORF27), SSH13a (ORF39), and SSH13b (ORF46) to gain a global perspective. Here, the three haplotypes do not share any individuals in terms of assigned *a priori* morphotypes according to the network analysis. These results confirm the reproducibility of methods and congruency of results [4,33].

The division of ORF39 and ORF46 (corresponding respectively to SSH13a and SSH13b) was not clear using sequence-based delineation methods. However, the phylogram, sl-FFR, and microsatellite analysis showed a significant and congruent delineation of the SSH13a and SSH13b haplotypes. SSH identification using 13 microsatellite markers resulted in the delineation of SSH13a and SSH13b, corresponding with our three ORF haplotypes and the three *a priori* morphotypes (Fig 3B). Microsatellite analysis resulted in three clusters for individuals of ORF27 (SSH09a); namely SSH09a-1, -2, and -3 [33]. These three clusters fell within the SSH09a delineation and showed no clear discerning colony morphology within SSH09a.
CR data (S1 Fig) were shorter in length and were congruent with the results of the ORF gene. The concatenated mtDNA (ORF+CR) resulted in five haplotypes (S2 Fig). This composite sequence set could not be compared to existing datasets due to the shorter CR length. Nevertheless, this dataset corroborates the ORF data. Thus, we only present the results from the ORF; CR and the mtDNA analysis in the supplementary information.

## Nuclear DNA phylogenetic analysis

The nuclear HSP70A data were congruent with the ITS2 dataset, even though only half the HSP70 gene was used (S2 Fig). The concatenated nDNA (ITS2+HSP70A) dataset resulted in 56 haplotypes (S1 Fig). These data were not included in our argument as they were mostly congruent with the primary data (the ORF and ITS2 gene, and the microsatellite data). These data were excluded as they were mostly not comparable, nor were they used consistently across previous studies. Thus, for the sake of continuity and reducing redundancy, we report these results as a supplementary addition to the study.

SSH09a was composed of 16 ITS2 haplotypes (two shared: one each with SSH13a and SSH13b). SSH13a was composed of six haplotypes (one shared with SSH09a). SSSH13b was composed of 10 haplotypes (one shared with SSH09a). Ten of the ITS haplotypes (Haplotypes ITS2-129 to ITS2-138) were novel. Four previously discovered haplotypes are present in South African populations: ITS2-002, ITS2-009, and ITS2-014 *sensu* [4]. The ITS2-014 haplotype was only associated with individuals of SSH09a for the first time. The ITS2-014 haplotype has only previously been reported [4] in association with ORF39 (SSH13a) and ORF47 (SSH13b).

ITS2 phylogenies were mostly congruent with the ORF data for the three major clades but with weaker support, more phylogenetic structure, and the mismatch assignment of some individuals to different clades. The tree topology can be attributed to a large number of haplotypes, the polymorphic nature of the nuclear ITS gene, and the possibility of hybridisation within the *Pocillopora* complex. Thus, due to its congruency with the ORF gene, these results are not presented here. The network analysis (Fig 4) revealed a complex of mostly homogeneity in colony form within the SSHs (based on ORF and microsatellite assignment tests) with two haplotypes consisting of more than one *a priori* morphotype (ITS2-134 and ITS2-142). Ten single locus co-occurring haplotypes were present (red dashed lines, Fig 4), two of which occur across PSH and SSH boundaries (solid red line, Fig 4). These co-occurring haplotypes delineate two pools of sl-FFR (Single locus fields for recombination); the SSH09a and SSH13b group, and the second pool composed of only SSH13a, SSH13a and SSH13b shared no sl-FFR.

## Coral morphology and species identification

Coral colonies were classified into three *a priori* morphotypes using previous descriptions [6] and these related to the morphological species descriptions listed in the two South African checklists [17,18]. Morphotypes 1, 2, and 3 were associated with the species delineation as discussed in the previous paragraphs, namely; species SSH09a, SSH13a, and SSH13b. Thus, each ORF haplotype and corresponding SSH could be matched with species of previous checklists; SSH09a related to *P. verrucosa*, SSH13a related to *P. verrucosa* and *P. damicornis*, and SSH13b related to *P. eydouxi*; according to the descriptions of previous South African checklists [17,18]. The three SSHs were then correlated to current accepted morpho-genetic species designations from contemporary literature [5,7,8,10,11,26,27,29,30,66] and a collated study including all ORF haplotypes [4]; SSH09a related to *P. eydouxi/meandrina*, SSH13a related to *P. verrucosa*, and SSH13b related to a distinct species, sometimes identified as *P. verrucosa*, but for which the name *P. villosa* has been proposed if no authority already exists [4].
<figure><img src="image_7.png" /><figcaption>Fig 4. Internal transcribed spacer 2 (ITS2) haplotype of *Pocillopora* diversity of South Africa in the Southwestern Indian Ocean (SWIO), constructed using minimum spanning. Each haplotype is sized according to the number of colonies and coloured according to the Secondary Species Hypothesis (SSH). Red dashed curved lines represent haplotypes found co-occurring in heterozygous individuals (solid lines indicate co-occurring haplotypes across SSHs). Haplotypes marked with an asterisk indicate newly discovered ITS2 haplotypes, and haplotypes in bold are previously reported in the literature.</figcaption></figure>

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

Thus, *P. verrucosa* previously described and listed in South African waters belong to SSH09a (ORF27) and appears to be congruent globally with the *P. eydouxi/meandrina* complex. *Pocillopora damicornis/verrucosa* previously described and listed in South African waters belongs to the *P. verrucosa* (SSH13a) species globally. Lastly, the *P. eydouxi* species described in South African waters belong to the morphologically and genetically distinct clade SSH13b, with a newly proposed name of *P. villosa* (if no authority already exists). These overall discrepancies between contemporary phylogenetics and morphological identification in the previous checklists [18,67] are evident, seem to be mostly incorrect but consistent, and are highlighted in Table 1.

Discussion

*Pocillopora* species from South Africa appear to belong to three commonly occurring species on the Maputaland reef complex (MRC), according to contemporary genetic delineation techniques; *Pocillopora eydouxi/meandrina* (SSH09a), *P. verrucosa* (SSH13a), and another genetically distinct species (SSH13b, and a proposed name of *P. villosa*). Morphological identification of South African *Pocillopora* seems to suffer from the same difficulties in identification as many regional studies from around the world. These designations are out of sync with contemporary phylogenetics which is now widely accepted. Previous South African checklists [18,67] and numerous studies employing the same taxonomic designations [16,68–72] may be affected. Each of the species identified in this study seems to be affected, and as such, we have tried to refrain from solely using previously accepted names, but also include
<table><caption>Table 1. Relating DNA analysis to morphology and current literature and checklists, of South African *Pocillopora*.</caption><thead><tr><th colspan="2">Literature</th><th colspan="4">This study</th><th colspan="2">Contemporary Literature</th></tr><tr><th colspan="2">South African checklists</th><th colspan="2">Morphology</th><th colspan="2">DNA</th><th colspan="2">Taxonomy</th></tr><tr><th>Riegl (1996)</th><th>Celliers and Schleyer (2003)</th><th>a priori colony morphotype</th><th>Open reading frame Haplotype</th><th>Primary Species Hypothesis</th><th>Secondary Species Hypothesis</th><th>Individual studies</th><th>Collated literature [4]</th></tr></thead><tbody><tr><td>P. verrucosa</td><td>Morphotype 1</td><td>[4]; ORF27[11]; A [5,9]; 1a [7]; em[26]; N/A [27]; N/A [8]; IIb</td><td>PSH09</td><td>SSH09a-1<br>SS09a-2<br>SSH09a-3</td><td>[11] eydouxi, meandrina [5,9] eydouxi, meandrina, verrucosa, damicornis [7] eydouxi, meandrina [26] meandrina [27] eydouxi, meandrina [8] eydouxi, meandrina</td><td>P. eydouxi, meandrina complex in the Indian Ocean.</td><td></td></tr><tr><td>P. verrucosa, or P. damicornis (noticeably in young colonies).</td><td>Morphotype 2</td><td>[4]; ORF39 (5, 10); 3c (27); NF</td><td>PSH13</td><td>SSH13a</td><td>[5,9] verrucosa (27) N/A</td><td>P. verrucosa of the Indian Ocean</td><td></td></tr><tr><td>P. eydouxi</td><td>Morphotype 3</td><td>[4]; ORF46 [5,9]; 3a [7]; γ [11]; N/A [27]; N/A</td><td>PSH13</td><td>SSH13b</td><td>[5,9] damicornis, verrucosa, meandrina [7] γ-verrucosa [11] molokensis [27] N/A</td><td>Genetically distinct subclade in the P. verrucosa complex (known as the "doom-trap"). [4] a new name proposed; P. villosa</td><td></td></tr></tbody></table>

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

the SSHs, which represent the genetic lineage of these species. Added to this, we do reconcile our genetic findings with previous checklists. This dynamic evolutionary environment suggests this marginal habitat-in terms of hermatypic coral assemblages-may harbour unique phylogenetic diversity and limited hybridisation amongst these species.

## Use of mitochondrial, nuclear, and microsatellite markers for *Pocillopora* in South Africa

Nuclear DNA data suggests a unique haplotype profile, but low diversity at the population level when compared with other lower latitude populations. Fifty-seven colonies of *Pocillopora* were sampled for species-level diversity in South Africa for the first time. The MRC marks the southwestern margin of hermatypic coral diversity in the WIO bioregion, where a paucity of research on coral biogeography exists. Recent insight [3] has highlighted the concordance of the ORF marker with holobiont, transcriptomic, SNP, and conspecificity data. Added to this, a collection of data now exists for the marker in the form of a collated list of haplotypes [4]. Thus, species assessment using the ORF marker was appropriate for delineation in South African *Pocillopora* and demonstrates the efficacy of this marker and suggests its use as a principal barcode marker. This data also corroborates and illustrates the reproducibility of our methods developed and employed by many studies over different regions [3,4,7,10,13,30,32-34,37].

The failure of the CR primers [11] to amplify the much longer CR gene across our dataset suggests mutations within the primer annealing sites. The CR exhibits phylogenetic congruency with the considerably shorter ORF gene. The ITS2 marker was easier to amplify than the HSP70A and B regions. Evidence of this was apparent by the high failure rate of the HSP70B marker in this study, and we were forced to use only half the gene region (HSP70A). Although the ITS2 gene has been reasoned to be a poor barcode marker, it may have merit in identifying
population structure, paternal lineages, possible backcrossing, incomplete lineage sorting, and
hybridisation events. Microsatellite data were able to resolve all species groupings with preci-
sion and congruency with the other marker sets. The combination of sequence-based species
delineation, microsatellite genotyping, and FFR methods across multiple markers resolved
reliable results that probably describe the species of *Pocillopora* in South Africa.

Taxonomic assessment of *Pocillopora* in South Africa

South African *Pocillopora* records are at odds with those described by contemporary phyloge-
netic methods (Table 1). Nevertheless, South African checklists [17,18] are largely consistent
with one another, allowing for easier identification of species and comparison with our phy-
logenies (Table 1). We identified three species in South Africa: *P. eydouxi/meandrina*
(SSH09a) *P. verrucosa* (SSH13a), and the distinct newly proposed *P. villosa* (SSH13b). All the
studied colonies assigned to characteristic SSH09a harbouring morphotype 1 (*P. meandrina/P.
eydouxi* morphotype; Table 1), and corresponded to *P. verrucosa* in South African checklists
[17,18]. Morphological similarities between the *P. meandrina/ P. eydouxi* complex and that of
*P. verrucosa* may explain the incorrect taxonomic designation (see previously described simi-
larities [2]). In addition, SSH09a colonies are phenotypically plastic which may blur species
boundaries that have been based solely on morphology explaining why both *P. meandrina* and
*P. eydouxi* names have been used. Nevertheless, the South African *P. eydouxi/meandrina*
(SSH09a) population consists of one haplotype, ORF 27, and is regarded as a distinct species
using our delineation methods. These species boundaries are supported by eight genetic stud-
ies and morphologically designated as *P. eydouxi* with a 75% frequency [3–5,7–9,11,29]. Mor-
phologically, SSH09a does fit conventional descriptions of *P. eydouxi* [2,6]. However, novel
research [3] using RAD-seq data found it is not possible to discern *P. meandrina* morphotypes
from *P. eydouxi* using the ORF marker. In addition, our microsatellite genotyping revealed the
presence of all three population clusters presently known to SSH09a (-1, -2, and -3, Table 1). It
is hypothesised that these three clusters possibly reflect three cryptic species existing in sym-
patry in the WIO [33]. Recently, a method using the Histone-3 region has been successful in
delineating the *P. eydouxi/meandrina* complex [3]. More investigation in this species complex
is needed using this method as our study was not able to correctly discern the morphologies
and genetics.

In a similar fashion to the *P. meandrina/eydouxi* (SSH09a) colonies, the MRC colonies
assigned to SSH13b harboured a single ORF haplotype (ORF46; [4]) and were all assigned to
morphotype 3 (Table 1). Although sequence-based delimitation methods failed to identify the
SSH13b clade as a species, phylogenetic and microsatellite data support this unique clade (Fig
3). Assignment tests [4] based on microsatellite data that this unique haplotype, with little
interspecies genetic divergence, exists in sympatry with SSH13a and SSH13c. Together with its
easily identifiable and unique colony morphology, and the argument that SSH13b should con-
stitute its own species; researchers [4] proposed the name *P. villosa* due to its particular and
recognizable velvety aspect. The SSH13b morphological description matches in some aspects
to the description [6] of *P. elegans* and *P. eydouxi*, and has been previously identified [17,18] as
*P. eydouxi* in South African checklists (Table 1). We believe the argument that SSH13a and
SSH13b represents one morphologically plastic species is unlikely. The evidence against this
possibility includes a unique genetic identity across all of the markers tested with unique
monophyletic and statistically robust clades, no evidence of shared genetic material (haplo-
types) or hybrid backcrosses. Thus, we conclude that *P. villosa* (SSH13b) should constitute its
own species.
*Pocillopora verrucosa* (SSH13a) colonies from the MRC also harboured a single ORF haplotype (ORF39) and they were all assigned to a unique morphotype (morphotype 2, i.e. *P. verrucosa* morphotype; Table 1), which were identified as *P. verrucosa* and *P. damicornis* by Riegl [17] and Schleyer and Celliers [18]. Contrary to SSH09a and SSH13b corresponding to a unique species in the various checklists, *P. verrucosa* (SSH13a) corresponds to two species: *P. verrucosa* and *P. damicornis* (*sensu lato*) in the South African checklists. Indeed, *P. damicornis* (*sensu lato*) is reported in both Scleractinian checklists in South Africa (Table 1). Yet, after a search for unique colony morphotypes, including many individuals that may resemble *P. damicornis*-like colonies, no colonies harbouring the corresponding ORF haplotypes were collected (haplotypes ORF01–ORF21, species PSH01–PSH05). It is possible that young colonies of other *Pocillopora* growing in marginal habitats (such as rock pools of intertidal zones) are mistaken for *P. damicornis*. Indeed, these colonies sometimes show finer branching where verrucae grade into branches and may resemble true *P. damicornis*. Interestingly though, researchers [69] report *P. damicornis* on Nine Mile Reef on the MRC, which possibly suggests healthy, phenotypically representative adults of *Pocillopora* sp. were mistaken as *P. damicornis* (*sensu lato*). Others [68] also report the coral’s presence in most tidal pools in the MRC and as far south as Isipingo, KwaZulu-Natal. Having sampled some of the exact localities of these studies, it is unlikely that we did not account for these morphotypes.

Incorrect morphological identification of small and finer branched colonies belonging to SSH13a (associated with the *P. verrucosa* morphotype) is the most likely case for misidentification of *Pocillopora* in South Africa. The results of our study illustrate the difficulty in identification of cryptic lineages without corroboration employing genetic analysis. It is also clear that certain instances of identification in various other studies were incorrect, owing to the varied morphological assignments in each haplotype, a probable artefact of misidentification [4].

## Genetic diversity and structure of South African *Pocillopora*

Analyses suggest low genetic diversity at the species and population resolutions, as evidenced by the low haplotypic diversity and lack of linkages in the haplotype network (Fig 3). South Africa harbours three *Pocillopora* species [possibly five, considering the three diverging clusters within SSH09a (4)] found in the WIO region. Nevertheless, when comparing haplotypic diversity, South Africa comprised only three of the sixteen haplotypes currently known in the WIO. This suggests a small founder population for each species, probably supplied to southern African reefs from the tropics.

South African populations were also demographically unique, exhibiting high levels of ITS2 haplotype endemism compared to the rest of the WIO (Fig 4). Indeed, distinct structuring in South African *P. verrucosa* compared to lower latitude populations has been found [73]. It is also likely that the same patterns exist for all three SSHs in South Africa given their similar life histories. The genetic demography of South African *Pocillopora* suggests that recruitment at population and generational resolutions is self-sustaining. Our results support these findings [73], reporting *P. verrucosa* in southern Mozambique as a separate cluster from South Africa. This suggests weak connectivity between southern Mozambique and South African reef communities over ecological time scales. Similar patterns of diversity have been reported for other scleractinians in South Africa, such as *Acropora tenuis* [24] and other tropical taxa that demonstrate a southward attenuation in diversity [25]. The physical conditions of the Maputo Bay are generally unsuitable for the colonisation of large hermatypic populations and is a likely altotropic barrier to genetic connectivity between central Mozambique and South Africa. This biogeographic region probably affects the demography of these populations by suppressing population connectivity, possibly explaining the observed population indices.
We also report that nuclear DNA reveals some hybridisation between the *P. eydouxi* and *P. villosa* lineages. *Pocillopora verrucosa* exhibits negligible hybridisation with *P. eydouxi* and *P. villosa*. There is considerable evidence for high levels of hybridisation in marginal marine environments (20–23) and in *Pocillopora* (72,73), leading us to expect greater levels of hybridisation than we observed. Causes for suppressed hybridisation are numerous, however, investigation of these mechanisms is beyond the scope of this study. Nevertheless, we can postulate mechanisms that may account for our observations: Firstly, pre- and post-zygotic barriers such as compatibility, genetic mules, and ecological fitness are common factors that determine the success of hybrids [19]. One-way introgressive hybridization of *Pocillopora* [74] may limit the possibility of numerous hybrid combinations. Secondly, ecological factors may play a role in inhibiting successful hybridisation. For example, *Pocillopora* communities on the MRC may resemble a functioning ecosystem more similar to those of lower latitude *Pocillopora* communities than anticipated [66]. Competition, larval development, niche space, or other factors may inhibit a suitable environment for sympatry of *Pocillopora* hybrids. Indeed, these are testable hypotheses and warrant further investigation.

## Conclusions

A recent drive to redress inconsistencies in *Pocillopora* classification has begun to align global congeners and taxonomic understanding. South Africa's Maputaland Reef Complex (MRC) corals are influenced by its position at the southwestern margin of the Indian Ocean. The MRC harbours unique populations of *Pocillopora*, consisting of three to five species (considering sub-structuring within SSH09a) co-occurring throughout the southwestern Indian ocean, namely: SSH09a corresponding to *P. meandrina*/*P. eydouxi* morphotype (commonly misidentified as *P. verrucosa*), SSH13a corresponding to *P. verrucosa* morphotype (sometimes correctly identified, but also commonly misidentified as *P. damicornis sensu lato*), and SSH13b corresponding to the *P. villosa* morphotype (almost always misidentified as *P. eydouxi*). The results of this study demonstrate the possibility that previously investigated *Pocillopora* species in South Africa may have been misidentified. Nevertheless, due to the precision and consistency of taxonomy of these local studies [68,69,72,73], any incorrect taxonomy is likely to have little effect on conclusions, except for taxonomy based checklists [18,67]. We also report a unique assemblage of nuclear diversity on South African reef complexes. This diversity represents a unique but possibly fragile assemblage of *Pocillopora*, owing to the low observed genetic diversity and minimal hybridisation compared to the rest of the WIO. Considering these populations and their isolation from lower latitude populations and their demographic subsidies, this warrants the MRC's protection and management. This study also demonstrates the universality of using the ORF marker for *Pocillopora* species identification and the importance of understanding the nature of species complexes at range margins.

## Supporting information

S1 Fig. Phylogenetic analysis of *Pocillopora* sp. of South Africa in the Southwestern Indian Ocean (SWIO), using the part of D-loop and COI of the Control Region (CR). (A) Ultrametric phylogram of the two D-loop haplotypes. Node values indicate Ultrametric probability/ Bayesian probability/Maximum Likelihood values. Scale bar represents Millions of Years Before present. (B) Haplotype network coloured according to Primary Species Hypothesis (PSH) delineation and sized according to the number of individuals. (TIF)
**S2 Fig. Phylogenetic analysis of *Pocillopora* sp. of South Africa in the Southwestern Indian Ocean (SWIO), using HSP70A.** (A) Ultrametric phylogram of the 28 HSP70A haplotypes. Node values indicate Ultrametric probability/ Bayesian probability/Maximum Likelihood values. Scale bar represents Millions of Years Before present. (B) Haplotype network coloured according to Primary Species Hypothesis (PSH) delineation and sized according to the number of individuals.

(TIF)

**S3 Fig. Phylogenetic analysis of *Pocillopora* sp. of South Africa in the Southwestern Indian Ocean (SWIO), using the mtDNA region: Open Reading Frame (ORF) and part of the D-loop of the Control Region (CR).** (A) Ultrametric phylogram of the concatenated haplotypes. Node values indicate Ultrametric probability/ Bayesian probability/Maximum Likelihood values. (B) Species delimitation methods using ABGD, ML PTP, B PTP, and GMYC methods using Primary Species Hypothesis (PSH) criteria, and Secondary Species Hypothesis (SSH) using thirteen microsatellite markers. C: Haplotype network coloured according to PSH delineation and sized according to the number of individuals. Scale bar represents Millions of Years Before present.

(TIF)

**S4 Fig. Phylogenetic analysis of *Pocillopora* sp. of South Africa in the Southwestern Indian Ocean (SWIO), using the nDNA region: ITS2 and HSP70A.** (A) Ultrametric phylogram of the concatenated haplotypes. Node values indicate Ultrametric probability/ Bayesian probability/Maximum Likelihood values. (B) Haplotype network coloured according to Primary Species Hypothesis (PSH) delineation and sized according to the number of individuals. Scale bar represents Millions of Years Before present.

(TIF)

**S1 Table. General morphological characters for of *Pocillopora* sp. of South Africa in the Southwestern Indian Ocean (SWIO).** Features of morphological characters presented here are strictly those that can be seen through general *in situ* and photographic observations.

(DOCX)

**Acknowledgments**

We are grateful to members of the Marine Evolution Lab (MAREVOL) that assisted in sampling: Kerry Etsebeth, Carolina Barausse, Kirsty Wiggill, and Jessica Chiazzari. We are also thankful to David Obura, and Louis Celliers who aided in identifying species in the field, and Louis Celliers for aiding in confirming previous information on species identification. Thanks go to Kerry Sink, Mari-Lise Franken and Prideel Madjiet and SANBI for logistical and funding support. We acknowledge the Plateforme Gentyane of the Institut National de la Recherche Agronomique (INRA, Clermont-Ferrand, France) for genotyping and technical support. Lastly, thanks to the three anonymous reviewers, their advice and comments have greatly improved the overall quality of this manuscript.

**Author Contributions**

**Conceptualization:** Brent Chiazzari, Hélène Magalon, Angus Macdonald.

**Data curation:** Brent Chiazzari.

**Formal analysis:** Brent Chiazzari, Pauline Gélin.

**Funding acquisition:** Brent Chiazzari, Angus Macdonald.
**Investigation:** Brent Chiazzari, Angus Macdonald.

**Methodology:** Brent Chiazzari, Hélène Magalon, Angus Macdonald.

**Project administration:** Brent Chiazzari, Angus Macdonald.

**Resources:** Hélène Magalon, Angus Macdonald.

**Supervision:** Angus Macdonald.

**Validation:** Brent Chiazzari.

**Visualization:** Brent Chiazzari.

**Writing – original draft:** Brent Chiazzari.

**Writing – review & editing:** Brent Chiazzari, Hélène Magalon, Angus Macdonald.

## References

1. Kitahara M V., Fukami H, Benzoni F, Danwei Huang. The New Systematics of Scleractinia: Integrating Molecular and Morphological Evidence. Cnidaria, Past, Present Future. 2016; 1–855. https://doi.org/10.1007/978-3-319-31305-4

2. Veron JEN, Pichon M. Scleractinia of Eastern Australia. Part I. Families *Thamnasteriidae*, *Astrocoeniidae*, *Pocilloporidae*. By Veron J.E.N. and Pichon M. Townsville: Queensland: Australian Institute of Marine Science, 1976. 86 pp., 166 figs. (Monograph Series Volume 1.) P. J Mar Biol Assoc United Kingdom. 2009/05/01. Cambridge University Press; 1978;58: 544. https://doi.org/10.1017/S0025315400028289

3. Johnston EC, Forsman ZH, Flot J-F, Schmidt-Roach S, Pinzón JH, Knapp ISS, et al. A genomic glance through the fog of plasticity and diversification in *Pocillopora*. Sci Rep. 2017; 7: 1–11. https://doi.org/10.1038/s41598-016-0028-x

4. Gélin P, Postaire B, Fauvelot C, Magalon H. Reevaluating species number, distribution and endemism of the coral genus *Pocillopora* Lamarck, 1816 using species delimitation methods and microsatellites. Mol Phylogenet Evol. 2017; 109: 430–446. https://doi.org/10.1016/j.ympev.2017.01.018 PMID: 28219759

5. Pinzón JHC, Sampayo EM, Cox EF, Chauka LJ, Chen CA, Voolstra CR, et al. Blind to morphology: Genetics identifies several widespread ecologically common species and few endemics among Indo-Pacific cauliflower corals (*Pocillopora*, Scleractinia). J Biogeogr. Wiley; 2013; 40: 1595–1608. https://doi.org/10.1111/jbi.12110

6. Veron JEN, Smith M. Corals of the world [Internet]. Australian Institute of Marine Science: 2000. https://doi.org/10.1126/Science.1061936

7. Schmidt-Roach S, Miller KJ, Lundgren P, Andreakis N. With eyes wide open: A revision of species within and closely related to the *Pocillopora damicornis* species complex (Scleractinia; *Pocilloporidae*) using morphology and genetics. Zool J Linn Soc. 2014; 170: 1–33. https://doi.org/10.1111/zoj.12092

8. Marti-Puig P, Forsman ZH, Haverkort-Yeh RD, Knapp IS, Maragos JE, Toonen RJ. Extreme phenotypic polymorphism in the coral genus *Pocillopora*; micro-morphology corresponds to mitochondrial groups, while colony morphology does not. Bull Mar Sci. 2014; 90: 1–22.

9. Pinzón JH, LaJeunesse TC. Species delimitation of common reef corals in the genus *Pocillopora* using nucleotide sequence phylogenies, population genetics and symbiosis ecology. Mol Ecol. Blackwell Publishing Ltd; 2011; 20: 311–325. https://doi.org/10.1111/j.1365-294X.2010.04939.x PMID: 21091563

10. Schmidt-Roach S, Lundgren P, Miller KJ, Gerlach G, Noreen AME, Andreakis N. Assessing hidden species diversity in the coral *Pocillopora damicornis* from Eastern Australia. Coral Reefs. 2013; 32: 161–172. https://doi.org/10.1007/s00338-012-0959-z

11. Flot J-F, Magalon H, Cruaud C, Couloux A, Tillier S. Patterns of genetic structure among Hawaiian corals of the genus *Pocillopora* yield clusters of individuals that are compatible with morphology. C R Biol. Elsevier; 2008; 331: 239–247. https://doi.org/10.1016/j.crvi.2007.12.003 PMID: 18280989

12. Tillier S, Flot J-F. Molecular phylogeny and systematics of the scleractinian coral genus *Pocillopora* in Hawaii. Proc 10th Int Coral Reef Symp. 2006; 24–29.

13. Soto D, De Palmas S, Ho MJ, Denis V, Chen CA. Spatial variation in the morphological traits of *Pocillopora verrucosa* along a depth gradient in Taiwan. PLoS One. Public Library of Science; 2018; 13: e0202586. https://doi.org/10.1371/journal.pone.0202586 PMID: 30118513
14. Pante E, Puillandre N, Viricel A, Arnaud-Haond S, Aurelle D, Castelin M, et al. Species are hypotheses: Avoid connectivity assessments based on pillars of sand. Mol Ecol. 2015; 24: 525–544. https://doi.org/10.1111/mec.13048 PMID: 25529046

15. Lutjeharms JRE. Three decades of research on the greater Agulhas Current. Ocean Sci. 2007; 3: 129–147. https://doi.org/10.5194/os-3-129-2007

16. Porter SN, Schleyer MH. Long-term dynamics of a high-latitude coral reef community at Sodwana Bay, South Africa. Coral Reefs. 2017; 36: 369–382. https://doi.org/10.1007/s00338-016-1531-z

17. Riegl BCF of SSAAC. Hermatypic Coral Fauna of Subtropical Southeast Africa: A Checklist! Pacific Sci. 1996; 50: 404–414. Available: http://nsuworks.nova.edu/occ_facarticles/337

18. Schleyer MH, Celliers L. Biodiversity on the marginal coral reefs of South Africa: What does the future hold? Verh Leiden. 2003; 345: 387–400.

19. Willis BL, Van Oppen MJH, Miller DJ, Vollmer S V, Ayre DJ. The Role of Hybridization in the Evolution of Reef Corals. Annu Rev Ecol Evol Syst. 2006; 37: 489–517. https://doi.org/10.1146/annurev.ecolsys.37.091305.110136

20. Hatta M, Fukami H, Wang W, Omori M, Shimoike K, Hayashibara T, et al. Reproductive and genetic evidence for a reticulate evolutionary history of mass spawning coral. Mol Biol Evol. 1999; 16: 1607–1613. https://doi.org/10.1093/oxfordjournals.molbev.a026073 PMID: 10555292

21. Miller KJ, Ayre DJ. The role of sexual and asexual reproduction in structuring high latitude populations of the reef coral *Pocillopora damicornis*. Heredity (Edinb). 2004; 92: 557–568. https://doi.org/10.1038/sj.hdy.6800459 PMID: 15100708

22. Budd AF, Pandolfi JM. Evolutionary Novelty Is Concentrated at the Edge of Coral Species Distributions. Science (80-). 2010; 328: 1558–1561. https://doi.org/10.1126/science.1188947 PMID: 20558718

23. Macdonald AHH, Schleyer MH, Lamb JM. *Acropora austera* connectivity in the south-western Indian Ocean assessed using nuclear intron sequence data. Mar Biol. 2011; 158: 613–621. https://doi.org/10.1007/s00227-010-1585-3

24. Chiazzari B, Macdonald A, Schleyer M. High-latitude connectivity of the scleractinian coral *Acropora tenuis* in the south-western Indian Ocean, identified using nuclear intron and mitochondrial sequence data. African J Mar Sci. 2013; 35: 233–241. https://doi.org/10.2989/1814232X.2013.795915

25. Teske PR, Bader S, Golla TR. Passive dispersal against an ocean current. Mar Ecol Prog Ser. 2015; 539: 153–163. https://doi.org/10.3354/meps11516

26. Paz DA, Alejandro G, Moreno A. Morphological variation and different branch modularity across contrasting Morphological variation and different branch modularity across contrasting flow conditions in dominant *Pocillopora* reef - building corals. 2015; https://doi.org/10.1007/s00442-014-3199-9 PMID: 25556295

27. Hsu C-M, de Palmas S, Kuo C-Y, Denis V, Chen CA. Identification of Scleractinian Coral Recruits Using Fluorescent Censusing and DNA Barcoding Techniques. PLoS One. Public Library of Science; 2014; 9: e1–13. https://doi.org/10.1371/journal.pone.0107366 PMID: 25211345

28. Baums IB, Devlin-durante M, Laing BAA, Feingold J, Smith T, Bruckner A, et al. Marginal coral populations: the densest known aggregation of *Pocillopora* in the Galápagos Archipelago is of asexual origin. 2014; 1: 1–11. https://doi.org/10.3389/fmars.2014.00059

29. Souter P. Hidden genetic diversity in a key model species of coral. Mar Biol. 2010; 157: 875–885.

30. Forsman ZH, Barshis DJ, Hunter CL, Toonen RJ. BMC Evolutionary Biology Shape-shifting corals: Molecular markers show morphology is evolutionarily plastic in *Porites*. 2009; https://doi.org/10.1186/1471-2148-9-45

31. Flot J-F, Tillier S. The mitochondrial genome of *Pocillopora* (Cnidaria: Scleractinia) contains two variable regions: The putative D-loop and a novel ORF of unknown function. Gene. 2007; 401: 80–87. https://doi.org/10.1016/j.gene.2007.07.006 PMID: 17716831

32. Torres AF, Ravago-Gotanco R. Rarity of the “common” coral *Pocillopora damicornis* in the western Philippine archipelago. Coral Reefs. 2018; 37: 1209–1216. https://doi.org/10.1007/s00338-018-1729-3

33. Gélin P, Fauvelot C, Bigot L, Baly J, Magalon H. From population connectivity to the art of striping Russian dolls: the lessons from *Pocillopora* corals. Ecol Evol. 2018; 8: 1411–1426. https://doi.org/10.1002/ece3.3747 PMID: 29375807

34. Johnston EC, Forsman ZH, Toonen RJ. A simple molecular technique for distinguishing species reveals frequent misidentification of Hawaiian corals in the genus *Pocillopora*. PeerJ. 2018; 6: e4355. https://doi.org/10.7717/peerj.4355 PMID: 29441239

35. Starger CJ, Yeoh SS, Dai CF, Baker AC, Desalle R. Ten polymorphic STR loci in the cosmopolitan reef coral, *Pocillopora damicornis*. Mol Ecol Resour. 2007; 8: 619–621.
36. Magalon H, Samadi S, Richard M, Adjeroud M, Veuille M. Development of coral and zooxanthella-specific microsatellites in three species of *Pocillopora* (Cnidaria, Scleractinia) from French Polynesia. Mol Ecol Notes. 2004; 4: 206–208.

37. Torda G, Schmidt-Roach S, Peplow LM, Lundgren P, van Oppen MJH. A rapid genetic assay for the identification of the most common *Pocillopora damicornis* genetic lineages on the great barrier reef. PLoS One. 2013; 8, e58447. https://doi.org/10.1371/journal.pone.0058447 PMID: 23505507

38. Gorospe KD, Karl SA. Genetic relatedness does not retain spatial pattern across multiple spatial scales: dispersal and colonization in the coral, *Pocillopora damicornis*. Mol Ecol. 2013; 22: 3721–3736. https://doi.org/10.1111/mec.12335 PMID: 23786173

39. Schuelke M. An economic method for the fluorescent labeling of PCR fragments. Nat Biotech. 2000; 18: 233–234.

40. Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic acids symposium series. 1999. pp. 95–98.

41. Stephens M, Smith NJ, Donnelly P. A New Statistical Method for Haplotype Reconstruction from Population Data. Am J Hum Genet. 2001; 68: 978–989. https://doi.org/10.1086/319501 PMID: 11254454

42. Librado P, Rozas J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009; 25: 1451–1452. https://doi.org/10.1093/bioinformatics/btp187 PMID: 19346325

43. Nylander JAA. MrModeltest v2. Evolutionary Biology Centre, Uppsala University; 2004.

44. Huelsenbeck JP, Ronquist F. MRBAYES: Bayesian inference of phylogeny. Bioinformatics. 2001; 17: 754–755. https://doi.org/10.1093/bioinformatics/17.8.754 PMID: 11524383

45. Bouckaert R, Heled J, Kühnert D, Vaughan T, Wu C-H, Xie D, et al. BEAST 2: A Software Platform for Bayesian Evolutionary Analysis. PLoS Comput Biol. 2014; 10. https://doi.org/10.1371/journal.pcbi.1003537 PMID: 24722319

46. Rambaut A, Drummond A, Xie D, Baele G, Suchard M. Tracer v1.7. 2018.

47. Wilgenbusch JC, Warren DL, Swofford DL. AWTY: A system for graphical exploration of MCMC convergence in Bayesian phylogenetic inference. [Internet]. 2004. Available: http://ceb.csit.fsu.edu/awty.

48. Rambaut A, Drummond AJ. TreeAnnotator v1.7.0. Available as part of the BEAST package at http://beast.bio.ed.ac.uk. 2013.

49. Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. Phylogenies; New algorithms and methods to estimate maximum-likelihood assessing the performance of PhyML 3.0. Syst Biol. 2010; 59: 307–321. https://doi.org/10.1093/sysbio/syq010 PMID: 20525638

50. Schweigger AF. Handbuch der Naturgeschichte der skelettlosen ungegliederten Thiere. Leipzig. 1820;

51. Lamarck JB. de Monet Histoire naturelle des animaux sans vertèbres: présentant les caractères génér-aux et particuliers de ces animaux, leur distribution, leurs classes, leurs familles, leur genres, et al cita-tion des principales espèces qui s'y rapportent. Ire Partie. 1816; 1–6: 1–343.

52. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mar Biol. 1990; 215: 403–410.

53. Erixon P, Svennblad B, Britton T, Oxelman B. Reliability of Bayesian Posterior Probabilities and Boot-strap Frequencies in Phylogenetics. Syst Biol. Society of Systematic Zoology; 2003; 52: 665–673. https://doi.org/10.1080/10635150390235485

54. Leigh J, Bryant D. PopART: Full-feature software for haplotype network construction. Methods Ecol ogical Evol. 2015; 6: 1110–1116.

55. Bandelt H, Forster P, Röhl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999; 16: 37–48. https://doi.org/10.1093/oxfordjournals.molbev.a026036 PMID: 10331250

56. Puillandre N, Lambert A, Brouillet S, Achaz G. ABGD, Automatic Barcode Gap Discovery for primary species delimitation. Mol Ecol. 2011;

57. Zhang J, Kapli P, Pavlidis P, Stamatakis A. A general species delimitation method with applications to phylogenetic placements. Bioinformatics. 2013; 29: 2869–2876. https://doi.org/10.1093/bioinformatics/btt499 PMID: 23990417

58. Pons J, Barraclough TG, Gomez-Zurita J, Cardoso A, Duran DP, Hazell S, et al. Sequence-based species delimitation for the DNA taxonomy of undescribed insects. Syst Biol. 2006; 55: 595–609. https://doi.org/10.1080/10635150600852011 PMID: 16967577

59. Flot J, Couloux A, Tillier S. Haplowebs as a graphical tool for delimiting species: a revival of Doyle 's "field for recombination" approach and its application to the coral genus *Pocillopora* in Clipperton. BMC Evol Biol. 2010; 10.

60. Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980; 16: 111–120. https://doi.org/10.1007/BF01731581 PMID: 7463489
61. Fujisawa T, Barraclough TG. Delimiting species using single-locus data and the Generalized Mixed Yule Coalescent approach: a revised method and evaluation on simulated data sets. *Syst Biol*. Oxford University Press; 2013; 62: 707–724.

62. Yule GU. A Mathematical Theory of Evolution, Based on the Conclusions of Dr. J. C. Willis, F.R.S. Philos Trans R Soc London Ser B, Contain Pap a Biol character. The Royal Society London; 1925; 213: 21–87.

63. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. *Genetics*. 2000; 155: 945–959. PMID: 10835412

64. Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. *Mol Ecol*. 2005; 14: 2611–2620. https://doi.org/10.1111/j.1365-294X.2005.02553.x PMID: 15969739

65. Clark K, Karsch-Mizrachi I, Lipman DJ, Ostell J, Sayers EW. GenBank. Nucleic Acids Res. 2015/11/20. Oxford University Press; 2016; 44: D67–D72. https://doi.org/10.1093/nar/gkv1276 PMID: 26590407

66. Schmidt-Roach S, Miller KJ, Woolsey E, Gerlach G, Baird AH. Broadcast spawning by *Pocillopora* species on the Great Barrier Reef. *PLoS One*. 2012; 7.

67. Riegl B. Hermatypic Coral Fauna of Subtropical Southeast Africa: A Checklist! *Pacific Sci*. 1996; 50: 404–414. Available: http://nsuworks.nova.edu/occ_facarticles

68. Smit K, Glassom D. Large fluctuations but constant mean temperatures allow corals to persist in intertidal rock pools on the east coast of South Africa Background. Smit Glas Helgol Mar Res. 2017; 71. https://doi.org/10.1186/s10152-017-0482-2

69. Glassom D, Celliers L, Schleyer MH. Coral recruitment patterns at Sodwana Bay, South Africa. *Coral Reefs*. 2006; 25: 485–492. https://doi.org/10.1007/s00338-006-0117-6

70. Schleyer MH, Celliers L. A Consideration of the Biodiversity and Future of Southern African Coral Reefs. *Coral Reefs*. 2000; 83–90.

71. Riegl B, Schleyer MH, Cook PJ, Branch GM. Structure of Africa's southernmost coral communities. *Bull Mar Sci*. 1995; 56: 676–691.

72. Celliers L, Schleyer MH. Coral community structure and risk assessment of high-latitude reefs at Sodwana Bay, South Africa. *Biodivers Conserv*. 2008; 17: 3097–3117. https://doi.org/10.1007/s10531-007-9271-6

73. Ridgway T, Riginos C, Davis J, Hoegh-guldberg O. Genetic connectivity patterns of *Pocillopora verrucosa* in southern African Marine Protected Areas. 2008; 354: 161–168. https://doi.org/10.3354/meps07245

74. Combosch DJ, Vollmer S V. Molecular Phylogenetics and Evolution Trans-Pacific RAD-Seq population genomics confirms introgressive hybridization in Eastern Pacific *Pocillopora* corals. *Mol Phylogenet Evol*. Elsevier Inc.; 2015;88. https://doi.org/10.1016/j.ympev.2015.01.005