Whole genome assembly and annotation of the King Angelfish (Holacanthus passer) gives insight into the evolution of marine fishes of the Tropical Eastern Pacific.

Holacanthus angelfishes are some of the most iconic marine fishes of the Tropical Eastern Pacific (TEP). However, very limited genomic resources currently exist for the genus. In this study we: (i) assembled and annotated the nuclear genome of the King Angelfish (Holacanthus passer), and (ii) examined the demographic history of H. passer in the TEP. We generated 43.8 Gb of ONT and 97.3 Gb Illumina reads representing 75× and 167× coverage, respectively. The final genome assembly size was 583 Mb with a contig N50 of 5.7 Mb, which captured 97.5% of the complete Actinoterygii Benchmarking Universal Single-Copy Orthologs (BUSCOs). Repetitive elements accounted for 5.09% of the genome, and 33,889 protein-coding genes were predicted, of which 22,984 were functionally annotated. Our demographic analysis suggests that population expansions of H. passer occurred prior to the last glacial maximum (LGM) and were more likely shaped by events associated with the closure of the Isthmus of Panama. This result is surprising, given that most rapid population expansions in both freshwater and marine organisms have been reported to occur globally after the LGM. Overall, this annotated genome assembly provides a novel molecular resource to study the evolution of Holacanthus angelfishes, while facilitating research into local adaptation, speciation, and introgression in marine fishes.

Invasion genomics of lionfish in the Mediterranean Sea.

The rate of biological invasions is steadily increasing, with major ecological and economic impacts accounting for billions of dollars in damage as a result. One spectacular example is the western Atlantic invasion by lionfishes. In the Mediterranean Sea, invasions from the Red Sea via the Suez Canal (termed Lessepsian invasions) comprise more than 100 fish species, including a recent invasion by lionfish. In light of the devastating effects of lionfish in the Caribbean Sea, understanding the dynamics of Mediterranean lionfish invasion is crucial. The Lessepsian lionfish invasion started in 2012, and rapidly spread to the central Mediterranean. Here, we used thousands of RAD seq genomic markers to study the population dynamics of this invasion. While we did not find a reduction in genetic diversity between source (Red Sea) and invasive (Mediterranean) populations (i.e., bottleneck effects), we found evidence of population structure within the invasive range in the Mediterranean Sea. We found that loci that are potentially under selection may play an important role in invasion success (in particular, genes involved in osmoregulation and fin spine sizes). Genomic approaches proved powerful in examining the ecological and evolutionary patterns of successful invaders and may be used as tools to understand and potentially mitigate future invasions.

Chromosome-level genome of the three-spot damselfish, Dascyllus trimaculatus.

Damselfishes (Family: Pomacentridae) are a group of ecologically important, primarily coral reef fishes that include over 400 species. Damselfishes have been used as model organisms to study recruitment (anemonefishes), the effects of ocean acidification (spiny damselfish), population structure, and speciation (Dascyllus). The genus Dascyllus includes a group of small-bodied species, and a complex of relatively larger bodied species, the Dascyllus trimaculatus species complex that is comprised of several species including D. trimaculatus itself. The three-spot damselfish, D. trimaculatus, is a widespread and common coral reef fish species found across the tropical Indo-Pacific. Here, we present the first-genome assembly of this species. This assembly contains 910 Mb, 90% of the bases are in 24 chromosome-scale scaffolds, and the Benchmarking Universal Single-Copy Orthologs score of the assembly is 97.9%. Our findings confirm previous reports of a karyotype of 2n = 47 in D. trimaculatus in which one parent contributes 24 chromosomes and the other 23. We find evidence that this karyotype is the result of a heterozygous Robertsonian fusion. We also find that the D. trimaculatus chromosomes are each homologous with single chromosomes of the closely related clownfish species, Amphiprion percula. This assembly will be a valuable resource in the population genomics and conservation of Damselfishes, and continued studies of the karyotypic diversity in this clade.

A high-quality reference genome of the kelp surfperch, Brachyistius frenatus (Embiotocidae), a wide-ranging Eastern Pacific reef fish with no pelagic larval stage.

The surfperches (family Embiotocidae) are a unique group of mostly marine fishes whose phylogenetic position within the Ovalentaria clade (Percomorpha) is still unresolved. As a result of their viviparity and lack of a dispersive larval stage, surfperches are an excellent model for the study of speciation, gene flow, and local adaptation in the ocean. They are also the target of an immensely popular recreational fishery. Very few high-quality molecular resources, however, are available for this group and only for a single species. Here, we describe a highly complete reference genome for the kelp surfperch, Brachyistius frenatus, assembled using a combination of short-read (Illumina, ~47× coverage) and long-read (Oxford Nanopore Technologies, ~27× coverage) sequencing. The 596 Mb assembly has a completeness level of 98.1% (BUSCO), a contig N50 of 2.6 Mb (n = 56), and a contig N90 of 406.6 kb (n = 293). Comparative analysis revealed a high level of synteny between B. frenatus and its close relative, Embiotoca jacksoni. This assembly will serve as a valuable molecular resource upon which future evolutionary dynamics research will build, such as the investigation of local adaptation and the genomic potential for climate adaptation in wild populations.

The complete mitochondrial genome of Dascyllus trimaculatus (Rüppell, 1829).

Damselfishes (family Pomacentridae) comprise approximately 400 species that play an important ecological role in temperate and coral reefs. Here, for the first time, we assemble and annotate the mitochondrial genome of Dascyllus trimaculatus, the three-spot dascyllus, a planktivorous damselfish that primarily recruits in anemones. The circular genome of D. trimaculatus is 16,967 bp in length and contains 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and a control region. Gene arrangement and codon usage is similar to reported mitochondrial genomes of other damselfish genera, and a phylogenetic analysis of a set of damselfish representatives is consistent with known evolutionary analyses.

Reference genome of the Woolly Sculpin, Clinocottus analis.

Sculpins (Family Cottidae) are generally cold-temperate intertidal reef fishes most commonly found in the North Pacific. As part of the California Conservation Genomics Project (CCGP), we sequenced the genome of the Woolly Sculpin, Clinocottus analis, to establish a genomic model for understanding phylogeographic structure of inshore marine taxa along the California coast. These patterns, in turn, should further inform the design of marine protected areas using dispersal models based on genomic data. The small genome of C. analis is typical of marine fishes at less than 1 Gb (genome size = 538 Mb), and our assembly is near-chromosome level (contig N50 = 9.1 Mb, scaffold N50 = 21 Mb, BUSCO completeness = 97.9%). Within the context of the CCGP, the Woolly Sculpin genome will be used as a reference for future whole-genome resequencing projects aimed at enhancing our knowledge of the population structure of the species, and efficacy of marine protected areas across the state.

Reference genome of the Monkeyface Prickleback, Cebidichthys violaceus.

Pricklebacks (Family Stichaeidae) are generally cold-temperate fishes most commonly found in the north Pacific. As part of the California Conservation Genomics Project (CCGP), we sequenced the genome of the Monkeyface Prickleback, Cebidichthys violaceus, to establish a genomic model for understanding phylogeographic patterns of marine organisms in California. These patterns, in turn, may inform the design of marine protected areas using dispersal models based on forthcoming population genomic data. The genome of C. violaceus is typical of many marine fishes at less than 1 Gb (genome size = 575.6 Mb), and our assembly is near-chromosome level (contig N50 = 1 Mb, scaffold N50 = 16.4 Mb, BUSCO completeness = 93.2%). Within the context of the CCGP, the genome will be used as a reference for future whole genome resequencing projects, enhancing our knowledge of the population structure of the species and more generally, the efficacy of marine protected areas as a primary conservation tool across California's marine ecosystems.

Reference Genome of the Black Surfperch, Embiotoca jacksoni (Embiotocidae, Perciformes), a California Kelp Forest Fish That Lacks a Pelagic Larval Stage.

Surfperches (Family Embiotocidae) are viviparous temperate reef fishes that brood their young. This life history trait translates into limited dispersal, strong population structure, and an unusually strong potential for local adaptation in a marine fish. As part of the California Conservation Genomics Project (CCGP), we sequenced the genome of the Black Surfperch, Embiotoca jacksoni, to establish a genomic model for understanding phylogeographic patterns of marine organisms in California. These patterns, in turn, may inform the design of marine protected areas using dispersal models based on genomic data. The genome of E. jacksoni is typical of marine fishes at less than 1Gb (genome size = 635 Mb), and our assembly is near-chromosome level (contig N50 = 6.5Mb, scaffold N50 = 15.5 Mb, BUSCO = 98.1%). Within the context of the CCGP, the genome will be used as a reference for future whole genome resequencing projects aimed at enhancing our knowledge of the population structure of the species, and efficacy of Marine Protected Areas across the state.

Reference Genome of the California Sheephead, Semicossyphus pulcher (Labridae, Perciformes), A Keystone Fish Predator in Kelp Forest Ecosystems.

Keystone species are known to play a critical role in kelp forest health, including the well-known killer whales, sea otter, sea urchin, kelp trophic cascade in the Aleutian Islands, Alaska, USA. In California, a major player in the regulation of sea urchin abundance, and in turn, the health of kelp forests ecosystems, is a large wrasse, the California Sheephead, Semicossyphus pulcher. We present a reference genome for this ecologically important species that will serve as a key resource for future conservation research of California's inshore marine environment utilizing genomic tools to address changes in life-history traits, dispersal, range shifts, and ecological interactions among members of the kelp forest ecological assemblages. Our genome assembly of S. pulcher has a total length of 0.794 Gb, which is similar to many other marine fishes. The assembly is largely contiguous (N50 = 31.9 Mb) and nearly complete (BUSCO single-copy core gene content = 98.1%). Within the context of the California Conservation Genomics Project (CCGP), the genome of S. pulcher will be used as an important reference resource for ongoing whole genome resequencing efforts of the species.

Upwelling-level acidification and pH/pCO2 variability moderate effects of ocean acidification on brain gene expression in the temperate surfperch, Embiotoca jacksoni.

Acidification-induced changes in neurological function have been documented in several tropical marine fishes. Here, we investigate whether similar patterns of neurological impacts are observed in a temperate Pacific fish that naturally experiences regular and often large shifts in environmental pH/pCO2 . In two laboratory experiments, we tested the effect of acidification, as well as pH/pCO2 variability, on gene expression in the brain tissue of a common temperate kelp forest/estuarine fish, Embiotoca jacksoni. Experiment 1 employed static pH treatments (target pH = 7.85/7.30), while Experiment 2 incorporated two variable treatments that oscillated around corresponding static treatments with the same mean (target pH = 7.85/7.70) in an eight-day cycle (amplitude ± 0.15). We found that patterns of global gene expression differed across pH level treatments. Additionally, we identified differential expression of specific genes and enrichment of specific gene sets (GSEA) in comparisons of static pH treatments and in comparisons of static and variable pH treatments of the same mean pH. Importantly, we found that pH/pCO2 variability decreased the number of differentially expressed genes detected between high and low pH treatments, and that interindividual variability in gene expression was greater in variable treatments than static treatments. These results provide important confirmation of neurological impacts of acidification in a temperate fish species and, critically, that natural environmental variability may mediate the impacts of ocean acidification.

The skeleton of Balanophyllia coral species suggests adaptive traits linked to the onset of mixotrophy.

The diversity in the skeletal features of coral species is an outcome of their evolution, distribution and habitat. Here, we explored, from macro- to nano-scale, the skeletal structural and compositional characteristics of three coral species belonging to the genus Balanophyllia having different trophic strategies. The goal is to address whether the onset of mixotrophy influenced the skeletal features of B. elegans, B. regia, and B. europaea. The macroscale data suggest that the presence of symbiotic algae in B. europaea can lead to a surplus of energy input that increases its growth rate and skeletal bulk density, leading to larger and denser corals compared to the azooxanthellate ones, B. regia and B. elegans. The symbiosis would also explain the higher intra-skeletal organic matrix (OM) content, which is constituted by macromolecules promoting the calcification, in B. europaea compared to the azooxanthellate species. The characterization of the soluble OM also revealed differences between B. europaea and the azooxanthellate species, which may be linked to diverse macromolecular machineries responsible for skeletal biosynthesis and final morphology. Differently, the crystallographic features were homogenous among species, suggesting that the basic building blocks of skeletons remained a conserved trait in these related species, regardless of the trophic strategy. These results show changes in skeletal phenotype that could be triggered by the onset of mixotrophy, as a consequence of the symbiotic association, displaying remarkable plasticity of coral skeletons which repeatedly allowed this coral group to adapt to a range of changing environments throughout its geological history.

Haplotype network branch diversity, a new metric combining genetic and topological diversity to compare the complexity of haplotype networks.

A common way of illustrating phylogeographic results is through the use of haplotype networks. While these networks help to visualize relationships between individuals, populations, and species, evolutionary studies often only quantitatively analyze genetic diversity among haplotypes and ignore other network properties. Here, we present a new metric, haplotype network branch diversity (HBd), as an easy way to quantifiably compare haplotype network complexity. Our metric builds off the logic of combining genetic and topological diversity to estimate complexity previously used by the published metric haplotype network diversity (HNd). However, unlike HNd which uses a combination of network features to produce complexity values that cannot be defined in probabilistic terms, thereby obscuring the values' implication for a sampled population, HBd uses frequencies of haplotype classes to incorporate topological information of networks, keeping the focus on the population and providing easy-to-interpret probabilistic values for randomly sampled individuals. The goal of this study is to introduce this more intuitive metric and provide an R script that allows researchers to calculate diversity and complexity indices from haplotype networks. A group of datasets, generated manually (model dataset) and based on published data (empirical dataset), were used to illustrate the behavior of HBd and both of its terms, haplotype diversity, and a new index called branch diversity. Results followed a predicted trend in both model and empirical datasets, from low metric values in simple networks to high values in complex networks. In short, the new combined metric joins genetic and topological diversity of haplotype networks, into a single complexity value. Based on our analysis, we recommend the use of HBd, as it makes direct comparisons of network complexity straightforward and provides probabilistic values that can readily discriminate situations that are difficult to resolve with available metrics.

Evolutionary origin of the Atlantic Cabo Verde nibbler (Girella stuebeli), a member of a primarily Pacific Ocean family of antitropical herbivorous reef fishes.

Nibblers (family Girellidae) are reef fishes that are mostly distributed in the Indo-Pacific, with one exception: Girella stuebeli, which is found in the Cabo Verde Archipelago, in the Atlantic Ocean. We capitalized on this unusual distribution to study the evolutionary history of the girellids, and determine the relationship between G. stuebeli and the remaining nibbler taxa. Based on thousands of genomic markers (RAD sequences), we identified the closest relatives of G. stuebeli as being a clade of three species endemic to the northwestern Pacific, restricted to the Sea of Japan and vicinity. This clade diverged from G. stuebeli approximately 2.2 Mya. Two alternative potential routes of migration may explain this affinity: a western route, from the Tropical Eastern Pacific and the Tropical Western Atlantic, and an eastern route via the Indian Ocean and Southern Africa. The geological history and oceanography of the regions combined with molecular data presented here, suggest that the eastern route of invasion (via the Indian Ocean and Southern Africa) is a more likely scenario.

Distinct patterns of hybridization across a suture zone in a coral reef fish (Dascyllus trimaculatus).

Hybrid zones are natural laboratories for investigating the dynamics of gene flow, reproductive isolation, and speciation. A predominant marine hybrid (or suture) zone encompasses Christmas Island (CHR) and Cocos (Keeling) Islands (CKE), where 15 different instances of interbreeding between closely related species from Indian and Pacific Oceans have been documented. Here, we report a case of hybridization between genetically differentiated Pacific and Indian Ocean lineages of the three-spot dascyllus, Dascyllus trimaculatus (Rüppell, 1829). Field observations indicate there are subtle color differences between Pacific and Indian Ocean lineages. Most importantly, population densities of color morphs and genetic analyses (mitochondrial DNA and SNPs obtained via RADSeq) suggest that the pattern of hybridization within the suture zone is not homogeneous. At CHR, both color morphs were present, mitochondrial haplotypes of both lineages were observed, and SNP analyses revealed both pure and hybrid genotypes. Meanwhile, in CKE, the Indian Ocean color morphs were prevalent, only Indian Ocean mitochondrial haplotypes were observed, and SNP analysis showed hybrid individuals with a large proportion (~80%) of their genotypes assigning to the Indian Ocean lineage. We conclude that CHR populations are currently receiving an influx of individuals from both ocean basins, with a greater influence from the Pacific Ocean. In contrast, geographically isolated CKE populations appear to be self-recruiting and with more influx of individuals from the Indian Ocean. Our research highlights how patterns of hybridization can be different at scales of hundreds of kilometers, due to geographic isolation and the history of interbreeding between lineages.

Patterns of Genomic Divergence and Signals of Selection in Sympatric and Allopatric Northeastern Pacific and Sea of Cortez Populations of the Sargo (Anisotremus davidsonii) and Longjaw Mudsucker (Gillichthys mirabilis).

Studying how isolation can impact population divergence and adaptation in co-distributed species can bring us closer to understanding how landscapes affect biodiversity. The Sargo, Anisotremus davidsonii (Haemulidae), and the Longjaw mudsucker, Gillichthys mirabilis (Gobiidae), offer a notable framework to study such mechanisms as their Pacific populations cross phylogeographic breaks at Point Conception, California, United States, and Punta Eugenia, Mexico, and are separated to those in the Sea of Cortez by the Baja California peninsula. Here, thousands of loci are genotyped from 48 Sargos and 73 mudsuckers using RADseq to characterize overall genomic divergence, and search for common patterns of putatively neutral and non-neutral structure based on outlier loci among populations with hypothesized different levels of isolation. We further search for parallels between population divergence and the total proportion of outliers, outlier FST distribution, and the proportion of outliers matching coding regions in GenBank. Statistically significant differentiation is seen across Point Conception in mudsucker (FST = 0.15), Punta Eugenia in Sargo (FST = 0.02), and on either side of the Baja California peninsula in both species (FST = 0.11 and 0.23, in Sargo and mudsucker, respectively). Each species shows structure using neutral and non-neutral loci. Finally, higher population divergence yields a more even distribution of outliers along their differentiation range but does not always translate into higher outlier proportions or higher rates in which outliers are matched to coding regions. If repeated in similar systems, observed genomic patterns might reveal speciation signatures in diverse networks of population isolation.

Phylogeography of the banded butterflyfish, Chaetodon striatus, indicates high connectivity between biogeographic provinces and ecosystems in the western Atlantic

Among the four butterflyfishes of the genus Chaetodon present in the western Atlantic, the banded butterflyfish Chaetodon striatus has the largest distribution range, spanning 44 degrees of latitude (from Massachusetts, USA to Santa Catarina, Brazil). Although the ecology of the banded butterflyfish has been well studied over its entire range, nothing is known about its phylogeography and how biogeographic barriers structure its populations. To assess the level of genetic connectivity among populations from distinct biogeographic provinces and environmental conditions, we collected samples from seven localities: Puerto Rico, in the Caribbean, and Tamandaré, Salvador, Abrolhos, Trindade Island, Arraial do Cabo and Florianópolis, in Brazil. One nuclear (rag 2) and two mitochondrial (control region and cyt b) molecular markers were sequenced. Our findings are consistent with a recent population expansion, around 30-120 thousand years ago, which was found for all populations. Haplotype network analyses point to the Caribbean as a refugium before the population expansion. Results show no geographic pattern of genetic diversity. Indeed, a lack of population structure was found and no isolation was observed across oceanographic barriers, as well as between coral and rocky reef ecosystems. Furthermore, no directionality in the migration pattern was found among populations. Since ecological and environmental characteristics are very diverse across such a vast geographic range, the lack of genetic differentiation suggests that C. striatus evolved ecological plasticity rather than local adaptation in the western Atlantic.(AU)

O peixe-borboleta listrado, Chaetodon striatus, possui a maior distribuição geográfica dentre as quatro espécies de peixes-borboleta do gênero Chaetodon presentes no Oceano Atlântico Ocidental, abrangendo 44° de latitude (entre Massachusetts, EUA até o sul do Brasil). A ecologia alimentar desta espécie é bastante conhecida, considerando a ampla distribuição, porém, pouco se sabe sobre a filogeografia e como as barreiras biogeográficas estruturam as populações. Para acessar a conectividade genética entre as populações de diferentes províncias biogeográficas e diferentes condições ambientais, foram coletadas amostras de sete localidades: Porto Rico, no Caribe, e Tamandaré, Salvador, Abrolhos, Ilha da Trindade, Arraial do Cabo e Florianópolis, no Brasil. Foram sequenciados um gene nuclear (rag 2) e dois genes mitocondriais (região controle e cit B). Para todas as populações, foi identificada uma expansão populacional recente, em torno de 30-120 mil anos atrás. A análise de rede de haplótipos sugere que o Caribe serviu como refúgio antes desta expansão populacional. Os resultados indicam que não há padrão geográfico de diversidade genética. Apesar da existência de barreiras oceanográficas e diferenças na constituição dos recifes (rochosos e coralíneos), não foi encontrada estruturação populacional. Também, não encontramos padrão na direção de migração entre as populações. Os resultados sugerem que C. striatus apresenta plasticidade ecológica, uma vez que não há diferenciação genética entre as populações que habitam ecossistemas tão diferentes ao longo da ampla distribuição no Atlântico Ocidental.(AU)

Alloparental care in the sea: Brood parasitism and adoption within and between two species of coral reef Altrichthys damselfish?

The evolution of parental care opens the door for the evolution of brood parasitic strategies that allow individuals to gain the benefits of parental care without paying the costs. Here we provide the first documentation for alloparental care in coral reef fish and we discuss why these patterns may reflect conspecific and interspecific brood parasitism. Species-specific barcodes revealed the existence of low levels (3.5% of all offspring) of mixed interspecific broods, mostly juvenile Amblyglyphidodon batunai and Pomacentrus smithi damselfish in Altrichthys broods. A separate analysis of conspecific parentage based on microsatellite markers revealed that mixed parentage broods are common in both species, and the genetic patterns are consistent with two different modes of conspecific brood parasitism, although further studies are required to determine the specific mechanisms responsible for these mixed parentage broods. While many broods had offspring from multiple parasites, in many cases a given brood contained only a single foreign offspring, perhaps a consequence of the movement of lone juveniles between nests. In other cases, broods contained large numbers of putative parasitic offspring from the same parents and we propose that these are more likely to be cases where parasitic adults laid a large number of eggs in the host nest than the result of movements of large numbers of offspring from a single brood after hatching. The evidence that these genetic patterns reflect adaptive brood parasitism, as well as possible costs and benefits of parasitism to hosts and parasites, are discussed.

Genetic diversity mirrors trophic ecology in coral reef fish feeding guilds.

Genetic diversity is essential for species persistence because it provides the raw material for evolution. For marine organisms, short pelagic larval duration (PLD) and small population size are characteristics generally assumed to associate with low genetic diversity. The ecological diversity of organisms may also affect genetic diversity with an expected corollary that more restricted habitat and dietary requirements could lead to a reduced genetic diversity because of pronounced genetic structuring. Here, we tested whether groups of species with narrower trophic niches displayed lower genetic diversity than those with broader niches. In order to test those predictions, we used different trophic guilds (i.e., groups of species having similar trophic habits) of coral reef damselfishes in Moorea (French Polynesia) for which we determined their genetic diversity using restriction site-associated DNA sequencing (RADseq) and their trophic ecology with stomach contents and stable isotope data. We found that pelagic feeders- the guild picking zooplankton in the water column- exhibited the lowest genetic diversity despite having the longest PLD and the largest population size. This guild had also the lowest variation in habitat characteristics and dietary composition compared to benthic feeders (i.e., those mainly grazing on algae) and the intermediate group (i.e., those feeding on zooplankton, filamentous algae and small benthic invertebrates). Our findings highlight the association between trophic ecology and genetic diversity that should be more commonly investigated in population genetics.

Taxonomic revisions within Embiotocidae (Teleostei, Perciformes) based on molecular phylogenetics.

Embiotocidae, a unique family within the Perciformes that has evolved a complex viviparous natural history, has lacked full resolution and strong support in several interspecific relationships until recently. Here we propose three taxonomic revisions within embiotocid surfperches based on recent molecular phylogenetic analyses that robustly resolve all interspecific relationship in the Eastern Pacific species: Hypsurus caryi (Agassiz, 1853) resurrected to its original name Embiotoca caryi Agassiz, 1853, Rhacochilus vacca (Girard, 1855) shifted into the genus Phanerodon Girard, 1854, and Hyperprosopon anale Agassiz, 1861 separated into the available genus Hypocritichthys Gill, 1862. The proposed changes would leave three previously paraphyletic groups monophyletic (Embiotoca, Hyperprosopon, and Phanerodon) and would maintain the current number of genera at 13.

Genomic islands of divergence in the Yellow Tang and the Brushtail Tang Surgeonfishes.

The current ease of obtaining thousands of molecular markers challenges the notion that full phylogenetic concordance, as proposed by phylogenetic species concepts, is a requirement for defining species delimitations. Indeed, the presence of genomic islands of divergence, which may be the cause, or in some cases the consequence, of speciation, precludes concordance. Here, we explore this issue using thousands of RAD markers on two sister species of surgeonfishes (Teleostei: Acanthuridae), Zebrasoma flavescens and Z. scopas, and several populations within each species. Species are readily distinguished based on their colors (solid yellow and solid brown, respectively), yet populations and species are neither distinguishable using mitochondrial markers (cytochrome c oxidase 1), nor using 5193 SNPs (pairwise Φst = 0.034). In contrast, when using outlier loci, some of them presumably under selection, species delimitations, and strong population structure follow recognized taxonomic positions (pairwise Φst = 0.326). Species and population delimitation differences based on neutral and selected markers are likely due to local adaptation, thus being consistent with the idea that these genomic islands of divergence arose as a consequence of isolation. These findings, which are not unique, raise the question of a potentially important pathway of divergence based on local adaptation that is only evident when looking at thousands of loci.