# Metabarcoding dietary analysis of coral dwelling predatory fish demonstrates the minor contribution of coral mutualists to their highly partitioned, generalist diet

Matthieu Leray, Christopher P. Meyer, Suzanne C. Mills

► To cite this version:

Matthieu Leray, Christopher P. Meyer, Suzanne C. Mills. Metabarcoding dietary analysis of coral dwelling predatory fish demonstrates the minor contribution of coral mutualists to their highly partitioned, generalist diet. PeerJ, 2015, 3, 10.7717/peerj.1047. http://dx.doi.org/10.7717/peerj.1047

HAL Id: hal-01174184

https://univ-perp.hal.science/hal-01174184

Submitted on 8 Jul 2015

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Metabarcoding dietary analysis of coral dwelling predatory fish demonstrates the minor contribution of coral mutualists to their highly partitioned, generalist diet

Matthieu Leray<sup>1,2,3</sup>, Christopher P. Meyer<sup>3</sup> and Suzanne C. Mills<sup>1,2</sup>

Submitted 7 April 2015

Accepted 2 June 2015

Published 25 June 2015

Corresponding author
Matthieu Leray,
leray.upmc@gmail.com

Academic editor
Mónica Medina

Additional Information and Declarations can be found on page 13

DOI 10.7717/peerj.1047

cc Copyright

2015 Leray et al.

Distributed under
Creative Commons CC-BY 4.0

OPEN ACCESS

## ABSTRACT

Understanding the role of predators in food webs can be challenging in highly diverse predator/prey systems composed of small cryptic species. DNA based dietary analysis can supplement predator removal experiments and provide high resolution for prey identification. Here we use a metabarcoding approach to provide initial insights into the diet and functional role of coral-dwelling predatory fish feeding on small invertebrates. Fish were collected in Moorea (French Polynesia) where the BIOCODE project has generated DNA barcodes for numerous coral associated invertebrate species. Pyrosequencing data revealed a total of 292 Operational Taxonomic Units (OTU) in the gut contents of the arc-eye hawkfish (*Paracirrhites arcatus*), the flame hawkfish (*Neocirrhites armatus*) and the coral croucher (*Caracanthus maculatus*). One hundred forty-nine (51%) of them had species-level matches in reference libraries (>98% similarity) while 76 additional OTUs (26%) could be identified to higher taxonomic levels. Decapods that have a mutualistic relationship with *Pocillopora* and are typically dominant among coral branches, represent a minor contribution of the predators' diets. Instead, predators mainly consumed transient species including pelagic taxa such as copepods, chaetognaths and siphonophores suggesting non random feeding behavior. We also identified prey species known to have direct negative interactions with stony corals, such as *Hapalocarcinus* sp, a gall crab considered a coral parasite, as well as species of vermetid snails known for their deleterious effects on coral growth. *Pocillopora* DNA accounted for 20.8% and 20.1% of total number of sequences in the guts of the flame hawkfish and coral croucher but it was not detected in the guts of the arc-eye hawkfish. Comparison of diets among the three fishes demonstrates remarkable partitioning with nearly 80% of prey items consumed by only one predator. Overall, the taxonomic resolution provided by the metabarcoding approach highlights a highly complex interaction web and demonstrates that levels of trophic partitioning among coral reef fishes have likely been underestimated. Therefore, we strongly encourage further empirical approaches to dietary studies prior to making assumptions of trophic equivalency in food web reconstruction.

<footnote n="1"><sup>1</sup> USR 3278 CRIOBE CNRS-EPHE-UPVD, CBETM de l'Université de Perpignan, Perpignan Cedex, France</footnote>

<footnote n="2"><sup>2</sup> Laboratoire d'Excellence "CORAIL"</footnote>

<footnote n="3"><sup>3</sup> Department of Invertebrate Zoology, National Museum of Natural History, Smithsonian Institution, Washington, D.C., USA</footnote>
**Subjects** Biodiversity, Ecology, Marine Biology, Molecular Biology, Zoology

**Keywords** Coral reefs, Food web, Interactions, Invertebrates, Mutualism, Exosymbionts, Pyrosequencing, COI

## INTRODUCTION

Anthropogenic stressors are impacting all ecosystems on Earth, causing both drastic changes in the structure of communities and a reduction in biodiversity (Wright, 2005; Hoegh-Guldberg & Bruno, 2010). Predators are among the most vulnerable trophic group, and have long been known to play a crucial role in stabilizing ecosystems by generating top-down forces and trophic cascades (Paine, 1966; Paine, 1969). Yet, because all predator species are not functionally equivalent, understanding how species partition their diet as well as their ecological role in food webs have become a major focus to help predict the consequences of their decline on ecosystem services (Harley, 2011).

A detailed knowledge of a predator's diet is a key element for deciphering its ecological function. Among the numerous techniques used in the literature to characterize a predator's diet, PCR-based molecular analysis of gut contents is among the most powerful because species-diagnostic DNA fragments can be detected even after several hours of digestion (Symondson, 2002). Moreover, the availability of versatile PCR primers targeting short hypervariable DNA regions combined with a high-throughput sequencing platform now offer the possibility to characterize the dietary breadth of any predator (Pompanon et al., 2012; Leray et al., 2013a). The ecological influence of a predator may then be inferred from its dietary selectivity as well as the traits and functional role of prey consumed (Chapman et al., 2013). On land, this tool is already proving invaluable for understanding the biological control potential of insect predators (Mollot et al., 2014) and the ecological effects of large herbivores (Kowalczyk et al., 2011) and carnivores (Shehzad et al., 2012). The use of high-throughput sequencing for understanding trophic links in marine systems has been more limited to date (Leray et al., 2013a).

On coral reefs, one of the most diverse and threatened of ecosystems, predatory fish feeding on benthic invertebrates are the dominant trophic category. They often dwell within the reef framework where they feed upon diverse communities of small cryptic species that are known to perform a variety of functions including direct positive or negative interactions with stony corals, the foundation species of the coral reef ecosystem (reviewed by Stella et al., 2011). Some invertebrate taxa promote the survival and growth of corals by slowing the progression of coral diseases (Pollock et al., 2012), protecting corals against corallivores (Glynn, 1980; Glynn, 1983; McKeon & Moore, 2014; Rouzé et al., 2014), removing sediments from their coral host (Stewart et al., 2006; Stier et al., 2012) and alleviating detrimental effects of coral competitors or parasites (Stier et al., 2010). Other invertebrates have deleterious effects on corals as they are known vectors of coral diseases (Sussman et al., 2003; Williams & Miller, 2005), are parasites of stony corals (Humes, 1985; Shima, Osenberg & Stier, 2010) or feed upon coral polyps (Turner, 1994; Rotjan & Lewis, 2008; Rawlinson et al., 2011) sometimes causing extensive and widespread coral mortality (Leray et al., 2012a; Kayal et al., 2012). As a consequence, the feeding behavior of these
predatory fish may have significant cascading effects on the dynamics of stony corals and subsequently the dynamics of the whole coral reef ecosystem, but it has proven challenging to understand their ecological role.

The flame hawkfish (*Neocirrhites armatus*), arc-eye hawkfish (*Paracirrhites arcatus*) and coral croucher (*Caracanthus maculatus*) are common predatory fish species on Indo-Pacific coral reefs. They co-occur among the branches of Pocilloporids (genus *Stylophora* and *Pocillopora*), one of the most important reef building corals, along with a wide diversity of invertebrates (*Patton, 1974*; *Coles, 1980*; *Odinetz, 1983*; *Stella, Jones & Pratchett, 2010*). These invertebrates include both coral mutualistic (family: Trapeziidae and some Alpheidae) and parasitic (family: Cryptochiridae) decapod species (*Simon-Blecher & Achituv, 1997*), which are potential prey for coral dwelling fish. A field manipulation of the two *Pocilloporid* obligate species, the flame hawkfish and the coral croucher (habitat specialists), highlighted that their presence among the branches of *Pocillopora eydouxi* reduced total abundance and diversity of decapod recruits by 34% and 20% respectively (*Stier & Leray, 2014*). These predators modified the composition and abundance of key mutualists (coral crabs, genus: *Trapezia*), whose benefits to *Pocillopora* are known to be both density- and diversity-dependent (*Stier et al., 2012*). Predator removal experiments have also shown that the presence of arc-eye hawkfish decreases the density of coral associated mutualist damselfish (*Holbrook, Schmitt & Brooks, 2011*). Preliminary molecular dietary analysis using traditional cloning showed the presence of coral mutualists in the gut contents of both hawkfish species (*Leray et al., 2013b*), but sampling and sequencing effort were too limited to understand their contribution to each species' diets.

In the present study, we use a high throughput sequencing approach targeting the mitochondrial Cytochrome c. Oxidase subunit I gene (COI) (also referred to as metabarcoding approach, *Taberlet et al., 2012*) to describe the dietary breadth of these predators. The study was conducted in Moorea, French Polynesia, where an extensive library of COI DNA barcodes, including all *Pocillopora* associated species, has been built by the BIOCODE project (*Leray et al., 2012b*). Implications of each predator's feeding behavior are further discussed in light of our findings.

# METHODS

## Predator and prey collections

Twenty-five adult specimens of each of the three predator fish species were speared after sunset, which corresponds to peak feeding time for all three species (M Leray, pers. obs., 2010), in the lagoon of the North shore of Moorea on the 8th, 10th and 15th of July 2010. We limited our collections to a single site (17°28′40S; 149°50′25W, Fig. 1) where coral populations had been little impacted by the recent outbreak of the corallivorous seastar, *Acanthaster planci* (*Adjeroud et al., 2009*; *Kayal et al., 2011*; *Rouzé et al., 2015*). Adults of the flame hawkfish and coral croucher always co-occurred among *Pocillopora* branches, whereas adult arc-eye hawkfish were occasionally present. Fish were individually preserved in cold 50% ethanol *in situ* after which their digestive track was dissected within 3 h and preserved in eppendorf tubes containing 80% ethanol.
<figure><img src="image_2.png" /><figcaption>**Figure 1** Map of the study location.</figcaption></figure>

## Laboratory protocol

The total content of the digestive track of each fish was dissected and used for total genomic DNA extraction using the QIAGEN DNeasy Blood & Tissue kit. Genomic DNA was then purified using the PowerClean DNA clean-up kit (MO BIO) to remove potential PCR inhibitors. We used a single set of versatile PCR primers (mlCOIintF/jgHCO2198, *Geller et al., 2013*; *Leray et al., 2013a*) known to perform well across the diversity of marine invertebrates, to amplify a 313bp region of the mitochondrial Cytochrome c. Oxidase subunit I (COI) region from each gut content sample. Moreover, this primer set was recently shown to provide reliable estimates of relative abundance for metabarcoding benthic samples (*Leray & Knowlton, 2015*). Because predator DNA co-amplification is known to impede prey detection (*Vestheim & Jarman, 2008*), predator-specific annealing blocking primers were included at ten times the concentration of versatile primers during PCR reactions as in *Leray et al. (2013a)*. All primer sequences are provided in Table 1. The PCR cocktail and touchdown temperature profile used in this study can both be found in *Leray et al. (2013a)*. Three PCR replications per sample were generated, pooled, gel excised to ensure that all primer dimers were screened away, purified using QIAGEN® MinElute columns and eluted in 12 µl of elution buffer. PCR product concentration was measured with the Qubit® Fluorometer (Invitrogen, Carlsbad, California, USA).

We pooled equimolar amounts of the combined amplicons per individual gut content for each predator species (e.g., 25 flame hawkfish gut content samples were pooled together) and 500 ng of PCR product was used per species for library preparation for Roche 454 FLX sequencing. Amplicons were end-repaired and dA-tailed using the NEBNext Quick DNA Sample Prep Reagent Set 2 chemistry (New England BioLabs, Ipswitch, Massachusetts, USA). We then performed a ligation of 454 Multiplex Identifiers (a total of three, each one containing a recognizable sequence tag) using the FLX Titanium Rapid Library MID Adaptors Kit (Roche, Basel Switzerland). Finally, the ligated PCR product of each sample was purified using Agencourt AMPure beads (Beckman Coulter
<table><caption>**Table 1 List of primers used in this study.**</caption><thead><tr><td>Primer label</td><td>Sequence (5'-3')</td><td>Reference</td></tr></thead><tbody><tr><td>mlCOIntF</td><td>GGWACWGGWTGAACWGTWTAYCCYCC</td><td>(Leray et al., 2013a)</td></tr><tr><td>jgHCO2198</td><td>TAIACYTCIGGRTGICCRAARAAYCA</td><td>(Geller et al., 2013)</td></tr><tr><td>Narmatus_Blocking</td><td>CAAAGAATCAAAACAGGTGTTGATAAAGA-C3</td><td>(Leray et al., 2013b)</td></tr><tr><td>Parcatus_Blocking</td><td>CAAAGAATCAGAACAGATGTTGGTAAAGA-C3</td><td>(Leray et al., 2013b)</td></tr><tr><td>Cmaculatus_Blocking</td><td>CAAAGAATCAGAATAGGTGTTGGTACAGA-C3</td><td>Herein</td></tr></tbody></table>

Genomics, Danvers, Massachusetts, USA), eluted in 40 µl of TE buffer, and the three samples pooled together for sequencing at the Duke Institute for Genome Sciences and Policy (Duke University, North Carolina, USA). Note that the three samples of the present study were combined with several other samples in the same 454 run.

## Analysis of FLX sequencing data

We followed a sequence data analysis pipeline optimized for analyzing large COI datasets. The pipeline detailed in *Leray et al. (2013a)* takes advantage of the coding properties of the barcoding region to discard all dubious sequences.

The initial step denoised flowgrams using Pyronoise (Quince et al., 2011) implemented in Mothur (Schloss et al., 2009). We then further quality filtered the dataset by removing any reads that met the following criteria: shorter than 200bp; more than two mismatches in the primers sequence; any ambiguous base calls (e.g., “N”); or with any homopolymer regions longer than 8bp. Remaining sequences were subsequently aligned to a high quality reference dataset (Moorea BIOCODE barcode library) based on amino acid translations using the option “enrichAlignment” in MACSE (Ranwez et al., 2011) and all sequences with any of the following were also discarded: stop codon; frame shift; insertion; or more than three deletions. Finally, potential chimeric sequences identified using UCHIME (Edgar et al., 2011) were removed to obtain a high quality sequence dataset for downstream analysis.

To evaluate prey richness and composition, sequences were clustered in Operational Taxonomic Units (OTUs) using a Bayesian algorithm implemented in CROP (Hao, Jiang & Chen, 2011). This program delineates OTUs based on the natural distribution of sequence dissimilarity in the data and within a range of sequence similarity values defined by the user. This approach performs better for clustering sequences obtained from environmental samples than a fixed dissimilarity cutoff (e.g., 5%) because they contain a diversity of phyla that differ in their rate of COI evolution. The lower and upper bound variance were set to 3 and 4 respectively (which corresponds to 6% and 8%) as they were shown to provide the best results for marine invertebrates (*Leray et al., 2013a; Leray et al., 2013b*). Following OTU delineation, a representative sequence per OTU was used for taxonomic identification using BLAST searches in the local BIOCODE database and in GENBANK. We considered that there was a species level match when sequence similarity was at least 98% (*Machida et al., 2009; Plaisance et al., 2009*). Whenever sequence similarity was lower
<figure><img src="image_3.png" /><figcaption>**Figure 2** Rarefaction curves to evaluate the completeness of the sequencing effort at describing the diversity of dietary items in the gut contents of three coral reef fish species.</figcaption></figure>

than 98%, we used a Bayesian approach implemented in the Statistical Assignment Package (SAP, Munch et al., 2008) to assign OTUs to a higher taxonomic group. SAP conducts assignments by building 10,000 unrooted phylogenetic trees from a collection of homologue sequences retrieved from a sequence database. It then calculates the probability that a query sequence belongs to a monophyletic group within that set of homologues. Here, we allowed SAP to retrieve 50 homologues from GENBANK with >70% sequence similarity to each query sequence (i.e., each OTU representative sequence) and accepted taxonomic assignments at an 80% posterior probability cutoff. Importantly, SAP can only assign sequences to taxonomic groups that are represented in the database, as is also the case with other assignment methods. Therefore, to minimize misidentification at lower taxonomic levels, we only report assignments to the phylum, class and order levels (Appendix S1).

## RESULTS

We obtained 69,663 reads of which 54,283 high quality reads were retained for downstream analysis (arc-eye hawkfish: 24,629; flame hawkfish: 13,536; coral croucher: 16,118). The Bayesian clustering algorithm delineated 292 OTUs in the gut contents of the three predatory fish species (Appendix S1). The number of dietary items was much lower in the gut contents of the coral croucher (64 taxa) than in both arc-eye (147 taxa) and flame hawkfish (149 taxa). BLAST searches provided high-resolution taxonomic assignments (>98% similarity) for 149 OTUs (51%) (Appendix S1) and the statistical assignment approach enabled the identification of 76 additional OTUs to a higher taxonomic level (26%). 67 OTUs (22.9%) remained unidentified (labeled as "Unidentified" in Appendix S1). None of the rarefaction curves reached a plateau (Fig. 2) which indicates that further sequencing effort would be necessary for a more exhaustive dietary analysis of these predatory fish.
<figure><img src="image_4.png" /><figcaption>**Figure 3** Proportion of identified OTUs in relation to the number of sequences they represent. Whenever OTU sequence similarity to a reference barcode was <98%, we used the Phylogenetic Bayesian assignment tool implemented in SAP to assign OTUs to a higher taxonomic group.</figcaption></figure>

The diversity of dietary items spanned 25 classes belonging to 17 phyla. Malacostraca was the dominant taxonomic prey group (36.7%, 21.5% and 43.7% for the arc-eye hawkfish, the flame hawkfish and the coral croucher, respectively). The arc-eye hawkfish also consumed numerous species of Actinopterygii (17.7% total OTUs) and Maxillopoda (10.9% total OTUs). A significant proportion of the flame hawkfish and coral croucher's diet was represented by Maxillopoda (12.1% and 6.2% total OTUs, respectively) and Gastropoda (9.4% and 4.7% total OTUs, respectively). Eighteen OTUs (28%) detected in the gut contents of the coral croucher remained unidentified. Direct matches to reference barcodes (>98% similarity) were more prevalent among Actinopterygii (94.1%), Malacostraca (74.1%) and Gastropoda (79.2%) compared to Maxillopoda (40%). Moreover, direct matches were more prevalent for OTUs represented by large numbers of sequences (Fig. 3). Almost nine of ten OTUs (86.7%) matched reference barcodes if they were common in the amplicon libraries (>1,000 sequences), whereas only a third (33.8%) of the single sequences matched a reference sequence. Probability of a match increased as the number of sequences increased (1: 33.8%; [2–9]: 43.7%; [10–99]: 52.2%; [100–999]: 70%; >1,000: 86.7%; Fig. 3).

Most Malacostraca OTUs were decapods (81.5%, 46.9% and 78.6% for the arc-eye hawkfish, the flame hawkfish and the coral croucher respectively— Appendix S1). All
three predatory fish fed upon *Pocillopora* obligate decapod species, but they represent a minor fraction of the total diversity of the prey they consumed (arc-eye hawkfish: 2%; flame hawkfish: 4%, coral croucher: 9.3%). Among them, we detected five coral crab species of the genus *Trapezia* that are mutualists of *Pocillopora* (*Trapezia bidentata*, *T. serenei*, *T. rufopunctata*, *T. areolata* and *T. spp*). These mutualists also represented a minor proportion of sequences in the gut contents of the arc-eye and flame hawkfish (proportion of total sequences: 5.6% and 2.4%; proportion of decapod sequences: 9.1% and 12.7%, respectively). *Pocillopora* mutualists represented a higher proportion of the coral croucher's diet with 15.3% of the total number of sequences and 47.9% of the total number of decapod sequences.

Additional trophic links involving non-decapod prey are of particular interest for understanding the effect of predators on coral and its associated communities. Predatory fish had fed upon coral associated planktivorous damselfishes of the family Pomacentridae (*Dascyllus flavicaudus*: 0.02%, 0% and 0.12%, *Chromis viridis*: 0.01%, 0.69% and 0% of total sequences in the diet of the arc-eye hawkfish, the flame hawkfish and the coral croucher, respectively) that benefit the growth of the coral host (*Holbrook et al., 2008*). Interestingly, Anthozoa were represented by two OTUs among which the host *Pocillopora* itself accounted for 20.8% and 20.1% of total number of sequences in the guts of flame hawkfish and coral croucher, but was completely absent from the gut of the arc-eye hawkfish. On the other hand, *Hapalocarcinus sp*, a gall crab considered a coral parasite, was recovered in the diet of both the arc-eye and flame hawkfish. Both hawkfish had also consumed vermetid snails known for their deleterious effects on coral growth (*Shima, Osenberg & Stier, 2010*). *Harpiliopsis beaupresii*, a caridean shrimp associated with *Pocillopora* but whose function is unknown, was also detected in the gut contents of the coral croucher. Almost 10 percent (8.3%) of the coral croucher's diet is composed of two snails (*Drupa ricinus* and *Pascula muricata*). Finally, predators had also consumed pelagic taxa including members of Maxillipoda, Chaetognatha and Hydrozoa (*Appendix S1*).

Prey species were remarkably partitioned among predators (Fig. 4). Almost eighty percent (79.5%) of prey species had been consumed by only one predator species (232 of 292). Eighteen percent (N = 52) were found in two predator diets and only eight prey species (>3%) had been ingested by all three predatory fish species analyzed. Of the shared components, the arc-eye hawkfish and the coral croucher had consumed 14 taxa in common among which six were Malacostraca. The arc-eye and flame hawkfish shared 29 prey taxa with a majority of Actinopterygii and Malacostraca. Prey sharing was lowest (nine OTUs; of which six were Malacostraca) between the two species that were always found co-occuring together in the coral host, the flame hawkfish and coral croucher. Analyses that included only prey OTUs consisting of >1% of either of the three species diets according to the relative abundance of reads demonstrate even greater partitioning (Fig. 4). Only six of the sixty-six prey items were shared at a proportion greater than 1% in any two fish species diets, and no prey species were shared among all three. Of the 66 prey items making up at least 1% of any diet, nine out of ten were consumed by only one predator.
<figure><img src="image_5.png" /><figcaption>**Figure 4** Dietary partitioning among the three predatory fish species. Left neighbor-joining phylogeny using LogDet distance model based on a constraint topology of major clades represents relationship among the 66 prey OTUs that comprise >1% of any one species diet. Thickness of linkages to right represents relative proportion of predatory diets. Six shared species are highlighted with circles. Fish images courtesy of D. Liittschwager. The 66 OTUs are highlighted in Appendix S1.</figcaption></figure>
DISCUSSION

Dietary analysis can be a powerful approach to gain insights into the ecological role of reef-dwelling predatory fish, but low taxonomic resolution in prey identification often obscures the complexity of trophic links (Longenecker, 2007). For example, the diet of the arc-eye hawkfish, flame hawkfish and coral croucher previously described from morphological identification of prey remains in gut contents was considered to be simply composed of small benthic crustaceans (class: Malacostraca) (Bachet, Zysman & Lefevre, 2007). Preliminary DNA analysis using traditional cloning revealed a breadth of prey species in the guts of the arc-eye and flame hawkfish, the majority of which were crustaceans 18 of 24 (75%) and 21 of 31 (68%) respectively (Leray et al., 2013b). This study highlights that a metabarcoding approach significantly increases the taxonomic scope by documenting an even broader taxonomic distribution of species consumed by hawkfishes (Appendix S1). The coral croucher diet also includes a wide spectrum of prey demonstrating that all three predatory fish feed broadly across community diversity. Our results highlight the importance of collecting empirical dietary data to understand processes of species coexistence in this high diversity marine ecosystem.

The ecological influence of a predator is contingent upon the prey it consumes. Their feeding behavior may induce cascading effects that will depend on the type of association that the prey they consume (or interfere with) have with keystone species. For example, in terrestrial ecosystems where up to 90% of flowering plant species use animal pollinators for reproduction (Bushmann & Nabhan, 1996), a predator's effect on plant reproductive success, growth and survival will depend on its relative consumption of pollinators and phytophageous insects (Dukas, 2005; Knight et al., 2006). Similarly, some coral reef dwelling predatory fish may either disrupt benefits to corals if they derive a significant proportion of their diet from coral mutualists or alternatively alleviate deleterious effects on corals if they consume coral parasites. Invertebrate communities occurring among the branches of live *Pocillopora* corals in Moorea or elsewhere in the Pacific are typically composed of a preponderance of decapod mutualists (>80% of diversity and abundance in live *Pocillopora*—see Patton, 1974; Coles, 1980; Odinetz, 1983; Stella, Jones & Pratchett, 2010; Leray et al., 2012a). Based on previous cloning studies (Leray et al., 2013b) only the arc-eye hawkfish consumed functionally important prey (*Trapezia tigrina*). With increased sequencing depth herein, we demonstrate that while many other mutualist decapod species do occur in the diets of the arc-eye hawkfish, flame hawkfish and coral croucher (5.6%, 2.4% and 15.3% of sequence abundance, respectively; Appendix S1), they represent a much smaller proportion of the diet than would be expected from their density in natural communities. Interestingly, we found evidence of the *Pocillopora* obligate pontoniid shrimp *Harpiliopsis beaupressi* but no detection of congeneric *H. depressa* and *H. spinigera* in the predators' gut contents, despite their very high abundance reported on head-size *Pocillopora* in Moorea (Leray et al., 2012a). It is also surprising not to discover *Alpheus lottini* in the diets of the three species although this is a common species found in all living *Pocillopora* observed and known to have beneficial effects on coral survivorship (Stier et al., 2012). Overall, our data indicate a non random pattern of prey consumption atypical of
an opportunistic feeding behavior (where prey would be consumed in proportion to their abundance—Heinlein, Stier & Steele, 2010) which suggests the outcome of coevolutionary dynamics between *Pocillopora* associated predator and prey.

Nevertheless, while our metabarcoding dietary analysis suggests limited predation pressure on mutualists, a four-month recruitment experiment conducted on the North shore of Moorea in 2009 showed a lower abundance of mutualists in corals where the coral croucher and the flame hawkfish occurred (Stier & Leray, 2014), a pattern that may be driven by non-consumptive effects of predators. For example, competent larvae may preferentially settle on corals where predators are absent. Regardless of the mechanisms, such predator effects have important implications for coral performance, because density and composition of mutualist assemblages are known to be important for the quality of the services provided to their host (Stier et al., 2012; Rouzé et al., 2014).

In addition, our metabarcoding analyses of gut contents revealed for the first time predation on a gall crab (*Hapalocarcinus sp*) and vermetid snails (genus: *Dendropoma*), which are considered detrimental to the coral host (*Simon-Blecher & Achituv, 1997; Shima, Osenberg & Stier, 2010*). Vermetid snails are particularly prevalent in Moorea where they can reduce coral growth by up to 81% and survival by up to 52% (*Shima, Osenberg & Stier, 2010*). Predation on parasites may compensate for the negative effects of the reduction in density of decapod mutualists in corals facing environmental stressors. We also recovered a significant proportion of sequences belonging to *Pocillopora* from the flame hawkfish and the coral croucher gut contents, which suggest that these predatory fish also feed on mucus released by their biogenic habitat. The absence of *Symbiodinium* COI sequences from our dataset also supports the consumption of mucus rather than coral polyps. Alternatively, *Pocillopora* DNA may have been sufficiently abundant and well preserved in the gut contents of mucus feeding prey (e.g., Trapeziidae) to be co-amplified (Harwood et al., 2001; Sheppard et al., 2005). Importantly though, *Pocillopora* was completely absent from the arc-eye hawkfish diet which also includes Trapeziid species, suggesting minimal secondary consumption or associated eDNA amplification. Overall, high-resolution dietary data are revealing a highly complex interaction web with very specialized functional roles played by each species. This highlights the shortcomings of the functional groups approach commonly used to evaluate redundancy and complementarity among coral reef species (Naeem & Wright, 2003; Micheli & Halpern, 2005).

Fine-scale spatial partitioning commonly occurs among coral reef fish species (Robert-
son & Lassig, 1980; Waldner & Robertson, 1980; Ebersole, 1985; Bouchonnavaro, 1986;
Munday, Jones & Caley, 1997; Depczynski & Bellwood, 2004; Gardiner & Jones, 2005) but the
extent of food partitioning remains controversial (Longenecker, 2007). In fact, most early
work investigating differences in diet among reef fish species showed high levels of diet
overlap (Hiatt & Strasburg, 1960; Randall, 1967; Hobson, 1974; Talbot, Russell & Anderson,
1978; Harmelin-Vivien, 1979; Anderson et al., 1981; Bouchonnavaro, 1986; Ross, 1986;
Depczynski & Bellwood, 2003; Kulbicki et al., 2005; Longenecker, 2007; Castellanos-Galindo
& Giraldo, 2008) which has led many to the conclusion that trophic partitioning was
not a mechanism promoting species coexistence on coral reefs. However, these studies,
which rely on morphological identification of semi-digested prey remains in gut contents grouped food items into broad categories therefore impeding accurate measures of partitioning (Longenecker, 2007). Alternative strategies such as field observations of feeding behavior (Pratchett, 2005; Pratchett, 2007; Pratchett & Berumen, 2008) or a combination of gut content and stable isotope analyses (Ho et al., 2007; Nagelkerken et al., 2009) helped describe dietary differences between closely related species, but generalizations about the importance of trophic partitioning for the maintenance of coral reef diversity remain difficult. In the present study, high-resolution molecular data highlight an unexpected level of dietary partitioning among the three study species. While both hawkfish species are from the same family (Cirrhitidae), they share only a single prey item at greater than 1% of either of their diets (Trapezia serenei). There is also a minor dietary overlap between the coral croucher (family Caracanthidae) and the flame hawkfish that were always found co-occurring in Pocillopora and are known to rarely venture outside the branching structure provided by their host coral (Hiatt & Strasburg, 1960; Stier & Leray, 2014). These results demonstrate that levels of trophic partitioning have likely been underestimated. We strongly encourage further empirical approaches to dietary studies prior to making assumptions of trophic equivalency in food web reconstruction (Leibold & McPeek, 2006).

The extent to which secondary prey co-amplification could lead to errors in food web analysis has not been quantified in marine systems (see Sheppard et al., 2005 for an example in a terrestrial system). In the present dataset, numerous prey species identified in fish gut contents are either grazers or detritivores (e.g., isopods, amphipods, decapods, ophiuroids and gastropods) and are therefore unlikely to consume each other. Some fish species detected in the gut contents are higher-level predators (e.g., *Caranx melampygus*) that could consume benthic grazers and detritivores as adults, but they were most likely fed upon at a younger developmental stage (egg, larva or juvenile) given the size of predators. Demospongiae, Ascidiacea and Gymnolaemata represented by few or a single sequence in the dataset were, however, possibly ingested unintentionally as secondary prey or epiphytes on the carapace of spider crabs (e.g., *Menaethius monoceros* and *Perinia tumida*). Parasites of prey (e.g., parasitic isopods of coral crabs of the genus *Trapezia*, Appendix S1) and parasites of a predator's digestive track (e.g., Trematoda and Cestoda) may also confound food web reconstructions and care should be taken to consider the targeted roles these fish predators have on various parasites. The recovery of secondary prey may artificially inflate dietary partitioning if those lower levels are also partitioned. However, we expect the amount of DNA that these secondary prey items represent in the guts of our target predators should be minor and highly digested in comparison to primary prey. A recent metabarcoding analysis of benthic samples (*Leray & Knowlton, 2015*) showed evidence of a correlation between amount of DNA and number of reads. Thus if secondary prey is quickly degraded, those taxa should be represented by one or few reads only. The present dataset shows minor dietary overlap both with and without rare OTUs (<1% of total OTUs, Fig. 4), further supporting our conclusions regarding the extent of trophic partitioning among all three fish species.
Importantly, our analysis shows that in-depth sequencing would enable a more comprehensive representation of trophic links in this multi-faceted ecosystem. Additional reads would provide more OTUs matching reference barcodes (in GENBANK, BOLD or BIOCODE) but also a higher proportion of unidentified OTUs represented by a single sequence (“singleton”, Fig. 3) that are likely to be either (1) small taxa underrepresented in DNA barcode libraries (Leray et al., 2013a), or (2) the product of sequencing artifacts despite our very stringent quality filtering based on amino-acid translation. Further barcoding initiatives aiming to catalogue small life forms (e.g., meiofauna) will be crucial to advance our understanding of food webs. Systematic removal of singletons may also be used as a conservative measure, although most of them likely represent valid taxa (Huse et al., 2010). As coral reef ecosystems decline worldwide, understanding the role of predator species in a dominant, yet largely understudied trophic category, is essential. Our study highlights the tremendous potential of metabarcoding as an approach to provide unprecedented taxonomic resolution in the diet of coral dwelling predatory fish. We encourage that further work should be conducted to understand the ecological role of reef dwelling fish and invertebrates.

## ACKNOWLEDGEMENTS

We thank Gustav Paulay, Arthur Anker, Joseph Poupin and the BIOCODE teams who collected both marine and terrestrial specimens, the “Centre de Recherche Insulaire et Observatoire de l’Environnement (CRIOBE) de Moorea” and the Richard B. Gump field station in Moorea for logistical support.

## ADDITIONAL INFORMATION AND DECLARATIONS

### Funding

Funding was provided by the Gordon and Betty Moore Foundation, France American Cultural Exchange program (FACE)/Partner University Fund (PUF), the Smithsonian Institution fellowship program and the Agence National de Recherche, ANR-11-JSV7-012-01 Live and Let Die. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

### Grant Disclosures

The following grant information was disclosed by the authors:
Gordon and Betty Moore Foundation.
France American Cultural Exchange program (FACE)/Partner University Fund (PUF).
Smithsonian Institution fellowship program and the Agence National de Recherche:
ANR-11-JSV7-012-01.

### Competing Interests

Matthieu Leray is a postdoctoral fellow at the National Museum of Natural History, Smithsonian Institution. Christopher P Meyer is a research zoologist at the Department of Invertebrate Zoology, National Museum of Natural History, Smithsonian Institution.
## Author Contributions

* Matthieu Leray conceived and designed the experiments, performed the experiments, analyzed the data, contributed reagents/materials/analysis tools, wrote the paper, prepared figures and/or tables, reviewed drafts of the paper.

* Christopher P. Meyer conceived and designed the experiments, contributed reagents/materials/analysis tools, prepared figures and/or tables, reviewed drafts of the paper.

* Suzanne C. Mills conceived and designed the experiments, contributed reagents/materials/analysis tools, reviewed drafts of the paper.

## Animal Ethics

The following information was supplied relating to ethical approvals (i.e., approving body and any reference numbers):

Approval was granted from our institutional animal ethics committee, le Centre National de la Recherche Scientifique (CNRS), for sacrificing and subsequently dissecting fish (Permit Number: 006725). None of the fish species are on the endangered species list and no specific authorization was required from the French Polynesian government for collection.

## DNA Deposition

The following information was supplied regarding the deposition of DNA sequences:

An alignment of all OTU representative sequences is provided in Appendix S2 and the denoised sequence dataset is deposited in the Dryad Repository DOI 10.5061/dryad.v0p71).

## Supplemental Information

Supplemental information for this article can be found online at http://dx.doi.org/10.7717/peerj.1047#supplemental-information.

## REFERENCES

* Adjeroud M, Michonneau F, Edmunds PJ, Chancerelle Y, de Loma TL, Penin L, Thibaut L, Vidal-Dupiol J, Salvat B, Galzin R. 2009. Recurrent disturbances, recovery trajectories, and resilience of coral assemblages on a South Central Pacific reef. *Coral Reefs* 28:775–780 DOI 10.1007/s00338-009-0515-7.

* Anderson GRV, Ehrlich AH, Ehrlich PR, Roughgarden JD, Russell BC, Talbot FH. 1981. The community structure of coral reef fishes. *American Naturalist* 117:476–495 DOI 10.1086/283729.

* Bachet P, Zysman T, Lefevre Y. 2007. Guide des poissons de Tahiti et de ses îles. In: *Au vent des îles, collection nature et environnement d'Oceanie*.

* Bouchonnavaro Y. 1986. Partitioning of food and space resources by chaetodontid fish on coral reefs. *Journal of Experimental Marine Biology and Ecology* 103:21–40 DOI 10.1016/0022-0981(86)90130-9.

* Bushmann S, Nabhan G. 1996. *The forgotten pollinators*. Washington, DC: Island Press.
**Castellanos-Galindo GA, Giraldo A. 2008.** Food resource use in a tropical eastern Pacific tidepool fish assemblage. *Marine Biology* 153:1023–1035 DOI 10.1007/s00227-007-0874-y.

**Chapman EG, Schmidt JM, Welch KD, Harwood JD. 2013.** Molecular evidence for dietary selectivity and pest suppression potential in an epigeal spider community in winter wheat. *Biological Control* 65:72–86 DOI 10.1016/j.biocontrol.2012.08.005.

**Coles SL. 1980.** Species diversity of decapods associated with living and dead reef coral *Pocillopora meandrina*. *Marine Ecology-Progress Series* 2:281–291 DOI 10.3354/meps002281.

**Depczynski M, Bellwood DR. 2003.** The role of cryptobenthic reef fishes in coral reef trophodynamics. *Marine Ecology-Progress Series* 256:183–191 DOI 10.3354/meps256183.

**Depczynski M, Bellwood DR. 2004.** Microhabitat utilisation patterns in cryptobenthic coral reef fish communities. *Marine Biology* 145:455–463 DOI 10.1007/s00227-004-1342-6.

**Dukas R. 2005.** Bumble bee predators reduce pollinator density and plant fitness. *Ecology* 86:1401–1406 DOI 10.1890/04-1663.

**Ebersole JP. 1985.** Niche separation of two damselfish species by aggregation and differential microhabitat utilization. *Ecology* 66:14–20 DOI 10.2307/1941302.

**Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. 2011.** UCHIME improves sensitivity and speed of chimera detection. *Bioinformatics* 27:2194–2200 DOI 10.1093/bioinformatics/btr381.

**Gardiner N, Jones G. 2005.** Habitat specialisation and overlap in a guild of coral reef cardinalfishes (Apogonidae). *Marine Ecology Progress Series* 305:163–175 DOI 10.3354/meps305163.

**Geller JB, Meyer CP, Parker M, Hawk H. 2013.** Redesign of PCR primers for mitochondrial Cytochrome c. oxidase subunit I for marine invertebrates and application in all-taxa biotic surveys. *Molecular Ecology Resources* 13:851–861 DOI 10.1111/1755-0998.12138.

**Glynn PW. 1980.** Defense by symbiotic crustacea of host corals elicited by chemical cues from predator. *Oecologia* 47:287–290 DOI 10.1007/BF00398518.

**Glynn PW. 1983.** Increased survivorship in corals harboring crustacean symbionts. *Marine Biology Letters* 4:105–111.

**Hao X, Jiang R, Chen T. 2011.** Clustering 16S rRNA for OTU prediction: a method of unsupervised Bayesian clustering. *Bioinformatics* 27:611–618 DOI 10.1093/bioinformatics/btq725.

**Harley CDG. 2011.** Climate change, keystone predation, and biodiversity loss. *Science* 334:1124–1127 DOI 10.1126/science.1210199.

**Harmelin-Vivien ML. 1979.** *Ichthyofaune des recifs coralliens de Tulear (Madagascar): ecologie et relations trophiques*. Aix-Marseilles, France: University of Aix-Marseilles.

**Harwood JD, Phillips SW, Sunderland KD, Symondson WOC. 2001.** Secondary predation: quantification of food chain errors in an aphid-spider-carabid system using monoclonal antibodies. *Molecular Ecology* 10:2049–2057 DOI 10.1046/j.0962-1083.2001.01349.x.

**Heinlein JM, Stier AC, Steele MA. 2010.** Predators reduce abundance and species richness of coral reef fish recruits via non-selective predation. *Coral Reefs* 29:527–532 DOI 10.1007/s00338-010-0592-7.

**Hiatt RW, Strasburg DW. 1960.** Ecological relationships of the fish fauna on coral reefs on the Marshall Islands. *Ecological Monographs* 30:66–127 DOI 10.2307/1942181.

**Ho CT, Kao SJ, Dai CF, Hsieh HL, Shiah FK, Jan RQ. 2007.** Dietary separation between two blennies and the Pacific gregory in northern Taiwan: evidence from stomach content and stable isotope analyses. *Marine Biology* 151:729–736 DOI 10.1007/s00227-006-0517-8.

**Hobson ES. 1974.** Feeding relationships of teleostean fishes on coral reefs in Kona, Hawaii. *Fishery Bulletin* 72:915–1031.
**Hoegh-Guldberg O, Bruno JF. 2010.** The impact of climate change on the world's marine ecosystems. *Science* **328**:1523-1528 DOI 10.1126/science.1189930.

**Holbrook SJ, Brooks AJ, Schmitt RJ, Stewart HL. 2008.** Effects of sheltering fish on growth of their host corals. *Marine Biology* **155**:521-530 DOI 10.1007/s00227-008-1051-7.

**Holbrook SJ, Schmitt RJ, Brooks AJ. 2011.** Indirect effects of species interactions on habitat provisioning. *Oecologia* **166**:739-749 DOI 10.1007/s00442-011-1912-5.

**Humes A. 1985.** A review of the Xarifiidae (Copepoda, Poecilostomatoida), parasites of scleractinian corals in the Indo-Pacific. *Bulletin of Marine Science* **36**:467-632.

**Huse SM, Welch DM, Morrison HG, Sogin ML. 2010.** Ironing out the wrinkles in the rare biosphere through improved OTU clustering. *Environmental Microbiology* **12**:1889-1898 DOI 10.1111/j.1462-2920.2010.02193.x.

**Kayal M, Lenihan HS, Pau C, Penin L, Adjeroud M. 2011.** Associational refuges among corals mediate impacts of a crown-of-thorns starfish *Acanthaster planci* outbreak. *Coral Reefs* **30**:827-837 DOI 10.1007/s00338-011-0763-1.

**Kayal M, Vercelloni J, Lison de Loma T, Bosserelle P, Chancerelle Y, Geoffroy S, Stievenart C, Michonneau F, Penin L, Planes S, Adjeroud M. 2012.** Predator crown-of-thorns starfish (*Acanthaster planci*) outbreak, mass mortality of corals, and cascading effects on reef fish and benthic communities. *PLoS ONE* **7**:e47363 DOI 10.1371/journal.pone.0047363.

**Knight TM, Chase JM, Hillebrand H, Holt RD. 2006.** Predation on mutualists can reduce the strength of trophic cascades. *Ecology Letters* **9**:1173-1178 DOI 10.1111/j.1461-0248.2006.00967.x.

**Kowalczyk R, Taberlet P, Coissac E, Valentini A, Miquel C, Kamiński T, Wójcik JM. 2011.** Influence of management practices on large herbivore diet—Case of European bison in Białowieża Primeval Forest (Poland). *Forest Ecology and Management* **261**:821-828 DOI 10.1016/j.foreco.2010.11.026.

**Kulbicki M, Bozec Y-M, Labrosse P, Letourneur Y, Mou-Tham G, Wantiez L. 2005.** Diet composition of carnivorous fishes from coral reef lagoons of New Caledonia. *Aquatic Living Resources* **18**:231-250 DOI 10.1051/alr:2005029.

**Leibold MA, McPeek MA. 2006.** Coexistence of the niche and neutral perspectives in community ecology. *Ecology* **87**:1399-1410 DOI 10.1890/0012-9658(2006)87[1399:COTNAN]2.0.CO;2.

**Leray M, Agudelo N, Mills SC, Meyer CP. 2013b.** Effectiveness of annealing blocking primers versus restriction enzymes for characterization of generalist diets: unexpected prey revealed in the gut contents of two coral reef fish species. *PLoS ONE* **8**:e58076 DOI 10.1371/journal.pone.0058076.

**Leray M, Beraud M, Anker A, Chancerelle Y, Mills S. 2012a.** *Acanthaster planci* outbreak: decline in coral health, coral size structure modification and consequences for obligate decapod assemblages. *PLoS ONE* **7**:e35456 DOI 10.1371/journal.pone.0035456.

**Leray M, Boehm JT, Mills SC, Meyer CP. 2012b.** Moorea BIOCODE barcode library as a tool for understanding predator-prey interactions: insights into the diet of common predatory coral reef fishes. *Coral Reefs* **31**:383-388 DOI 10.1007/s00338-011-0845-0.

**Leray M, Knowlton N. 2015.** DNA barcoding and metabarcoding of standardized samples reveal patterns of marine benthic diversity. *Proceedings of the National Academy of Sciences of the United States of America* **112**:2076-2081 DOI 10.1073/pnas.1424997112.

**Leray M, Yang JY, Meyer CP, Mills SC, Agudelo N, Ranwez V, Boehm JT, Machida RJ. 2013a.** A new versatile primer set targeting a short fragment of the mitochondrial COI region for metabarcoding metazoan diversity: application for characterizing coral reef fish gut contents. *Frontiers in Zoology* **10**:34 DOI 10.1186/1742-9994-10-34.
**Longenecker K. 2007.** Devil in the details: high-resolution dietary analysis contradicts a basic assumption of reef-fish diversity models. *Copeia* **3**:543–555
DOI 10.1643/0045-8511(2007)2007[543:DITDHD]2.0.CO;2.

**Machida RJ, Hashiguchi Y, Nishida M, Nishida S. 2009.** Zooplankton diversity analysis through single-gene sequencing of a community sample. *BMC Genomics* **10**:438
DOI 10.1186/1471-2164-10-438.

**McKeon CS, Moore JM. 2014.** Species and size diversity in protective services offered by coral guard-crabs. *PeerJ* **2**:e574 DOI 10.7717/peerj.574.

**Micheli F, Halpern BS. 2005.** Low functional redundancy in coastal marine assemblages. *Ecology Letters* **8**:391–400 DOI 10.1111/j.1461-0248.2005.00731.x.

**Mollot G, Duyck P-F, Lefeuvre P, Lescourret F, Martin J-F, Piry S, Canard E, Tixier P. 2014.** Cover cropping alters the diet of arthropods in a banana plantation: a metabarcoding approach. *PLoS ONE* **9**:e93740 DOI 10.1371/journal.pone.0093740.

**Munch K, Boomsma W, Huelsenbeck JP, Willerslev E, Nielsen R. 2008.** Statistical assignment of DNA sequences using Bayesian phylogenetics. *Systematic Biology* **57**:750–757
DOI 10.1080/10635150802422316.

**Munday PL, Jones GP, Caley MJ. 1997.** Habitat specialisation and the distribution and abundance of coral-dwelling gobies. *Marine Ecology-Progress Series* **152**:227–239 DOI 10.3354/meps152227.

**Naeem S, Wright JP. 2003.** Disentangling biodiversity effects on ecosystem functioning: deriving solutions to a seemingly insurmountable problem. *Ecology Letters* **6**:567–579
DOI 10.1046/j.1461-0248.2003.00471.x.

**Nagelkerken I, Van der Velde G, Wartenbergh SLJ, Nugues MM, Pratchett MS. 2009.** Cryptic dietary components reduce dietary overlap among sympatric butterflyfishes (Chaetodontidae). *Journal of Fish Biology* **75**:1123–1143 DOI 10.1111/j.1095-8649.2009.02303.x.

**Odinetz MO. 1983.** *Ecologie et structure des peuplements de crustaces decapodes associes aux coraux du genre Pocillopora en Polynesie Francaise et en Micronesie (Guam)*. Paris: Universite Pierre et Marie Curie.

**Paine RT. 1966.** Food web complexity and species diversity. *American Naturalist* **100**:65–75
DOI 10.1086/282400.

**Paine RT. 1969.** A note on trophic complexity and community stability. *American Naturalist* **103**:91–93 DOI 10.1086/282586.

**Patton WK. 1974.** Community structure among the animals inhabiting the coral *Pocillopora damicornis* at Heron Island, Australia. In: Vernberg WB, ed. *Symbiosis in the sea*. Columbia: University of South Carolina Press.

**Plaisance L, Knowlton N, Paulay G, Meyer C. 2009.** Reef-associated crustacean fauna: biodiversity estimates using semi-quantitative sampling and DNA barcoding. *Coral Reefs* **28**:977–986
DOI 10.1007/s00338-009-0543-3.

**Pollock FJ, Katz SM, Bourne DG, Willis BL. 2012.** *Cymo melanodactylus* crabs slow progression of white syndrome lesions on corals. *Coral Reefs* **32**:43–48 DOI 10.1007/s00338-012-0978-9.

**Pompanon F, Deagle BE, Symondson WOC, Brown DS, Jarman SN, Taberlet P. 2012.** Who is eating what: diet assessment using next generation sequencing. *Molecular Ecology* **21**:1931–1950
DOI 10.1111/j.1365-294X.2011.05403.x.

**Pratchett MS. 2005.** Dietary overlap among coral-feeding butterflyfishes (Chaetodontidae) at Lizard Island, northern Great Barrier Reef. *Marine Biology* **148**:373–382
DOI 10.1007/s00227-005-0084-4.
**Pratchett MS. 2007.** Dietary selection by coral-feeding butterflyfishes (Chaetodontidae) on the Great Barrier Reef, Australia. *Raffles Bulletin of Zoology* **14**:171–176.

**Pratchett MS, Berumen ML. 2008.** Interspecific variation in distributions and diets of coral reef butterflyfishes (Teleostei: Chaetodontidae). *Journal of Fish Biology* **73**:1730–1747 DOI 10.1111/j.1095-8649.2008.02062.x.

**Quince C, Lanzen A, Davenport RJ, Turnbaugh PJ. 2011.** Removing noise from pyrosequenced amplicons. *BMC Bioinformatics* **12**:38 DOI 10.1186/1471-2105-12-38.

**Randall JE. 1967.** Food habits of reef fishes of the West Indies. *Studies in Tropical Oceanography* **5**:665–847.

**Ranwez V, Harispe S, Delsuc F, Douzery EJP. 2011.** MACSE: multiple Alignment of Coding SEquences accounting for frameshifts and stop codons. *PLoS ONE* **6**:e22594 DOI 10.1371/journal.pone.0022594.

**Rawlinson KA, Gillis JA, Billings RE, Borneman EH. 2011.** Taxonomy and life history of the *Acropora*-eating flatworm *Amakusaplana acroporae nov. sp.* (Polycladida: Prothiostomidae). *Coral Reefs* **30**:693–705 DOI 10.1007/s00338-011-0745-3.

**Robertson DR, Lassig B. 1980.** Spatial distribution patterns and coexistence of a group of territorial damselfishes from the Great Barrier reef. *Bulletin of Marine Science* **30**:187–203.

**Ross ST. 1986.** Resource partitioning in fish assemblages—a review of field studies. *Copeia* **2**:352–388 DOI 10.2307/1444996.

**Rotjan R, Lewis S. 2008.** Impact of coral predators on tropical reefs. *Marine Ecology Progress Series* **367**:73–91 DOI 10.3354/meps07531.

**Rouzé H, Lecellier G, Langlade MJ, Planes S, Berteaux-Lecellier V. 2015.** Fringing reefs exposed to different levels of eutrophication and sedimentation can support similar benthic communities. *Marine Pollution Bulletin* **92**:212–221 DOI 10.1016/j.marpolbul.2014.12.016.

**Rouzé H, Lecellier G, Mills SC, Planes S, Berteaux-Lecellier V, Stewart H. 2014.** Juvenile *Trapezia spp.* crabs can increase juvenile host coral survival by protection from predation. *Marine Ecology Progress Series* **515**:151–159 DOI 10.3354/meps10970.

**Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, Sahl JW, Stres B, Thallinger GG, Van Horn DJ, Weber CF. 2009.** Introducing Mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. *Applied and Environmental Microbiology* **75**:7537–7541 DOI 10.1128/AEM.01541-09.

**Shehzad W, McCarthy TM, Pompanon F, Purevjav L, Coissac E, Riaz T, Taberlet P. 2012.** Prey preference of snow leopard (*Panthera uncia*) in South Gobi, Mongolia. *PLoS ONE* **7**:e32104 DOI 10.1371/journal.pone.0032104.

**Sheppard SK, Bell J, Sunderland KD, Fenlon J, Skervin D, Symondson WOC. 2005.** Detection of secondary predation by PCR analyses of the gut contents of invertebrate generalist predators. *Molecular Ecology* **14**:4461–4468 DOI 10.1111/j.1365-294X.2005.02742.x.

**Shima JS, Osenberg CW, Stier AC. 2010.** The vermetid gastropod *Dendropoma maximum* reduces coral growth and survival. *Biology Letters* **6**:815–818 DOI 10.1098/rsbl.2010.0291.

**Simon-Blecher N, Achituv Y. 1997.** Relationship between the coral pit crab *Cryptochirus coralliodytes* Heller and its host coral. *Journal of Experimental Marine Biology and Ecology* **215**:93–102 DOI 10.1016/S0022-0981(97)00002-6.

**Stella JS, Jones GP, Pratchett MS. 2010.** Variation in the structure of epifaunal invertebrate assemblages among coral hosts. *Coral Reefs* **29**:957–973 DOI 10.1007/s00338-010-0648-8.
**Stella JS, Pratchett MS, Hutchings PA, Jones GP. 2011.** Coral-associated invertebrates: diversity, ecological importance and vulnerability to disturbance. *Oceanography and Marine Biology: An Annual Review* **49**:43–104 DOI 10.1201/b11009-3.

**Stewart HL, Holbrook SJ, Schmitt RJ, Brooks AJ. 2006.** Symbiotic crabs maintain coral health by clearing sediments. *Coral Reefs* **25**:609–615 DOI 10.1007/s00338-006-0132-7.

**Stier AC, Gil MA, McKeon CS, Lemer S, Leray M, Mills SC, Osenberg CW. 2012.** Housekeeping mutualisms: do more symbionts facilitate host performance? *PLoS ONE* **7**:e32079 DOI 10.1371/journal.pone.0032079.

**Stier AC, Leray M. 2014.** Predators alter community organization of coral reef cryptofauna and reduce abundance of coral mutualists. *Coral Reefs* **33**:181–191 DOI 10.1007/s00338-013-1077-2.

**Stier AC, McKeon CS, Osenberg CW, Shima JS. 2010.** Guard crabs alleviate deleterious effects of vermetid snails on a branching coral. *Coral Reefs* **29**:1019–1022 DOI 10.1007/s00338-010-0663-9.

**Sussman M, Loya Y, Fine M, Rosenberg E. 2003.** The marine fireworm *Hermodice carunculata* is a winter reservoir and spring-summer vector for the coral-bleaching pathogen *Vibrio shiloi*. *Environmental Microbiology* **5**:250–255 DOI 10.1046/j.1462-2920.2003.00424.x.

**Symondson WOC. 2002.** Molecular identification of prey in predator diets. *Molecular Ecology* **11**:627–641 DOI 10.1046/j.1365-294X.2002.01471.x.

**Taberlet P, Coissac E, Pompanon F, Brochmann C, Willerslev E. 2012.** Towards next-generation biodiversity assessment using DNA metabarcoding. *Molecular Ecology* **21**:2045–2050 DOI 10.1111/j.1365-294X.2012.05470.x.

**Talbot FH, Russell BC, Anderson GRV. 1978.** Coral reef fish communities—unstable, high diversity systems. *Ecological Monographs* **48**:425–440 DOI 10.2307/2937241.

**Turner S. 1994.** The biology and population outbreaks of the corallivorous gastropod *Drupella* on Indo-Pacific reefs. *Oceanography and Marine Biology* **32**:461–530.

**Vestheim H, Jarman SN. 2008.** Blocking primers to enhance PCR amplification of rare sequences in mixed samples—a case study on prey DNA in Antarctic krill stomachs. *Frontiers in Zoology* **5**:12 DOI 10.1186/1742-9994-5-12.

**Waldner RE, Robertson DR. 1980.** Patterns of habitat partitioning by eight species of territorial Carribbean damselfishes (Pisces: Pomacentridae). *Bulletin of Marine Science* **30**:171–186.

**Williams D, Miller M. 2005.** Coral disease outbreak: pattern, prevalence and transmission in *Acropora cervicornis*. *Marine Ecology Progress Series* **301**:119–128 DOI 10.3354/meps301119.

**Wright SJ. 2005.** Tropical forests in a changing environment. *Trends in Ecology & Evolution* **20**:553–560 DOI 10.1016/j.tree.2005.07.009.