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

Genetic structure, UV-vision, wing coloration and size coincide with colour polymorphism in Fabriciana adippe butterflies

Colour polymorphisms have long served as model systems in evolutionary studies and continue to inform about processes involved in the origin and dynamics of biodiversity. Modern sequencing tools allow for evaluating whether phenotypic differences between morphs reflect genetic differentiation rather than developmental plasticity, and for investigating whether polymorphisms represent intermediate stages of diversification towards speciation. We investigated phenotypic and genetic differentiation between two colour morphs of the butterfly Fabriciana adippe using a combination of ddRAD-sequencing and comparisons of body size, colour patterns and optical properties of bright wing spots. The silvery-spotted adippe form had larger and darker wings and reflected UV light, while the yellow cleodoxa form displayed more green scales and reflected very little UV, showcasing that they constitute distinct and alternative integrated phenotypes. Genomic analyses revealed genetic structuring according to source population, and to colour morph, suggesting that the phenotypic differentiation reflects evolutionary modifications. We report 17 outlier loci associated with colour morph, including ultraviolet-sensitive visual pigment (UVRh1), which is associated with intraspecific communication and mate choice in butterflies. Together with the demonstration that the wings of the adippe (but essentially not the cleodoxa) morph reflect UV light, that UV reflectance is higher in females than males and that morphs differ in wing size, this suggests that these colour morphs might represent genetically integrated phenotypes, possibly adapted to different microhabitats. We propose that non-random mating might contribute to the differentiation and maintenance of the polymorphism.

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.

Ecological speciation by sympatric host shifts in a clade of herbivorous sea slugs, with introgression and localized mitochondrial capture between species

Host shifting in insect-plant systems was historically important to the development of ecological speciation theory, yet surprisingly few studies have examined whether host shifting drives diversification of marine herbivores. When small-bodied consumers feed and also mate on a preferred host, disruptive selection can split a population into host races despite gene flow. Support for host shifts is notably lacking for invertebrates associated with macroalgae, where the scale of dispersal by planktonic larvae often far exceeds the grain of host patchiness, and adults are typically less specialized than terrestrial herbivores. Here, we present a candidate example of ecological speciation in a clade of sea slugs that primarily consume green algae in the genus Caulerpa, including highly invasive species. Ancestral character state reconstructions supported 'sea grapes' (C. racemosa, C. lentillifera) as the ancestral host for a tropical radiation of 12 Elysia spp., with one shift onto alternative Caulerpa spp. in the Indo-Pacific. A Caribbean radiation of three species included symaptric host shifts to Rhipocephalus brevicaulis in the ancestor of E. pratensis Ortea & Espinosa, 1996, and to C. prolifera in E. hamanni Krug, Vendetti & Valdes 2016, plus a niche expansion to a range of Caulerpa spp. in E. subornata Verrill, 1901. All three species are broadly sympatric across the Caribbean but are host-partitioned at a fine grain, and distinct by morphology and at nuclear loci. However, non-recombining mtDNA revealed a history of gene flow between E. pratensis and E. subornata: COI haplotypes from E. subornata were 10.4% divergent from E. pratensis haplotypes from four sites, but closely related to all E. pratensis haplotypes sampled from six Bahamian islands, indicating historical introgression and localized "mitochondrial capture." Disruptive selective likely fueled divergence and adaptation to distinct host environments, indicating ecological speciation may be an under-appreciated driver of diversification for marine herbivores as well as epibionts and other resource specialists.

The interaction of resource use and gene flow on the phenotypic divergence of benthic and pelagic morphs of Icelandic Arctic charr (Salvelinus alpinus)

Conceptual models of adaptive divergence and ecological speciation in sympatry predict differential resource use, phenotype-environment correlations, and reduced gene flow among diverging phenotypes. While these predictions have been assessed in past studies, connections among them have rarely been assessed collectively. We examined relationships among phenotypic, ecological, and genetic variation in Arctic charr (Salvelinus alpinus) from six Icelandic localities that have undergone varying degrees of divergence into sympatric benthic and pelagic morphs. We characterized morphological variation with geometric morphometrics, tested for differential resource use between morphs using stable isotopes, and inferred the amount of gene flow from single nucleotide polymorphisms. Analysis of stable isotopic signatures indicated that sympatric morphs showed similar difference in resource use across populations, likely arising from the common utilization of niche space within each population. Carbon isotopic signature was also a significant predictor of individual variation in body shape and size, suggesting that variation in benthic and pelagic resource use is associated with phenotypic variation. The estimated percentage of hybrids between sympatric morphs varied across populations (from 0% to 15.6%) but the majority of fish had genotypes (ancestry coefficients) characteristic of pure morphs. Despite evidence of reduced gene flow between sympatric morphs, we did not detect the expected negative relationship between divergence in resource use and gene flow. Three lakes showed the expected pattern, but morphs in the fourth showed no detectable hybridization and had relatively low differences in resource use between them. This coupled with the finding that resource use and genetic differentiation had differential effects on body shape variation across populations suggests that reproductive isolation maintains phenotypic divergence between benthic and pelagic morphs when the effects of resource use are relatively low. Our ability to assess relationships between phenotype, ecology, and genetics deepens our understanding of the processes underlying adaptive divergence in sympatry.

Species Radiations in the Sea: What the Flock?

Species flocks are proliferations of closely-related species, usually after colonization of depauperate habitat. These radiations are abundant on oceanic islands and in ancient freshwater lakes, but rare in marine habitats. This contrast is well documented in the Hawaiian Archipelago, where terrestrial examples include the speciose silverswords (sunflower family Asteraceae), Drosophila fruit flies, and honeycreepers (passerine birds), all derived from one or a few ancestral lineages. The marine fauna of Hawai'i is also the product of rare colonization events, but these colonizations usually yield only one species. Dispersal ability is key to understanding this evolutionary inequity. While terrestrial fauna rarely colonize between oceanic islands, marine fauna with pelagic larvae can make this leap in every generation. An informative exception is the marine fauna that lack a pelagic larval stage. These low-dispersal species emulate a "terrestrial" mode of reproduction (brooding, viviparity, crawl-away larvae), yielding marine species flocks in scattered locations around the world. Elsewhere, aquatic species flocks are concentrated in specific geographic settings, including the ancient lakes of Baikal (Siberia) and Tanganyika (eastern Africa), and Antarctica. These locations host multiple species flocks across a broad taxonomic spectrum, indicating a unifying evolutionary phenomenon. Hence marine species flocks can be singular cases that arise due to restricted dispersal or other intrinsic features, or they can be geographically clustered, promoted by extrinsic ecological circumstances. Here, we review and contrast intrinsic cases of species flocks in individual taxa, and extrinsic cases of geological/ecological opportunity, to elucidate the processes of species radiations.

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.

Host-symbiont coevolution, cryptic structure, and bleaching susceptibility, in a coral species complex (Scleractinia; Poritidae)

The 'species' is a key concept for conservation and evolutionary biology, yet the lines between population and species-level variation are often blurred, especially for corals. The 'Porites lobata species complex' consists of branching and mounding corals that form reefs across the Pacific. We used reduced representation meta-genomic sequencing to examine genetic relationships within this species complex and to identify candidate loci associated with colony morphology, cryptic genetic structure, and apparent bleaching susceptibility. We compared existing Porites data with bleached and unbleached colonies of the branching coral P. compressa collected in Kāne'ohe Bay Hawai'i during the 2015 coral bleaching event. Loci that mapped to coral, symbiont, and microbial references revealed genetic structure consistent with recent host-symbiont co-evolution. Cryptic genetic clades were resolved that previous work has associated with distance from shore, but no genetic structure was associated with bleaching. We identified many candidate loci associated with morphospecies, including candidate host and symbiont loci with fixed differences between branching and mounding corals. We also found many loci associated with cryptic genetic structure, yet relatively few loci associated with bleaching. Recent host-symbiont co-evolution and rapid diversification suggests that variation and therefore the capacity of these corals to adapt may be underappreciated.

Genetic Evidence for a Mixed Composition of the Genus Myoxocephalus (Cottoidei: Cottidae) Necessitates Generic Realignment

Sculpin fishes belonging to the family Cottidae represent a large and complex group, inhabiting a wide range of freshwater, brackish-water, and marine environments. Numerous studies based on analysis of their morphology and genetic makeup frequently provided controversial results. In the present work, we sequenced complete mitochondrial (mt) genomes and fragments of nuclear ribosomal DNA (rDNA) of the fourhorn sculpin Myoxocephalus quadricornis and some related cottids to increase the power of phylogenetic and taxonomic analyses of this complex fish group. A comparison of the My. quadricornis mt genomes obtained by us with other complete mt genomes available in GenBank has revealed a surprisingly low divergence (3.06 ± 0.12%) with Megalocottus platycephalus and, at the same time, a significantly higher divergence (7.89 ± 0.16%) with the species of the genus Myoxocephalus. Correspondingly, phylogenetic analyses have shown that My. quadricornis is clustered with Me. platycephalus but not with the Myoxocephalus species. Completely consistent patterns of divergence and tree topologies have been obtained based on nuclear rDNA. Thus, the multi-gene data in the present work indicates obvious contradictions in the relationships between the Myoxocephalus and Megalocottus species studied. An extensive phylogenetic analysis has provided evidence for a closer affinity of My. quadricornis with the species of the genus Megalocottus than with the species of the genus Myoxocephalus. A recombination analysis, along with the additional GenBank data, excludes introgression and/or incorrect taxonomic identification as the possible causative factors responsible for the observed closer affinity between the two species from different genera. The above facts necessitate realignment of the genera Myoxocephalus and Megalocottus. The genetic data supports the two recognized genera, Myoxocephalus and Megalocottus, but suggests changing their compositions through transferring My. quadricornis to the genus Megalocottus. The results of the present study resolve the relationships within a complex group of sculpin fishes and show a promising approach to phylogenetic systematics (as a key organizing principle in biodiversity research) for a better understanding of the taxonomy and evolution of fishes and for supplying relevant information to address various fish biodiversity conservation and management issues.

Does color matter? Molecular and ecological divergence in four sympatric color morphs of a coral reef fish

Non-sex-linked color polymorphism is common in animals and can be maintained in populations via balancing selection or, when under diversifying selection, can promote divergence. Despite their potential importance in ecological interactions and the evolution of biodiversity, their function and the mechanisms by which these polymorphisms are maintained are still poorly understood. Here, we combine field observations with life history and molecular data to compare four sympatric color morphs of the coral reef fish Paracirrhites forsteri (family Cirrhitidae) in the central Red Sea. Our findings verify that the color morphs are not sex-limited, inhabit the same reefs, and do not show clear signs of avoidance or aggression among them. A barcoding approach based on 1,276 bp of mitochondrial DNA could not differentiate the color morphs. However, when 36,769 SNPs were considered, we found low but significant population structure. Focusing on 1,121 F ST outliers, we recovered distinct population clusters that corresponded to shifts in allele frequencies with each color morph harboring unique alleles. Genetic divergence at these outlier loci is accompanied by differences in growth and marginal variation in microhabitat preference. Together, life history and molecular analysis suggest subtle divergence between the color morphs in this population, the causes for which remain elusive.

Decoding diversity in a coral reef fish species complex with restricted range using metagenomic sequencing of gut contents

AIM

Identification of the processes that generate and maintain species diversity within the same region can provide insight into biogeographic patterns at broader spatiotemporal scales. Hawkfishes in the genus Paracirrhites are a unique taxon to explore with respect to niche differentiation, exhibiting diagnostic differences in coloration, and an apparent center of distribution outside of the Indo-Malay-Philippine (IMP) biodiversity hotspot for coral reef fishes. Our aim is to use next-generation sequencing methods to leverage samples of a taxon at their center of maximum diversity to explore phylogenetic relationships and a possible mechanism of coexistence.

LOCATION

Flint Island, Southern Line Islands, Republic of Kiribati.

METHODS

A comprehensive review of museum records, the primary literature, and unpublished field survey records was undertaken to determine ranges for four "arc-eye" hawkfish species in the Paracirrhites species complex and a potential hybrid. Fish from four Paracirrhites species were collected from Flint Island in the Southern Line Islands, Republic of Kiribati. Hindgut contents were sequenced, and subsequent metagenomic analyses were used to assess the phylogenetic relatedness of the host fish, the microbiome community structure, and prey remains for each species.

RESULTS

Phylogenetic analyses conducted with recovered mitochondrial genomes revealed clustering of P. bicolor with P. arcatus and P. xanthus with P. nisus, which were unexpected on the basis of previous morphological work in this species complex. Differences in taxonomic composition of gut microbial communities and presumed prey remains indicate likely separation of foraging niches.

MAIN CONCLUSIONS

Our findings point toward previously unidentified relationships in this cryptic species complex at its proposed center of distribution. The three species endemic to the Polynesian province (P. nisus, P. xanthus, and P. bicolor) cluster separately from the more broadly distributed P. arcatus on the basis of relative abundance of metazoan sequences in the gut (presumed prey remains). Discordance between gut microbial communities and phylogeny of the host fish further reinforce the hypothesis of niche separation.

Patterns of morphological diversification in giant Berberomeloe blister beetles (Coleoptera: Meloidae) reveal an unexpected taxonomic diversity concordant with mtDNA phylogenetic structure

Delimiting species boundaries is a complex challenge usually hindered by overlooked morphological diversification or misinterpretation of geographically structured phenotypic variability. Independent molecular data are extremely useful to characterize and understand such morphological diversity. Morphological and molecular variability of the non-phoretic and apterous, widely distributed, giant blister beetles of the genus Berberomeloe, were investigated within and between lineages across most of the distributional range of the genus. We used two mtDNA gene fragments to characterize genetic variability and to produce a time-calibrated phylogeny of the genus. Our results reveal several mitochondrial lineages, allopatrically, parapatrically and sympatrically distributed. Most clades are not distinguishable between each other based on morphometrics. However, no morphometric overlap is observed between two closely related clades, one of them occurring in sympatry with a distantly congeneric species (B. insignis), suggesting that sympatry could trigger morphological diversification. Although most species share a morphometric space, they can be morphologically identified by a combination of easily observed characteristic qualitative features. Based on the concordance between mtDNA clades and morphological units, we describe six new species of Berberomeloe (B. castuo sp. nov., B. comunero sp. nov., B. indalo sp. nov, B. yebli sp. nov., B. payoyo sp. nov. and B. tenebrosus sp. nov.), revalidate two taxa (B. maculifrons comb. nov. and B. laevigatus comb. nov.) and redefine B. majalis.

The combined use of raw and phylogenetically independent methods of outlier detection uncovers genome-wide dynamics of local adaptation in a lizard

Local adaptation is a dynamic process by which different allele combinations are selected in different populations at different times, and whose genetic signature can be inferred by genome-wide outlier analyses. We combined gene flow estimates with two methods of outlier detection, one of them independent of population coancestry (CIOA) and the other one not (ROA), to identify genetic variants favored when ecology promotes phenotypic convergence. We analyzed genotyping-by-sequencing data from five populations of a lizard distributed over an environmentally heterogeneous range that has been changing since the split of eastern and western lineages ca. 3 mya. Overall, western lizards inhabit forest habitat and are unstriped, whereas eastern ones inhabit shrublands and are striped. However, one population (Lerma) has unstriped phenotype despite its eastern ancestry. The analysis of 73,291 SNPs confirmed the east-west division and identified nonoverlapping sets of outliers (12 identified by ROA and 9 by CIOA). ROA revealed ancestral adaptive variation in the uncovered outliers that were subject to divergent selection and differently fixed for eastern and western populations at the extremes of the environmental gradient. Interestingly, such variation was maintained in Lerma, where we found high levels of heterozygosity for ROA outliers, whereas CIOA uncovered innovative variants that were selected only there. Overall, it seems that both the maintenance of ancestral variation and asymmetric migration have counterbalanced adaptive lineage splitting in our model species. This scenario, which is likely promoted by a changing and heterogeneous environment, could hamper ecological speciation of locally adapted populations despite strong genetic structure between lineages.

Drift versus selection as drivers of phenotypic divergence at small spatial scales: The case of Belgjarskógur threespine stickleback

Divergence in phenotypic traits is facilitated by a combination of natural selection, phenotypic plasticity, gene flow, and genetic drift, whereby the role of drift is expected to be particularly important in small and isolated populations. Separating the components of phenotypic divergence is notoriously difficult, particularly for multivariate phenotypes. Here, we assessed phenotypic divergence of threespine stickleback (Gasterosteus aculeatus) across 19 semi-interconnected ponds within a small geographic region (~7.5 km2) using comparisons of multivariate phenotypic divergence (PST), neutral genetic (FST), and environmental (EST) variation. We found phenotypic divergence across the ponds in a suite of functionally relevant phenotypic traits, including feeding, defense, and swimming traits, and body shape (geometric morphometric). Comparisons of PSTs with FSTs suggest that phenotypic divergence is predominantly driven by neutral processes or stabilizing selection, whereas phenotypic divergence in defensive traits is in accordance with divergent selection. Comparisons of population pairwise PSTs with ESTs suggest that phenotypic divergence in swimming traits is correlated with prey availability, whereas there were no clear associations between phenotypic divergence and environmental difference in the other phenotypic groups. Overall, our results suggest that phenotypic divergence of these small populations at small geographic scales is largely driven by neutral processes (gene flow, drift), although environmental determinants (natural selection or phenotypic plasticity) may play a role.

Multilocus molecular systematics of the circumtropical reef-fish genus Abudefduf (Pomacentridae): history, geography and ecology of speciation

We investigated a pantropical sub-family and genus of damselfishes, the sergeant-majors (Pomacentridae: Abudefdufinae: Abudefduf), to identify the tempo and mechanisms of speciation in the lineage. We examined sequence capture data from 500 loci and 20 species, with multiple individuals sampled from across the geographic ranges of widespread species. Utilizing a maximum likelihood framework, as well as a time-calibrated Bayesian phylogeny, the following key questions are addressed: What is the historical tempo of speciation? What are the relative contributions of vicariant, peripatric and parapatric speciation to sergeant-major diversity? How is speciation related to major variation in trophic ecology? The approximately 20 species of sergeant-majors fall into three main lineages. The ancestral condition appears to be benthivory, which is predominant in two lineages comprising six species. The remaining species of sergeant-majors, of which there are at least 15, fall within a clade composed entirely of planktivores. This clade is sister to a benthivore clade that included one species, Abudefduf notatus, in transition to planktivory. Most speciation of sergeant-majors, which appeared ∼24 million years ago, occurred in the last 10 million years. Present distributional patterns indicate vicariant speciation precipitated by the closure of land barriers between both sides of the Atlantic and the Pacific, and the emergence of land between the Indian and Pacific Oceans. Within this backdrop, frequent oscillations in sea level over the last 10 million years also appear to have generated conditions suitable for both peripatric and vicariant speciation, and most speciation within the genus appears linked to these changes in sea level. Diversification within the genus has been concentrated in planktivorous seargeant-majors rather than benthivores. The root cause is unclear, but does not appear to be related to differences in dispersal potential, which is greater in the planktivorous species, due to the ability of their post-larval juveniles to raft with floating debris. This elevated speciation rate in planktivores and their propensity to form local endemics may reflect relaxation of selective pressures (e.g., on crypticity) that limit speciation in benthivorous sergeant-majors. Finally, our data allow us to clarify relationships of geminate sergeant-major species, indicating that there are subdivisions within the Atlantic for both benthivore and planktivore geminate pairs that may have misled previous studies.