Abstract
We analyzed the spatial distribution patterns of a white syndrome (WS) outbreak affecting Pseudodiploria strigosa colonies in the northern Mexican Caribbean during 2018–2019. The purpose of the study was to describe the outbreak progression in a single species and determine if this WS incidence is related to the nearest diseased neighbor distance. Two separated sites with different P. strigosa colonial densities (Bocana: 0.08 col/m2; Picudas: 0.2 col/m2) were selected in similar habitats of the same reef complex. P. strigosa colonies within the survey sites were mapped, and their status was recorded (healthy, diseased, or dead) in sequential surveys until colonies died or the study terminated (306 days). Spatial distribution modes were assessed using Ripley’s K function. The spatial colony distribution was random in one site (Bocana) and clustered in the other (Picudas). However, the WS disease incidence per survey was randomly distributed in both sites throughout the observation period of the outbreak, suggesting that WS transmission at small spatial scales was independent of the colony distribution pattern and from the nearest diseased colonies. Survival probability since WS onset in surveyed colonies was different: 0% at Bocana and 14% at Picudas by April 2019. But, eventually, all diseased colonies died in both sites. WS outbreak timing was different at the two sites: Initial prevalence 8% at the Bocana site vs. 44% at Picudas site. Distribution of time to disease onset shown multimodality, with modes varying from 17 to 184 days and wide main modes amplitude suggest a highly variable resistance to the WS. Disease incidence was not abated during winter surveys. Differences between sites in the WS disease outbreak distribution and progression suggest that colony condition, environmental quality, and perhaps several transmission events played an essential role in the complex outbreak dynamics at the local spatial scale of our study.
Introduction
In recent years outbreaks of a lethal white syndrome (WS) affecting many species in several Caribbean reefs have been reported. The nature, timing, and regional spread of the outbreak in the Caribbean region (), strongly suggest a similar outbreak in the NE Mexican Caribbean () as well as that previously reported for Florida (). In recent publications, these outbreaks have been named “stony coral tissue loss disease” (SCTLD; ), although in our surveys signs also correspond to the white plague type II (WP-II) description (). The spreading spatial patterns have suggested a primary infectious causal agent and diminishing tissue loss rates in affected colonies treated with antibiotics reinforce the idea that a bacterial infection is involved in SCTDL syndrome development (). Nevertheless, so far, although several putative bacterial agents have been proposed (; ; ; ; Thomé et al., 2021) no causality has been demonstrated, only association, which is not evidenced enough (; ). Moreover, given the above, it seems that opportunistic pathogens may play a confusing role to the extent that consider that not just bacteria, but viruses or even protozoa may play a causal role in the SCTLD. In contrast, , ponder the possibility that there may be no infectious agent in the primary etiology of the SCTLD.
Disease transmission of infectious diseases may occur by direct contact between diseased and susceptible organisms or by indirect transmission by vector-borne (mechanical or biological) and in the case of the marine environment by water-borne transmission (). In one host one parasite system, the diseased triad (host, pathogen, and environment) can be equally applied to disease transmission, instead of the classic one pathogen–one host point of view (). And a more complex panorama arises for disease transmission in highly diverse ecological systems, more so if subject to continuous environmental change (Webster et al., 2017; ). It is considered that a primary infection outbreak tends to show a spatial dependence in diseased corals (; Thrusfield, 2016). However, argue that a clustered distribution of diseased colonies follows a random one and is determined by direct colony contact, a secondary transmission process. Zvuloni et al. (2009) found a clustered distribution of Black Band disease affected colonies at very small spatial scales in dense multispecies coral assemblages. In contrast, report no clustered distribution for SCTLD in multispecies assemblages at relatively small scales in Florida, whereas, report significant clustering at different large spatial scales for SCTLD, also along the Florida reef tract. Webster et al. (2017) point out that in complex ecosystems many simultaneous transmission modes may coexist, further confounding spatial patterns interpretation.
However, our present insight into outbreak transmission patterns in the reef environment is based on multispecies coral assemblages that may not adequately describe the single species response. To improve our understanding of the spatial spread of this WS outbreak at the species level, we studied the temporal incidence variation along with the outbreak in two populations of the highly susceptible scleractinian coral Pseudodiploria strigosa (; ; ) in northern Mexican Caribbean reefs. To characterize the incidence patterns we analyze the WS incidence spatial distribution patterns at several intervals during the outbreak to find out to what extent colony proximity was a determining factor in the spatial distribution incidence patterns of this WS in P. strigosa. We also contrast survival probabilities and the temporal patterns of apparently healthy colonies to become diseased among the relatively nearby two study sites, to further describe the dynamics of the outbreak.
We decided to maintain the hypernym of WS for this outbreak because: (a) the uncertainty of identifying a coral disease simply by visual signs in the field (Work and Aeby, 2006; Thrusfield, 2016); (b) that many “old” Caribbean coral diseases, such as the plagues as originally described are still extant; and (c) that in our surveys we observed both SCTLD and WP-II specific disease signs (SCTLD: tissue sloughing and multifocal lesions. WP-II: basal initiation and unifocal rapid in-colony spread).
Materials and Methods
Study Site
The study was conducted in two back reef zone sites of Puerto Morelos Reef in the northeastern Yucatan peninsula, Mexico: Picudas (20°53′02′′N and 86°50′54′′W) and Bocana (20°52′29′′N and 86°51′06′′W (Figure 1). At these sites P. strigosa colonies were abundant and the WS outbreak was beginning. Both sites were deployed in back reef coral grounds behind the reef crest, over a flat calcareous pavement 4–5 m deep and nearby reef inlets, a distance of 1,200 m separates the two sites. Picudas site has a denser and more diverse reef biota assemblage and the local wave regime is more pronounced due to a lower reef crest and nearby a wide reef opening; the site is rarely visited. Bocana site is at the lagoon-ward end of an active channel through which reef lagoon waters are flushed out to the sea, and is well protected by a shallower reef crest. It is a favorite site both by snorkelers and scuba divers.
FIGURE 1
Colony Mapping
Colonies at both sites were mapped following standard topographic procedures, but underwater. The spatial position of each P. strigosa colony was obtained by estimating the distance and bearing of each colony from one or more fixed central points. Bearings to each colony were measured with an alidade table (0.25 degrees precision), oriented to the magnetic north with a compass, and anchored to prevent its movement. A stretched measuring tape running from the alidade center to the same colony provides distances. Conversion to Cartesian coordinates was carried out by trigonometry, and also plotted in an X/Y plane to allow spatial visualization. The extension of the study area at each site was determined by underwater horizontal visibility when mapping, axes in the order of 15 m, and thus surveyed areas were similar: Bocana 453 m2 (Figure 2A), and Picudas 435 m2 (Figure 3A).
FIGURE 2
FIGURE 3
The monitoring of P. strigosa colonies began in August 2018. Successive surveys were carried out with an irregular frequency (biweekly to monthly), until April 2019 in Bocana and June 2019 in Picudas, comprising the duration of the outbreak at each site. On each survey, the condition of each colony was recorded as healthy (apparently healthy), diseased (showing WS signs), or dead. All colonies were photographed at each survey to record the signs present and the relative progression of the disease.
Spatial Distribution Analysis
The spatial distribution patterns of P. strigosa colonies at each site, and later the spatial incidence distribution (newly diseased colonies) per survey were calculated using Ripley’s K function (). Ripley’s K function computes the expected number of colonies within an increasing radius from an arbitrary colony (; ). It is a cumulative function quantifying the type of spatial pattern (clustered, disperse, and random) along the spatial dimensions of the surveyed area. K(r) was calculated as:
where A is the study area, n is the number of colonies in this area, dij is the distance between two colonies i and j. The variable indicator Ir(dij) takes the value of 1 if the distance among points i and j is less than r and 0 otherwise; Wij is introduced to avoid bias due to border effects, as it weights the points in function to their distance to the border. For easier visualization, we normalize K(r) with Besag’s L(r) transformation (). We use the theoretical L-function for the homogeneous Poisson point process (complete spatial randomness, CSR) as a benchmark that separates spatial clustering from spatial regularity (); If a value of = 0, the observed pattern is random, positive values indicates clustering and negative values indicates overdispersion.
To find whether the spatial patterns of colonies (Figures 2B, 3B) were significantly different from a random distribution we applied a CSR test based on Monte Carlo procedures (). Then was compared with a 95% confidence global envelope of null distribution and if L(r) falls into the envelope the spatial pattern is completely random with 95% of confidence (p < 0.05).
In spatial epidemiology, the distribution pattern of a disease is known as a “marked” pattern because it is the result of two processes: one is determined by the spatial distribution of colonies of the study population and the second is determined by the spatial distribution of the status labeled colonies (diseased, healthy, or dead; ; ). To test for colony spatial pattern dependence among newly diseased colonies we applied a random labeling (RL) test, to each survey on both sites. The null hypothesis is that the label of each colony is randomly allocated, regardless of other points. H0 was evaluated through the D(r) statistic which calculates the difference of the newly diseased colonies per site per survey (Ldis(r)) and the surveyed population at each site (Lpop(r)): D(r)=Ldis(r)-Lpop(r). If H0 is true, then D(r) = 0. Significant deviations were evaluated through a 95% confidence global envelope for null distribution constructed from random permutations of the labels of the complete data. If D(r) falls into the envelope the null is accepted, implying that newly diseased colonies are randomly distributed among the P. strigosa studied population at each site. The hypothesis is rejected when a value falls out of the confidence envelope. Positive significant values (above envelope) suggest that new colonies with WS signs are more likely to be found close to diseased colonies than would be expected if randomly allocated: a clustered distribution of newly diseased colonies. Negative values (under envelope) indicate that cases are more likely to be found close to a colony without WS signs than would expect under RL (; ).
Survival and Time to Disease Onset Analysis
Colony condition status was obtained at intervals, not daily, thus it is known that disease or death onset for a given colony occurred within a given interval, but the date is unknown. Colony time to WS signs appearance or death once diseased (time to event) was calculated for each survey, as the number of days that each colony remains in a given condition (healthy, diseased, or dead). The sum of the colony-days per interval without change in colony condition, plus half the preceding interval length when condition change, provided the total colony-days for a colony in a given condition. This procedure assigns a midpoint to the number of colony-days uncertainty within that interval (actuarial method in survival analyses; ). This type of data is common in epidemiological studies and several statistical techniques have been developed to analyze it (). To estimate the survival probabilities, we apply a non-parametric maximum likelihood estimator (NPMLE) for the distribution of the censored data, to generalize the Kaplan–Meir survival estimates for interval-censored data (). Then the time to an event (disease, death) addresses the estimated probability at each survival interval, which is an average probability that takes into account censored data.
The temporal distribution of colonies to develop WS signs at each site was analyzed utilizing empirical cumulative distributions and kernel density estimators (KDE). Comparisons among sites were performed with Kolmogorov–Smirnov bootstrapped tests designed for discrete data with ties, given that the original KS test requires continuous data (). Observing that the distribution of time to disease onset at each site was not uniform or unimodal, a multimodal analysis was carried out (R multimode package, ). Potential modes were identified with KDE curves generated with the standard nrd0 and plug-in bandwidths () for each site. The magnitude of each local maxima displayed by the KDE plot was considered to decide the number of apparently relevant modes. This number was the required input to obtain the critical bandwidth value of the KD’s as to identify the number and location of numerically valid modes.
All analyses were carried using R (). Survival probabilities were estimated with interval package (); CSR and RL tests, and L(r) plots were made with the packages spatstat () and GET (; ). For uniformity tests (uniftest package, ), bootstrapped KS comparisons (matching package, ) and for modal analysis (multimode package, ).
Results
Outbreak Development
While both reef study sites had a similar survey area (≈440 m2), P. strigosa colony density at Picudas site was more than twice that of Bocana (0.20 cols/m2 vs. 0.08 cols/m2). Notwithstanding, at the beginning of the study WS prevalence at Picudas site was 8.5%, much lower than at Bocana: 44%. At the first survey some colonies were already dead probably attributable to the WS outbreak as recolonization of the denuded skeleton was at its initial stages and there were no other evident causes of death. Coincidental with the initial prevalence levels registered at each site, recent colony mortality was lower in Picudas than in Bocana (three and seven dead colonies, respectively).
The survival probability of colonies since the onset of WS signs was different among sites. For the Bocana site average colony survival probability since WS signs onset was 20% at the K–M estimated interval starting around day 158 (by end of January 2019; Figure 4) and close to 0% around day 214 (April 2019). At Picudas site average colony survival probability since WS signs onset was 14% at the K–M estimated survival interval starting around day 172 (by mid January, 2019) and 13% at the interval starting around day 225 (beginning of April 25). As shown by the proportion of healthy colonies surviving till the end of the study in Figure 5. While the changes in the slope of both K–M curves are relatively different, there was no significant difference in the Kaplan–Meir survival curves among sites (asymptotic log-rank two-sample test: Z = 0.7506, p-value = 0.46), most likely because both slopes approach zero survival. As modeled, all diseased colonies eventually died. Furthermore, no recovery of a diseased colony was registered at either site.
FIGURE 4
FIGURE 5
The time to disease onset data distribution at each site showed no continuous or regular pattern. Assuming the possibility of constant development of disease signs amongst survey intervals a Neyman–Barton uniform distribution test () was applied. Both data sets showed no uniformity (p < 0.00), indicate non-constant signs of development at either site. A Kolmogorov–Smirnov bootstrapped test to allow for discrete data and presence of ties, suggests that empirical cumulative distribution indeed of both sites are marginally different (p = 0.048) in terms of absolute maximum distance among the corresponding functions. However, analyzing the site date using KDEs, it was apparent that time to disease onset patterns was different among sites (Figure 6).
FIGURE 6
While Picudas KDE strongly suggests multimodality, Bocana site data suggest a relatively week bimodality. Testing for unimodality in the underlying density by bootstrapping, it was found that both sites had multimodal distributions (Bocana: p = 0.014; Picudas: p = 0.002). After selecting the appropriate critical bandwidth (Bocana: 19.5757; Picudas: 16.2998), the number and location of local maxima were found to be two for Bocana (modes at 17 and 151 days from study initiation), and three for Picudas (modes at 70, 129, and 184 days). The amplitude of the main modes is quite large suggesting that colony resistance is rather variable at both sites. Aside from the variable colony resistance, other striking findings are the lag in disease outbreak among sites and the discontinuous time to disease onset per survey (as opposed to constant or progressively decreasing trends), as shown by multimodality in the distribution of time to disease data. More so, as different multimodal patterns exist between the sites.
Spatial and Temporal Patterns of the WS Outbreak
No colonies were in direct contact at the Bocana site, but at the Picudas site, there were seven pairs. In six of these pairs, the two colonies appeared healthy at the first survey. The healthy-healthy initial pairs eventually developed signs in one of the colonies and later in the other, except for one couple remaining healthy until the end of the study. In no instance, both colonies in a pair developed WS signs within the same survey (Table 1). The mean time for the apparently healthy colony of a pair to become diseased was 44 days. This period was significantly shorter than that among independent colonies (mean = 83.4 days; Welch t test = t: 3.59, p = 0.006).
TABLE 1
| Pair | Size | Tissue in contact | Time to onset (days) |
| 1 | L–M | AH–AH | 0 |
| Dis–AH | 67 | ||
| Dis–Dis | 31.5 | ||
| 2 | S–S | AH–AH | 0 |
| Dead–Dis | 24 | ||
| 3 | L–M | AH–AH | 0 |
| Dead–Dis | 38 | ||
| 4 | L–S | AH–AH | 0 |
| Dis–AH | 7 | ||
| Dead–AH | 21 | ||
| Dead–Dis | 74.5 | ||
| 5 | M–L | AH–AH | 0 |
| Dead–AH | >107 | ||
| 6 | S–M | AH–AH | 0 |
| Dead–AH | 33 | ||
| Dead–Dis | 35 | ||
| 7 | M–L | Dead–Dis | 0 |
Development of colony condition through time for the contact colony pairs at Picudas site.
Tissue in contact: Apparently Healthy (AH), Diseased (Dis), and Dead. Colony size: Large (L), Medium (M), and Small (S).
The spatial distribution pattern of P. strigosa colonies surveyed at the Bocana site was random all along the distance scale of the Bessag’s–Ripley’s K diagram (1–6 m; CSR test p = 0.53, Figure 2B), whereas at the Picudas site the colony distribution were clustered all along the distance scale (1–5 m; CSR test p = 0.001, Figure 3B). The incidence distribution patterns per survey were adjusted because although in both study sites newly diseased colonies were registered at every survey, in several instances less than 15 newly diseased colonies were registered per survey. Performing Ripley’s K analysis with less than 15 data, result in biased results due to insufficient sample size to estimate the variance adequately (, ). To obtain a sufficient sample size, survey data of successive samplings were pooled until 15 data points were obtained. Pooling data sets are commonly applied to K function analyses (). For instance, combined spatial data, we pooled time data. For the Bocana site, two sets of pooled newly diseased colonies were constructed, for the periods of August–September 2018 and October 2018–April 2019. These periods coincide with the hottest and coldest seasons. At Picudas site, we obtained six adequate (≥15 data points each) periods for the spatial analysis: August 22–October 15, October 26, November 14, and November 27 of 2018, January 29, and March 18–June 24 of 2019. With these data sets, we obtained L(r) functions of newly diseased colonies (Supplementary Figures 1, 2 for Bocana and Picudas sites, respectively), and compute the difference between incidence and population patterns at each site (Ldis(r) – Lpop(r)), to test the RL hypothesis (Figures 7A,B, 8A–F for Bocana and Picudas sites, respectively). The tests show that the WS incidence per period analyzed was randomly distributed except in two cases: the second period at the Bocana site (October 2018–April 2019) indicating an overdispersed pattern, marginally significant at larger distance scales (4.5–6 m, RL test p = 0.05, Figure 7B), perhaps due to a space increase between diseased colonies and scarce healthy colonies. The second exception is for the summer survey at Picudas site (August–September 2018) where incidence distribution shows a marginally significant clustering at 2.5–4.5 m (RL test p = 0.05, Figure 8A) when the abundance of apparently healthy colonies was still high.
FIGURE 7
FIGURE 8
Discussion
A remarkable finding of this study was the different dynamics in the WS outbreak in P. strigosa surveyed populations of the study sites, notwithstanding the two sites are relatively close to each other and in similar environmental settings. Although site surveys started within a week of each other the initially recorded WS prevalence at the Bocana site (44%) was five times greater than at Picudas sites (8.5%), even though P. strigosa colony density at Picudas was more than twice the one at Bocana. This first finding shows that the WS outbreak started at Picudas site, about the same time that our survey started, but earlier at the Bocana site and/or that P. strigosa colonies at Bocana site were more susceptible to the WS. If the outbreak should have occurred simultaneously at both sites it would have resulted in relatively similar survival probabilities, if not lower at Picudas, due to its higher colony density; unless P. strigosa colonies resistance to the WS in Picudas was higher than at Bocana.
A possible outbreak starting time for the Bocana site can be roughly estimated by finding the median time to disease onset, since our observations started, and subtracting that period from our initial survey dates. A median time of apparently healthy colonies to become diseased of 74 days was obtained (pooling data from both sites), with a bootstrapped 95% confidence interval of 60–90 days. The 2–3 months range comprised by the confidence interval roughly corresponds with the probable time frame of the outbreak in the northern Mexican Caribbean as reported by . On the other hand, the median time to disease onset in Bocana colonies since observations started was of 97 days, suggesting that the outbreak developed earlier at that site.
The WS outbreak progression in terms of time to disease onset in apparently healthy colonies was not constant or progressively decreasing at neither site. Trends that may be expected by continuous successful transmission and decreasing host abundance, respectively, if nearest neighbor distances were important. Instead, a multimodal distribution of time to disease onset data was found for both sites, implying a wavy frequency in the WS sign’s appearance. Furthermore, the modal location for each site was different for Bocana (modes at 17 and 151 days from study initiation), and three for Picudas (modes at 70, 129, and 184 days). While at Bocana both modes had a similar percentage of the data (Figure 6), in Picudas the first mode at 70 days has the most of data, and the other two were quite smaller and with similar values. The wavy frequency in the time to disease onset data, suggests different events of successful transmission whereas the large amplitude of the main modes is probably related to large variability in colony resistance to the WS. In this late case, amplitudes of main modes at both sites are similarly wide and in that sense at both sites colony resistance is variable. As the outbreak in our study sites progressed, many colonies died, and live colony density was greatly reduced, thus increasing average nearest neighbor distances, and still, the WS incidence continued unabated, even in colder winter months.
Incidence spatial patterns of the WS outbreak, throughout the study, don’t reflect the surveyed population spatial patterns. Also, a transition of incidence spatial patterns was observed, indicating that the dynamics of this WS outbreak are complex even at small spatial scales. The initial Picudas incidence patterns showed clustering suggesting a nearest neighbor relation in the transmission of the disease; and if so, probably driven by the relatively high colony density at the site. Later surveys, however, showed a transition to random distribution patterns, although many apparently healthy colonies remained in the surveyed area. On the other hand, at the Bocana site the incidence distribution pattern was random in the summer period, but in the winter period transition toward an overdispersed pattern; probably an effect of a reduced population of susceptible colonies in the study area. reported that Aspergillosis in the sea fan Gorgonia ventalina had different distribution patterns as the outbreak developed. At low prevalence levels the distribution pattern of diseased colonies was random, but a transition to clustered distribution at high prevalence levels; due to direct contact after the initial long-distance water-borne driven outbreak. An apparent lack of colony density dependence for the spread of coral diseases has also been documented in Florida reefs for outbreaks of WP-II (), of a WS (; ), and SCTDL (). In our study, lack of spatial dependence seems to extend to very small spatial scales as no simultaneous (not observed in the same survey) disease WS signs development were observed in colonies that grew in direct contact. Nevertheless, although few colony pairs were registered, the shorter period of the apparently healthy colony to become diseased in diseased-healthy pairs suggests that direct contact enhances colony susceptibility to this WS. More so, as there were no cases where a WS lesion appeared in the contact area of a colony pair.
On the other hand, given the density independence of the incidence events, multimodal distributions may be related to periodic transport events of the causal agent or necessary additional factors to trigger the WS disease () at the spatial scale of the study areas. It may also imply different transmission events and mechanisms for this WS outbreak (). infer that the WS they observed in Florida was water-borne transmitted due to the low density of coral colonies. shared this point of view after assessing the rates of spread of SCTLD in the Florida reef tract. , modeling the SCTLD in the Florida reef tract, consider that the disease is likely to be transported by currents as a neutrally buoyant material. Zvuloni et al. (2009) claimed that even at very short-inter colony distances and high colony densities waterborne transmission of a Black Band disease was the main mechanism for a clustered distribution of diseased colonies, in the absence of biological vectors. If there is a primary agent of the disease and is waterborne transported (; ; ) given the short distance between our sites (1,200 m) and high water connectivity in the study area, probably the WS outbreak should have started more or less simultaneously at both study sites.
reported that white plagues are transmitted effectively by Coralliophilia abbreviata, a corallivorous snail, and also through water even when filtered. Williams and Miller (2005) reported that white band disease in Acropora cervicornis was transmitted both by direct colony contact and also by C. abbreviate. In the present study, we observed the fireworm Hermodice carunculata, feeding on tissue remnants of the WS lesion of a P. strigosa colony at Picudas site in one of the surveys; but not on the adjacent, apparently healthy coral tissue (Figure 9). There were no sightings of fireworms or other coral predators such as C. abbreviata feeding on healthy P. strigosa colonies. report active biological vectoring for SCTLD by butterflyfishes in Florida, the first biological vector report to our knowledge. Passive biological vectors had also been suggested, as the macroalgae Halimeda opuntia () or sediment biofilm (). These reports point out the possibility of several co-occurring transmission mechanisms (), during the disease outbreak in our study sites. And the observed particular outcomes may have depended on minor variations in local conditions and colony stress history, as well on the epigenetics of both hosts and potential primary and secondary pathogens present through the outbreak (Webster et al., 2017; Zaneveld et al., 2017).
FIGURE 9
The importance of changing ecosystems effects in coral disease etiology has been stressed by many authors, from the ecological (i.e., ; ; ) to the microbiological context (i.e., ; Zaneveld et al., 2016; Webster et al., 2017; ; Vega-Thurber et al., 2020). Models were the role of opportunistic pathogens, and stress-driven reduced host resistance should perhaps be the norm, rather than the exception (). Such interactions are likely to result in greater susceptibility to disease as a result of unbalanced coral microbiomes (Zaneveld et al., 2017), leading both to host depressed immune response and reduce overall performance (). In addition, environmental degradation can lead to the emergence of new strains of presently purported pathogens or even new ones (; Vega-Thurber et al., 2009). For instance, all colonies in our study sites developed the WS but at different times and with varying survival periods, strongly suggest that individual colony condition played a role in susceptibility to the disease and its progression within a colony. But also opens the possibility that different diseases are co-occurring in the WS outbreak. If consideration that there may be no primary infectious etiology in the SCTLD holds, then multiple secondary bacterial pathogens associated with the SCTLD (; ; ; ; ) and for a more general WS as well (Thomé et al., 2021), may play an important role in these diseases development. As it is well known that external stressors may compromise host immunity (; ) and/or produce an imbalance in the structure and diversity of the coral microbiome, contributing to the proliferation of opportunistic pathogens (; ; Woodley et al., 2016; Sweet and Bulling, 2017; ).
Two main drivers of coral diseases so far considered are increased seawater temperatures and poor water quality, both individually and synergistically (Wang et al., 2018). Coral diseases seem to bear a positive relationship with high seawater temperatures, but not with low seawater temperatures (; ; Zvuloni et al., 2009; ; ; ). However, at our sites, despite the decrease of apparently healthy colonies, the proportion of newly diseased colonies per interval between summer and winter surveys showed no significant difference. Thomé et al. (2021) observed that WS lesion progression was not hinder during winter at Picudas site and observed newly WS diseased colonies in relatively nearby Cozumel Island in winter. If this WS causality is water-borne, the winter strong north winds may contribute to this density-independent pattern by transporting fomites () resulting in waves of effective local transmissions; despite the lower seawater temperatures.
Global and regional stressful factors further confound local factors’ relevance, for instance, , considered that a previous thermal bleaching event was of importance in triggering a WS Florida outbreak, a posteriori identify as SCTLD (). At our sites, no previous strong bleaching occurred but a potential stressful factor affecting these reefs could be the Sargassum spp. brown tides (Van Tussenbroek et al., 2017). The decomposition of thousands of tons of these algae on the shore deteriorates water quality, increases nutrient input reduces transparency, pH, and oxygen levels, and increases ammonia, phosphorus, and concentration of toxic metals (Van Tussenbroek et al., 2017; ). These brown tides were associated with the mortality of seagrass beds and fauna (Van Tussenbroek et al., 2017; ). However, the potential effects of brown tides and drifting Sargasso mats, are coastal broad, and is difficult to imagine a differential effect at the small spatial scale of our study sites. The same holds regarding the relatively low coastal waters quality, pre-Sargasso events, in the Puerto Morelos coastal waters due to explosive coastal development and inadequate water treatment practices (; ; Van Tussenbroek et al., 2017; ).
The difference in direct anthropogenic effects between both sites is on the level of recreational activities. Bocana site receives many snorkelers and scuba divers, whereas the Picudas site is rarely visited. The risks of coral reef tourism-related activities for coral biota are well known (; Zaneveld et al., 2016; Spalding et al., 2017) but are difficult to demonstrate. For instance, some level of stressful effects from snorkeling and scuba diving on the Bocana coral assemblages may be suspected, but directly linking it to the particular WS outbreak dynamics by association is no demonstration.
We have described in detail the dynamics of the outbreak events at our study sites and shown differences in the dynamics of this WS outbreak at the study sites including a lag in disease outbreak, discontinuous events of the time to disease onset, lack of colony density dependence for the spread of the WS, unhampered incidence in winter months and extremely high morbidity. All indications of a complex and powerful phenomenon, even at the small spatial scales of this study. However, explaining the observed patterns is pending. In this context, further investigations with not-so-aggressive diseases may provide a better natural model to understand the processes behind such complex patterns. Meanwhile, we must recognize that these recent outbreaks are but a continuation of the lethal coral epizootics that had decimated Caribbean reefs, during the last 40 years. And doing so, as has been proposed in many other works, preserving a good environmental quality at the local reef scale may effectively contribute to reducing the susceptibility of corals to WS, but also any other extant or future diseases.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
EG-U performed the samplings, analyzed the data, prepared the figures and tables, authored, reviewed the drafts of the manuscript, and approved the final draft. EJ-D conceived, designed and performed the samplings, analyzed the data, authored, or reviewed drafts of the manuscript, and approved the final draft. Both authors contributed to the article and approved the submitted version.
Funding
This work was internally supported by the ICML-UNAM. There was no additional external funding received for this study. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Acknowledgments
We are grateful to the ICMyL and UNAM for providing financial and logistic support to EJ-D; the CONACyT for a scholarship for EG-U; J. Rivera Ortega and R. Rodríguez Martínez for their support in monitoring and samples collection; the Meteorological and Oceanographic Monitoring Academic Service of Puerto Morelos for SST data; and F. Negrete for boat and diving support. We specially thank the two reviewers whose comments greatly improved the quality and scope of this work. Part of the content of this work appears in EG-U’s master’s thesis ().
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2021.669171/full#supplementary-material
References
1
AebyG. S.UshijimaB.CampbellJ. E.JonesS.WilliamsG. J.MeyerJ. L.et al (2019). Pathogenesis of a tissue loss disease affecting multiple species of corals along the florida reef tract.Front. Mar. Sci.6:678. 10.3389/fmars.2019.00678
2
AGRRA (2019). Coral Disease Outbreak. Stony Coral Tissue Loss Disease. Available online at: http://www.agrra.org/coral-disease-outbreak/(accessed January 2020).
3
Álvarez-FilipL.Estrada-SaldívarN.Pérez-CervantesE.Molina-HernándezA.Gonzalez-BarriosF. J. (2019). A rapid spread of the stony coral tissue loss disease outbreak in the Mexican Caribbean.PeerJ7:e8069. 10.7717/peerj.8069
4
Ameijeiras-AlonsoJ.CrujeirasR. M.Rodríguez-CasalA. (2021). multimode: an R package for mode assessment.J. Stat. Softw.971–32. 10.18637/jss.v097.i09
5
AntonovicsJ. (2017). Transmission dynamics: critical questions and challenges.Phil. Trans. R. Soc. B372:20160087. 10.1098/rstb.2016.0087
6
BaddeleyA.RubakE.TurnerR. (2016). Spatial Point Patterns Methodology and Applications with R.Boca Raton, FL: CRC Press.
7
BaddeleyA.TurnerR.RubakE. (2018). Spatial Point Pattern Analysis, Model-Fitting, Simulation, Test: Package ‘spatstat’. Available online at: http://www.spatstat.org(accessed October 2018).
8
BeedenR.MaynardJ. A.MarshalP. A.HeronS. F.WillisB. L. (2012). A framework for responding to coral disease outbreaks that facilitates adaptive management.Environ. Manag.491–13. 10.1007/s00267-011-9770-9
9
BenichouJ.PaltaM. (2005). “Rates, risks, measures of association and impact,” in Handbook of Epidemiology, edsAhrensW. W.PigeotI. (Berlin: Spinger-Verlag), 89–156. 10.1007/3-540-26577-5_3
10
BesagJ. (1977). Discussion of Dr Ripley’s paper.J. R. Stat. Soc. Ser. B39193–195.
11
BivandR. S.PebesmaE.Gómez-RubioV. (2013). Applied Spatial Data Analysis with R, 2nd Edn. New York, NY: Springer.
12
BrooksD. R.BoegerW. A. (2019). Climate change and emerging infectious diseases: evolutionary complexity in action.Curr. Opin. Syst. Biol.1375–81. 10.1016/j.coisb.2018.11.001
13
BrunoJ. F.CôtéI. M.TothL. T. (2019). Climate change, coral loss, and the curious case of the parrotfish paradigm: why don’t marine protected areas improve reef resilience?Annu. Rev. Mar. Sci.11307–334. 10.1146/annurev-marine-010318-095300
14
BrunoJ. F.SeligE. R. (2007). Regional decline of coral cover in the Indo-Pacific: timing, extent, and subregional comparisons.PLoS One2:e711. 10.1371/journal.pone.0000711
15
CarruthersT. J. B.van TussenbroekB. I.DennisonW. C. (2005). Influence of submarine springs and wastewater on nutrient dynamics of Caribbean seagrass meadows.Estuar. Coast. Shelf Sci.64191–199. 10.1016/j.ecss.2005.01.015
16
ClemensE.BrandtM. E. (2015). Multiple mechanisms of transmission of the Caribbean coral disease white plague.Coral Reefs341179–1188. 10.1007/s00338-015-1327-6
17
DiggleP. J.ChetwyndA. G. (1991). Second-order analysis of spatial clustering for inhomogeneous populations.Biometrics471155–1163. 10.2307/2532668
18
DobbelaereT.MullerE. M.GramerL. J.HolsteinD. M.HanertE. (2020). Coupled epidemio-hydrodynamic modeling to understand the spread of a deadly coral disease in Florida.Front. Mar. Sci.7:591881. 10.3389/fmars.2020.591881
19
EllnerS. P.JonesL. E.MydlarzL. D.HarvellC. D. (2007). Within-host disease ecology in the sea fan gorgonia ventalina: modeling the spatial immunodynamics of a coral-pathogen interaction.Am. Nat.170143–161.
20
EvansS. A. (1976). Causation and disease: the henle-koch postulates revisited.Yale J. Biol. Med.49175–195.
21
FayM. P. (2020). Package ‘interval’. Weighted Logrank Tests and NPMLE for Interval Censored Data. Available online at: https://cran.r-project.org/web/packages/interval/interval.pdf(accessed April 2021).
22
FayM. P.ShawP. A. (2010). Exact and asymptotic weighted logrank tests for interval censored data: the interval package.J. Stat. Softw.36:i02.
23
FKNMS (Florida Keys National Marine Sanctuary) (2018). Florida Reef Track Coral Disease Outbreak. Case Definition: Stony Coral Tissue Loss Disease (SCTLD). Available online at: https://nmsfloridakeys.blob.core.windows.net/floridakeys-prod/media/docs/20181002-stony-coral-tissue-loss-disease-case-definition.pdf(accessed October 2, 2018).
24
GlaslB.HerndlG. J.FradeP. R. (2016). The microbiome of coral surface mucus has a key role in mediating holobiont health and survival upon disturbance.ISME J.102280–2292. 10.1038/ismej.2016.9
25
Guzmán-UrietaE. O. (2020). Patrones Espaciales de un Síndrome Blanco en el Coral Pseudodiploria Strigosa del Arrecife de Puerto Morelos, Quintana Roo. Ph. D. thesis. Mexico City: Universidad Nacional Autónoma de México.
26
HarvellC. D.Jordan-DahlgrenE.MerkelS.RosenbergE.RaymundoL.GarrietS.et al (2007). Coral disease, environmental drivers, and the balance between coral and microbial associates.Oceanography20172–195. 10.5670/oceanog.2007.91
27
Hernández-TerronesL.Rebolledo-VieyraM.Merino-IbarraM.SotoM.Le-CossecA.Monroy-RíosE. (2011). Groundwater pollution in a karstic region (NE Yucatan): baseline nutrient content and flux to coastal ecosystems.Water Air Soil Pollut.218517–528. 10.1007/s11270-010-0664-x
28
HeronS. F.WillisB. L.SkirvingW. J.EakinM. C.PageC. A.MillerI. R. (2010). Summer hot snaps and winter conditions: modeling white syndrome outbreaks on great barrier reef corals.PLoS One5:e12210. 10.1371/journal.pone.0012210
29
HillA. B. (1965). The environment and disease: association or causation?Proc. R. Soc. Med.58295–300. 10.1177/003591576505800503
30
Hoegh-GuldbergO. (1999). Climate change, coral bleaching and the future of the world’s coral reefs.Mar. Freshwater Res.50839–866.
31
IwanowiczD. D.SchillW. B.WoodleyC. M.BrucknerA.NeelyK.BriggsK. M. (2020). Exploring the stony coral tissue loss disease bacterial pathobiome.Preprint bioRxiv [Preprint]. 10.1101/2020.05.27.120469
32
JollesA. E.SulivanP.AlkerA. P.HarvellD. (2002). Disease transmission of aspergillosis in sea fans: inferring process from spatial pattern.Ecology832373–2378. 10.1890/0012-9658(2002)083[2373:dtoais]2.0.co;2
33
JonesR. J.BowyerJ.Hoegh-GuldbergO.BlackallL. (2004). Dynamics of a temperature-related coral disease outbreak.Mar. Ecol. Prog. Ser.28163–77. 10.3354/meps281063
34
LambJ. B.TrueD. J.PiromvaragornS.WillisB. L. (2014). Scuba diving damage and intensity of tourist activities increases coral disease prevalence.Biol. Conserv.17888–96. 10.1016/j.biocon.2014.06.027
35
LandsbergJ. H.KiryuY.PetersE. C.WilsonP. W.PerryN.WatersY.et al (2020). Stony coral tissue loss disease in florida is associated with disruption of host–zooxanthellae physiology.Front. Mar. Sci.7:576013. 10.3389/fmars.2020.576013
36
LesserM. P.BythellJ. C.GatesR. D.JohnstoneR. W.Hoegh-GuldbergO. (2007). Are infectious diseases really killing corals? Alternative interpretations of the experimental and ecological data.J. Exp. Mar. Biol. Ecol.34636–44. 10.1016/j.jembe.2007.02.015
37
MelnikM.PusevR. (2015). uniftest: Tests for Uniformity. R package version 1.1.
38
MeyerJ. L.Castellanos-GellJ.AebyG. S.HäseC. C.UshijimaB.PaulV. J. (2019). Microbial community shifts associated with the ongoing stony coral tissue loss disease outbreak on the florida reef track.Front. Microbiol.10:2244. 10.3389/fmicb.2019.02244
39
MullerE. M.SartorC.AlcarazN. I.van WoesikR. (2020). Spatial epidemiology of the stony-coral-tissue-loss disease in Florida.Front. Mar. Sci.7:163. 10.3389/fmars.2020.00163
40
MullerE. M.van WoesikR. (2012). Caribbean coral diseases: primary transmission or secondary infection?Glob. Chang. Biol.183529–3535. 10.1111/gcb.12019
41
MyllymäkiM.MrkvickaT. (2020). GET: global envelopes in R.Preprint arXiv [Preprint]. arXiv:1911.06583 [stat.ME],
42
MyllymäkiM.MrkvickaT.GrabarnikP.SeijoH.HahnU. (2017). Global envelope test for spatial processes.J. R. Stat. Soc. B79381–404. 10.1111/rssb.12172
43
NeelyK. L.MacaulayK. A.HowerE. K.DoblerM. A. (2020). Effectiveness of topical antibiotics in treating corals affected by stony coral tissue loss disease.PeerJ8:e9289. 10.7717/peerj.9289
44
NoonanK. R.ChildressM. J. (2020). Association of butterflyfishes and stony coral tissue loss disease in the Florida Keys.Coral Reefs391581–1590. 10.1007/s00338-020-01986-8
45
NuguesM.SmithG.van HooidonkR.SeabraM. I.BakR. P. M. (2004). Algal contact as a trigger for coral disease.Ecol. Lett.7919–923. 10.1111/j.1461-0248.2004.00651.x
46
PaulV. J.UshijimaB.AebyG. (2019). Studies of the Ecology and Microbiology of Florida’s Coral Tissue Loss Diseases.Miami, FL: Florida DEP, 1–33.
47
PrechtW. F.GintertB. E.RobbartM. L.FuraR.van WoesikR. (2016). Unprecedented disease-related coral mortality in southeastern Florida.Sci. Rep.6:31374. 10.1038/srep31374
48
R Core Team (2019). R: A Language and Environment for Statistical Computing.Vienna: R Foundation for Statistical Computing.
49
RichardsonL. (1996). Motility patterns of Phormidium corallyticum and Beggiatoa spp. associated with black band disease of corals.Microb. Ecol.32323–335.
50
RichardsonL. L.GoldbergW. M.KutaK. G.AronsonR. B.SmithG. W.RitchieK. B. (1998). Florida’s mystery coral-killer identified.Nature392557–558. 10.1038/33302
51
RipleyB. (1977). Modelling spatial patterns.J. R. Stat. Soc. B30172–192.
52
RippeJ. P.KriefallN. G.DaviesS. W.CastilloK. D. (2019). Differential disease incidence and mortality of inner and outer reef corals of the upper Florida Keys in association with a white syndrome outbreak.Bull. Mar. Sci.95305–316. 10.5343/bms.2018.0034
53
Rodríguez-MartínezR. E.Medina-ValmasedaA. E.BlanchonP.Monroy-VelázquezL. V.Almazán-BecerrilA.Delgado-PechB.et al (2019). Faunal mortality associated with massive beaching and decomposition of pelagic Sargassum.Mar. Pollut. Bull.146201–205. 10.1016/j.marpolbul.2019.06.015
54
Rodríguez-MartínezR. E.RoyP. D.Torrescano-ValleN.Cabanillas-TeránN.Carrillo-DomínguezS.Collado-VidesL.et al (2020). Element concentrations in pelagic Sargassum along the Mexican Caribbean coast in 2018-2019.PeerJ8:e8667. 10.7717/peerj.8667
55
Ruiz-MorenoD.WillisB. L.PageA. C.WeilE.CróquerA.Vargas-ÁngelB.et al (2012). Global coral disease prevalence associated with sea temperature anomalies and local factors.Dis. Aquat. Org.100249–261. 10.3354/dao02488
56
SekhonJ. S. (2011). Multivariate and propensity score matching software with automated balance optimization: the matching package for R.J. Stat. Softw.421–52.
57
SharpW. C.SheaC. P.MaxwellK. E.MullerE. M.HuntJ. H. (2020). Evaluating the small-scale epidemiology of the stony-coral-tissue-loss-disease in the middle Florida Keys.PLoS One15:e0241871. 10.1371/journal.pone.0241871
58
SheatherS. J.JonesM. C. (1991). A reliable data-based bandwidth selection method for kernel density estimation.J. R. Stat. Soc. B.53683–690.
59
ShoreA.CaldwellJ. M. (2019). Modes of coral disease transmission: how do diseases spread between individuals and among populations?Mar. Biol.166:45. 10.1007/s00227-019-3490-8
60
SokolowS. (2009). Effects of a changing climate on the dynamics of coral infectious disease: a review of the evidence.Dis. Aquat. Organ.875–18. 10.3354/dao02099
61
SpaldingM. D.BurkeL.WoodS. A.AshpoleJ. (2017). Mapping the global value and distribution of coral reef tourism.Mar. Policy82104–113. 10.1016/j.marpol.2017.05.014
62
SweetM. J.BullingM. T. (2017). On the importance of the microbiome and pathobiome in coral health and disease.Front. Mar. Sci.4:9. 10.3389/fmars.2017.00009
63
ThoméP. E.Rivera-OrtegaJ.Rodríguez-VillalobosJ. C.Cerqueda-GarcíaD.Guzmán-UrietaE. O.García-MaldonadoJ. Q.et al (2021). Local dynamics of a white syndrome outbreak and changes in the microbial community structure associated with colonies of the scleractinian brain coral Pseudodiploria strigosa.PeerJ9:e10695. 10.7717/peerj.10695
64
ThrusfieldM. (2016). “Epidemiology,” in Diseases of Coral, edsWoodleyC. M.DownsC. A.BrucknerA. W.PorterJ. W.GallowayS. B. (Hoboken, NJ: John Wiley & Sons, Inc), 28–51.
65
Van TussenbroekB. I.Hernández-AranaH. A.Rodríguez-MartínezR. E.Espinoza-AvalosJ.Canizales-FloresH. M.González-GodoyC. E.et al (2017). Severe impacts of brown tides caused by Sargassum spp. on nearshore Caribbean seagrass communities.Mar. Pollut. Bull.122272–281. 10.1016/j.marpolbul.2017.06.057
66
Vega-ThurberR.MydlarzL. D.BrandtM.HarvellD.WeilE.RaymundoL.et al (2020). Deciphering coral disease dynamics: integrating host, microbiome, and the changing environment.Front. Ecol. Evol.8:575927. 10.3389/fevo.2020.575927
67
Vega-ThurberR.Willner-HallD.Rodríguez-MuellerB.DesnuesC.EdwardsR. A.AnglyF.et al (2009). Metagenomic analysis of stressed coral holobionts.Environ. Microbiol.112148–2163. 10.1111/j.1462-2920.2009.01935.x
68
WangL.ShantzA. A.PayetJ. P.SharptonT. J.FosterA.BurkepileD. E.et al (2018). Corals and their microbiomes are differentially affected by exposure to elevated nutrients and a natural thermal anomaly.Front. Mar. Sci.5:101. 10.3389/fmars.2018.00101
69
WebsterJ. P.BorlaseA.RudgeJ. W. (2017). Who acquires infection from whom and how? Disentangling multi-host and multi-mode transmission dynamics in the ‘elimination’ era.Phil. Trans. R. Soc. B372:20160091. 10.1098/rstb.2016.0091
70
WilliamsD. E.MillerM. W. (2005). Coral disease outbreak: pattern, prevalence and transmission in Acropora cervicornis.Mar. Ecol. Prog. Ser.301119–128. 10.3354/meps301119
71
WoodleyC. M.DownsC. A.BrucknerA. W.PorterJ. W.GallowayS. B. (2016). Diseases of Coral.Hoboken, NJ: John Wiley & Sons Inc.
72
WorkT. M.AebyG. S. (2006). Systematically describing gross lesions in corals.Dis. Aquat. Org.70155–160. 10.3354/dao070155
73
ZaneveldJ.McMindsR.Vega ThurberR. (2017). Stress and stability: applying the Anna Karenina principle to animal microbiomes.Nat. Microbiol.2:17121. 10.1038/nmicrobiol.2017.121
74
ZaneveldJ. R.BurkepileD. E.ShantzA. A.PritchardC. E.McMindsR.PayetJ. P.et al (2016). Overfishing and nutrient pollution interact with temperature to disrupt coral reefs down to microbial scales.Nat. Commun.7:11833. 10.1038/ncomms11833
75
ZvuloniA.Artzy-RandrupY.StoneL.Kramarsky-WinterE.BarkanR.LoyaY. (2009). Spatio-temporal transmission patterns of black-band disease in a coral community.PLoS One4:e4993. 10.1371/journal.pone.0004993
Summary
Keywords
SCTLD, white plague, coral disease, incidence, transmission
Citation
Guzmán-Urieta EO and Jordán-Dahlgren E (2021) Spatial Patterns of a Lethal White Syndrome Outbreak in Pseudodiploria strigosa. Front. Mar. Sci. 8:669171. doi: 10.3389/fmars.2021.669171
Received
18 February 2021
Accepted
03 August 2021
Published
03 September 2021
Volume
8 - 2021
Edited by
William F. Precht, Dial Cordy and Associates, Inc., United States
Reviewed by
Greta Smith Aeby, Qatar University, Qatar; Brian Keith Walker, Nova Southeastern University, United States
Updates
Copyright
© 2021 Guzmán-Urieta and Jordán-Dahlgren.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Edgar O. Guzmán-Urieta, edomguzman@outlook.com
This article was submitted to Coral Reef Research, a section of the journal Frontiers in Marine Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.