Introduction

Coral reefs are diverse ecosystems that directly support numerous human populations providing valuable sources of dietary protein and income from fisheries and tourism (Hicks et al. 2019). Unfortunately, the world’s coral reefs are in severe decline (Pandolfi et al. 2003, 2011) with an estimated 50% decline in coral cover from the Indo-Pacific (Bruno and Selig 2007; Death et al. 2012) and a 40% loss from Caribbean reefs (Gardner et al. 2003; Jackson et al. 2014). Coral reefs are frequently situated along coastlines and are subject to increasing interconnections of human populations and their associated impacts (Williams et al. 2019). Overfishing which removes critical herbivores and land-based pollution with resultant sedimentation, excess nutrients and other pollutants are cited as major local stressors contributing heavily to coral decline (Pandolfi 2003; Hughes et al. 2010). Compounding these threats are increases in sea surface temperatures, changing weather patterns and ocean acidification associated with climate change (Pandolfi et al. 2011; Hughes et al. 2017). Coral disease has been linked to various stressors, including temperature stress, anthropogenic exploitation and pollution (Harvell et al. 2007), and disease outbreaks have occurred in numerous regions (Aeby 2005; Williams and Miller 2005; Williams et al. 2011a; Aeby et al. 2016; Precht et al. 2016). Despite this importance, we are just beginning to understand the ecology and pathogenesis of coral diseases. Numerous questions remain about the causative agents of coral diseases, their prevalence, epidemiology and ecological impacts. While coral disease prevalence has been quantified in many regions across the Indo-Pacific (Willis et al. 2004; Aeby et al. 2008, 2011a, 2015; Haapkylä et al. 2009; Vargas-Ángel 2009; Williams et al. 2011b), we know little about the levels and types of coral disease on reefs in many other regions, including around the Arabian Peninsula (Riegl et al. 2012).

The Persian/Arabian Gulf (hereafter called the Persian Gulf) and Arabian Sea host fringing and patch reefs with over 25 coral genera documented (Coles 2003). These coral species, typical of Indo-Pacific fauna, represent a hardy subset of corals that are able to exist under demanding environmental conditions such as extreme temperature and salinity fluctuations (Riegl et al. 2012). These coral reefs are now under constant threat from multiple human impacts including coastal construction, dredging activity and effluent discharge from seawater desalinization plants (Sheppard et al. 2010; Burt et al. 2014) and, recently, disease (Coles 1994; Riegl 2002; Riegl et al. 2012).

The first step in understanding threats from disease is to develop baseline information on the types, distribution and prevalence of coral diseases on reefs. Here, we conducted baseline coral disease surveys, along patch and fringing reefs, of the northeastern Arabian Peninsula including within the Persian Gulf, Strait of Hormuz and Oman Sea. These three regions represent a gradient of environmental conditions presenting a natural experiment to examine the effect of stressors on coral health. The Persian Gulf is a semi-enclosed marginal sea, separated from the Indian Ocean by the narrow Strait of Hormuz (Sheppard 1993). Its shallow depth and restricted water exchange with the wider Indian Ocean result in large seasonal fluctuations in ocean temperatures, 15 °C–36 °C (John et al. 1990) and chronic hypersalinity (> 40 ppt) (Vaughan and Al-Mansoori 2019). These extremes in temperature and salinity create stressful conditions that control the distribution and diversity of coral reef species (Bauman et al. 2013). In contrast, the Strait of Hormuz and Oman Sea have wider connections with the larger Indian Ocean, and so, coral reefs exist in less extreme environments with temperatures ranging from 23 to 30 °C although there are regional exceptions associated with seasonal monsoons and upwelling events (Claereboudt 2019) and salinities averaging 37 ppt (Pous et al. 2004).

The aim of our study was to determine the types and prevalence of coral diseases on coral reefs in the southern Persian Gulf, Strait of Hormuz and Oman Sea. We document taxonomic differences in disease susceptibility and explore environmental drivers of regional variation in disease prevalence. Such information will provide critical data that are essential to help guide reef management efforts in a region where such data are currently deficient. Baseline disease data give managers the capacity to track the health of their reefs through time and can focus monitoring efforts on sites or corals that are more prone to disease.

Materials and methods

Disease surveys

Coral community structure and disease prevalence were recorded along the northeastern Arabian Peninsula at 17 sites having different environmental conditions. Surveys were conducted at eight sites in the Persian Gulf, five sites in the Strait of Hormuz and four sites in Oman Sea, which represent a gradient of salinity and temperature extremes and potential human impacts (Fig. 1; Supplementary Table 1). Additionally, reefs closer to shore, within these regions, may be more strongly impacted by anthropogenic activities such as coastline modification and dredging, as well as discharge from desalinization, wastewater treatment and industrial facilities (Sheppard et al. 2010). At each site, teams of divers surveyed two replicate 25 m belt transects recording the diversity and abundance of coral genera and the prevalence of coral disease. Corals were identified to the genus level along 25 × 1 m belts and substrate characteristics were documented by line-intercept method with the substratum underlying the tape measure recorded at 50-cm intervals.

Fig. 1
figure 1

Map of the Arabian Peninsula showing sites surveyed for coral diseases in the southern Persian Gulf, Strait of Hormuz and the Gulf of Oman. Site GPS coordinates and depths are provided in Supplementary Table S1

Coral lesions were assessed along wider 25 × 2 m belts. Gross lesions were classified into three broad categories including tissue loss, discoloration and growth anomaly, and nomenclature for lesions was based on the host genus affected and the lesion type (e.g., Porites growth anomaly; Work and Aeby 2006). Tissue loss lesions were further classified as known cause (e.g., predation) or unknown cause (e.g., infectious or non-infectious disease processes) based on the lesion size, shape, presence of predators, knowledge of what common predation marks look like and evidence of lesion progression based on degree of algal colonization onto the bare coral skeleton. Only lesions that were thought to be due to disease processes were included in the analyses. Paired normal and lesion tissues were sampled for histopathologic examination and examined for the presence of organisms that could be contributing to gross lesions (Work et al. 2012). Tissue loss lesions not found associated with obvious microorganisms were termed ‘white syndrome’ indicating a tissue loss disease of unknown etiology. If ciliates or cyanobacteria were found microscopically invading coral tissue in the lesion, then they were termed ciliate or cyanobacterial infections.

Geospatial Environmental Predictors

To examine the association between disease prevalence and environmental variables, survey data were combined with spatially explicit environmental data to develop statistical models (Williams et al. 2010; Aeby et al. 2011a). To model different anthropogenic stressors on coral reefs, we used existing 2013 data for three metrics of human impacts on marine environments from the Ocean Health Index (OHI) (Halpern et al. 2015): coastal human population, organic pollution and ocean pollution (Table 1, Supplementary Table 1). The OHI creates a spatially explicit stress index (1 km resolution) for each metric of ocean impact by synthesizing remotely sensed data. The impact of human population was based on satellite estimates of human population size within 10 km of each coast-adjacent ocean grid cell. Organic pollution impact was based on plume models of pesticide use and land use. Ocean pollution impact was derived from pollutant plumes originating in ports and number of shipping tracks from Automatic Identification System and Volunteer Observing System data. For each of these metrics, environmental predictor values at each georeferenced survey location were found using ArcGIS 10.5. Proximity to shore was estimated by measuring the Euclidean distance (m) of each survey site from the coastline using the Near function in ArcGIS. Long-term temperature-based metrics using satellite-derived sea surface temperature (SST) were calculated, and the metric chosen to represent temperature stress was temperature range for each site. This metric was calculated by taking the mean of the annual temperature range each year over a 5-yr time period leading up to the survey date. The annual values from 5 km daily temperature values were calculated from the National Oceanic and Atmospheric Administration (NOAA) Coral Reef Watch ‘CoralTemp’ SST version 3.1 (https://coralreefwatch.noaa.gov/product/5km/index_5km_sst.php). Maximum degree heating week (DHW) values were used and derived from monthly maximum composites of NOAA Coral Reef Watch DHWs, a measure of cumulative heat stress over a 12-week rolling window and measured at 5 km spatial resolution; the monthly composites were from June, July and August 2016 (Table 1).

Table 1 Predictor variables used to model environmental drivers of coral disease along the Arabian Peninsula

Statistical analyses

Quantifying disease prevalence and coral composition (univariate analyses)

Time constraints underwater prevented enumeration of all coral colonies within the wider belt transects surveyed for disease. Hence, the prevalence of lesions was calculated by extrapolating colony counts within the 25 × 1 m transect to the wider 25 × 2 m disease survey area and by using this as the denominator of prevalence calculations (e.g., (number of colonies with lesions/total number of estimated colonies) × 100). The frequency of disease occurrence (FOC) was defined as the number of sites having corals with lesions divided by total number of sites surveyed. FOC reflects the spatial distribution of diseases on reefs. Data were not normally distributed, even with transformation, so nonparametric analyses were used. For each ocean region, a Chi square test for equality of distributions examined differences in the distribution of the number of diseased versus healthy colonies among the scleractinian genera affected by disease. All calculations for disease prevalence or FOC were by coral genera (e.g., prevalence of Porites yellow-banded tissue loss disease was calculated as the total number of infected Porites colonies divided by the estimated number of Porites colonies surveyed, multiplied by 100). A Kruskal–Wallis test examined regional differences in coral cover, number of coral genera and colony densities. Nonparametric statistics were performed using JMP Pro 13 statistical software (SAS Institute Inc., Buckinghamshire, UK).

Examining coral and disease assemblage patterns (multivariate analyses)

Differences in coral communities and disease assemblages among regions were investigated using multivariate community analyses. We tested for differences in coral host community patterns (proportional abundance of coral genera) and coral disease assemblage patterns (disease prevalence) for the fixed factor ‘Region’ (three levels: Persian Gulf, Strait of Hormuz and Oman Sea) using permutational multivariate analysis of variance (PERMANOVA) (Anderson 2001). The raw percentage coral cover values were used to construct a Bray–Curtis similarity matrix, as the values were directly comparable across regions, and models were based on 9999 random permutations of the raw data under a reduced model and type III (partial) sums of squares. For the disease assemblage data, we used a zero-adjusted Bray–Curtis similarity matrix due to the relatively large number of zero prevalence records for some of the diseases within the sites (Clarke et al. 2006). Three sites that contained non-true zero prevalence values (i.e., the particular coral host was not recorded at that site) or denuded (all zero) disease prevalence values were removed from the analysis (Musandam Island, Delma West, Delma north), leaving 14 sites in total. Where the global model was significant, we calculated PERMANOVA pairwise comparisons between factor levels (i.e., regions) and visualized the results using a constrained canonical analysis of principal coordinates (CAP) (Anderson and Willis 2003). To identify the core coral genera and disease types driving differences among factor levels within the CAP plot, we calculated Pearson’s correlation values of the CAP ordination axes with the original coral and disease variables. Coral genera and disease types with strong correlations (defined as ≥ 0.6 in this study) were then overlaid as a bi-plot to visualize the direction and strength of the correlation.

Modeling disease–environment associations (multivariate analyses)

To identify correlations between overall disease prevalence (univariate response) and disease assemblages (multivariate response) with environmental and human drivers that could affect disease processes, we used a permutational distance-based multiple regression technique (DISTLM) (McArdle and Anderson 2001). DISTLM carries out a partitioning of variation in a data set described by a resemblance matrix according to a regression (or multiple regression) model. Predictor variables can be categorical or continuous and the technique makes no prior assumptions about the nature of the response variable distribution, meaning that normality does not have to be satisfied. Predictor variables included coral cover, depth, measures of potential heat stress (degree heating weeks, SST range) and potential human stress (nearby human population size, organic pollution, ocean pollution, proximity to shore) (Table 1). The model was based on 9999 random permutations of the raw data and a Euclidean distance matrix. The predictors were first investigated for collinearity by calculating their pairwise Pearson’s correlation values (r = < 0.6 in all cases), and therefore, all eight predictors were included in the model fitting exercise. During the model fitting process, all possible candidate models (i.e., all possible unique combinations of the predictors) were fitted. Importantly, when more than one predictor was included, its ability to capture variation in the response variable was conditional on the performance of other predictors included in the model. The final optimal model (balancing performance against model complexity) was identified using Akaike’s information criterion (Akaike 1973) with a second-order bias correction applied (AICc) to help account for the relatively low sample sizes (Hurvich and Tsai 1989; Burnham and Anderson 2004). This was chosen as the method to identify the most parsimonious model as it adds a ‘penalty’ for increases in the number of predictor variables. We report the top four models for both response variables (overall disease prevalence and disease assemblage structure) to highlight the change in AICc values with changing model complexity.

Results

Coral reef characteristics and coral community structure

Across the 17 sites we surveyed along the Arabian Peninsula, overall average live coral cover was 51.2% (range = 13–83%) and average macroalgae cover was 0.1% (range = 0–2%). No significant differences among regions were found for colony densities (Kruskal–Wallis, X2 = 0.2, df = 2, p = 0.89) or coral cover (Kruskal–Wallis, X2 = 5.3, df = 2, p = 0.07), but the number of coral genera differed significantly among regions (Kruskal–Wallis, X2 = 8.8, df = 2, p = 0.01) with sites in the Persian Gulf having the lowest number of coral genera compared to sites within the Strait of Hormuz or Oman Sea (Table 2). We identified 25 coral genera within transects (Table 3) and overall (all surveys combined) Porites was the most common coral genus (avg. = 32.8% of the coral community), followed by Acropora (avg. = 14.9%), Faviidae (Favites and Dipsastraea) (avg. = 9.7%), Platygyra (avg. = 9.6%) and Montipora (avg. = 7.6%). Coral community structure was significantly different among regions (PERMAVOVA, pseudo-F2,16 = 2.88, p = 0.006) and was mainly driven by a higher abundance of Cyphastrea at sites within the Persian Gulf compared to a higher abundance of several genera at sites within the Strait of Hormuz (Oxypora, Stylophora, Montipora, Echinopora, Pocillopora, Pavona) and higher abundances of Favites and Leptastrea at sites within the Oman Sea (Fig. 2).

Table 2 Comparison of substrate characteristics among regions surveyed in 2016
Table 3 Coral community composition of three regions surveyed in August, 2016
Fig. 2
figure 2

Canonical analysis of principal coordinates (CAP) ordination quantifying the relative similarity in coral community structure across three regions of the Arabian Peninsula (n = 17 sites). The closer the sites are in multivariate space, the more similar their coral communities are, with those coral genera exerting the strongest influence on spatial variations in community structure overlaid as a bi-plot (the length of the vector is proportional to the strength of the correlation). Squared canonical correlations equal 0.937 and 0.746 for CAP1 and CAP2, respectively

Types of diseases, occurrence and disease assemblage patterns among regions

Across the northeastern Arabian Peninsula, we found 13 coral diseases. These included tissue loss diseases of unknown etiology (white syndromes) in Porites, Platygyra, Dipsastraea, Cyphastrea, Acropora and Goniopora; growth anomalies in Porites, Platygyra and Dipsastraea (formerly Favia); black band disease (BBD) in Platygyra, Dipsastraea, Acropora, Echinopora and Pavona; bleached patches in Porites and Goniopora and yellow-banded tissue loss disease in Porites (Figs. 3, 4). Disease prevalence among sites ranged from 0 to 4.6% with sites within the Persian Gulf having the highest average prevalence compared to the other two regions (Kruskal–Wallis, X2 = 6.1, df = 2, p = 0.048) (Fig. 5). Frequency of occurrence (FOC) differed among diseases and regions (Table 4) with the three most widespread diseases, which were found in all regions, being Platygyra growth anomalies found at 52.9% of the sites (all regions combined), Acropora white syndrome (47.1% of the sites) and Porites bleached patches (35.3% of the sites). Prevalence differed among diseases and regions (Table 5) with the highest average disease prevalence being Echinopora BBD (2.8%) within the Strait of Hormuz, followed by Acropora white syndrome in the Persian Gulf (2.5%). Disease assemblages differed significantly between regions (PERMANOVA, pseudo-F2, 13 = 3.1, P = 0.004), with sites in the southern Persian Gulf having significantly different disease assemblages than sites in the Strait of Hormuz (t = 1.6, P = 0.018) and Oman Sea (t = 2.1, P = 0.004), primarily driven by all types of disease being elevated in the Persian Gulf as compared to the other two regions.

Fig. 3
figure 3

Representative photos of coral diseases found during surveys at sites within the Persian Gulf, Strait of Hormuz and Oman Sea. aPorites yellow-banded tissue loss, bPorites white syndrome, cPorites growth anomaly, dPorites bleached patch, ePlatygyra white syndrome, fPlatygyra growth anomaly

Fig. 4
figure 4

Representative photographs of coral diseases found during surveys at sites within the Persian Gulf, Strait of Hormuz and Oman Sea. a Black band disease on Acropora, bAcropora white syndrome, cCyphastrea white syndrome, d Faviidae (Dipsastraea) growth anomaly, e Faviidae (Dipsastraea) white syndrome, fGoniopora bleached patch, gGoniopora white syndrome

Fig. 5
figure 5

Differences in average disease prevalence (%) among three regions surveyed along the Arabian Peninsula in August 2016. Data reflect mean and standard error

Table 4 Frequency of occurrence (FOC) of different coral diseases within three ocean regions surveyed along the Arabian Peninsula in August 2016
Table 5 Average prevalence (%) of different coral diseases within three regions surveyed along the Arabian Peninsula in August, 2016

Lacking any historical information on normal disease levels in these regions, we defined a site as having a localized disease outbreak if prevalence was higher than the highest average regional prevalence we found, which was within the Persian Gulf (2.1%). Based upon this criterium, we identified eight of 17 sites that had localized outbreaks of a specific disease with prevalence ranging from 2.3 to 21.1%. The Persian Gulf had the highest number of sites with outbreaks (5 of 8 sites) which included an outbreak of Porites yellow-banded tissue loss (prevalence = 6.9%) at Saadiyat, Platygyra white syndrome (prevalence = 10.2%) at Ras Ghanada and three sites (Coral 85, Sir Bu Nair North, Sir Bu Nair Northwest) with outbreaks of Acropora white syndrome: prevalence = 2.3%, 6.6% and 6.8%. Two sites in the Strait of Hormuz had BBD outbreaks, Eagle Bay had 3.9% BBD prevalence on Echinopora and Abu Rasheed Island had 7.1% BBD prevalence on Echinopora and 3.3% BBD prevalence on Pavona. In the Oman Sea, one site (Dibba) had an outbreak of black band disease (BBD) on Pavona (21.1%).

Differences in disease prevalence among coral genera and regions

Twenty-five coral genera were found within survey transects, but only seven showed signs of disease. Disease prevalence differed among the coral genera within each ocean region (Persian Gulf X2 = 261.4, df = 6, p < 0.01; Strait of Hormuz X2 = 58.4, df = 8, p < 0.01; and Oman Sea X2 = 276, df = 7, p < 0.01). Coral genera that were affected more commonly than expected based on their abundance in the community included Acropora and Platygyra at sites within the Persian Gulf, Echinopora, Pavona, Platygyra and Porites at sites within the Strait of Hormuz and Pavona at sites within the Oman Sea (Fig. 6). Consistent among all three regions, coral genera that showed less disease than expected (disease resistant) included Cyphastrea and Faviidae (Favites and Dipsastraea) (Fig. 6).

Fig. 6
figure 6

Relationship between overall abundance of coral genera within transects and overall proportion of corals with disease signs (lesions) along the Arabian Peninsula. Top graph shows data (all surveys combined) from the Persian Gulf, middle graph from Strait of Hormuz and bottom graph from Oman Sea. Faviidae includes Favites and Dipsastraea

Disease–environment associations

Of all the possible models containing unique combinations of the eight predictor variables, the optimal model, for predicting spatial variation in disease prevalence, explained 82.2% of the variability across sites and contained three predictors: range in mean annual sea surface temperature (59.4% of the variation), organic pollution (13.6%) and proximity to shore (9.1%) (Table 6). A positive relationship was found between disease prevalence and each of the three predictors (Fig. 7). However, it should be noted that the positive relationship between disease prevalence and organic pollution is skewed by a single site (Saadiyat) so that relationship is less clear. Disease assemblage patterns were best explained by two predictors: range in mean annual sea surface temperature (27.9% of the variability) and proximity to shore (21.1%). These two variables explained 49.0% of the overall variation in disease assemblages across sites (Table 6).

Table 6 Distance-based permutational multiple regression analyses (DISTLM) relating overall disease prevalence and disease assemblages to eight predictor variables (depth, coral cover, organic pollution, ocean pollution, distance to shore, degree heating weeks and average range of SST)
Fig. 7
figure 7

Scattergrams with trend lines showing the relationship between disease prevalence and a mean annual SST range, b modeled organic pollution, and c proximity to shore. These three variables make up the optimal model explaining spatial variation in disease prevalence at survey sites within the Persian Gulf, Strait of Hormuz and Oman Sea

Discussion

We surveyed 17 reefs along the northeastern Arabian coast and found high coral cover (40 to 65%) and diversity (25 coral genera). Thirteen diseases were found affecting seven coral genera present within transects. Diseases included tissue loss diseases of unknown etiology (white syndromes) in six genera, black band disease in five genera, growth anomalies in three genera and yellow-banded tissue loss (Arabic yellow band) in one genus.

Common and unique diseases impact Arabian corals

Tissue loss disease has been reported in a number of coral genera from the Persian Gulf (Riegl 2002; Riegl et al. 2012; Alidoost Salimi et al. 2017) and from Oman Sea (Coles 1994), but our study extends the list of affected genera to include Cyphastrea and Goniopora. Tissue loss is a common manifestation of disease in corals and has been found worldwide either affecting a single coral genus or species at a site (Raymundo et al. 2005; Williams and Miller 2005; Aeby et al. 2011b; Williams et al. 2011b) or affecting multiple coral genera (Willis et al. 2004; Precht et al. 2016). The underlying etiologies of tissue loss diseases differ with multiple bacterial pathogens implicated (Patterson et al. 2002; Sussman et al. 2008; Ushijima et al. 2012, 2014, 2016; Beurmann et al. 2017), chimeric parasites (Work et al. 2011) as well as a variety of other organisms (Work et al. 2012). Tissue loss lesions in a single coral species can also have different etiologies. For example, three different pathogenic bacterial species can cause tissue loss disease in the reef coral, Montipora capitata, in Kaneohe Bay, Oahu. The pathogen, Vibrio owensii, creates chronic tissue loss lesions (Ushijima et al. 2012), and both Vibrio coralliilyticus and Pseudoalteromonas piratica produce acute tissue loss lesions, but under different environmental conditions (Ushijima et al. 2014; Beurmann et al. 2017). As such, we consider tissue loss diseases in this region as non-specific host responses to a variety of agents.

Bleached patches and growth anomalies on corals along the Arabian Peninsula present biogeographic range extensions for both types of lesions. Bleached patches have been reported in multiple coral genera (Lawrence et al. 2014; Aeby et al.2015), and while their etiologies are unknown, viruses have been implicated for Porites on the Great Barrier Reef (Lawrence et al. 2014). Growth anomalies have also been reported in multiple coral genera from other regions (Sutherland et al. 2004) and Acropora from Oman (Coles and Seapy 1998), but this is the first report of growth anomalies in Porites, Platygyra and Faviidae (Favites and Dipsastraea) in this region. Growth anomalies in Porites are widespread across the Indo-Pacific (Aeby et al. 2011c), but there are limited reports of growth anomalies in Faviidae (Sutherland et al. 2004) and only two reports of growth anomalies in Platygyra: one from the Great Barrier Reef (Loya et al. 1984) and one from New Caledonia (Aeby et al. 2015).

The one disease that appears to be unique to the Persian Gulf is Porites yellow-banded tissue loss. This disease presents with subacute progressive tissue loss, but the area between denuded skeleton and healthy tissue is a yellow mucus band as opposed to a bare, white skeleton as found with typical white syndromes. Similar disease signs on Porites and other coral genera have been reported in the Persian Gulf and were termed Arabic yellow-banded disease (Riegl 2002; Riegl et al. 2012) and Porites peeling tissue loss (Alidoost Salimi et al. 2017), which distinguishes them from the standard ‘white syndrome’ classification. Under the dissecting microscope, the yellow band was found to be overlying bare skeleton and was composed of a viscous material with zooxanthellae present, suggesting the yellow band is derived from coral tissue. In the field, we also frequently observed red filamentous cyanobacteria in the region between the yellow band and healthy tissue. It is not known what role, if any, the cyanobacteria play in the disease process, but the lesion appearance and the ecology of the disease presented very differently than black band disease.

Some diseases were widespread such as Platygyra growth anomalies, Porites bleached patches and Acropora white syndrome, which were found in all three regions examined. Other diseases had a smaller geographic spread such as Porites white syndrome found only at sites in the Persian Gulf, Porites yellow-banded tissue loss which was found mainly in the Persian Gulf with one site in the Strait of Hormuz and Porites growth anomalies which occurred exclusively in the Strait of Hormuz and Oman Sea. Although Porites is a common coral in all three regions, the types of lesions that developed on Porites in these regions differ (e.g., growth anomalies in Strait of Hormuz and Oman Sea versus tissue loss lesions in the Persian Gulf). The emergence of coral disease occurs from a complex interplay between the host, causative agent and environment (Work and Aeby 2006), and so, it is not surprising that the environmental conditions occurring within the different regions are facilitating the development of different types of diseases. Environmental conditions are a known component affecting coral health (Harvell et al. 2007) and this has been found in numerous studies. For example, Porites is a common coral genus across the Indo-Pacific, but the prevalence of growth anomalies in this genus is strongly associated with human presence (Aeby et al. 2011c).

Taxonomic variation in disease susceptibility

Corals possess constituent immunity that contribute to their capacity to resist disease, and investment in immunity differs among coral taxa resulting in certain species or families being more susceptible to disease than others (Palmer et al. 2010). Consistent differences in disease susceptibility across families or species have been found in numerous regions (Willis et al. 2004; Aeby et al. 2008, 2011b; Vargas-Ángel 2009), and we also found differences in disease prevalence among coral genera with only seven out of 25 genera within transects showing any signs of disease. Coral genera within our surveys with higher disease prevalence than expected based on their abundance in the field (higher disease susceptibility) included Acropora, Pavona, Platygyra and Echinopora. Acropora is reported as especially vulnerable to disease in many regions of the Pacific (Willis et al. 2004; Aeby et al. 2008; Vargas-Ángel 2009; Myers and Raymundo 2009; Page et al. 2009; Haapkylä et al. 2009; Aeby et al. 2011c; Weil et al. 2012), Indian Ocean (Thinesh et al. 2011; Montano et al. 2012) as well as the Western Atlantic (Santavy et al. 2001; Porter et al. 2001). Reports of localized outbreaks of disease in Acropora have also been reported from across the Pacific (Willis et al. 2004; Aeby 2005; Aeby et al. 2011b; Williams et al. 2011a), and disease has almost extirpated acroporids from Florida and some regions in the Caribbean (Gladfelter 1982; Aronson and Precht 2001; Patterson et al. 2002). Our study confirms that Acropora within the Persian Gulf are also vulnerable to diseases.

Discovering relatively high disease prevalence in Pavona, Platygyra and Echinopora is unusual and is in contrast to the low disease levels usually found in these genera in other regions (Aeby et al. 2008; Myers and Raymundo 2009; Vargas-Ángel 2009; Thinesh et al. 2011; Montano et al. 2012; Aeby et al. 2015). However, in our study Echinopora and Pavona were both affected exclusively by black band disease (BBD). Black band disease affects numerous coral genera worldwide and is more common during the summer months (Sutherland et al. 2004) which is when we surveyed the reefs. So, it may be more accurate to say these two genera were more susceptible to BBD than other coral genera, but not necessarily to disease in general, consistent with other studies. In contrast, Platygyra was affected by BBD, growth anomalies and white syndromes and high disease levels in this coral genus may represent a difference in disease susceptibility in corals along the northeastern Arabian coast as compared to other regions. However, there are many areas across the world that lack baseline disease surveys and that need to be addressed before potential differences in disease susceptibility among coral genera can be understood.

Regional patterns in coral disease prevalence are characterized by a high frequency of localized outbreaks

Our study found an overall prevalence of 1.1% which is within the range reported for other regions. Disease prevalence of less than 1% has been reported (Aeby et al. 2008, 2015; Haapkylä et al. 2009; Vargas-Ángel 2009) as well as over 7% (Willis et al. 2004; Weil et al. 2012). Disease prevalence clearly varies between regions and points to the importance of developing an understanding of what ‘normal’ disease levels are for your local reefs. Having baseline disease data and building disease assessment into monitoring programs will give managers the capacity to recognize and respond quickly to any changes in disease levels. Although we lack historical information on disease levels in this region, compared to baseline field surveys in other regions (Willis et al. 2004; Aeby et al. 2008; Haapkylä et al. 2009; Williams et al. 2011b; Vargas-Ángel 2009), we found a high number of localized outbreaks of specific diseases with almost half of our total survey sites having outbreak levels of disease. Disease is a natural component of all populations, but disease outbreaks indicate a shift in the host–pathogen–environment triad of disease causation. Temperature stress has been associated with a number of disease outbreaks (Miller et al. 2009; Bruno et al. 2007; Brandt and McManus 2009), and our field surveys were conducted in August, corresponding to peak summer temperatures in the region and that might, in part, help explain the high number of outbreaks. Three of our sites had BBD outbreaks and two sites had white syndrome outbreaks, two diseases which are known to occur more commonly in the warmer summer months (Kuta and Richardson 1996; Willis et al. 2004; Bruno et al. 2007; Sato et al. 2009; Williams et al. 2011a). Alternatively, perhaps localized disease outbreaks regularly occur on reefs in these regions due to other environmental stressors. Longer-term disease studies, including surveys during different seasons, are needed to understand disease dynamics on northeastern Arabian coral reefs.

Persian Gulf corals are particularly vulnerable to coral disease impacts

Coral reefs within the Persian Gulf have to contend with large fluctuations in ocean temperatures and salinities, leading to a lower diversity of coral species (John et al. 1990; Vaughan and Al-Mansoori 2019). However, many coral species live under these challenging environmental conditions and their ability to tolerate extreme temperatures has been linked to thermal adaptations by both the corals and their photosynthetic symbiotic algae (Howells et al. 2016a; Smith et al. 2017; Kirk et al. 2018). Although Persian Gulf corals may be able to tolerate high ocean temperatures, our findings of a high number of outbreaks on coral reefs within the southern Persian Gulf suggest that their abilities to fight diseases may be compromised. For some coral species, hosting thermally tolerant symbiont taxa comes at a cost such as reduced growth (Little 2004; Jones and Berkelmans 2010), lipid stores and egg sizes (Jones and Berkelmans 2011) or increased disease susceptibility (Stat et al. 2008; Rouzé et al. 2016). Shore-Maggio et al. (2018) found that for the reef coral, Montipora capitata, in Hawaii, the trade-off for corals was between disease and bleaching susceptibility. Coral colonies hosting the thermally tolerant variants of Durusdinium were more susceptible to tissue loss diseases compared to colonies hosting Cladocopium symbionts. It is possible that similar trade-offs between bleaching resistance and disease susceptibility are occurring in some Persian Gulf coral species and further research in this area would be of interest.

Conversely, perhaps the higher disease levels found in Persian Gulf corals are a result of chronic poor environmental conditions. Chronic temperature and salinity stress both exert significant energetic costs on corals, with reduced growth, reproduction and survivorship the outcome (Coles and Jokiel 1992; Porter et al. 1999). Corals in the southern Persian Gulf are also exposed to significant sedimentation (Riegl et al. 2012) which negatively affects coral health (Fabricius 2005). The energetic costs of chronic multiple environmental stressors do not come without trade-offs. Bauman et al. (2013) found that the size structure of coral populations in the southern Persian Gulf (Dubai and Abu Dhabi) was significantly smaller and their size structured skewed compared to conspecifics in the Strait of Hormuz (Musandam). Additionally, Howells et al. (2016b, 2018) found that Persian Gulf corals tend to have lower rates of calcification and fecundity compared to conspecifics in the Gulf of Oman. Corals also mount an immune response when faced with environmental challenges which can be energetically costly (Palmer 2018). In laboratory experiments, Sheridan et al. (2014) found that acute exposure to sterile sediments induced a significant immune response in the reef coral, M. patula, and resulted in a reduction in storage lipids. In addition to chronic challenges, coral reefs within the southern Persian Gulf also experienced major bleaching events in 1996, 1998, 2002, 2010, 2012 and 2017 (Riegl 1999, 2002; Riegl et al. 2011; Burt et al. 2019). Bleaching events results in direct mortality and short- and long-term negative effects on growth and reproduction (Baird and Marshall 2002; Mendes and Woodley 2002; Howells et al. 2016b) and higher susceptibility to tissue loss diseases (Miller et al. 2009; Mydlarz et al. 2009; Heron et al. 2010; Muller et al. 2018). For corals in the Persian Gulf, energy depletion resulting from repeated bleaching events and chronic multiple stressors could result in a weakening of the immune system, reducing the capacity of corals to resist infections and thus help explain the higher disease prevalence we found in this region.

Corals in the Persian Gulf live in an extreme environment which has led to some interesting adaptations and findings about how corals respond to thermal stress (Burt et al. 2014; Vaughan and Burt 2016). Thermal stress not only causes bleaching, but it is known to drive a variety of coral diseases, and thus, this region provides a platform to study disease at the edge of corals’ thermal range. Our study provides an important baseline of information on diseases that occur in this region, but much more research is needed to understand how corals and their pathogens interact under these extreme conditions. This may provide critical insights into the future of coral reefs facing the continued impacts of human stressors and global climate change.