Gene flow throughout the evolutionary history of a colour polymorphic and generalist clownfish.

Even seemingly homogeneous on the surface, the oceans display high environmental heterogeneity across space and time. Indeed, different soft barriers structure the marine environment, which offers an appealing opportunity to study various evolutionary processes such as population differentiation and speciation. Here, we focus on Amphiprion clarkii (Actinopterygii; Perciformes), the most widespread of clownfishes that exhibits the highest colour polymorphism. Clownfishes can only disperse during a short pelagic larval phase before their sedentary adult lifestyle, which might limit connectivity among populations, thus facilitating speciation events. Consequently, the taxonomic status of A. clarkii has been under debate. We used whole-genome resequencing data of 67 A. clarkii specimens spread across the Indian and Pacific Oceans to characterize the species' population structure, demographic history and colour polymorphism. We found that A. clarkii spread from the Indo-Pacific Ocean to the Pacific and Indian Oceans following a stepping-stone dispersal and that gene flow was pervasive throughout its demographic history. Interestingly, colour patterns differed noticeably among the Indonesian populations and the two populations at the extreme of the sampling distribution (i.e. Maldives and New Caledonia), which exhibited more comparable colour patterns despite their geographic and genetic distances. Our study emphasizes how whole-genome studies can uncover the intricate evolutionary past of wide-ranging species with diverse phenotypes, shedding light on the complex nature of the species concept paradigm.

Diverse interactions and ecosystem engineering can stabilize community assembly.

The complexity of an ecological community can be distilled into a network, where diverse interactions connect species in a web of dependencies. Species interact directly with each other and indirectly through environmental effects, however to our knowledge the role of these ecosystem engineers has not been considered in ecological network models. Here we explore the dynamics of ecosystem assembly, where species colonization and extinction depends on the constraints imposed by trophic, service, and engineering dependencies. We show that our assembly model reproduces many key features of ecological systems, such as the role of generalists during assembly, realistic maximum trophic levels, and increased nestedness with mutualistic interactions. We find that ecosystem engineering has large and nonlinear effects on extinction rates. While small numbers of engineers reduce stability by increasing primary extinctions, larger numbers of engineers increase stability by reducing primary extinctions and extinction cascade magnitude. Our results suggest that ecological engineers may enhance community diversity while increasing persistence by facilitating colonization and limiting competitive exclusion.

Environmental DNA survey captures patterns of fish and invertebrate diversity across a tropical seascape.

Accurate, rapid, and comprehensive biodiversity assessments are critical for investigating ecological processes and supporting conservation efforts. Environmental DNA (eDNA) surveys show promise as a way to effectively characterize fine-scale patterns of community composition. We tested whether a single PCR survey of eDNA in seawater using a broad metazoan primer could identify differences in community composition between five adjacent habitats at 19 sites across a tropical Caribbean bay in Panama. We paired this effort with visual fish surveys to compare methods for a conspicuous taxonomic group. eDNA revealed a tremendous diversity of animals (8,586 operational taxonomic units), including many small taxa that would be undetected in traditional in situ surveys. Fish comprised only 0.07% of the taxa detected by a broad COI primer, yet included 43 species not observed in the visual survey. eDNA revealed significant differences in fish and invertebrate community composition across adjacent habitats and areas of the bay driven in part by taxa known to be habitat-specialists or tolerant to wave action. Our results demonstrate the ability of broad eDNA surveys to identify biodiversity patterns in the ocean.

The indirect paths to cascading effects of extinctions in mutualistic networks.

Biodiversity loss is a hallmark of our times, but predicting its consequences is challenging. Ecological interactions form complex networks with multiple direct and indirect paths through which the impacts of an extinction may propagate. Here we show that accounting for these multiple paths connecting species is necessary to predict how extinctions affect the integrity of ecological networks. Using an approach initially developed for the study of information flow, we estimate indirect effects in plant-pollinator networks and find that even those species with several direct interactions may have much of their influence over others through long indirect paths. Next, we perform extinction simulations in those networks and show that although traditional connectivity metrics fail in the prediction of coextinction patterns, accounting for indirect interaction paths allows predicting species' vulnerability to the cascading effects of an extinction event. Embracing the structural complexity of ecological systems contributes towards a more predictive ecology, which is of paramount importance amid the current biodiversity crisis.

Revealing biases in the sampling of ecological interaction networks.

The structure of ecological interactions is commonly understood through analyses of interaction networks. However, these analyses may be sensitive to sampling biases with respect to both the interactors (the nodes of the network) and interactions (the links between nodes), because the detectability of species and their interactions is highly heterogeneous. These ecological and statistical issues directly affect ecologists' abilities to accurately construct ecological networks. However, statistical biases introduced by sampling are difficult to quantify in the absence of full knowledge of the underlying ecological network's structure. To explore properties of large-scale ecological networks, we developed the software EcoNetGen, which constructs and samples networks with predetermined topologies. These networks may represent a wide variety of communities that vary in size and types of ecological interactions. We sampled these networks with different mathematical sampling designs that correspond to methods used in field observations. The observed networks generated by each sampling process were then analyzed with respect to the number of components, size of components and other network metrics. We show that the sampling effort needed to estimate underlying network properties depends strongly both on the sampling design and on the underlying network topology. In particular, networks with random or scale-free modules require more complete sampling to reveal their structure, compared to networks whose modules are nested or bipartite. Overall, modules with nested structure were the easiest to detect, regardless of the sampling design used. Sampling a network starting with any species that had a high degree (e.g., abundant generalist species) was consistently found to be the most accurate strategy to estimate network structure. Because high-degree species tend to be generalists, abundant in natural communities relative to specialists, and connected to each other, sampling by degree may therefore be common but unintentional in empirical sampling of networks. Conversely, sampling according to module (representing different interaction types or taxa) results in a rather complete view of certain modules, but fails to provide a complete picture of the underlying network. To reduce biases introduced by sampling methods, we recommend that these findings be incorporated into field design considerations for projects aiming to characterize large species interaction networks.

Spatial distribution of environmental DNA in a nearshore marine habitat.

In the face of increasing threats to biodiversity, the advancement of methods for surveying biological communities is a major priority for ecologists. Recent advances in molecular biological technologies have made it possible to detect and sequence DNA from environmental samples (environmental DNA or eDNA); however, eDNA techniques have not yet seen widespread adoption as a routine method for biological surveillance primarily due to gaps in our understanding of the dynamics of eDNA in space and time. In order to identify the effective spatial scale of this approach in a dynamic marine environment, we collected marine surface water samples from transects ranging from the intertidal zone to four kilometers from shore. Using PCR primers that target a diverse assemblage of metazoans, we amplified a region of mitochondrial 16S rDNA from the samples and sequenced the products on an Illumina platform in order to detect communities and quantify their spatial patterns using a variety of statistical tools. We find evidence for multiple, discrete eDNA communities in this habitat, and show that these communities decrease in similarity as they become further apart. Offshore communities tend to be richer but less even than those inshore, though diversity was not spatially autocorrelated. Taxon-specific relative abundance coincided with our expectations of spatial distribution in taxa lacking a microscopic, pelagic life-history stage, though most of the taxa detected do not meet these criteria. Finally, we use carefully replicated laboratory procedures to show that laboratory treatments were remarkably similar in most cases, while allowing us to detect a faulty replicate, emphasizing the importance of replication to metabarcoding studies. While there is much work to be done before eDNA techniques can be confidently deployed as a standard method for ecological monitoring, this study serves as a first analysis of diversity at the fine spatial scales relevant to marine ecologists and confirms the promise of eDNA in dynamic environments.

Genetic signatures of ecological diversity along an urbanization gradient.

Despite decades of work in environmental science and ecology, estimating human influences on ecosystems remains challenging. This is partly due to complex chains of causation among ecosystem elements, exacerbated by the difficulty of collecting biological data at sufficient spatial, temporal, and taxonomic scales. Here, we demonstrate the utility of environmental DNA (eDNA) for quantifying associations between human land use and changes in an adjacent ecosystem. We analyze metazoan eDNA sequences from water sampled in nearshore marine eelgrass communities and assess the relationship between these ecological communities and the degree of urbanization in the surrounding watershed. Counter to conventional wisdom, we find strongly increasing richness and decreasing beta diversity with greater urbanization, and similar trends in the diversity of life histories with urbanization. We also find evidence that urbanization influences nearshore communities at local (hundreds of meters) rather than regional (tens of km) scales. Given that different survey methods sample different components of an ecosystem, we then discuss the advantages of eDNA-which we use here to detect hundreds of taxa simultaneously-as a complement to traditional ecological sampling, particularly in the context of broad ecological assessments where exhaustive manual sampling is impractical. Genetic data are a powerful means of uncovering human-ecosystem interactions that might otherwise remain hidden; nevertheless, no sampling method reveals the whole of a biological community.

Indexed PCR Primers Induce Template-Specific Bias in Large-Scale DNA Sequencing Studies.

Massively parallel sequencing is rapidly emerging as an efficient way to quantify biodiversity at all levels, from genetic variation and expression to ecological community assemblage. However, the number of reads produced per sequencing run far exceeds the number required per sample for many applications, compelling researchers to sequence multiple samples per run in order to maximize efficiency. For studies that include a PCR step, this can be accomplished using primers that include an index sequence allowing sample origin to be determined after sequencing. The use of indexed primers assumes they behave no differently than standard primers; however, we found that indexed primers cause substantial template sequence-specific bias, resulting in radically different profiles of the same environmental sample. Likely the outcome of differential amplification efficiency due to primer-template mismatch, two indexed primer sets spuriously change the inferred sequence abundance from the same DNA extraction by up to 77.1%. We demonstrate that a double PCR approach alleviates these effects in applications where indexed primers are necessary.

Assessing vertebrate biodiversity in a kelp forest ecosystem using environmental DNA.

Preserving biodiversity is a global challenge requiring data on species' distribution and abundance over large geographic and temporal scales. However, traditional methods to survey mobile species' distribution and abundance in marine environments are often inefficient, environmentally destructive, or resource-intensive. Metabarcoding of environmental DNA (eDNA) offers a new means to assess biodiversity and on much larger scales, but adoption of this approach for surveying whole animal communities in large, dynamic aquatic systems has been slowed by significant unknowns surrounding error rates of detection and relevant spatial resolution of eDNA surveys. Here, we report the results of a 2.5 km eDNA transect surveying the vertebrate fauna present along a gradation of diverse marine habitats associated with a kelp forest ecosystem. Using PCR primers that target the mitochondrial 12S rRNA gene of marine fishes and mammals, we generated eDNA sequence data and compared it to simultaneous visual dive surveys. We find spatial concordance between individual species' eDNA and visual survey trends, and that eDNA is able to distinguish vertebrate community assemblages from habitats separated by as little as ~60 m. eDNA reliably detected vertebrates with low false-negative error rates (1/12 taxa) when compared to the surveys, and revealed cryptic species known to occupy the habitats but overlooked by visual methods. This study also presents an explicit accounting of false negatives and positives in metabarcoding data, which illustrate the influence of gene marker selection, replication, contamination, biases impacting eDNA count data and ecology of target species on eDNA detection rates in an open ecosystem.

Variation in partner benefits in a shrimp-sea anemone symbiosis.

Symbiotic interactions, where two species occur in close physical proximity for the majority of the participants' lifespans, may constrain the fitness of one or both of the participants. Host choice could result in lineage divergence in symbionts if fitness benefits vary across the interaction with hosts. Symbiotic interactions are common in the marine environment, particularly in the most diverse marine ecosystems: coral reefs. However, the variation in symbiotic interactions that may drive diversification is poorly understood in marine systems. We measured the fecundity of the symbiotic shrimp Periclimenes yucatanicus on two anemone hosts on coral reefs in Panama, and found that while fecundity varies among host species, this variation is explained largely by host size, not species. This suggests that shrimp on larger hosts may have higher fitness regardless of host species, which in turn could drive selection for host choice, a proposed driver of diversification in this group.

Indirect effects of a key ecosystem engineer alter survival and growth of foundation coral species.

Stegastes nigricans, a "farmerfish" that cultivates algal turf and defends territories from grazers and other intruders, can affect coral indirectly due to increased competition with farmed algal turf and/or reduced predation resulting from territorial aggression directed at corallivores. To investigate the indirect effects of this key ecosystem engineer on coral mortality and growth, we transplanted caged and exposed fragments of four coral species to patch reefs in French Polynesia on which we manipulated the presence of S. nigricans and turf, and to reefs naturally devoid of S. nigricans. Reef access was two to four times higher for herbivorous fishes, and two times higher for corallivorous fishes, when S. nigricans was removed, indicating that reef access is reduced for two important guilds of fishes when S. nigricans is present. Stegastes' territoriality indirectly benefited delicate acroporids (Montipora floweri and Acropora striata), yielding a twofold to fivefold reduction in skeletal loss due to lower predation frequencies in the presence of S. nigricans. Three corals, A. striata, M. floweri, and especially Porites australiensis, suffered mortality due to overgrowth significantly more frequently in the presence of farmed turf, but Pocillopora verrucosa did not. Algal abundance predicted the frequency of overgrowth for only A. striata and P. australiensis. M. floweri were more likely to be overgrown when exposed (uncaged) in the presence of S. nigricans, suggesting an interaction modification, in this case that initial predation increased susceptibility to competition with turf. In this community, the presence of S. nigricans may increase algal overgrowth of massive Porites by facilitating its turf competitors and simultaneously reduce predation of branching corals through territorial exclusion of corallivores. These indirect interactions may underlie previously documented community transitions from disturbance-resistant massive coral to recovering branching corals within S. nigricans territories.