Tree of life metabarcoding can serve as a biotic benchmark for shifting baselines in urbanized estuaries.

Urbanization of estuaries drastically changed existing shorelines and bathymetric contours, in turn modifying habitat for marine foundational species that host critical biodiversity. And yet we lack approaches to characterize a significant fraction of the biota that inhabit these ecosystems on time scales that align with rates of urbanization. Environmental DNA (or eDNA) metabarcoding that combines multiple assays targeting a broad range of taxonomic groups can provide a solution, but we need to determine whether the biological communities it detects ally with different habitats in these changing aquatic environments. In this study, we tested whether tree of life metabarcoding (ToL-metabarcoding) data extracted from filtered seawater samples correlated with four known geomorphic habitat zones across a heavily urbanized estuary (Sydney Harbour, Australia). Using this method, we substantially expanded our knowledge on the composition and spatial distribution of marine biodiversity across the tree of life in Sydney Harbour, particularly for organisms where existing records are sparse. Excluding terrestrial DNA inputs, we identified significant effects of both distance from the mouth of Sydney Harbour and geomorphic zone on biological community structure in the ToL-metabarcoding dataset (entire community), as well as in each of the taxonomic subgroups that we considered (fish, macroinvertebrates, algae and aquatic plants, bacteria). This effect appeared to be driven by taxa as a collective versus a few individual taxa, with each taxon explaining no more than 0.62% of the variation between geomorphic zones. Similarly, taxonomic richness was significantly higher within geomorphic zones with large sample sizes, but also decreased by 1% with each additional kilometer from the estuary mouth, a result consistent with a reduction in tidal inputs and available habitat in upper catchments. Based on these results, we suggest that ToL-metabarcoding can be used to benchmark biological monitoring in other urbanized estuaries globally, and in Sydney Harbour at future time points based on detection of bioindicators across the tree of life. We also suggest that robust biotic snapshots can be archived following extensive curation of taxonomic assignments that incorporates ecological affinities, supported by records from relevant and regional biodiversity repositories.

Phylogenomic Analysis of Concatenated Ultraconserved Elements Reveals the Recent Evolutionary Radiation of the Fairy Wrasses (Teleostei: Labridae: Cirrhilabrus).

The fairy wrasses (genus Cirrhilabrus) are among the most successful of the extant wrasse lineages (Teleostei: Labridae), with their 61 species accounting for nearly 10$\%$ of the family. Although species complexes within the genus have been diagnosed on the basis of coloration patterns and synapomorphies, attempts to resolve evolutionary relationships among these groups using molecular and morphological data have largely been unsuccessful. Here, we use a phylogenomic approach with a data set comprising 991 ultraconserved elements (UCEs) and mitochondrial COI to uncover the evolutionary history and patterns of temporal and spatial diversification of the fairy wrasses. Our analyses of phylogenetic signal suggest that most gene-tree incongruence is caused by estimation error, leading to poor resolution in a summary-coalescent analysis of the data. In contrast, analyses of concatenated sequences are able to resolve the major relationships of Cirrhilabrus. We determine the placements of species that were previously regarded as incertae sedis and find evidence for the nesting of Conniella, an unusual, monotypic genus, within Cirrhilabrus. Our relaxed-clock dating analysis indicates that the major divergences within the genus occurred around the Miocene-Pliocene boundary, followed by extensive cladogenesis of species complexes in the Pliocene-Pleistocene. Biogeographic reconstruction suggests that the fairy wrasses emerged within the Coral Triangle, with episodic fluctuations of sea levels during glacial cycles coinciding with shallow divergence events but providing few opportunities for more widespread dispersal. Our study demonstrates both the resolving power and limitations of UCEs across shallow timescales where there is substantial estimation error in individual gene trees.[Biogeography; concatenation; gene genealogy interrogation; gene trees; molecular dating; summary coalescent; UCEs.].

Hyporthodus griseofasciatus (Perciformes: Epinephelidae), a new species of deep-water grouper from the west coast of Australia.

A new species of deep-water epinephelid fish is described from the west coast of Australia based on 14 specimens, 99-595 mm standard length. Hyporthodus griseofasciatus sp. nov. is endemic to Western Australia from Barrow Island to Two Peoples Bay in depths of 76-470 m. It has a series of eight grey bands alternating with eight brown bands along the body and the soft dorsal, soft anal and caudal fin margins are pale cream to white. It is distinguished from its nearest congener, H. ergastularius, by the presence of a star-like pattern of radiating lines on the head versus an overall brownish colour in the latter as well as significant differences in the quantitative analyses of 25 morphological characters. The two species have allopatric distributions on either side of the Australian continent. H. griseofasciatus is distinguished from H. octofasciatus by several grey bands being distinctly narrower than other grey bands (vs. all grey bands subequal in the latter) and the presence of broad white margins on the dorsal, caudal and anal fins (vs. narrow or absent in the latter). Some scale counts appear to also differ. Analyses of mitochondrial cytochrome oxidase subunit 1 sequences revealed reciprocally monophyletic clades with fixed differences and genetic distances typical of recently diverged species of fishes.

Angels in disguise: sympatric hybridization in the marine angelfishes is widespread and occurs between deeply divergent lineages.

Hybridization events are not uncommon in marine environments where physical barriers are attenuated. Studies of coral reef taxa have suggested that hybridization predominantly occurs between parapatric species distributed along biogeographic suture zones. By contrast, little is known about the extent of sympatric hybridization on coral reefs, despite the large amount of biogeographic overlap shared by many coral reef species. Here, we investigate if the propensity for hybridization along suture zones represents a general phenomenon among coral reef fishes, by focusing on the marine angelfishes (family Pomacanthidae). Although hybridization has been reported for this family, it has not been thoroughly surveyed, with more recent hybridization studies focusing instead on closely related species from a population genetics perspective. We provide a comprehensive survey of hybridization among the Pomacanthidae, characterize the upper limits of genetic divergences between hybridizing species and investigate the occurrence of sympatric hybridization within this group. We report the occurrence of hybridization involving 42 species (48% of the family) from all but one genus of the Pomacanthidae. Our results indicate that the marine angelfishes are among the groups of coral reef fishes with the highest incidences of hybridization, not only between sympatric species, but also between deeply divergent lineages.

eDNA metabarcoding survey reveals fine-scale coral reef community variation across a remote, tropical island ecosystem.

Environmental DNA (eDNA) metabarcoding, a technique for retrieving multispecies DNA from environmental samples, can detect a diverse array of marine species from filtered seawater samples. There is a growing potential to integrate eDNA alongside existing monitoring methods in order to establish or improve the assessment of species diversity. Remote island reefs are increasingly vulnerable to climate-related threats and as such there is a pressing need for cost-effective whole-ecosystem surveying to baseline biodiversity, study assemblage changes and ultimately develop sustainable management plans. We investigated the utility of eDNA metabarcoding as a high-resolution, multitrophic biomonitoring tool at the Cocos (Keeling) Islands, Australia (CKI)-a remote tropical coral reef atoll situated within the eastern Indian Ocean. Metabarcoding assays targeting the mitochondrial 16S rRNA and CO1 genes, as well as the 18S rRNA nuclear gene, were applied to 252 surface seawater samples collected from 42 sites within a 140 km2 area. Our assays successfully detected a wide range of bony fish and elasmobranchs (244 taxa), crustaceans (88), molluscs (37) and echinoderms (7). Assemblage composition varied significantly between sites, reflecting habitat partitioning across the island ecosystem and demonstrating the localisation of eDNA signals, despite extensive tidal and oceanic movements. In addition, we document putative new occurrence records for 46 taxa and compare the efficiency of our eDNA approach to visual survey techniques at CKI. Our study demonstrates the utility of a multimarker metabarcoding approach in capturing multitrophic biodiversity across an entire coral reef atoll and sets an important baseline for ongoing monitoring and management.

Combined use of eDNA metabarcoding and video surveillance for the assessment of fish biodiversity.

Monitoring communities of fish is important for the management and sustainability of fisheries and marine ecosystems. Baited remote underwater video systems (BRUVs) are among the most effective nondestructive techniques for sampling bony fishes and elasmobranchs (sharks, rays, and skates). However, BRUVs sample visually conspicuous biota; hence, some taxa are undersampled or not recorded at all. We compared the diversity of fishes characterized using BRUVs with diversity detected via environmental DNA (eDNA) metabarcoding. We sampled seawater and captured BRUVs imagery at 48 locales that included reef and seagrass beds inside and outside a marine reserve (Jurien Bay in Western Australia). Eighty-two fish genera from 13 orders were detected, and the community of fishes described using eDNA and BRUVs combined yielded >30% more generic richness than when either method was used alone. Rather than detecting a homogenous genetic signature, the eDNA assemblages mirrored the BRUVs' spatial explicitness; differentiation of taxa between seagrass and reef was clear despite the relatively small geographical scale of the study site (∼35 km2 ). Taxa that were not sampled by one approach, due to limitations and biases intrinsic to the method, were often detected with the other. Therefore, using BRUVs and eDNA in concert provides a more holistic view of vertebrate marine communities across habitats. Both methods are noninvasive, which enhances their potential for widespread implementation in the surveillance of marine ecosystems.

Comparative phylogeography of the ocean planet.

Understanding how geography, oceanography, and climate have ultimately shaped marine biodiversity requires aligning the distributions of genetic diversity across multiple taxa. Here, we examine phylogeographic partitions in the sea against a backdrop of biogeographic provinces defined by taxonomy, endemism, and species composition. The taxonomic identities used to define biogeographic provinces are routinely accompanied by diagnostic genetic differences between sister species, indicating interspecific concordance between biogeography and phylogeography. In cases where individual species are distributed across two or more biogeographic provinces, shifts in genotype frequencies often align with biogeographic boundaries, providing intraspecific concordance between biogeography and phylogeography. Here, we provide examples of comparative phylogeography from (i) tropical seas that host the highest marine biodiversity, (ii) temperate seas with high productivity but volatile coastlines, (iii) migratory marine fauna, and (iv) plankton that are the most abundant eukaryotes on earth. Tropical and temperate zones both show impacts of glacial cycles, the former primarily through changing sea levels, and the latter through coastal habitat disruption. The general concordance between biogeography and phylogeography indicates that the population-level genetic divergences observed between provinces are a starting point for macroevolutionary divergences between species. However, isolation between provinces does not account for all marine biodiversity; the remainder arises through alternative pathways, such as ecological speciation and parapatric (semiisolated) divergences within provinces and biodiversity hotspots.

Risky business for a juvenile marine predator? Testing the influence of foraging strategies on size and growth rate under natural conditions.

Mechanisms driving selection of body size and growth rate in wild marine vertebrates are poorly understood, thus limiting knowledge of their fitness costs at ecological, physiological and genetic scales. Here, we indirectly tested whether selection for size-related traits of juvenile sharks that inhabit a nursery hosting two dichotomous habitats, protected mangroves (low predation risk) and exposed seagrass beds (high predation risk), is influenced by their foraging behaviour. Juvenile sharks displayed a continuum of foraging strategies between mangrove and seagrass areas, with some individuals preferentially feeding in one habitat over another. Foraging habitat was correlated with growth rate, whereby slower growing, smaller individuals fed predominantly in sheltered mangroves, whereas larger, faster growing animals fed over exposed seagrass. Concomitantly, tracked juveniles undertook variable movement behaviours across both the low and high predation risk habitat. These data provide supporting evidence for the hypothesis that directional selection favouring smaller size and slower growth rate, both heritable traits in this shark population, may be driven by variability in foraging behaviour and predation risk. Such evolutionary pathways may be critical to adaptation within predator-driven marine ecosystems.

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.

Regal phylogeography: Range-wide survey of the marine angelfish Pygoplites diacanthus reveals evolutionary partitions between the Red Sea, Indian Ocean, and Pacific Ocean.

The regal angelfish (Pygoplites diacanthus; family Pomacanthidae) occurs on reefs from the Red Sea to the central Pacific, with an Indian Ocean/Rea Sea color morph distinct from a Pacific Ocean morph. To assess population differentiation and evaluate the possibility of cryptic evolutionary partitions in this monotypic genus, we surveyed mtDNA cytochrome b and two nuclear introns (S7 and RAG2) in 547 individuals from 15 locations. Phylogeographic analyses revealed four mtDNA lineages (d=0.006-0.015) corresponding to the Pacific Ocean, the Red Sea, and two admixed lineages in the Indian Ocean, a pattern consistent with known biogeographic barriers. Christmas Island in the eastern Indian Ocean had both Indian and Pacific lineages. Both S7 and RAG2 showed strong population-level differentiation between the Red Sea, Indian Ocean, and Pacific Ocean (ΦST=0.066-0.512). The only consistent population sub-structure within these three regions was at the Society Islands (French Polynesia), where surrounding oceanographic conditions may reinforce isolation. Coalescence analyses indicate the Pacific (1.7Ma) as the oldest extant lineage followed by the Red Sea lineage (1.4Ma). Results from a median-joining network suggest radiations of two lineages from the Red Sea that currently occupy the Indian Ocean (0.7-0.9Ma). Persistence of a Red Sea lineage through Pleistocene glacial cycles suggests a long-term refuge in this region. The affiliation of Pacific and Red Sea populations, apparent in cytochrome b and S7 (but equivocal in RAG2) raises the hypothesis that the Indian Ocean was recolonized from the Red Sea, possibly more than once. Assessing the genetic architecture of this widespread monotypic genus reveals cryptic evolutionary diversity that merits subspecific recognition. We recommend P.d. diacanthus and P.d. flavescens for the Pacific and Indian Ocean/Red Sea forms.

Surgeons and suture zones: Hybridization among four surgeonfish species in the Indo-Pacific with variable evolutionary outcomes.

Closely related species can provide valuable insights into evolutionary processes through comparison of their ecology, geographic distribution and the history recorded in their genomes. In the Indo-Pacific, many reef fishes are divided into sister species that come into secondary contact at biogeographic borders, most prominently where Indian Ocean and Pacific Ocean faunas meet. It is unclear whether hybridization in this contact zone represents incomplete speciation, secondary contact, an evolutionary dead-end (for hybrids) or some combination of the above. To address these issues, we conducted comprehensive surveys of two widely-distributed surgeonfish species, Acanthurus leucosternon (N=141) and A. nigricans (N=412), with mtDNA cytochrome b sequences and ten microsatellite loci. These surgeonfishes are found primarily in the Indian and Pacific Oceans, respectively, but overlap at the Christmas and Cocos-Keeling Islands hybrid zone in the eastern Indian Ocean. We also sampled the two other Pacific members of this species complex, A. achilles (N=54) and A. japonicus (N=49), which are known to hybridize with A. nigricans where their ranges overlap. Our results indicate separation between the four species that range from the recent Pleistocene to late Pliocene (235,000-2.25million years ago). The Pacific A. achilles is the most divergent (and possibly ancestral) species with mtDNA dcorr≈0.04, whereas the other two Pacific species (A. japonicus and A. nigricans) are distinguishable only at a population or subspecies level (ΦST=0.6533, P<0.001). Little population structure was observed within species, with evidence of recent population expansion across all four geographic ranges. We detected sharing of mtDNA haplotypes between species and extensive hybridization based on microsatellites, consistent with later generation hybrids but also the effects of allele homoplasy. Despite extensive introgression, 98% of specimens had concordance between mtDNA lineage and species identification based on external morphology, indicating that species integrity may not be eroding. The A. nigricans complex demonstrates a range of outcomes from incomplete speciation to secondary contact to decreasing hybridization with increasing evolutionary depth.

Angelfishes, Paper Tigers, and the Devilish Taxonomy of the Centropyge flavissima Complex.

The pygmy angelfishes (genus Centropyge) provide numerous examples of discordance between color morphology, taxonomy, and evolutionary genetic lineages. This discordance is especially evident in the Centropyge flavissima complex, which includes three primary color morphs, three previously recognized species (C. flavissima, Centropyge eibli, and Centropyge vrolikii) and three distinct mitochondrial (mtDNA) lineages that do not align with species designations. Our previous research showed that the putative C. flavissima arose independently in the Indian and Pacific Oceans, and the three mtDNA lineages align with geography rather than species assignments. Here, we add 157 specimens to the previous dataset of 291 specimens, spread across a greater geographic range, to pinpoint the distribution of mtDNA lineages and color morphs. We found that the mtDNA lineages show remarkably strong geographic boundaries corresponding to the Indian Ocean, Central-West Pacific, and Central-South Pacific. We also test the validity of the "Black Tiger Centropyge" in the C. flavissima species complex, a taxonomic oddity that is restricted to shoals and atolls off the coast of northwestern Australia, and the newly named Centropyge cocosensis assigned to the C. flavissima lineage in the Indian Ocean. We conclude that the Black Tiger Centropyge is not a valid species but rather an intermediate between sympatric color morphs that correspond to the putative species C. eibli and C. vrolikii Our greater sampling efforts also do not support the genetic distinctiveness of C. cocosensis given shared mtDNA haplotypes with the sympatric C. eibli and C. vrolikii, but instead we find conflicting lines of evidence concerning the taxonomy of this group. We urge caution and taxonomic restraint until the true nature of this species complex can be revealed.

Seascape genetics along environmental gradients in the Arabian Peninsula: insights from ddRAD sequencing of anemonefishes.

Understanding the processes that shape patterns of genetic structure across space is a central aim of landscape genetics. However, it remains unclear how geographical features and environmental variables shape gene flow, particularly for marine species in large complex seascapes. Here, we evaluated the genomic composition of the two-band anemonefish Amphiprion bicinctus across its entire geographical range in the Red Sea and Gulf of Aden, as well as its close relative, Amphiprion omanensis endemic to the southern coast of Oman. Both the Red Sea and the Arabian Sea are complex and environmentally heterogeneous marine systems that provide an ideal scenario to address these questions. Our findings confirm the presence of two genetic clusters previously reported for A. bicinctus in the Red Sea. Genetic structure analyses suggest a complex seascape configuration, with evidence of both isolation by distance (IBD) and isolation by environment (IBE). In addition to IBD and IBE, genetic structure among sites was best explained when two barriers to gene flow were also accounted for. One of these coincides with a strong oligotrophic-eutrophic gradient at around 16-20˚N in the Red Sea. The other agrees with a historical bathymetric barrier at the straight of Bab al Mandab. Finally, these data support the presence of interspecific hybrids at an intermediate suture zone at Socotra and indicate complex patterns of genomic admixture in the Gulf of Aden with evidence of introgression between species. Our findings highlight the power of recent genomic approaches to resolve subtle patterns of gene flow in marine seascapes.

Two decades of genetic profiling yields first evidence of natal philopatry and long-term fidelity to parturition sites in sharks.

Sharks are a globally threatened group of marine fishes that often breed in their natal region of origin. There has even been speculation that female sharks return to their exact birthplace to breed ('natal philopatry'), which would have important conservation implications. Genetic profiling of lemon sharks (Negaprion brevirostris) from 20 consecutive cohorts (1993-2012) at Bimini, Bahamas, showed that certain females faithfully gave birth at this site for nearly two decades. At least six females born in the 1993-1997 cohorts returned to give birth 14-17 years later, providing the first direct evidence of natal philopatry in the chondrichthyans. Long-term fidelity to specific nursery sites coupled with natal philopatry highlights the merits of emerging spatial and local conservation efforts for these threatened predators.

Evolution of pygmy angelfishes: recent divergences, introgression, and the usefulness of color in taxonomy.

The pygmy angelfishes (genus Centropyge, family Pomacanthidae) are brightly colored species that occupy reef habitats in every tropical ocean. Some species are rarely observed because they occur below conventional scuba depths. Their striking coloration can command thousands of U.S. dollars in the aquarium trade, and closely related species are often distinguished only by coloration. These factors have impeded phylogenetic resolution, and every phylogeographic survey to date has reported discordance between coloration, taxonomy, and genetic partitions. Here we report a phylogenetic survey of 29 of the 34 recognized species (N=94 plus 23 outgroups), based on two mtDNA and three nuclear loci, totaling 2272 bp. The resulting ML and Baysian trees are highly concordant and indicate that the genus Centropyge is paraphyletic, consistent with a previous analysis of the family Pomacanthidae. Two recognized genera (Apolemichthys and Genicanthus) nest within Centropyge, and two subgenera (Xiphypops and Paracentropyge) comprise monophyletic lineages that should be elevated to genus level. Based on an age estimate of 38 Ma for the family Pomacanthidae, Centropyge diverged from the closest extant genus Pygoplites about 33 Ma, three deep lineages within Centropyge diverged about 18-28 Ma, and four species complexes diverged 3-12 Ma. However, in 11 of 13 cases, putative species in these complexes are indistinguishable based on morphology and genetics, being defined solely by coloration. These cases indicate either emerging species or excessive taxonomic splitting based on brightly colored variants.

The spotted parrotfish genome provides insights into the evolution of a coral reef dietary specialist (Teleostei: Labridae: Scarini: Cetoscarus ocellatus).

With over 600 valid species, the wrasses (family Labridae) are among the largest and most successful families of the marine teleosts. They feature prominently on coral reefs where they are known not only for their impressive diversity in colouration and form but also for their functional specialisation and ability to occupy a wide variety of trophic guilds. Among the wrasses, the parrotfishes (tribe Scarini) display some of the most dramatic examples of trophic specialisation. Using abrasion-resistant biomineralized teeth, parrotfishes are able to mechanically extract protein-rich micro-photoautotrophs growing in and among reef carbonate material, a dietary niche that is inaccessible to most other teleost fishes. This ability to exploit an otherwise untapped trophic resource is thought to have played a role in the diversification and evolutionary success of the parrotfishes. In order to better understand the key evolutionary innovations leading to the success of these dietary specialists, we sequenced and analysed the genome of a representative species, the spotted parrotfish (Cetoscarus ocellatus). We find significant expansion, selection and duplications within several detoxification gene families and a novel poly-glutamine expansion in the enamel protein ameloblastin, and we consider their evolutionary implications. Our genome provides a useful resource for comparative genomic studies investigating the evolutionary history of this highly specialised teleostean radiation.

The spatial patterns of directional phenotypic selection.

Local adaptation, adaptive population divergence and speciation are often expected to result from populations evolving in response to spatial variation in selection. Yet, we lack a comprehensive understanding of the major features that characterise the spatial patterns of selection, namely the extent of variation among populations in the strength and direction of selection. Here, we analyse a data set of spatially replicated studies of directional phenotypic selection from natural populations. The data set includes 60 studies, consisting of 3937 estimates of selection across an average of five populations. We performed meta-analyses to explore features characterising spatial variation in directional selection. We found that selection tends to vary mainly in strength and less in direction among populations. Although differences in the direction of selection occur among populations they do so where selection is often weakest, which may limit the potential for ongoing adaptive population divergence. Overall, we also found that spatial variation in selection appears comparable to temporal (annual) variation in selection within populations; however, several deficiencies in available data currently complicate this comparison. We discuss future research needs to further advance our understanding of spatial variation in selection.

Gut microbial communities of hybridising pygmy angelfishes reflect species boundaries.

Hybridisation and introgression of eukaryotic genomes can generate new species or subsume existing ones, with direct and indirect consequences for biodiversity. An understudied component of these evolutionary forces is their potentially rapid effect on host gut microbiomes, and whether these pliable microcosms may serve as early biological indicators of speciation. We address this hypothesis in a field study of angelfishes (genus Centropyge), which have one of the highest prevalence of hybridisation within coral reef fish. In our study region of the Eastern Indian Ocean, the parent fish species and their hybrids cohabit and display no differences in their diet, behaviour, and reproduction, often interbreeding in mixed harems. Despite this ecological overlap, we show that microbiomes of the parent species are significantly different from each other in form and function based on total community composition, supporting the division of parents into distinct species, despite the confounding effects of introgression acting to homogenize parent species identity at other molecular markers. The microbiome of hybrid individuals, on the other hand, are not significantly different to each of the parents, instead harbouring an intermediate community composition. These findings suggest that shifts in gut microbiomes may be an early indicator of speciation in hybridising species.

Gut content metabarcoding of specialized feeders is not a replacement for environmental DNA assays of seawater in reef environments.

In tropical marine ecosystems, the coral-based diet of benthic-feeding reef fishes provides a window into the composition and health of coral reefs. In this study, for the first time, we compare multi-assay metabarcoding sequences of environmental DNA (eDNA) isolated from seawater and partially digested gut items from an obligate corallivore butterflyfish (Chaetodon lunulatus) resident to coral reef sites in the South China Sea. We specifically tested the proportional and statistical overlap of the different approaches (seawater vs gut content metabarcoding) in characterizing eukaryotic community composition on coral reefs. Based on 18S and ITS2 sequence data, which differed in their taxonomic sensitivity, we found that gut content detections were only partially representative of the eukaryotic communities detected in the seawater based on low levels of taxonomic overlap (3 to 21%) and significant differences between the sampling approaches. Overall, our results indicate that dietary metabarcoding of specialized feeders can be complimentary to, but is no replacement for, more comprehensive environmental DNA assays of reef environments that might include the processing of different substrates (seawater, sediment, plankton) or traditional observational surveys. These molecular assays, in tandem, might be best suited to highly productive but cryptic oceanic environments (kelp forests, seagrass meadows) that contain an abundance of organisms that are often small, epiphytic, symbiotic, or cryptic.

Phylogeography of two closely related Indo-Pacific butterflyfishes reveals divergent evolutionary histories and discordant results from mtDNA and microsatellites.

Marine biogeographic barriers can have unpredictable consequences, even among closely related species. To resolve phylogeographic patterns for Indo-Pacific reef fauna, we conducted range-wide surveys of sister species, the scrawled butterflyfish (Chaetodon meyeri; N = 134) and the ornate butterflyfish (Chaetodon ornatissimus; N = 296), using mitochondrial DNA cytochrome b sequences and 10 microsatellite loci. The former is distributed primarily in the Indian Ocean but also extends to the Line Islands in the Central Pacific, whereas the latter is distributed primarily in the Central-West Pacific (including Hawaii and French Polynesia) but extends to the eastern margin of the Indian Ocean. Analyses of molecular variance and Bayesian STRUCTURE results revealed 1 range-wide group for C. meyeri and 3 groups for C. ornatissimus: 1) eastern Indian Ocean and western Pacific, 2) Central Pacific, and 3) Hawaii. Estimates of the last population expansion were much more recent for C. meyeri (61 500 to 95 000 years) versus C. ornatissimus (184 700 to 286 300 years). Despite similarities in ecology, morphology, life history, and a broadly overlapping distribution, these sister species have divergent patterns of dispersal and corresponding evolutionary history. The mtDNA and microsatellite markers did not provide concordant results within 1 of our study species (C. meyeri), or in 7 out of 12 other cases of marine fishes in the published literature. This discordance renews caution in relying on one or a few markers for reconstructing historical demography.