Seascape genomics reveals fine‐scale patterns of dispersal for a reef fish along the ecologically divergent coast of Northwestern Australia

Understanding the drivers of dispersal among populations is a central topic in marine ecology and fundamental for spatially explicit management of marine resources. The extensive coast of Northwestern Australia provides an emerging frontier for implementing new genomic tools to comparatively identify patterns of dispersal across diverse and extreme environmental conditions. Here, we focused on the stripey snapper (Lutjanus carponotatus), which is important to recreational, charter‐based and customary fishers throughout the Indo‐West Pacific. We collected 1,016 L. carponotatus samples at 51 locations in the coastal waters of Northwestern Australia ranging from the Northern Territory to Shark Bay and adopted a genotype‐by‐sequencing approach to test whether realized connectivity (via larval dispersal) was related to extreme gradients in coastal hydrodynamics. Hydrodynamic simulations using CONNIE and a more detailed treatment in the Kimberley Bioregion provided null models for comparison. Based on 4,402 polymorphic single nucleotide polymorphism loci shared across all individuals, we demonstrated significant genetic subdivision between the Shark Bay Bioregion in the south and all locations within the remaining, more northern bioregions. More importantly, we identified a zone of admixture spanning a distance of 180 km at the border of the Kimberley and Canning bioregions, including the Buccaneer Archipelago and adjacent waters, which collectively experiences the largest tropical tidal range and some of the fastest tidal currents in the world. Further testing of the generality of this admixture zone in other shallow water species across broader geographic ranges will be critical for our understanding of the population dynamics and genetic structure of marine taxa in our tropical oceans.     

Spatial patterns of self‐recruitment of a coral reef fish in relation to island‐scale retention mechanisms

Oceanographic features influence the transport and delivery of marine larvae, and physical retention mechanisms, such as eddies, can enhance self‐recruitment (i.e. the return of larvae to their natal population). Knowledge of exact locations of hatching (origin) and settlement (arrival) of larvae of reef animals provides a means to compare observed patterns of self‐recruitment ‘connectivity’ with those expected from water circulation patterns. Using parentage inference based on multiple sampling years in Moorea, French Polynesia, we describe spatial and temporal variation in self‐recruitment of the anemonefish Amphiprion chrysopterus, evaluate the consistency of net dispersal distances of self‐recruits against the null expectation of passive particle dispersal and test the hypothesis that larvae originating in certain reef habitats (lagoons and passes) would be retained and thus more likely to self‐recruit than those originating on the outer (fore) reef. Estimates of known self‐recruitment were consistent across the sampling years (~25–27% of sampled recruits). For most (88%) of these self‐recruits, the net distance between hatching and settlement locations was within the maximum dispersal distance expected for a neutrally buoyant passive particle based on the longest duration of the larval dispersive phase and the average direction and speed of current flow around Moorea. Furthermore, a parent of a given body size on the outer (fore) reef of Moorea was less likely to produce self‐recruits than those in passes. Our findings show that even a simple dispersal model based on net average flow and direction of alongshore currents can provide insight into landscape‐scale retention patterns of reef fishes.     

Carrion fly‐derived DNA metabarcoding is an effective tool for mammal surveys: Evidence from a known tropical mammal community

Metabarcoding of vertebrate DNA derived from carrion flies has been proposed as a promising tool for biodiversity monitoring. To evaluate its efficacy, we conducted metabarcoding surveys of carrion flies on Barro Colorado Island (BCI), Panama, which has a well‐known mammal community, and compared our results against diurnal transect counts and camera trapping. We collected 1,084 flies in 29 sampling days, conducted metabarcoding with mammal‐specific (16S) and vertebrate‐specific (12S) primers, and sequenced amplicons on Illumina MiSeq. For taxonomic assignment, we compared blast with the new program protax , and we found that protax improved species identifications. We detected 20 mammal, four bird, and one lizard species from carrion fly metabarcoding, all but one of which are known from BCI. Fly metabarcoding detected more mammal species than concurrent transect counts (29 sampling days, 13 species) and concurrent camera trapping (84 sampling days, 17 species), and detected 67% of the number of mammal species documented by 8 years of transect counts and camera trapping combined, although fly metabarcoding missed several abundant species. This study demonstrates that carrion fly metabarcoding is a powerful tool for mammal biodiversity surveys and has the potential to detect a broader range of species than more commonly used methods.     

Evidence of stable genetic structure across a remote island archipelago through self‐recruitment in a widely dispersed coral reef fish

We used microsatellite markers to assess the population genetic structure of the scribbled rabbitfish Siganus spinus in the western Pacific. This species is a culturally important food fish in the Mariana Archipelago and subject to high fishing pressure. Our primary hypothesis was to test whether the individuals resident in the southern Mariana Island chain were genetically distinct and hence should be managed as discrete stocks. In addition to spatial sampling of adults, newly‐settled individuals were sampled on Guam over four recruitment events to assess the temporal stability of the observed spatial patterns, and evidence of self‐recruitment. We found significant genetic structure in S. spinus across the western Pacific, with Bayesian analyses revealing three genetically distinct clusters: the southern Mariana Islands, east Micronesia, and the west Pacific; with the southern Mariana Islands being more strongly differentiated from the rest of the region. Analyses of temporal samples from Guam indicated the southern Mariana cluster was stable over time, with no genetic differentiation between adults versus recruits, or between samples collected across four separate recruitment events spanning 11 months. Subsequent assignment tests indicated seven recruits had self‐recruited from within the Southern Mariana Islands population. Our results confirm the relative isolation of the southern Mariana Islands population and highlight how local processes can act to isolate populations that, by virtue of their broad‐scale distribution, have been subject to traditionally high gene flows. Our results add to a growing consensus that self‐recruitment is a highly significant influence on the population dynamics of tropical reef fish.     

The community and ecosystem consequences of intraspecific diversity: a meta‐analysis

Understanding the relationships between biodiversity and ecosystem functioning has major implications. Biodiversity–ecosystem functioning relationships are generally investigated at the interspecific level, although intraspecific diversity (i.e. within‐species diversity) is increasingly perceived as an important ecological facet of biodiversity. Here, we provide a quantitative and integrative synthesis testing, across diverse plant and animal species, whether intraspecific diversity is a major driver of community dynamics and ecosystem functioning. We specifically tested (i) whether the number of genotypes/phenotypes (i.e. intraspecific richness) or the specific identity of genotypes/phenotypes (i.e. intraspecific variation) in populations modulate the structure of communities and the functioning of ecosystems, (ii) whether the ecological effects of intraspecific richness and variation are strong in magnitude, and (iii) whether these effects vary among taxonomic groups and ecological responses. We found a non‐linear relationship between intraspecific richness and community and ecosystem dynamics that follows a saturating curve shape, as observed for biodiversity–function relationships measured at the interspecific level. Importantly, intraspecific richness modulated ecological dynamics with a magnitude that was equal to that previously reported for interspecific richness. Our results further confirm, based on a database containing more than 50 species, that intraspecific variation also has substantial effects on ecological dynamics. We demonstrated that the effects of intraspecific variation are twice as high as expected by chance, and that they might have been underestimated previously. Finally, we found that the ecological effects of intraspecific variation are not homogeneous and are actually stronger when intraspecific variation is manipulated in primary producers than in consumer species, and when they are measured at the ecosystem rather than at the community level. Overall, we demonstrated that the two facets of intraspecific diversity (richness and variation) can both strongly affect community and ecosystem dynamics, which reveals the pivotal role of within‐species biodiversity for understanding ecological dynamics.     

Whole‐genome patterns of linkage disequilibrium across flycatcher populations clarify the causes and consequences of fine‐scale recombination rate variation in birds

Recombination rate is heterogeneous across the genome of various species and so are genetic diversity and differentiation as a consequence of linked selection. However, we still lack a clear picture of the underlying mechanisms for regulating recombination. Here we estimated fine‐scale population recombination rate based on the patterns of linkage disequilibrium across the genomes of multiple populations of two closely related flycatcher species (Ficedula albicollis and F. hypoleuca). This revealed an overall conservation of the recombination landscape between these species at the scale of 200 kb, but we also identified differences in the local rate of recombination despite their recent divergence (<1 million years). Genetic diversity and differentiation were associated with recombination rate in a lineage‐specific manner, indicating differences in the extent of linked selection between species. We detected 400–3,085 recombination hotspots per population. Location of hotspots was conserved between species, but the intensity of hotspot activity varied between species. Recombination hotspots were primarily associated with CpG islands (CGIs), regardless of whether CGIs were at promoter regions or away from genes. Recombination hotspots were also associated with specific transposable elements (TEs), but this association appears indirect due to shared preferences of the transposition machinery and the recombination machinery for accessible open chromatin regions. Our results suggest that CGIs are a major determinant of the localization of recombination hotspots, and we propose that both the distribution of TEs and fine‐scale variation in recombination rate may be associated with the evolution of the epigenetic landscape.     

Long‐distance dispersal and local retention of larvae as mechanisms of recruitment in an island population of a coral reef fish

Abstract Although recruitment of pelagic larvae is a fundamental and well‐documented process in the dynamics of benthic marine populations, identifying the sources of recruitment, or the degree to which populations are connected via dispersal of larvae, has remained elusive for most marine taxa. In this study we used natural environmental markers (trace elements) recorded in fish otoliths (ear stones) as tags of natal origin. Specifically, we used the otolith core and edge chemistries of a locally endemic wrasse (Coris bulbifrons) from Lord Howe Island (LHI), Australia, and a widely distributed species (Coris picta) from three potential mainland source regions, to determine the likely sources of recruitment to C. picta populations on LHI. The use of a local endemic species, which is by definition self‐recruiting, is a novel approach for ground‐truthing the dispersal history of non‐endemic coral reef fish. Discriminant function analyses were able to separate LHI from mainland fish, using both edge and core signatures, with a high degree of accuracy, suggesting at least some of the C. picta collected on LHI were of local origin. This result was corroborated when half of the C. bulbifrons and LHI C. picta were introduced as unknowns into a discriminant function analysis using the remaining C. bulbifrons, LHI C. picta, and the mainland C. picta as a training data set. Overall, our findings suggest that both long distance dispersal and local retention are important sources of recruitment to populations of C. picta on LHI and that otolith chemistry of endemic species could be a useful benchmark for determining the prevalence of self‐recruitment in insular populations of other widespread species.     

Effects of diversity loss on ecosystem function across trophic levels and ecosystems: A test in a detritus‐based tropical food web

Abstract We investigated the effects of biodiversity loss across trophic levels and across ecosystems (terrestrial to aquatic) on ecosystem function, in a detritus‐based tropical food web. Diversities of consumers (stream shredders) and resources (leaf litter) were experimentally manipulated by varying the number of species from 3 to 1, using different species combinations, and the effects on leaf breakdown rates were examined. In single‐species shredder treatments, leaf diversity loss affected breakdown rates, but the effect depended on the identity of the leaves remaining in the system: they increased when the most preferred leaf species remained, but decreased when this species was lost (leaf preferences were the same for all shredders). In multi‐species shredder assemblages, breakdown rates were lower than expected from single‐species treatments, suggesting an important role of interspecific competition. This pattern was also evident when oneleaf species was available but not with higher leaf diversity, suggesting that lowered leaf diversity promotes competitive interactions among shredders. The influence of diversity and identity of species across trophic levels and ecosystems on stream functioning points to complex interactions that may well be reflected in other types of ecosystem.     

Cycles of coral reef ‘turn‐on’, rapid growth and ‘turn‐off’ over the past 8500 years: a context for understanding modern ecological states and trajectories

Human activities threaten reef ecosystems globally, forcing ecological change at rates and scales regarded as unprecedented in the Holocene. These changes are so profound that a cessation of reef accretion (reef ‘turn‐off’) and net erosion of reef structures is argued by many as the ultimate and imminent trajectory. Here, we use a regional scale reef growth dataset, based on 76 core records (constrained by 211 radiometric dates) from 22 reefs along and across the inner‐shelf of the Great Barrier Reef, Australia, to examine the timing of different phases of reef initiation (‘turn‐on’), growth and ‘turn‐off’ during the Holocene. This dataset delineates two temporally discrete episodes of reef‐building over the last 8500 years: the first associated with the Holocene transgression‐early highstand period [～8.5–5.5 k calibrated years bp (cal ybp )]; the second since ～2.3 k cal ybp . During both periods, reefs accreted rapidly to sea level before entering late evolutionary states – states naturally characterized by reduced coral cover and low accretion potential – and a clear hiatus occurs between these reef‐building episodes for which no records of reef initiation exist. These transitions mimic those projected under current environmental disturbance regimes, but have been driven entirely by natural forcing factors. Our results demonstrate that, even through the late Holocene, reef health and growth has fluctuated through cycles independent of anthropogenic forcing. Consequently, degraded reef states cannot de facto be considered to automatically reflect increased anthropogenic stress. Indeed, in many cases degraded or nonaccreting reef communities may reflect past reef growth histories (as dictated by reef growth–sea level interactions) as much as contemporary environmental change. Recognizing when changes in reef condition reflect these natural ‘turn‐on’– growth –‘turn‐off’ cycles and how they interact with on‐going human disturbance is critical for effective coral reef management and for understanding future reef ecological trajectories.     

Evaluating the effectiveness of teeth and dorsal fin spines for non‐lethal age estimation of a tropical reef fish,coral trout Plectropomus leopardus

This study investigated whether teeth and dorsal fin spines could be used as non‐lethal methods of age estimation for a vulnerable and highly valued tropical fisheries species, coral trout Plectropomus leopardus. Age estimation of individuals from 2 to 9 years old revealed that dorsal spines represent an accurate ageing method (90% agreement with otoliths) that was more precise [average per cent error (APE) = 4·1, coefficient of variation (c.v .) = 5·8%] than otoliths (APE = 6·2, c.v . = 8·7%). Of the three methods for age estimation (otoliths, dorsal spines and teeth), spines were the most time and cost efficient. An aquarium‐based study also found that removing a dorsal spine or tooth did not affect survivorship or growth of P. leopardus. No annuli were visible in teeth despite taking transverse and longitudinal sections throughout the tooth and trialling several different laboratory methods. Although teeth may not be suitable for estimating age of P. leopardus, dorsal spines appear to be an acceptably accurate, precise and efficient method for non‐lethal ageing of individuals from 2 to 9 years old in this tropical species.     

RAD genotyping reveals fine‐scale genetic structuring and provides powerful population assignment in a widely distributed marine species,the American lobster (Homarus americanus)

Deciphering genetic structure and inferring connectivity in marine species have been challenging due to weak genetic differentiation and limited resolution offered by traditional genotypic methods. The main goal of this study was to assess how a population genomics framework could help delineate the genetic structure of the American lobster (Homarus americanus) throughout much of the species’ range and increase the assignment success of individuals to their location of origin. We genotyped 10 156 filtered SNPs using RAD sequencing to delineate genetic structure and perform population assignment for 586 American lobsters collected in 17 locations distributed across a large portion of the species’ natural distribution range. Our results revealed the existence of a hierarchical genetic structure, first separating lobsters from the northern and southern part of the range (F_CT = 0.0011; P‐value = 0.0002) and then revealing a total of 11 genetically distinguishable populations (mean F_ST = 0.00185; CI: 0.0007–0.0021, P‐value < 0.0002), providing strong evidence for weak, albeit fine‐scale population structuring within each region. A resampling procedure showed that assignment success was highest with a subset of 3000 SNPs having the highest F_ST. Applying Anderson's (Molecular Ecology Resources, 2010, 10, 701) method to avoid ‘high‐grading bias’, 94.2% and 80.8% of individuals were correctly assigned to their region and location of origin, respectively. Lastly, we showed that assignment success was positively associated with sample size. These results demonstrate that using a large number of SNPs improves fine‐scale population structure delineation and population assignment success in a context of weak genetic structure. We discuss the implications of these findings for the conservation and management of highly connected marine species, particularly regarding the geographic scale of demographic independence.