Mixed Patterns of Intergenerational DNA Methylation Inheritance in Acropora.

For sessile organisms at high risk from climate change, phenotypic plasticity can be critical to rapid acclimation. Epigenetic markers like DNA methylation are hypothesized as mediators of plasticity; methylation is associated with the regulation of gene expression, can change in response to ecological cues, and is a proposed basis for the inheritance of acquired traits. Within reef-building corals, gene-body methylation (gbM) can change in response to ecological stressors. If coral DNA methylation is transmissible across generations, this could potentially facilitate rapid acclimation to environmental change. We investigated methylation heritability in Acropora, a stony reef-building coral. Two Acropora millepora and two Acropora selago adults were crossed, producing eight offspring crosses (four hybrid, two of each species). We used whole-genome bisulfite sequencing to identify methylated loci and allele-specific alignments to quantify per-locus inheritance. If methylation is heritable, differential methylation (DM) between the parents should equal DM between paired offspring alleles at a given locus. We found a mixture of heritable and nonheritable loci, with heritable portions ranging from 44% to 90% among crosses. gBM was more heritable than intergenic methylation, and most loci had a consistent degree of heritability between crosses (i.e. the deviation between parental and offspring DM were of similar magnitude and direction). Our results provide evidence that coral methylation can be inherited but that heritability is heterogenous throughout the genome. Future investigations into this heterogeneity and its phenotypic implications will be important to understanding the potential capability of intergenerational environmental acclimation in reef building corals.

Should I stay or should I go? Coral bleaching from the symbionts' perspective.

Coral bleaching, the stress-induced breakdown of coral-algal symbiosis, threatens reefs globally. Paradoxically, despite adverse fitness effects, corals bleach annually, even outside of abnormal temperatures. This generally occurs shortly after the once-per-year mass coral spawning. Here, we propose a hypothesis linking annual coral bleaching and the transmission of symbionts to the next generation of coral hosts. We developed a dynamic model with two symbiont growth strategies, and found that high sexual recruitment and low adult coral survivorship and growth favour bleaching susceptibility, while the reverse promotes bleaching resilience. Otherwise, unexplained trends in the Indo-Pacific align with our hypothesis, where reefs and coral taxa exhibiting higher recruitment are more bleaching susceptible. The results from our model caution against interpreting potential shifts towards more bleaching-resistant symbionts as evidence of climate adaptation-we predict such a shift could also occur in declining systems experiencing low recruitment rates, a common scenario on today's reefs.

Evidence of sweepstakes reproductive success in a broadcast-spawning coral and its implications for coral metapopulation persistence.

Processes governing genetic diversity and adaptive potential in reef-building corals are of interest both for fundamental evolutionary biology and for reef conservation. Here, we investigated the possibility of "sweepstakes reproductive success" (SRS) in a broadcast spawning coral, Acropora hyacinthus, at Yap Island, Micronesia. SRS is an extreme yearly variation in the number of surviving offspring among parents. It is predicted to generate genetically differentiated, low-genetic-diversity recruit cohorts, containing close kin individuals. We have tested these predictions by comparing genetic composition of size classes (adults and juveniles) at several sites on the island of Yap. We did see the genome-wide dip in genetic diversity in juveniles compared to adults at two of the four sites; however, both adults and juveniles varied in genetic diversity across sites, and there was no detectable genetic structure among juveniles, which does not conform to the classical SRS scenario. Yet, we have identified a pair of juvenile siblings at the site where juveniles had the lowest genetic diversity compared to adults, an observation that is hard to explain without invoking SRS. While further support for SRS is needed to fully settle the issue, we show that incorporating SRS into the Indo-West Pacific coral metapopulation adaptation model had surprisingly little effect on mean rates of coral cover decline during warming. Still, SRS notably increases year-to-year variation in coral cover throughout the simulation.

Gene expression underlying variation in the survival skills of red drum larvae (Sciaenops ocellatus).

Mortality rates of marine fish larvae are incredibly high and can determine year-class strength. The major causes of larval mortality are predation and starvation, and the performance of larvae in survival skills that can mitigate this mortality (predator evasion, foraging) varies among individuals and cohorts, but the causes of the variation are not known. Transcriptomics can link gene expression variation to phenotypic variation at the whole-system level to investigate the molecular basis of behavioural variation. We used tag-based RNA-sequencing to examine the molecular basis of variation in predator evasion and routine swimming (trait related to foraging efficiency) in the larval red drum, Sciaenops ocellatus. We looked for functional gene networks in which interindividual variation would explain variation in larval behavioural performance. We identified co-expressed gene groups ("modules") associated with predator evasion traits and found enrichment of motor, neural and energy metabolism pathways. These functional associations and pattern of correlations between modules and traits suggest that energy availability and allocation were responsible for the magnitude of startle responses, while differential neural and motor activation were associated with differences in response latency.

Changes in gene body methylation do not correlate with changes in gene expression in Anthozoa or Hexapoda.

BACKGROUND: As human activity alters the planet, there is a pressing need to understand how organisms adapt to environmental change. Of growing interest in this area is the role of epigenetic modifications, such as DNA methylation, in tailoring gene expression to fit novel conditions. Here, we reanalyzed nine invertebrate (Anthozoa and Hexapoda) datasets to validate a key prediction of this hypothesis: changes in DNA methylation in response to some condition correlate with changes in gene expression. RESULTS: In accord with previous observations, baseline levels of gene body methylation (GBM) positively correlated with transcription, and negatively correlated with transcriptional variation between conditions. Correlations between changes in GBM and transcription, however, were negligible. There was also no consistent negative correlation between methylation and transcription at the level of gene body methylation class (either highly- or lowly-methylated), anticipated under the previously described "seesaw hypothesis". CONCLUSION: Our results do not support the direct involvement of GBM in regulating dynamic transcriptional responses in invertebrates. If changes in DNA methylation regulate invertebrate transcription, the mechanism must involve additional factors or regulatory influences.

Millennia-old coral holobiont DNA provides insight into future adaptive trajectories.

Ancient DNA (aDNA) has been applied to evolutionary questions across a wide variety of taxa. Here, for the first time, we utilized aDNA from millennia-old fossil coral fragments to gain new insights into a rapidly declining western Atlantic reef ecosystem. We sampled four Acropora palmata fragments (dated 4215 BCE to 1099 CE) obtained from two Florida Keys reef cores. From these samples, we established that it is possible both to sequence aDNA from reef cores and place the data in the context of modern-day genetic variation. We recovered varying amounts of nuclear DNA exhibiting the characteristic signatures of aDNA from the A. palmata fragments. To describe the holobiont sensu lato, which plays a crucial role in reef health, we utilized metagenome-assembled genomes as a reference to identify a large additional proportion of ancient microbial DNA from the samples. The samples shared many common microbes with modern-day coral holobionts from the same region, suggesting remarkable holobiont stability over time. Despite efforts, we were unable to recover ancient Symbiodiniaceae reads from the samples. Comparing the ancient A. palmata data to whole-genome sequencing data from living acroporids, we found that while slightly distinct, ancient samples were most closely related to individuals of their own species. Together, these results provide a proof-of-principle showing that it is possible to carry out direct analysis of coral holobiont change over time, which lays a foundation for studying the impacts of environmental stress and evolutionary constraints.

Coral-bleaching responses to climate change across biological scales.

The global impacts of climate change are evident in every marine ecosystem. On coral reefs, mass coral bleaching and mortality have emerged as ubiquitous responses to ocean warming, yet one of the greatest challenges of this epiphenomenon is linking information across scientific disciplines and spatial and temporal scales. Here we review some of the seminal and recent coral-bleaching discoveries from an ecological, physiological, and molecular perspective. We also evaluate which data and processes can improve predictive models and provide a conceptual framework that integrates measurements across biological scales. Taking an integrative approach across biological and spatial scales, using for example hierarchical models to estimate major coral-reef processes, will not only rapidly advance coral-reef science but will also provide necessary information to guide decision-making and conservation efforts. To conserve reefs, we encourage implementing mesoscale sanctuaries (thousands of km2 ) that transcend national boundaries. Such networks of protected reefs will provide reef connectivity, through larval dispersal that transverse thermal environments, and genotypic repositories that may become essential units of selection for environmentally diverse locations. Together, multinational networks may be the best chance corals have to persist through climate change, while humanity struggles to reduce emissions of greenhouse gases to net zero.

Comparative neurotranscriptomics reveal widespread species differences associated with bonding.

BACKGROUND: Pair bonding with a reproductive partner is rare among mammals but is an important feature of human social behavior. Decades of research on monogamous prairie voles (Microtus ochrogaster), along with comparative studies using the related non-bonding meadow vole (M. pennsylvanicus), have revealed many of the neural and molecular mechanisms necessary for pair-bond formation in that species. However, these studies have largely focused on just a few neuromodulatory systems. To test the hypothesis that neural gene expression differences underlie differential capacities to bond, we performed RNA-sequencing on tissue from three brain regions important for bonding and other social behaviors across bond-forming prairie voles and non-bonding meadow voles. We examined gene expression in the amygdala, hypothalamus, and combined ventral pallidum/nucleus accumbens in virgins and at three time points after mating to understand species differences in gene expression at baseline, in response to mating, and during bond formation. RESULTS: We first identified species and brain region as the factors most strongly associated with gene expression in our samples. Next, we found gene categories related to cell structure, translation, and metabolism that differed in expression across species in virgins, as well as categories associated with cell structure, synaptic and neuroendocrine signaling, and transcription and translation that varied among the focal regions in our study. Additionally, we identified genes that were differentially expressed across species after mating in each of our regions of interest. These include genes involved in regulating transcription, neuron structure, and synaptic plasticity. Finally, we identified modules of co-regulated genes that were strongly correlated with brain region in both species, and modules that were correlated with post-mating time points in prairie voles but not meadow voles. CONCLUSIONS: These results reinforce the importance of pre-mating differences that confer the ability to form pair bonds in prairie voles but not promiscuous species such as meadow voles. Gene ontology analysis supports the hypothesis that pair-bond formation involves transcriptional regulation, and changes in neuronal structure. Together, our results expand knowledge of the genes involved in the pair bonding process and open new avenues of research in the molecular mechanisms of bond formation.

Fantastic Beasts and How To Sequence Them: Ecological Genomics for Obscure Model Organisms.

The application of genomic approaches to 'obscure model organisms' (OMOs), meaning species with no prior genomic resources, enables increasingly sophisticated studies of the genomic basis of evolution, acclimatization, and adaptation in real ecological contexts. I consider here ecological questions that can be addressed using OMOs, and indicate optimal sequencing and data-handling solutions for each case. With this I hope to promote the diversity of OMO-based projects that would capitalize on the peculiarities of the natural history of OMOs and could feasibly be completed within the scope of a single PhD thesis.

Shuffling between Cladocopium and Durusdinium extensively modifies the physiology of each symbiont without stressing the coral host.

As sea surface temperatures increase, many coral species that used to harbour symbionts of the genus Cladocopium have become colonized with the thermally tolerant genus, Durusdinium. Here, we asked how gene expression in the symbionts of one genus changes depending on the abundance of another symbiont genus within the same coral host, and what effect this interaction has on the host. Symbiont gene expression was overwhelmingly driven by whether the genus was the minority or the majority within the host, which affected 79% (Durusdinium) and 96% (Cladocopium) of all genes. Particularly strong effects in both genera were observed for photosynthesis components (upregulated in the minority state) and proteins putatively associated with cell motility (upregulated in the majority state). Importantly, there was no distinct gene expression signature associated with the mixed symbiosis state when both genera were represented in comparable proportions within the host, which could lead to more intense competition. The mixed symbiosis was also not associated with elevated host stress: in fact, after heat treatment, stress signatures were the lowest in mixed-symbiosis corals compared to both Cladocopium- and Durusdinium-dominated corals. In conclusion, during shuffling between Cladocopium and Durusdinium both symbiont genera go through extensive and largely reciprocal physiological transitions, but there is no evidence of intensifying antagonistic interactions that are detrimental to the host. Unless the mixed-symbiosis corals in this study are not representative of the typical transition between Cladocopium and Durusdinium, the process of shuffling from one symbiont genus to another appears to be cost-free for the coral host, and even appears to be associated with lower stress susceptibility. This raises optimism for the future corals, which will probably have to rely on symbiont shuffling more and more to withstand environmental challenges.

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.

How mitonuclear discordance and geographic variation have confounded species boundaries in a widely studied snake.

As DNA sequencing technologies and methods for delimiting species with genomic data become more accessible and numerous, researchers have more tools than ever to investigate questions in systematics and phylogeography. However, easy access to sophisticated computational tools is not without its drawbacks. Choosing the right approach for one's question can be challenging when presented with multitudinous options, some of which fail to distinguish between species and intraspecific population structure. Here, we employ a methodology that emphasizes intensive geographic sampling, particularly at contact zones between populations, with a focus on differentiating intraspecific genetic clusters from species in the Pantherophis guttatus complex, a group of North American ratsnakes. Using a mitochondrial marker as well as ddRADseq data, we find evidence of mitonuclear discordance which has contributed to historical confusion about the relationships within this group. Additionally, we identify geographically and genetically structured populations within the species Pantherophis emoryi that are congruent with previously described morphological variation. Importantly, we find that these structured populations within P. emoryi are highly admixed throughout the range of the species and show no evidence of any reproductive isolation. Our data support a revision of the taxonomy of this group, and we recognize two species within the complex and three subspecies within P. emoryi. This study illustrates the importance of thorough sampling of contact zones and consideration of gene flow when delimiting species in widespread complexes containing parapatric lineages.

Potential and limits for rapid genetic adaptation to warming in a Great Barrier Reef coral.

Can genetic adaptation in reef-building corals keep pace with the current rate of sea surface warming? Here we combine population genomics, biophysical modeling, and evolutionary simulations to predict future adaptation of the common coral Acropora millepora on the Great Barrier Reef (GBR). Genomics-derived migration rates were high (0.1-1% of immigrants per generation across half the latitudinal range of the GBR) and closely matched the biophysical model of larval dispersal. Both genetic and biophysical models indicated the prevalence of southward migration along the GBR that would facilitate the spread of heat-tolerant alleles to higher latitudes as the climate warms. We developed an individual-based metapopulation model of polygenic adaptation and parameterized it with population sizes and migration rates derived from the genomic analysis. We find that high migration rates do not disrupt local thermal adaptation, and that the resulting standing genetic variation should be sufficient to fuel rapid region-wide adaptation of A. millepora populations to gradual warming over the next 20-50 coral generations (100-250 years). Further adaptation based on novel mutations might also be possible, but this depends on the currently unknown genetic parameters underlying coral thermal tolerance and the rate of warming realized. Despite this capacity for adaptation, our model predicts that coral populations would become increasingly sensitive to random thermal fluctuations such as ENSO cycles or heat waves, which corresponds well with the recent increase in frequency of catastrophic coral bleaching events.

Role of gene body methylation in acclimatization and adaptation in a basal metazoan.

Gene body methylation (GBM) has been hypothesized to modulate responses to environmental change, including transgenerational plasticity, but the evidence thus far has been lacking. Here we show that coral fragments reciprocally transplanted between two distant reefs respond predominantly by increase or decrease in genome-wide GBM disparity: The range of methylation levels between lowly and highly methylated genes becomes either wider or narrower. Remarkably, at a broad functional level this simple adjustment correlated very well with gene expression change, reflecting a shifting balance between expressions of environmentally responsive and housekeeping genes. In our experiment, corals in a lower-quality habitat up-regulated genes involved in environmental responses, while corals in a higher-quality habitat invested more in housekeeping genes. Transplanted fragments showing closer GBM match to local corals attained higher fitness characteristics, which supports GBM's role in acclimatization. Fixed differences in GBM between populations did not align with plastic GBM changes and were mostly observed in genes with elevated FST, which suggests that they arose predominantly through genetic divergence. However, we cannot completely rule out transgenerational inheritance of acquired GBM states.

CRISPR/Cas9-mediated genome editing in a reef-building coral.

Reef-building corals are critically important species that are threatened by anthropogenic stresses including climate change. In attempts to understand corals' responses to stress and other aspects of their biology, numerous genomic and transcriptomic studies have been performed, generating a variety of hypotheses about the roles of particular genes and molecular pathways. However, it has not generally been possible to test these hypotheses rigorously because of the lack of genetic tools for corals. Here, we demonstrate efficient genome editing using the CRISPR/Cas9 system in the coral Acropora millepora We targeted the genes encoding fibroblast growth factor 1a (FGF1a), green fluorescent protein (GFP), and red fluorescent protein (RFP). After microinjecting CRISPR/Cas9 ribonucleoprotein complexes into fertilized eggs, we detected induced mutations in the targeted genes using changes in restriction-fragment length, Sanger sequencing, and high-throughput Illumina sequencing. We observed mutations in ∼50% of individuals screened, and the proportions of wild-type and various mutant gene copies in these individuals indicated that mutation induction continued for at least several cell cycles after injection. Although multiple paralogous genes encoding green fluorescent proteins are present in A. millepora, appropriate design of the guide RNA allowed us to induce mutations simultaneously in more than one paralog. Because A. millepora larvae can be induced to settle and begin colony formation in the laboratory, CRISPR/Cas9-based gene editing should allow rigorous tests of gene function in both larval and adult corals.

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.

Cladocopium community divergence in two Acropora coral hosts across multiple spatial scales.

Many broadly-dispersing corals acquire their algal symbionts (Symbiodiniaceae) "horizontally" from their environment upon recruitment. Horizontal transmission could promote coral fitness across diverse environments provided that corals can associate with divergent algae across their range and that these symbionts exhibit reduced dispersal potential. Here we quantified community divergence of Cladocopium algal symbionts in two coral host species (Acropora hyacinthus, Acropora digitifera) across two spatial scales (reefs on the same island, and between islands) across the Micronesian archipelago using microsatellites. We find that both hosts associated with a variety of multilocus genotypes (MLG) within two genetically distinct Cladocopium lineages (C40, C21), confirming that Acropora coral hosts associate with a range of Cladocopium symbionts across this region. Both C40 and C21 included multiple asexual lineages bearing identical MLGs, many of which spanned host species, reef sites within islands, and even different islands. Both C40 and C21 exhibited moderate host specialization and divergence across islands. In addition, within every island, algal symbiont communities were significantly clustered by both host species and reef site, highlighting that coral-associated Cladocopium communities are structured across small spatial scales and within hosts on the same reef. This is in stark contrast to their coral hosts, which never exhibited significant genetic divergence between reefs on the same island. These results support the view that horizontal transmission could improve local fitness for broadly dispersing Acropora coral species.

Estimating the potential for coral adaptation to global warming across the Indo-West Pacific.

The potential of reef-building corals to adapt to increasing sea-surface temperatures is often debated but has rarely been comprehensively modeled on a region-wide scale. We used individual-based simulations to model adaptation to warming in a coral metapopulation comprising 680 reefs and representing the whole of the Central Indo-West Pacific. Encouragingly, some reefs-most notably Vietnam, Japan, Taiwan, New Caledonia and the southern half of the Great Barrier Reef-exhibited high capacity for adaptation and, in our model, maintained coral cover even under a rapid "business-as-usual" warming scenario throughout the modeled period (200 years). Higher resilience of these reefs was observed under all tested parameter settings except the models prohibiting selection and/or migration during warming. At the same time, the majority of reefs in the region tended to collapse within the first 100 years of warming. The adaptive potential (odds of maintaining high coral cover) of a given reef could be predicted based on two metrics: the reef's present-day temperature, and the proportion of recruits immigrating from warmer locations. The latter metric explains the most variation in adaptive potential, and significantly correlates with actual coral cover changes observed throughout the region between the 1970s and the early 2000s. These findings will help prioritize coral conservation efforts and plan assisted gene flow interventions to boost the adaptive potential of specific coral populations.

Positive genetic associations among fitness traits support evolvability of a reef-building coral under multiple stressors.

Climate change threatens organisms in a variety of interactive ways that requires simultaneous adaptation of multiple traits. Predicting evolutionary responses requires an understanding of the potential for interactions among stressors and the genetic variance and covariance among fitness-related traits that may reinforce or constrain an adaptive response. Here we investigate the capacity of Acropora millepora, a reef-building coral, to adapt to multiple environmental stressors: rising sea surface temperature, ocean acidification, and increased prevalence of infectious diseases. We measured growth rates (weight gain), coral color (a proxy for Symbiodiniaceae density), and survival, in addition to nine physiological indicators of coral and algal health in 40 coral genets exposed to each of these three stressors singly and combined. Individual stressors resulted in predicted responses (e.g., corals developed lesions after bacterial challenge and bleached under thermal stress). However, corals did not suffer substantially more when all three stressors were combined. Nor were trade-offs observed between tolerances to different stressors; instead, individuals performing well under one stressor also tended to perform well under every other stressor. An analysis of genetic correlations between traits revealed positive covariances, suggesting that selection to multiple stressors will reinforce rather than constrain the simultaneous evolution of traits related to holobiont health (e.g., weight gain and algal density). These findings support the potential for rapid coral adaptation under climate change and emphasize the importance of accounting for corals' adaptive capacity when predicting the future of coral reefs.

Role of host genetics and heat-tolerant algal symbionts in sustaining populations of the endangered coral Orbicella faveolata in the Florida Keys with ocean warming.

Identifying which factors lead to coral bleaching resistance is a priority given the global decline of coral reefs with ocean warming. During the second year of back-to-back bleaching events in the Florida Keys in 2014 and 2015, we characterized key environmental and biological factors associated with bleaching resilience in the threatened reef-building coral Orbicella faveolata. Ten reefs (five inshore, five offshore, 179 corals total) were sampled during bleaching (September 2015) and recovery (May 2016). Corals were genotyped with 2bRAD and profiled for algal symbiont abundance and type. O. faveolata at the inshore sites, despite higher temperatures, demonstrated significantly higher bleaching resistance and better recovery compared to offshore. The thermotolerant Durusdinium trenchii (formerly Symbiondinium trenchii) was the dominant endosymbiont type region-wide during initial (78.0% of corals sampled) and final (77.2%) sampling; >90% of the nonbleached corals were dominated by D. trenchii. 2bRAD host genotyping found no genetic structure among reefs, but inshore sites showed a high level of clonality. While none of the measured environmental parameters were correlated with bleaching, 71% of variation in bleaching resistance and 73% of variation in the proportion of D. trenchii was attributable to differences between genets, highlighting the leading role of genetics in shaping natural bleaching patterns. Notably, D. trenchii was rarely dominant in O. faveolata from the Florida Keys in previous studies, even during bleaching. The region-wide high abundance of D. trenchii was likely driven by repeated bleaching associated with the two warmest years on record for the Florida Keys (2014 and 2015). On inshore reefs in the Upper Florida Keys, O. faveolata was most abundant, had the highest bleaching resistance, and contained the most corals dominated by D. trenchii, illustrating a causal link between heat tolerance and ecosystem resilience with global change.