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

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

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 (F_ST = 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 F_ST 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.     

Population genetic analyses of Hypoplectrus coral reef fishes provide evidence that local processes are operating during the early stages of marine adaptive radiations

Large-scale, spatially explicit models of adaptive radiation suggest that the spatial genetic structure within a species sampled early in the evolutionary history of an adaptive radiation might be higher than the genetic differentiation between different species formed during the same radiation over all locations. Here we test this hypothesis with a spatial population genetic analysis of Hypoplectrus coral reef fishes (Serranidae), one of the few potential cases of a recent adaptive radiation documented in the marine realm. Microsatellite analyses of Hypoplectrus puella (barred hamlet) and Hypoplectrus nigricans (black hamlet) from Belize, Panama and Barbados validate the population genetic predictions at the regional scale for H. nigricans despite the potential for high levels of gene flow between populations resulting from the 3-week planktonic larval phase of Hypoplectrus . The results are different for H. puella , which is characterized by significantly lower levels of spatial genetic structure than H. nigricans . An extensive field survey of Hypoplectrus population densities complemented by individual-based simulations shows that the higher abundance and more continuous distribution of H. puella could account for the reduced spatial genetic structure within this species. The genetic and demographic data are also consistent with the hypothesis that H. puella might represent the ancestral form of the Hypoplectrus radiation, and that H. nigricans might have evolved repeatedly from H. puella through ecological speciation. Altogether, spatial genetic analysis within and between Hypoplectrus species indicate that local processes can operate at a regional scale within recent marine adaptive radiations.     

Model‐based demographic inference of introgression history in European whitefish species pairs'

Parallel phenotypic differentiation is generally attributed to parallel adaptive divergence as an evolutionary response to similar environmental contrasts. Such parallelism may actually originate from several evolutionary scenarios ranging from repeated parallel divergence caused by divergent selection to a unique divergence event followed by gene flow. Reconstructing the evolutionary history underlying parallel phenotypic differentiation is thus fundamental to understand the relative contribution of demography and selection on genomic divergence during speciation. In this study, we investigate the divergence history of replicate European whitefish (Coregonus lavaretus), limnetic and benthic species pairs from two lakes in Norway and two lakes in Switzerland. Demographic models accounting for semi‐permeability and linked selection were fitted to the unfolded joint allele frequency spectrum built from genome‐wide SNPs and compared to each other in each species pair. We found strong support for a model of asymmetrical post‐glacial secondary contact between glacial lineages in all four lakes. Moreover, our results suggest that heterogeneous genomic differentiation has been shaped by the joint action of linked selection accelerating lineage sorting during allopatry, and heterogeneous migration eroding divergence at different rates along the genome following secondary contact. Our analyses reveal how the interplay between demography, selection and historical contingency has influenced the levels of diversity observed in previous whitefish phylogeographic studies. This study thus provides new insights into the historical demographic and selective processes that shaped the divergence associated with ecological speciation in European whitefish.     

Population genomics of local adaptation versus speciation in coral reef fishes (Hypoplectrus spp,Serranidae)

Are the population genomic patterns underlying local adaptation and the early stages of speciation similar? Addressing this question requires a system in which (i) local adaptation and the early stages of speciation can be clearly identified and distinguished, (ii) the amount of genetic divergence driven by the two processes is similar, and (iii) comparisons can be repeated both taxonomically (for local adaptation) and geographically (for speciation). Here, we report just such a situation in the hamlets (Hypoplectrus spp), brightly colored reef fishes from the wider Caribbean. Close to 100,000 SNPs genotyped in 126 individuals from three sympatric species sampled in three repeated populations provide genome‐wide levels of divergence that are comparable among allopatric populations (F_st estimate = 0.0042) and sympatric species (F_st estimate = 0.0038). Population genetic, clustering, and phylogenetic analyses reveal very similar patterns for local adaptation and speciation, with a large fraction of the genome undifferentiated (F_st estimate ≈ 0), a very small proportion of F_st outlier loci (0.05–0.07%), and remarkably few repeated outliers (1–3). Nevertheless, different loci appear to be involved in the two processes in Hypoplectrus, with only 7% of the most differentiated SNPs and outliers shared between populations and species comparisons. In particular, a tropomyosin (Tpm4) and a previously identified hox (HoxCa) locus emerge as candidate loci (repeated outliers) for local adaptation and speciation, respectively. We conclude that marine populations may be locally adapted notwithstanding shallow levels of genetic divergence, and that from a population genomic perspective, this process does not appear to differ fundamentally from the early stages of speciation.     

Increased CO2 stimulates reproduction in a coral reef fish

Ocean acidification is predicted to negatively impact the reproduction of many marine species, either by reducing fertilization success or diverting energy from reproductive effort. While recent studies have demonstrated how ocean acidification will affect larval and juvenile fishes, little is known about how increasing partial pressure of carbon dioxide (pCO₂) and decreasing pH might affect reproduction in adult fishes. We investigated the effects of near‐future levels of pCO₂ on the reproductive performance of the cinnamon anemonefish, Amphiprion melanopus, from the Great Barrier Reef, Australia. Breeding pairs were held under three CO₂ treatments [Current‐day Control (430 μatm), Moderate (584 μatm) and High (1032 μatm)] for a 9‐month period that included the summer breeding season. Unexpectedly, increased CO₂ dramatically stimulated breeding activity in this species of fish. Over twice as many pairs bred in the Moderate (67% of pairs) and High (55%) compared to the Control (27%) CO₂ treatment. Pairs in the High CO₂ group produced double the number of clutches per pair and 67% more eggs per clutch compared to the Moderate and Control groups. As a result, reproductive output in the High group was 82% higher than that in the Control group and 50% higher than that in the Moderate group. Despite the increase in reproductive activity, there was no difference in adult body condition among the three treatment groups. There was no significant difference in hatchling length between the treatment groups, but larvae from the High CO₂ group had smaller yolks than Controls. This study provides the first evidence of the potential effects of ocean acidification on key reproductive attributes of marine fishes and, contrary to expectations, demonstrates an initially stimulatory (hormetic) effect in response to increased pCO₂. However, any long‐term consequences of increased reproductive effort on individuals or populations remain to be determined.     

Colonisation of toxic environments drives predictable life‐history evolution in livebearing fishes (Poeciliidae)

New World livebearing fishes (family Poeciliidae) have repeatedly colonised toxic, hydrogen sulphide‐rich waters across their natural distribution. Physiological considerations and life‐history theory predict that these adverse conditions should favour the evolution of larger offspring. Here, we examined nine poeciliid species that independently colonised toxic environments, and show that these fishes have indeed repeatedly evolved much larger offspring size at birth in sulphidic waters, thus uncovering a widespread pattern of predictable evolution. However, a second pattern, only indirectly predicted by theory, proved additionally common: a reduction in the number of offspring carried per clutch (i.e. lower fecundity). Our analyses reveal that this secondary pattern represents a mere consequence of a classic life‐history trade‐off combined with strong selection on offspring size alone. With such strong natural selection in extreme environments, extremophile organisms may commonly exhibit multivariate phenotypic shifts even though not all diverging traits necessarily represent adaptations to the extreme conditions.     

Urbanization drives contemporary evolution in stream fish

Human activities reduce biodiversity but may also drive diversification by modifying selection. Urbanization alters stream hydrology by increasing peak water velocities, which should in turn alter selection on the body morphology of aquatic species. Here, we show how urbanization can generate evolutionary divergence in the body morphology of two species of stream fish, western blacknose dace (Rhinichthys obtusus) and creek chub (Semotilus atromaculatus). We predicted that fish should evolve more streamlined body shapes within urbanized streams. We found that in urban streams, dace consistently exhibited more streamlined bodies while chub consistently showed deeper bodies. Comparing modern creek chub populations with historical museum collections spanning 50 years, we found that creek chub (1) rapidly became deeper bodied in streams that experienced increasing urbanization over time, (2) had already achieved deepened bodies 50 years ago in streams that were then already urban (and showed no additional deepening over time), and (3) remained relatively shallow bodied in streams that stayed rural over time. By raising creek chub from five populations under common conditions in the laboratory, we found that morphological differences largely reflected genetically based differences, not velocity–induced phenotypic plasticity. We suggest that urbanization can drive rapid, adaptive evolutionary responses to disturbance, and that these responses may vary unpredictably in different species.     

Effects of cyanobacteria,fish kairomones,and the presence of ostracods on the demography of Simocephalus vetulus (Cladocera)

Simocephalus vetulus is a large (2.0–4.0 mm at maturity) cladoceran often found in the littoral region of lakes and ponds, and capable of moderate growth rates even on poor‐quality cyanobacterial diets. It frequently co‐occurs with fishes and similar sized ostracods such as Heterocypris incongruens, but little is known of its response to fish kairomones or its interactions with potential competitors. We studied the demographic responses of S. vetulus fed the green alga Scenedesmus acutus, Microcystis cf. aeruginosa strain A, Microcystis cf. aeruginosa strain B, or Limnothrix sp. Experiments were conducted separately and together in the presence of Heterocypris incongruens and cichlid fish (Oreochromis) kairomones. A diet of Limnothrix sp. resulted in the lowest population growth rate (0.21±0.023 d^?1), while on diets of S. acutus or Microcystis, population growth was higher (0.30±0.009 d^?1). The presence of ostracods resulted in significantly higher growth rates of S. vetulus fed Limnothrix (0.33±0.01 d^?1), but not Microcystis or S. acutus. Regardless of the diet, the presence of fish kairomones resulted in significantly higher growth rates as compared with controls, particularly when ostracods were also present. Coexistence with ostracods may be beneficial to S. vetulus, particularly when food quality is poor.     

Reticulate evolution and sea‐level fluctuations together drove species diversification of slipper orchids (Paphiopedilum) in South‐East Asia

South‐East Asia covers four of the world's biodiversity hotspots, showing high species diversity and endemism. Owing to the successive expansion and contraction of distribution and the fragmentation by geographical barriers, the tropical flora greatly diversified in this region during the Tertiary, but the evolutionary tempo and mode of species diversity remain poorly investigated. Paphiopedilum, the largest genus of slipper orchids comprising nearly 100 species, is mainly distributed in South‐East Asia, providing an ideal system for exploring how plant species diversity was shaped in this region. Here, we investigated the evolutionary history of this genus with eight cpDNA regions and four low‐copy nuclear genes. Discordance between gene trees and network analysis indicates that reticulate evolution occurred in the genus. Ancestral area reconstruction suggests that vicariance and long‐distance dispersal together led to its current distribution. Diversification rate variation was detected and strongly correlated with the species diversity in subg. Paphiopedilum (~80 species). The shift of speciation rate in subg. Paphiopedilum was coincident with sea‐level fluctuations in the late Cenozoic, which could have provided ecological opportunities for speciation and created bridges or barriers for gene flow. Moreover, some other factors (e.g. sympatric distribution, incomplete reproductive barriers and clonal propagation) might also be advantageous for the formation and reproduction of hybrid species. In conclusion, our study suggests that the interplay of reticulate evolution and sea‐level fluctuations has promoted the diversification of the genus Paphiopedilum and sheds light into the evolution of Orchidaceae and the historical processes of plant species diversification in South‐East Asia.     

Genomic data reveal ancient microendemism in forest scorpions across the California Floristic Province

The California Floristic Province (CFP) in western North America is a globally significant biodiversity hotspot. Elucidating patterns of endemism and the historical drivers of this diversity has been an important challenge of comparative phylogeography for over two decades. We generated phylogenomic data using ddRADseq to examine genetic structure in Uroctonus forest scorpions, an ecologically restricted and dispersal‐limited organism widely distributed across the CFP north to the Columbia River. We coupled our genetic data with species distribution models (SDMs) to determine climatically suitable areas for Uroctonus both now and during the Last Glacial Maximum. Based on our analyses, Uroctonus is composed of two major genetic groups that likely diverged over 2 million years ago. Each of these groups itself contains numerous genetic groups that reveal a pattern of vicariance and microendemism across the CFP. Migration rates among these populations are low. SDMs suggest forest scorpion habitat has remained relatively stable over the last 21 000 years, consistent with the genetic data. Our results suggest tectonic plate rafting, mountain uplift, river drainage formation and climate‐induced habitat fragmentation have all likely played a role in the diversification of Uroctonus. The intricate pattern of genetic fragmentation revealed across a temporal continuum highlights the potential of low‐dispersing species to shed light on small‐scale patterns of biodiversity and the underlying processes that have generated this diversity in biodiversity hotspots.     

Molecular phylogeny of the diatom genera Amphora and Halamphora (Bacillariophyta) with a focus on morphological and ecological evolution

The taxonomic history of the diatom genus Amphora is one of a broad early morphological concept resulting in the inclusion of a diversity of taxa, followed by an extended period of revision and refinement. The introduction of molecular systematics has increased the pace of revision and has largely resolved the relationships between the major lineages, indicating homoplasy in the evolution of amphoroid symmetry. Within the two largest monophyletic lineages, the genus Halamphora and the now taxonomically refined genus Amphora, the intrageneric morphological and ecological relationships have yet to be explored within a phylogenetic framework. Critical among this is whether the range of morphological features exhibited within these diverse genera are reflective of evolutionary groupings or, as with many previously studied amphoroid features, are nonhomologous when examined phylogenetically. Presented here is a four‐marker molecular phylogeny that includes 31 taxa from the genus Amphora and 77 taxa from the genus Halamphora collected from fresh, brackish, and salt waters from coastal and inland habitats of the United States and Japan. These phylogenies illustrate complex patterns in the evolution of frustule morphology and ecology within the genera and the implications of this on the taxonomy, classification, and organization of the genera are discussed.     

Global warming may disproportionately affect larger adults in a predatory coral reef fish

Global warming is expected to reduce body sizes of ectothermic animals. Although the underlying mechanisms of size reductions remain poorly understood, effects appear stronger at latitudinal extremes (poles and tropics) and in aquatic rather than terrestrial systems. To shed light on this phenomenon, we examined the size dependence of critical thermal maxima (CTmax) and aerobic metabolism in a commercially important tropical reef fish, the leopard coral grouper (Plectropomus leopardus) following acclimation to current‐day (28.5 °C) vs. projected end‐of‐century (33 °C) summer temperatures for the northern Great Barrier Reef (GBR). CTmax declined from 38.3 to 37.5 °C with increasing body mass in adult fish (0.45–2.82 kg), indicating that larger individuals are more thermally sensitive than smaller conspecifics. This may be explained by a restricted capacity for large fish to increase mass‐specific maximum metabolic rate (MMR) at 33 °C compared with 28.5 °C. Indeed, temperature influenced the relationship between metabolism and body mass (0.02–2.38 kg), whereby the scaling exponent for MMR increased from 0.74 ± 0.02 at 28.5 °C to 0.79 ± 0.01 at 33 °C, and the corresponding exponents for standard metabolic rate (SMR) were 0.75 ± 0.04 and 0.80 ± 0.03. The increase in metabolic scaling exponents at higher temperatures suggests that energy budgets may be disproportionately impacted in larger fish and contribute to reduced maximum adult size. Such climate‐induced reductions in body size would have important ramifications for fisheries productivity, but are also likely to have knock‐on effects for trophodynamics and functioning of ecosystems.     

Evolution of putative barrier loci at an intermediate stage of speciation with gene flow in campions (Silene)

Understanding the origin of new species is a central goal in evolutionary biology. Diverging lineages often evolve highly heterogeneous patterns of genetic differentiation; however, the underlying mechanisms are not well understood. We investigated evolutionary processes governing genetic differentiation between the hybridizing campions Silene dioica (L.) Clairv. and S. latifolia Poiret. Demographic modelling indicated that the two species diverged with gene flow. The best‐supported scenario with heterogeneity in both migration rate and effective population size suggested that a small proportion of the loci evolved without gene flow. Differentiation (F _ST) and sequence divergence (d _XY) were correlated and both tended to peak in the middle of most linkage groups, consistent with reduced gene flow at highly differentiated loci. Highly differentiated loci further exhibited signatures of selection. In between‐species population pairs, isolation by distance was stronger for genomic regions with low between‐species differentiation than for highly differentiated regions that may contain barrier loci. Moreover, differentiation landscapes within and between species were only weakly correlated, suggesting that linked selection due to shared recombination and gene density landscapes is not the dominant determinant of genetic differentiation in these lineages. Instead, our results suggest that divergent selection shaped the genomic landscape of differentiation between the two Silene species, consistent with predictions for speciation in the face of gene flow.     

Glass in the water: Molecular phylogenetics and evolution of Indian glassy perchlets (Teleostei: Ambassidae)

Glassfishes of the family Ambassidae, comprising around 50 species, are distributed in the Indo‐West Pacific where they inhabit marine, estuarine, and freshwater ecosystems. We investigated for the first time the molecular phylogenetic and evolutionary relationships of this group using a combined dataset of mitochondrial and nuclear genes, particularly focusing on the taxa occurring in the Indian subcontinent. Results revealed that marine and freshwater genera of Ambassidae diverged during the Paleocene (~62 mya). The enigmatic monotypic genus Chanda is nested within the larger clade currently recognized as Parambassis, indicating its paraphyly. Based on cleared and stained osteological preparations and phylogenetic placement of Chanda nama, we hypothesize that the elongated and protruding lower jaw is an autapomorphic character that might have evolved for the lepidophagous habit of the species. The southern Indian species of Parambassis, Parambassis dayi, and Parambassis thomassi, which formed a monophyletic group, probably diverged from other species of Parambassis and Chanda nama around the Eocene (~42 mya) and can potentially be recognized as a distinct genus in view of the apomorphic characters such as the presence of serration on the ventral fringe of interopercle, densely serrated palatine and ectopterygoid, and the presence of more than 30 serrations along the lower preopercle and the posterior edge. Our analysis provides new insights into the evolution and phylogenetic relationships of glassy perchlets, including detailed relationships among the Indian species within this family.     

The genomic bases of morphological divergence and reproductive isolation driven by ecological speciation in Senecio (Asteraceae)

Ecological speciation, driven by adaptation to contrasting environments, provides an attractive opportunity to study the formation of distinct species, and the role of selection and genomic divergence in this process. Here, we focus on a particularly clear‐cut case of ecological speciation to reveal the genomic bases of reproductive isolation and morphological differences between closely related Senecio species, whose recent divergence within the last ~200 000 years was likely driven by the uplift of Mt. Etna (Sicily). These species form a hybrid zone, yet remain morphologically and ecologically distinct, despite active gene exchange. Here, we report a high‐density genetic map of the Senecio genome and map hybrid breakdown to one large and several small quantitative trait loci (QTL). Loci under diversifying selection cluster in three 5 cM regions which are characterized by a significant increase in relative (F_ST), but not absolute (d_XY), interspecific differentiation. They also correspond to some of the regions of greatest marker density, possibly corresponding to ‘cold‐spots’ of recombination, such as centromeres or chromosomal inversions. Morphological QTL for leaf and floral traits overlap these clusters. We also detected three genomic regions with significant transmission ratio distortion (TRD), possibly indicating accumulation of intrinsic genetic incompatibilities between these recently diverged species. One of the TRD regions overlapped with a cluster of high species differentiation, and another overlaps the large QTL for hybrid breakdown, indicating that divergence of these species may have occurred due to a complex interplay of ecological divergence and accumulation of intrinsic genetic incompatibilities.     

The accumulation of reproductive isolation in early stages of divergence supports a role for sexual selection

Models of speciation by sexual selection propose that male–female coevolution leads to the rapid evolution of behavioural reproductive isolation. Here, we compare the strength of behavioural isolation to ecological isolation, gametic incompatibility and hybrid inviability in a group of dichromatic stream fishes. In addition, we examine whether any of these individual barriers, or a combined measure of total isolation, is predicted by body shape differences, male colour differences, environmental differences or genetic distance. Behavioural isolation reaches the highest values of any barrier and is significantly greater than ecological isolation. No individual reproductive barrier is associated with any of the predictor variables. However, marginally significant relationships between male colour and body shape differences with ecological and behavioural isolation are discussed. Differences in male colour and body shape predict total reproductive isolation between species; hierarchical partitioning of these two variables' effects suggests a stronger role for male colour differences. Together, these results suggest an important role for divergent sexual selection in darter speciation but raise new questions about the mechanisms of sexual selection at play and the role of male nuptial ornaments.