Phenotypic plasticity or speciation? A case from a clonal marine organism
BMC Evolutionary Biology volume 8, Article number: 47 (2008)
Clonal marine organisms exhibit high levels of morphological variation. Morphological differences may be a response to environmental factors but also they can be attributed to accumulated genetic differences due to disruption of gene flow among populations. In this study, we examined the extensive morphological variation (of 14 characters) in natural populations observed in the gorgonian Eunicea flexuosa, a widely distributed Caribbean octocoral. Eco-phenotypic and genetic effects were evaluated by reciprocal transplants of colonies inhabiting opposite ends of the depth gradient and analysis of population genetics of mitochondrial and nuclear genes, respectively.
Significant differences (P < 0.001) in 14 morphological traits were found among colonies inhabiting 12 locations distributed in seven reefs in southwest Puerto Rico. Results from principal component analysis indicated the presence of two groups based on depth distribution, suggesting the presence of two discrete morphotypes (i.e. shallow type < 5 m and deep type > 17 m). A discriminant function analysis based on a priori univariate and multivariate analyses (which separated the colonies in morphotypes) correctly classified 93% of the colonies for each environment. Light, water motion and sediment transport might influence the distribution of the two morphotypes. Reaction norms of morphological characters of colonies reciprocally transplanted showed gradual significant changes through the 15 months of transplantation. Sclerites of shallow water colonies became larger when transplanted to deeper environments and vice versa, but neither of the two transplanted groups overlapped with the residents' morphology. Genetic analysis of mitochondrial and nuclear genes suggested that such discrete morphology and non-overlapping phenotypic plasticity is correlated with the presence of two independent evolutionary lineages. The distribution of the lineages is non-random and may be related to adaptational responses of each lineage to the environmental demands of each habitat.
The extensive distribution and ample morphological variation of Eunicea flexuosa corresponds to two distinct genetic lineages with narrower distributions and more rigid phenotypic plasticity than the original description. The accepted description sensu Bayer (1961) of E. flexuosa is a complex of at least two distinct genetic lineages, adapted to different habitats and do not exchange genetic material despite living in sympatry. The present study highlights the importance of correctly defining species, because the unknowingly use of species complexes can overestimate geographical distribution, population abundance, and physiological tolerance.
The phenotype is considered the product of inherited genetic information and its interaction with the environment. Thus, differences in the phenotype can be explained by variations in environmental conditions, but also such phenotypic differences may reflect accumulated genetic variation due to disruption of gene flow between populations, and their subsequent speciation into biological species.
First, phenotypic plasticity enhances the survival and reproductive success of individuals by contributing to their ability to cope with environmental changes and to potentially adapt to new niches. Plasticity is an emergent property of the genotype and therefore also susceptible to natural selection . The change of the plastic response is often continuous, when the trait under analysis is subjected to an environmental gradient suspected to induce changes . The spectrum of phenotypes due to the environmental change describes the norms of reaction [2, 3]. Among the metazoans that exhibit the most extensive phenotypic plasticity are the marine modular species.
Phenotypic plasticity has been studied in algae , sponges , barnacles , gastropods [7, 8], bryozoans  and anthozoans [10–16]. This plasticity provides organisms with the ability to generate the fittest phenotype suiting local conditions. Morphology is then acquired through development under the current environment and can be changed in the next generation, if conditions are modified. Strong environmental gradients in the sea (e.g. light, water flow, sediment transport) may restrict the distribution of individuals to habitats, representing opposite ends of the gradient, where each phenotype is adapted [17, 18]. Furthermore, the fitness of the phenotypes varies along the environmental gradient . Disruptive selection may enhance the success of the two phenotypes at the opposite ends of the gradient by ecologically favoring each phenotype in its more suitable environment and by increasing genetic divergence. In this case, organisms settle and suffer high mortalities in non-optimal environments. Disruptive selection may be an influential evolutionary force leading to two disparate phenotypes by the existence of non-random mating related to habitat utilization .
In the absence of local adaptation, the high dispersal potential of marine propagules usually results in genetic homogeneity over large distances [20–22]. However, allopatric speciation is possible mainly because changes in oceanographic conditions, the emergence of land masses , and disconnection of populations by lower sea levels . As gene flow is disrupted by a geographic barrier, populations become isolated and diverge due to genetic drift. After genetic divergence has been acquired through generations of genetic drift and restricted gene flow, secondary contact can be achieved when the two new lineages attain similar geographic distributions .
Apart from allopatric divergence, sympatric divergence is also plausible. Speciation has occurred in spawning organisms with larvae capable of long dispersal [26–28] and genetic differences have been detected in sympatric populations [29, 30]. Ecological specializations to different habitats [31, 32], variable symbiotic relationships related to habitat distribution  and unsynchronized gamete release [34–36] may prevent organisms to reproduce randomly in sympatry, leading to a rapid evolution of mating systems [37–39] and eventually to speciation. It is not surprising that sibling species in the sea are more common than previously thought . Species with novel gene combinations can also be formed sympatrically through hybridization, an important process of diversification in marine and terrestrial systems [41–43].
In octocorals, phenotypic plasticity along environmental gradients or habitats is not uncommon [12, 14, 29, 44–47]. Octocorals are relatively abundant and visually dominant in low relief hard ground habitats with preference for high water motion areas [48, 49]. Light, water motion and sediment transport are determining factors in the distribution of gorgonians . These abiotic factors may induce morphological adjustments in broadly distributed species to optimize fitness under suboptimal conditions. Colonies of Eunicea flexuosa (Lamouroux 1821), in shallow forereef areas are susceptible to high water motion and are generally taller; grow in a single plane with thicker branches and bigger calices. In contrast, colonies in deeper environments are exposed to low water motion and less light. There, the colonies exhibit multiplane growth, are smaller with fewer terminal branches than their shallow counterparts with smaller and more sparse calices . Sclerite plasticity has also been correlated with differences in water motion and light [12, 14]. Smaller and thinner clubs and spindles are present in high water motion environments (i.e., forereef areas), providing a stronger structure and support to the colony . The high morphological variability in E. flexuosa could be due to phenotypic plasticity, genetic differentiation or a combined effect. In this study, first, the morphological variation of 14 traits of E. flexuosa was evaluated in seven reefs (from protected to exposed areas to water motion) at two depths (< 5 m and >17 m) and the correlation of morphology with light, water motion and sediment patterns was inferred. Second, environmental and/or genetic factors were studied to define the morphological variation of E. flexuosa. Reciprocal transplants of colonies inhabiting the opposite ends of the depth gradient were used to infer patterns, magnitude and direction of the phenotypic response. Gene genealogies of the mitochondrial gene msh1 and the nuclear gene 18S were used to elucidate possible genetic-phenotypic interactions. Lastly, allopatric and sympatric divergence (through ecological differentiation) and hybridization was considered as possible evolutionary processes that produce and maintain the morphological and genetic variation found in E. flexuosa.
Natural Variability of Morphological Traits
Mean values per trait and habitat are shown in Fig. 1. Measurements of morphological characters of most colonies in a given location fell within a narrow range around the mean. However, there were few colonies that did not match the mean population values. These atypical colonies were closer morphologically to colonies inhabiting the opposite depth habitat than to their neighbours and showed that in a given location, colonies under similar environmental conditions develop different morphologies. All traits varied significantly between reefs and among depths (Fig. 1). Values for CL, CW, CAL, CAW CD, SA, BT, CH, PD, TBN, and TBM decreased as depth increased, whereas values for SL, SW and ID increased with increasing depth (Fig. 1).
Depth differences were also recovered in the principal component analysis. The first principal component (PC1) explained 64% of the variance among measurements and was characterized by high negative loadings for SL, SW and ID and high positive loadings for CW, CAW, SA and BT (Table 1). PC2 explained 16% of the variance and was characterized mainly by positively weighted CL, CW and CAL and negatively by CD. Two-way ANOVA on the loading values of the first three principal components showed significant differences in reef and depth in PC1 (Table 2) and PC2 loadings showed significant differences in reef. Since 80% of the variation was explained by the first two principal components, environmental gradients related to depth, water motion, light and sediment transport could probably explain a portion of the morphological differences. As these co-vary with depth, it was not possible to quantify their individual effects on the morphological variation of E. flexuosa.
Results of the discriminant function analysis based on Wilk's λ (depth = 0.19, F = 56.29, P < 0.0001; reef λ = 0.31, F = 3.60, P < 0.0001 and zone λ = 0.21, F = 15.76, P < 0.0001) showed significant differences in colonies of different depths, reefs and zones. Depth differences as expected from the univariate and PCA analyses were clearly defined (Fig. 2) and the correct classification percentage was high (93%). The colonies in shallow habitats are concentrated on the right side of the canonical plot, while colonies in deep areas are mostly in the left side. There were a few deep water colonies misclassified and embedded within the shallow ones and vice versa (Fig. 2). Canonical plots of the DFAs of reef and zone showed high percentages of misclassification, 39% and 41%, however, a tendency from protected and deep to exposed and shallow was observed (data not shown).
The results suggest the presence of two morphotypes (Fig. 3). The first morphotype fits well with Bayer's description; adult colonies exhibit bushy-like shape and branch profusely in a single plane. Microscopically, size of leaf clubs, structural spindles and fused capstans are of ~200 μm, ~1000 μm and ~200 μm, respectively. The two morphotypes found in two depth habitats (hereafter deep and shallow indicating the depth of habitat) are exposed to different water motion, light, and sediment transport regimes. However, each morphotype was recorded in very low numbers at the opposite habitat. The shallow type was found more frequently in deep habitats than the deep type in shallow areas. Also, colonies sampled in the outer shelf reefs were morphologically closer to those in shallow areas of inshore and mid reefs; however the presence of both morphotypes was higher in such environments and each type was represented by at least six colonies (15 colonies total). In Weinberg, the presence of deep type colonies was associated with sand channels, areas characterized by high sediment transport. Along these channels deep type colonies were more frequently encountered.
On average, transplanted branches grew 1.94 ± 0.34 (1 SD) cm during the 15 months of the study. The new tissue deposited at the tip of the branches was enough to allow sampling of sclerites developed under novel conditions, thus the analysis excluded premature sclerites present at the very tip of the colony.
A two-way ANOVA test on linear growth values revealed a significant difference across depths (df = 1, F = 12.15, P = 0.001). Colonies in shallow environments regardless of population source (residents or transplanted) grew almost twice as fast as their deep counterparts. Population source and population X depth interaction were not significant (df = 1, 1; F = 0.067, 0.10; P = 0.80, 0.75; respectively). Of the 90 initial colonies transplanted, 59 were recovered for sclerite analysis. The mortality was not independent among groups (X2 = 10.449, df = 3, P < 0.025). The control colonies either from deep to deep, or shallow to shallow had higher survivorship (93% and 80%, respectively) than the transplanted ones. Mortality was highest (57% ex 17 of 30 colonies) in colonies transplanted from shallow to deep areas. Most of these colonies died by the 6th month; however neither detachment nor presence of predators was noticeable during the experiment. Competition with other reef organisms was also not evident. In most cases the tissue started to peel away until the entire axis was exposed. Colonies transplanted from deep to shallow areas had high survivorship (77% ex 23 of 30 colonies), nonetheless bleaching was recorded in 14 colonies (46%) and 23% of those died.
After 15 months of the experiment, measurements of sclerites for shallow and deep colonies were within the range suggested by the natural variability analyses (Fig. 4). However, spindle width and length of the reciprocally transplanted colonies gradually became similar to those of the resident colonies but never overlapped, colonies transplanted from shallow to deep areas tended to increase in sclerite size while spindles of colonies from deep to shallow became smaller. CL, CW, CAL and CAW showed similar trends but exhibited lower percent variation within groups, therefore resulting in non-significance between groups comparisons (i.e., DD, SS, SD and DS) (Fig. 4). Some colonies (i.e., SD16, SD18, SD19, DS7, DS8, DS14, DS 15, DS17, DS22, DS 23 and DD2) differed significantly from their source population either shallow or deep. Such atypical colonies did not show the gradual change in spicule size over time, observed within their transplanted group. Two-way repetitive measurement ANOVAs (Table 3) also showed significant differences among populations (DD, SS, SD and DS) for most traits (except CL and CAL) and the time X population interaction in SL and SW.
All colonies found in either shallow or deep areas were classified as deep habitat or shallow habitat. Morphology was used to define the second group of populations (deep or shallow types). The two classification schemes are not identical, as the atypical colonies found during the morphological analysis could be correctly classified. Also, each reef was treated as one population so that among reef comparisons could be established.
A total of 130 sequences of msh1 (723 bp) resulted in 10 distinct haplotypes, with three to nine haplotypes per population. The numbers of segregating sites were similar between populations, 9 and 8 sites were observed for colonies inhabiting shallow and deep environments, respectively.
The shallow type possessed 6 haplotypes and the deep type contained 4 haplotypes, among which the maximum difference was 6 substitutions. There were 2 and 3 segregating sites in the shallow and deep types, respectively. The most common haplotypes within each of shallow and deep types were represented by 54 and 30 individuals, respectively. The values of the nucleotide diversity indices (π, θ) for the pooled data were 0.0039 and 0.0024, respectively. π ranged from 0.0007 (shallow type) to 0.0039 (Romero reef), whereas θ was lowest in deep type (0.0009) and highest in Media Luna reef (0.0028). Most of the haplotypes were singletons; this mutation pattern was more common in colonies of shallow habitats. Fu's Fs test for the shallow and deep type revealed a significant departure from equilibrium only for the shallow type (-3.05038, P = 0.037). The excess of rare mutations observed in the shallow type is consistent with population expansion or purifying selection. Tajima's D tests were not significant.
The within shallow and deep type divergence varied from 0% to 0.3%, as estimated with the Kimura-2-parameter model . Divergence between shallow and deep types ranged from 0.61% to 1.07%. The shallow type was more closely related (0.61% to 0.76%) to Plexaura homomalla than the deep type (0.92% to 1.23%).
A 251 bp fragment of 18S was sequenced from 143 colonies of E. flexuosa. Most of the sequences (90) were from Media Luna reef (colonies used in the transplant experiment), while 32 sequences were from Romero and an additional 21 from Culebra. Among all sequences, 15 haplotypes were discernable with colonies in shallow habitats presenting the highest haplotype diversity (h = 11). Among reefs, Media Luna had the highest number of haplotypes (h = 11). The values of the nucleotide diversity indices (π, θ) for the pooled data were 0.0088 and 0.0079, respectively. π ranged from 0.0014 (deep type) to 0.0093 (Romero reef), whereas the lowest (0.0044) and the highest θ (0.0089) was also observed in the same data partitions. The within shallow type genetic variability ranged from 0% to 0.97%, and the within deep type variability ranged from 0 to 0.48%, as estimated by the Kimura-2-parameter model. Divergence between shallow and deep types was similar to the msh1 divergence and ranged from 0.97% to 1.95%. Fu's Fs test for the data partitions revealed a significant departure from equilibrium for the shallow type (-3.7390, P = 0.0459) and the deep type (-4.4450, P = 0.0036). The excess of rare mutations observed in the ribosomal gene of the shallow type may be explained by population expansion or purifying selection. Tajima's D was significant only for the deep type (-1.6119, P = 0.0224). Departure from equilibrium was consistently observed across the two genes only for the shallow type group. The significant departure from neutrality recorded in the deep type by the 18S gene and not the msh1 gene, may be attributed to the fact that the two genes belong to different genomes and are under different selection and stochastic processes.
AMOVA tests showed significant differentiation among populations (Table 4), regardless whether the populations were divided per habitat (shallow or deep) or by morphotype (shallow type and deep type). Nonetheless, population differentiation was maximized when the latter assignment was used. A comparison among reefs (Media Luna and Romero) yielded non-significant FST values, suggesting an extensive gene flow among these two reefs. Also, a comparison of all shallow habitats (Romero shallow, Media Luna shallow and Culebra) suggested homogeneity among populations. The genetic homogeneity within La Parguera is not surprising because of the short distance between the reefs. Even when including colonies from Culebra, which is 100 km northeast from La Parguera, genetic homogeneity was still observed. However, when comparisons of reefs lacking the deep habitats (e.g., in Culebra only shallow areas were sampled), FST values were significant, suggesting population subdivision. As inferred from the AMOVAs the haplotype network analysis showed similar patterns (Fig. 5). Two groups were extracted from the haplotype network analysis, the first consists of the vast majority of sequences that have the shallow morphotype and the second group contains most of the deep morphotype.
Gene genealogies were constructed in PAUP using ML with the HKY and Jukes Cantor as most suitable substitution models for msh1 and 18S, respectively. Analysis performed using neighbour joining and parsimony yielded similar patterns. Also, the topology of msh1 and 18S was similar, however the 18S analysis recovered only one of the clades, as 18S was less variable than msh1.
The genealogy divided the individuals in two clades (Fig. 6), which are highly indicative of the habitat origin (shallow or deep). In the bottom clade, 62 of the colonies are from deep areas; however nine colonies fell within this clade. These nine colonies are the atypical colonies found during the morphological analysis (e.g., Media Luna shallow), where comparisons of morphological characters of these atypical colonies did not match the overall population mean, despite living in the same habitat. Similar patterns are displayed for the top clade, which represents 68 colonies sampled in shallow areas and three atypical colonies found in deep areas. Furthermore, these atypical colonies in each clade are those colonies that were misclassified (7%) by the discriminant function analysis. Table 5 shows a summary of the number of individuals within each genetic lineage divided by habitat and morphotype (shallow and deep types). The atypical colonies within each habitat were later shown by the DNA analysis (Fig. 6) to be phylogenetically closer to their opposite habitat congeners, reinforcing the morphological differences.
Eunicea flexuosa is divided into two discrete morphological forms in southwest Puerto Rico. The shallow type is pervasive in shallow areas, but a few colonies of the shallow type can be found in deep habitats. The second morphotype (deep type) could be described as small colonies, with fewer terminal branches, more sparsely polyps, thinner branches and bigger spindles. The deep type is found at the muddy bed at the base of forereef of inside and mid reefs and largely confined to deep areas with low water motion patterns, high sediment transport and lower light levels. The observed variation in E. flexuosa resembles findings on phenotypic plasticity of other modular organisms. Mechanical stimuli such as wind speed or water motion induce morphological adjustments in both terrestrial and marine plants [4, 51]. Other marine modular taxa such as the bryozoan Membranipora membranacea adapt to water flow to maximize food capture , the demosponge Halichondria panacea develop stiffer branches in high energy habitats  and the coral Madracis mirabilis exhibit sparser and thicker branches in high flow areas . Light has also been associated with phenotypic changes in plants , algae  and corals . Other non-mechanical stimuli such as presence of predators in the bryozoan M. membranacea can induce plastic responses (i.e. defensive spines) .
The discrete morphological distribution found in E. flexuosa is in concordance with previous studies of octocorals, where colonies inhabiting deep forereef areas were thinner with more sparse and fewer calices and have bigger spindles than the colonies in shallow habitats in either back or forereef areas in the Florida Keys . Spindles at the branch tip were significantly smaller and presumably underdeveloped in E. flexuosa, while in Briareum asbestinum, West  reported bigger spindles at the tips. The discrepancy of the findings may be related to different ecological pressures on the two species, different function of spindles in each species (i.e., B. asbestinum lacks a central axis) or may depict difference responses constrained by their separate phylogenetic history. Moreover in terms of habitat differences, smaller spindles at high water motion areas are thought to increase stiffness of the colony to avoid breakage. Grigg , West et al.  and Kim et al.  reported that sclerite size is negatively correlated with water motion for B. asbestinum and E. flexuosa. However, the opposite observation has been made in Eunicella singularis . Colony size, calice diameter and polyp density may compensate water flow changes and enhance respiration, feeding or structure as has been reported for scleractinian corals [10, 57, 58]. Furthermore, the overall discrepancy in colony development (more and bigger terminal branches) may be interpreted as a response to the water drag forces acting upon the colony, yet colonies in shallow areas are exposed to higher water flow and vice versa.
In E. flexuosa, the differences of the two morphotypes are related to depth profiles, water motion and sediment transport. However, as water motion, sediment transport and light co-vary along depth profiles, the study could not distinguish the individual influences. Furthermore, the deep habitats of the inside and mid reefs differ significantly from those of the outer reefs. In the inside and mid reefs, loose sediments move along the slope and are deposited at the base of the forereef, where deep colonies were sampled. Therefore, settling larvae in such areas are subject to higher sediment deposition and probably have developed an adaptational response to such environmental condition. Similarly, colonies in the two outer reefs, El Hoyo and Weinberg, developed in the edges of sand patches were different than those spread along the hard ground bottom. In La Parguera, the octocoral community observed at the sand bed adjacent to the forereef (~20 m), is less diverse than the community of hard ground habitats. In such suboptimal habitats, Pseudopterogorgia spp., Plexaurella spp. and some Muricea spp. (except Muricia muricata) are the dominant species. These species are probably better adapted to overcome sediment suffocation and develop. Pseudopterogorgia spp. and Muricea spp. avoid sediment burial by rapidly reaching a safe size due to the lack of thick central axis which allows for faster growth. In Plexaurella spp., the presence of big polyps can efficiently remove sediment from the colony. The deep type of E. flexuosa is a thin colony with a small central axis and few polyps, which might grow fast enough to overcome sediment overload.
The morphology of E. flexuosa can also be affected by light and other factors associated with depth such as food resources, presence of predators and hydrostatic pressure. Nonetheless, water motion covaries with depth and therefore both are related. Light is likely a major factor related to depth and affects anthozoan morphology, as presence of zooxanthellae is tightly related to carbonate uptake and energy supply for the host. Light levels influence the phenotype of anthozoans due to changes in zooxanthellae concentration and decrease of calcification rates [14, 55, 59, 60]. Moreover, plate-like colonies most likely represent an evolutionary response to compensate for low levels of light by increasing area and zooxanthellae concentrations.
The two morphologies of E. flexuosa associated with different habitats showed some degree of phenotypic plasticity in sclerite characters (especially spindles), which showed a clear tendency to increase and decrease when transplanted to deep and shallow areas, respectively. Nonetheless, after the 15 months of the experiment in neither case the spindles of transplanted colonies became similar in size as the residents' spindles. The results suggest that either there was not enough time for the colonies to produce new tissue under the novel conditions or there is a non-environmental factor accounting for the rest of the morphological variation. There is less evidence for the former; colonies grew on average ~2 cm during the transplant experiment, and since the analysis was performed using tissue 1 cm from the branch tip, there was more than 0.5 cm tissue to be analyzed. Also, the reaction norm graphs showed a stabilization of the curve at 10 and 15 months, suggesting that the maximum of variation was reached.
Gorgonian corals resemble plants in many ways and morphologically they are modular organisms where integration/disintegration can occur in response to environmental gradients . If phenotypic plasticity in gorgonians is expressed at the modular level (polyp and branch), the transplant results may be interpreted as the sum of all historical (branches before transplantation) and experimental (branches during transplantation) interactive and antagonistic changes between modules. The observed changes in sclerite size which are changes affecting individual modules suggest there is some plasticity at the modular level. The plastic response of modules has cascade effects on the colony level. For example, when colonies are transplanted from deep to shallow areas, sclerites become smaller and denser enhancing the flexibility of the colony.
Other transplant experiments with modular organisms (e.g., the sponge Halichondria panacea) exhibit variable responses to environmental stimuli . When colonies of H. panacea were placed in high energy areas, they became stiffer, a response which is in agreement with our observations on E. flexuosa. Alternatively, plasticity may be most likely promoted during the early life stages rather than in mature colonies. Other stochastic environmental factors such as the 2005 Caribbean wide bleaching event may have affected the results.
The possibility that no environmental factors account for the majority of morphological variation is supported by the genetic data. Analysis of the msh1 and 18S genes resulted in two major lineages associated with depth. Colonies transplanted from shallow to deep habitats suffered about 50% mortality; those transplanted from deep to shallow suffered 20% mortality. Since survivorship is a fundamental attribute of fitness then increase in mortality would indicate lower fitness. Therefore, each morphotype is better adapted to the deep or shallow areas. The results from the transplant experiments may represent experimental error or a drastic response of the colonies to novel environmental conditions. Previous studies have shown that adult colonies (> 20 cm) in natural populations are stable and have a normal survivorship rates above 90% , suggesting that such high mortalities are natural responses rather than sampling error. Although the two genetic lineages are associated with the two habitats, there is a noticeable response of both shallow and deep type to environmental stimuli. As previous studies have reported, there is a positive correlation of spindles to increase in size as depth increases and vice versa [12, 14]. Clubs and capstans slightly changed with depth but the tendency was not consistent.
In Puerto Rico, the morphological divergence found in colonies of E. flexuosa is genetically based. Gene genealogies, haplotype networks and AMOVA analysis of both nuclear and mitochondrial genes suggested that such discrete morphological distribution is correlated with the presence of two distinct lineages, distributed non-randomly in shallow and deep environments. However, either of the two lineages can infrequently be found in both depth habitats. The genetic break found in nuclear and mitochondrial DNA suggests that gene flow ceased a long time ago and divergence may have led to speciation. Fixed differences in both nuclear and mitochondrial genes are comparable to those reported between species of octocorals [63, 64]. Therefore, the current description sensu Bayer  of E. flexuosa is a complex of at least two distinct genetic lineages, adapted to different habitats and that do not exchange genetic material despite living in sympatry. The extensive distribution and ample morphological variation corresponds to two distinct genetic lineages with narrower distributions and more rigid phenotypic plasticity. The observed genetic pattern may have resulted from 1) secondary contact after populations diverged in allopatry and reproductive incompatibility developed, 2) by divergence with gene flow through ecological specialization in sympatry or 3) by the poorly understood process of hybridization in anthozoan evolution.
Processes of Divergence
The planulae of most anthozoan broadcasters such as E. flexuosa are capable of staying days to weeks in the water column producing genetically homogeneous populations across large geographic areas [21, 22, 65]. Despite the high potential for dispersion and population connectivity, allopatric speciation is probably the most common mode of speciation in marine environments [23, 24, 66, 67]. Recent genetic studies [68–71] have contradicted earlier assumptions of population homogeneity in Caribbean populations of marine taxa [72, 73]. Even though allopatric distribution could have caused the divergence of the lineages, there is no recent evidence of geological processes that may have altered the patterns of Caribbean circulation. However, distinct sympatric lineages have been uncovered in other Caribbean invertebrates .
On the other hand, sympatric speciation by ecological differentiation  and disruption of gene flow in proximate populations is plausible. The two genetic lineages of E. flexuosa are found at the opposite ends of a depth gradient. Diversifying selection may favor the two phenotypes at the extremes of the depth gradient, preventing gene flow through assortative mating and eventually leading to new species . These ecological specializations to depth habitats have been pointed out in earlier reviews . Such ecological differences in niche utilization may be reinforced by dissimilar characteristics associated with the habitats they occupy. Symbiotic relationships to host  or to environment , differential timing of gamete release due to depth related differences [34, 35] may have provided different resources to populations at different habitats and eventually prevent random mating. As a consequence, rapid evolution of mating systems may have been favoured [37–39]. In Caribbean corals, diversifying selection has been proposed in at least two species: Montastraea annularis  and Favia fragum . It is likely that the genetic differences reported by Brazeau and Harvell  in the gorgonian B. asbestinum could have arisen through the same mechanism. Although, disruptive selection seems to explain the divergence of the two lineages of E. flexuosa, hybridization is also another plausible mechanism to cause diversification in marine taxa.
Hybridization is a common phenomenon in plants and the rise of new lineages due to reticulations is often reported in flowering plants [41, 78]. In the marine environment, there are instances of hybridization in angelfishes , cichlids , blue mussels  and corals [82–84], suggesting that hybridization is an important evolutionary mechanism for speciation. Nonetheless, it is often assumed that discrepancies in gene phylogenies or FST statistics from different molecular markers is interpreted as incomplete lineage sorting, rather than reticulations as it is interpreted often in plants [41, 78]. The phenomenon of reticulate evolution may have great influence in the speciation of marine species especially those living in sympatry with high potential for hybridization (e.g., spawners).
Veron  has provided a theoretical framework to consider reticulate evolution as an important factor of coral evolution. Direct measures of chromosome differences established in Acropora , genetic surveys of the nuclear genome of corals [83, 84, 86] and direct crosses of gametes  have shown introgression in natural populations. In E. flexuosa, the two uncovered lineages may have arisen by hybridization between the common form of E. flexuosa with another Eunicea or Plexaura species. Furthermore, hybrid fitness may increase over parent fitness in novel environment or in extreme habitats. Hybrids tend to explore novel habitats avoiding introgression and competition with their parents [87, 88]. Therefore, it is likely that the lineage related to deep muddy areas is of hybrid origin and is better adapted to such conditions than the parental species.
The accepted description sensu Bayer  of E. flexuosa is a complex of at least two distinct genetic lineages, adapted to different habitats and that do not exchange genetic material despite living in sympatry. The present study highlights the importance of correctly defining species, because the unknowingly use of species complexes can overestimate geographical distribution, population abundance, and physiological tolerance. The non-random distribution of both morphotypes can yield misleading population genetic inferences in the absence of adequate taxonomy. Consequently, decisions based on these estimates will have repercussions in conservation programs . Detailed studies of the mechanism by which anthozoans achieve assortative mating and become reproductively isolated would give us insights in the speciation process. Also, cross fertilization experiments, genetic assessment of shared alleles through genetic markers and karyotyping may shed light on the speciation process via hybridization. Reticulations are common in plants, a group that resembles most of the ecological aspects (bet hedging strategies, modular organization, philopatric recruitment, etc.) that govern marine modular organisms.
Formerly known as Plexaura flexuosa, the species has been recently placed into the genus Eunicea based on molecular data and the size ratio of spindles and clubs . Eunicea flexuosa is an octocoral cnidarian forming colonies of ~1 m in height found in coral reefs and tropical rocky walls [46, 49, 90]. Eunicea flexuosa is relatively abundant in low relief hard ground habitats with preference for high water motion areas [48, 49]. The distribution of the species spans through several environmental gradients (e.g., depth, light, water motion and sediment transport)  and display an unusual amount of morphological variability. Eunicea flexuosa is a gonochoric gorgonian that reproduces sexually by spawning gametes . Asexual reproduction through fragmentation can also take place when loose branches spread to the surroundings (< 10 m) and new colonies develop from clone-mate propagules but is not as common as in P. kuna . The colony is produced by the asexual budding of its polyps, generating a branching dichotomous morphology, arranged in an axial skeleton of fused sclerites of calcium carbonate . Adult colonies exhibited bush-like shape and branch profusely in a single plane. Microscopically, the apertures (calices) present an inconspicuous lower lip with an unarmed collaret. Sclerites are arranged in three layers. The axial sheath is made up of fused capstans usually purple and ≤ 200 μm in length. The external layer contains leaf clubs of ~200 μm with 3 or 4 serrate folia and structural spindles (~2000 μm in length) are disposed in a mid layer . The last two features distinguish the species from the other plexaurids. For both P. homomalla and P. nina the leaf clubs and spindles are smaller (150 μm and 700 μm, respectively), spindles are also more slender and calyces usually exerted in P. nina (Bayer 1961). Contrary to P. homomalla, colonies of E. flexuosa in the field usually tend to be branched in one plane but not in a net-like shape with calyces as a lower lip. All colonies sampled at different depths and reefs fell within the original description of Eunicea flexuosa (sensu Bayer 1961, modified by Grajales et al. 2007). None of the colonies exhibited prominent and well developed calyces, typical of other Eunicea species.
The study was carried out in La Parguera southwest Puerto Rico (Table 6) during September 2004 to August 2006. The hydrography of the area has been studied and described elsewhere . Reefs off La Parguera are exposed to wave action, generated by the easterly trade winds. The shallow (< 10 m) benthic community is visually dominated by octocorals, forming colourful gardens, especially in shallow reef platforms with high water motion and wave energy. The current abundance of gorgonians may have been enhanced by the die off of the herbivorous sea urchin Diadema antillarum . La Parguera is composed by three prominent reef formations, located parallel to the shoreline: 1) inshore reefs which are more protected to waves and currents, but subjected to higher sedimentation rates, direct contact with sewage discharges and lower scleractinian cover, 2) mid reefs and 3) outer reefs which are fully exposed to wave energy, with higher coral cover and significant bottom relief.
Gorgonian colonies were sampled in seven reefs for morphological measurements (Table 6). Three inshore, protected reefs (Conservas, Pelotas and Romeo); two mid-shelf reefs exposed to wave action (Media Luna and Turrumote) and two outer-shelf reefs (El Hoyo and Weinberg) were included. In each site, morphological variability of E. flexuosa was assessed at two depths (shallow < 5 m and deep > 17 m), except in the outer reefs El Hoyo and Weinberg, which are at depth 23–27 m. Therefore, in most reefs there are two depths, except the outer reefs. In each location, 15 colonies were sampled (n = 30/reef).
Natural Variability of Morphological Traits
Three macro-morphological traits: colony height (CH), branch thickness (BT) and branch development (BD) were measured in 180 colonies, representing fifteen colonies per location (7 reefs per 2 depths, except El Hoyo and Weinberg which represent deep habitats). Branch development was assessed as the total number of terminal branches (Fig. 7b) and the average length of 10 haphazardly chosen branches from branch tip to the first node in terminal branches (Fig. 7b).
Nine microscopic characters from the same 180 colonies were measured at 2 cm off the tip to avoid measuring underdeveloped sclerites because new tissue is generated at the tips and contains higher number of smaller sclerites, especially in the spindle width and length . Two cm off the branch tip was far enough to avoid smaller, underdeveloped sclerites, but within the range in which the new tissue would be generated under the novel environmental conditions during the reciprocal transplant experiment. Polyp density (PD), calice diameter (CD), inter-calice distance (ID) and width and length of external clubs (CW and CL) mid-layer spindles (SW and SL) and axial sheath capstans (CAW and CAL) were selected for the analysis (Fig. 7). Polyp density was estimated by counting the number of polyps/cm2 and standardized with branch thickness, assuming cylindrical shape of the branches. In these micro-morphological traits, 20 measurements (except polyp density) were performed in randomly selected calices in CD and ID or sclerites in CW, CL, SW, SL, CAW and CAL; representing 300 measurements per depth, 600 per reef, 3,600 per character and 28,800 in total, excluding polyp density measures. Octocoral branches were collected by clipping off a 5 cm section at the branch tip, slightly bleached with Clorox (5%) to remove some tissue, rinsed in distilled water, and dried. The slight bleaching kept the colony shape intact, avoiding colony dissolution while allowing clear observations of the calices. For sclerite analysis 1 cm section at 2 cm from the tip was collected and dissolved with Clorox (5%), following Bayer's protocol . For clubs and capstans the samples were taken by placing the spicules in slides and random samples were obtained by moving blindly the slide and measuring all sclerites in each new visual field until 20 sclerites were measured. Spindle, calice and polyp were analyzed by photographing the characters with an Olympus BX-51 compound microscope. All measurements were carried out using photographs (calibrated with a slide of 10 μm accuracy) taken with an Olympus C-5050 camera system attached to an Olympus SZH-10 stereo microscope. Analysis of the photos was performed using SigmaScan (SPSS Inc.).
The extent of variation among habitats (i.e., PS, PD) in morphological traits (14) was analyzed by individual One-way ANOVAs. Principal component analysis (PCA) was applied to test the association of the traits and visually examine overall trends, the analysis was done with both raw and standardized (in a 0 to 1 scale) data. The PCA scores of the first three principal components were used to test for differences among groups and factor interactions in a two-way ANOVA. The two factors were degree of protection (zone) and depth (< 5 m and < 17 m). Reefs were nested within zones and depth was factorial with reefs. Kolmogorov-Smirnov and Levene's test were used to test for normality and homogeneity of variance, respectively . Discriminant function analysis (DFA) of 14 morphological characters (CH, BT, BD, PD, CD, ID, CW, CL, SW, SL, CAW and CAL, Fig. 7) was used to explore if colonies could be assigned independently to one of the two resulting morphotypes and replicate the same predetermined groups in the sampled populations (i.e., inner zone, mid zone, shallow and deep). Wilk's λ was used in the DFA to test for multivariate differences among groups . Individual DFA's were employed for depth, location and zone.
If reciprocally transplanted colonies become similar in sclerite size to the neighbouring colonies, the environment is likely to have an effect on the phenotype. As gorgonian growth is slow (< 3 cm a year; Yoshioka unpub. data) only sclerite related traits were evaluated. We examined the subapical parts of branches of transplanted gorgonians because this is where new spicules and skeleton material are generated. By avoiding the sampling of tips no underdeveloped sclerites were counted .
To test for eco-phenotypic responses of the species to different depth habitats, thirty branches (< 30 cm) of different colonies separated by >10 m from each other, from each depth were reciprocally transplanted in Media Luna reef. Additionally, 15 branches from each of the two depths were auto-transplanted (transplanted to the same depth) to serve as controls. In each depth, 45 branches were transplanted per depth. Fifteen of them were residents and 30 were introduced from either the shallow or the deep habitat, for a total of 90 colonies in the experiment. Each branch was glued into a short piece of plastic PVC pipe with marine hydraulic cement. The cement did not affect the octocoral growth, as live tissue quickly overgrew the dried cement; the same material has been successfully used in restoration of gorgonian populations (Yoshioka et al. unpub. data). Each of the plastic pipes was attached to a 1.5 × 1 m cement panel. Fifteen branches were randomly assigned in each cement panel (allowing enough distance to avoid conspecific aggressions) and three panels were placed per depth (45 colonies per depth) separated by 5 m from each other. Each branch was assigned a code and individually tagged on the panels. The same code was also used for genetic analysis. The experiments were monitored every month to remove fouling organisms and sediment from the panels. Three times during the 15 months of the study (once every five months) surveys were carried out and linear growth and survivorship recorded. As the transplanted colonies were introduced to the novel habitat, the survivorship and linear growth were used to evaluate the overall colony response in each treatment. Linear growth was measured from base to tip of the colony and survivorship was recorded as live or dead. Branch thickness and polyp density could also reflect phenotypic plasticity but these parameters were not measured because detectable variation requires a longer experimental time due to the slow growth (~3 cm y-1) of E. flexuosa.
Linear growth of the transplanted colonies was assessed by subtracting the initial length from the final length. A two-way ANOVA was performed to test for differences in population (colonies from shallow and deep areas) and depth (colonies grew in shallow and deep environments) and their interactions. Both factors were considered fixed effects. Chi-square analysis on the number of dead colonies was used to test for survivorship independence among groups. To check for environmental effects on the groups, morphological trends were visually examined by graphing the reaction norm through time per group per trait. Repetitive measurement ANOVAs were used to test the effect of population and time and their interaction. ANOVA assumptions were evaluated as previously described. All statistical analyses were performed in JMP ver 5.0.1, InfoStat ver. 2004, and SigmaStat.
Portions of msh1 and the 18S genes were sequenced to test for genetic differentiation among populations of E. flexuosa living in different depths and reefs (Media Luna, Romero and Isla Culebra) in Puerto Rico. Msh1 is a mitochondrial gene unique to octocorals that codes for a DNA mismatch repair protein and is not present in other Cnidaria . Msh1 and the nuclear 18S provide enough resolution to discriminate between closely related species of octocorals [63, 64]. Segments of msh1 and 18S genes of all transplanted colonies in the eco-phenotypic experiment were examined, 90 colonies in total (45 individuals for each depth habitat). Additional samples from Romero and Culebra were included in the analysis to test the reproducibility of the preliminary results found in Media Luna. Thirty one colonies from Romero and 21 colonies from Culebra were included for the msh1 analysis. Thirty two colonies from Romero and 21 from Culebra were included for the 18S analysis. A 3 to 5 cm colony fragment was brought to the laboratory for immediate DNA extraction or stored in 95% ethanol for subsequent work. The PureGene DNA isolation kit (Gentra) was used for DNA extraction. The PCR amplifications was performed in an Eppendorf MasterCycler with the same cycling conditions for both genes, consisted of an initial denaturation at 95°C for 3 min, followed by a touch-down routine of annealing of 10 cycles at 43°C of 45 sec and 25 cycles at 48°C of 45 sec; denaturation at 94°C for 45 sec and elongation at 72°C for 5 min. The primers used to amplify 18S (A18S – 5'-GATCGAACGGTTTAGTGAGG-3' and ITS-4 – 5'-TCCTCCGCTTATTGATATGC-3') and msh1 were developed by Takabayashi  and France and Hoover , respectively. Sequencing reactions were prepared with a DYEnamic ET Terminator Cycle Sequencing Kit (GE) and loaded in a MEGABase 96 lane Sequencer for capillary electrophoresis.
DNA sequencing trace files were imported into the Phrap/Phred/Consed programs  for base calling, quality assessment, assembly and visualization. Mutations were verified in both the forward and reverse direction. Sequences were then aligned in MacClade  and compared with BLAST to publicly available sequences of closely related gorgonians. The haplotype (h) and the nucleotide diversity (π), number of segregating sites (S), the Watterson's estimator (θ w ) were evaluated according to Nei  as implemented in DNAsp 4.0 . Tajima's D  and Fu's Fs  were used in ARLEQUIN  to test for deviations from neutrality. A parsimony haplotype network was constructed for the msh1 sequences using the Templeton et al.  algorithm as implemented in TCS version 1.21 . The parsimony network was constructed with confidence level set at 95%. Analysis of molecular variation [AMOVA, ] of among reefs, between habitats and between morphotypes was performed in ARLEQUIN. The AMOVAs were performed with 10,000 permutations by using conventional F-statistics with haplotype frequencies.
Gene genealogies were constructed for msh1 and 18S using the maximum likelihood (ML) method in PAUP* . Hierarchical likelihood ratio tests in MODELTEST 3.6  suggested that the HKY  and JC models were the best substitution models for msh1 and 18S, respectively. For the heuristic searches in ML, data were bootstrapped 100 times, and sequences were added randomly ten times. Phylogenetic relationships were also constructed using neighbor-joining and maximum parsimony with 1000 bootstrap replicates as implemented in PAUP*. Given the uncertain phylogenetic position of E. flexuosa and its conspecifics , several outgroups were used, including P. homomalla, Muricea muricata and Eunicea spp. Regardless the outgroup, identical topologies were obtained. All sequences have been deposited in GenBank (Accession numbers EF659469-EF659598 and EF659599- EF659741)
Pigliucci M: Evolution of phenotypic plasticity: where are we going now? . Trends Ecol Evol. 2005, 20 (9): 481-486. 10.1016/j.tree.2005.06.001.
Futuyma DJ: Evolutionary Biology. 1998, Sunderland, MA. , Sinauer
Schlichting CD: The evolution of phenotypic plasticity in plants. Annu Rev Ecol Syst. 1986, 17: 667-693. 10.1146/annurev.es.17.110186.003315.
Fowler-Walker MJ, Wernberg T, Connell SD: Differences in kelp morphology between wave sheltered and exposed localities: morphologically plastic or fixed traits?. Mar Biol. 2006, 148: 755-767. 10.1007/s00227-005-0125-z.
Palumbi SR: Tactics of acclimation: morphological changes of sponges in an unpredictable environment. Science. 1984, 225: 1478-1480. 10.1126/science.225.4669.1478.
Marchinko KB: Dramatic phenotypic plasticity in barnacle legs (Balanus glandula Darwin): Magnitude, age-dependence and speed of response. Evolution. 2003, 57: 1281-1290.
Trussell GC: Phenotypic plasticity in an intertidal snail: the role of a common crab predator. Evolution. 1996, 50: 448-454. 10.2307/2410815.
Vermeij GJ: Phenotypic evolution in a poorly dispersing snail after arrival of a predator. Nature. 1982, 299: 349-350. 10.1038/299349a0.
Yoshioka PM: Predator-Induced Polymorphism in the Bryozoan Membranipora membranacea. J Exp Mar Biol Ecol. 1982, 61: 233-242. 10.1016/0022-0981(82)90071-5.
Bruno JF, Edmunds PJ: Clonal variation for phenotypic plasticity in the coral Madracis mirabilis. Ecology. 1997, 78: 2177-2190.
Gleason DF: The adaptive significance of morphological plasticity in the reef coral Porites astreoides. Am Zool. 1992, 32: 92D-
Kim E, Lasker HR, Coffroth MA, Kim K: Morphological and genetic variation across reef habitats in a broadcast-spawning octocoral. Hydrobiologia. 2004, 530/531: 423–432-10.1007/s10750-004-2646-8.
Takabayashi M, Hoegh-Guldberg O: Ecological and physiological differences between two colour morphs of the coral Pocillopora damicornis. Mar Biol. 1995, 123: 705-714. 10.1007/BF00349113.
West JM, Harvell CD, Walls AM: Morphological plasticity in a gorgonian coral (Briareum asbestinum) over a depth cline. Mar Ecol Prog Ser. 1993, 94: 61-69. 10.3354/meps094061.
Willis BL, Ayre DJ: Asexual reproduction and genetic determination of growth form in the coral Pavona cactus: biochemical genetic and immunogenetic evidence. Oecologia. 1985, 65: 516-525. 10.1007/BF00379666.
Grajales A, Aguilar C, Sánchez JA: Phylogenetic reconstruction using secondary structures of Internal Transcribed Spacer 2 (ITS2, rDNA): finding the molecular and morphological gap in Caribbean gorgonian corals. BMC Evolutionary Biology. 2007, 7: 90-10.1186/1471-2148-7-90.
Doebeli M, Dieckmann U: Speciation along environmental gradients. Nature. 2003, 421: 259-264. 10.1038/nature01274.
Hilbish TJ: Demographic and temporal structure of an allele frequency cline in the mussel Mytilus edulis. Mar Biol. 1985, 86: 163-171. 10.1007/BF00399023.
Doebeli M, Dieckmann U, Metz JAJ, Tautz D: What we have also learned: adaptive speciation is theoretically plausible. Evolution. 2005, 59: 691-695.
Grosberg RK, Cunningham CW: Genetic structure in the sea: from populations to communities. Marine Community Ecology. Edited by: M. D. Bertness, S. Gaines, Hay ME. 2001, Sunderland, MA , Sinauer, 61-84.
McFadden CS, R. K. Grosberg, B. B. Cameron, D. P. Karlton, D. Secord: Genetic relationships within and between solitary and clonal forms of the sea anemone Anthopleura elegantissima revisited: Evidence for the existence of two species. Mar Biol. 1997, 128: 127-139. 10.1007/s002270050076.
Lessios HA, Robertson DR: Crossing the impassable: genetic connections in 20 reef fishes across the Eastern Pacific Barrier. Proc R Soc Lond B. 2006, 273: 2201-2208. 10.1098/rspb.2006.3543.
Knowlton N, Mills DEK: The systematic importance of color and color Pattern: Evidence for complex of sibling species of snapping shrimp (Caridea: Alpheidae: Alpheus) from the Caribbean and pacific coasts of Panama. Proc San Diego Soc Nat Hist. 1992, 18: 1-5.
Barber PH, S. R. Palumbi, M. V. Erdmann, Moosa MK: Sharp genetic breaks among populations of Haptosquilla pulchella (Stomatopoda) indicate limits to larval transport: patterns, causes, and consequences. 2002, 11: 659-674.
Jiggins CD, Mallet J: Bimodal hybrid zones and speciation. Trends Ecol Evol. 2000, 15: 250-255. 10.1016/S0169-5347(00)01873-5.
Knowlton N, Weil E, Weigt LA, Guzman HM: Sibling species in Montastraea annularis, coral bleaching, and the coral climate record. Science. 1992, 255: 330 333-10.1126/science.255.5042.330.
Weil E, Knowlton N: A multi-character analysis of the Caribbean coral Montastraea annularis (Ellis and Solander, 1786) and its two sibling species, M. faveolata (Ellis and Solander, 1786) and M. franksi (Gregory, 1895). Bull Mar Sci. 1994, 55: 151-175.
Palumbi SR, Metz E: Strong reproductive isolation in closely related tropical sea urchins (genus Echinometra). Mol Biol Evol. 1991, 8: 227-239.
Brazeau DA, Harvell CD: Genetic structure of local populations and divergence between growth forms in a clonal invertebrate, the Caribbean octocoral Briareum asbestinum. Mar Biol. 1994, 119: 53-60. 10.1007/BF00350106.
Carlon DB, Budd AF: Incipient speciation across a depth gradient in a scleractinian coral?. Evolution. 2002, 56: 2227-2242.
Duffy JE: Species boundaries, specialization, and the radiation of sponge-dwelling alpheid shrimp. Biol J Linn Soc. 1996, 58: 307-324.
Knowlton N, Jackson JBC: New taxonomy and niche partitioning on coral reefs: Jack of all trades or master of some?. Trends Ecol Evol. 1994, 9: 7-9. 10.1016/0169-5347(94)90224-0.
Santos SR, T.L. Shearer, A.R. Hannes, Coffroth. MA: Fine-scale diversity and specificity in the most prevalent lineage of symbiotic dinoflagellates (Symbiodinium, Dinophyceae) of the Caribbean. Mol Ecol. 2004, 13: 459-469. 10.1046/j.1365-294X.2003.02058.x.
Knowlton N, Mate JL, Guzman HM, Rowan R, Jara J: Direct evidence for reproductive isolation among the three species of the Montastraea annularis complex in Central America (Panama and Honduras). Mar Biol. 1997, 127: 705-711. 10.1007/s002270050061.
Lessios HA: Possible prezygotic reproductive isolation in sea urchins separated by the Isthmus of Panama. Evolution. 1984, 35: 1144-1148. 10.2307/2408446.
Levitan DR, Fukami H, Jara J, Kline D, McGovern TA, McGhee KM, Swanson CA, Knowlton N: Mechanisms of reproductive isolation among sympatric broadcast-spawning corals of the Montastraea annularis complex. Evolution. 2004, 58: 308-323.
Hellberg ME, Moy GW, Vacquier VD: Positive selection and propeptide repeats promote rapid Interspecific divergence of a gastropod sperm protein. Mol Biol Evol. 2000, 17 (3): 458-466.
Palumbi SR: Species formation and the evolution of gamete recognition loci. Endless Forms: Species and Speciation. Edited by: Howard DJ, Berlocher SH. 1998, NY , Oxford Univ. Press, 271-278.
Vacquier VD, Swanson WJ, Lee YH: Positive Darwinian selection on two homologous fertilization proteins: what is the selective pressure driving their divergence?. J Mol Evol. 1997, 44: 15-22. 10.1007/PL00000049.
Knowlton N: Sibling species in the sea. Annu Rev Ecol Syst. 1993, 24: 189-216. 10.1146/annurev.es.24.110193.001201.
Arnold ML: Natural Hybridization and Evolution. Oxford Series in Ecology and Evolution. 1997, New York , Oxford Univ. Press
Rieseberg LH, Fossen CV, Desrochers AM: Hybrid speciation accompanied by genomic reorganization in wild comparative mapping between Arabidopsis sunflowers. Nature. 1995, 375: 313-316. 10.1038/375313a0.
Veron JEN: Corals in space and time. Sydney. 1995, Sydney , Univ. of New South Wales Press, 321-
Kim K, Lasker HR: Flow-mediated competition among suspension feeding gorgonians. J Exp Mar Biol Ecol. 1997, 215: 49-64. 10.1016/S0022-0981(97)00015-4.
Skoufas G: Comparative biometry of Eunicella singularis (gorgonian) sclerites at East Mediterranean Sea (North Aegean Sea, Greece) . Mar Biol. 2006, 149: 1365-1370. 10.1007/s00227-006-0314-4.
Bayer FM: The Shallow-Water Octocorallia of the West Indian Region. A Manual for marine biologists. 1961, The Hague , Martinus Nijhoff, 373 pp-
West JM: Plasticity in the sclerites of a gorgonian coral: Tests of water motion, light level, and damage cues. Biol Bull. 1997, 192: 279-289. 10.2307/1542721.
Kinzie RA: The zonation of West Indian gorgonians. Bul Mar Sci. 1973, 23: 93-155.
Yoshioka PM, Yoshioka BB: Effects of water motion, topographic relief and sediment transport on the distribution of shallow-water gorgonian community. Mar Ecol Prog Ser. 1989, 54: 257-264. 10.3354/meps054257.
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. 10.1007/BF01731581.
Pigliucci M: Touchy and bushy: phenotypic plasticity and integration in response to wind stimulation in Arabidopsis thaliana. Int J Plant Sci. 2002, 163: 399-408. 10.1086/339158.
Okamura BJ, Partridge C: Suspension Feeding Adaptations to Extreme Flow Environments in a Marine Bryozoan. Biol Bull. 1999, 196: 205-215. 10.2307/1542566.
Callahan H, Pigliucci M: Shade-induced plasticity and its ecological significance in wild populations of Arabidopsis thaliana. Ecology. 2002, 83: 1965-1980.
Monro K, Poore AGB: Light quantity and quality induce shade-avoiding plasticity in a marine macroalga. J Evol Biol. 2005 , 18 : 426-435. 10.1111/j.1420-9101.2004.00826.x.
Graus RR, Macintyre. IG: Light control of growth forms in colonial reef corals: computer simulation. Science. 1976, 193: 895-897. 10.1126/science.193.4256.895.
Grigg RW: Orientation and growth form of sea fans. Limn Oceanogr. 1972, 17: 185-192.
Lesser MP, Weis VM, Patterson MP, Jokiel PL: Effects of morphology and water motion on carbon delivery and productivity in the reef coral, Pocillopora damicornis (Linnaeus): diffusion barriers, inorganic carbon limitation, and biochemical plasticity. J Exp Mar Biol Ecol. 1994, 178: 153-179. 10.1016/0022-0981(94)90034-5.
Sebens KP, Witting J, Helmuth B: Effects of water flow and branch spacing on particle capture by the reef coral Madracis mirabilis (Duchassaing and Michelotti). J Exp Mar Biol Ecol. 1997, 211: 1-28. 10.1016/S0022-0981(96)02636-6.
Meroz E, Brickner I, Loya Y, Peretzman-Shemer A, Ilan M: The effect of gravity on coral morphology. Proc R Soc Lond B. 2002, 269: 717-720. 10.1098/rspb.2001.1924.
Stambler N, Dubinsky Z: Corals as light collectors: an integrating sphere approach. Coral Reefs. 2005, 24: 1-9. 10.1007/s00338-004-0452-4.
de Kroon H, H.Huber, J.F.Stuefer, van Groenendael JM: A modular concept of phenotypic plasticity in plants. New Phytologist. 2005, 166: 73-82. 10.1111/j.1469-8137.2004.01310.x.
Yoshioka PM, Yoshioka BB: A comparison of the survivorship and growth of shallow-water gorgonian species of Puerto Rico. Mar Ecol Prog Ser. 1991, 69: 253-260. 10.3354/meps069253.
McFadden CS, France SC, Sánchez JA, Alderslade P: A molecular phylogenetic analysis of the Octocorallia (Cnidaria: Anthozoa) based on mitochondrial protein-coding sequences. Mol Phyl Evol. 2006, 41: 513-527. 10.1016/j.ympev.2006.06.010.
Sanchez JA, McFadden CS, France SC, Lasker HR: Molecular phylogenetic analyses of shallow-water Caribbean octocorals. Mar Biol. 2003, 142: 975-987.
Hellberg ME: Dependence of gene flow on geographic distance in two solitary corals with different larval dispersal capabilities. Evolution. 1996, 50: 1167-1175. 10.2307/2410657.
Avise JC: Molecular Markers, Natural History and Evolution. 2004, Sunderland MA , Sinauer, 684 p-2nd
Mayr E: Systematics and the Origin of Species. 1942, New York , Columbia Univ. Press
Baums IB, Miller MW, Hellberg ME: Regionally isolated populations of an imperiled Caribbean coral, Acropora palmata. Molecular Ecology. 2005, 14: 1377-1390. 10.1111/j.1365-294X.2005.02489.x.
Gutiérrez-Rodríguez C, Lasker HR: Microsatellite variation reveals high levels of genetic variability and population structure in the gorgonian coral Pseudopterogorgia elisabethae across the Bahamas. Mol Ecol. 2004, 13: 2211-2221. 10.1111/j.1365-294X.2004.02247.x.
Taylor MS, Hellberg ME: Marine radiations at small geographic scales: speciation in Neotropical reef gobies (Elacatinus). Evolution. 2005, 59: 374-385.
Vollmer SV, Palumbi SR: Restricted gene flow in the Caribbean staghorn coral Acropora cervicornis: Implications for the recovery of endangered reefs. J Heredity. 2007, 98: 40-50. 10.1093/jhered/esl057.
Briggs JC: Operation of zoogeographic barriers. Syst Zool. 1973, 23: 248-256. 10.2307/2412136.
Shulman MJ, Bermingham E: Early life histories, ocean currents, and the population genetics of Caribbean reef fishes. Evolution. 1995, 49: 897-910. 10.2307/2410412.
Zardus JD, Hadfield MG: Multiple origins and incursions of the Atlantic barnacle Chthamalus proteus in the Pacific. Mol Ecol. 2005, 14 (12): 3719–3733-10.1111/j.1365-294X.2005.02701.x.
Palumbi SR: Marine speciation on a small planet. Trends Ecol Evol. 1992, 7: 114-117. 10.1016/0169-5347(92)90144-Z.
Diekmann OE, Olsen JL, Stam WT, Bak RPM: Genetic variation within Symbiodinium clade B from the coral genus Madracis in the Caribbean (Netherlands Antilles). Coral Reefs. 2003, 22: 29-33.
Rowan R, Knowlton N: Intraspecific diversity and ecological zonation in coral-algal symbiosis. Proc Natl Acad Sci. 1995, 92: 2850-2853. 10.1073/pnas.92.7.2850.
Rieseberg LH, Carter R, Zona S: Molecular tests of the hypothesized hybrid origin of two diploid Helianthus species (Asteraceae). Evolution. 1990, 44: 1498-1511. 10.2307/2409332.
Pyle RL, Randall JE: A review of hybridization in marine angelfishes (Perciformes: Pomacanthidae) . Environ Biol Fishes. 1994, 41: 127-145.
Schelly R, Salzburger W, Koblmuller S, Duftner N, Sturmbauer C: Phylogenetic relationships of the lamprologine cichlid genus Lepidiolamprologus (Teleostei: Perciformes) based on mitochondrial and nuclear sequences, suggesting introgressive hybridization. Mol Phyl Evol. 2006, 38: 426-438. 10.1016/j.ympev.2005.04.023.
Gilg MR, Hilbish TJ: Patterns of larval dispersal and their effect on the maintenance of a blue mussel hybrid zone in Southwest England. Evolution. 2003, 57: 1061-1077.
Hatta M, Fukami H, Wang WQ, Omori M, Shimoike K, Hayashibara T, Ina Y, Sugiyama T: Reproductive and genetic evidence for a reticulate evolutionary history of mass-spawning corals. Mol Biol Evol. 1999, 16: 1607-1613.
van Oppen MJH, Willis BL, van Vugt HWJA, Miller DJ: Examination of species boundaries in the Acropora cervicornis group (Scleractinia, Cnidaria) using nuclear DNA sequence analyses. Mol Ecol. 2000, 9: 1363-1373. 10.1046/j.1365-294x.2000.01010.x.
Vollmer SV, Palumbi SR: Hybridization and the evolution of reef coral diversity. Science. 2002, 296: 2023-2025. 10.1126/science.1069524.
Kenyon JC: Models of reticulate evolution in the coral genus Acropora based on chromosome numbers: parallels with plants. Evolution. 1997, 51: 756-767. 10.2307/2411152.
Willis BL, van Oppen MJH, Miller DJ, V. VS, Ayre DJ: The role of hybridization in the evolution of reef corals. Annu Rev Ecol Syst. 2006, 37: 489-517. 10.1146/annurev.ecolsys.37.091305.110136.
Burke JM, Rieseberg LH: The fitness effects of transgenic disease resistance in wild sunflowers. Science. 2003, 300: 1250-10.1126/science.1084960.
Gompert Z, Fordyce JA, Forister ML, Shapiro AM, Nice CC: Homoploid hybrid speciation in an extreme habitat. Science. 2006, 314: 1923-1925. 10.1126/science.1135875.
Sánchez JA, Aguilar C, Dorado D, Manrique N: Phenotypic plasticity and morphological integration in a marine modular invertebrate. BMC Evolutionary Biology. 2007, 7: 122-10.1186/1471-2148-7-122.
Lasker HR, Coffroth MA: Octocoral distribution at Carrie Bow Cay, Belize. Mar Ecol Prog Ser. 1983, 13: 21-28. 10.3354/meps013021.
Beiring EA, Lasker HR: Egg production by colonies of a gorgonian coral. Mar Ecol Prog Ser. 2000, 196: 169-177. 10.3354/meps196169.
Lasker HR: Asexual reproduction, fragmentation, and skeletal morphology of a plexaurid gorgonian. Mar Ecol Prog Ser. 1984, 19: 261-268. 10.3354/meps019261.
Garcia JR, Schmitt C, Heberer G, Winter A: La Parguera Puerto Rico UNESCO. CARICOMP- Caribbean coral reef, sea grass and mangrove sites. 1998, Paris , UNESCO, 347-
Prada C: Phenotypic and genetic variability in the octocoral Plexaura flexuosa. Master's thesis. Department of Marine Sciences. 2007, Puerto Rico , University of Puerto Rico, Mayaguez. , M.S.:
Sokal RR, Rohlf FJ: Biometry. 3rd ed. 1995, New York , W. H. and Freeman
Quinn GP, Keough MJ: Experimental design and data analysis for biologists. 2002, Cambridge, UK , Cambridge Univ. Press
Pont-Kingdon GA, Okada NA, Macfarlane JL, Beagley CT, Wolstenholme DR, Cavalier-Smith T, Clark-Walker GD: A coral mitochondrial mutS gene. Nature. 1995, 375: 109-111. 10.1038/375109b0.
Takabayashi M, Carter DA, Loh WKW, Hoegh-Guldberg O: A coral specific primer for PCR amplification of the internal transcribed spacer region in ribosomal DNA. Mol Ecol. 1998, 7: 925-931. 10.1111/j.1365-294X.1998.00368.x.
France SC, Hoover LL: Analysis of variation in mitochondrial DNA sequences (ND3, ND4L, MSH) among Octocorallia (=Alcyonaria) (Cnidaria: Anthozoa). Bull Biol Soc Wash. 2001, 10: 110-118.
Gordon D: Viewing and editing assembled sequences using consed. Current Protocols in Bioinformatics. Edited by: Baxevanis AD DDB. 2004, New York , John Wiley & Co, 11.12.11-11.12.43.
Maddison WP, Maddison DR: MacClade version 3: Analysis of phylogeny and character evolution. (book) + 900K (computer program). 1992, Sunderland MA , Sinauer, 398 pp-
Nei M: Molecular evolutionary genetics. 1987, New York , Columbia Univ. Press
Rozas J, Sandez-DelBurrio JC, Messeguer X, Rozas. R: DnaSP: DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359.
Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 105: 437-460.
Fu YX: Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997, 147 (2): 915-925.
Schneider S, Roessli D, Excoffier L: ARLEQUIN. A software for population genetic data analysis, Version 2.00. 2000, Geneva, Switzerland , Genetics and Biometry Laboratory, Department of Anthropology, Univ. of Geneva
Templeton AR, Crandall KA, Sing CF: A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation. Genetics. 1992, 132: 619-633.
Clement M, Posada D, Crandall K: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9(10): 1657-1660. 10.1046/j.1365-294x.2000.01020.x.
Excoffier L, Smouse PE, Quattro JM: Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial restriction data. Genetics. 1992, 131: 479-491.
Swofford DL: PAUP*: Phylogenetic analysis using parsimony*, version 4.0b10. 2002, Sunderland, MA , Sinauer
Posada D, Crandall KA: Modeltest: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.
Hasegawa M, Kishino H, Yano T: Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol. 1985, 22: 160-174. 10.1007/BF02101694.
McGehee A: Comparisons of water motion in coral reef by measuring corrosion rates of dissimilar metals. Carib J Sci. 1998, 34: 286-297.
We thank Ernesto Weil, Aldo Croquer, Joe Staton, Raul Machiavelli, Hernan Torres-Pratts and two anonymous reviewers for their comments. Alex Mercado, Diana Beltrán, and Elizabeth McLean helped in the field. Prada thanks Cielo Montoya for her unconditional support. This work was presented as partial fulfilment of the requirements for Prada's Masters Thesis at the Department of Marine Sciences of UPRM. Funds for this project came from Sigma Xi research-in-aid grant (G200610111157185072) and Sea Grant Seed Money awarded to Prada and Schizas, a "Restoration of Gorgonian populations" Grant (535480) from the Sea Grant program of the UPRM, and the Coral Reef Ecosystem Studies from NOAA to Yoshioka and a Caribbean Coral Reef Institute grant awarded to Schizas. SEMs were made possible through the help of José Almodovar. DNA sequencing was performed in the Sequencing and Genotyping facility, University of Puerto Rico-Río Piedras which is supported in part by NCRR AABRE Grant (P20 RR16470), NIH-SCORE Grant (S06GM08102), University of Puerto Rico Biology Department, NSFCREST Grant (0206200).
CP conceived the project, collected the data, carried out the analyses and wrote the manuscript. PY and NVS advised on analyses, NVS provided samples and lab materials. PY and NVS participated in discussions. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Prada, C., Schizas, N.V. & Yoshioka, P.M. Phenotypic plasticity or speciation? A case from a clonal marine organism. BMC Evol Biol 8, 47 (2008). https://doi.org/10.1186/1471-2148-8-47
- Phenotypic Plasticity
- Discriminant Function Analysis
- Outer Reef
- Shallow Habitat
- Deep Habitat