Evaluating genomic divergence and parallelism in replicate ecomorphs from young and old cichlid adaptive radiations

Genomic architecture underlying morphological and physiological adaptation to high elevation in a songbird

Organisms often acquire physiological and morphological modifications to conquer ecological challenges when colonizing new environments which lead to their adaptive evolution. However, deciphering the genomic mechanism of ecological adaptation is difficult because ecological environments are often too complex for straightforward interpretation. Thus, we examined the adaptation of a widespread songbird-the rufous-capped babbler (Cyanoderma ruficeps)-to a relatively simple system: distinct environments across elevational gradients on the mountainous island of Taiwan. We focused on the genomic sequences of 43 birds from five populations to show that the Taiwan group split from its sister group in mainland China around 1-2 million years ago (Ma) and colonized the montane habitats of Taiwan at least twice around 0.03-0.22 Ma. The montane and lowland Taiwan populations diverged with gene flow between them, suggesting strong selection associated with different elevations. We found that the montane babblers had smaller beaks than the lowland ones, consistent with Allen's rule, and identified candidate genes-COL9A1 and SOX11-underlying the beak size changes. We also found that altitudinally divergent mutations were mostly located in noncoding regions and tended to accumulate in chromosomal inversions and autosomes. The altitudinally divergent mutations might regulate genes related to haematopoietic, metabolic, immune, auditory and vision functions, as well as cerebrum morphology and plumage development. The results reveal the genomic bases of morphological and physiological adaptation in this species to the low temperature, hypoxia and high UV light environment at high elevation. These findings improve our understanding of how ecological adaptation drives population divergence from the perspective of genomic architecture.

Phylogenomics of trophically diverse cichlids disentangles processes driving adaptive radiation and repeated trophic transitions

Cichlid fishes of the tribe Tropheini are a striking case of adaptive radiation, exemplifying multiple trophic transitions between herbivory and carnivory occurring in sympatry with other established cichlid lineages. Tropheini evolved highly specialized eco-morphologies to exploit similar trophic niches in different ways repeatedly and rapidly. To better understand the evolutionary history and trophic adaptations of this lineage, we generated a dataset of 532 targeted loci from 21 out of the 22 described Tropheini species. We resolved the Tropheini into seven monophyletic genera and discovered one to be polyphyletic. The polyphyletic genus, Petrochromis, represents three convergent origins of the algae grazing trophic specialization. This repeated evolution of grazing may have been facilitated by adaptive introgression as we found evidence for gene flow among algae grazing genera. We also found evidence of gene flow among algae browsing genera, but gene flow was restricted between herbivorous and carnivorous genera. Furthermore, we observed no evidence supporting a hybrid origin of this radiation. Our molecular evolutionary analyses suggest that opsin genes likely evolved in response to selection pressures associated with trophic ecology in the Tropheini. We found surprisingly little evidence of positive selection in coding regions of jaw-shaping genes in this trophically diverse lineage. This suggests low degrees of freedom for further change in these genes, and possibly a larger role for regulatory variation in driving jaw adaptations. Our study emphasizes Tropheini cichlids as an important model for studying the evolution of trophic specialization and its role in speciation.

Somewhere I belong: phylogeny and morphological evolution in a species-rich lineage of ectoparasitic flatworms infecting cichlid fishes

A substantial portion of biodiversity has evolved through adaptive radiation. However, the effects of explosive speciation on species interactions remain poorly understood. Metazoan parasites infecting radiating host lineages could improve our knowledge because of their intimate host relationships. Yet limited molecular, phenotypic and ecological data discourage multivariate analyses of evolutionary patterns and encourage the use of discrete characters. Here, we assemble new molecular, morphological and host range data widely inferred from a species-rich lineage of parasites (Cichlidogyrus, Platyhelminthes: Monogenea) infecting cichlid fishes to address data scarcity. We infer a multimarker (28S/18S rDNA, ITS1, COI mtDNA) phylogeny of 58 of 137 species and characterize major lineages through synapomorphies inferred from mapping morphological characters. We predict the phylogenetic position of species without DNA data through shared character states, a morphological phylogenetic analysis, and a classification analysis with support vector machines. Based on these predictions and a cluster analysis, we assess the systematic informativeness of continuous characters, search for continuous equivalents for discrete characters, and suggest new characters for morphological traits not analysed to date. We also model the attachment/reproductive organ and host range evolution using the data for 136 of 137 described species and multivariate phylogenetic comparative methods (PCMs). We show that discrete characters not only can mask phylogenetic signals, but also are key for characterizing species groups. Regarding the attachment organ morphology, a divergent evolutionary regime for at least one lineage was detected and a limited morphological variation indicates host and environmental parameters affecting its evolution. However, moderate success in predicting phylogenetic positions, and a low systematic informativeness and high multicollinearity of morphological characters call for a revaluation of characters included in species characterizations.

The repeated evolution of stripe patterns is correlated with body morphology in the adaptive radiations of East African cichlid fishes

Color patterns are often linked to the behavioral and morphological characteristics of an animal, contributing to the effectiveness of such patterns as antipredatory strategies. Species-rich adaptive radiations, such as the freshwater fish family Cichlidae, provide an exciting opportunity to study trait correlations at a macroevolutionary scale. Cichlids are also well known for their diversity and repeated evolution of color patterns and body morphology. To study the evolutionary dynamics between color patterns and body morphology, we used an extensive dataset of 461 species. A phylogenetic supertree of these species shows that stripe patterns evolved ~70 times independently and were lost again ~30 times. Moreover, stripe patterns show strong signs of correlated evolution with body elongation, suggesting that the stripes' effectiveness as antipredatory strategy might differ depending on the body shape. Using pedigree-based analyses, we show that stripes and body elongation segregate independently, indicating that the two traits are not genetically linked. Their correlation in nature is therefore likely maintained by correlational selection. Lastly, by performing a mate preference assay using a striped CRISPR-Cas9 mutant of a nonstriped species, we show that females do not differentiate between striped CRISPR mutant males and nonstriped wild-type males, suggesting that these patterns might be less important for species recognition and mate choice. In summary, our study suggests that the massive rates of repeated evolution of stripe patterns are shaped by correlational selection with body elongation, but not by sexual selection.

Effects of hybridization on pelvic morphology: A macaque model

Ancient DNA analyses have shown that interbreeding between hominin taxa occurred multiple times. Although admixture is often reflected in skeletal phenotype, the relationship between the two remains poorly understood, hampering interpretation of the hominin fossil record. Direct study of this relationship is often impossible due to the paucity of hominin fossils and difficulties retrieving ancient genetic material. Here, we use a sample of known ancestry hybrids between two closely related nonhuman primate taxa (Indian and Chinese Macaca mulatta) to investigate the effect of admixture on skeletal morphology. We focus on pelvic shape, which has potential fitness implications in hybrids, as mismatches between maternal pelvic and fetal cranial morphology are often fatal to mother and offspring. As the pelvis is also one of the skeletal regions that differs most between Homo sapiens and Neanderthals, investigating the pelvic consequences of interbreeding could be informative regarding the viability of their hybrids. We find that the effect of admixture in M. mulatta is small and proportional to the relatively small morphological difference between the parent taxa. Sexual dimorphism appears to be the main determinant of pelvic shape in M. mulatta. The lack of difference in pelvic shape between Chinese and Indian M. mulatta is in contrast to that between Neanderthals and H. sapiens, despite a similar split time (in generations) between the hybridizing pairs. Greater phenotypic divergence between hominins may relate to adaptations to disparate environments but may also highlight how the unique degree of cultural buffering in hominins allowed for greater neutral divergence. In contrast to some previous work identifying extreme morphologies in first- and second-generation hybrids, here the relationship between pelvic shape and admixture is linear. This linearity may be because most sampled animals have a multigenerational admixture history or because of relatively high constraints on the pelvis compared with other skeletal regions.

Manipulation of the Tyrosinase gene permits improved CRISPR/Cas editing and neural imaging in cichlid fish

Direct tests of gene function have historically been performed in a limited number of model organisms. The CRISPR/Cas system is species-agnostic, offering the ability to manipulate genes in a range of models, enabling insights into evolution, development, and physiology. Astatotilapia burtoni, a cichlid fish from the rivers and shoreline around Lake Tanganyika, has been extensively studied in the laboratory to understand evolution and the neural control of behavior. Here we develop protocols for the creation of CRISPR-edited cichlids and create a broadly useful mutant line. By manipulating the Tyrosinase gene, which is necessary for eumelanin pigment production, we describe a fast and reliable approach to quantify and optimize gene editing efficiency. Tyrosinase mutants also remove a major obstruction to imaging, enabling visualization of subdermal structures and fluorophores in situ. These protocols will facilitate broad application of CRISPR/Cas9 to studies of cichlids as well as other non-traditional model aquatic species.

Intercontinental genomic parallelism in multiple three-spined stickleback adaptive radiations

Parallelism, the evolution of similar traits in populations diversifying in similar conditions, provides strong evidence of adaptation by natural selection. Many studies of parallelism focus on comparisons of different ecotypes or contrasting environments, defined a priori, which could upwardly bias the apparent prevalence of parallelism. Here, we estimated genomic parallelism associated with components of environmental and phenotypic variation at an intercontinental scale across four freshwater adaptive radiations (Alaska, British Columbia, Iceland and Scotland) of the three-spined stickleback (Gasterosteus aculeatus). We combined large-scale biological sampling and phenotyping with restriction site associated DNA sequencing (RAD-Seq) data from 73 freshwater lake populations and four marine ones (1,380 fish) to associate genome-wide allele frequencies with continuous distributions of environmental and phenotypic variation. Our three main findings demonstrate that (1) quantitative variation in phenotypes and environments can predict genomic parallelism; (2) genomic parallelism at the early stages of adaptive radiations, even at large geographic scales, is founded on standing variation; and (3) similar environments are a better predictor of genome-wide parallelism than similar phenotypes. Overall, this study validates the importance and predictive power of major phenotypic and environmental factors likely to influence the emergence of common patterns of genomic divergence, providing a clearer picture than analyses of dichotomous phenotypes and environments.

Different Sources of Allelic Variation Drove Repeated Color Pattern Divergence in Cichlid Fishes

The adaptive radiations of East African cichlid fish in the Great Lakes Victoria, Malawi, and Tanganyika are well known for their diversity and repeatedly evolved phenotypes. Convergent evolution of melanic horizontal stripes has been linked to a single locus harboring the gene agouti-related peptide 2 (agrp2). However, where and when the causal variants underlying this trait evolved and how they drove phenotypic divergence remained unknown. To test the alternative hypotheses of standing genetic variation versus de novo mutations (independently originating in each radiation), we searched for shared signals of genomic divergence at the agrp2 locus. Although we discovered similar signatures of differentiation at the locus level, the haplotypes associated with stripe patterns are surprisingly different. In Lake Malawi, the highest associated alleles are located within and close to the 5' untranslated region of agrp2 and likely evolved through recent de novo mutations. In the younger Lake Victoria radiation, stripes are associated with two intronic regions overlapping with a previously reported cis-regulatory interval. The origin of these segregating haplotypes predates the Lake Victoria radiation because they are also found in more basal riverine and Lake Kivu species. This suggests that both segregating haplotypes were present as standing genetic variation at the onset of the Lake Victoria adaptive radiation with its more than 500 species and drove phenotypic divergence within the species flock. Therefore, both new (Lake Malawi) and ancient (Lake Victoria) allelic variation at the same locus fueled rapid and convergent phenotypic evolution.

Genetic Variation and Hybridization in Evolutionary Radiations of Cichlid Fishes

Evolutionary radiations are responsible for much of the variation in biodiversity across taxa. Cichlid fishes are well known for spectacular evolutionary radiations, as they have repeatedly evolved into large and phenotypically diverse arrays of species. Cichlid genomes carry signatures of past events and, at the same time, are the substrate for ongoing evolution. We survey genome-wide data and the available literature covering 438 cichlid populations (412 species) across multiple radiations to synthesize information about patterns and sharing of genetic variation. Nucleotide diversity within species is low in cichlids, with 92% of surveyed populations having less diversity than the median value found in other vertebrates. Divergence within radiations is also low, and a large proportion of variation is shared among species due to incomplete lineage sorting and widespread hybridization. Population genetics therefore provides a suitable conceptual framework for evolutionary genomic studies of cichlid radiations. We focus in detail on the roles of hybridization in shaping the patterns of genetic variation and in promoting cichlid diversification.

Ecological speciation in European whitefish is driven by a large-gaped predator

Lake-dwelling fish that form species pairs/flocks characterized by body size divergence are important model systems for speciation research. Although several sources of divergent selection have been identified in these systems, their importance for driving the speciation process remains elusive. A major problem is that in retrospect, we cannot distinguish selection pressures that initiated divergence from those acting later in the process. To address this issue, we studied the initial stages of speciation in European whitefish (Coregonus lavaretus) using data from 358 populations of varying age (26-10,000 years). We find that whitefish speciation is driven by a large-growing predator, the northern pike (Esox lucius). Pike initiates divergence by causing a largely plastic differentiation into benthic giants and pelagic dwarfs: ecotypes that will subsequently develop partial reproductive isolation and heritable differences in gill raker number. Using an eco-evolutionary model, we demonstrate how pike's habitat specificity and large gape size are critical for imposing a between-habitat trade-off, causing prey to mature in a safer place or at a safer size. Thereby, we propose a novel mechanism for how predators may cause dwarf/giant speciation in lake-dwelling fish species.

Ancestral Hybridization Facilitated Species Diversification in the Lake Malawi Cichlid Fish Adaptive Radiation

The adaptive radiation of cichlid fishes in East African Lake Malawi encompasses over 500 species that are believed to have evolved within the last 800,000 years from a common founder population. It has been proposed that hybridization between ancestral lineages can provide the genetic raw material to fuel such exceptionally high diversification rates, and evidence for this has recently been presented for the Lake Victoria region cichlid superflock. Here, we report that Lake Malawi cichlid genomes also show evidence of hybridization between two lineages that split 3-4 Ma, today represented by Lake Victoria cichlids and the riverine Astatotilapia sp. "ruaha blue." The two ancestries in Malawi cichlid genomes are present in large blocks of several kilobases, but there is little variation in this pattern between Malawi cichlid species, suggesting that the large-scale mosaic structure of the genomes was largely established prior to the radiation. Nevertheless, tens of thousands of polymorphic variants apparently derived from the hybridization are interspersed in the genomes. These loci show a striking excess of differentiation across ecological subgroups in the Lake Malawi cichlid assemblage, and parental alleles sort differentially into benthic and pelagic Malawi cichlid lineages, consistent with strong differential selection on these loci during species divergence. Furthermore, these loci are enriched for genes involved in immune response and vision, including opsin genes previously identified as important for speciation. Our results reinforce the role of ancestral hybridization in explosive diversification by demonstrating its significance in one of the largest recent vertebrate adaptive radiations.

Evolutionary Dynamics of Structural Variation at a Key Locus for Color Pattern Diversification in Cichlid Fishes

Color patterns in African cichlid fishes vary spectacularly. Although phylogenetic analysis showed already 30 years ago that many color patterns evolved repeatedly in these adaptive radiations, only recently have we begun to understand the genomic basis of color variation. Horizontal stripe patterns evolved and were lost several times independently across the adaptive radiations of Lake Victoria, Malawi, and Tanganyika and regulatory evolution of agouti-related peptide 2 (agrp2/asip2b) has been linked to this phenotypically labile trait. Here, we asked whether the agrp2 locus exhibits particular characteristics that facilitate divergence in color patterns. Based on comparative genomic analyses, we discovered several recent duplications, insertions, and deletions. Interestingly, one of these events resulted in a tandem duplication of the last exon of agrp2. The duplication likely precedes the East African radiations that started 8-12 Ma, is not fixed within any of the radiations, and is found to vary even within some species. Moreover, we also observed variation in copy number (two to five copies) and secondary loss of the duplication, illustrating a surprising dynamic at this locus that possibly promoted functional divergence of agrp2. Our work suggests that such instances of exon duplications are a neglected mechanism potentially involved in the repeated evolution and diversification that deserves more attention.

Genomics of rapid ecological divergence and parallel adaptation in four tidal marsh sparrows

Theory suggests that different taxa having colonized a similar, challenging environment will show parallel or lineage-specific adaptations to shared selection pressures, but empirical examples of parallel evolution in independent taxa are exceedingly rare. We employed comparative genomics to identify parallel and lineage-specific responses to selection within and among four species of North American sparrows that represent four independent, post-Pleistocene colonization events by an ancestral, upland subspecies and a derived salt marsh specialist. We identified multiple cases of parallel adaptation in these independent comparisons following salt marsh colonization, including selection of 12 candidate genes linked to osmoregulation. In addition to detecting shared genetic targets of selection across multiple comparisons, we found many novel, species-specific signatures of selection, including evidence of selection of loci associated with both physiological and behavioral mechanisms of osmoregulation. Demographic reconstructions of all four species highlighted their recent divergence and small effective population sizes, as expected given their rapid radiation into saline environments. Our results highlight the interplay of both shared and lineage-specific selection pressures in the colonization of a biotically and abiotically challenging habitat and confirm theoretical expectations that steep environmental clines can drive repeated and rapid evolutionary diversification in birds.

Recent Secondary Contacts, Linked Selection, and Variable Recombination Rates Shape Genomic Diversity in the Model Species Anolis carolinensis

Gaining a better understanding on how selection and neutral processes affect genomic diversity is essential to gain better insights into the mechanisms driving adaptation and speciation. However, the evolutionary processes affecting variation at a genomic scale have not been investigated in most vertebrate lineages. Here, we present the first population genomics survey using whole genome resequencing in the green anole (Anolis carolinensis). Anoles have been intensively studied to understand mechanisms underlying adaptation and speciation. The green anole in particular is an important model to study genome evolution. We quantified how demography, recombination, and selection have led to the current genetic diversity of the green anole by using whole-genome resequencing of five genetic clusters covering the entire species range. The differentiation of green anole's populations is consistent with a northward expansion from South Florida followed by genetic isolation and subsequent gene flow among adjacent genetic clusters. Dispersal out-of-Florida was accompanied by a drastic population bottleneck followed by a rapid population expansion. This event was accompanied by male-biased dispersal and/or selective sweeps on the X chromosome. We show that the interaction between linked selection and recombination is the main contributor to the genomic landscape of differentiation in the anole genome.

Local ancestry analysis reveals genomic convergence in extremophile fishes

The molecular basis of convergent phenotypes is often unknown. However, convergence at a genomic level is predicted when there are large population sizes, gene flow among diverging lineages or strong genetic constraints. We used whole-genome resequencing to investigate genomic convergence in fishes ( Poecilia spp.) that have repeatedly colonized hydrogen sulfide (H₂S)-rich environments in Mexico. We identified genomic similarities in both single nucleotide polymorphisms (SNPs) and structural variants (SVs) among independently derived sulfide spring populations, with approximately 1.2% of the genome being shared among sulfidic ecotypes. We compared these convergent genomic regions to candidate genes for H₂S adaptation identified from transcriptomic analyses and found that a significant proportion of these candidate genes (8%) were also in regions where sulfidic individuals had similar SNPs, while only 1.7% were in regions where sulfidic individuals had similar SVs. Those candidate genes included genes involved in sulfide detoxification, the electron transport chain (the main toxicity target of H₂S) and other processes putatively important for adaptation to sulfidic environments. Regional genomic similarity across independent populations exposed to the same source of selection is consistent with selection on standing variation or introgression of adaptive alleles across divergent lineages. However, combined with previous analyses, our data also support that adaptive changes in mitochondrially encoded subunits arose independently via selection on de novo mutations. Pressing questions remain on what conditions ultimately facilitate the independent rise of adaptive alleles at the same loci in separate populations, and thus, the degree to which evolution is repeatable or predictable. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Genomic insights into the vulnerability of sympatric whitefish species flocks

The erosion of habitat heterogeneity can reduce species diversity directly but can also lead to the loss of distinctiveness of sympatric species through speciation reversal. We know little about changes in genomic differentiation during the early stages of these processes, which can be mediated by anthropogenic perturbation. Here, we analyse three sympatric whitefish species (Coregonus spp) sampled across two neighbouring and connected Swiss pre-alpine lakes, which have been differentially affected by anthropogenic eutrophication. Our data set comprises 16,173 loci genotyped across 138 whitefish using restriction-site associated DNA sequencing (RADseq). Our analysis suggests that in each of the two lakes, the population of a different, but ecologically similar, whitefish species declined following a recent period of eutrophication. Genomic signatures consistent with hybridization are more pronounced in the more severely impacted lake. Comparisons between sympatric pairs of whitefish species with contrasting ecology, where one is shallow benthic and the other one more profundal pelagic, reveal genomic differentiation that is largely correlated along the genome, while differentiation is uncorrelated between pairs of allopatric provenance with similar ecology. We identify four genomic loci that provide evidence of parallel divergent adaptation between the shallow benthic species and the two different more profundal species. Functional annotations available for two of those loci are consistent with divergent ecological adaptation. Our genomic analysis indicates the action of divergent natural selection between sympatric whitefish species in pre-alpine lakes and reveals the vulnerability of these species to anthropogenic alterations of the environment and associated adaptive landscape.

(Non)Parallel Evolution

Parallel evolution across replicate populations has provided evolutionary biologists with iconic examples of adaptation. When multiple populations colonize seemingly similar habitats, they may evolve similar genes, traits, or functions. Yet, replicated evolution in nature or in the laboratory often yields inconsistent outcomes: Some replicate populations evolve along highly similar trajectories, whereas other replicate populations evolve to different extents or in distinct directions. To understand these heterogeneous outcomes, biologists are increasingly treating parallel evolution not as a binary phenomenon but rather as a quantitative continuum ranging from parallel to nonparallel. By measuring replicate populations’ positions along this (non)parallel continuum, we can test hypotheses about evolutionary and ecological factors that influence the extent of repeatable evolution. We review evidence regarding the manifestation of (non)parallel evolution in the laboratory, in natural populations, and in applied contexts such as cancer. We enumerate the many genetic, ecological, and evolutionary processes that contribute to variation in the extent of parallel evolution.

Phylogenomics uncovers early hybridization and adaptive loci shaping the radiation of Lake Tanganyika cichlid fishes

Lake Tanganyika is the oldest and phenotypically most diverse of the three East African cichlid fish adaptive radiations. It is also the cradle for the younger parallel haplochromine cichlid radiations in Lakes Malawi and Victoria. Despite its evolutionary significance, the relationships among the main Lake Tanganyika lineages remained unresolved, as did the general timescale of cichlid evolution. Here, we disentangle the deep phylogenetic structure of the Lake Tanganyika radiation using anchored phylogenomics and uncover hybridization at its base, as well as early in the haplochromine radiation. This suggests that hybridization might have facilitated these speciation bursts. Time-calibrated trees support that the radiation of Tanganyika cichlids coincided with lake formation and that Gondwanan vicariance concurred with the earliest splits in the cichlid family tree. Genes linked to key innovations show signals of introgression or positive selection following colonization of lake habitats and species' dietary adaptations are revealed as major drivers of colour vision evolution. These findings shed light onto the processes shaping the evolution of adaptive radiations.

Genomic divergence between Spanish Littorina saxatilis ecotypes unravels limited admixture and extensive parallelism associated with population history

The rough periwinkle, Littorina saxatilis, is a model system for studying parallel ecological speciation in microparapatry. Phenotypically parallel wave-adapted and crab-adapted ecotypes that hybridize within the middle shore are replicated along the northwestern coast of Spain and have likely arisen from two separate glacial refugia. We tested whether greater geographic separation corresponding to reduced opportunity for contemporary or historical gene flow between parallel ecotypes resulted in less parallel genomic divergence. We sequenced double-digested restriction-associated DNA (ddRAD) libraries from individual snails from upper, mid, and low intertidal levels of three separate sites colonized from two separate refugia. Outlier analysis of 4256 SNP markers identified 34.4% sharing of divergent loci between two geographically close sites; however, these sites each shared only 9.9%-15.1% of their divergent loci with a third more-distant site. STRUCTURE analysis revealed that genotypes from only three of 166 phenotypically intermediate mid-shore individuals appeared to result from recent hybridization, suggesting that hybrids cannot be reliably identified using shell traits. Hierarchical AMOVA indicated that the primary source of genomic differentiation was geographic separation, but also revealed greater similarity of the same ecotype across the two geographically close sites than previously estimated with dominant markers. These results from a model system for ecological speciation suggest that genomic parallelism is affected by the opportunity for historical or contemporary gene flow between populations.

Experimental evidence for rapid genomic adaptation to a new niche in an adaptive radiation

A substantial part of biodiversity is thought to have arisen from adaptive radiations in which one lineage rapidly diversified into multiple lineages adapted to many different niches. However, selection and drift reduce genetic variation during adaptation to new niches and may thus prevent or slow down further niche shifts. We tested whether rapid adaptation is still possible from a highly derived ecotype in the adaptive radiation of threespine stickleback on the Haida Gwaii archipelago, Western Canada. In a 19-years selection experiment, we let giant stickleback from a large blackwater lake evolve in a small clearwater pond without vertebrate predators. 56 whole genomes from the experiment and 26 natural populations revealed that adaptive genomic change was rapid in many small genomic regions and encompassed 75% of the adaptive genomic change between 12,000 years old ecotypes. Adaptive genomic change was as fast as phenotypic change in defence and trophic morphology and both were largely parallel between the short-term selection experiment and long-term natural adaptive radiation. Our results show that functionally relevant standing genetic variation can persist in derived adaptive radiation members, allowing adaptive radiations to unfold very rapidly.

A genomic map of clinal variation across the European rabbit hybrid zone

Speciation is a process proceeding from weak to complete reproductive isolation. In this continuum, naturally hybridizing taxa provide a promising avenue for revealing the genetic changes associated with the incipient stages of speciation. To identify such changes between two subspecies of rabbits that display partial reproductive isolation, we studied patterns of allele frequency change across their hybrid zone using whole-genome sequencing. To connect levels and patterns of genetic differentiation with phenotypic manifestations of subfertility in hybrid rabbits, we further investigated patterns of gene expression in testis. Geographic cline analysis revealed 253 regions characterized by steep changes in allele frequency across their natural region of contact. This catalog of regions is likely to be enriched for loci implicated in reproductive barriers and yielded several insights into the evolution of hybrid dysfunction in rabbits: (i) incomplete reproductive isolation is likely governed by the effects of many loci, (ii) protein-protein interaction analysis suggest that genes within these loci interact more than expected by chance, (iii) regulatory variation is likely the primary driver of incompatibilities, and (iv) large chromosomal rearrangements appear not to be a major mechanism underlying incompatibilities or promoting isolation in the face of gene flow. We detected extensive misregulation of gene expression in testis of hybrid males, but not a statistical overrepresentation of differentially expressed genes in candidate regions. Our results also did not support an X chromosome-wide disruption of expression as observed in mice and cats, suggesting variation in the mechanistic basis of hybrid male reduced fertility among mammals.

Geographic isolation facilitates the evolution of reproductive isolation and morphological divergence

Geographic isolation is known to contribute to divergent evolution, resulting in unique phenotypes. Oftentimes morphologically distinct populations are found to be interfertile while reproductive isolation is found to exist within nominal morphological species revealing the existence of cryptic species. These disparities can be difficult to predict or explain especially when they do not reflect an inferred history of common ancestry which suggests that environmental factors affect the nature of ecological divergence. A series of laboratory experiments and observational studies were used to address what role biogeographic factors may play in the ecological divergence of Hyalella amphipods. It was found that geographic isolation plays a key role in the evolution of reproductive isolation and divergent morphology and that divergence cannot be explained by molecular genetic variation.

Molecular ecology studies of species radiations: current research gaps, opportunities and challenges

Understanding the drivers and limits of species radiations is a crucial goal of evolutionary genetics and molecular ecology, yet research on this topic has been hampered by the notorious difficulty of connecting micro- and macroevolutionary approaches to studying the drivers of diversification. To chart the current research gaps, opportunities and challenges of molecular ecology approaches to studying radiations, we examine the literature in the journal Molecular Ecology and revisit recent high-profile examples of evolutionary genomic research on radiations. We find that available studies of radiations are highly unevenly distributed among taxa, with many ecologically important and species-rich organismal groups remaining severely understudied, including arthropods, plants and fungi. Most studies employed molecular methods suitable over either short or long evolutionary time scales, such as microsatellites or restriction site-associated DNA sequencing (RAD-seq) in the former case and conventional amplicon sequencing of organellar DNA in the latter. The potential of molecular ecology studies to address and resolve patterns and processes around the species level in radiating groups of taxa is currently limited primarily by sample size and a dearth of information on radiating nuclear genomes as opposed to organellar ones. Based on our literature survey and personal experience, we suggest possible ways forward in the coming years. We touch on the potential and current limitations of whole-genome sequencing (WGS) in studies of radiations. We suggest that WGS and targeted ('capture') resequencing emerge as the methods of choice for scaling up the sampling of populations, species and genomes, including currently understudied organismal groups and the genes or regulatory elements expected to matter most to species radiations.

Gene flow, ancient polymorphism, and ecological adaptation shape the genomic landscape of divergence among Darwin's finches

Genomic comparisons of closely related species have identified “islands” of locally elevated sequence divergence. Genomic islands may contain functional variants involved in local adaptation or reproductive isolation and may therefore play an important role in the speciation process. However, genomic islands can also arise through evolutionary processes unrelated to speciation, and examination of their properties can illuminate how new species evolve. Here, we performed scans for regions of high relative divergence (F_ST) in 12 species pairs of Darwin's finches at different genetic distances. In each pair, we identify genomic islands that are, on average, elevated in both relative divergence (F_ST) and absolute divergence (d_XY). This signal indicates that haplotypes within these genomic regions became isolated from each other earlier than the rest of the genome. Interestingly, similar numbers of genomic islands of elevated d_XY are observed in sympatric and allopatric species pairs, suggesting that recent gene flow is not a major factor in their formation. We find that two of the most pronounced genomic islands contain the ALX1 and HMGA2 loci, which are associated with variation in beak shape and size, respectively, suggesting that they are involved in ecological adaptation. A subset of genomic island regions, including these loci, appears to represent anciently diverged haplotypes that evolved early during the radiation of Darwin's finches. Comparative genomics data indicate that these loci, and genomic islands in general, have exceptionally low recombination rates, which may play a role in their establishment.

Ancient hybridization fuels rapid cichlid fish adaptive radiations

Understanding why some evolutionary lineages generate exceptionally high species diversity is an important goal in evolutionary biology. Haplochromine cichlid fishes of Africa's Lake Victoria region encompass >700 diverse species that all evolved in the last 150,000 years. How this 'Lake Victoria Region Superflock' could evolve on such rapid timescales is an enduring question. Here, we demonstrate that hybridization between two divergent lineages facilitated this process by providing genetic variation that subsequently became recombined and sorted into many new species. Notably, the hybridization event generated exceptional allelic variation at an opsin gene known to be involved in adaptation and speciation. More generally, differentiation between new species is accentuated around variants that were fixed differences between the parental lineages, and that now appear in many new combinations in the radiation species. We conclude that hybridization between divergent lineages, when coincident with ecological opportunity, may facilitate rapid and extensive adaptive radiation.

Cichlids: A Host of Opportunities for Evolutionary Parasitology

Thanks to high species diversity and a broad range of speciation mechanisms, cichlid fishes represent a textbook model in evolutionary biology. They are also of substantial economic value. Despite this importance, cichlid parasites remain understudied, although some are more diverse than their hosts. They may offer important insights into cichlid evolution and the evolution of host-parasite interactions. We review five major lines of research conducted on cichlid parasites so far: the study of parasite diversity and speciation; the role of parasites in cichlid diversification; the evolutionary ecology of host specificity; historical biogeography; and biological invasions. We call for more research in these areas and suggest approaches to valorise the potential that cichlid parasites hold for the study of evolutionary parasitology.

Rapid and high throughput molecular identification of diverse mosquito species by high resolution melting analysis

Mosquitoes are a diverse group of invertebrates, with members that are among the most important vectors of diseases. The correct identification of mosquitoes is paramount to the control of the diseases that they transmit. However, morphological techniques depend on the quality of the specimen and often unavailable taxonomic expertise, which may still not be able to distinguish mosquitoes among species complexes (sibling and cryptic species). High resolution melting (HRM) analyses, a closed-tube, post-polymerase chain reaction (PCR) method used to identify variations in nucleic acid sequences, has been used to differentiate species within the Anopheles gambiae and Culex pipiens complexes. We validated the use of PCR-HRM analyses to differentiate species within Anopheles and within each of six genera of culicine mosquitoes, comparing primers targeting cytochrome b ( cyt b), NADH dehydrogenase subunit 1 (ND1), intergenic spacer region (IGS) and cytochrome c oxidase subunit 1 ( COI) gene regions. HRM analyses of amplicons from all the six primer pairs successfully differentiated two or more mosquito species within one or more genera ( Aedes ( Ae. vittatus from Ae. metallicus), Culex ( Cx. tenagius from Cx. antennatus, Cx. neavei from Cx. duttoni, cryptic Cx. pipiens species), Anopheles ( An. gambiae s.s. from An. arabiensis) and Mansonia ( Ma. africana from Ma. uniformis)) based on their HRM profiles. However, PCR-HRM could not distinguish between species within Aedeomyia ( Ad. africana and Ad. furfurea), Mimomyia ( Mi. hispida and Mi. splendens) and Coquillettidia ( Cq. aurites, Cq. chrysosoma, Cq. fuscopennata, Cq. metallica, Cq. microannulatus, Cq. pseudoconopas and Cq. versicolor) genera using any of the primers. The IGS and COI barcode region primers gave the best and most definitive separation of mosquito species among anopheline and culicine mosquito genera, respectively, while the other markers may serve to confirm identifications of closely related sub-species. This approach can be employed for rapid identification of mosquitoes.