Ecological speciation in marine v. freshwater fishes

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.

ddRAD‐Seq reveals evolutionary insights into population differentiation and the cryptic phylogeography of Hyporhamphus intermedius in Mainland China

Species differentiation and local adaptation in heterogeneous environments have attracted much attention, although little is known about the mechanisms involved. Hyporhamphus intermedius is an anadromous, brackish‐water halfbeak that is widely distributed in coastal areas and hyperdiverse freshwater systems in China, making it an interesting model for research on phylogeography and local adaptation. Here, 156 individuals were sampled at eight sites from heterogeneous aquatic habitats to examine environmental and genetic contributions to phenotypic divergence. Using double‐digest restriction‐site‐associated DNA sequencing (ddRAD‐Seq) in the specimens from the different watersheds, 5498 single nucleotide polymorphisms (SNPs) were found among populations, with obvious population differentiation. We find that present‐day Mainland China populations are structured into distinct genetic clusters stretching from southern and northern ancestries, mirroring geography. Following a transplant event in Plateau Lakes, there were virtually no variations of genetic diversity occurred in two populations, despite the fact two main splits were unveiled in the demographic history. Additionally, dorsal, and anal fin traits varied widely between the southern group and the others, which highlighted previously unrecognized lineages. We then explore genotype–phenotype‐environment associations and predict candidate loci. Subgroup ranges appeared to correspond to geographic regions with heterogeneous hydrological factors, indicating that these features are likely important drivers of diversification. Accordingly, we conclude that genetic and phenotypic polymorphism and a moderate amount of genetic differentiation occurred, which might be ascribed to population subdivision, and the impact of abiotic factors.

Diversification Rate is Associated with Rate of Molecular Evolution in Ray-Finned Fish (Actinopterygii)

Understanding the factors that drive diversification of taxa across the tree of life is a key focus of macroevolutionary research. While the effects of life history, ecology, climate and geography on diversity have been studied for many taxa, the relationship between molecular evolution and diversification has received less attention. However, correlations between rates of molecular evolution and diversification rate have been detected in a range of taxa, including reptiles, plants and birds. A correlation between rates of molecular evolution and diversification rate is a prediction of several evolutionary theories, including the evolutionary speed hypothesis which links variation in mutation rates to differences in speciation rates. If it is widespread, such correlations could also have significant practical impacts, if they are not adequately accounted for in phylogenetic inference of evolutionary rates and timescales. Ray-finned fish (Actinopterygii) offer a prime target to test for this relationship due to their extreme variation in clade size suggesting a wide range of diversification rates. We employ both a sister-pairs approach and a whole-tree approach to test for correlations between substitution rate and net diversification. We also collect life history and ecological trait data and account for potential confounding factors including body size, latitude, max depth and reef association. We find evidence to support a relationship between diversification and synonymous rates of nuclear evolution across two published backbone phylogenies, as well as weak evidence for a relationship between mitochondrial nonsynonymous rates and diversification at the genus level.

Global diversification and evolutionary history of onchidiid slugs (Gastropoda, Pulmonata)

Many marine species are specialized to specific parts of a habitat. In a mangrove forest, for instance, species may be restricted to the mud surface, the roots and trunks of mangrove trees, or rotting logs, which can be regarded as distinct microhabitats. Shifts to new microhabitats may be an important driver of sympatric speciation. However, the evolutionary history of these shifts is still poorly understood in most groups of marine organisms, because it requires a well-supported phylogeny with relatively complete taxon sampling. Onchidiid slugs are an ideal case study for the evolutionary history of habitat and microhabitat shifts because onchidiid species are specialized to different tidal zones and microhabitats in mangrove forests and rocky shores, and the taxonomy of the family in the Indo-West Pacific has been recently revised in a series of monographs. Here, DNA sequences for onchidiid species from the North and East Pacific, the Caribbean, and the Atlantic are used to reconstruct phylogenetic relationships among Onchidella species, and are combined with new data for Indo-West Pacific species to reconstruct a global phylogeny of the family. The phylogenetic relationships of onchidiid slugs are reconstructed based on three mitochondrial markers (COI, 12S, 16S) and three nuclear markers (28S, ITS2, H3) and nearly complete taxon sampling (all 13 genera and 62 of the 67 species). The highly-supported phylogeny presented here suggests that ancestral onchidiids most likely lived in the rocky intertidal, and that a lineage restricted to the tropical Indo-West Pacific colonized new habitats, including mudflats, mangrove forests, and high-elevation rainforests. Many onchidiid species in the Indo-West Pacific diverged during the Miocene, around the same time that a high diversity of mangrove plants appears in the fossil record, while divergence among Onchidella species occurred earlier, likely beginning in the Eocene. It is demonstrated that ecological specialization to microhabitats underlies the divergence between onchidiid genera, as well as the diversification through sympatric speciation in the genera Wallaconchis and Platevindex. The geographic distributions of onchidiid species also indicate that allopatric speciation played a key role in the diversification of several genera, especially Onchidella and Peronia. The evolutionary history of several morphological traits (penial gland, rectal gland, dorsal eyes, intestinal loops) is examined in relation to habitat and microhabitat evolutionary transitions and that the rectal gland of onchidiids is an adaptation to high intertidal and terrestrial habitats.

Genomic Evidence for Speciation with Gene Flow in Broadcast Spawning Marine Invertebrates

How early stages of speciation in free-spawning marine invertebrates proceed is poorly understood. The Western Pacific abalones, Haliotis discus, H. madaka, and H. gigantea, occur in sympatry with shared breeding season and are capable of producing viable F1 hybrids in spite of being ecologically differentiated. Population genomic analyses revealed that although the three species are genetically distinct, there is evidence for historical and ongoing gene flow among these species. Evidence from demographic modeling suggests that reproductive isolation among the three species started to build in allopatry and has proceeded with gene flow, possibly driven by ecological selection. We identified 27 differentiation islands between the closely related H. discus and H. madaka characterized by high FST and dA, but not high dXY values, as well as high genetic diversity in one H. madaka population. These genomic signatures suggest differentiation driven by recent ecological divergent selection in presence of gene flow outside of the genomic islands of differentiation. The differentiation islands showed low polymorphism in H. gigantea, and both high FST, dXY, and dA values between H. discus and H. gigantea, as well as between H. madaka and H. gigantea. Collectively, the Western Pacific abalones appear to occupy the early stages speciation continuum, and the differentiation islands associated with ecological divergence among the abalones do not appear to have acted as barrier loci to gene flow in the younger divergences but appear to do so in older divergences.

Cryptic species and grey zone speciation of the Barbodes binotatus complex (Teleostei, Cyprinidae) in Sundaland

Morphology-based taxonomy of freshwater fish is effective when there are representative specimens covering large regions. However, in Sundaland, where the presence of cryptic species is high, the technique has its limitations. This is compounded by uncritical descriptions of holotypes in old literature. We demonstrate the problem using Barbodes binotatus first described from an ink drawing. Several species in the Barbodes genus of Sundaland exhibit morphological similarity to B. binotatus. We applied new DNA sequences of 16S, cytochrome c oxidase subunit I (COI), cytochrome b (Cytb) and recombination-activating gene 1 (RAG1), and pigmentation markers to clarify species complex boundaries in the Malay Peninsula, namely B. aff. binotatus "Malay Peninsula", Barbodes cf. banksi and Barbodes rhombeus. Results suggest B. binotatus-like specimens in the Malay Peninsula are B. rhombeus based on a threshold of 3% COI genetic divergence. B. aff. binotatus recorded in Sumatra, Borneo and the Philippines are likely valid but undescribed species. However, if the 2% COI threshold is applied, some populations in the northern Malay Peninsula would qualify as new and undescribed species. The implications of the 2% threshold and the likelihood of "grey zone" incipient populations are discussed. We further found a rapid visual method, not reported previously, to delineate B. aff. binotatus and B. cf. banksi, but it requires further validation. Additionally, we offer fresh perspectives by discussing the roles of biological species concept, morphological species concept, genetic species concept and mate recognition concept in the B. binotatus complex. Our findings reinforce the standpoint that species delineation is not entirely a binary process, but there is a spectrum to consider, especially in biogeography intersection regions.

Rapid adaptive radiation in a hillstream cyprinid fish in the East African White Nile River basin

Adaptive radiation of freshwater fishes was long thought to be possible only in lacustrine environments. Recently, several studies have shown that riverine and stream environments also provide the ecological opportunity for adaptive radiation. In this study, we report on a riverine adaptive radiation of six ecomorphs of cyprinid hillstream fishes of the genus Garra in a river located in the Ethiopian Highlands in East Africa. Garra are predominantly highly specialized algae-scrapers with a wide distribution ranging from Southeast Asia to West Africa. However, adaptive phenotypic diversification in mouth type, sucking disc morphology, gut length and body shape have probably been found among these ecomorphs in a single Ethiopian river. Moreover, we found two novel phenotypes of Garra ("thick-lipped" and "predatory") that had not been discovered before in this species-rich genus (>160 species). Mitochondrial and genome-wide data suggest monophyletic, intrabasin evolution of Garra phenotypic diversity with signatures of gene flow from other local populations. Although sympatric ecomorphs are genetically distinct and can be considered to being young species as suggested by genome-wide single nucleotide polymorphism data, mitochondrial DNA was unable to identify any genetic structure suggesting recent and rapid speciation events. Some data suggest a hybrid origin of the novel "thick-lipped" ecomorph. Here we highlight how, driven by ecological opportunity, an ancestral trophically highly specialized lineage is likely to have rapidly radiated in a riverine environment promoted by the evolution of novel feeding strategies.

Speciation in marine environments: Diving under the surface

Marine environments are inhabited by a broad representation of the tree of life, yet our understanding of speciation in marine ecosystems is extremely limited compared with terrestrial and freshwater environments. Developing a more comprehensive picture of speciation in marine environments requires that we 'dive under the surface' by studying a wider range of taxa and ecosystems is necessary for a more comprehensive picture of speciation. Although studying marine evolutionary processes is often challenging, recent technological advances in different fields, from maritime engineering to genomics, are making it increasingly possible to study speciation of marine life forms across diverse ecosystems and taxa. Motivated by recent research in the field, including the 14 contributions in this issue, we highlight and discuss six axes of research that we think will deepen our understanding of speciation in the marine realm: (a) study a broader range of marine environments and organisms; (b) identify the reproductive barriers driving speciation between marine taxa; (c) understand the role of different genomic architectures underlying reproductive isolation; (d) infer the evolutionary history of divergence using model-based approaches; (e) study patterns of hybridization and introgression between marine taxa; and (f) implement highly interdisciplinary, collaborative research programmes. In outlining these goals, we hope to inspire researchers to continue filling this critical knowledge gap surrounding the origins of marine biodiversity.

Speciation across the Tree of Life

Much of what we know about speciation comes from detailed studies of well-known model systems. Although there have been several important syntheses on speciation, few (if any) have explicitly compared speciation among major groups across the Tree of Life. Here, we synthesize and compare what is known about key aspects of speciation across taxa, including bacteria, protists, fungi, plants, and major animal groups. We focus on three main questions. Is allopatric speciation predominant across groups? How common is ecological divergence of sister species (a requirement for ecological speciation), and on what niche axes do species diverge in each group? What are the reproductive isolating barriers in each group? Our review suggests the following patterns. (i) Based on our survey and projected species numbers, the most frequent speciation process across the Tree of Life may be co-speciation between endosymbiotic bacteria and their insect hosts. (ii) Allopatric speciation appears to be present in all major groups, and may be the most common mode in both animals and plants, based on non-overlapping ranges of sister species. (iii) Full sympatry of sister species is also widespread, and may be more common in fungi than allopatry. (iv) Full sympatry of sister species is more common in some marine animals than in terrestrial and freshwater ones. (v) Ecological divergence of sister species is widespread in all groups, including ~70% of surveyed species pairs of plants and insects. (vi) Major axes of ecological divergence involve species interactions (e.g. host-switching) and habitat divergence. (vii) Prezygotic isolation appears to be generally more widespread and important than postzygotic isolation. (viii) Rates of diversification (and presumably speciation) are strikingly different across groups, with the fastest rates in plants, and successively slower rates in animals, fungi, and protists, with the slowest rates in prokaryotes. Overall, our study represents an initial step towards understanding general patterns in speciation across all organisms.

Genomic landscape of geographically structured colour polymorphism in a temperate marine fish

The study of phenotypic variation patterns among populations is fundamental to elucidate the drivers of evolutionary processes. Empirical evidence that supports ongoing genetic divergence associated with phenotypic variation remains very limited for marine species where larval dispersal is a common homogenizing force. We present a genome‐wide analysis of a marine fish, Labrus bergylta, comprising 144 samples distributed from Norway to Spain, a large geographical area that harbours a gradient of phenotypic differentiation. We analysed 39,602 biallelic single nucleotide polymorphisms and found a clear latitudinal gradient of genomic differentiation strongly correlated with the variation in phenotypic morph frequencies observed across the North Atlantic. We also detected a strong association between the latitude and the number of loci that appear to be under divergent selection, which increased with differences in coloration but not with overall genetic differentiation. Our results demonstrate that strong reproductive isolation is occurring between sympatric colour morphs of L. bergylta found at the southern areas and provide important new insights into the genomic changes shaping early stages of differentiation that might precede speciation with gene flow.

Patterns of genomic divergence and introgression between Japanese stickleback species with overlapping breeding habitats

With only a few absolute geographic barriers in marine environments, the factors maintaining reproductive isolation among marine organisms remain elusive. However, spatial structuring in breeding habitat can contribute to reproductive isolation. This is particularly important for marine organisms that migrate to use fresh- or brackish water environments to breed. The Japanese Gasterosteus stickleback species, the Pacific Ocean three-spined stickleback (G. aculeatus) and the Japan Sea stickleback (G. nipponicus) overwinter in the sea, but migrate to rivers for spawning. Although they co-occur at several locations across the Japanese islands, they are reproductively isolated. Our previous studies in Bekanbeushi River showed that the Japan Sea stickleback spawns in the estuary, while the Pacific Ocean stickleback mainly spawns further upstream in freshwater. Overall genomic divergence was very high with many interspersed regions of introgression. Here, we investigated genomic divergence and introgression between the sympatric species in the much shorter Tokotan River, where they share spawning sites. The levels of genome-wide divergence were reduced and introgression was increased, suggesting that habitat isolation substantially contributes to a reduction in gene flow. We also found that genomic regions of introgression were largely shared between the two systems. Furthermore, some regions of introgression were located near loci with a heterozygote advantage for juvenile survival. Taken together, introgression may be partially driven by adaptation in this system. Although, the two species remain clearly genetically differentiated. Regions with low recombination rates showed especially low introgression. Speciation reversal is therefore likely prevented by barriers other than habitat isolation.

In the rivers: Multiple adaptive radiations of cyprinid fishes (Labeobarbus) in Ethiopian Highlands

Multiple repeated patterns of adaptive radiation were revealed in cyprinid fish inhabiting the compact geographic region of the Ethiopian Highlands. We found four independently evolved radiations in the evolutionary hexaploid (2n = 150) Labeobarbus lineage based on matrilineal relationships of >800 individuals. Each radiation displayed similar patterns of mouth phenotype diversification, and included ecomorphs/species of the generalized, lipped, scraping (one or two), and large-mouthed (one to three) types. All radiations were detected in geographically isolated rivers, and originated from different ancestral populations. This is the first documented case in which numerous parallel radiations of fishes occurred-via different ways-in a riverine environment. Some radiations are very recent and monophyletic, while others are older and include ecomorphs that originated in separate mini flocks and later combined into one. The diversification bursts among Ethiopian Labeobarbus were detected in the mid-upper reaches of rivers (1050-1550 m above sea level), which likely offer ecological opportunities that include diverse habitats yet poor fauna (i.e. lower competition and relaxed selection). This promising example of parallel evolution of adaptive radiation warrants further investigation.

The evolution of ecological specialization across the range of a broadly distributed marine species

Ecological specialization is an important engine of evolutionary change and adaptive radiation, but empirical evidence of local adaptation in marine environments is rare, a pattern that has been attributed to the high dispersal ability of marine taxa and limited geographic barriers to gene flow. The broad-nosed pipefish, Syngnathus typhle, is one of the most broadly distributed syngnathid species and shows pronounced variation in cranial morphology across its range, a factor that may contribute to its success in colonizing new environments. We quantified variation in cranial morphology across the species range using geometric morphometrics, and tested for evidence of trophic specialization by comparing individual-level dietary composition with the community of prey available at each site. Although the diets of juvenile pipefish from each site were qualitatively similar, ontogenetic shifts in dietary composition resulted in adult populations with distinctive diets consistent with their divergent cranial morphology. Morphological differences found in nature are maintained under common garden conditions, indicating that trophic specialization in S. typhle is a heritable trait subject to selection. Our data highlight the potential for ecological specialization in response to spatially variable selection pressures in broadly distributed marine species.

Genomic signatures of host-associated divergence and adaptation in a coral-eating snail, Coralliophila violacea (Kiener, 1836)

The fluid nature of the ocean, combined with planktonic dispersal of marine larvae, lowers physical barriers to gene flow. However, divergence can still occur despite gene flow if strong selection acts on populations occupying different ecological niches. Here, we examined the population genomics of an ectoparasitic snail, Coralliophila violacea (Kiener 1836), that specializes on Porites corals in the Indo-Pacific. Previous genetic analyses revealed two sympatric lineages associated with different coral hosts. In this study, we examined the mechanisms promoting and maintaining the snails' adaptation to their coral hosts. Genome-wide single nucleotide polymorphism (SNP) data from type II restriction site-associated DNA (2b-RAD) sequencing revealed two differentiated clusters of C. violacea that were largely concordant with coral host, consistent with previous genetic results. However, the presence of some admixed genotypes indicates gene flow from one lineage to the other. Combined, these results suggest that differentiation between host-associated lineages of C. violacea is occurring in the face of ongoing gene flow, requiring strong selection. Indeed, 2.7% of all SNP loci were outlier loci (73/2,718), indicative of divergence with gene flow, driven by adaptation of each C. violacea lineage to their specific coral hosts.

Comparative transcriptomics of sympatric species of coral reef fishes (genus: Haemulon)

Background

Coral reefs are major hotspots of diversity for marine fishes, yet there is still ongoing debate on the mechanisms that promote divergence in these rich ecosystems. Our understanding of how diversity originates in this environment could be enhanced by investigating the evolutionary dynamics of closely related fishes with overlapping ranges. Here, we focus on grunts of the genus Haemulon, a group of coral reef fishes with 15 species in the Western Atlantic, 11 of which are syntopic.

Methods

Wild fish samples from three sympatric species of the Caribbean: Haemulon flavolineatum, H. carbonarium and H. macrostomum, were collected while SCUBA diving. RNA was extracted from livers, and the transcriptomes were assembled and annotated to investigate positive selection (Pairwise d _N/d _S) and patterns of gene expression between the three species.

Results

Pairwise d _N/d _S analyses showed evidence of positive selection for genes associated with immune response, cranial morphology and formation of the anterior-posterior axis. Analyses of gene expression revealed that despite their sympatric distribution, H. macrostomum showed upregulation of oxidation-reduction machinery, while there was evidence for activation of immune response in H. carbonarium.

Discussion

Overall, our analyses suggest closely related grunts show important differences in genes associated with body shape and feeding morphology, a result in-line with previous morphological studies in the group. Further, despite their overlapping distribution they interact with their environment in distinct fashions. This is the largest compendium of genomic information for grunts thus far, representing a valuable resource for future studies in this unique group of coral reef fishes.

An exhaustive phylogeny of the combtooth blenny genus Salaria (Pisces, Blenniidae) shows introgressive hybridization and lack of reciprocal mtDNA monophyly between the marine species Salaria basilisca and Salaria pavo

A comprehensive phylogeny of the genus Salaria based on mitochondrial and nuclear markers grouped the extant species of the genus in well-characterised marine and freshwater clades, thus rejecting the hypothesis of a polytypic origin of the freshwater Salaria populations and supporting the occurrence of a single invasion event of the inland waters by the genus. Based on both mitochondrial and nuclear DNA datasets, the Salaria species of the freshwater clade proved to be vicariant taxa originating from a common ancestor which could possibly spread throughout the circum-Mediterranean inland waters during the late Miocene Messinian salinity crisis, then experiencing a process of allopatric differentiation after the re-flooding of the Mediterranean basin. Within the marine clade, although the nuDNA datasets showed the existence of well-supported subclades in accordance to the morphological identification of the studied specimens, one of the two subclades obtained in the phylogenetic tree based on the mtDNA dataset included both S. basilisca and S. pavo specimens, thus failing to find the two species as reciprocally monophyletic. Such a mito-nuclear discordance is here ascribed to multiple mtDNA unidirectional introgression events from S. basilisca to S. pavo, and the molecular diversity pattern of the marine Salaria species is here ascribed to a Pleistocene speciation event nowadays partly concealed by the occurrence of introgressive hybridization phenomena between the two taxa. Our results urge for prudence when implementing DNA barcoding approaches since, in the presence of mito-nuclear discordance phenomena, single-marker mtDNA-only analyses might lead to significant misidentifications.

Insights into the Genomics of Clownfish Adaptive Radiation: Genetic Basis of the Mutualism with Sea Anemones

Clownfishes are an iconic group of coral reef fishes, especially known for their mutualism with sea anemones. This mutualism is particularly interesting as it likely acted as the key innovation that triggered clownfish adaptive radiation. Indeed, after the acquisition of the mutualism, clownfishes diversified into multiple ecological niches linked with host and habitat use. However, despite the importance of this mutualism, the genetic mechanisms allowing clownfishes to interact with sea anemones are still unclear. Here, we used a comparative genomics and molecular evolutionary analyses to investigate the genetic basis of clownfish mutualism with sea anemones. We assembled and annotated the genome of nine clownfish species and one closely related outgroup. Orthologous genes inferred between these species and additional publicly available teleost genomes resulted in almost 16,000 genes that were tested for positively selected substitutions potentially involved in the adaptation of clownfishes to live in sea anemones. We identified 17 genes with a signal of positive selection at the origin of clownfish radiation. Two of them (Versican core protein and Protein O-GlcNAse) show particularly interesting functions associated with N-acetylated sugars, which are known to be involved in sea anemone discharge of toxins. This study provides the first insights into the genetic mechanisms of clownfish mutualism with sea anemones. Indeed, we identified the first candidate genes likely to be associated with clownfish protection form sea anemones, and thus the evolution of their mutualism. Additionally, the genomic resources acquired represent a valuable resource for further investigation of the genomic basis of clownfish adaptive radiation.

Origins of the Greenland shark (Somniosus microcephalus): Impacts of ice-olation and introgression

Herein, we use genetic data from 277 sleeper sharks to perform coalescent‐based modeling to test the hypothesis of early Quaternary emergence of the Greenland shark (Somniosus microcephalus) from ancestral sleeper sharks in the Canadian Arctic‐Subarctic region. Our results show that morphologically cryptic somniosids S. microcephalus and Somniosus pacificus can be genetically distinguished using combined mitochondrial and nuclear DNA markers. Our data confirm the presence of genetically admixed individuals in the Canadian Arctic and sub‐Arctic, and temperate Eastern Atlantic regions, suggesting introgressive hybridization upon secondary contact following the initial species divergence. Conservative substitution rates fitted to an Isolation with Migration (IM) model indicate a likely species divergence time of 2.34 Ma, using the mitochondrial sequence DNA, which in conjunction with the geographic distribution of admixtures and Pacific signatures likely indicates speciation associated with processes other than the closing of the Isthmus of Panama. This time span coincides with further planetary cooling in the early Quaternary period followed by the onset of oscillating glacial‐interglacial cycles. We propose that the initial S. microcephalus–S. pacificus split, and subsequent hybridization events, were likely associated with the onset of Pleistocene glacial oscillations, whereby fluctuating sea levels constrained connectivity among Arctic oceanic basins, Arctic marginal seas, and the North Atlantic Ocean. Our data demonstrates support for the evolutionary consequences of oscillatory vicariance via transient oceanic isolation with subsequent secondary contact associated with fluctuating sea levels throughout the Quaternary period—which may serve as a model for the origins of Arctic marine fauna on a broad taxonomic scale.

Evolution of Hemoglobin Genes in Codfishes Influenced by Ocean Depth

Understanding the genetic basis of adaptation is one of the main enigmas of evolutionary biology. Among vertebrates, hemoglobin has been well documented as a key trait for adaptation to different environments. Here, we investigate the role of hemoglobins in adaptation to ocean depth in the diverse teleost order Gadiformes, with species distributed at a wide range of depths varying in temperature, hydrostatic pressure and oxygen levels. Using genomic data we characterized the full hemoglobin (Hb) gene repertoire for subset of species within this lineage. We discovered a correlation between expanded numbers of Hb genes and ocean depth, with the highest numbers in species occupying shallower, epipelagic regions. Moreover, we demonstrate that the Hb genes have functionally diverged through diversifying selection. Our results suggest that the more variable environment in shallower water has led to selection for a larger Hb gene repertoire and that Hbs have a key role in adaptive processes in marine environments.

Extraordinarily rapid speciation in a marine fish

Divergent selection may initiate ecological speciation extremely rapidly. How often and at what pace ecological speciation proceeds to yield strong reproductive isolation is more uncertain. Here, we document a case of extraordinarily rapid speciation associated with ecological selection in the postglacial Baltic Sea. European flounders (Platichthys flesus) in the Baltic exhibit two contrasting reproductive behaviors: pelagic and demersal spawning. Demersal spawning enables flounders to thrive in the low salinity of the Northern Baltic, where eggs cannot achieve neutral buoyancy. We show that demersal and pelagic flounders are a species pair arising from a recent event of speciation. Despite having a parapatric distribution with extensive overlap, the two species are reciprocally monophyletic and show strongly bimodal genotypic clustering and no evidence of contemporary migration, suggesting strong reproductive isolation. Divergence across the genome is weak but shows strong signatures of selection, a pattern suggestive of a recent ecological speciation event. We propose that spawning behavior in Baltic flounders is the trait under ecologically based selection causing reproductive isolation, directly implicating a process of ecological speciation. We evaluated different possible evolutionary scenarios under the approximate Bayesian computation framework and estimate that the speciation process started in allopatry ∼2,400 generations ago, following the colonization of the Baltic by the demersal lineage. This is faster than most known cases of ecological speciation and represents the most rapid event of speciation ever reported for any marine vertebrate.

Assessing pre- and post-zygotic barriers between North Atlantic eels (Anguilla anguilla and A. rostrata)

Elucidating barriers to gene flow is important for understanding the dynamics of speciation. Here we investigate pre- and post-zygotic mechanisms acting between the two hybridizing species of Atlantic eels: Anguilla anguilla and A. rostrata. Temporally varying hybridization was examined by analyzing 85 species-diagnostic single-nucleotide polymorphisms (SNPs; F_ST ⩾0.95) in eel larvae sampled in the spawning region in the Sargasso Sea in 2007 (N=92) and 2014 (N=460). We further investigated whether genotypes at these SNPs were nonrandomly distributed in post-F1 hybrids, indicating selection. Finally, we sequenced the mitochondrial ATP6 and nuclear ATP5c1 genes in 19 hybrids, identified using SNP and restriction site associated DNA (RAD) sequencing data, to test a previously proposed hypothesis of cytonuclear incompatibility leading to adenosine triphosphate (ATP) synthase dysfunction and selection against hybrids. No F1 hybrids but only later backcrosses were observed in the Sargasso Sea in 2007 and 2014. This suggests that interbreeding between the two species only occurs in some years, possibly controlled by environmental conditions at the spawning grounds, or that interbreeding has diminished through time as a result of a declining number of spawners. Moreover, potential selection was found at the nuclear and the cytonuclear levels. Nonetheless, one glass eel individual showed a mismatch, involving an American ATP6 haplotype and European ATP5c1 alleles. This contradicted the presence of cytonuclear incompatibility but may be explained by that (1) cytonuclear incompatibility is incomplete, (2) selection acts at a later life stage or (3) other genes are important for protein function. In total, the study demonstrates the utility of genomic data when examining pre- and post-zyotic barriers in natural hybrids.

Decrease of Population Divergence in Eurasian Perch (Perca fluviatilis) in Browning Waters: Role of Fatty Acids and Foraging Efficiency

Due to altered biogeochemical processes related to climate change, highly colored dissolved organic carbon (DOC) from terrestrial sources will lead to a water “brownification” in many freshwater systems of the Northern Hemisphere. This will create deteriorated visual conditions that have been found to affect habitat-specific morphological variations in Eurasian perch (Perca fluviatilis) in a previous study. So far, potential drivers and ultimate causes of these findings have not been identified. We conducted a field study to investigate the connection between morphological divergence and polyunsaturated fatty acid (PUFA) composition of perch from six lakes across a gradient of DOC concentration. We expected a decrease in the prevalence of PUFAs, which are important for perch growth and divergence with increasing DOC concentrations, due to the restructuring effects of DOC on aquatic food webs. In general, rate of morphological divergence in perch decreased with increasing DOC concentrations. Proportions of specific PUFAs (22:6n-3, 18:3n-3, 20:5n-3, and 20:4n-6) identified to primarily contribute to overall differences between perch caught in clear and brown-water lakes tended to be connected to overall decline of morphological divergence. However, no overall significant relationship was found, indicating no severe limitation of essential fatty acids for perch inhabiting brown water lakes. We further broaden our approach by conducting a laboratory experiment on foraging efficiency of perch. Therefore, we induced pelagic and littoral phenotypes by differences in habitat-structure and feeding mode and recorded attack rate in a feeding experiment. Generally, fish were less efficient in foraging on littoral prey (Ephemeroptera) when visual conditions were degraded by brown water color. We concluded that browning water may have a strong effect on the forager’s ability to find particular food resources, resulting in the reduced development of evolutionary traits, such as habitat- specific morphological divergence.

Genomic atolls of differentiation in coral reef fishes (Hypoplectrus spp., Serranidae)

Because the vast majority of species are well diverged, relatively little is known about the genomic architecture of speciation during the early stages of divergence. Species within recent evolutionary radiations are often minimally diverged from a genomic perspective, and therefore provide rare opportunities to address this question. Here, we leverage the hamlet radiation (Hypoplectrus spp., brightly coloured reef fishes from the tropical western Atlantic) to characterize genomic divergence during the early stages of speciation. Transect surveys and spawning observations in Belize, Honduras and Panama confirm that sympatric barred (H. puella), black (H. nigricans) and butter (H. unicolor) hamlets are phenotypically distinct and reproductively isolated, although hybrid spawnings and individuals with intermediate phenotypes are seen on rare occasions. A survey of approximately 100 000 restriction site-associated SNPs in 126 samples from the three species across the three replicate populations reveals extremely slight genomewide divergence among species (FST  = 0.0038), indicating that ecomorphological differences and functional reproductive isolation are maintained in sympatry in a backdrop of extraordinary genomic similarity. Nonetheless, a very small proportion of SNPs (0.05% on average) are identified as FST outliers among sympatric species. Remarkably, a single SNP is identified as an outlier in repeated populations for the same species pair. A minicontig assembled de novo around this SNP falls into the genomic region containing the HoxCa10 and HoxCa11 genes in 10 teleost species, suggesting an important role for Hox gene evolution in this radiation. This finding, if confirmed, would provide a better understanding of the links between micro- and macroevolutionary processes.

Delphinid systematics and biogeography with a focus on the current genus Lagenorhynchus: multiple pathways for antitropical and trans-oceanic radiation

The six species currently classified within the genus Lagenorhynchus exhibit a pattern of antitropical distribution common among marine taxa. In spite of their morphological similarities they are now considered an artificial grouping, and include both recent and the oldest representatives of the Delphinidae radiation. They are, therefore, a good model for studying questions about the evolutionary processes that have driven dolphin speciation, dispersion and distribution. Here we used two different approaches. First we constructed a multigenic phylogeny with a minimum amount of missing data (based on 9 genes, 11,030bp, using the 6 species of the genus and their closest relatives) to infer their relationships. Second, we built a supermatrix phylogeny (based on 33 species and 27 genes) to test the effect of taxon sampling on the phylogeny of the genus, to provide inference on biogeographic history, and provide inference on the main events shaping the dispersion and radiation of delphinids. Our analyses suggested an early evolutionary history of marine dolphins in the North Atlantic Ocean and revealed multiple pathways of migration and radiation, probably guided by paleoceanographic changes during the Miocene and Pliocene. L. acutus and L. albirostris likely shared a common ancestor that arose in the North Atlantic around the Middle Miocene, predating the radiation of subfamilies Delphininae, Globicephalinae and Lissodelphininae.

Population structure of guppies in north-eastern Venezuela, the area of putative incipient speciation

BACKGROUND

Geographic barriers to gene flow and divergence among populations in sexual traits are two important causes of genetic isolation which may lead to speciation. Genetic isolation may be facilitated if these two mechanisms act synergistically. The guppy from the Cumaná region (within the Cariaco drainage) of eastern Venezuela has been previously described as a case of incipient speciation driven by sexual selection, significantly differentiated in sexual colouration and body shape from the common guppy, Poecilia reticulata. The latter occurs widely in northern Venezuela, including the south-eastern side of Cordillera de la Costa, where it inhabits streams belonging to the San Juan drainage. Here, we present molecular and morphological analyses of differentiation among guppy populations in the Cariaco and San Juan drainages. Our analyses are based on a 953 bp long mtDNA fragment, a set of 15 microsatellites (519 fish from 20 populations), and four phenotypic traits.

RESULTS

Both microsatellite and mtDNA data showed that guppies inhabiting the two drainages are characterised by a significant genetic differentiation, but a higher proportion of the genetic variance was distributed among populations within regions. Most guppies in the Cariaco drainage had mtDNA from a distinct lineage, but we also found evidence for widespread introgression of mtDNA from the San Juan drainage into the Cariaco drainage. Phenotypically, populations in the two regions differed significantly only in the number of black crescents. Phenotypic clustering did not support existence of two distinct groupings, but indicated a degree of distinctiveness of Central Cumaná (CC) population. However, CC population showed little differentiation at the neutral markers from the proximate populations within the Cariaco drainage.

CONCLUSIONS

Our findings are consistent with only partial genetic isolation between the two geographic regions and indicate that the geographic barrier of Cordillera de la Costa has not played an important role in strengthening the incomplete pre-zygotic reproductive barrier between Cumaná and common guppy. Significant phenotypic differentiation between genetically similar (in terms of neutral variation) populations suggests that mate choice can maintain divergence at sexually selected traits despite gene flow. However, neither genetic nor phenotypic clustering supported delineation of two species within the region.

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.

Speciation in fishes

The field of speciation has seen much renewed interest in the past few years, with theoretical and empirical advances that have moved it from a descriptive field to a predictive and testable one. The goal of this review is to provide a general background on research on speciation as it pertains to fishes. Three major components to the question are first discussed: the spatial, ecological and sexual factors that influence speciation mechanisms. We then move to the latest developments in the field of speciation genomics. Affordable and rapidly available, massively parallel sequencing data allow speciation studies to converge into a single comprehensive line of investigation, where the focus has shifted to the search for speciation genes and genomic islands of speciation. We argue that fish present a very diverse array of scenarios, making them an ideal model to study speciation processes.

Demographic variation between colour patterns in a temperate protogynous hermaphrodite, the ballan wrasse Labrus bergylta

Fish populations are often treated as homogeneous units in typical fishery management, thereby tacitly ignoring potential intraspecific variation which can lead to imprecise management rules. However, intraspecific variation in life-history traits is widespread and related to a variety of factors. We investigated the comparative age-based demography of the two main colour patterns of Labrus bergylta (plain and spotted, which coexist in sympatry), a commercially valuable resource in the NE Atlantic. Individuals were aged based on otolith readings after validating the annual periodicity of annuli deposition. The relationships between the otolith weight and fish age and between otolith length and fish length were strong but differed between colour patterns. The fit of the growth models to the age and length data resulted in divergent growth curves between colour morphotypes and between sexes. Males and spotted individuals attained larger mean asymptotic sizes (L_inf) than females and plain individuals, respectively, but converged to them more slowly (smaller k). Estimates of mortality based on catch curves from two independent datasets provided a global total mortality (Z) of 0.35 yr^–1, although Z was larger in plain and female individuals. Overall, the results of this research have direct implications for management of L. bergylta and, as a precautionary measure, we recommend considering both colour patterns as two different management units.

Triton's trident: cryptic Neogene divergences in a marine clam (Lasaea australis) correspond to Australia's three temperate biogeographic provinces

The southern coast of Australia is composed of three distinct biogeographic provinces distinguished primarily by intertidal community composition. Several ecological mechanisms have been proposed to explain their formation and persistence, but no consensus has been reached. The marine clam Lasaea australis is arguably the most common bivalve on southern Australian rocky shores and occurs in all three provinces. Here, we tested if this species exhibits cryptic genetic structuring corresponding to the provinces and if so, what mechanisms potentially drove its divergence. Variation in two mitochondrial genes (16S and COIII) and one nuclear gene (ITS2) was assayed to test for genetic structuring and to reconstruct the clam's phylogenetic history. Our results showed that L. australis is comprised of three cryptic mitochondrial clades, each corresponding almost perfectly to one of the three biogeographic provinces. Divergence time estimates place their cladogenesis in the Neogene. The trident-like topology and Neogene time frame of L. australis cladogenesis are incongruent with Quaternary vicariance predictions: a two-clade topology produced by Pleistocene Bass Strait land bridge formation. We hypothesize that the interaction of the Middle Miocene Climate Transition with the specific geography of the southern coastline of Australia was the primary cladogenic driver in this clam lineage. Additional in-depth studies of the endemic southern Australian marine biota across all three provinces are needed to establish the generality of this proposed older framework for regional cladogenesis.

Mutualism with sea anemones triggered the adaptive radiation of clownfishes

Background

Adaptive radiation is the process by which a single ancestral species diversifies into many descendants adapted to exploit a wide range of habitats. The appearance of ecological opportunities, or the colonisation or adaptation to novel ecological resources, has been documented to promote adaptive radiation in many classic examples. Mutualistic interactions allow species to access resources untapped by competitors, but evidence shows that the effect of mutualism on species diversification can greatly vary among mutualistic systems. Here, we test whether the development of obligate mutualism with sea anemones allowed the clownfishes to radiate adaptively across the Indian and western Pacific oceans reef habitats.

Results

We show that clownfishes morphological characters are linked with ecological niches associated with the sea anemones. This pattern is consistent with the ecological speciation hypothesis. Furthermore, the clownfishes show an increase in the rate of species diversification as well as rate of morphological evolution compared to their closest relatives without anemone mutualistic associations.

Conclusions

The effect of mutualism on species diversification has only been studied in a limited number of groups. We present a case of adaptive radiation where mutualistic interaction is the likely key innovation, providing new insights into the mechanisms involved in the buildup of biodiversity. Due to a lack of barriers to dispersal, ecological speciation is rare in marine environments. Particular life-history characteristics of clownfishes likely reinforced reproductive isolation between populations, allowing rapid species diversification.

Multilocus phylogeny and recent rapid radiation of the viviparous sea snakes (Elapidae: Hydrophiinae)

The viviparous sea snakes (Hydrophiinae: Hydrophiini) comprise a young but morphologically and ecologically diverse clade distributed throughout the Indo-Pacific. Despite presenting a very promising model for marine diversification studies, many relationships among the 62 species and 16 genera in Hydrophiini remain unresolved. Here, we extend previous taxonomic and genomic sampling for Hydrophiini using three mitochondrial fragments and five nuclear loci for multiple individuals of 39 species in 15 genera. Our results highlight many of the impediments to inferring phylogenies in recent rapid radiations, including low variation at all five nuclear markers, and conflicting relationships supported by mitochondrial and nuclear trees. However, concatenated Bayesian and likelihood analyses, and a multilocus coalescent tree, recovered concordant support for primary clades and several previously unresolved inter-specific groupings. The Aipysurus group is monophyletic, with egg-eating specialists forming separate, early-diverging lineages. All three monotypic semi-aquatic genera (Ephalophis, Parahydrophis and Hydrelaps) are robustly placed as early diverging lineages along the branch leading to the Hydrophis group, with Ephalophis recovered as sister to Parahydrophis. The molecular phylogeny implies extensive evolutionary convergence in feeding adaptations within the Hydrophis group, especially the repeated evolution of small-headed (microcephalic) forms. Microcephalophis (Hydrophis) gracilis is robustly recovered as a relatively distant sister lineage to all other sampled Hydrophis group species, here termed the 'core Hydrophis group'. Within the 'core Hydrophis group', Hydrophis is recovered as broadly paraphyletic, with several other genera nested within it (Pelamis, Enhydrina, Astrotia, Thalassophina, Acalyptophis, Kerilia, Lapemis, Disteira). Instead of erecting multiple new genera, we recommend dismantling the latter (mostly monotypic) genera and recognising a single genus, Hydrophis Latreille 1802, for the core Hydrophis group. Estimated divergence times suggest that all Hydrophiini last shared a common ancestor ∼6million years ago, but that the majority of extant lineages diversified over the last ∼3.5million years. The core Hydrophis group is a young and rapidly speciating clade, with 26 sampled species and 9 genera and dated at only ∼1.5-3million years old.

Ecological partitioning and diversity in tropical planktonic foraminifera

BACKGROUND

Ecological processes are increasingly being viewed as an important mode of diversification in the marine environment, where the high dispersal potential of pelagic organisms, and a lack of absolute barriers to gene flow may limit the occurrence of allopatric speciation through vicariance. Here we focus on the potential role of ecological partitioning in the diversification of a widely distributed group of marine protists, the planktonic foraminifera. Sampling was conducted in the tropical Arabian Sea, during the southwest (summer) monsoon, when pronounced environmental conditions result in a strong disparity in temperature, salinity and productivity between distinct northern and southern water masses.

RESULTS

We uncovered extensive genetic diversity within the Arabian Sea planktonic foraminifera, identifying 13 morphospecies, represented by 20 distinct SSU rRNA genetic types. Several morphospecies/genetic types displayed non-random biogeographical distributions, partitioning between the northern and southern water masses, giving a strong indication of independent ecological adaptations.

CONCLUSIONS

We propose sea-surface primary productivity as the main factor driving the geographical segregation of Arabian Sea planktonic foraminifera, during the SW monsoon, with variations in symbiotic associations possibly playing a role in the specific ecological adaptations observed. Our findings suggest that ecological partitioning could be contributing to the high levels of 'cryptic' genetic diversity observed within the planktonic foraminifera, and support the view that ecological processes may play a key role in the diversification of marine pelagic organisms.

The genetic signature of recent speciation in manta rays (Manta alfredi and M. birostris)

Manta rays have been taxonomically revised as two species, Manta alfredi and M. birostris, on the basis of morphological and meristic data, yet the two species occur in extensive mosaic sympatry. We analysed the genetic signatures of the species boundary using a portion of the nuclear RAG1 (681 base pairs), mitochondrial CO1 (574 bp) and ND5 genes (1188 bp). The assay with CO1 sequences, widely used in DNA barcoding, failed to distinguish the two species. The two species were clearly distinguishable, however, with no shared RAG1 or ND5 haplotypes. The species were reciprocally monophyletic for RAG1, but paraphyletic for ND5 sequences. Qualitative evidence and statistical inferences using the 'Isolation-with-Migration models' indicated that these results were better explained with post-divergence gene flow in the recent past rather than incomplete lineage sorting with zero gene flow since speciation. An estimate of divergence time was less than 0.5 Ma with an upper confidence limit of within 1 Ma. Recent speciation of highly mobile species in the marine environment is of great interest, as it suggests that speciation may have occurred in the absence of long-term physical barriers to gene flow. We propose that the ecologically driven forces such as habitat choice played a significant role in speciation in manta rays.

Population connectivity: dam migration mitigations and contemporary site fidelity in arctic char

BACKGROUND

Animal feeding and spawning migrations may be limited by physical barriers and behavioral interactions. Dam constructions (e.g. hydropower) commonly include gateways for fish migrations to sustain ecological connectivity. Relative genetic impacts of fish passage devices versus natural processes (e.g. hybrid inferiority) are, however, rarely studied. We examined genetic (i.e. microsatellite) population connectivity of highly migrating lake-dwelling Arctic char (Salvelinus alpinus), introduced 20 generations ago, across and within two subalpine lakes separated by a dam with a subterranean tunnel and spill gates after 7 generations. Due to water flow regime, the time window for fish migration is highly restricted.

RESULTS

Char populations, with similar genetic structuring and diversity observed across and within lakes, were admixed across the dam with fishways during feeding. For spawning, however, statistically significant, but very low population differentiation (θ; 0.002 - 0.013) was found in nine out of ten reproductive site comparisons, reflecting interactions between extensive migration (mean first generation (F0) = 10.8%) and initial site fidelity. Simulations indicated that genetic drift among relatively small effective populations (mean N(e) = 62) may have caused the observed contemporary differentiation. Novel Bayesian analyses indicated mean contributions of 71% F0 population hybrids in spawning populations, of which 76% had maternal or paternal native origin.

CONCLUSIONS

Ecological connectivity between lakes separated by a dam has been retained through construction of fishways for feeding migration. Considerable survival and homing to ancestral spawning sites in hybrid progeny was documented. Population differentiation despite preceding admixture is likely caused by contemporary reduced reproductive fitness of population hybrids. The study documents the beginning stages of population divergence among spatial aggregations with recent common ancestry.

A transcriptomic scan for positively selected genes in two closely related marine fishes: Sebastes caurinus and S. rastrelliger

Comparative genomic analyses can provide valuable insight into functional evolutionary divergence among closely related species. Here we employ a comparative evolutionary analysis of expressed sequence tags (ESTs) from two closely related species of marine fishes (genus Sebastes--rockfish). Sebastes is a highly diverse group of marine fishes that inhabit a wide array of marine habitats and the study of this group can provide insights into speciation in the marine environment. ESTs were developed for S. caurinus (23,668 from brain, kidney, and spleen tissues) and S. rastrelliger (11,207 from brain and pituitary tissues). Following assembly we were able to identify, with high confidence, 257 orthologous sequence pairs between the two species through a reciprocal best hit blast search. An analysis of functional divergence between orthologs revealed that 19.46% had Ka/Ks values greater than 0.5 and 8.17% had Ka/Ks values greater than one, identifying a large pool of candidate genes to further study adaptive divergence in the group. Genes with elevated Ka/Ks values belonged to the following functional categories: immune function, metabolism, longevity, and reproductive behavior, indicating that adaptive divergence in these functional groups may be important in the diversification of this group of fishes. This study provides the ground work to better understand the molecular evolution of genes involved in a radiation of marine fishes.

Speciation along a depth gradient in a marine adaptive radiation

Oceans are home to much of the world's biodiversity, but we know little about the processes driving speciation in marine ecosystems with few geographical barriers to gene flow. Ecological speciation resulting from divergent natural selection between ecological niches can occur in the face of gene flow. Sister species in the young and ecologically diverse rockfish genus Sebastes coexist in the northeast Pacific, implying that speciation may not require geographical isolation. Here, I use a novel phylogenetic comparative analysis to show that rockfish speciation is instead associated with divergence in habitat depth and depth-associated morphology, consistent with models of parapatric speciation. Using the same analysis, I find no support for alternative hypotheses that speciation involves divergence in diet or life history, or that speciation involves geographic isolation by latitude. These findings support the hypothesis that rockfishes undergo ecological speciation on an environmental gradient.