RADseq population genomics confirms divergence across closely related species in blue coral (Heliopora coerulea).

BACKGROUND: Heliopora coerulea, the blue coral, is the octocoral characterized by its blue skeleton. Recently, two Heliopora species were delimited by DNA markers: HC-A and HC-B. To clarify the genomic divergence of these Heliopora species (HC-A and HC-B) from sympatric and allopatric populations in Okinawa, Japan, we used a high throughput reduced representation genomic DNA sequencing approach (ezRAD). RESULTS: We found 6742 biallelic SNPs shared among all target populations, which successfully distinguished the HC-A and HC-B species in both the sympatric and allopatric populations, with no evidence of hybridization between the two. In addition, we detected 410 fixed SNPs linking functional gene differences, including heat resilience and reproductive timing, between HC-A and HC-B. CONCLUSIONS: We confirmed clear genomic divergence between Heliopora species and found possible genes related to stress-responses and reproduction, which may shed light on the speciation process and ecological divergence of coral species.

Rare coral under the genomic microscope: timing and relationships among Hawaiian Montipora.

BACKGROUND: Evolutionary patterns of scleractinian (stony) corals are difficult to infer given the existence of few diagnostic characters and pervasive phenotypic plasticity. A previous study of Hawaiian Montipora (Scleractinia: Acroporidae) based on five partial mitochondrial and two nuclear genes revealed the existence of a species complex, grouping one of the rarest known species (M. dilatata, which is listed as Endangered by the International Union for Conservation of Nature - IUCN) with widespread corals of very different colony growth forms (M. flabellata and M. cf. turgescens). These previous results could result from a lack of resolution due to a limited number of markers, compositional heterogeneity or reflect biological processes such as incomplete lineage sorting (ILS) or introgression. RESULTS: All 13 mitochondrial protein-coding genes from 55 scleractinians (14 lineages from this study) were used to evaluate if a recent origin of the M. dilatata species complex or rate heterogeneity could be compromising phylogenetic inference. Rate heterogeneity detected in the mitochondrial data set seems to have no significant impacts on the phylogenies but clearly affects age estimates. Dating analyses show different estimations for the speciation of M. dilatata species complex depending on whether taking compositional heterogeneity into account (0.8 [0.05-2.6] Myr) or assuming rate homogeneity (0.4 [0.14-0.75] Myr). Genomic data also provided evidence of introgression among all analysed samples of the complex. RADseq data indicated that M. capitata colour morphs may have a genetic basis. CONCLUSIONS: Despite the volume of data (over 60,000 SNPs), phylogenetic relationships within the M. dilatata species complex remain unresolved most likely due to a recent origin and ongoing introgression. Species delimitation with genomic data is not concordant with the current taxonomy, which does not reflect the true diversity of this group. Nominal species within the complex are either undergoing a speciation process or represent ecomorphs exhibiting phenotypic polymorphisms.

Divergence times in demosponges (Porifera): first insights from new mitogenomes and the inclusion of fossils in a birth-death clock model.

BACKGROUND: Approximately 80% of all described extant sponge species belong to the class Demospongiae. Yet, despite their diversity and importance, accurate divergence times are still unknown for most demosponge clades. The estimation of demosponge divergence time is key to answering fundamental questions on the origin of Demospongiae, their diversification and historical biogeography. Molecular sequence data alone is not informative on an absolute time scale, and therefore needs to be "calibrated" with additional data such as fossils. Here, we calibrate the molecular data with the fossilized birth-death model, which compared to strict node dating, allows for the inclusion of young and old fossils in the analysis of divergence time. We use desma-bearing sponges, a diverse group of demosponges that form rigid skeletons and have a rich and continuous fossil record dating back to the Cambrian (~500 Ma), to date the demosponge radiation and constrain the timing of key evolutionary events, like the transition from marine to freshwater habitats. To infer a dated phylogeny of Demospongiae we assembled the mitochondrial genomes of six desma-bearing demosponges from reduced-representation genomic libraries. The total dataset included 33 complete demosponge mitochondrial genomes and 30 fossils. RESULTS: Our study supports a Neoproterozoic origin of Demospongiae. Novel age estimates for the split of freshwater and marine sponges dating back to the Carboniferous and the previously assumed recent (~18 Ma) diversification of freshwater sponges is supported. Moreover, we provide detailed age estimates for a possible diversification of Tetractinellidae (~315 Ma), the Astrophorina (~240 Ma), the Spirophorina (~120 Ma) and the family Corallistidae (~188 Ma) all of which are considered as key groups for dating the Demospongiae due to their extraordinary rich and continuous fossil history. CONCLUSION: This study provides novel insights into the evolution of Demospongiae. Observed discrepancies of our dated phylogeny with their putative first fossil appearance dates are discussed for selected sponge groups. For instance, a Carboniferous origin of the order Tetractinellida seems to be too late, compared to their first appearance in the fossil record in the Middle Cambrian. This would imply that Paleozoic spicule forms are not homologous to post-Paleozoic forms.

Complex signatures of genomic variation of two non-model marine species in a homogeneous environment.

BACKGROUND: Genomic tools are increasingly being used on non-model organisms to provide insights into population structure and variability, including signals of selection. However, most studies are carried out in regions with distinct environmental gradients or across large geographical areas, in which local adaptation is expected to occur. Therefore, the focus of this study is to characterize genomic variation and selective signals over short geographic areas within a largely homogeneous region. To assess adaptive signals between microhabitats within the rocky shore, we compared genomic variation between the Cape urchin (Parechinus angulosus), which is a low to mid-shore species, and the Granular limpet (Scutellastra granularis), a high shore specialist. RESULTS: Using pooled restriction site associated DNA (RAD) sequencing, we described patterns of genomic variation and identified outlier loci in both species. We found relatively low numbers of outlier SNPs within each species, and identified outlier genes associated with different selective pressures than those previously identified in studies conducted over larger environmental gradients. The number of population-specific outlier loci differed between species, likely owing to differential selective pressures within the intertidal environment. Interestingly, the outlier loci were highly differentiated within the two northernmost populations for both species, suggesting that unique evolutionary forces are acting on marine invertebrates within this region. CONCLUSIONS: Our study provides a background for comparative genomic studies focused on non-model species, as well as a baseline for the adaptive potential of marine invertebrates along the South African west coast. We also discuss the caveats associated with Pool-seq and potential biases of sequencing coverage on downstream genomic metrics. The findings provide evidence of species-specific selective pressures within a homogeneous environment, and suggest that selective forces acting on small scales are just as crucial to acknowledge as those acting on larger scales. As a whole, our findings imply that future population genomic studies should expand from focusing on model organisms and/or studying heterogeneous regions to better understand the evolutionary processes shaping current and future biodiversity patterns, particularly when used in a comparative phylogeographic context.

Ocean warming and acidification have complex interactive effects on the dynamics of a marine fungal disease.

Diseases threaten the structure and function of marine ecosystems and are contributing to the global decline of coral reefs. We currently lack an understanding of how climate change stressors, such as ocean acidification (OA) and warming, may simultaneously affect coral reef disease dynamics, particularly diseases threatening key reef-building organisms, for example crustose coralline algae (CCA). Here, we use coralline fungal disease (CFD), a previously described CCA disease from the Pacific, to examine these simultaneous effects using both field observations and experimental manipulations. We identify the associated fungus as belonging to the subphylum Ustilaginomycetes and show linear lesion expansion rates on individual hosts can reach 6.5 mm per day. Further, we demonstrate for the first time, to our knowledge, that ocean-warming events could increase the frequency of CFD outbreaks on coral reefs, but that OA-induced lowering of pH may ameliorate outbreaks by slowing lesion expansion rates on individual hosts. Lowered pH may still reduce overall host survivorship, however, by reducing calcification and facilitating fungal bio-erosion. Such complex, interactive effects between simultaneous extrinsic environmental stressors on disease dynamics are important to consider if we are to accurately predict the response of coral reef communities to future climate change.

Coral micro-fragmentation assays for optimizing active reef restoration efforts.

The global decline of coral reefs has driven considerable interest in active coral restoration. Despite their importance and dominance on mature reefs, relatively few coral restoration projects use slower growth forms like massive and encrusting coral species. Micro-fragmentation can increase coral cover by orders of magnitude faster than natural growth, which now allows cultivation of slow growing massive forms and shows promise and flexibility for active reef restoration. However, the major causes of variation in growth and survival of outplanted colonies remain poorly understood. Here, we report simple outplanting assays to aid in active reef restoration of slower growing species and increase the likelihood of restoration success. We used two different micro-fragmentation assays. Pyramid assays were used to examine variation associated with fragment size (ranging from ≈1-9 cm2), nursery residence time (for both in-situ and ex-situ nurseries), and 2D vs. 3D measurements of growth. Block assays were used to examine spatial variation among individual performance at outplanting sites in the field. We found 2D and 3D measurements correlated well, so measured survivorship and growth using top-down planar images for two of the main Hawaiian reef building corals, the plating Montipora capitata and the massive Porites compressa. Pyramid assays housed and outplanted from the in-situ nursery showed no effect of residence time or size on overall survivorship or growth for either species. Results from the ex-situ nursery, however, varied by species, with P. compressa again showing no effect of nursery residence time or size on survivorship or growth. In contrast, nursery culture resulted in improved survivorship of small M. capitata fragments, but net growth showed a weak positive effect of nursery time for medium fragments. Small fragments still suffered higher mortality than either medium or large fragments. Due to their lower mortality, medium fragments (≈3 cm2) appear to be the best compromise between growth and survivorship for outplanting. Likewise, given weak positive gains relative to the investment, our results suggest that it could be more cost-effective to simply outplant medium fragments as soon as possible, without intermediate culture in a nursery. Furthermore, the block assay revealed significant differences in survivorship and growth among sites for individuals of both species, emphasizing the importance of considering spatial variation in coral survival and growth following outplanting. These results highlight the value of using short-term micro-fragmentation assays prior to outplanting to assess size, and location specific performance, optimizing the efficiency of active reef restoration activities and maximizing the probability of success for active coral restoration projects.

Gross and microscopic morphology of lesions in Cnidaria from Palmyra Atoll, Central Pacific.

We conducted gross and microscopic characterizations of lesions in Cnidaria from Palmyra Atoll, Central Pacific. We found growth anomalies (GA) to be the most commonly encountered lesion. Cases of discoloration and tissue loss were rare. GAs had a focal or multi-focal distribution and were predominantly nodular, exophytic, and umbonate. In scleractinians, the majority of GAs manifested as hyperplasia of the basal body wall (52% of cases), with an associated absence or reduction of polyp structure (mesenteries and filaments, actinopharynx and tentacles), and depletion of zooxanthellae in the gastrodermis of the upper body wall. In the soft corals Sinularia sp. and Lobophytum sp., GAs exclusively manifested as prominent hyperplasia of the coenenchyme with an increased density of solenia. In contrast to scleractinians, soft coral GAs displayed an inflammatory and necrotizing component with marked edema of the mesoglea, accompanied by infiltrates of variably-sized granular amoebocytes. Fungi, algae, sponges, and Crustacea were present in some scleractinian GAs, but absent in soft coral GAs. Fragmentation of tissues was a common finding in Acropora acuminata and Montipora cf. dilatata colonies with tissue loss, although no obvious causative agents were seen. Discoloration in the zoanthid, Palythoa tuberculosa, was found to be the result of necrosis, while in Lobophytum sp. discoloration was the result of zooxanthellar depletion (bleaching). Soft corals with discoloration or tissue loss showed a marked inflammatory response, however no obvious causative organisms were seen. Lesions that appeared similar at the gross level were revealed to be distinct by microscopy, emphasizing the importance of histopathology.

Spatial and temporal patterns of scleractinian coral, soft coral, and zoanthid disease on a remote, near-pristine coral reef (Palmyra Atoll, central Pacific).

There is an urgent need for accurate baselines of coral disease prevalence across our oceans in order for sudden or unnatural changes to be recognized. Palmyra Atoll allows us to study disease dynamics under near-pristine, functionally intact conditions. We examined disease prevalence among all known species of scleractinian coral, soft coral and zoanthid (Palythoa) at a variety of coral reef habitats at Palmyra over a 2 yr period. In 2008, overall disease prevalence across the atoll was low (0.33%), but higher on the shallower backreef (0.88%) and reef terrace (0.80%) than on the deeper forereef (0.09%). Scleractinian coral disease prevalence was higher (0.30%) than were soft coral and zoanthid disease (0.03% combined). Growth anomalies (GAs) were the most commonly encountered lesions, with scleractinian species in the genera Astreopora (2.12%), Acropora (1.30%), and Montipora (0.98%) showing the highest prevalence atoll-wide. Discoloration necrosis (DN) was most prevalent in the zoanthid Palythoa tuberculosa (1.18%), although the soft coral Sinulana and Montipora also had a prevalence of 0.44 and 0.01%, respectively. Overall disease prevalence within permanently marked transects increased from 0.65% in 2008 to 0.79% in 2009. Palythoa DN contributed most to this increased prevalence, which coincided with rising temperatures during the 2009 El Niño. GAs on the majority of susceptible genera at Palmyra increased in number over time, and led to tissue death. Host distribution and environmental factors (e.g., temperature) appear to be important for determining spatiotemporal patterns of disease at Palmyra. More sophisticated analyses are required to tease apart the likely inter-correlated proximate drivers of disease occurrence on remote, near-pristine reefs.

Genomic divergence and differential gene expression between crustacean ecotypes across a marine thermal gradient.

Environmental gradients between marine biogeographical provinces separate distinct faunal communities. However, the absence of absolute dispersal barriers allows numerous species to occur on both sides of such boundaries. While the regional populations of such widespread species are often morphologically indistinguishable from each other, genetic evidence suggests that they represent unique ecotypes, and likely even cryptic species, that may be uniquely adapted to their local environment. Here, we explored genomic divergence in four sympatric southern African decapod crustaceans whose ranges span the boundary between the cool-temperate west coast (south-eastern Atlantic) and the warm-temperate south coast (south-western Indian Ocean) near the southern tip of the African continent. Using genome-wide data, we found that all four species comprise distinct west coast and south coast ecotypes, with molecular dating suggesting divergence during the Pleistocene. Transcriptomic data from the hepatopancreas of twelve specimens of one of these species, the mudprawn Upogebia africana, which were exposed to either 10 °C or 20 °C, showed a clear difference in gene expression profiles between the west- and south coast ecotypes. This difference was particularly clear at 10 °C, where individuals from the south coast experienced a 'transcriptomic shock'. This low temperature is more typical of the west coast during upwelling events, and the physiological stress experienced by the south coast ecotype under such conditions may explain its absence from that coastline. Our results shed new light on the processes involved in driving genomic divergence and incipient speciation along coastlines with porous dispersal barriers.

Evolutionary genomics of endangered Hawaiian tree snails (Achatinellidae: Achatinellinae) for conservation of adaptive capacity.

Phylogenomic studies can provide insights into speciation, adaptation, and extinction, while providing a roadmap for conservation. Hawaiian tree snails are a model system for an adaptive radiation facing an extinction crisis. In the last 5 years, nearly all populations of Hawaiian tree snails across the 30 remaining species in the subfamily Achatinellinae (Achatinellidae) have declined from hundreds or thousands in the wild down to undetectable levels. Nearly 100 species historically occurred across dramatic environmental gradients on five of the Hawaiian Islands, but habitat loss, overcollection, and predation by invasive species have decimated populations. As such, this system offers the opportunity to integrate efforts to conserve evolutionary potential into conservation planning for a rapidly declining subfamily. Here, we used genome-wide, restriction-site associated DNA sequencing (RADseq), along with mitochondrial genome reconstruction, to resolve evolutionary relationships to inform conservation efforts. Phylogenetic analysis of nearly 400k genome-wide SNPs from 59 populations and 25 species across six genera in the family Achatinellidae, was generally concordant with taxonomy, geography, and mtDNA with several notable exceptions; mtDNA was unable to resolve some deeper nodes (e.g., the monophyly of Achatinella), while SNP data did not resolve as many shallow nodes. Both phylogenetic and coalescent analysis revealed deep divergences between populations within Achatinella mustelina that were consistent with species-level differences. Given cryptic species-level divergence within populations that are geographically proximate, they are at higher risk of extirpation from invasive predators and climate change than previously assumed. This study clarifies evolutionary relationships within this model system for adaptive radiation, forming the basis for conservation strategies such as translocation, captive rearing, and hybridization trials to prevent the loss of capacity to adapt to rapidly changing environmental conditions.

Unexpectedly high genetic diversity in a rare and endangered seabird in the Hawaiian Archipelago.

Seabirds in the order of Procellariiformes have one of the highest proportions of threatened species of any avian order. Species undergoing recovery may be predicted to have a genetic signature of a bottleneck, low genetic diversity, or higher rates of inbreeding. The Hawaiian Band-rumped Storm Petrel ('Ake'ake; Hydrobates castro), a long-lived philopatric seabird, suffered massive population declines resulting in its listing under the Endangered Species Act in 2016 as federally Endangered. We used high-throughput sequencing to assess patterns of genetic diversity and potential for inbreeding in remaining populations in the Hawaiian Islands. We compared a total of 24 individuals, including both historical and modern samples, collected from breeding colonies or downed individuals found on the islands of Kaua'i, O'ahu, Maui, and the Big Island of Hawai'i. Genetic analyses revealed little differentiation between breeding colonies on Kaua'i and the Big Island colonies. Although small sample sizes limit inferences regarding other island colonies, downed individuals from O'ahu and Maui did not assign to known breeding colonies, suggesting the existence of an additional distinct breeding population. The maintenance of genetic diversity in future generations is an important consideration for conservation management. This study provides a baseline of population structure for the remaining nesting colonies that could inform potential translocations of the Endangered H. castro.

Atlantia, a new genus of Dendrophylliidae (Cnidaria, Anthozoa, Scleractinia) from the eastern Atlantic.

Atlantia is described as a new genus pertaining to the family Dendrophylliidae (Anthozoa, Scleractinia) based on specimens from Cape Verde, eastern Atlantic. This taxon was first recognized as Enallopsammia micranthus and later described as a new species, Tubastraea caboverdiana, which then changed the status of the genus Tubastraea as native to the Atlantic Ocean. Here, based on morphological and molecular analyses, we compare fresh material of T. caboverdiana to other dendrophylliid genera and describe it as a new genus named Atlantia in order to better accommodate this species. Evolutionary reconstruction based on two mitochondrial and one nuclear marker for 67 dendrophylliids and one poritid species recovered A. caboverdiana as an isolated clade not related to Tubastraea and more closely related to Dendrophyllia cornigera and Leptopsammia pruvoti. Atlantia differs from Tubastraea by having a phaceloid to dendroid growth form with new corallites budding at an acute angle from the theca of a parent corallite. The genus also has normally arranged septa (not Portualès Plan), poorly developed columella, and a shallow-water distribution all supporting the classification as a new genus. Our results corroborate the monophyly of the genus Tubastraea and reiterate the Atlantic non-indigenous status for the genus. In the light of the results presented herein, we recommend an extensive review of shallow-water dendrophylliids from the Eastern Atlantic.

A genomic glance through the fog of plasticity and diversification in Pocillopora.

Scleractinian corals of the genus Pocillopora (Lamarck, 1816) are notoriously difficult to identify morphologically with considerable debate on the degree to which phenotypic plasticity, introgressive hybridization and incomplete lineage sorting obscure well-defined taxonomic lineages. Here, we used RAD-seq to resolve the phylogenetic relationships among seven species of Pocillopora represented by 15 coral holobiont metagenomic libraries. We found strong concordance between the coral holobiont datasets, reads that mapped to the Pocillopora damicornis (Linnaeus, 1758) transcriptome, nearly complete mitochondrial genomes, 430 unlinked high-quality SNPs shared across all Pocillopora taxa, and a conspecificity matrix of the holobiont dataset. These datasets also show strong concordance with previously published clustering of the mitochondrial clades based on the mtDNA open reading frame (ORF). We resolve seven clear monophyletic groups, with no evidence for introgressive hybridization among any but the most recently derived sister species. In contrast, ribosomal and histone datasets, which are most commonly used in coral phylogenies to date, were less informative and contradictory to these other datasets. These data indicate that extant Pocillopora species diversified from a common ancestral lineage within the last ~3 million years. Key to this evolutionary success story may be the high phenotypic plasticity exhibited by Pocillopora species.

The complete mitochondrial genome of Achatinella sowerbyana (Gastropoda: Pulmonata: Stylommatophora: Achatinellidae).

In this study, we report the complete mitochondrial genome sequence of Achatinella sowerbyana, an endangered Hawaiian tree snail. The mitogenome is 15,374 bp in length and has a base composition of A (36.4%), T (42.7%), C (9.1%), and G (11.8%). Similar to other pulmonates, it contains 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. The gene order is the same as its sister species, A. mustelina. To our knowledge, this is the second mitochondrial genome sequenced within the superfamily Achatinelloidea, and will aid in the examination of the genus Achatinella, part of a spectacular radiation in the Hawaiian Islands.

Cryptic species obscure introduction pathway of the blue Caribbean sponge (Haliclona (Soestella) caerulea), (order: Haplosclerida) to Palmyra Atoll, Central Pacific.

Cryptic species are widespread across the phylum Porifera making the identification of non-indigenous species difficult, an issue not easily resolved by the use of morphological characteristics. The widespread order Haplosclerida is a prime example due to limited and plastic morphological features. Here, we study the reported introduction of Haliclona (Soestella) caerulea from the Caribbean to Palmyra Atoll via Hawai'i using morphological characteristics and genetic analyses based on one nuclear (18s rDNA) and three mitochondrial (COI, the barcoding COI extension (COI ext.) and rnl rDNA) markers. Despite no clear division in lengths of the oxea spicules between the samples, both mtDNA and nDNA phylogenetic trees supported similar topologies resolving two distinct clades. Across the two clades, the concatenated mtDNA tree resolved twelve subclades, with the COI ext. yielding most of the variability between the samples. Low sequence divergence values (0.68%) between two of the subclades indicate that the same species is likely to occur at Palmyra, Hawai'i and the Caribbean, supporting the hypothesis that H. caerulea was introduced to Palmyra from the Caribbean, although whether species came directly from the Caribbean to Palmyra or from Hawai'i remains unresolved. Conversely, the pattern of highly divergent cryptic species supports the notion that traditionally used spicule measurements are taxonomically unreliable in this group. This study illustrates how understanding the scale of within- as opposed to between-species level genetic variation is critical for interpreting biogeographic patterns and inferring the origins of introduced organisms.

Restriction of sponges to an atoll lagoon as a result of reduced environmental quality.

The lagoon at Palmyra Atoll in the central Pacific was subject to major military modifications during WWII and now the dominant fauna on the lagoon's hard substrate are sponges, not corals. In this study, we quantified the physical and biological factors explaining the variation in sponge distribution patterns across 11 sites to determine the potential for the sponges in the lagoon at Palmyra to invade the surrounding reef systems. Significant differences in sponge assemblages were found among all but three sites. For all the models we examined the strongest environmental relationships were found for variables related to sedimentation/turbidity and food/habitat availability. Our findings suggest that the sponges in Palmyra's lagoon are likely to be restricted to this habitat type where they are associated with conditions resulting from the earlier heavy disturbance and are unlikely to spread to the outer reef environments unless there is a dramatic decline in environmental quality.

Proximate environmental drivers of coral communities at Palmyra Atoll: establishing baselines prior to removing a WWII military causeway.

A management proposal aims to partly remove a WWII military causeway at Palmyra Atoll to improve lagoon water circulation and alleviate sedimentation stress on the southeast backreef, an area of high coral cover and diversity. This action could result in a shift in sedimentation across reef sites. To provide management advice, we quantified the proximate environmental factors driving scleractinian coral cover and community patterns at Palmyra. The proportion of fine sedimentation was the optimal predictor of coral cover and changes in community structure, explaining 23.7% and 24.7% of the variation between sites, respectively. Scleractinian coral cover was negatively correlated with increases in fine sedimentation. Removing the causeway could negatively affect the Montipora corals that dominate the western reef terrace, as this genus was negatively correlated with levels of fine sedimentation. The tolerance limits of corals, and sediment re-distribution patterns, should be determined prior to complete removal of the causeway.

Modeling patterns of coral bleaching at a remote Central Pacific atoll.

A mild bleaching event (9.2% prevalence) at Palmyra Atoll occurred in response to the 2009 ENSO, when mean water temperature reached 29.8-30.1 degrees C. Prevalence among both abundant and sparse taxa varied with no clear pattern in susceptibility relating to coral morphology. Seven taxon-specific models showed that turbidity exacerbated while prior exposure to higher background temperatures alleviated bleaching, with these predictors explaining an average 16.3% and 11.5% variation in prevalence patterns, respectively. Positive associations occurred between bleaching prevalence and both immediate temperature during the bleaching event (average 8.4% variation explained) and increased sand cover (average 3.7%). Despite these associations, mean unexplained variation in prevalence equalled 59%. Lower bleaching prevalence in areas experiencing higher background temperatures suggests acclimation to temperature stress among several coral genera, while WWII modifications may still be impacting the reefs via shoreline sediment re-distribution and increased turbidity, exacerbating coral bleaching susceptibility during periods of high temperature stress.