Scleractinian corals (hereafter referred to as corals) are under threat in the Anthropocene due to the synergistic effects of stressors, natural and man-made 1. As we progress into uncertainty, including the likelihood of seawater temperatures breaching a 1.5 ºC threshold 2 and frequent and prolonged seawater temperature anomalies below and above threshold limits for coral survival 1,3–5, there is increased awareness and pressure for coral and coral reef conservation worldwide 7–12. While the conservation of corals is a top priority, planning for it can be confusing and overwhelming due to diverse species composition. With 1698 coral species in the order Scleractinia 13 distributed from the tropics to high-latitude non-reefal communities 14, working toward their protection is very difficult as they are morphologically plastic and challenging to identify genetically due to the lack of efficient molecular markers 15,16. Nevertheless, recent studies on coral taxonomy using integrative approaches (morphology, reproduction, and molecular) have gradually advanced our understanding of coral species 17–21.
The most dominant group, the genus Acropora, with approximately 400 nominal species being abundant in tropical and subtropical reefs and high-latitudinal non-reefal communities in the three ocean provinces, represents a typical showcase for demonstrating integrative approaches to understanding coral taxonomy 14,22. Traditional methods for describing species that include the use of morphological characteristics such as growth form, corallite, and corallum structure22,23, have not been effective in delineating species within Acropora as its morphological traits are often unstable and subject to variation, either interspecifically or interspecifically 24. For instance, skeletal characteristics such as size and shape of corallites can differ dramatically due to genetics, environmental factors 25, or even the location of corallites within the same colony14,22. 23 divided Acropora into 14 “species groups” based on skeletal morphology to assist in classification and subsequent research. 22 recognized 19 Indo-Pacific and one Caribbean species group within the genus. Nevertheless, only 122 of approximately 400 nominal species were considered valid in the latest revision of the genus 26, and the phylogenetic implications of Acropora species groups remain challenging 20,21,27–29.
Inferring Acropora species phylogeny using conventional single-locus markers has shown intertwined patterns of mitochondrial and nuclear gene trees since the late 20th Century 24,30–33. This has been widely interpreted as evidence supporting either ongoing hybridization between coral species 34,35 or incomplete lineage sorting due to recent diversification of species groups as well as in species with large effective population sizes 36. In addition, the incorrect identification of specimens, which is highly likely in morphologically diverse groups such as Acropora, could also result in polyphyletic patterns observed in gene trees 37,38. Nevertheless, cryptic lineages of Acropora morphospecies complexes that have a broad geographic distribution in the Indo-Pacific region have recently been identified by multi-locus (microsatellites) and single-nucleotide polymorphism (SNP) genetic data in combination with the genotypic cluster species definition 20,21,28,29,39–42. 39 use DNA sequence data from 12 genomic regions to clarify that A. hyacinthus collections from the Indo-Pacific region comprise at least four cryptic lineages. Two of the A. hyacinthus lineages have an unusual distributional pattern in East Asia, with one dominating marginal locations in Japan and Taiwan and the other dominating reefs of the Ryukyus Archipelago, Japan 40. Using microsatellites and single nucleotide polymorphisms (SNPs) derived from 2b-RAD sequencing, one lineage of A. hyacinthus identified in temperate Japan showed distinct bottleneck pressures, including higher clonality, increased linkage disequilibrium, and lower genetic diversity compared to subtropical populations 41,42. Combining morphological examination, genetic data of conventional single-locus markers, SNPs derived from the ultraconserved elements (UCEs) of hexacorallian genomes 27,43, and breeding trials successfully delineate distinct species boundaries among A. bifurcata, A. cytherea, and A. hyacinthus 28. The application of UCEs to examine species phylogeny in the A. tenuis clade derived from Cowman et al. 27 show this clade contains over 11 distinct lineages, only four of which correspond to then-accepted species based on morphological and geographical evidence, and allowed the description of two new species, A. rongi and A. tenuissima 20. Population structure and principal component analyses with SNPs (>60,000) indicate that A. cf bifurcata, A. cf cytherea, A. cf hyacinthus, and A. cf subulata are genetically distinct and do not show signs of introgression 21.
In this study, we apply molecular phylogenetic and morphological approaches to delineate species boundaries for Acropora divaricata 44 and A. solitaryensis 23 collected from coral ecosystems in East Asia. Based on morphology, both species have previously been categorized into the divaricata species group 14,22. The holotype of A. divaricata was described by Dana in 1846 from Fiji with nariform radial corallites characterized by large, open calices and possessing a reticulate coenosteum, despite variations in the density and arrangement of spinules (Table 1). At the colony level, A. divaricata exhibits an open caespitose-corymbose branching pattern, forming branching colonies with tapering branches that can curve and anastomose to create a network within the colony 22. Acropora solitaryensis is a relatively new species described by 23 with a similar growth form with respect to branching patterns, corallite, and corallum structures as A. divaricata, but with a tendency to fuse into solid plates along basal branches 22,23,45–47. Acropora solitaryensis is abundant in high-latitude coral ecosystems, including the low-latitude Flinders Reef (Moreton Bay), Middleton Reef, and Solitary Islands in Australia. However, specimens of A. solitaryensis has also been reported and collected from the tropical Murray Reef, Martha Ridgeway Reef, and Palam Islands in the Great Barrier Reef 23. This unusual geographic distribution and its morphological variabilities leads to a suspicion that A. solitaryensis might be readily divisible into five geographic subspecies that are widely separated spatially and environmentally 23. On the other hand, morphological plasticity also results in the sympatric and continuous occurrences of colony morphs, including arborescent (AR), solid plate (PL), and intermediate (IM) forms of A. solitaryensis, mainly in high-latitude coral ecosystems 14,22,46–48.
Both A. divaricata and A. solitaryensis have broad geographic distributions in the Pacific and Indian Oceans, with the former extending more to the Red Sea and the Persian/ Arabian Gulf 14,22,49,50. While 22 confirms the occurrence of A. solitaryensis in the central and west Indian Ocean, 14 confines its distribution to Sumatra and Indonesia in the east Indian Ocean. Both species are recorded from the island chain of East Asia, including Taiwan, Ryukyus Archipelago, and mainland Japan, with A. divaricata more in tropical and subtropical coral reefs and A. solitaryensis distributed further into high-latitude non-reefal coral ecosystems 14,22,46,48,51. Molecular phylogeny and cross-fertilization experiments was used to examine the relationships of the three morphs, namely arborescent (AR), intermediate (IM), and solid plate (PL), of A. solitarynensis in the high-latitude non-reefal region in Japan 47. Their results showed that AR is clearly distinct from PL, suggesting an absence of gene flow between the morphs. In their cross-fertilization experiments, gametic compatibility between AR and PL was extremely low, suggesting pre-zygotic isolation of these morphs. They concluded that AR and IM forms are variations of A. solitaryensis, whereas the PL form may be an undescribed species. Interestingly, the AR and IM types formed a monophyletic group in the molecular phylogenetic trees with samples of A. divaricata collected from the subtropical islands Ishigaki and Miyako of the Ryukyus Archipelago 47. In contrast, applying morphometric, molecular phylogenetic, and cross-fertilization experiments to examine the cryptic boundary of two morphs of A. divaricata, “slender” and “robust,” in the central Ryukyus, show that although inter-morphotype gamete compatibility was high in the year of overlapping spawning seasons for these two distinct morphotypes, population genetics analyses and molecular phylogenetic analysis show that they are genetically distinct and rarely hybridize 52. The molecular phylogenetic tree grouped the “slender” form into the clade containing mitochondrial control region DNA sequences of A. divaricata and the AR of A. solitaryensis. On the other hand, the “robust” form was grouped with the PL of A. solitaryensis from the high-latitude non-reefal region in Japan 47,52.
We hypothesized that the different growth forms of these two Acropora species collected from different geographical localities in the island groups of East Asia might represent two distinct lineages of the A. divaricata species group, one being A. divaricata and the other A. solitaryensis. To test this hypothesis, we collected A. solitaryensis from the holotype locality in the north Solitary Island, Australia, and extended the sampling of both species to the waters off Taiwan, the largest continental island in East Asia with distinct development of tropical reefs to the southeast and subtropical non-reefal coral communities to the northeast of mainland Taiwan and Penghu Archipelago in the Taiwan Strait 53 and high-latitude coral communities in Shikoku, Japan 54. By applying multi-locus molecular phylogenetic and morphological analyses, we confirm the species status of A. divaricata and A. solitaryensis in the island chain of East Asia. While both species have sympatric occurrences in the subtropical reefs of the Ryukyus Archipelago and high-latitude non-reefal coral communities in mainland Japan, they have discrete distributions following a boundary separating tropical coral reefs and subtropical non-reefal coral communities in Taiwan.