Polygenic evolution drives species divergence and climate adaptation in corals

In-depth single-cell transcriptomic exploration of the regenerative dynamics in stony coral

Coral reef ecosystems face escalating threats from anthropogenic global climate challenges, leading to frequent bleaching events. A key issue in coral transplantation is the inability of fragments to rapidly grow to sizes that can resist environmental pressures. The observation of accelerated growth during the early stages of coral regeneration provides new insights for addressing this challenge. To investigate the underlying molecular mechanisms, we study the fast-growing stony coral Acropora muricata. Using single-cell RNA sequencing, bulk RNA sequencing, and high-resolution micro-computed tomography, we identify a critical regeneration phase around 2-4 weeks post-injury. Single-cell transcriptome analysis reveals 11 function-specific cell clusters. Pseudotime analysis indicates epidermal cell differentiation into calicoblasts. Bulk RNA-seq results highlight a temporal limitation in coral's rapid regeneration. Through integrated multi-omics analysis, this study emphasizes the importance of a comprehensive understanding of coral regeneration, providing insights beyond fundamental knowledge and offering potential protective strategies to promote coral growth.

Genetic Basis of Reproductive Isolation in Torrey Pine (Pinus torreyana Parry): Insights From Hybridization and Adaptation

Tree species are often locally adapted to their environments, but the extent to which environmental adaptation contributes to incipient speciation is unclear. One of the rarest pines in the world, Torrey pine (Pinus torreyana Parry), persists naturally across one island and one mainland population in southern California. The two populations are morphologically and genetically differentiated but experience some connectivity, making it an ideal system for assessing the evolution of reproductive isolation. Previous work has found evidence of heterosis in F1 mainland-island hybrids, suggesting genetic rescue could be beneficial in the absence of reproductive barriers. Using ddRADseq and GWAS for a common garden experiment of island, mainland, and F1 individuals, we identified candidate loci for environmentally driven reproductive isolation, their function, and their relationship to fitness proxies. By simulating neutral evolution and admixture between the two populations, we identified loci that exhibited reduced heterozygosity in the F1s, evidence of selection against admixture. SNPs with reduced F1 heterozygosity were enriched for growth and pollination functions, suggesting genetic variants that could be involved in the evolution of reproductive barriers between populations. One locus with reduced F1 heterozygosity exhibited strong associations with growth and reproductive fitness proxies in the common garden, with the mainland allele conferring increased fitness. If this locus experiences divergent selection in the two natural populations, it could promote their reproductive isolation. Finally, although hybridization largely reduced allele fixation in the F1s initially, indicating heterosis is likely due to the masking of deleterious alleles, the emergence of reproductive isolation between populations may diminish the longer-term benefits of genetic rescue in F2 or advanced-generation hybrids. As Torrey pine is a candidate for interpopulation genetic rescue, caution is warranted where longer-term gene flow between diverged populations may result in reduced fitness if barriers have evolved.

A Genetic Variant of Delta-9 Desaturase Is Associated With Latitudinal Adaptation in a Coral from the Great Barrier Reef

Coral populations across the Great Barrier Reef (GBR) could rapidly adapt to the warming climate if they have standing genetic variation for thermal tolerance. Here, we describe a locus likely involved in latitudinal adaptation of Acropora millepora. This locus shows a steep latitudinal gradient of derived allele frequency increasing at higher latitudes, and harbours a cluster of eight tandemly repeated Δ9-desaturase genes adjacent to a region in the genome where a hard selective sweep likely occurred. In colonies reciprocally transplanted across 4.5° of latitude, the expression of Δ9-desaturase is upregulated at the high-latitude reef. Furthermore, corals from the low-latitude reef bearing the derived Δ9-desaturase allele express the gene more and grow faster than their peers when transplanted to the high-latitude reef. In other organisms ranging from bacteria to fish, Δ9-desaturase is upregulated under cold conditions to adjust membrane fluidity by introducing double bonds into fatty acid chains of membrane lipids. It is therefore plausible that the signal of latitudinal adaptation at the Δ9-desaturase locus is due to its involvement in adaptation to cooler temperatures at higher latitudes.

Cryptic coral diversity is associated with symbioses, physiology, and response to thermal challenge

Coral persistence in the Anthropocene depends on interactions among holobiont partners (coral animals and microbial symbionts) and their environment. Cryptic coral lineages-genetically distinct yet morphologically similar groups-are critically important as they often exhibit functional diversity relevant to thermal tolerance. In addition, environmental parameters such as thermal variability may promote tolerance, but how variability interacts with holobiont partners to shape responses to thermal challenge remains unclear. Here, we identified three cryptic lineages of Siderastrea siderea in Bocas del Toro, Panamá that differ in distributions across inshore and offshore reefs, microbial associations, phenotypic traits of holobiont partners (i.e., phenomes), and skeleton morphologies. A thermal variability experiment failed to increase thermal tolerance, but subsequent thermal challenge and recovery revealed that one lineage maintained elevated energetic reserves, photochemical efficiency, and growth. Last, coral cores highlighted that this lineage also exhibited greater growth historically. Functional variation among cryptic lineages highlights their importance in predicting coral reef responses to climate change.

Population Genomics Reveals Local Adaptation Related to Temperature Variation in Two Stream Frog Species: Implications for Vulnerability to Climate Warming

Identifying populations at highest risk from climate change is a critical component of conservation efforts. However, vulnerability assessments are usually applied at the species level, even though intraspecific variation in exposure, sensitivity and adaptive capacity play a crucial role in determining vulnerability. Genomic data can inform intraspecific vulnerability by identifying signatures of local adaptation that reflect population-level variation in sensitivity and adaptive capacity. Here, we address the question of local adaptation to temperature and the genetic basis of thermal tolerance in two stream frogs (Ascaphus truei and A. montanus). Building on previous physiological and temperature data, we used whole-genome resequencing of tadpoles from four sites spanning temperature gradients in each species to test for signatures of local adaptation. To support these analyses, we developed the first annotated reference genome for A. truei. We then expanded the geographic scope of our analysis using targeted capture at an additional 11 sites per species. We found evidence of local adaptation to temperature based on physiological and genomic data in A. montanus and genomic data in A. truei, suggesting similar levels of sensitivity (i.e., susceptibility) among populations regardless of stream temperature. However, invariant thermal tolerances across temperatures in A. truei suggest that populations occupying warmer streams may be most sensitive. We identified high levels of evolutionary potential in both species based on genomic and physiological data. While further integration of these data is needed to comprehensively evaluate spatial variation in vulnerability, this work illustrates the value of genomics in identifying spatial patterns of climate change vulnerability.

Evolutionary adaptation under climate change: Aedes sp. demonstrates potential to adapt to warming

Climate warming is expected to shift the distributions of mosquitoes and mosquito-borne diseases, promoting expansions at cool range edges and contractions at warm range edges. However, whether mosquito populations could maintain their warm edges through evolutionary adaptation remains unknown. Here, we investigate the potential for thermal adaptation in Aedes sierrensis, a congener of the major disease vector species that experiences large thermal gradients in its native range, by assaying tolerance to prolonged and acute heat exposure, and its genetic basis in a diverse, field-derived population. We found pervasive evidence of heritable genetic variation in mosquito heat tolerance, and phenotypic trade-offs in tolerance to prolonged versus acute heat exposure. Further, we found genomic variation associated with prolonged heat tolerance was clustered in several regions of the genome, suggesting the presence of larger structural variants such as chromosomal inversions. A simple evolutionary model based on our data estimates that the maximum rate of evolutionary adaptation in mosquito heat tolerance will exceed the projected rate of climate warming, implying the potential for mosquitoes to track warming via genetic adaptation.

Climate adaptive loci revealed by seascape genomics correlate with phenotypic variation in heat tolerance of the coral Acropora millepora

One of the main challenges in coral reef conservation and restoration is the identification of coral populations resilient under global warming. Seascape genomics is a powerful tool to uncover genetic markers potentially involved in heat tolerance among large populations without prior information on phenotypes. Here, we aimed to provide first insights on the role of candidate heat associated loci identified using seascape genomics in driving the phenotypic response of Acropora millepora from New Caledonia to thermal stress. We subjected 7 colonies to a long-term ex-situ heat stress assay (4 °C above the maximum monthly mean) and investigated their physiological response along with their Symbiodiniaceae communities and genotypes. Despite sharing similar thermal histories and associated symbionts, these conspecific individuals differed greatly in their tolerance to heat stress. More importantly, the clustering of individuals based on their genotype at heat-associated loci matched the phenotypic variation in heat tolerance. Colonies that sustained on average lower mortality, higher Symbiodiniaceae/chlorophyll concentrations and photosynthetic efficiency under prolonged heat stress were also the closest based on their genotypes, although the low sample size prevented testing loci predictive accuracy. Together these preliminary results support the relevance of coupling seascape genomics and long-term heat stress experiments in the future, to evaluate the effect size of candidate heat associated loci and pave the way for genomic predictive models of corals heat tolerance.

Evolutionary adaptation under climate change: Aedes sp. demonstrates potential to adapt to warming

Climate warming is expected to shift the distributions of mosquitoes and mosquito-borne diseases, facilitating expansions at cool range edges and contractions at warm range edges. However, whether mosquito populations could maintain their warm edges through evolutionary adaptation remains unknown. Here, we investigate the potential for thermal adaptation in Aedes sierrensis, a congener of the major disease vector species that experiences large thermal gradients in its native range, by assaying tolerance to prolonged and acute heat exposure, and its genetic basis in a diverse, field-derived population. We found pervasive evidence of heritable genetic variation in acute heat tolerance, which phenotypically trades off with tolerance to prolonged heat exposure. A simple evolutionary model based on our data shows that the estimated maximum rate of evolutionary adaptation in mosquito heat tolerance typically exceeds that of projected climate warming under idealized conditions. Our findings indicate that natural mosquito populations may have the potential to track projected warming via genetic adaptation. Prior climate-based projections may thus underestimate the range of mosquito and mosquito-borne disease distributions under future climate conditions.

Complex parental effects impact variation in larval thermal tolerance in a vertically transmitting coral

Coral populations must be able to adapt to changing environmental conditions for coral reefs to persist under climate change. The adaptive potential of these organisms is difficult to forecast due to complex interactions between the host animal, dinoflagellate symbionts and the environment. Here we created 26 larval families from six Montipora capitata colonies from a single reef, showing significant, heritable variation in thermal tolerance. Our results indicate that 9.1% of larvae are expected to exhibit four times the thermal tolerance of the general population. Differences in larval thermotolerance were driven mainly by maternal contributions, but we found no evidence that these effects were driven by symbiont identity despite vertical transmission from the dam. We also document no evidence of reproductive incompatibility attributable to symbiont identity. These data demonstrate significant genetic variation within this population which provides the raw material upon which natural selection can act.

Incorporating adaptive genomic variation into predictive models for invasion risk assessment

Global climate change is expected to accelerate biological invasions, necessitating accurate risk forecasting and management strategies. However, current invasion risk assessments often overlook adaptive genomic variation, which plays a significant role in the persistence and expansion of invasive populations. Here we used Molgula manhattensis, a highly invasive ascidian, as a model to assess its invasion risks along Chinese coasts under climate change. Through population genomics analyses, we identified two genetic clusters, the north and south clusters, based on geographic distributions. To predict invasion risks, we employed the gradient forest and species distribution models to calculate genomic offset and species habitat suitability, respectively. These approaches yielded distinct predictions: the gradient forest model suggested a greater genomic offset to future climatic conditions for the north cluster (i.e., lower invasion risks), while the species distribution model indicated higher future habitat suitability for the same cluster (i.e, higher invasion risks). By integrating these models, we found that the south cluster exhibited minor genome-niche disruptions in the future, indicating higher invasion risks. Our study highlights the complementary roles of genomic offset and habitat suitability in assessing invasion risks under climate change. Moreover, incorporating adaptive genomic variation into predictive models can significantly enhance future invasion risk predictions and enable effective management strategies for biological invasions in the future.

Finding genes and pathways that underlie coral adaptation

Mass coral bleaching is one of the clearest threats of climate change to the persistence of marine biodiversity. Despite the negative impacts of bleaching on coral health and survival, some corals may be able to rapidly adapt to warming ocean temperatures. Thus, a significant focus in coral research is identifying the genes and pathways underlying coral heat adaptation. Here, we review state-of-the-art methods that may enable the discovery of heat-adaptive loci in corals and identify four main knowledge gaps. To fill these gaps, we describe an experimental approach combining seascape genomics with CRISPR/Cas9 gene editing to discover and validate heat-adaptive loci. Finally, we discuss how information on adaptive genotypes could be used in coral reef conservation and management strategies.

(Epi)genomic adaptation driven by fine geographical scale environmental heterogeneity after recent biological invasions

Elucidating processes and mechanisms involved in rapid local adaptation to varied environments is a poorly understood but crucial component in management of invasive species. Recent studies have proposed that genetic and epigenetic variation could both contribute to ecological adaptation, yet it remains unclear on the interplay between these two components underpinning rapid adaptation in wild animal populations. To assess their respective contributions to local adaptation, we explored epigenomic and genomic responses to environmental heterogeneity in eight recently colonized ascidian (Ciona intestinalis) populations at a relatively fine geographical scale. Based on MethylRADseq data, we detected strong patterns of local environment-driven DNA methylation divergence among populations, significant epigenetic isolation by environment (IBE), and a large number of local environment-associated epigenetic loci. Meanwhile, multiple genetic analyses based on single nucleotide polymorphisms (SNPs) showed genomic footprints of divergent selection. In addition, for five genetically similar populations, we detected significant methylation divergence and local environment-driven methylation patterns, indicating the strong effects of local environments on epigenetic variation. From a functional perspective, a majority of functional genes, Gene Ontology (GO) terms, and biological pathways were largely specific to one of these two types of variation, suggesting partial independence between epigenetic and genetic adaptation. The methylation quantitative trait loci (mQTL) analysis showed that the genetic variation explained only 18.67% of methylation variation, further confirming the autonomous relationship between these two types of variation. Altogether, we highlight the complementary interplay of genetic and epigenetic variation involved in local adaptation, which may jointly promote populations' rapid adaptive capacity and successful invasions in different environments. The findings here provide valuable insights into interactions between invaders and local environments to allow invasive species to rapidly spread, thus contributing to better prediction of invasion success and development of management strategies.

Exploring coral speciation: Multiple sympatric Stylophora pistillata taxa along a divergence continuum on the Great Barrier Reef

Understanding how biodiversity originates and is maintained are fundamental challenge in evolutionary biology. Speciation is a continuous process and progression along this continuum depends on the interplay between evolutionary forces driving divergence and forces promoting genetic homogenisation. Coral reefs are broadly connected yet highly heterogeneous ecosystems, and divergence with gene flow at small spatial scales might therefore be common. Genomic studies are increasingly revealing the existence of closely related and sympatric taxa within taxonomic coral species, but the extent to which these taxa might still be exchanging genes and sharing environmental niches is unclear. In this study, we sampled extensively across diverse habitats at multiple reefs of the Great Barrier Reef (GBR) and comprehensively examined genome-wide diversity and divergence histories within and among taxa of the Stylophora pistillata species complex. S. pistillata is one of the most abundant and well-studied coral species, yet we discovered five distinct taxa, with wide geographic ranges and extensive sympatry. Demographic modelling showed that speciation events have occurred with gene flow and that taxa are at different stages along a divergence continuum. We found significant correlations between genetic divergence and specific environmental variables, suggesting that niche partitioning may have played a role in speciation and that S. pistillata taxa might be differentially adapted to different environments. Conservation actions rely on estimates of species richness, population sizes and species ranges, which are biased if divergent taxa are lumped together. As coral reefs are rapidly degrading due to climate change, our study highlights the importance of recognising evolutionarily distinct and differentially adapted coral taxa to improve conservation and restoration efforts aiming at protecting coral genetic diversity.

Coral Reef Population Genomics in an Age of Global Change

Coral reefs are both exceptionally biodiverse and threatened by climate change and other human activities. Here, we review population genomic processes in coral reef taxa and their importance for understanding responses to global change. Many taxa on coral reefs are characterized by weak genetic drift, extensive gene flow, and strong selection from complex biotic and abiotic environments, which together present a fascinating test of microevolutionary theory. Selection, gene flow, and hybridization have played and will continue to play an important role in the adaptation or extinction of coral reef taxa in the face of rapid environmental change, but research remains exceptionally limited compared to the urgent needs. Critical areas for future investigation include understanding evolutionary potential and the mechanisms of local adaptation, developing historical baselines, and building greater research capacity in the countries where most reef diversity is concentrated.

Marine heatwaves threaten cryptic coral diversity and erode associations among coevolving partners

Climate change-amplified marine heatwaves can drive extensive mortality in foundation species. However, a paucity of longitudinal genomic datasets has impeded understanding of how these rapid selection events alter cryptic genetic structure. Heatwave impacts may be exacerbated in species that engage in obligate symbioses, where the genetics of multiple coevolving taxa may be affected. Here, we tracked the symbiotic associations of reef-building corals for 6 years through a prolonged heatwave, including known survivorship for 79 of 315 colonies. Coral genetics strongly predicted survival of the ubiquitous coral, Porites (massive growth form), with variable survival (15 to 61%) across three morphologically indistinguishable-but genetically distinct-lineages. The heatwave also disrupted strong associations between these coral lineages and their algal symbionts (family Symbiodiniaceae), with symbiotic turnover in some colonies, resulting in reduced specificity across lineages. These results highlight how heatwaves can threaten cryptic genotypes and decouple otherwise tightly coevolved relationships between hosts and symbionts.

Moving beyond heritability in the search for coral adaptive potential

Global environmental change is happening at unprecedented rates. Coral reefs are among the ecosystems most threatened by global change. For wild populations to persist, they must adapt. Knowledge shortfalls about corals' complex ecological and evolutionary dynamics, however, stymie predictions about potential adaptation to future conditions. Here, we review adaptation through the lens of quantitative genetics. We argue that coral adaptation studies can benefit greatly from "wild" quantitative genetic methods, where traits are studied in wild populations undergoing natural selection, genomic relationship matrices can replace breeding experiments, and analyses can be extended to examine genetic constraints among traits. In addition, individuals with advantageous genotypes for anticipated future conditions can be identified. Finally, genomic genotyping supports simultaneous consideration of how genetic diversity is arrayed across geographic and environmental distances, providing greater context for predictions of phenotypic evolution at a metapopulation scale.

Cryptic diversity and spatial genetic variation in the coral Acropora tenuis and its endosymbionts across the Great Barrier Reef

Genomic studies are uncovering extensive cryptic diversity within reef-building corals, suggesting that evolutionarily and ecologically relevant diversity is highly underestimated in the very organisms that structure coral reefs. Furthermore, endosymbiotic algae within coral host species can confer adaptive responses to environmental stress and may represent additional axes of coral genetic variation that are not constrained by taxonomic divergence of the cnidarian host. Here, we examine genetic variation in a common and widespread, reef-building coral, Acropora tenuis, and its associated endosymbiotic algae along the entire expanse of the Great Barrier Reef (GBR). We use SNPs derived from genome-wide sequencing to characterize the cnidarian coral host and organelles from zooxanthellate endosymbionts (genus Cladocopium). We discover three distinct and sympatric genetic clusters of coral hosts, whose distributions appear associated with latitude and inshore-offshore reef position. Demographic modelling suggests that the divergence history of the three distinct host taxa ranges from 0.5 to 1.5 million years ago, preceding the GBR's formation, and has been characterized by low-to-moderate ongoing inter-taxon gene flow, consistent with occasional hybridization and introgression typifying coral evolution. Despite this differentiation in the cnidarian host, A. tenuis taxa share a common symbiont pool, dominated by the genus Cladocopium (Clade C). Cladocopium plastid diversity is not strongly associated with host identity but varies with reef location relative to shore: inshore colonies contain lower symbiont diversity on average but have greater differences between colonies as compared with symbiont communities from offshore colonies. Spatial genetic patterns of symbiont communities could reflect local selective pressures maintaining coral holobiont differentiation across an inshore-offshore environmental gradient. The strong influence of environment (but not host identity) on symbiont community composition supports the notion that symbiont community composition responds to habitat and may assist in the adaptation of corals to future environmental change.

A decade of population genetics studies of scleractinian corals: A systematic review

Coral reefs are the most diverse marine ecosystems. However, coral cover has decreased worldwide due to natural disturbances, climate change, and local anthropogenic drivers. In recent decades, various genetic methods and molecular markers have been developed to assess genetic diversity, structure, and connectivity in different coral species to determine the vulnerability of their populations. This review aims to identify population genetic studies of scleractinian corals in the last decade (2010-2020), and the techniques and molecular markers used. Bibliometric analysis was conducted to identify journals and authors working in this field. We then calculated the number of genetic studies by species and ecoregion based on data obtained from 178 studies found in Scopus and Web of Science. Coral Reefs and Molecular Ecology were the main journals published population genetics studies, and microsatellites are the most widely used molecular markers. The Caribbean, Australian Barrier Reef, and South Kuroshio in Japan are among the ecoregions with the most population genetics data. In contrast, we found limited information about the Coral Triangle, a region with the highest biodiversity and key to coral reef conservation. Notably, only 117 (out of 1500 described) scleractinian coral species have genetic studies. This review emphasizes which coral species have been studied and highlights remaining gaps and locations where such data is critical for coral conservation.

Divergence in interspecific and intersubspecific gene expression between two closely related horseshoe bats (Rhinolophus)

Closely related species have been used as representative systems to investigate the genetic mechanisms involved in the early stages of species differentiation. Previous studies have indicated that variation in gene expression might be a sensitive indicator of initial species divergence, although the role of expression divergence, and especially that associated with phenotypic variation remained relatively undefined. For three organs (cochlea, brain, and liver) from two closely related bat species (Rhinolophus siamensis and R. episcopus), the interspecific and intersubspecific gene expression profiles were compared using transcriptomics in this study. Striking organ specificity of expression was observed, and expression profiles exhibited similarities between cochlea and brain tissues. Numerous differentially expressed genes (DEGs) were identified for each organ in the interspecific comparison (cochlea/brain/liver: 1,069/647/692) and intersubspecific comparison (608/528/368). Functional enrichment analysis indicated vital variation in expression related to the immune system, ion activities, neuronal function, and multisensory system regulation in both comparisons. DEGs relevant to the variation in echolocation calls (RF) were found, and some of them were involved in the pivotal patterns of expression variation. The regulation of immune, ion channel, neural activity, and sophisticated sensory functions at the expression level might be key mechanisms in the early species divergence of bats, and the expression variation related to acoustical signal could have played a crucial part. This study expands our knowledge of gene expression and patterns of variation for three key organs to echolocation at both the interspecific and intersubspecific levels. Further, the framework described here provides insight into the genetic basis of phenotypic variation during the incipient stage of species differentiation.

Climatic similarity and genomic background shape the extent of parallel adaptation in Timema stick insects

Evolution can repeat itself, resulting in parallel adaptations in independent lineages occupying similar environments. Moreover, parallel evolution sometimes, but not always, uses the same genes. Two main hypotheses have been put forth to explain the probability and extent of parallel evolution. First, parallel evolution is more likely when shared ecologies result in similar patterns of natural selection in different taxa. Second, parallelism is more likely when genomes are similar because of shared standing variation and similar mutational effects in closely related genomes. Here we combine ecological, genomic, experimental and phenotypic data with Bayesian modelling and randomization tests to quantify the degree of parallelism and its relationship with ecology and genetics. Our results show that the extent to which genomic regions associated with climate are parallel among species of Timema stick insects is shaped collectively by shared ecology and genomic background. Specifically, the extent of genomic parallelism decays with divergence in climatic conditions (that is, habitat or ecological similarity) and genomic similarity. Moreover, we find that climate-associated loci are likely subject to selection in a field experiment, overlap with genetic regions associated with cuticular hydrocarbon traits and are not strongly shaped by introgression between species. Our findings shed light on when evolution is most expected to repeat itself.

The immediate costs and long-term benefits of assisted gene flow in large populations

With the genetic health of many plant and animal populations deteriorating due to climate change outpacing adaptation, interventions, such as assisted gene flow (AGF), may provide genetic variation necessary for populations to adapt to climate change. We ran genetic simulations to mimic different AGF scenarios in large populations and measured their outcomes on population-level fitness to determine circumstances in which it is worthwhile to perform AGF. In the absence of inbreeding depression, AGF was beneficial within a few generations only when introduced genotypes had much higher fitness than local individuals and traits affecting fitness were controlled by a few genes of large effect. AGF was harmful over short periods (e.g., first ∼10-20 generations) if there was strong outbreeding depression or introduced deleterious genetic variation. When the adaptive trait was controlled by many loci of small effect, the benefits of AGF took over 10 generations to realize-potentially too long for most climate-related management scenarios. The genomic integrity of the recipient population typically remained intact following AGF; the amount of genetic material from the donor population usually constituted no more of the recipient population's genome than the fraction of the population introduced. Significant genomic turnover (e.g., >50% replacement) only occurred when the selective advantage of the adaptive trait and translocation fraction were extremely high. Our results will be useful when adaptive management is used to maintain the genetic health and productivity of large populations under climate change.

Intrapopulation adaptive variance supports thermal tolerance in a reef-building coral

Coral holobionts are multi-species assemblages, which adds significant complexity to genotype-phenotype connections underlying ecologically important traits like coral bleaching. Small scale heterogeneity in bleaching is ubiquitous in the absence of strong environmental gradients, which provides adaptive variance needed for the long-term persistence of coral reefs. We used RAD-seq, qPCR and LC-MS/MS metabolomics to characterize host genomic variation, symbiont community and biochemical correlates in two bleaching phenotypes of the vertically transmitting coral Montipora capitata. Phenotype was driven by symbiosis state and host genetic variance. We documented 5 gene ontologies that were significantly associated with both the binary bleaching phenotype and symbiont composition, representing functions that confer a phenotype via host-symbiont interactions. We bred these corals and show that symbiont communities were broadly conserved in bulk-crosses, resulting in significantly higher survivorship under temperature stress in juveniles, but not larvae, from tolerant parents. Using a select and re-sequence approach, we document numerous gene ontologies selected by heat stress, some of which (cell signaling, antioxidant activity, pH regulation) have unique selection dynamics in larvae from thermally tolerant parents. These data show that vertically transmitting corals may have an adaptive advantage under climate change if host and symbiont variance interact to influence bleaching phenotype.

Genetic approaches for increasing fitness in endangered species

The global rate of wildlife extinctions is accelerating, and the persistence of many species requires conservation breeding programs. A central paradigm of these programs is to preserve the genetic diversity of the founder populations. However, this may preserve original characteristics that make them vulnerable to extinction. We introduce targeted genetic intervention (TGI) as an alternative approach that promotes traits that enable species to persist in the face of threats by changing the incidence of alleles that impact on fitness. The TGI toolkit includes methods with established efficacy in model organisms and agriculture but are largely untried for conservation, such as synthetic biology and artificial selection. We explore TGI approaches as a species-restoration tool for intractable threats including infectious disease and climate change.

Genomic analysis of distinct bleaching tolerances among cryptic coral species

Reef-building coral species are experiencing an unprecedented decline owing to increasing frequency and intensity of marine heatwaves and associated bleaching-induced mortality. Closely related species from the Acropora hyacinthus species complex differ in heat tolerance and in their association with heat-tolerant symbionts. We used low-coverage full genome sequencing of 114 colonies monitored across the 2015 bleaching event in American Samoa to determine the genetic differences among four cryptic species (termed HA, HC, HD and HE) that have diverged in these species traits. Cryptic species differed strongly at thousands of single nucleotide polymorphisms across the genome which are enriched for amino acid changes in the bleaching-resistant species HE. In addition, HE also showed two particularly divergent regions with strong signals of differentiation. One approximately 220 kb locus, HES1, contained the majority of fixed differences in HE. A second locus, HES2, was fixed in HE but polymorphic in the other cryptic species. Surprisingly, non-HE individuals with HE-like haplotypes at HES2 were more likely to bleach. At both loci, HE showed particular sequence similarity to a congener, Acropora millepora. Overall, resilience to bleaching during the third global bleaching event was strongly structured by host cryptic species, buoyed by differences in symbiont associations between these species.

Environmental specialization and cryptic genetic divergence in two massive coral species from the Florida Keys Reef Tract

Broadcast-spawning coral species have wide geographical ranges spanning strong environmental gradients, but it is unclear how much spatially varying selection these gradients actually impose. Strong divergent selection might present a considerable barrier for demographic exchange between disparate reef habitats. We investigated whether the cross-shelf gradient is associated with spatially varying selection in two common coral species, Montastraea cavernosa and Siderastrea siderea, in the Florida Keys. To this end, we generated a de novo genome assembly for M. cavernosa and used 2bRAD to genotype 20 juveniles and 20 adults of both species from each of the three reef zones to identify signatures of selection occurring within a single generation. Unexpectedly, each species was found to be composed of four genetically distinct lineages, with gene flow between them still ongoing but highly reduced in 13.0%-54.7% of the genome. Each species includes two sympatric lineages that are only found in the deep (20 m) habitat, while the other lineages are found almost exclusively on the shallower reefs (3-10 m). The two "shallow" lineages of M. cavernosa are also specialized for either nearshore or offshore: comparison between adult and juvenile cohorts indicates that cross-shelf migrants are more than twice as likely to die before reaching adulthood than local recruits. S. siderea and M. cavernosa are among the most ecologically successful species on the Florida Keys Reef Tract, and this work offers important insight into the genomic background of divergent selection and environmental specialization that may in part explain their resilience and broad environmental range.

Morphological stasis masks ecologically divergent coral species on tropical reefs

Coral reefs are the epitome of species diversity, yet the number of described scleractinian coral species, the framework-builders of coral reefs, remains moderate by comparison. DNA sequencing studies are rapidly challenging this notion by exposing a wealth of undescribed diversity, but the evolutionary and ecological significance of this diversity remains largely unclear. Here, we present an annotated genome for one of the most ubiquitous corals in the Indo-Pacific (Pachyseris speciosa) and uncover, through a comprehensive genomic and phenotypic assessment, that it comprises morphologically indistinguishable but ecologically divergent lineages. Demographic modeling based on whole-genome resequencing indicated that morphological crypsis (across micro- and macromorphological traits) was due to ancient morphological stasis rather than recent divergence. Although the lineages occur sympatrically across shallow and mesophotic habitats, extensive genotyping using a rapid molecular assay revealed differentiation of their ecological distributions. Leveraging "common garden" conditions facilitated by the overlapping distributions, we assessed physiological and quantitative skeletal traits and demonstrated concurrent phenotypic differentiation. Lastly, spawning observations of genotyped colonies highlighted the potential role of temporal reproductive isolation in the limited admixture, with consistent genomic signatures in genes related to morphogenesis and reproduction. Overall, our findings demonstrate the presence of ecologically and phenotypically divergent coral species without substantial morphological differentiation and provide new leads into the potential mechanisms facilitating such divergence. More broadly, they indicate that our current taxonomic framework for reef-building corals may be scratching the surface of the ecologically relevant diversity on coral reefs, consequently limiting our ability to protect or restore this diversity effectively.

Signals of selection beyond bottlenecks between exotic populations of the bull-headed dung beetle, Onthophagus taurus

Colonization of new environments can lead to population bottlenecks and rapid phenotypic evolution that could be due to neutral and selective processes. Exotic populations of the bull-headed dung beetle (Onthophagus taurus) have differentiated in opposite directions from native beetles in male horn-to-body size allometry and female fecundity. Here we test for genetic and transcriptional differences among two exotic and one native O. taurus populations after three generations in common garden conditions. We sequenced RNA from 24 individuals for each of the three populations including both sexes, and spanning four developmental stages for the two exotic, differentiated populations. Identifying 270,400 high-quality single nucleotide polymorphisms, we revealed a strong signal of genetic differentiation between the three populations, and evidence of recent bottlenecks within and an excess of outlier loci between exotic populations. Differences in gene expression between populations were greatest in prepupae and early adult life stages, stages during which differences in male horn development and female fecundity manifest. Finally, genes differentially expressed between exotic populations also had greater genetic differentiation and performed functions related to chitin biosynthesis and nutrient sensing, possibly underlying allometry and fecundity trait divergences. Our results suggest that beyond bottlenecks, recent introductions have led to genetic and transcriptional differences in genes correlated with observed phenotypic differences.

Genome-wide diversity and habitat underlie fine-scale phenotypic differentiation in the rainbow darter (Etheostoma caeruleum)

Adaptation to environmental change requires that populations harbor the necessary genetic variation to respond to selection. However, dispersal-limited species with fragmented populations and reduced genetic diversity may lack this variation and are at an increased risk of local extinction. In freshwater fish species, environmental change in the form of increased stream temperatures places many cold-water species at-risk. We present a study of rainbow darters (Etheostoma caeruleum) in which we evaluated the importance of genetic variation on adaptive potential and determined responses to extreme thermal stress. We compared fine-scale patterns of morphological and thermal tolerance differentiation across eight sites, including a unique lake habitat. We also inferred contemporary population structure using genomic data and characterized the relationship between individual genetic diversity and stress tolerance. We found site-specific variation in thermal tolerance that generally matched local conditions and morphological differences associated with lake-stream divergence. We detected patterns of population structure on a highly local spatial scale that could not be explained by isolation by distance or stream connectivity. Finally, we showed that individual thermal tolerance was positively correlated with genetic variation, suggesting that sites with increased genetic diversity may be better at tolerating novel stress. Our results highlight the importance of considering intraspecific variation in understanding population vulnerability and stress response.

Population genetics of the coral Acropora millepora: Toward genomic prediction of bleaching

Although reef-building corals are declining worldwide, responses to bleaching vary within and across species and are partly heritable. Toward predicting bleaching response from genomic data, we generated a chromosome-scale genome assembly for the coral Acropora millepora We obtained whole-genome sequences for 237 phenotyped samples collected at 12 reefs along the Great Barrier Reef, among which we inferred little population structure. Scanning the genome for evidence of local adaptation, we detected signatures of long-term balancing selection in the heat-shock co-chaperone sacsin We conducted a genome-wide association study of visual bleaching score for 213 samples, incorporating the polygenic score derived from it into a predictive model for bleaching in the wild. These results set the stage for genomics-based approaches in conservation strategies.

Meta-analysis of the coral environmental stress response: Acropora corals show opposing responses depending on stress intensity

As climate change progresses, reef-building corals must contend more often with suboptimal conditions, motivating a need to understand coral stress response. Here, we test the hypothesis that there is a stereotyped transcriptional response that corals enact under all stressful conditions, functionally characterized by downregulation of growth, and activation of cell death, response to reactive oxygen species, immunity, and protein folding and degradation. We analyse RNA-seq and Tag-Seq data from 14 previously published studies and supplement them with four new experiments involving different stressors, totaling over 600 gene expression profiles from the genus Acropora. Contrary to expectations, we found not one, but two distinct types of response. The type A response was observed under all kinds of high-intensity stress, was correlated between independent projects and was functionally consistent with the hypothesized stereotyped response. The consistent correlation between projects, irrespective of stress type, supports the type A response as the general coral environmental stress response (ESR), a blanket solution to severely stressful conditions. The distinct type B response was observed under lower intensity stress and was more variable among studies. Unexpectedly, at the level of individual genes and functional categories, the type B response was broadly opposite the type A response. Finally, taking advantage of the breadth of the data set, we present contextual annotations for previously unannotated genes based on consistent stress-induced differences across independent projects.

Genetic Characterization of Cupped Oyster Resources in Europe Using Informative Single Nucleotide Polymorphism (SNP) Panels

The Pacific oyster, Crassostrea gigas, was voluntarily introduced from Japan and British Columbia into Europe in the early 1970s, mainly to replace the Portuguese oyster, Crassostrea angulata, in the French shellfish industry, following a severe disease outbreak. Since then, the two species have been in contact in southern Europe and, therefore, have the potential to exchange genes. Recent evolutionary genomic works have provided empirical evidence that C. gigas and C. angulata exhibit partial reproductive isolation. Although hybridization occurs in nature, the rate of interspecific gene flow varies across the genome, resulting in highly heterogeneous genome divergence. Taking this biological property into account is important to characterize genetic ancestry and population structure in oysters. Here, we identified a subset of ancestry-informative makers from the most differentiated regions of the genome using existing genomic resources. We developed two different panels in order to (i) easily differentiate C. gigas and C. angulata, and (ii) describe the genetic diversity and structure of the cupped oyster with a particular focus on French Atlantic populations. Our results confirm high genetic homogeneity among Pacific cupped oyster populations in France and reveal several cases of introgressions between Portuguese and Japanese oysters in France and Portugal.

Local adaptation across a complex bioclimatic landscape in two montane bumble bee species

Understanding evolutionary responses to variation in temperature and precipitation across species ranges is of fundamental interest given ongoing climate change. The importance of temperature and precipitation for multiple aspects of bumble bee (Bombus) biology, combined with large geographic ranges that expose populations to diverse environmental pressures, make these insects well-suited for studying local adaptation. Here, we analyzed genome-wide sequence data from two widespread bumble bees, Bombus vosnesenskii and Bombus vancouverensis, using multiple environmental association analysis methods to investigate climate adaptation across latitude and altitude. The strongest signatures of selection were observed in B. vancouverensis, but despite unique responses between species for most loci, we detected several shared responses. Genes relating to neural and neuromuscular function and ion transport were especially evident with respect to temperature variables, while genes relating to cuticle formation, tracheal and respiratory system development, and homeostasis were associated with precipitation variables. Our data thus suggest that adaptive responses for tolerating abiotic variation are likely to be complex, but that several parallels among species can emerge even for these complex traits and landscapes. Results provide the framework for future work into mechanisms of thermal and desiccation tolerance in bumble bees and a set of genomic targets that might be monitored for future conservation efforts.

Resilience in reef-building corals: The ecological and evolutionary importance of the host response to thermal stress

Coral reefs are under extreme threat due to a number of stressors, but temperature increases due to changing climate are the most severe. Rising ocean temperatures coupled with local extremes lead to extensive bleaching, where the coral-algal symbiosis breaks down and corals may die, compromising the structure and function of reefs. Although the symbiotic nature of the coral colony has historically been a focus of research on coral resilience, the host itself is a foundational component in the response to thermal stress. Fixed effects in the coral host set trait baselines through evolutionary processes, acting on many loci of small effect to create mosaics of thermal tolerance across latitudes and individual coral reefs. These genomic differences can be strongly heritable, producing wide variation among clones of different genotypes or families of a specific larval cross. Phenotypic plasticity is overlaid on these baselines and a growing body of knowledge demonstrates the potential for acclimatization of reef-building corals through a variety of mechanisms that promote resilience and stress tolerance. The long-term persistence of coral reefs will require many of these mechanisms to adjust to warmer temperatures within a generation, bridging the gap to reproductive events that allow recombination of standing diversity and adaptive change. Business-as-usual climate scenarios will probably lead to the loss of some coral populations or species in the future, so the interaction between intragenerational effects and evolutionary pressure is critical for the survival of reefs.

Unknown to Known: Advancing Knowledge of Coral Gene Function

Given the catastrophic changes befalling coral reefs, understanding coral gene function is essential to advance reef conservation. This has proved challenging due to the paucity of genomic data and genetic tools available for corals. Recently, CRISPR/Cas9 gene editing was applied to these species; however, a major bottleneck is the identification and prioritization of candidate genes for manipulation. This issue is exacerbated by the many unknown ('dark') coral genes that may play key roles in the stress response. We review the use of gene coexpression networks that incorporate both known and unknown genes to identify targets for reverse genetic analysis. This approach also provides a framework for the annotation of dark genes in established interaction networks to improve our fundamental knowledge of coral gene function.

Regional heterogeneity impacts gene expression in the subarctic zooplankter Neocalanus flemingeri in the northern Gulf of Alaska

Marine pelagic species are being increasingly challenged by environmental change. Their ability to persist will depend on their capacity for physiological acclimatization. Little is known about limits of physiological plasticity in key species at the base of the food web. Here we investigate the capacity for acclimatization in the copepod Neocalanus flemingeri, which inhabits the Gulf of Alaska, a heterogeneous and highly seasonal environment. RNA-Seq analysis of field-collected pre-adults identified large regional differences in expression of genes involved in metabolic and developmental processes and response to stressors. We found that lipid synthesis genes were up-regulated in individuals from Prince William Sound and down-regulated in the Gulf of Alaska. Up-regulation of lipid catabolic genes in offshore individuals suggests they are experiencing nutritional deficits. The expression differences demonstrate physiological plasticity in response to a steep gradient in food availability. Our transcriptional analysis reveals mechanisms of acclimatization that likely contribute to the observed resilience of this population.

Genomic signatures of sympatric speciation with historical and contemporary gene flow in a tropical anthozoan (Hexacorallia: Actiniaria)

Sympatric diversification is recognized to have played an important role in the evolution of biodiversity. However, an in situ sympatric origin for codistributed taxa is difficult to demonstrate because different evolutionary processes can lead to similar biogeographic outcomes, especially in ecosystems that can readily facilitate secondary contact due to a lack of hard barriers to dispersal. Here we use a genomic (ddRADseq), model-based approach to delimit a species complex of tropical sea anemones that are codistributed on coral reefs throughout the Tropical Western Atlantic. We use coalescent simulations in fastsimcoal2 and ordinary differential equations in Moments to test competing diversification scenarios that span the allopatric-sympatric continuum. Our results suggest that the corkscrew sea anemone Bartholomea annulata is a cryptic species complex whose members are codistributed throughout their range. Simulation and model selection analyses from both approaches suggest these lineages experienced historical and contemporary gene flow, supporting a sympatric origin, but an alternative secondary contact model receives appreciable model support in fastsimcoal2. Leveraging the genome of the closely related Exaiptasia diaphana, we identify five loci under divergent selection between cryptic B. annulata lineages that fall within mRNA transcripts or CDS regions. Our study provides a rare empirical, genomic example of sympatric speciation in a tropical anthozoan and the first range-wide molecular study of a tropical sea anemone, underscoring that anemone diversity is under-described in the tropics, and highlighting the need for additional systematic studies into these ecologically and economically important species.

WhoGEM: an admixture-based prediction machine accurately predicts quantitative functional traits in plants

The explosive growth of genomic data provides an opportunity to make increased use of sequence variations for phenotype prediction. We have developed a prediction machine for quantitative phenotypes (WhoGEM) that overcomes some of the bottlenecks limiting the current methods. We demonstrated its performance by predicting quantitative disease resistance and quantitative functional traits in the wild model plant species, Medicago truncatula, using geographical locations as covariates for admixture analysis. The method's prediction reliability equals or outperforms all existing algorithms for quantitative phenotype prediction. WhoGEM analysis produces evidence that variation in genome admixture proportions explains most of the phenotypic variation for quantitative phenotypes.

Evolutionary genomics of gypsy moth populations sampled along a latitudinal gradient

The European gypsy moth (Lymantria dispar L.) was first introduced to Massachusetts in 1869 and within 150 years has spread throughout eastern North America. This large-scale invasion across a heterogeneous landscape allows examination of the genetic signatures of adaptation potentially associated with rapid geographical spread. We tested the hypothesis that spatially divergent natural selection has driven observed changes in three developmental traits that were measured in a common garden for 165 adult moths sampled from six populations across a latitudinal gradient covering the entirety of the range. We generated genotype data for 91,468 single nucleotide polymorphisms based on double digest restriction-site associated DNA sequencing and used these data to discover genome-wide associations for each trait, as well as to test for signatures of selection on the discovered architectures. Genetic structure across the introduced range of gypsy moth was low in magnitude (F_ST  = 0.069), with signatures of bottlenecks and spatial expansion apparent in the rare portion of the allele frequency spectrum. Results from applications of Bayesian sparse linear mixed models were consistent with the presumed polygenic architectures of each trait. Further analyses indicated spatially divergent natural selection acting on larval development time and pupal mass, with the linkage disequilibrium component of this test acting as the main driver of observed patterns. The populations most important for these signals were two range-edge populations established less than 30 generations ago. We discuss the importance of rapid polygenic adaptation to the ability of non-native species to invade novel environments.

Tolerance traits related to climate change resilience are independent and polygenic

The resilience of organisms to climate change through adaptive evolution is dependent on the extent of genetically based variation in key phenotypic traits and the nature of genetic associations between them. For aquatic animals, upper thermal tolerance and hypoxia tolerance are likely to be a important determinants of sensitivity to climate change. To determine the genetic basis of these traits and to detect associations between them, we compared naturally occurring populations of two subspecies of Atlantic killifish, Fundulus heteroclitus, that differ in both thermal and hypoxia tolerance. Multilocus association mapping demonstrated that 47 and 35 single nucleotide polymorphisms (SNPs) explained 43.4% and 51.9% of variation in thermal and hypoxia tolerance, respectively, suggesting that genetic mechanisms underlie a substantial proportion of variation in each trait. However, no explanatory SNPs were shared between traits, and upper thermal tolerance varied approximately linearly with latitude, whereas hypoxia tolerance exhibited a steep phenotypic break across the contact zone between the subspecies. These results suggest that upper thermal tolerance and hypoxia tolerance are neither phenotypically correlated nor genetically associated, and thus that rates of adaptive change in these traits can be independently fine-tuned by natural selection. This modularity of important traits can underpin the evolvability of organisms to complex future environmental change.

Analysis of Genome-Wide Differentiation between Native and Introduced Populations of the Cupped Oysters Crassostrea gigas and Crassostrea angulata

The Pacific cupped oyster is genetically subdivided into two sister taxa, Crassostrea gigas and Crassostrea angulata, which are in contact in the north-western Pacific. The nature and origin of their genetic and taxonomic differentiation remains controversial due the lack of known reproductive barriers and the high degree of morphologic similarity. In particular, whether the presence of ecological and/or intrinsic isolating mechanisms contributes to species divergence is unknown. The recent co-introduction of both taxa into Europe offers a unique opportunity to test how genetic differentiation is maintained under new environmental and demographic conditions. We generated a pseudochromosome assembly of the Pacific oyster genome using a combination of BAC-end sequencing and scaffold anchoring to a new high-density linkage map. We characterized genome-wide differentiation between C. angulata and C. gigas in both their native and introduced ranges, and showed that gene flow between species has been facilitated by their recent co-introductions in Europe. Nevertheless, patterns of genomic divergence between species remain highly similar in Asia and Europe, suggesting that the environmental transition caused by the co-introduction of the two species did not affect the genomic architecture of their partial reproductive isolation. Increased genetic differentiation was preferentially found in regions of low recombination. Using historical demographic inference, we show that the heterogeneity of differentiation across the genome is well explained by a scenario whereby recent gene flow has eroded past differentiation at different rates across the genome after a period of geographical isolation. Our results thus support the view that low-recombining regions help in maintaining intrinsic genetic differences between the two species.