Use of 16S ribosomal DNA for delineation of marine bacterioplankton species

All of the marine bacterioplankton-derived 16S ribosomal DNA sequences previously deposited in GenBank were reanalyzed to determine the number of bacterial species in the oceanic surface waters. These sequences have been entered into the database since 1990. The rate of new additions reached a peak in 1999 and subsequently leveled off, suggesting that much of the marine microbial species richness has been sampled. When the GenBank sequences were dereplicated by using 97% similarity as a cutoff, 1,117 unique ribotypes were found. Of the unique sequences, 609 came from uncultured environmental clones and 508 came from cultured bacteria. We conclude that the apparent bacterioplankton species richness is relatively low.     

Status of the microbial census.

Over the past 20 years, more than 78,000 16S rRNA gene sequences have been deposited in GenBank and the Ribosomal Database Project, making the 16S rRNA gene the most widely studied gene for reconstructing bacterial phylogeny. While there is a general appreciation that these sequences are largely unique and derived from diverse species of bacteria, there has not been a quantitative attempt to describe the extent of sequencing efforts to date. We constructed rarefaction curves for each bacterial phylum and for the entire bacterial domain to assess the current state of sampling and the relative taxonomic richness of each phylum. This analysis quantifies the general sense among microbiologists that we are a long way from a complete census of the bacteria on Earth. Moreover, the analysis indicates that current sampling strategies might not be the most effective ones to describe novel diversity because there remain numerous phyla that are globally distributed yet poorly sampled. Based on the current level of sampling, it is not possible to estimate the total number of bacterial species on Earth, but the minimum species richness is 35,498. Considering previous global species richness estimates of 10(7) to 10(9), we are certain that this estimate will increase with additional sequencing efforts. The data support previous calls for extensive surveys of multiple chemically disparate environments and of specific phylogenetic groups to advance the census most rapidly.     

Status of the Microbial Census

Over the past 20 years, more than 78,000 16S rRNA gene sequences have been deposited in GenBank and the Ribosomal Database Project, making the 16S rRNA gene the most widely studied gene for reconstructing bacterial phylogeny. While there is a general appreciation that these sequences are largely unique and derived from diverse species of bacteria, there has not been a quantitative attempt to describe the extent of sequencing efforts to date. We constructed rarefaction curves for each bacterial phylum and for the entire bacterial domain to assess the current state of sampling and the relative taxonomic richness of each phylum. This analysis quantifies the general sense among microbiologists that we are a long way from a complete census of the bacteria on Earth. Moreover, the analysis indicates that current sampling strategies might not be the most effective ones to describe novel diversity because there remain numerous phyla that are globally distributed yet poorly sampled. Based on the current level of sampling, it is not possible to estimate the total number of bacterial species on Earth, but the minimum species richness is 35,498. Considering previous global species richness estimates of 10⁷ to 10⁹, we are certain that this estimate will increase with additional sequencing efforts. The data support previous calls for extensive surveys of multiple chemically disparate environments and of specific phylogenetic groups to advance the census most rapidly.     

Diversity of basidiomycetes in michigan agricultural soils

We analyzed the communities of soil basidiomycetes in agroecosystems that differ in tillage history at the Kellogg Biological Station Long-Term Ecological Research site near Battle Creek, Michigan. The approach combined soil DNA extraction through a bead-beating method modified to increase recovery of fungal DNA, PCR amplification with basidiomycete-specific primers, cloning and restriction fragment length polymorphism screening of mixed PCR products, and sequencing of unique clones. Much greater diversity was detected than was anticipated in this habitat on the basis of culture-based methods or surveys of fruiting bodies. With "species" defined as organisms yielding PCR products with > or =99% identity in the 5' 650 bases of the nuclear large-subunit ribosomal DNA, 241 "species" were detected among 409 unique basidiomycete sequences recovered. Almost all major clades of basidiomycetes from basidiomycetous yeasts and other heterobasidiomycetes through polypores and euagarics (gilled mushrooms and relatives) were represented, with a majority from the latter clade. Only 24 of 241 "species" had 99% or greater sequence similarity to named reference sequences in GenBank, and several clades with multiple "species" could not be identified at the genus level by phylogenetic comparisons with named sequences. The total estimated "species" richness for this 11.2-ha site was 367 "species" of basidiomycetes. Since >99% of the study area has not been sampled, the accuracy of our diversity estimate is uncertain. Replication in time and space is required to detect additional diversity and the underlying community structure.     

Identifying gastropod spawn from DNA barcodes: possible but not yet practicable

Identifying life stages of species with complex life histories is problematic as species are often only known and/or described from a single stage. DNA barcoding has been touted as an important tool for linking life-history stages of the same species. To test the current efficacy of DNA barcodes for identifying unknown mollusk life stages, 24 marine gastropod egg capsules were collected off the Philippines in deep water and sequenced for partial fragments of the COI, 16S and 12S mitochondrial genes. Two egg capsules of known shallow-water Mediterranean species were used to calibrate the method. These sequences were compared to those available in GenBank and the Barcode of Life Database (BOLD). Using COI sequences alone, only a single Mediterranean egg capsule was identified to species, and a single Philippine egg capsule was identified tentatively to genus; all other COI sequences recovered matches between 76% and 90% with sequences from BOLD and GenBank. Similarity-based identification using all three markers confirmed the Mediterranean specimens' identifications. A phylogenetic approach was also implemented to confirm similarity-based identifications and provide a higher-taxonomic identification when species-level identifications were not possible. Comparison of available GenBank sequences to the diversity curve of a well-sampled coral reef habitat in New Caledonia highlights the poor taxonomic coverage achieved at present in existing genetic databases, emphasizing the need to develop DNA barcoding projects for megadiverse and often taxonomically challenging groups such as mollusks, to fully realize its potential as an identification and discovery tool.     

Taxonomic and geographic bias in the genetic study of helminth parasites

The use of genetic information is now fundamental in parasite taxonomy and systematics, for resolving parasite phylogenies, discovering cryptic species, and elucidating patterns of gene flow among parasite populations. The accumulation of available gene sequences per geographical area or per parasite taxonomic group is likely proportional to species richness, but not without some biases. Certain areas and certain taxonomic groups receive more research effort than others, possibly causing a deficit in the relative number of parasite species being characterized genetically in some areas or taxonomic groups. Here, we use data on the number of parasite records per country or helminth family from the London Natural History Museum host-parasite database, and matching data on the number of gene sequences available from the National Center for Biotechnology Information (NCBI) GenBank database, to determine how available gene sequences scale with species richness across countries or parasitic helminth families. Our quantitative analysis identified countries/regions of the world and helminth families that have received the most effort in genetic research. More importantly, it allowed us to generate lists (based on residuals from the statistical model) of the 20 countries/regions and the 20 helminth families with the largest deficit in available gene sequences relative to their helminth species richness. We propose these lists as useful guides toward future allocation of effort to maximise advances in parasite biodiscovery, systematics and population structure.     

The ‘Expansion–Contraction’ model of Pleistocene biogeography: rocky shores suffer a sea change?

Approximately 20 000 years ago the last glacial maximum (LGM) radically altered the distributions of many Northern Hemisphere terrestrial organisms. Fewer studies describing the biogeographic responses of marine species to the LGM have been conducted, but existing genetic data from coastal marine species indicate that fewer taxa show clear signatures of post-LGM recolonization. We have assembled a mitochondrial DNA (mtDNA) data set for 14 co-distributed northeastern Pacific rocky-shore species from four phyla by combining new sequences from ten species with previously published sequences from eight species. Nuclear sequences from four species were retrieved from GenBank, plus we gathered new elongation factor 1-α sequences from the barnacle Balanus glandula . Results from demographic analyses of mtDNA for five (36%) species ( Evasterias troschelii, Pisaster ochraceus, Littorina sitkana, L. scutulata, Xiphister mucosus ) were consistent with large population expansions occurring near the LGM, a pattern expected if these species recently recolonized the region. However, seven (50%) species ( Mytilus trossulus, M. californianus, B. glandula, S. cariosus, Patiria miniata, Katharina tunicata , X. atropurpureus ) exhibited histories consistent with long-term stability in effective population size, a pattern indicative of regional persistence during the LGM. Two species of Nucella with significant mtDNA genetic structure showed spatially variable demographic histories. Multilocus analyses for five species were largely consistent with mtDNA: the majority of multilocus interpopulation divergence times significantly exceeded the LGM. Our results indicate that the LGM did not extirpate the majority of species in the northeastern Pacific; instead, regional persistence during the LGM appears a common biogeographic history for rocky-shore organisms in this region.     

Barcoding marine nematodes: an improved set of nematode 18S rRNA primers to overcome eukaryotic co-interference

Nematodes form an important component of many benthic marine ecosystems and DNA barcoding approaches could provide an insight into nematode community composition from different environments globally. We have amplified nematode 18S rRNA sequences using standard nematode18S rRNA primers from environmental DNA extracted from intertidal sediment collected from New Jersey coast, USA to test whether the published marine nematode 18S rRNA sequences from GenBank and EMBL databases can effectively assign unknown nematode sequences into genus or species level. Most of the sequenced clones showed some degree of identities with published marine nematode 18S rRNA sequences. However, relatively very few of the sequences could be assigned even to genus level based on sequence assignment rule. In addition, other eukaryotic 18S rRNA sequences were found to be co-amplified with commonly used nematode 18S rRNA primers. We found that the majority of the current nematode 18S rRNA primers will co-amplify other eukaryotes if environmental DNA is the target template. We therefore designed a new set of nematode 18S rRNA primers and evaluated them using environmental DNA in intertidal sediment from the New Jersey coast. In total, 40 clones were screened and subsequently sequenced and all the sequences showed varying degree of identities with published nematode 18S rRNA sequences from GenBank and EMBL databases, and no obvious eukaryotic co-amplicons were detected with new primers. Only 13 out of 40 clones amplified with the new primer set showed 100% identity to published Daptonema and Metachromadora 18S rRNA sequences. The current molecular databases for nematodes are dominated by sequences from NW Europe and need to be more extensively populated with new full length 18S rRNA nematode sequences collected from different biogeographic locations. The new primers developed in this study, in combination with an updated nematode 18S rRNA sequence database, would help us to better investigate and understand the diversity and community composition of free-living marine nematodes based on DNA barcoding approaches during biodiversity or biomonitoring surveys on a global-scale.     

Diversity of Siphonaria Sowerby I, 1823 (Gastropoda,Siphonariidae) in the Seychelles Bank and beyond

Molecular phylogenetic studies have shown that the characters of the reduced shell of the false limpets of the genus Siphonaria Sowerby I, 1823 are highly variable and often insufficient for species delimitation. The taxonomy and distribution of Siphonaria in the Indian Ocean are poorly known. We sampled Siphonaria in the Seychelles Bank to check the occurrence of recorded species using DNA sequences and to study the paths through which Siphonaria species have colonised the Seychelles Bank by reconstructing their phylogenetic relationships. Analyses of a dataset comprising 16 S rRNA gene sequences of 33 specimens from the Seychelles Bank and 300 additional Siphonaria sequences from other regions from GenBank with various methods for species delimitation resulted in 19–102 primary species hypotheses. Assemble Species by Automatic Partitioning provided a conservative estimate of the species number (42) in which several indisputable species were lumped. The results of Automatic Barcode Gap Discovery depended strongly on the assumed prior maximum intraspecific divergence, whereas the tree-based methods Generalised Mixed Yule Coalescent and Poisson Tree Processes resulted in high overestimates. The specimens from the Seychelles Bank represent three clades, belonging to the Siphonaria ‘atra’ group, the Siphonaria ‘normalis’ group and a possibly undescribed species recorded previously only from Hainan. At least two of the three species recorded from the Seychelles Bank came from the east, i.e., from the Coral Triangle in the Indo-Australian Archipelago, the region with the highest marine biodiversity worldwide. A major transport mechanism across the Indian Ocean was probably the South Equatorial Current.     

DNA Barcoding for Species Assignment: The Case of Mediterranean Marine Fishes

Background

DNA barcoding enhances the prospects for species-level identifications globally using a standardized and authenticated DNA-based approach. Reference libraries comprising validated DNA barcodes (COI) constitute robust datasets for testing query sequences, providing considerable utility to identify marine fish and other organisms. Here we test the feasibility of using DNA barcoding to assign species to tissue samples from fish collected in the central Mediterranean Sea, a major contributor to the European marine ichthyofaunal diversity.

Methodology/Principal Findings

A dataset of 1278 DNA barcodes, representing 218 marine fish species, was used to test the utility of DNA barcodes to assign species from query sequences. We tested query sequences against 1) a reference library of ranked DNA barcodes from the neighbouring North East Atlantic, and 2) the public databases BOLD and GenBank. In the first case, a reference library comprising DNA barcodes with reliability grades for 146 fish species was used as diagnostic dataset to screen 486 query DNA sequences from fish specimens collected in the central basin of the Mediterranean Sea. Of all query sequences suitable for comparisons 98% were unambiguously confirmed through complete match with reference DNA barcodes. In the second case, it was possible to assign species to 83% (BOLD-IDS) and 72% (GenBank) of the sequences from the Mediterranean. Relatively high intraspecific genetic distances were found in 7 species (2.2%–18.74%), most of them of high commercial relevance, suggesting possible cryptic species.

Conclusion/Significance

We emphasize the discriminatory power of COI barcodes and their application to cases requiring species level resolution starting from query sequences. Results highlight the value of public reference libraries of reliability grade-annotated DNA barcodes, to identify species from different geographical origins. The ability to assign species with high precision from DNA samples of disparate quality and origin has major utility in several fields, from fisheries and conservation programs to control of fish products authenticity.     

Use of SSCP to improve the efficiency of microsatellite identification from microsatellite‐enriched libraries

An inefficient aspect of marker identification from microsatellite‐enriched libraries is the proportion of clones with identical sequences. This can substantially increase the number of clones that need to be sequenced in order to identify a sufficient number of microsatellite loci. We propose the use of single‐stranded conformation polymorphism (SSCP) analysis to identify unique clones prior to sequencing. We used this approach prior to sequencing from microsatellite‐enriched libraries for three marine invertebrate species and were able to obtain a given number of unique clone sequences for only 28% of the sequencing effort that would have been required without SSCP screening.     

Large-scale species-richness gradients in the Atlantic Ocean

The increase in species richness from the poles to the Equator has been observed in numerous terrestrial and aquatic taxa. A number of different hypotheses have been put forward as explanations for this trend, e.g. area and energy availability. However, whether these hypotheses apply to large spatial scales in marine environments remains unclear. The present study shows a clear latitudinal gradient from high to low latitude (from 80 degrees N to 70 degrees S) in marine species richness for 6643 species (fishes and invertebrates) in 10 different taxa dwelling in benthic and pelagic habitats on both sides of the Atlantic. The patterns in benthic taxa are strongly influenced by coastal hydrographic processes, with marked peaks and troughs, and consequently the gradients are not symmetric along both Atlantic sides. Pelagic taxa show a plateau-shaped distribution and the influence from coastal events on gradients could not be demonstrated. The relationships between species richness and different environmental factors indicate that area size does not explain the latitudinal pattern in benthic species richness on a large spatial scale. Sea-surface temperature (positive relationship) is the best predictor of this pattern for benthic species, and nitrate concentration (negative relationship) is the best predictor for pelagic species. The results call into question the existence of a single primary cause that would explain the pattern in marine species richness on a large spatial scale.     

Fungal community analysis by large-scale sequencing of environmental samples

Fungi are an important and diverse component of soil communities, but these communities have proven difficult to study in conventional biotic surveys. We evaluated soil fungal diversity at two sites in a temperate forest using direct isolation of small-subunit and internal transcribed spacer (ITS) rRNA genes by PCR and high-throughput sequencing of cloned fragments. We identified 412 sequence types from 863 fungal ITS sequences, as well as 112 ITS sequences from other eukaryotic microorganisms. Equal proportions of Basidiomycota and Ascomycota sequences were present in both the ITS and small-subunit libraries, while members of other fungal phyla were recovered at much lower frequencies. Many sequences closely matched sequences from mycorrhizal, plant-pathogenic, and saprophytic fungi. Compositional differences were observed among samples from different soil depths, with mycorrhizal species predominating deeper in the soil profile and saprophytic species predominating in the litter layer. Richness was consistently lowest in the deepest soil horizon samples. Comparable levels of fungal richness have been observed following traditional specimen-based collecting and culturing surveys, but only after much more extensive sampling. The high rate at which new sequence types were recovered even after sampling 863 fungal ITS sequences and the dominance of fungi in our libraries relative to other eukaryotes suggest that the abundance and diversity of fungi in forest soils may be much higher than previously hypothesized. All sequences were deposited in GenBank, with accession numbers AY 969316 to AY 970290 for the ITS sequences and AY 969135 to AY 969315 for the SSU sequences.     

Characterization of Cryptocaryon irritans,a parasite isolated from marine fishes in Taiwan

The ciliated protozoan parasite Cryptocaryon irritans infecting marine fishes in Taiwan is described. Developmental characteristics and sequences of the ribosomal DNA regions such as part of 18 S, the entire first internal transcribed spacer, and part of 5.8 S of various Taiwan isolates of C. irritans were investigated. A total of 5 isolates was obtained from different fish-host species and localities, the majority from cultured fish species. C. irritans from Taiwan is able to shift its developmental characteristics, i.e. from non-adherent to adherent tomonts, from individualistic to aggregate-forming tomonts, from infection of the gills only to infection of the gills and body. Thus, it is not possible to classify strains of C. irritans on the basis of these parameters. Premature tomonts that developed from dead fishes were able to produce theronts that could infect fish host. Isolates from Pingtung and the USA had identical nucleotide sequences while an isolate from Malaysia was identical to an Israel isolate. Percentage variation among pairs of Taiwan isolates showed a higher degree of variation than isolate sequences listed in GenBank. Sequence analysis revealed highly aberrant isolates in Taiwan, and a phylogenetic tree distinguished a marine and a low-salinity variant. C. irritans from marine fishes in Taiwan, therefore, display some characteristics not previously reported. Since manipulation of salinity in brackishwater ponds and marine cage sites is not feasible, there is a need to develop new strategies for the control and prevention of cryptocaryoniasis.     

Analysis of messages expressed by Echinostoma paraensei miracidia and sporocysts, obtained by random EST sequencing

A lambdaZAP Express cDNA library was constructed with mRNA obtained from immature miracidia within eggs, hatched miracidia, and sporocysts of Echinostoma paraensei. This cDNA library was amplified and 213 expressed sequence tag (EST) sequences (averaging 466 nucleotides in length) were obtained. The mean percentage of unresolved bases within the EST sequences was 0.4%, ranging from 0 to 4.6%. The 213 ESTs represent 151 unique messages. BLAST (version 2.0.8) analysis disclosed that 64 unique E. paraensei messages (42.4%) had significant similarities (BLAST score < or =e-5), at deduced amino acid or nucleotide levels, with known sequences in the nonredundant GenBank databases or the dbEST database (NCBI). The remainder, 57.6% of the unique EST-encoded messages, scored nonsignificant hits. Most of the E. paraensei messages that could be assigned a cellular role based on sequence similarities were involved in gene/protein expression. Several ESTs scored highest similarities with sequences obtained from trematode species. A total of 22,560 nucleotides present in open reading frames from ESTs that aligned with known sequences was used to determine codon usage for E. paraensei. Analysis of a subset of eight ESTs that contained full-length open reading frames did not reveal a bias in codon usage. Also, EST sequences were found to contain 3' untranslated regions with an average length of 69.9 +/- 88.4 nucleotides (n = 46). The EST sequences were submitted to GenBank/dbEST, adding to the 51 available Echinostoma-derived sequences, to provide reference information for both phylogenetic analysis and study of general trematode biology.     

Isolation and diversity of planctomycetes from the sponge Niphates sp., seawater, and sediment of Moreton Bay, Australia

Planctomycetes are ubiquitous in marine environment and were reported to occur in association with multicellular eukaryotic organisms such as marine macroalgae and invertebrates. Here, we investigate planctomycetes associated with the marine sponge Niphates sp. from the sub-tropical Australian coast by assessing their diversity using culture-dependent and -independent approaches based on the 16S rRNA gene. The culture-dependent approach resulted in the isolation of a large collection of diverse planctomycetes including some novel lineages of Planctomycetes from the sponge as well as sediment and seawater of Moreton Bay where this sponge occurs. The characterization of these novel planctomycetes revealed that cells of one unique strain do not possess condensed nucleoids, a phenotype distinct from other planctomycetes. In addition, a culture-independent clone library approach identified unique planctomycete 16S rRNA gene sequences closely related to other sponge-derived sequences. The analysis of tissue of the sponge Niphates sp. showed that the mesohyl of the sponge is almost devoid of microbial cells, indicating this species is in the group of ‘low microbial abundant’ (LMA) sponges. The unique planctomycete 16S rRNA gene sequences identified in this study were phylogenetically closely related to sequences from LMA sponges in other published studies. This study has revealed new insights into the diversity of planctomycetes in the marine environment and the association of planctomycetes with marine sponges.     

Molecular detection of fungal communities in the Hawaiian marine sponges Suberites zeteki and Mycale armata

Symbiotic microbes play a variety of fundamental roles in the health and habitat ranges of their hosts. While prokaryotes in marine sponges have been broadly characterized, the diversity of sponge-inhabiting fungi has barely been explored using molecular approaches. Fungi are an important component of many marine and terrestrial ecosystems, and they may be an ecologically significant group in sponge-microbe interactions. This study tested the feasibility of using existing fungal primers for molecular analysis of sponge-associated fungal communities. None of the eight selected primer pairs yielded satisfactory results in fungal rRNA gene or internal transcribed spacer (ITS) clone library constructions. However, 3 of 10 denaturing gradient gel electrophoresis (DGGE) primer sets, which were designed to preferentially amplify fungal rRNA gene or ITS regions from terrestrial environmental samples, were successfully amplified from fungal targets in marine sponges. DGGE analysis indicated that fungal communities differ among different sponge species (Suberites zeteki and Mycale armata) and also vary between sponges and seawater. Sequence analysis of DGGE bands identified 23 and 21 fungal species from each of the two sponge species S. zeteki and M. armata, respectively. These species were representatives of 11 taxonomic orders and belonged to the phyla of Ascomycota (seven orders) and Basidiomycota (four orders). Five of these taxonomic orders (Malasseziales, Corticiales, Polyporales, Agaricales, and Dothideomycetes et Chaetothyriomcetes incertae sedis) have now been identified for the first time in marine sponges. Seven and six fungal species from S. zeteki and M. armata, respectively, are potentially new species because of their low sequence identity (< or =98%) with their references in GenBank. Phylogenetic analysis indicated sponge-derived sequences were clustered into "marine fungus clades" with those from other marine habitats. This is the first report of molecular analysis of fungal communities in marine sponges, adding depth and dimension to our understanding of sponge-associated microbial communities.     

Analysis of bacterial communities in sponges and coral inhabiting Red Sea,using barcoded 454 pyrosequencing

Microbial communities are linked with marine sponge are diverse in their structure and function. Our understanding of the sponge-associated microbial diversity is limited especially from Red Sea in Saudi Arabia where few species of sponges have been studied. Here we used pyrosequencing to study two marine sponges and coral species sampled from Obhur region from Red sea in Jeddah. A total of 168 operational taxonomic units (OTUs) were identified from Haliclona caerulea, Stylissa carteri and Rhytisma fulvum. Taxonomic identification of tag sequences of 16S ribosomal RNA revealed 6 different bacterial phyla and 9 different classes. A proportion of unclassified reads were was also observed in sponges and coral sample. We found diverse bacterial communities associated with two sponges and a coral sample. Diversity and richness estimates based on OUTs revealed that sponge H. caerulea had significantly high bacterial diversity. The identified OTUs showed unique clustering in three sponge samples as revealed by Principal coordinate analysis (PCoA). Proteobacteria (88–95%) was dominant phyla alonwith Bacteroidetes, Planctomycetes, Cyanobacteria, Firmicutes and Nitrospirae. Seventeen different genera were identified where genus Pseudoalteromonas was dominant in all three samples. This is first study to assess bacterial communities of sponge and coral sample that have never been studied before to unravel their microbial communities using 454-pyrosequencing method.     

Evidence of artisanal fishing impacts and depth refuge in assemblages of Fijian reef fish

Protection from fishing generally results in an increase in the abundance and biomass of species targeted by fisheries within marine reserve boundaries. Natural refuges such as depth may also protect such species, yet few studies in the Indo Pacific have investigated the effects of depth concomitant with marine reserves. We studied the effects of artisanal fishing and depth on reef fish assemblages in the Kubulau District of Vanua Levu Island, Fiji, using baited remote underwater stereo-video systems. Video samples were collected from shallow (5–8 m) and deep (25–30 m) sites inside and outside of a large old marine reserve (60.6 km², 13 years old) and a small new marine reserve (4.25 km², 4 years old). Species richness tended to be greater in the shallow waters of the large old reserve when compared to fished areas. In the deeper waters, species richness appeared to be comparable. The difference in shallow waters was driven by species targeted by fisheries, indicative of a depth refuge effect. In contrast, differences in the abundance composition of the fish assemblage existed between protected and fished areas for deep sites, but not shallow. Fish species targeted by local fisheries were 89% more abundant inside the large old reserve than surrounding fished areas, while non-targeted species were comparable. We observed no difference in the species richness or abundance of species targeted by fisheries inside and outside of the small new reserve. This study suggests that artisanal fishing impacts on the abundance and species richness of coral reef fish assemblages and effects of protection are more apparent with large reserves that have been established for a long period of time. Observed effects of protection also vary with depth, highlighting the importance of explicitly incorporating multiple depth strata in studies of marine reserves.