2bRAD reveals fine-scale genetic structuring among populations within the Mediterranean zoanthid Parazoanthus axinellae (Schmidt, 1862)

The zoanthid Parazoanthus axinellae (Schmidt, 1862) is a widespread coral species in the Mediterranean coralligenous assemblages where two morphotypes are found: Slender and Stocky, differing in size, color, and preferred substrate. Due to these marked differences, Slender and Stocky morphotypes were hypothesized to be two species. Here, we used 2bRAD to obtain genome‐wide genotyped single nucleotide polymorphisms (SNPs) to investigate the genetic differentiation between Slender and Stocky morphs, as well as their population structure. A total of 101 specimens of P. axinellae were sampled and genotyped from eight locations along the Italian coastline. In four locations, samples of the two morphotypes were collected in sympatry. 2bRAD genome-wide SNPs were used to assess the genetic divergence between the two morphotypes (1319 SNPs), and population connectivity patterns within Slender (1926 SNPs) and Stocky (1871 SNPs) morphotypes. Marked and consistent differentiation was detected between Slender and Stocky morphotypes. The widely distributed Slender morphotype showed higher population mixing patterns, while populations of the Stocky morphotype exhibited a stronger genetic structure at a regional scale. The strong genetic differentiation observed between P. axinellae Slender and Stocky morphotypes provides additional evidence that these morphs could be attributed to different species, although further morphological and ecological studies are required to validate this hypothesis. Our study highlights the importance of resolving phylogenetic and taxonomic disparities within taxonomically problematic groups, such as the P. axinellae species complex, when performing genetic connectivity studies for management and conservation purposes.