skip to main content

Title: Speciation‐by‐depth on coral reefs: Sympatric divergence with gene flow or cryptic transient isolation?

The distributions of many sister species in the sea overlap geographically but are partitioned along depth gradients. The genetic changes leading to depth segregation may evolve in geographic isolation as a prerequisite to coexistence or may emerge during primary divergence leading to new species. These alternatives can now be distinguished via the power endowed by the thousands of scorable loci provided by second‐generation sequence data. Here, we revisit the case of two depth‐segregated, genetically isolated ecotypes of the nominal Caribbean candelabrum coralEunicea flexuosa. Previous analyses based on a handful of markers could not distinguish between models of genetic exchange after a period of isolation (consistent with secondary contact) and divergence with gene flow (consistent with primary divergence). Analyses of the history of isolation, genetic exchange and population size based on 15,640 new SNP markers derived from RNAseq data best support models where divergence began 800K BP and include epochs of divergence with gene flow, but with an intermediate period of transient isolation. Results also supported the previous conclusion that recent exchange between the ecotypes occurs asymmetrically from the Shallow lineage to the Deep. Parallel analyses of data from two other corals with depth‐segregated populations (Agaricia fragilisandPocillopora damicornis) suggest divergence leading to depth‐segregated populations may begin with a period of symmetric exchange, but that an epoch of population isolation precedes more complete isolation marked by asymmetric introgression. Thus, while divergence‐with‐gene flow may account for much of the differentiation that separates closely related, depth‐segregated species, it remains to be seen whether any critical steps in the speciation process only occur when populations are isolated.

more » « less
Author(s) / Creator(s):
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Journal of Evolutionary Biology
Medium: X Size: p. 128-137
["p. 128-137"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Context

    Processes that shape genomic and ecological divergence can reveal important evolutionary dynamics to inform the conservation of threatened species.Fontaineais a genus of rainforest shrubs and small trees including critically endangered and threatened species restricted to narrow, but complex geographic and ecological regions. Several species ofFontaineaare subject to spatially explicit conditions and experience limited intra-specific gene flow, likely generating genetic differentiation and local adaptation.


    Here, we explored the genetic and ecological mechanisms underlying patterns of diversification in two, closely related threatenedFontaineaspecies. Our aim was to compare spatial patterns of genetic variation between the vulnerableFontainea australis(Southern Fontainea) and critically endangeredF. oraria(Coastal Fontainea), endemic to the heterogeneous subtropical region of central, eastern Australia, where large-scale clearing has severely reduced rainforest habitat to a fraction (< 1%) of its pre-European settlement extent.


    We used a set of 10,000 reduced-representation markers to infer genetic relationships and the drivers of spatial genetic variation across the two species. In addition, we employed a combination of univariate and multivariate genome-environment association analysis using a set of topo-climatic variables to explore potential patterns of local adaptation as a factor impacting genomic divergence.


    Our study revealed that Coastal Fontainea have a close genetic relationship with Southern Fontainea. We showed that isolation by distance has played a key role in their genetic variation, indicating that vicariance can explain the spatial genetic distribution of the two species. Genotype-environment analyses showed a strong association with temperature and topographic features, suggesting adaptation to localised thermal environments. We used a multivariate redundancy analysis to identify a range of putatively adapted loci associated with local environmental conditions.


    Divergent selection at the local-habitat scale as a result of dispersal limitations and environmental heterogeneity (including physical barriers) are likely contributors to adaptive divergence between the twoFontaineaspecies. Our findings have presented evidence to indicate that Southern and Coastal Fontainea were comprised of distinct genetic groups and ecotypes, that together may form a single species continuum, with further phenotype research suggested to confirm the current species boundaries. Proactive conservation actions, including assisted migration to enhance the resilience of populations lacking stress-tolerant single nucleotide polymorphisms (SNPs) may be required to secure the long-term future of both taxa. This is especially vital for the critically endangered Coastal Fontainea given projections of habitat decline for the species under future climate scenarios.

    more » « less
  2. Background

    Vestimentiferan tubeworms are some of the most recognizable fauna found at deep-sea cold seeps, isolated environments where hydrocarbon rich fluids fuel biological communities. Several studies have investigated tubeworm population structure; however, much is still unknown about larval dispersal patterns at Gulf of Mexico (GoM) seeps. As such, researchers have applied microsatellite markers as a measure for documenting the transport of vestimentiferan individuals. In the present study, we investigate the utility of microsatellites to be cross-amplified within the escarpiid clade of seep vestimentiferans, by determining if loci originally developed forEscarpiaspp. could be amplified in the GoM seep tubeworm,Seepiophila jonesi. Additionally, we determine if cross-amplified loci can reliably uncover the same signatures of high gene flow seen in a previous investigation ofS. jonesi.


    Seventy-sevenS. jonesiindividuals were collected from eight seep sites across the upper Louisiana slope (<1,000 m) in the GoM. Forty-eight microsatellite loci that were originally developed forEscarpia laminata(18 loci) andEscarpia southwardae(30 loci) were tested to determine if they were homologous and polymorphic inS. jonesi. Loci found to be both polymorphic and of high quality were used to test for significant population structuring inS. jonesi.


    Microsatellite pre-screening identified 13 (27%) of theEscarpialoci were homologous and polymorphic inS. jonesi, revealing that microsatellites can be amplified within the escarpiid clade of vestimentiferans. Our findings uncovered low levels of heterozygosity and a lack of genetic differentiation amongstS. jonesifrom various sites and regions, in line with previous investigations that employed species-specific polymorphic loci onS. jonesiindividuals retrieved from both the same and different seep sites. The lack of genetic structure identified from these populations supports the presence of significant gene flow via larval dispersal in mixed oceanic currents.


    The ability to develop “universal” microsatellites reduces the costs associated with these analyses and allows researchers to track and investigate a wider array of taxa, which is particularly useful for organisms living at inaccessible locations such as the deep sea. Our study highlights that non-species specific microsatellites can be amplified across large evolutionary distances and still yield similar findings as species-specific loci. Further, these results show thatS. jonesicollected from various localities in the GoM represents a single panmictic population, suggesting that dispersal of lecithotrophic larvae by deep sea currents is sufficient to homogenize populations. These data are consistent with the high levels of gene flow seen inEscarpiaspp., which advocates that differences in microhabitats of seep localities lead to variation in biogeography of separate species.

    more » « less
  3. Abstract

    Although the impact of Pleistocene glacial cycles on the diversification of the tropical biota was once dismissed, increasing evidence suggests that Pleistocene climatic fluctuations greatly affected the distribution and population divergence of tropical organisms. Landscape genomic analyses coupled with paleoclimatic distribution models provide a powerful way to understand the consequences of past climate changes on the present‐day tropical biota. Using genome‐wide SNP data and mitochondrial DNA, combined with projections of the species distribution across the late Quaternary until the present, we evaluate the effect of paleoclimatic shifts on the genetic structure and population differentiation ofHypsiboas lundii, a treefrog endemic to the South American Cerrado savanna. Our results show a recent and strong genetic divergence inH. lundiiacross the Cerrado landscape, yielding four genetic clusters that do not seem congruent with any current physical barrier to gene flow. Isolation by distance (IBD) explains some of the population differentiation, but we also find strong support for past climate changes promoting range shifts and structuring populations even in the presence of IBD. Post‐Pleistocene population persistence in four main areas of historical stable climate in the Cerrado seems to have played a major role establishing the present genetic structure of this treefrog. This pattern is consistent with a model of reduced gene flow in areas with high climatic instability promoting isolation of populations, defined here as “isolation by instability,” highlighting the effects of Pleistocene climatic fluctuations structuring populations in tropical savannas.

    more » « less
  4. Abstract

    Divergence is often ephemeral, and populations that diverge in response to regional topographic and climatic factors may not remain reproductively isolated when they come into secondary contact. We investigated the geographical structure and evolutionary history of population divergence withinSceloporus occidentalis(western fence lizard), a habitat generalist with a broad distribution that spans the major biogeographical regions of Western North America. We used double digest RAD sequencing to infer population structure, phylogeny and demography. Population genetic structure is hierarchical and geographically structured with evidence for gene flow between biogeographical regions. Consistent with the isolation–expansion model of divergence during Quaternary glacial–interglacial cycles, gene flow and secondary contact are supported as important processes explaining the demographic histories of populations. Although populations may have diverged as they spread northward in a ring‐like manner around the Sierra Nevada and southern Cascade Ranges, there is strong evidence for gene flow among populations at the northern terminus of the ring. We propose the concept of an “ephemeral ring species” and contrastS. occidentaliswith the classic North American ring species,Ensatina eschscholtzii. Contrary to expectations of lower genetic diversity at northern latitudes following post‐Quaternary‐glaciation expansion, the ephemeral nature of divergence inS. occidentalishas produced centres of high genetic diversity for different reasons in the south (long‐term stability) vs. the north (secondary contact).

    more » « less
  5. Abstract

    Effective management of threatened and exploited species requires an understanding of both the genetic connectivity among populations and local adaptation. The Olympia oyster (Ostrea lurida), patchily distributed from Baja California to the central coast of Canada, has a long history of population declines due to anthropogenic stressors. For such coastal marine species, population structure could follow a continuous isolation‐by‐distance model, contain regional blocks of genetic similarity separated by barriers to gene flow, or be consistent with a null model of no population structure. To distinguish between these hypotheses inO. lurida, 13,424 single nucleotide polymorphisms (SNPs) were used to characterize rangewide population structure, genetic connectivity, and adaptive divergence. Samples were collected across the species range on the west coast of North America, from southern California to Vancouver Island. A conservative approach for detecting putative loci under selection identified 235SNPs across 129GBSloci, which were functionally annotated and analyzed separately from the remaining neutral loci. While strong population structure was observed on a regional scale in both neutral and outlier markers, neutral markers had greater power to detect fine‐scale structure. Geographic regions of reduced gene flow aligned with known marine biogeographic barriers, such as Cape Mendocino, Monterey Bay, and the currents around Cape Flattery. The outlier loci identified as under putative selection included genes involved in developmental regulation, sensory information processing, energy metabolism, immune response, and muscle contraction. These loci are excellent candidates for future research and may provide targets for genetic monitoring programs. Beyond specific applications for restoration and management of the Olympia oyster, this study lends to the growing body of evidence for both population structure and adaptive differentiation across a range of marine species exhibiting the potential for panmixia. Computational notebooks are available to facilitate reproducibility and future open‐sourced research on the population structure ofO. lurida.

    more » « less