Red Grouper (Epinephelus morio) Shape Faunal Communities via Multiple Ecological Pathways

Organisms that modify the availability of abiotic resources for other species can alter the structure and function of ecological communities through multiple pathways. In Florida Bay, red grouper (Epinephelus morio) engineer habitats by excavating sediment and detritus from karst solution holes and are also predators that consume a variety of benthic crustaceans and fish, some of which colonize engineered habitats. The effect of red grouper on these communities is complex as colonizing species interact with red grouper in different ways, including both direct (e.g., predator–prey) and indirect interactions. Here, I present the results of an experiment designed to test the direct effects of red grouper on faunal communities associated with Florida Bay solution holes by excluding red grouper from solution holes for four weeks. Red grouper presence generally had positive effects on the abundance, richness, and diversity of faunal communities associated with engineered habitats. Few strong interactions were observed between red grouper and colonizing species, mainly juvenile coral reef fishes. These results suggest that by acting as both a predator and habitat engineer, red grouper shape unique communities, distinct from those of surrounding areas, and influence the composition of communities associated with manipulated habitats.


Introduction
Habitat-modifying species generally have positive effects on the diversity of the biotic communities that are associated with those habitats [1,2]. This is true for those species that create habitat as a function of their presence (foundation species or autogenic ecosystem engineers), and for species that modify habitats mechanically (allogenic ecosystem engineers) [3,4]. Often, the net positive effect of habitat engineers on local species diversity accrues through changes in habitat availability, where engineers create or modify habitat, allowing more species to persist in a given area [1]. Because species diversity is strongly associated with community and ecosystem function, engineers that modify the availability of habitat for other species play critical roles in maintaining both the diversity and function of ecosystems [5][6][7]. The mechanisms by which habitat engineers affect local species diversity are mainly indirect: The engineer species modifies a habitat for its own purposes and in the process indirectly alters the availability and types of resources for other species. Many engineer species, both terrestrial and marine, are known to modify habitats by digging pits or burrows which are used by other species. On land, gopher tortoises (Gopherus polyphemus), native to xeric habitats in the southeastern US, dig burrows that function as shelter for at least 60 species of vertebrate and 300 species of invertebrate animals [8,9]. In the ocean, golden tilefish (Lopholatilus chamaeleonticeps) have been documented excavating burrows that are colonized by other fishes and crustaceans [10].
Conversely, predatory fish generally have negative effects on the local species diversity of biotic communities [11][12][13][14]. These effects can be particularly strong on young fishes that have recently settled to benthic habitats following their larval planktonic stage, usually referred to as recruits, and The effort to disentangle the various pathways by which species colonizing solution holes interact with red grouper will require a combination of observational and manipulative experiments. The primary goal of this study was to investigate the direct effects of red grouper on the faunal communities associated with excavated solution holes in Florida Bay by measuring the effects of red grouper presence at the community, functional group, and individual species levels. By assuming that the effects of red grouper as habitat engineer versus predator on solution hole faunal communities will occur over different time scales, a short-term manipulative experiment allowed me to isolate the direct effects of red grouper predation and agonistic behaviors on communities from the indirect effects that occur via changes in habitat availability. Red grouper were experimentally excluded from some solution holes to estimate the direct effects of red grouper presence on biotic communities over a period of days to weeks. These results were then compared to those from a multi-year observational study of solution hole habitats that was conducted concurrently and that estimated the cumulative effects of red grouper habitat manipulation on solution-hole-associated communities [29]. By comparing the results of the short-term experiment with a longer-term observational study, my goal is to more fully describe the role of red grouper in shaping these communities and to measure the various interactions that occur between red grouper and species that colonize modified solution holes.

Experimental Design
During the summers of 2011 and 2012, the effect of red grouper presence was tested on faunal communities associated with solution holes located in southwestern Florida Bay (Figure 1). A total of 24 solution holes were used across the two years. All solution holes used in the study were in water less than 4 m in depth (range = 2.4 to 3.7 m depth). This part of Florida Bay is characterized by patches of seagrasses, mainly Thalassia testudinum and Syringodium filiforme, interspersed with areas of exposed limestone hardbottom that is colonized by sponges and corals and pockmarked by solution holes. The solution holes surveyed for this study varied in areal size from 1.69 to 6.99 m 2 (mean solution hole area ± standard error (SE) = 4.44 ± 0.37). Most holes had multiple openings with one or few deeper parts that were excavated by the red grouper. The maximum excavated depth of the solution holes surveyed for this study ranged from 26 to 77 cm deep (mean excavated depth ± SE = 45.8 ± 2.70). Prior to the start of the experiment in each year, a total of 35 solution holes were surveyed to determine occupancy by red grouper: In 2011, 26 of the 35 holes were occupied by a red grouper; in 2012, red grouper occupied 24 of the 35 holes. I randomly selected 14 sites with red grouper in 2011 and 18 sites with red grouper in 2012 to be used as experimental treatment sites and assigned half of each to one of two groups: red grouper present (control treatment), and red grouper removed (exclusion treatment). Only 14 solution holes were used in 2011 due to concerns about the amount of time needed to survey all holes; additional holes were added in 2012 to maximize the experimental sample size given our surveying ability. The total sample size for each group (control and exclusion) summed across both years was 16. In 2011 the experiment was initiated during the last week of May; in 2012 the experiment was initiated during the second week in June. Red grouper that were occupying holes chosen for the exclusion treatment were captured using hook-and-line fishing gear or fish traps and were measured, tagged, and released in empty solution holes located more than 5 km from the experimental site. All solution holes were checked at least once per every 48 h to ensure that no other red groupers had moved into the hole; if a new red grouper was present, it was captured, measured, tagged, and transported to a new site as above. Divers trained in fish and invertebrate identification surveyed the faunal communities associated with solution holes once per week for four weeks. At each solution hole, a team of two scuba divers conducted a census of all fishes and motile invertebrates found within, in the watercolumn immediately above, and in a 1 m band around each hole. Fishes were counted using the pointcount method, where each diver remained stationary and counted all observable individuals (mean visual survey duration = 12 min; range = 5 to 18 min). Cryptic fish species and motile macroinvertebrates were counted by examining inside all areas of the solution hole after the pointcount survey was completed. Both divers counted all species present and the abundance of each species was determined as the maximum number of individuals observed by a single diver. Total species richness was determined as the sum of all species observed by both divers. All organisms were identified to species with two exceptions: juvenile grunts (Haemulon spp.) estimated to be < 5 cm total length (TL) and peppermint shrimp (Lysmata spp.). For both grunts and peppermint shrimp, multiple closely-related species appear very similar so that positive identification at the species level was impossible during diver surveys. At least six species of larger juvenile and sub-adult grunts were observed during the study, all of which are virtually identical to each other until they grow larger than ~5 cm TL, after which they develop distinguishing coloration. The peppermint shrimp species complex is composed of at least five closely-related species that may co-occur in Florida Bay, all of which have color patterns that were indistinguishable while scuba diving [41]. A complete set of visual surveys at all control and treatment holes generally took two days, where half of the holes were visited on one day and the other half on the next day. Red grouper occupancy at a given solution hole was determined by the positive identification of an individual within a solution hole during a Divers trained in fish and invertebrate identification surveyed the faunal communities associated with solution holes once per week for four weeks. At each solution hole, a team of two scuba divers conducted a census of all fishes and motile invertebrates found within, in the water-column immediately above, and in a 1 m band around each hole. Fishes were counted using the point-count method, where each diver remained stationary and counted all observable individuals (mean visual survey duration = 12 min; range = 5 to 18 min). Cryptic fish species and motile macroinvertebrates were counted by examining inside all areas of the solution hole after the point-count survey was completed. Both divers counted all species present and the abundance of each species was determined as the maximum number of individuals observed by a single diver. Total species richness was determined as the sum of all species observed by both divers. All organisms were identified to species with two exceptions: juvenile grunts (Haemulon spp.) estimated to be < 5 cm total length (TL) and peppermint shrimp (Lysmata spp.). For both grunts and peppermint shrimp, multiple closely-related species appear very similar so that positive identification at the species level was impossible during diver surveys. At least six species of larger juvenile and sub-adult grunts were observed during the study, all of which are virtually identical to each other until they grow larger than~5 cm TL, after which they develop distinguishing coloration. The peppermint shrimp species complex is composed of at least five closely-related species that may co-occur in Florida Bay, all of which have color patterns that were indistinguishable while scuba diving [41]. A complete set of visual surveys at all control and treatment holes generally took two days, where half of the holes were visited on one day and the other half on the next day. Red grouper occupancy at a given solution hole was determined by the positive identification of an individual within a solution hole during a diver survey, and both divers took care to closely observe all parts of the solution hole to determine if a red grouper was present during the survey.

Statistical Analysis
In addition to the total abundance of all fishes and motile invertebrates present at each solution hole, I calculated diversity metrics based on Hill's numbers for each community observation using the equation: where p i is the relative proportion of the community made up by species i [42]. When evaluated for integer values for a of 0, 1, and 2, H a reduces to species richness, the antilog of the Shannon-Weiner index, and the reciprocal of Simpson's index, respectfully. Generally, as a increases from 0 to 2, the Hill's index gives greater weight to more abundant species. I also calculated evenness using the equation: which is the ratio of the Simpson's reciprocal and the Shannon-Weiner antilog. Evenness represents a measure of the distribution of the abundances of each species in a community, and Hill initially proposed this version of evenness because it does not include species richness (H 0 ), thus is relatively insensitive to sample size [43]. This version of evenness converges to 1 when all species are equally abundant, and smaller values indicate more uneven communities. I included Hill's evenness here because the specific ways that red grouper may alter faunal communities are complex, so including multiple diversity indices should give enough variety of metrics to evaluate changes in faunal communities that could be attributed to red grouper. Because red grouper presence was the treatment variable, they were not included as a species when calculating community metrics. I tested for differences in each of the five biotic community metrics-total faunal abundance, Hill's numbers 0 to 2, and Hill's evenness-based on red grouper presence and time from experiment start using two-way repeated measures analysis of variance (ANOVA). Species abundances were square-root transformed to conform to assumptions of normality for ANOVA tests. I tested for differences in community structure attributable to red grouper presence with permutational multivariate analysis of variance (PERMANOVA), which is similar to traditional multivariate analysis of variance but does not require that data conform to multivariate normal distributions [19]. To quantify the contribution of individual species in driving differences in the structure of communities with and without red grouper, a similarity percentage analysis (SIMPER) was performed using the abundance of each species present at the end of the experiment. To visualize patterns of community structure between solution hole communities at the start and end of the experiment, I used non-metric multi-dimensional scaling (NMDS; [44]) to ordinate the data. All NMDS ordinations were performed using raw abundance data and the default options set by the "metaMDS" function of the "vegan" R package [45].
The strength of species interactions between red grouper and solution-hole-associated fauna was measured with Paine's index (PI; [46]). This interaction index calculates the per capita interaction strength for all species using the equation: for each species i, where N is the mean prey abundance in the presence of a predator, D is prey abundance in the absence of a predator, and Y is a measure of predator abundance [47]. Since red grouper are generally solitary and were never observed together in holes during this experiment, Y = 1 for all replicates. Paine's index was initially derived to quantify the effect of a consumer when the absence of said consumer resulted in a monoculture of the competitive dominant prey but has since been used to quantify the distribution of per capita effects among community members [47]. I used it here to determine the distribution of interaction strengths within the community and to identify which species have strong interactions with red grouper. I followed the methods described in Paine [46] where PI values for each species are estimated using experimental replicates for N and a mean value for D, then constructing bootstrapped confidence intervals, given that the relatively small sample size used here may underestimate the variability in D estimate. Bootstrapped 95% confidence intervals were calculated using the "resample" package in R [48]. Finally, to estimate the effect of red grouper presence on different functional groups of organisms, I calculated standardized effect sizes using Hedge's g with the equation where m with is the mean group abundance with red grouper, m without is the mean group abundance without red grouper, and s pooled is the pooled within-group variance [49,50]. The use of standardized effect scores has previously been used to compare across studies that use different methodologies or that study different organisms [50,51]. For this study, I wanted a way to compare the effect of red grouper on the various functional groups present in experimental communities, and to compare between these experimental results and those from a multi-year observational study that was conducted concurrently at alternate sites in Florida Bay (see [29]). Faunal abundance data were not normally distributed, so bootstrapping was used to nonparametrically estimate effect sizes and confidence intervals. Functional groups were constructed for diet and habitat preference. Diet groups were based on reported diets from Fishbase and included two groups of predators-piscivores and invertivores (fishes known to consume mostly motile invertebrates)-plus benthivores (fishes known to consume mostly benthic prey), planktivores, and herbivores ( [52]; see Table A1). Habitat preferences were based upon where fish species were observed relative to solution holes: Species that were always found inside holes ("demersal"), species found milling about above or adjacent to holes ("water-column"), or species that only visited holes temporarily and were not observed milling about or are known to have a more pelagic life history ("transient"). Motile invertebrates were grouped based on whether they are known prey of red grouper or not (see Table A2). I also included a group of "cleaners"-organisms known to consume ectoparasites from other species-which included both fish and invertebrate species. Effect sizes and bootstrapped 95% confidence intervals were calculated using the "bootES" package in R [53].

Results
A total of 51 species were documented during the two years of red grouper exclusion experiments: 42 species of fish and nine species of motile macro-invertebrates (see Appendix A). Based on the repeated measures ANOVA analysis using both years (2011 and 2012) combined, red grouper presence had a significant positive effect on total abundance (F 1,143 = 58.9; p = 0.003; Figure 2   Differences in community structure before and after the four-week exclusion experiments were visualized with NMDS plots, comparing communities found at control solution holes with red grouper to those found at exclusion holes without red grouper ( Figure 4). The NMDS ordination of the faunal communities at the start of the experiment shows clear overlap between all solution holes prior to treatment assignments, which was confirmed by PERMANOVA analysis (Pseudo-F1,24 = 0.737, p = 0.599). Although the NMDS plot of faunal communities after four weeks appears to show separation between control and exclusion communities, PERMANOVA analysis results did not support a significant effect of red grouper presence on faunal community structure (Pseudo-F1,30 = 1.32, p = 0.233). However, year did have a significant effect on community structure (Pseudo-F1,30 = 8.39, p = 0.001), suggesting differences between communities between the two study years. SIMPER analysis results revealed the five most abundant species-small juvenile grunts, white grunts (Haemulon plumierii), Caribbean spiny lobster, French grunts (Haemulon flavolineatum), and Pederson's cleaner shrimp (Ancylomenes pedersoni) accounted for 84.2% of the difference in community structure (Table 1). Of the top 20 species ranked by SIMPER contribution, only four species were less abundant in holes with red grouper compared to their abundance in the exclusion treatment holes-channel clinging crabs (Mithrax spinosissimus), gray angelfish (Pomacanthus arcuatus), spotted cleaner shrimp  (Simpson; B), and Hill's evenness (C) of faunal communities associated with solution holes in Florida Bay that were occupied by red grouper (control treatment, "CON"; n = 16) versus holes where the red grouper was experimentally removed (removal treatment, "REM"; n = 16) Error bars are ± 1 SE.
Differences in community structure before and after the four-week exclusion experiments were visualized with NMDS plots, comparing communities found at control solution holes with red grouper to those found at exclusion holes without red grouper ( Figure 4). The NMDS ordination of the faunal communities at the start of the experiment shows clear overlap between all solution holes prior to treatment assignments, which was confirmed by PERMANOVA analysis (Pseudo-F 1,24 = 0.737, p = 0.599). Although the NMDS plot of faunal communities after four weeks appears to show separation between control and exclusion communities, PERMANOVA analysis results did not support a significant effect of red grouper presence on faunal community structure (Pseudo-F 1,30 = 1.32, p = 0.233). However, year did have a significant effect on community structure (Pseudo-F 1,30 = 8.39, p = 0.001), suggesting differences between communities between the two study years. SIMPER analysis results revealed the five most abundant species-small juvenile grunts, white grunts (Haemulon plumierii), Caribbean spiny lobster, French grunts (Haemulon flavolineatum), and Pederson's cleaner shrimp (Ancylomenes pedersoni) accounted for 84.2% of the difference in community structure (Table 1). Of the top 20 species ranked by SIMPER contribution, only four species were less abundant in holes with red grouper compared to their abundance in the exclusion treatment holes-channel clinging crabs (Mithrax spinosissimus), gray angelfish (Pomacanthus arcuatus), spotted cleaner shrimp (Periclimenes yucatanicus), and sand perch (Diplectrum formosum)-the other 16 species were all more abundant in holes with red grouper.       The distribution of interaction strengths between red grouper and associated species varied based on the calculations of Paine's index (PI) values for each species. Of the 51 species observed during the experiment, 28 species had non-zero PI values and 23 species had interaction strengths equal to zero, indicating that their abundance was not different between treatment groups ( Figure 5). Bootstrapped means and their corresponding standard error values indicated that seven species or species complexes had PI values with 95% confidence intervals that did not include zero ( Table 2) (Table 3); the 95% confidence intervals of effect size on all other functional groups included zero, suggesting that red grouper did not have a significant effect on the abundance of these functional groups.
based on the calculations of Paine's index (PI) values for each species. Of the 51 species observed during the experiment, 28 species had non-zero PI values and 23 species had interaction strengths equal to zero, indicating that their abundance was not different between treatment groups ( Figure 5). Bootstrapped means and their corresponding standard error values indicated that seven species or species complexes had PI values with 95% confidence intervals that did not include zero ( Table 2): juvenile white grunts (PI ± SE: 4.97 ± 2.38), porkfish (Anisotremus virginicus; 0.36 ± 0.16), gray angelfish (−0.17 ± 0.063), hogfish (−0.033 ± 0.015), peppermint shrimp (−0.178 ± 0.069), sand perch (−0.211 ± 0.091), and spotted cleaner shrimp (−0.244 ± 0.079). Standardized effect size analysis of functional groups indicated significant positive effects of red grouper presence on the abundance of only the invertivore fish group (Table 3); the 95% confidence intervals of effect size on all other functional groups included zero, suggesting that red grouper did not have a significant effect on the abundance of these functional groups.      [29]. Functional group classifications for the fishes were based on reported diet information or by the location of individuals in relation to solution holes as observed during diver surveys. Functional group classifications for the invertebrates were based on reports of species known to be prey for red grouper. The functional group "cleaners" includes both fish and invertebrate species known to consume ectoparasites. Values in bold are statistically significant from zero. Across the two years of experiments, 21 red grouper colonized one of the exclusion holes following the removal of the initial individual (N 2011 = 8; N 2012 = 13). On average, empty holes were recolonized within a week (mean days to replacement ± SE = 6.62 ± 1.14). Red grouper that moved into empty solution holes were significantly smaller than the individual that was removed immediately prior, by approximately 4 cm (mean difference ± SE = −3.91 ± 1.39; T 20 = −2.82, p = 0.011; Figure 6). Of the 21 total re-colonization events, nine were multiple re-colonization events where the removed red grouper was either the third or fourth individual to be removed from the hole during the experiment. Including the initial removals at the start of the experiments, the mean size of red grouper removed from solution holes was 45.4 cm TL (± SE = 1.25); the mean size of red grouper in control treatment holes was 48.3 cm TL (± SE = 1.75). Across the two years of experiments, 21 red grouper colonized one of the exclusion holes following the removal of the initial individual (N2011 = 8; N2012 = 13). On average, empty holes were recolonized within a week (mean days to replacement ± SE = 6.62 ± 1.14). Red grouper that moved into empty solution holes were significantly smaller than the individual that was removed immediately prior, by approximately 4 cm (mean difference ± SE = -3.91 ± 1.39; T20 = -2.82, p = 0.011; Figure 6). Of the 21 total re-colonization events, nine were multiple re-colonization events where the removed red grouper was either the third or fourth individual to be removed from the hole during the experiment. Including the initial removals at the start of the experiments, the mean size of red grouper removed from solution holes was 45.4 cm TL (± SE = 1.25); the mean size of red grouper in control treatment holes was 48.3 cm TL (± SE = 1.75).

Discussion
Evaluating multiple community-level metrics indicated that red grouper presence resulted in higher abundance, species richness, and diversity of faunal communities associated with solution holes. These results generally match the predicted effects of a habitat engineer-communities with the engineer had more members and were more diverse compared to those where the engineer was excluded [1,2]. These results also confirm those of previous observational studies of Florida Bay solution holes that also found more abundant and species-rich communities in solution holes with red grouper [28,29]. While these patterns are consistent with the predicted effects of habitat engineers, they are contrary to the general effects expected of piscine predators that have been shown to reduce the species richness of prey communities through both selective and non-selective predation [12,15]. Instead, increased abundance and richness in the presence of a predator more closely resembles the effects of keystone predators that enhance diversity by limiting herbivory [22,25]. Red grouper have multiple roles in the ecosystem, acting as both habitat-engineer and predator in these communities and so they affect solution hole communities in multiple ways. By assuming that the indirect effects of red grouper on communities will manifest via changes in habitat availability (i.e., engineering effects) and thus over longer time scales (months to years), the results presented here from short-term (days to weeks) exclusion experiments isolate the direct effects of red grouper predation and agonistic behaviors on communities. Even at the short-time scale of these experiments, red grouper had positive effects on multiple community metrics, suggesting that their role as habitat engineer may be stronger than their effect as a predator.
At the community level, red grouper presence had a positive effect on abundance and richness but variable effects on diversity, as shown by the results of enhanced diversity on H 1 (Shannon-Weiner) but no effect on H 2 (Simpson's) diversity, and depressed evenness (see Figure 3). Analyzing the same data with different Hill's indices allows for the separation of the relative contributions of species within a single mathematical framework [54]. As the value of a increases (within Ha), more weight is given to abundant species, so here the effect of red grouper presence on faunal community diversity diminishes as abundant species are given more weight [42]. Similarly, Hill's evenness (E) is the ratio of abundant species to common species and here was used as a metric to compare the structure of species assemblages across solution holes [55]. Higher species richness and lower evenness could both occur if both rare and common species increased relative to abundant ones in communities with red grouper. Again, this result runs counter to the predictions for piscine predators, which tend to reduce species richness while increasing the evenness of prey communities [12].
Somewhat surprisingly, red grouper did not show strong effects on the abundance of their prey, as measured at either the functional group and individual species levels. The prey group was dominated by Caribbean spiny lobsters, which were abundant in solution holes irrespective of the presence of red grouper. Red grouper are well-known predators of lobsters, and thus a strong, direct interaction with red grouper is not entirely unexpected [35,36]. However, the experimental evidence described here suggests this interaction is positive: spiny lobsters increased when in the presence of the red grouper, a result that is at odds with an assumed predator-prey dynamic. Spiny lobsters are primarily nocturnal and seek out crevice shelters during the day, and solution holes in Florida Bay are known to host lobsters in high abundance relative to surrounding habitats [34,56]. Through excavation, red grouper likely increase the availability of suitable habitat for spiny lobster, and the positive interaction measured here suggests that habitat availability may be more important than the increased risk of mortality that comes with cohabitating with a predator. Here, the results of this experiment included all spiny lobsters together, irrespective of their size, which could mask asymmetry in predation effects across lobster size. In 2013, I conducted additional experiments to test how the interaction between red grouper and Caribbean spiny lobsters changed across lobster ontogeny. These experiments revealed that predation by red grouper was strongest on small juvenile lobsters and decreased with increasing lobster size [57]. While the cumulative interaction between red grouper and all spiny lobsters was positive, this result was driven by the increased abundance of large juvenile and adult lobsters that overshadowed the decrease in small juveniles due to predation.
Effects of short-term red grouper exclusion at the functional group level were negligible. The only significant effect of red grouper on any functional group was on the invertebrate predators (see Table 3). However, the absolute number of invertivores encountered during the experiments totaled just three individuals across both years: one goliath grouper (Epinephelus itajara) and one green moray (Gymnothorax funebris) in 2011, and one green moray in 2012, all observed in control treatment holes with red grouper. No invertebrate predators were encountered in exclusion treatment holes in either year. Although goliath groupers and green morays could potentially have significant effects on the abundance of invertebrate prey in solution holes, given their rarity and the relatively small effect size of red grouper on all invertebrate predators measured during the observational study (see Table 3), it is unlikely that red grouper have ecologically significant effects on the abundance or presence of other invertebrate predators in this system. The fact that no other effects were detected among any of the other functional groups suggests that the exclusion experiment may have been too short to detect any significant direct effects of red grouper at the functional group level, and potentially for individual species as well.
Alternatively, the strong direct effects of red grouper could be limited to a narrow range of species, regardless of the duration of experiments. While red grouper presence generally enhanced the abundance, richness, and diversity of communities, the effects of red grouper presence on individual species were highly variable. Structural and abiotic changes that result from habitat engineering will have variable consequences for specific members of associated biotic communities, as the effects of engineer presence will depend on both species-specific responses to modified habitats and on species-specific interactions with the engineer [58]. The results of the exclusion experiment show that the presence or absence of red grouper had no effect on the abundance of most of the fish and invertebrate species present in these communities but had a strong effect on the abundance of just a few species. In fact, only two species had significant positive interactions with red grouper-white grunts and porkfish (see Table 2). Both white grunts and porkfish are reef-associated species and were mainly observed as juveniles in solution holes; most of the white grunts and porkfish counted here were all smaller than 15 cm TL. Time from the start of the exclusion experiment was a significant factor for explaining patterns of both abundance and evenness, likely driven by pulses of reef fish recruitment to solution holes. Fish recruitment in Florida Bay and the adjacent reef tract tends to increase between May and November, with peak recruitment occurring between June and August [59,60]. In both 2011 and 2012, the exclusion experiment was conducted during June and July, well within the prime settlement window for fish recruitment. The growth rate of juvenile grunts was estimated previously at 0.25 to 0.75 mm per day [59], thus a 1 cm long grunt recruit that settled to a solution hole at the start of the experiment would initially be counted as a juvenile grunt recruit (a combined group of all six grunt species due to the difficulty in visually determining the species of individuals < 5 cm TL) but later catalogued by species anywhere from six to 16 days after settlement. The growth of juvenile grunt recruits from the species complex into the individual species categories could explain the significance of time as a factor in describing patterns of species richness and evenness, as the grunt species complex differentiated from the complex into the individual species across the duration of the experiments.
The effect of red grouper on the abundance of juvenile reef fishes was particularly large according to multiple metrics analyzed here and in the observational study that was conducted concurrently at an alternate set of solution holes [29]. In addition to the significant positive interactions measured experimentally between red grouper and white grunts and porkfish, the two functional groups that included juvenile grunts-planktivores, water column fishes-both showed large effects due to red grouper presence in the observational study (see Table 3). The SIMPER analysis of abundance data found that the juvenile grunt species complex was the most influential group for influencing community structure, followed by white grunts, Caribbean spiny lobster, and French grunts (see Table 1). Together, juvenile grunts, white grunts, and French grunts accounted for 56.4% of the difference in communities with and without red grouper. Adult grunts were very rarely encountered during diver surveys; most grunts were small (< 15 cm TL) and were observed milling about above and around solution holes in mixed schools with the occasional juvenile parrotfish. Juvenile grunts have not been documented as part of red grouper diets, and red grouper had negligible effects on the abundance of piscivores in the community, based on the standardized effect size on piscivore group abundance. These results suggest that red grouper may enhance juvenile reef fish abundance via a behaviorally-mediated indirect interaction (BMII) mediated through changes in the behavior of resident and transient predators.
In 2013, this hypothesis was tested by transplanting invasive lionfish (Pterois volitans and P. miles) into solution holes with and without red grouper [33]. After six weeks there were significantly more juvenile reef fish when red grouper and lionfish were both present compared to lionfish alone. In addition, analysis of lionfish stomach contents revealed that lionfish switched from a fish-dominated diet in the absence of red grouper to a crustacean-dominated diet in the presence of red grouper. These results suggest that territorial behaviors by red grouper toward other large-bodied fishes that live in or visit solution holes may disrupt the predatory ability of these piscivores and result in enhanced juvenile reef fish survival and abundance.
Five species had small but significantly negative interactions with red grouper: hogfish, sand perch, peppermint shrimp, spotted cleaner shrimp, and gray angelfish ( Table 2). Both hogfish and sand perch are transient benthivores, which may avoid solution holes occupied by a red grouper or be deterred by agonistic territorial displays. Peppermint shrimp and spotted cleaner shrimp both belong to the "cleaners" functional group, which was numerically dominated by a third species-the Pederson's cleaner shrimp. Negative responses by the rarer members of this group could be the result of competition for space among cleaners, as both Pederson's and spotted cleaner shrimp are both strongly associated with anemones, a potentially limited resource around solution holes [61]. Angelfishes are primarily considered reef-associated species, though both juveniles and adults of all four species commonly found in south Florida were regularly encountered in Florida Bay solution holes [29]. Adult angelfishes primarily consume sponges, but juveniles also consume algae and cleaning behavior by juveniles of at least one species has been reported [52,62]. In terms of individual species interactions, three of the four species had negative interactions with red grouper, though only the interaction between gray angelfish and red grouper was significant. Adult angelfishes are relatively large bodied, and though they were most often observed milling about above and around solution holes, they would also commonly dart into holes as divers approached, thus the negative interactions observed between angelfishes and red grouper could result from competition for space within solution holes. However, although both juveniles (< 10 cm TL) and adults were observed in solution holes, they were counted together in a single category. Juvenile angelfish may benefit from red grouper presence through the same mechanisms as juvenile grunts, and further investigation into possible ontogenetic variations in this interaction would clarify the nature of this interaction.
While the community-level effects of red grouper on the abundance and diversity of solution hole communities are clear, their effects on the larger Florida Bay ecosystem remain untested. Solution holes make up a relatively small portion of the benthic habitats in Florida Bay (~3%) but are important diurnal refuges for lobsters and other invertebrates and may function similarly for juvenile reef fish before they move onto adjacent reef habitats as adults [34,56,63]. Comparing the cumulative (i.e., direct plus indirect) effects of red grouper described in Ellis et al. [29] with the direct effects of red grouper presence on solution hole faunal communities described here suggests that there are relatively few negative direct effects of red grouper presence. Instead, red grouper presence had consistently positive effects on the abundance and species diversity of faunal communities over short (days to weeks) and long (multi-year) time scales, likely attributable to increased habitat availability and BMIIs mediated by territorial behaviors. Estimating the effects of red grouper presence in a hierarchical framework (e.g., community, functional group, individual species) by using both observational and experimental methods allowed me to describe the cumulative effects of red grouper as both a habitat engineer and predator along with the direct effects that red grouper had on individual species. These results