Speciation in fishes

Historical biogeography of the Mugil cephalus species complex and its rapid global colonization

Our understanding of speciation processes in marine environments remains very limited and the role of different reproductive barriers are still debated. While physical barriers were considered important drivers causing reproductive isolation, recent studies highlight the importance of climatic and hydrological changes creating unsuitable habitat conditions as factors promoting population isolation. Although speciation in marine fishes has been investigated from different perspectives, these studies often have a limited geographical extant. Therefore, data on speciation within widely distributed species are largely lacking. Species complexes offer valuable opportunities to study the initial stages of speciation. Herein we study speciation within the Mugil cephalus species complex (MCSC) which presents a unique opportunity due to its circumglobal distribution. We used a whole-genome shotgun analysis approach to identify SNPs among the 16 species within the MCSC. We inferred the phylogenetic relationships within the species complex followed by a time-calibration analysis. Subsequently, we estimated the ancestral ranges within the species complex to explore their biogeographical history. Herein, we present a fully resolved and well-supported phylogeny of the MCSC. Its origin is dated at around 3.79 Ma after which two main clades emerged: one comprising all West Atlantic and East Pacific species and the other all East Atlantic and Indo-Pacific species. Rapid dispersal following an initial founder colonization from the West to the East Atlantic led to the population of all major realms worldwide in less than 2 Myr. Physical and climatic barriers heavily impacted the ancestral distribution ranges within the MCSC and triggered the onset of speciation.

Morphological evidence supports splitting of species in the North Atlantic Sebastes spp. complex

The redfishes (genus Sebastes) are long-lived, commercial species in the North Atlantic. Excessive harvest through decades has led to a decline in the mature population of golden redfish (Sebastes norvegicus) in Norwegian waters, which is currently considered severely depleted. Accumulating genetic evidence suggests a more complex structure within this genus in the North Atlantic, which has recently inspired the hypotheses of cryptic species within S. norvegicus. Despite apparent genetic divergence between two types, they have yet to be verified morphologically. The morphology of genetically assigned fishes from Norwegian and Greenland waters was investigated using traditional morphometric methods, applying Linear Discriminant Analysis and Random Forest classification procedures to identify and evaluate the performance of descriptive characters. Combined with non-parametric meristic analysis, the results show that features such as beak length and eye diameter provide sufficient discrimination between the proposed cryptic species as well as separating them from the sympatric species S. mentella and S. viviparus. These findings support the presence of an additional redfish species in the North Atlantic, distinguishable both by morphological and genetic characters. This needs to be taken into consideration in future monitoring and management strategies for North Atlantic redfish.

A new discovery of cryptic species in barred sand burrower Limnichthys fasciatus (Pisces: Creediidae) from the western Pacific, with evolutionary perspectives of anti-equatorial species

Although the barred sand burrower, Limnichthys fasciatus, is widely distributed throughout the western Pacific, including Japan, Taiwan, and Australia, its morphology and genetics are poorly known. We discovered four cryptic species of Limnichthys from the western Pacific based on mtDNA cytochrome c oxidase subunit I (COI) and 16S ribosomal RNA (rRNA) sequences. Genetic distances showed remarkably large differences (12.7%-24.3% in COI and 7.9%-19.6% in 16S rRNA) between true L. fasciatus (type locality: southeastern Australia) and the others. A relaxed clock model with optimized selected substitution models showed that their deep divergence began in the middle Miocene epoch and subsequently diverged into the current cryptic species in the Plio-Pleistocene. A eurythermal common ancestor may have evolved independently in each region due to geographical events and paleoclimatic fluctuations, which made it possible for L. fasciatus complex to be an anti-equatorial species. Despite their deep genetic divergence, they showed marked phenotypic similarity, suggesting that they have experienced similar selective pressures related to their specific behavior.

Genome-wide species delimitation analyses of a silverside fish species complex in central Mexico indicate taxonomic over-splitting

BACKGROUND

Delimiting species across a speciation continuum is a complex task, as the process of species origin is not generally instantaneous. The use of genome-wide data provides unprecedented resolution to address convoluted species delimitation cases, often unraveling cryptic diversity. However, because genome-wide approaches based on the multispecies coalescent model are known to confound population structure with species boundaries, often resulting in taxonomic over-splitting, it has become increasingly evident that species delimitation research must consider multiple lines of evidence. In this study, we used phylogenomic, population genomic, and coalescent-based species delimitation approaches, and examined those in light of morphological and ecological information, to investigate species numbers and boundaries comprising the Chirostoma "humboltianum group" (family Atherinidae). The humboltianum group is a taxonomically controversial species complex where previous morphological and mitochondrial studies produced conflicting species delimitation outcomes. We generated ddRADseq data for 77 individuals representing the nine nominal species in the group, spanning their distribution range in the central Mexican plateau.

RESULTS

Our results conflict with the morphospecies and ecological delimitation hypotheses, identifying four independently evolving lineages organized in three geographically cohesive clades: (i) chapalae and sphyraena groups in Lake Chapala, (ii) estor group in Lakes Pátzcuaro and Zirahuén, and (iii) humboltianum sensu stricto group in Lake Zacapu and Lerma river system.

CONCLUSIONS

Overall, our study provides an atypical example where genome-wide analyses delineate fewer species than previously recognized on the basis of morphology. It also highlights the influence of the geological history of the Chapala-Lerma hydrological system in driving allopatric speciation in the humboltianum group.

Speciation genomics and the role of depth in the divergence of rockfishes (Sebastes) revealed through Pool-seq analysis of enriched sequences

Speciation in the marine environment is challenged by the wide geographic distribution of many taxa and potential for high rates of gene flow through larval dispersal mechanisms. Depth has recently been proposed as a potential driver of ecological divergence in fishes, and yet it is unclear how adaptation along these gradients' shapes genomic divergence. The genus Sebastes contains numerous species pairs that are depth-segregated and can provide a better understanding of the mode and tempo of genomic diversification. Here, we present exome data on two species pairs of rockfishes that are depth-segregated and have different degrees of divergence: S. chlorostictus-S. rosenblatti and S. crocotulus-S. miniatus. We were able to reliably identify "islands of divergence" in the species pair with more recent divergence (S. chlorostictus-S. rosenblatti) and discovered a number of genes associated with neurosensory function, suggesting a role for this pathway in the early speciation process. We also reconstructed demographic histories of divergence and found the best supported model was isolation followed by asymmetric secondary contact for both species pairs. These results suggest past ecological/geographic isolation followed by asymmetric secondary contact of deep to shallow species. Our results provide another example of using rockfish as a model for studying speciation and support the role of depth as an important mechanism for diversification in the marine environment.

Neutral and outlier single nucleotide polymorphisms disentangle the evolutionary history of a coastal Solanaceae species

Speciation begins with the isolation of some individuals or subpopulations due to drivers promoting a diverging genetic distribution. Such isolation may occur, followed by different processes and pressures. Isolation-by-distance (IBD), isolation-by-adaptation (IBA), and isolation-by-colonization (IBC) have been recognized as the main divergence patterns. Still, it is not easy to distinguish which one is the main pattern as each one may act at different points in time or even simultaneously. Using an extensive genome coverage from a Petunia species complex with coastal and inland distribution and multiple analytical approaches on population genomics and phylogeography, we showed a complex interplay between neutral and selective forces acting on the divergence process. We found 18,887 SNPs potentially neutral and 924 potentially under selection (outlier) loci. All analyses pointed that each subspecies displays its own genetic component and evolutionary history. We suggested plausible ecologic drivers for such divergence in a southernmost South Atlantic coastal plain in Brazil and Uruguay and identified a connection between adaptation and environment heterogeneity.

Genetic relationships between sympatric and allopatric Coregonus ciscoes in North and Central Europe

BACKGROUND

Sympatric speciation along ecological gradients has been studied repeatedly, in particular in freshwater fishes. Rapid post-glacial ecological divergence has resulted in numerous endemic species or ecologically distinct populations in lakes of the temperate zones. Here, we focus on the Baltic cisco (Coregonus albula) complex, to study the genetic similarity among two pairs of sympatric autumn- and spring-spawning populations from post-glacial German Lakes Stechlin and Breiter Luzin. For comparison, we included a similar pair of sympatric populations from the Swedish Lake Fegen. We wanted to explore potential genetic similarities between the three sympatric cisco population pairs in the three lakes, to evaluate whether the pairs may have emerged independently in the three lakes, or whether two different species may have colonized all three lakes independently. Furthermore, we considered allopatric C. albula populations from three Polish, three Finnish, and four Swedish locations, added one Siberian population of the sister species C. sardinella and a Swedish C. maraena (whitefish) population. By genotyping nine microsatellite markers in 655 individuals from these 18 populations, we wanted to elucidate how strongly the cisco populations differ across a larger geographical area within Europe. Finally, we compared the genetic differences between the spring- and autumn-spawning populations of ciscoes in the two German lakes to infer the potentially deteriorating effect of strong anthropogenic pressure on the lakes.

RESULTS

Dendrogram, Principal Coordinate Analysis and admixture analysis all indicated strong correspondence between population differentiation and geographical location for most cisco populations in Europe, including the Siberian population of C. sardinella. However, populations from some Swedish lakes deviated from this general pattern, by showing a distinct genetic structure. We found evidence for independent evolution of the three sympatric population pairs, because the populations co-occurring in the same lake were always most closely related. However, genetic differentiation was weak in the two German population pairs, but strong in the Swedish Lake Fegen, indicating that the weak differentiation in the German pairs reported earlier has eroded further.

CONCLUSIONS

Our results suggest that the genetic differentiation at neutral genetic markers among populations of the Baltic cisco complex has evolved (and is maintained) by random genetic drift in isolated populations. However, earlier studies on the Swedish populations combining mitochondrial DNA and microsatellite data indicate that also post-glacial immigration from separate glacial refugia has shaped the present genetic population structure. The low neutral differentiation of the German sympatric pairs in contrast to the Swedish pair suggests that recent anthropogenic effects on the lakes in Germany may put the endemic spring-spawners at risk to extinction.

Genomic data and multi-species demographic modelling uncover past hybridization between currently allopatric freshwater species

Evidence for ancient interspecific gene flow through hybridization has been reported in many animal and plant taxa based on genetic markers. The study of genomic patterns of closely related species with allopatric distributions allows the assessment of the relative importance of vicariant isolating events and past gene flow. Here, we investigated the role of gene flow in the evolutionary history of four closely related freshwater fish species with currently allopatric distributions in western Iberian rivers-Squalius carolitertii, S. pyrenaicus, S. torgalensis and S. aradensis-using a population genomics dataset of 23,562 SNPs from 48 individuals, obtained through genotyping by sequencing (GBS). We uncovered a species tree with two well-differentiated clades: (i) S. carolitertii and S. pyrenaicus; and (ii) S. torgalensis and S. aradensis. By using D-statistics and demographic modelling based on the site frequency spectrum, comparing alternative demographic scenarios of hybrid origin, secondary contact and isolation, we found that the S. pyrenaicus North lineage is likely the result of an ancient hybridization event between S. carolitertii (contributing ~84%) and S. pyrenaicus South lineage (contributing ~16%), consistent with a hybrid speciation scenario. Furthermore, in the hybrid lineage, we identify outlier loci potentially affected by selection favouring genes from each parental lineage at different genomic regions. Our results suggest that ancient hybridization can affect speciation and that freshwater fish species currently in allopatry are useful to study these processes.

Adding New Pieces to the Puzzle of Karyotype Evolution in Harttia (Siluriformes, Loricariidae): Investigation of Amazonian Species

A remarkable morphological diversity and karyotype variability can be observed in the Neotropical armored catfish genus Harttia. These fishes offer a useful model to explore both the evolution of karyotypes and sex chromosomes, since many species possess male-heterogametic sex chromosome systems and a high rate of karyotype repatterning. Based on the karyotype organization, the chromosomal distribution of several repetitive DNA classes, and the rough estimates of genomic divergences at the intraspecific and interspecific levels via Comparative Genomic Hybridization, we identified shared diploid chromosome numbers (2n = 54) but different karyotype compositions in H. dissidens (20m + 26sm + 8a) and Harttia sp. 3 (16m + 18sm + 14st + 6a), and different 2n in H. guianensis (2n = 58; 20m + 26sm + 2st + 10a). All species further displayed similar patterns of chromosomal distribution concerning constitutive heterochromatin, 18S ribosomal DNA (rDNA) sites, and most of the surveyed microsatellite motifs. Furthermore, differences in the distribution of 5S rDNA sites and a subset of microsatellite sequences were identified. Heteromorphic sex chromosomes were lacking in H. dissidens and H. guianensis at the scale of our analysis. However, one single chromosome pair in Harttia sp. 3 males presented a remarkable accumulation of male genome-derived probe after CGH, pointing to a tentative region of early sex chromosome differentiation. Thus, our data support already previously outlined evidence that Harttia is a vital model for the investigation of teleost karyotype and sex chromosome dynamics.

A Review on Fish Sensory Systems and Amazon Water Types With Implications to Biodiversity

The Amazon has the highest richness of freshwater organisms in the world, which has led to a multitude of hypotheses on the mechanisms that generated this biodiversity. However, most of these hypotheses focus on the spatial distance of populations, a framework that fails to provide an explicit mechanism of speciation. Ecological conditions in Amazon freshwaters can be strikingly distinct, as it has been recognized since Alfred Russel Wallace’s categorization into black, white, and blue (= clear) waters. Water types reflect differences in turbidity, dissolved organic matter, electrical conductivity, pH, amount of nutrients and lighting environment, characteristics that directly affect the sensory abilities of aquatic organisms. Since natural selection drives evolution of sensory systems to function optimally according to environmental conditions, the sensory systems of Amazon freshwater organisms are expected to vary according to their environment. When differences in sensory systems affect chances of interbreeding between populations, local adaptations may result in speciation. Here, we briefly present the limnologic characteristics of Amazonian water types and how they are expected to influence photo-, chemical-, mechano-, and electro-reception of aquatic organisms, focusing on fish. We put forward that the effect of different water types on the adaptation of sensory systems is an important mechanism that contributed to the evolution of fish diversity. We point toward underexplored research perspectives on how divergent selection may act on sensory systems and thus contribute to the origin and maintenance of the biodiversity of Amazon aquatic environments.

Hydrography rather than lip morphology better explains the evolutionary relationship between Gymnogeophagus labiatus and G. lacustris in Southern Brazil (Cichlidae: Geophagini)

ABSTRACT Gymnogeophagus labiatus and G. lacustris have been long recognized as sister species exhibiting different ecological requirements. Gymnogeophagus labiatus occurs in rock bottom rivers in the hydrographic basins of Patos Lagoon (HBP) and Tramandaí River (HBT), while G. lacustris is exclusive from sand bottom coastal lagoons of the HBT. In this study, we used molecular markers, morphological measurements and data from nuptial male coloration to investigate the evolutionary relationship between these species in each hydrographic basin. We found, for all data sets, a closer relationship between G. labiatus and G. lacustris from the HBT than between G. labiatus populations from HBT and HBP. In particular, lip area had a large intraspecific plasticity, being uninformative to diagnose G. lacustris from G. labiatus. Molecular clock-based estimates suggest a recent divergence between species in the HBT (17,000 years ago), but not between G. labiatus from HBP and HBT (3.6 millions of years ago). Finally, we also found a divergent G. labiatus genetic lineage from the Camaquã River, in the HBP. These results show that the current taxonomy of G. labiatus and G. lacustris does not properly represent evolutionary lineages in these species.

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.

Speciation history of European (Anguilla anguilla) and American eel (A. rostrata), analysed using genomic data

Speciation in the ocean could differ from terrestrial environments due to fewer barriers to gene flow. Hence, sympatric speciation might be common, with American and European eel being candidates for exemplifying this. They show disjunct continental distributions on both sides of the Atlantic, but spawn in overlapping regions of the Sargasso Sea from where juveniles are advected to North American, European and North African coasts. Hybridization and introgression are known to occur, with hybrids almost exclusively observed in Iceland. Different speciation scenarios have been suggested, involving either vicariance or sympatric ecological speciation. Using RAD sequencing and whole-genome sequencing data from parental species and F1 hybrids, we analysed speciation history based on the joint allele frequency spectrum (JAFS) and pairwise sequentially Markovian coalescent (PSMC) plots. JAFS supported a model involving a split without gene flow 150,000-160,000 generations ago, followed by secondary contact 87,000-92,000 generations ago, with 64% of the genome experiencing restricted gene flow. This supports vicariance rather than sympatric speciation, likely associated with Pleistocene glaciation cycles and ocean current changes. Whole-genome PSMC analysis of F1 hybrids from Iceland suggested divergence 200,000 generations ago and indicated subsequent gene flow rather than strict isolation. Finally, simulations showed that results from both approaches (JAFS and PSMC) were congruent. Hence, there is strong evidence against sympatric speciation in North Atlantic eels. These results reiterate the need for careful consideration of cases of possible sympatric speciation, as even in seemingly barrier-free oceanic environments palaeoceanographic factors may have promoted vicariance and allopatric speciation.

Natural hybridization between divergent lineages in a selfing hermaphroditic fish

By definition, mating between individuals is infrequent in highly selfing organisms, and so too, therefore, hybridization should be rare between genetically divergent lineages in predominantly self-fertilizing species. Notwithstanding these expectations, here we report a remarkable case of natural hybridization between highly diverged phylogeographic lineages of the mangrove rivulus, a small killifish that reproduces predominantly by self-fertilization and typically is found as highly homozygous lines in most parts of its extensive geographical range. Two distinctive genetic lineages (Kryptolebias marmoratus and a 'Central clade' closely related to K. hermaphroditus) previously were not known in sympatry, but were found by us to co-occur on San Salvador, Bahamas. Genetic analyses of a mitochondrial and multiple nuclear markers determined the direction of a cross producing a hybrid fish. Furthermore, we show that this hybrid individual was viable, as it successfully reproduced by self-fertilization for two generations. Additional sampling of this population will be necessary to determine if backcrossing of hybrids to the parental lineages occurs in nature and to analyse whether such backcross progeny are viable. Application of the biological species concept (BSC) is traditionally difficult in clonally reproducing organisms. Our results show that although mangrove rivulus fish are mostly highly selfing in nature (resulting in isogenic, effectively clonal and homozygous progeny), classification within this taxonomic complex need not be incompatible with the BSC.

Genomic, ecological, and morphological approaches to investigating species limits: A case study in modern taxonomy from Tropical Eastern Pacific surgeonfishes

A wide variety of species are distinguished by slight color variations. However, molecular analyses have repeatedly demonstrated that coloration does not always correspond to distinct evolutionary histories between closely related groups, suggesting that this trait is labile and can be misleading for species identification. In the present study, we analyze the evolutionary history of sister species of Prionurus surgeonfishes in the Tropical Eastern Pacific (TEP), which are distinguished by the presence or absence of dark spots on their body. We examined the species limits in this system using comparative specimen-based approaches, a mitochondrial gene (COI), more than 800 nuclear loci (Ultraconserved Elements), and abiotic niche comparisons. The results indicate there is a complete overlap of meristic counts and morphometric measurements between the two species. Further, we detected multiple individuals with intermediate spotting patterns suggesting that coloration is not diagnostic. Mitochondrial data recovered a single main haplotype shared between the species and all locations resulting in a complete lack of structure (ΦST = 0). Genomic analyses also suggest low levels of genetic differentiation (F ST = 0.013), and no alternatively fixed SNPs were detected between the two phenotypes. Furthermore, niche comparisons could not reject niche equivalency or similarity between the species. These results suggest that these two phenotypes are conspecific and widely distributed in the TEP. Here, we recognize Prionurus punctatus Gill 1862 as a junior subjective synonym of P. laticlavius (Valenciennes 1846). The underlying causes of phenotypic variation in this species are unknown. However, this system gives insight into general evolutionary dynamics within the TEP.

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.

Flickers of speciation: Sympatric colour morphs of the arc-eye hawkfish, Paracirrhites arcatus, reveal key elements of divergence with gene flow

One of the primary challenges of evolutionary research is to identify ecological factors that favour reproductive isolation. Therefore, studying partially isolated taxa has the potential to provide novel insight into the mechanisms of evolutionary divergence. Our study utilizes an adaptive colour polymorphism in the arc-eye hawkfish (Paracirrhites arcatus) to explore the evolution of reproductive barriers in the absence of geographic isolation. Dark and light morphs are ecologically partitioned into basaltic and coral microhabitats a few metres apart. To test whether ecological barriers have reduced gene flow among dark and light phenotypes, we evaluated genetic variation at 30 microsatellite loci and a nuclear exon (Mc1r) associated with melanistic coloration. We report low, but significant microsatellite differentiation among colour morphs and stronger divergence in the coding region of Mc1r indicating signatures of selection. Critically, we observed greater genetic divergence between colour morphs on the same reefs than that between the same morphs in different geographic locations. We hypothesize that adaptation to the contrasting microhabitats is overriding gene flow and is responsible for the partial reproductive isolation observed between sympatric colour morphs. Combined with complementary studies of hawkfish ecology and behaviour, these genetic results indicate an ecological barrier to gene flow initiated by habitat selection and enhanced by assortative mating. Hence, the arc-eye hawkfish fulfil theoretical expectations for the earliest phase of speciation with gene flow.

Widespread hybridization and bidirectional introgression in sympatric species of coral reef fish

Coral reefs are highly diverse ecosystems, where numerous closely related species often coexist. How new species arise and are maintained in these high geneflow environments have been long-standing conundrums. Hybridization and patterns of introgression between sympatric species provide a unique insight into the mechanisms of speciation and the maintenance of species boundaries. In this study, we investigate the extent of hybridization between two closely related species of coral reef fish: the common coral trout (Plectropomus leopardus) and the bar-cheek coral trout (Plectropomus maculatus). Using a complementary set of 25 microsatellite loci, we distinguish pure genotype classes from first- and later-generation hybrids, identifying 124 interspecific hybrids from a collection of 2,991 coral trout sampled in inshore and mid-shelf reefs of the southern Great Barrier Reef. Hybrids were ubiquitous among reefs, fertile and spanned multiple generations suggesting both ecological and evolutionary processes are acting to maintain species barriers. We elaborate on these finding to investigate the extent of genomic introgression and admixture from 2,271 SNP loci recovered from a ddRAD library of pure and hybrid individuals. An analysis of genomic clines on recovered loci indicates that 261 SNP loci deviate from a model of neutral introgression, of which 132 indicate a pattern of introgression consistent with selection favouring both hybrid and parental genotypes. Our findings indicate genome-wide, bidirectional introgression between two sympatric species of coral reef fishes and provide further support to a growing body of evidence for the role of hybridization in the evolution of coral reef fishes.

DNA Barcoding Marine Biodiversity: Steps from Mere Cataloguing to Giving Reasons for Biological Differences

DNA barcoding has become a useful tool in many contexts and has opened up a completely new avenue for taxonomy. DNA barcoding has its widest application in biodiversity and ecological research to detect and describe diversity whenever morphological discrimination is difficult or impossible (e.g., in the case of species lacking diagnostic characters, early life stages, or cryptic species). In this chapter, we outline the utility of including physiological parameters as part of species description in publicly available databases that catalog taxonomic information resulting from barcoding projects. Cryptic species or different life stages of a species often differ in their physiological traits. Thus, if physiological aspects were included in species definitions, the presently cryptic species could be distinguished. We furthermore give suggestions for physiological information that should be included in a species description and describe potential applications of DNA barcoding for research with physiological components.

Cryptic lineage divergence in marine environments: genetic differentiation at multiple spatial and temporal scales in the widespread intertidal goby Gobiosoma bosc

The adaptive radiation of the seven-spined gobies (Gobiidae: Gobiosomatini) represents a classic example of how ecological specialization and larval retention can drive speciation through local adaptation. However, geographically widespread and phenotypically uniform species also do occur within Gobiosomatini. This lack of phenotypic variation across large geographic areas could be due to recent colonization, widespread gene flow, or stabilizing selection acting across environmental gradients. We use a phylogeographic approach to test these alternative hypotheses in the naked goby Gobiosoma bosc, a widespread and phenotypically invariable intertidal fish found along the Atlantic Coast of North America. Using DNA sequence from 218 individuals sampled at 15 localities, we document marked intraspecific genetic structure in mitochondrial and nuclear genes at three main geographic scales: (i) between Gulf of Mexico and Atlantic Coast, (ii) between the west coast of the Florida peninsula and adjacent Gulf of Mexico across the Apalachicola Bay, and (iii) at local scales of a few hundred kilometers. Clades on either side of Florida diverged about 8 million years ago, whereas some populations along the East Cost show divergent phylogroups that have differentiated within the last 200,000 years. The absence of noticeable phenotypic or ecological differentiation among lineages suggests the role of stabilizing selection on ancestral phenotypes, together with isolation in allopatry due to reduced dispersal and restricted gene flow, as the most likely explanation for their divergence. Haplotype phylogenies and spatial patterns of genetic diversity reveal frequent population bottlenecks followed by rapid population growth, particularly along the Gulf of Mexico. The magnitude of the genetic divergence among intraspecific lineages suggests the existence of cryptic species within Gobiosoma and indicates that modes of speciation can vary among lineages within Gobiidae.

Genetic Diversity and Hybridisation between Native and Introduced Salmonidae Fishes in a Swedish Alpine Lake

Understanding the processes underlying diversification can aid in formulating appropriate conservation management plans that help maintain the evolutionary potential of taxa, particularly under human-induced activities and climate change. Here we assessed the microsatellite genetic diversity and structure of three salmonid species, two native (Arctic charr, Salvelinus alpinus and brown trout, Salmo trutta) and one introduced (brook charr, Salvelinus fontinalis), from an alpine lake in sub-arctic Sweden, Lake Ånn. The genetic diversity of the three species was similar and sufficiently high from a conservation genetics perspective: corrected total heterozygosity, H’_T = 0.54, 0.66, 0.60 and allelic richness, A_R = 4.93, 5.53 and 5.26 for Arctic charr, brown trout and brook charr, respectively. There were indications of elevated inbreeding coefficients in brown trout (G_IS = 0.144) and brook charr (G_IS = 0.129) although sibling relationships were likely a confounding factor, as a high proportion of siblings were observed in all species within and among sampling locations. Overall genetic structure differed between species, Fst = 0.01, 0.02 and 0.04 in Arctic charr, brown trout and brook charr respectively, and there was differentiation at only a few specific locations. There was clear evidence of hybridisation between the native Arctic charr and the introduced brook charr, with 6% of individuals being hybrids, all of which were sampled in tributary streams. The ecological and evolutionary consequences of the observed hybridisation are priorities for further research and the conservation of the evolutionary potential of native salmonid species.

Phylogeny and historical demography of endemic fishes in Lake Biwa: the ancient lake as a promoter of evolution and diversification of freshwater fishes in western Japan

To elucidate the origins of the endemic fish of Lake Biwa, an ancient lake in Japan, and the role of the lake in the diversification of freshwater fish in western Japan, we established a molecular phylogenetic framework with an absolute time scale and inferred the historical demography of a large set of fish species in and around the lake. We used mtDNA sequences obtained from a total of 190 specimens, including 11 endemic species of Lake Biwa and their related species, for phylogenetic analyses with divergence time estimations and from a total of 2319 specimens of 42 species (including 14 endemics) occurring in the lake for population genetic analyses. Phylogenetic analysis suggested that some of the endemic species diverged from their closest relatives earlier (1.3–13.0 Ma) than the period in which the present environmental characteristics of the lake started to develop (ca. 0.4 Ma), whereas others diverged more recently (after 0.4 Ma). In contrast, historical demographic parameters suggested that almost all species, including endemic and nonendemic ones, expanded their populations after the development of the present lake environment. In phylogeographic analyses, common or very close haplotypes of some species were obtained from Lake Biwa and other regions of western Japan. The phylogenetic and historical demographic evidence suggests that there was a time lag between phylogenetic divergence and population establishment and that phenotypic adaptation of some endemic species to the limnetic environment occurred much later than the divergences of those endemic lineages. Population structure and phylogeographic patterns suggest that Lake Biwa has functioned not only as the center of adaptive evolution but also as a reservoir for fish diversity in western Japan.

Assessing models of speciation under different biogeographic scenarios; an empirical study using multi-locus and RNA-seq analyses

Evolutionary biology often seeks to decipher the drivers of speciation, and much debate persists over the relative importance of isolation and gene flow in the formation of new species. Genetic studies of closely related species can assess if gene flow was present during speciation, because signatures of past introgression often persist in the genome. We test hypotheses on which mechanisms of speciation drove diversity among three distinct lineages of desert tortoise in the genus Gopherus. These lineages offer a powerful system to study speciation, because different biogeographic patterns (physical vs. ecological segregation) are observed at opposing ends of their distributions. We use 82 samples collected from 38 sites, representing the entire species' distribution and generate sequence data for mtDNA and four nuclear loci. A multilocus phylogenetic analysis in *BEAST estimates the species tree. RNA‐seq data yield 20,126 synonymous variants from 7665 contigs from two individuals of each of the three lineages. Analyses of these data using the demographic inference package ∂a∂i serve to test the null hypothesis of no gene flow during divergence. The best‐fit demographic model for the three taxa is concordant with the *BEAST species tree, and the ∂a∂i analysis does not indicate gene flow among any of the three lineages during their divergence. These analyses suggest that divergence among the lineages occurred in the absence of gene flow and in this scenario the genetic signature of ecological isolation (parapatric model) cannot be differentiated from geographic isolation (allopatric model).

High levels of interspecific gene flow in an endemic cichlid fish adaptive radiation from an extreme lake environment

Studying recent adaptive radiations in isolated insular systems avoids complicating causal events and thus may offer clearer insight into mechanisms generating biological diversity. Here, we investigate evolutionary relationships and genomic differentiation within the recent radiation of Alcolapia cichlid fish that exhibit extensive phenotypic diversification, and which are confined to the extreme soda lakes Magadi and Natron in East Africa. We generated an extensive RAD data set of 96 individuals from multiple sampling sites and found evidence for genetic admixture between species within Lake Natron, with the highest levels of admixture between sympatric populations of the most recently diverged species. Despite considerable environmental separation, populations within Lake Natron do not exhibit isolation by distance, indicating panmixia within the lake, although individuals within lineages clustered by population in phylogenomic analysis. Our results indicate exceptionally low genetic differentiation across the radiation despite considerable phenotypic trophic variation, supporting previous findings from smaller data sets; however, with the increased power of densely sampled SNPs, we identify genomic peaks of differentiation (FST outliers) between Alcolapia species. While evidence of ongoing gene flow and interspecies hybridization in certain populations suggests that Alcolapia species are incompletely reproductively isolated, the identification of outlier SNPs under diversifying selection indicates the radiation is undergoing adaptive divergence.

High levels of genetic connectivity among populations of yellowtail snapper, Ocyurus chrysurus (Lutjanidae-Perciformes), in the western South Atlantic revealed through multilocus analysis

In the present study, five loci (mitochondrial and nuclear) were sequenced to determine the genetic diversity, population structure, and demographic history of populations of the yellowtail snapper, Ocyurus chrysurus, found along the coast of the western South Atlantic. O. chrysurus is a lutjanid species that is commonly associated with coral reefs and exhibits an ample geographic distribution, and it can therefore be considered a good model for the investigation of phylogeographic patterns and genetic connectivity in marine environments. The results reflected a marked congruence between the mitochondrial and nuclear markers as well as intense gene flow among the analyzed populations, which represent a single genetic stock along the entire coast of Brazil between the states of Pará and Espírito Santo. Our data also showed high levels of genetic diversity in the species (mainly mtDNA), as well a major historic population expansion, which most likely coincided with the sea level oscillations at the end of the Pleistocene. In addition, this species is intensively exploited by commercial fisheries, and data on the genetic structure of its populations will be essential for the development of effective conservation and management plans.

Population genomics of marine fishes: next-generation prospects and challenges

Over the past few years, technological advances have facilitated giant leaps forward in our ability to generate genome-wide molecular data, offering exciting opportunities for gaining new insights into the ecology and evolution of species where genomic information is still limited. Marine fishes are valuable organisms for advancing our understanding of evolution on historical and contemporary time scales, and here we highlight areas in which research on these species is likely to be particularly important in the near future. These include possibilities for gaining insights into processes on ecological time scales, identifying genomic signatures associated with population divergence under gene flow, and determining the genetic basis of phenotypic traits. We also consider future challenges pertaining to the implementation of genome-wide coverage through next-generation sequencing and genotyping methods in marine fishes. Complications associated with fast decay of linkage disequilibrium, as expected for species with large effective population sizes, and the possibility that adaptation is associated with both soft selective sweeps and polygenic selection, leaving complex genomic signatures in natural populations, are likely to challenge future studies. However, the combination of high genome coverage and new statistical developments offers promising solutions. Thus, the next generation of studies is likely to truly facilitate the transition from population genetics to population genomics in marine fishes. This transition will advance our understanding of basic evolutionary processes and will offer new possibilities for conservation and management of valuable marine resources.

Genetic structure of Scomber japonicus (Perciformes: Scombridae) along the coast of China revealed by complete mitochondrial cytochrome b sequences

The phylogeography history and contemporary agents of selection for many marine fisheries, characterized by widespread species distributions in the face of significant harvest, remains poorly understood. Chub mackerel (Scomber japonicus) are a widespread species in the Indo-Pacific and represent one of the top five commercially fished species in the world, yet their phylogeographic history remains unknown. We characterized the genetic diversity, structure and demographic history of S. japonicus throughout adjacent Chinese seas (from the Yellow Sea to the South China Sea). Using 220 individuals from 11 sites, we inferred 55 distinct haplotypes from complete mitochondrial cytochrome b gene sequences. Haplotype diversity ranged from 0.505 to 0.967 and nucleotide diversity ranged from 0.00056 to 0.01042. Genetic differentiation (F_st) statistics suggested that the highest level of differentiation existed between the SanYa and SanSha localities (F_st = 0.86977), while the lowest levels of differentiation occurred between the DongGang and ShiDao localities (F_st∼ 0). Kimura's genetic distances ranged from 0.001 to 0.011 within and from 0.001 to 0.018 between populations. Hierarchical analysis of molecular variance, Neighbor-joining and median-joining network analyses identified significant phylogeographic structure with two localities (SanYa, Hainan of the South China Sea and LianYunGang, Jiangsu of the East China Sea) explaining most of the genetic variation observed, while the remaining populations were poorly differentiated.

Combining morphology and genetics in resolving taxonomy--a systematic revision of spined loaches (Genus Cobitis; Cypriniformes, Actinopterygii) in the Adriatic watershed

Taxonomic investigation of spined loaches from Dalmatia and Herzegovina was conducted on specimens from 14 localities. The results of the detailed morphological investigations were combined with genetic data (based on one mitochondrial and two nuclear genes) in order to resolve the taxonomic status of each Cobitis population. Among the investigated features of external morphology, the appearance of spots on the caudal fin base turned out to have the greatest diagnostic value. Furthermore, the number of branched fin rays enabled the discrimination of several species. No morphometric character alone could ensure determination of any Cobitis species. Nevertheless, groups of populations that are more similar in their body shapes correspond to mitochondrial phylogenetic lineages. Based on molecular genetic markers, Dalmatian and Herzegovinian spined loaches form independent lineages inside the Adriatic phylogenetic group. Mitochondrial DNA phylogenetic reconstruction revealed six monophyletic lineages, corresponding to six species distributed in the investigated area. The population distributed in Mostarsko blato karstic field in Bosnia and Herzegovina is described as a new species based on a unique combination of morphological characters: a single triangular Canestrini scale; usually 5¹/₂ branched anal fin rays, 6¹/₂ branched dorsal fin rays, 14 branched caudal fin rays; no spots in the surface pigmentation layer on the caudal fin base; scales on the body very small.

Speciation and demographic history of Atlantic eels (Anguilla anguilla and A. rostrata) revealed by mitogenome sequencing

Processes leading to speciation in oceanic environments without obvious physical barriers remain poorly known. European and American eel (Anguilla anguilla and A. rostrata) spawn in partial sympatry in the Sargasso Sea. Larvae are advected by the Gulf Stream and other currents towards the European/North African and North American coasts, respectively. We analyzed 104 mitogenomes from the two species along with mitogenomes of other Anguilla and outgroup species. We estimated divergence time between the two species to identify major events involved in speciation. We also considered two previously stated hypotheses: one where the ancestral species was present in only one continent but was advected across the Atlantic by ocean current changes and another where population declines during Pleistocene glaciations led to increasing vicariance, facilitating speciation. Divergence time was estimated to ∼3.38 Mya, coinciding with the closure of the Panama Gateway that led to reinforcement of the Gulf Stream. This could have advected larvae towards European/North African coasts, in which case American eel would be expected to be the ancestral species. This scenario could, however, not be unequivocally confirmed by analyses of dN/dS, nucleotide diversity and effective population size estimates. Extended bayesian skyline plots showed fluctuations of effective population sizes and declines during glaciations, and thus also lending support to the importance of vicariance during speciation. There was evidence for positive selection at the ATP6 and possibly ND5 genes, indicating a role in speciation. The findings suggest an important role of ocean current changes in speciation of marine organisms.

Phylogeography of the California sheephead, Semicossyphus pulcher: the role of deep reefs as stepping stones and pathways to antitropicality

In the past decade, the study of dispersal of marine organisms has shifted from focusing predominantly on the larval stage to a recent interest in adult movement. Antitropical distributions provide a unique system to assess vagility and dispersal. In this study, we have focused on an antitropical wrasse genus, Semicossyphus, which includes the California sheephead, S. pulcher, and Darwin's sheephead, S. darwini. Using a phylogenetic approach based on mitochondrial and nuclear markers, and a population genetic approach based on mitochondrial control region sequences and 10 microsatellite loci, we compared the phylogenetic relationships of these two species, as well as the population genetic characteristics within S. pulcher. While S. pulcher and S. darwini are found in the temperate eastern Pacific regions of the northern and southern hemispheres, respectively, their genetic divergence was very small (estimated to have occurred between 200 and 600 kya). Within S. pulcher, genetic structuring was generally weak, especially along mainland California, but showed weak differentiation between Sea of Cortez and California, and between mainland California and Channel Islands. We highlight the congruence of weak genetic differentiation both within and between species and discuss possible causes for maintenance of high gene flow. In particular, we argue that deep and cooler water refugia are used as stepping stones to connect distant populations, resulting in low levels of genetic differentiation.