Three problems in the genetics of speciation by selection

Floral phenotypic divergence and genomic insights in an Ophrys orchid: unraveling early speciation processes

Adaptive radiation in Ophrys orchids leads to complex floral phenotypes that vary in scent, color and shape. Using a novel pipeline to quantify these phenotypes, we investigated trait divergence at early stages of speciation in six populations of Ophrys aveyronensis experiencing recent allopatry. By integrating different genetic/genomic techniques, we investigated: variation and integration of floral components (scent, color and shape); phenotypes and genomic regions under divergent selection; and the genomic bases of trait variation. We identified a large genomic island of divergence, likely associated with phenotypic variation in particular in floral odor. We detected potential divergent selection on macular color, while stabilizing selection was suspected on floral morphology and for several volatile olfactive compounds. We also identified candidate genes involved in anthocyanin and in steroid biosynthesis pathways associated with standing genetic variation in color and odor. This study sheds light on early differentiation in Ophrys, revealing patterns that often become invisible over time, that is the geographic mosaic of traits under selection and the early appearance of strong genomic divergence. It also supports a crucial genomic region for future investigation and highlights the value of a multifaceted approach in unraveling speciation within taxa with large genomes.

Genomic Architecture Underlying the Striking Colour Variation in the Presence of Gene Flow for the Guinan Toad-Headed Lizard

How divergence occurs between closely related organisms in the absence of geographic barriers to gene flow stands as one of the long-standing questions in evolutionary biology. Previous studies suggested that the interplay between selection, gene flow and recombination strongly affected the process of divergence with gene flow. However, the extent to which these forces interact to drive divergence remains largely ambiguous. Guinan toad-headed lizards (Phrynocephalus guinanensis) in the Mugetan Desert exhibit striking colour differences from lizards outside the desert and provide an excellent model to address this question. Through extensive sampling and whole genome sequencing, we obtained genotypes for 191 samples from 14 populations inside and outside the desert. Despite the colour differences, continuous and asymmetric gene flow was detected across the desert border. More importantly, 273 highly diverged regions (HDRs) were identified between them, accounting only for 0.47% of the genome but widely distributed across 20 (out of the total 24) chromosomes. Strong signatures of selection were identified in HDRs, and local recombination rates were repressed. Furthermore, five HDRs exhibited significantly higher divergence, which contained key genes associated with crucial functions in animal coloration, including pteridine and melanocyte pigmentation. Genes related to retinal cells and steroid hormones were identified in other HDRs, which might have also contributed to the formation of colour variation in the presence of gene flow. This study provided novel insights into the understanding of the evolutionary mechanisms of genetic divergence in the presence of gene flow.

Phenological and morphological variations of Oryza rufipogon and O. nivara in Sri Lanka and their evolutionary implications

Phenological and morphological variation are widely viewed as a pivotal driver of ecological adaptation and speciation. Here, we investigate variation patterns of flowering phenology and morphological traits within and between O. rufipogon and O. nivara populations in Sri Lanka by incorporating the in situ observation in natural habitats and manipulative experiments in the common gardens. First, we observed varying degrees of phenological variation under different temporal and spatial conditions, suggesting that flowering phenology of two Oryza species varied depending on both environments and management practices. Particularly, the Sri Lankan O. nivara exhibits high plasticity in flowering phenology, implying that O. nivara might not be an annual in the strict sense. Second, the observation that flowering time of the two species overlapped suggests that the primary factor to maintain the species divergence in Sri Lanka may not be flowering time but rather environments. Third, our selection analysis suggests that interspecific divergence in the traits related to reproduction and habitat preference is adaptive and most likely driven by natural selection. Together, our case study on the Sri Lankan O. rufipogon and O. nivara enhances the understanding of the roles of phenotypic plasticity and environmental factors in the processes of adaptation and speciation.

The Genetic Architecture of Local Adaptation and Reproductive Character Displacement in Scutiger boulengeri Complex (Anura: Megophryidae)

Speciation is a continuous process driven by barriers to gene flow. Based on genome-wide SNPs (single nucleotide polymorphisms) of 190 toads from 31 sampling sites of Scutiger boulengeri complex, we found evidence for monophyly which represented a continuous speciation process of at least six lineages in S. boulengeri, which radiated and exhibited varying degrees of divergence and gene flow. The SNP-based phylogenetic tree was largely discordant with the multilocus mitochondrial tree (i.e., S. mammatus and S. glandulatus nested in the lineages of S. boulengeri) published before. The Min Mountains (MM) and Qinghai-Tibet Plateau (QTP) lineages differ fundamentally in habitat (i.e., elevation) and morphology (i.e., SVL), we detected signatures of potential high-altitude and cold adaptation genes in QTP (vs. MM). We found the evidence of reproductive trait disparity (i.e., SVL and nuptial pads) is key to promoting sympatric rather than allopatric species pairs. In addition, we identified selection signals for genes related to sympatric character displacement, genes linked to obesity-related traits, nuptial spines morphology and enlarged chest nuptial pads in S. mammatus (vs. QTP group of S. boulengeri). Our study provided new insight and paradigm for a varied speciation pattern from local adaptation of allopatry to sympatric character displacement in the S. boulengeri complex.

The Role of Hybridization in Species Formation and Persistence

Hybridization, or interbreeding between different taxa, was traditionally considered to be rare and to have a largely detrimental impact on biodiversity, sometimes leading to the breakdown of reproductive isolation and even to the reversal of speciation. However, modern genomic and analytical methods have shown that hybridization is common in some of the most diverse clades across the tree of life, sometimes leading to rapid increase of phenotypic variability, to introgression of adaptive alleles, to the formation of hybrid species, and even to entire species radiations. In this review, we identify consensus among diverse research programs to show how the field has progressed. Hybridization is a multifaceted evolutionary process that can strongly influence species formation and facilitate adaptation and persistence of species in a rapidly changing world. Progress on testing this hypothesis will require cooperation among different subdisciplines.

Tangled banks, braided rivers, and complex hierarchies: beyond microevolution and macroevolution

Ever since the Modern Synthesis, a debate about the relationship between microevolution and macroevolution has persisted-specifically, whether they are equivalent, distinct, or explain one another. How one answers these questions has become shorthand for a much broader set of theoretical debates in evolutionary biology. Here, we examine microevolution and macroevolution in the context of the vast proliferation of data, knowledge, and theory since the advent of the Modern Synthesis. We suggest that traditional views on microevolution and macroevolution are too binary and reductive given current empirical and theoretical advances in biology. For example, patterns and processes are interconnected at various temporal and spatial scales and among hierarchical entities, rather than defining micro- or macro-domains. Further, biological entities have variably fuzzy boundaries, resulting in complex evolutionary processes that influence macroevolution occuring at both micro- and macro-levels. In addition, conceptual advances in phylodynamics have yet to be fully integrated with contemporary macroevolutionary approaches. Finally, holding microevolution and macroevolution as distinct domains thwarts synthesis and collaboration on important research questions. Instead, we propose that the focal entities and processes considered by evolutionary studies be contextualized within the complexity of the multidimensional, multimodal, multilevel phylogenetic system.

Jumping through hoops: Structural rearrangements and accelerated mutation rates on Dendrodorididae (Mollusca: Nudibranchia) mitogenomes rumble their evolution

The systematics of the family Dendrodorididae, with only three valid genera, is a challenge for integrative taxonomists. Its members lack hard structures for morphological comparisons and their mitochondrial and nuclear markers provide contradictory phylogenetic signals, making phylogenetic reconstructions difficult. This molecular discordance has been hypothesized to be the result of nuclear pseudogenes or exogenous contamination. However, these hypotheses have not been tested. Here, we assembled the first genome drafts of seven Dendrodorididae species to investigate the evolutionary processes of this family. Two of the mitogenomes displayed an identical structural rearrangement involving the translocation of three coding genes and five tRNAs, described for the first time in nudibranchs. In addition, we found particularly high dN and dN/dS values and multiple insertions and deletions on the mitochondrial genes of smooth Dendrodoris. In contrast, nuclear single-copy ortholog genes showed no such mutational differences. Models of protein structures from mitochondrial genes are conserved, suggesting conserved functionality. Phylogenies using mitogenomic and nuclear data showed that species with rearranged mitogenomes form a clade, although Dendrodorididae relationships remained unresolved. The present study provides novel evidence for accelerated mutation rates in the mitogenomes of Dendrodorididae, which presumably have implications on respiratory adaptation, and highlights the importance of using genomic data to unveil rare evolutionary processes, crucial for correctly inferring phylogenies.

Ice age-driven range shifts of diploids and expanding autotetraploids of Biscutella laevigata within a conserved niche

Early studies of the textbook mixed-ploidy system Biscutella laevigata highlighted diploids restricted to never-glaciated lowlands and tetraploids at high elevations across the European Alps, promoting the hypothesis that whole-genome duplication (WGD) is advantageous under environmental changes. Here we addressed long-held hypotheses on the role of hybridisation at the origin of the tetraploids, their single vs multiple origins, and whether a shift in climatic niche accompanied WGD. Climatic niche modelling together with spatial genetics and coalescent modelling based on ddRAD-seq genotyping of 17 diploid and 19 tetraploid populations was used to revisit the evolution of this species complex in space and time. Diploids differentiated into four genetic lineages corresponding to allopatric glacial refugia at the onset of the last ice age, whereas tetraploids displaying tetrasomic inheritance formed a uniform group that originated from southern diploids before the last glacial maximum. Derived from diploids occurring at high elevation, autotetraploids likely inherited their adaptation to high elevation rather than having evolved it through or after WGD. They further presented considerable postglacial expansion across the Alps and underwent admixture with diploids. Although the underpinnings of the successful expansion of autotetraploids remain elusive, differentiation in B. laevigata was chiefly driven by the glacial history of the Alps.

Assessment of Population Genetic Diversity of Medicinal Meconopsis integrifolia (Maxim.) Franch. Using Newly Developed SSR Markers

Meconopsis integrifolia is an endangered Tibetan medicinal plant with significant medicinal and ornamental value. Understanding its genetic diversity and structure is crucial for its sustainable utilization and effective conservation. Here, we develop a set of SSR markers based on transcriptome data to analyze the genetic diversity and structure of 185 individuals from 16 populations of M. integrifolia. The results indicate that M. integrifolia exhibits relatively high genetic diversity at the species level (the percentage of polymorphic bands PPB = 91.67%, Nei's genetic diversity index He = 0.2989, Shannon's information index I = 0.4514) but limited genetic variation within populations (PPB = 12.08%, He = 0.0399, I = 0.0610). The genetic differentiation among populations is relatively high (the coefficient of gene differentiation GST = 0.6902), and AMOVA analysis indicates that 63.39% of the total variation occurs among populations. This suggests that maintaining a limited number of populations is insufficient to preserve the overall diversity of M. integrifolia. Different populations are categorized into four representative subclusters, but they do not cluster strictly according to geographical distribution. Limited gene flow (Nm = 0.2244) is likely the main reason for the high differentiation among these populations. Limited seed and pollen dispersal abilities, along with habitat fragmentation, may explain the restricted gene flow among populations, highlighting the necessity of conserving as many populations in the wild as possible.

The Ecology of Hybrid Incompatibilities

Ecologically mediated selection against hybrids, caused by hybrid phenotypes fitting poorly into available niches, is typically viewed as distinct from selection caused by epistatic Dobzhansky-Muller hybrid incompatibilities. Here, we show how selection against transgressive phenotypes in hybrids manifests as incompatibility. After outlining our logic, we summarize current approaches for studying ecology-based selection on hybrids. We then quantitatively review QTL-mapping studies and find traits differing between parent taxa are typically polygenic. Next, we describe how verbal models of selection on hybrids translate to phenotypic and genetic fitness landscapes, highlighting emerging approaches for detecting polygenic incompatibilities. Finally, in a synthesis of published data, we report that trait transgression-and thus possibly extrinsic hybrid incompatibility in hybrids-escalates with the phenotypic divergence between parents. We discuss conceptual implications and conclude that studying the ecological basis of hybrid incompatibility will facilitate new discoveries about mechanisms of speciation.

Chromosomal inversions harbour excess mutational load in the coral, Acropora kenti, on the Great Barrier Reef

The future survival of coral reefs in the Anthropocene depends on the capacity of corals to adapt as oceans warm and extreme weather events become more frequent. Targeted interventions designed to assist evolutionary processes in corals require a comprehensive understanding of the distribution and structure of standing variation, however, efforts to map genomic variation in corals have so far focussed almost exclusively on SNPs, overlooking structural variants that have been shown to drive adaptive processes in other taxa. Here, we show that the reef-building coral, Acropora kenti, harbours at least five large, highly polymorphic structural variants, all of which exhibit signatures of strongly suppressed recombination in heterokaryotypes, a feature commonly associated with chromosomal inversions. Based on their high minor allele frequency, uniform distribution across habitats and elevated genetic load, we propose that these inversions in A. kenti are likely to be under balancing selection. An excess of SNPs with high impact on protein-coding genes within these loci elevates their importance both as potential targets for adaptive selection and as contributors to genetic decline if coral populations become fragmented or inbred in future.

RIDGE, a tool tailored to detect gene flow barriers across species pairs

Characterizing the processes underlying reproductive isolation between diverging lineages is central to understanding speciation. Here, we present RIDGE-Reproductive Isolation Detection using Genomic polymorphisms-a tool tailored for quantifying gene flow barrier proportion and identifying the relevant genomic regions. RIDGE relies on an Approximate Bayesian Computation with a model-averaging approach to accommodate diverse scenarios of lineage divergence. It captures heterogeneity in effective migration rate along the genome while accounting for variation in linked selection and recombination. The barrier detection test relies on numerous summary statistics to compute a Bayes factor, offering a robust statistical framework that facilitates cross-species comparisons. Simulations revealed RIDGE's efficiency in capturing signals of ongoing migration. Model averaging proved particularly valuable in scenarios of high model uncertainty where no migration or migration homogeneity can be wrongly assumed, typically for recent divergence times <0.1 2Ne generations. Applying RIDGE to four published crow data sets, we first validated our tool by identifying a well-known large genomic region associated with mate choice patterns. Second, while we identified a significant overlap of outlier loci using RIDGE and traditional genomic scans, our results suggest that a substantial portion of previously identified outliers are likely false positives. Outlier detection relies on allele differentiation, relative measures of divergence and the count of shared polymorphisms and fixed differences. Our analyses also highlight the value of incorporating multiple summary statistics including our newly developed outlier ones that can be useful in challenging detection conditions.

Structural Variants and Speciation: Multiple Processes at Play

Research on the genomic architecture of speciation has increasingly revealed the importance of structural variants (SVs) that affect the presence, abundance, position, and/or direction of a nucleotide sequence. SVs include large chromosomal rearrangements such as fusion/fissions and inversions and translocations, as well as smaller variants such as duplications, insertions, and deletions (CNVs). Although we have ample evidence that SVs play a key role in speciation, the underlying mechanisms differ depending on the type and length of the SV, as well as the ecological, demographic, and historical context. We review predictions and empirical evidence for classic processes such as underdominance due to meiotic aberrations and the coupling effect of recombination suppression before exploring how recent sequencing methodologies illuminate the prevalence and diversity of SVs. We discuss specific properties of SVs and their impact throughout the genome, highlighting that multiple processes are at play, and possibly interacting, in the relationship between SVs and speciation.

Genetic architecture of ecological divergence between Oryza rufipogon and Oryza nivara

Ecological divergence due to habitat difference plays a prominent role in the formation of new species, but the genetic architecture during ecological speciation and the mechanism underlying phenotypic divergence remain less understood. Two wild ancestors of rice (Oryza rufipogon and Oryza nivara) are a progenitor-derivative species pair with ecological divergence and provide a unique system for studying ecological adaptation/speciation. Here, we constructed a high-resolution linkage map and conducted a quantitative trait locus (QTL) analysis of 19 phenotypic traits using an F2 population generated from a cross between the two Oryza species. We identified 113 QTLs associated with interspecific divergence of 16 quantitative traits, with effect sizes ranging from 1.61% to 34.1% in terms of the percentage of variation explained (PVE). The distribution of effect sizes of QTLs followed a negative exponential, suggesting that a few genes of large effect and many genes of small effect were responsible for the phenotypic divergence. We observed 18 clusters of QTLs (QTL hotspots) on 11 chromosomes, significantly more than that expected by chance, demonstrating the importance of coinheritance of loci/genes in ecological adaptation/speciation. Analysis of effect direction and v-test statistics revealed that interspecific differentiation of most traits was driven by divergent natural selection, supporting the argument that ecological adaptation/speciation would proceed rapidly under coordinated selection on multiple traits. Our findings provide new insights into the understanding of genetic architecture of ecological adaptation and speciation in plants and help effective manipulation of specific genes or gene cluster in rice breeding.

Divergent Selection in Low Recombination Regions Shapes the Genomic Islands in Two Incipient Shorebird Species

Speciation in the face of gene flow is usually associated with a heterogeneous genomic landscape of divergence in nascent species pairs. However, multiple factors, such as divergent selection and local recombination rate variation, can influence the formation of these genomic islands. Examination of the genomic landscapes of species pairs that are still in the early stages of speciation provides an insight into this conundrum. In this study, population genomic analyses were undertaken using a wide range of sampling and whole-genome resequencing data from 96 unrelated individuals of Kentish plover (Charadrius alexandrinus) and white-faced plover (Charadrius dealbatus). We suggest that the two species exhibit varying levels of population admixture along the Chinese coast and on the Taiwan Island. Genome-wide analyses for introgression indicate that ancient introgression had occurred in Taiwan population, and gene flow is still ongoing in mainland coastal populations. Furthermore, we identified a few genomic regions with significant levels of interspecific differentiation and local recombination suppression, which contain several genes potentially associated with disease resistance, coloration, and regulation of plumage molting and thus may be relevant to the phenotypic and ecological divergence of the two nascent species. Overall, our findings suggest that divergent selection in low recombination regions may be a main force in shaping the genomic islands in two incipient shorebird species.

The era of the ARG: An introduction to ancestral recombination graphs and their significance in empirical evolutionary genomics

In the presence of recombination, the evolutionary relationships between a set of sampled genomes cannot be described by a single genealogical tree. Instead, the genomes are related by a complex, interwoven collection of genealogies formalized in a structure called an ancestral recombination graph (ARG). An ARG extensively encodes the ancestry of the genome(s) and thus is replete with valuable information for addressing diverse questions in evolutionary biology. Despite its potential utility, technological and methodological limitations, along with a lack of approachable literature, have severely restricted awareness and application of ARGs in evolution research. Excitingly, recent progress in ARG reconstruction and simulation have made ARG-based approaches feasible for many questions and systems. In this review, we provide an accessible introduction and exploration of ARGs, survey recent methodological breakthroughs, and describe the potential for ARGs to further existing goals and open avenues of inquiry that were previously inaccessible in evolutionary genomics. Through this discussion, we aim to more widely disseminate the promise of ARGs in evolutionary genomics and encourage the broader development and adoption of ARG-based inference.

A novel tetra-primer ARMS-PCR approach for the molecular karyotyping of chromosomal inversion 2Ru in the main malaria vectors Anopheles gambiae and Anopheles coluzzii

BACKGROUND

Chromosomal inversion polymorphisms have been associated with adaptive behavioral, physiological, morphological and life history traits in the two main Afrotropical malaria vectors, Anopheles coluzzii and Anopheles gambiae. The understanding of the adaptive value of chromosomal inversion systems is constrained by the feasibility of cytological karyotyping. In recent years in silico and molecular approaches have been developed for the genotyping of most widespread inversions (2La, 2Rb and 2Rc). The 2Ru inversion, spanning roughly 8% of chromosome 2R, is commonly polymorphic in West African populations of An. coluzzii and An. gambiae and shows clear increases in frequency with increasing rainfall seasonally and geographically. The aim of this work was to overcome the constraints of currently available cytological and high-throughput molecular assays by developing a simple PCR assay for genotyping the 2Ru inversion in individual specimens of both mosquito species.

METHODS

We designed tetra-primer amplification refractory mutation system (ARMS)-PCR assays based on five tag single-nucleotide polymorphisms (SNPs) previously shown to be strongly correlated with 2Ru inversion orientation. The most promising assay was validated against laboratory and field samples of An. coluzzii and An. gambiae karyotyped either cytogenetically or molecularly using a genotyping-in-thousands by sequencing (GT-seq) high-throughput approach that employs targeted sequencing of multiplexed PCR amplicons.

RESULTS

A successful assay was designed based on the tag SNP at position 2R, 31710303, which is highly predictive of the 2Ru genotype. The assay, which requires only one PCR, and no additional post-PCR processing other than electrophoresis, produced a clear banding pattern for 98.5% of the 454 specimens tested, which is a 96.7% agreement with established karyotyping methods. Sequences were obtained for nine of the An. coluzzii specimens manifesting 2Ru genotype discrepancies with GT-seq. Possible sources of these discordances are discussed.

CONCLUSIONS

The tetra-primer ARMS-PCR assay represents an accurate, streamlined and cost-effective method for the molecular karyotyping of the 2Ru inversion in An. coluzzii and An. gambiae. Together with approaches already available for the other common polymorphic inversions, 2La, 2Rb and 2Rc, this assay will allow investigations of the adaptive value of the complex set of inversion systems observed in the two major malaria vectors in the Afrotropical region.

The era of the ARG: an empiricist's guide to ancestral recombination graphs

In the presence of recombination, the evolutionary relationships between a set of sampled genomes cannot be described by a single genealogical tree. Instead, the genomes are related by a complex, interwoven collection of genealogies formalized in a structure called an ancestral recombination graph (ARG). An ARG extensively encodes the ancestry of the genome(s) and thus is replete with valuable information for addressing diverse questions in evolutionary biology. Despite its potential utility, technological and methodological limitations, along with a lack of approachable literature, have severely restricted awareness and application of ARGs in empirical evolution research. Excitingly, recent progress in ARG reconstruction and simulation have made ARG-based approaches feasible for many questions and systems. In this review, we provide an accessible introduction and exploration of ARGs, survey recent methodological breakthroughs, and describe the potential for ARGs to further existing goals and open avenues of inquiry that were previously inaccessible in evolutionary genomics. Through this discussion, we aim to more widely disseminate the promise of ARGs in evolutionary genomics and encourage the broader development and adoption of ARG-based inference.

Mechanisms of Intrinsic Postzygotic Isolation: From Traditional Genic and Chromosomal Views to Genomic and Epigenetic Perspectives

Intrinsic postzygotic isolation typically appears as reduced viability or fertility of interspecific hybrids caused by genetic incompatibilities between diverged parental genomes. Dobzhansky-Muller interactions among individual genes, and chromosomal rearrangements causing problems with chromosome synapsis and recombination in meiosis, have both long been considered as major mechanisms behind intrinsic postzygotic isolation. Recent research has, however, suggested that the genetic basis of intrinsic postzygotic isolation can be more complex and involves, for example, overall divergence of the DNA sequence or epigenetic changes. Here, we review the mechanisms of intrinsic postzygotic isolation from genic, chromosomal, genomic, and epigenetic perspectives across diverse taxa. We provide empirical evidence for these mechanisms, discuss their importance in the speciation process, and highlight questions that remain unanswered.

Genomics of plant speciation

Studies of plants have been instrumental for revealing how new species originate. For several decades, botanical research has complemented and, in some cases, challenged concepts on speciation developed via the study of other organisms while also revealing additional ways in which species can form. Now, the ability to sequence genomes at an unprecedented pace and scale has allowed biologists to settle decades-long debates and tackle other emerging challenges in speciation research. Here, we review these recent genome-enabled developments in plant speciation. We discuss complications related to identification of reproductive isolation (RI) loci using analyses of the landscape of genomic divergence and highlight the important role that structural variants have in speciation, as increasingly revealed by new sequencing technologies. Further, we review how genomics has advanced what we know of some routes to new species formation, like hybridization or whole-genome duplication, while casting doubt on others, like population bottlenecks and genetic drift. While genomics can fast-track identification of genes and mutations that confer RI, we emphasize that follow-up molecular and field experiments remain critical. Nonetheless, genomics has clarified the outsized role of ancient variants rather than new mutations, particularly early during speciation. We conclude by highlighting promising avenues of future study. These include expanding what we know so far about the role of epigenetic and structural changes during speciation, broadening the scope and taxonomic breadth of plant speciation genomics studies, and synthesizing information from extensive genomic data that have already been generated by the plant speciation community.

Ancient standing genetic variation facilitated the adaptive radiation of Lake Victoria cichlids

Cichlid fishes are textbook examples of explosive speciation and adaptive radiation, providing a great opportunity to understand how the genomic substrate yields extraordinary species diversity. Recently, we performed comparative genomic analyses of three Lake Victoria cichlids to reveal the genomic substrates underlying their rapid speciation and adaptation. We found that long divergent haplotypes derived from large-scale standing genetic variation, which originated before the adaptive radiation of Lake Victoria cichlids, may have contributed to their rapid diversification. In addition, the present study on genomic data from other East African cichlids suggested the reuse of alleles that may have originated in the ancestral lineages of Lake Tanganyika cichlids during cichlid evolution. Therefore, our results highlight that the primary factor that could drive repeated adaptive radiation across East African cichlids was allelic reuse from standing genetic variation to adapt to their own specific environment. In this report, we summarize the main results and discuss the evolutionary mechanisms of cichlids, based on our latest findings.

Ecology, the pace-of-life, epistatic selection and the maintenance of genetic variation in life-history genes

Evolutionary genetics has long struggled with understanding how functional genes under selection remain polymorphic in natural populations. Taking as a starting point that natural selection is ultimately a manifestation of ecological processes, we spotlight an underemphasized and potentially ubiquitous ecological effect that may have fundamental effects on the maintenance of genetic variation. Negative frequency dependency is a well-established emergent property of density dependence in ecology, because the relative profitability of different modes of exploiting or utilizing limiting resources tends to be inversely proportional to their frequency in a population. We suggest that this may often generate negative frequency-dependent selection (NFDS) on major effect loci that affect rate-dependent physiological processes, such as metabolic rate, that are phenotypically manifested as polymorphism in pace-of-life syndromes. When such a locus under NFDS shows stable intermediate frequency polymorphism, this should generate epistatic selection potentially involving large numbers of loci with more minor effects on life-history (LH) traits. When alternative alleles at such loci show sign epistasis with a major effect locus, this associative NFDS will promote the maintenance of polygenic variation in LH genes. We provide examples of the kind of major effect loci that could be involved and suggest empirical avenues that may better inform us on the importance and reach of this process.

Extensive introgression among strongylocentrotid sea urchins revealed by phylogenomics

Gametic isolation is thought to play an important role in the evolution of reproductive isolation in broadcast-spawning marine invertebrates. However, it is unclear whether gametic isolation commonly evolves early in the speciation process or only accumulates after other reproductive barriers are already in place. It is also unknown whether gametic isolation is an effective barrier to introgression following speciation. Here, we used whole-genome sequencing data and multiple complementary phylogenomic approaches to test whether the well-documented gametic incompatibilities among the strongylocentrotid sea urchins have limited introgression. We quantified phylogenetic discordance, inferred reticulate phylogenetic networks, and applied the Δ statistic using gene tree topologies reconstructed from multiple sequence alignments of protein-coding single-copy orthologs. In addition, we conducted ABBA-BABA tests on genome-wide single nucleotide variants and reconstructed a phylogeny of mitochondrial genomes. Our results revealed strong mito-nuclear discordance and considerable nonrandom gene tree discordance that cannot be explained by incomplete lineage sorting alone. Eight of the nine species examined demonstrated a history of introgression with at least one other species or ancestral lineage, indicating that introgression was common during the diversification of the strongylocentrotid urchins. There was strong support for introgression between four extant species pairs (Strongylocentrotus pallidus ⇔ S. droebachiensis, S. intermedius ⇔ S. pallidus, S. purpuratus ⇔ S. fragilis, and Mesocentrotus franciscanus ⇔ Pseudocentrotus depressus) and additional evidence for introgression on internal branches of the phylogeny. Our results suggest that the existing gametic incompatibilities among the strongylocentrotid urchin species have not been a complete barrier to hybridization and introgression following speciation. Their continued divergence in the face of widespread introgression indicates that other reproductive isolating barriers likely exist and may have been more critical in establishing reproductive isolation early in speciation.

Conceptual and empirical bridges between micro- and macroevolution

Explaining broad molecular, phenotypic and species biodiversity patterns necessitates a unifying framework spanning multiple evolutionary scales. Here we argue that although substantial effort has been made to reconcile microevolution and macroevolution, much work remains to identify the links between biological processes at play. We highlight four major questions of evolutionary biology whose solutions require conceptual bridges between micro and macroevolution. We review potential avenues for future research to establish how mechanisms at one scale (drift, mutation, migration, selection) translate to processes at the other scale (speciation, extinction, biogeographic dispersal) and vice versa. We propose ways in which current comparative methods to infer molecular evolution, phenotypic evolution and species diversification could be improved to specifically address these questions. We conclude that researchers are in a better position than ever before to build a synthesis to understand how microevolutionary dynamics unfold over millions of years.

Genetic architecture of a pollinator shift and its fate in secondary hybrid zones of two Petunia species

BACKGROUND

Theory suggests that the genetic architecture of traits under divergent natural selection influences how easily reproductive barriers evolve and are maintained between species. Divergently selected traits with a simple genetic architecture (few loci with major phenotypic effects) should facilitate the establishment and maintenance of reproductive isolation between species that are still connected by some gene flow. While empirical support for this idea appears to be mixed, most studies test the influence of trait architectures on reproductive isolation only indirectly. Petunia plant species are, in part, reproductively isolated by their different pollinators. To investigate the genetic causes and consequences of this ecological isolation, we deciphered the genetic architecture of three floral pollination syndrome traits in naturally occurring hybrids between the widespread Petunia axillaris and the highly endemic and endangered P. exserta.

RESULTS

Using population genetics, Bayesian linear mixed modelling and genome-wide association studies, we found that the three pollination syndrome traits vary in genetic architecture. Few genome regions explain a majority of the variation in flavonol content (defining UV floral colour) and strongly predict the trait value in hybrids irrespective of interspecific admixture in the rest of their genomes. In contrast, variation in pistil exsertion and anthocyanin content (defining visible floral colour) is controlled by many genome-wide loci. Opposite to flavonol content, the genome-wide proportion of admixture between the two species predicts trait values in their hybrids. Finally, the genome regions strongly associated with the traits do not show extreme divergence between individuals representing the two species, suggesting that divergent selection on these genome regions is relatively weak within their contact zones.

CONCLUSIONS

Among the traits analysed, those with a more complex genetic architecture are best maintained in association with the species upon their secondary contact. We propose that this maintained genotype-phenotype association is a coincidental consequence of the complex genetic architectures of these traits: some of their many underlying small-effect loci are likely to be coincidentally linked with the actual barrier loci keeping these species partially isolated upon secondary contact. Hence, the genetic architecture of a trait seems to matter for the outcome of hybridization not only then when the trait itself is under selection.

Standing genetic variation and chromosome differences drove rapid ecotype formation in a major malaria mosquito

Species distributed across heterogeneous environments often evolve locally adapted ecotypes, but understanding of the genetic mechanisms involved in their formation and maintenance in the face of gene flow is incomplete. In Burkina Faso, the major African malaria mosquito Anopheles funestus comprises two strictly sympatric and morphologically indistinguishable yet karyotypically differentiated forms reported to differ in ecology and behavior. However, knowledge of the genetic basis and environmental determinants of An. funestus diversification was impeded by lack of modern genomic resources. Here, we applied deep whole-genome sequencing and analysis to test the hypothesis that these two forms are ecotypes differentially adapted to breeding in natural swamps versus irrigated rice fields. We demonstrate genome-wide differentiation despite extensive microsympatry, synchronicity, and ongoing hybridization. Demographic inference supports a split only ~1,300 y ago, closely following the massive expansion of domesticated African rice cultivation ~1,850 y ago. Regions of highest divergence, concentrated in chromosomal inversions, were under selection during lineage splitting, consistent with local adaptation. The origin of nearly all variations implicated in adaptation, including chromosomal inversions, substantially predates the ecotype split, suggesting that rapid adaptation was fueled mainly by standing genetic variation. Sharp inversion frequency differences likely facilitated adaptive divergence between ecotypes by suppressing recombination between opposing chromosomal orientations of the two ecotypes, while permitting free recombination within the structurally monomorphic rice ecotype. Our results align with growing evidence from diverse taxa that rapid ecological diversification can arise from evolutionarily old structural genetic variants that modify genetic recombination.

Incompatibility and Interchangeability in Molecular Evolution

There is remarkable variation in the rate at which genetic incompatibilities in molecular interactions accumulate. In some cases, minor changes-even single-nucleotide substitutions-create major incompatibilities when hybridization forces new variants to function in a novel genetic background from an isolated population. In other cases, genes or even entire functional pathways can be horizontally transferred between anciently divergent evolutionary lineages that span the tree of life with little evidence of incompatibilities. In this review, we explore whether there are general principles that can explain why certain genes are prone to incompatibilities while others maintain interchangeability. We summarize evidence pointing to four genetic features that may contribute to greater resistance to functional replacement: (1) function in multisubunit enzyme complexes and protein-protein interactions, (2) sensitivity to changes in gene dosage, (3) rapid rate of sequence evolution, and (4) overall importance to cell viability, which creates sensitivity to small perturbations in molecular function. We discuss the relative levels of support for these different hypotheses and lay out future directions that may help explain the striking contrasts in patterns of incompatibility and interchangeability throughout the history of molecular evolution.

Sympatric or micro-allopatric speciation in a glacial lake? Genomic islands support neither

Apparent cases of sympatric speciation may actually be due to micro-allopatric or micro-parapatric speciation. One way to distinguish between these models is to examine the existence and nature of genomic islands of divergence, wherein divergent DNA segments are interspersed with low-divergence segments. Such islands should be rare or absent under micro-allopatric speciation but common in cases of speciation with gene flow. Sympatric divergence of endemic fishes is known from isolated saline, crater, postglacial, and ancient lakes. Two morphologically distinct cyprinid fishes, Gymnocypris eckloni scoliostomus (GS) and G. eckloni eckloni (GE), in a small glacial lake on the Qinghai-Tibet Plateau, Lake Sunmcuo, match the biogeographic criteria of sympatric speciation. In this study, we examined genome-wide variation in 46 individuals from these two groups. The divergence time between the GS and GE lineages was estimated to be 20-60 Kya. We identified 54 large genomic islands (≥100 kb) of speciation, which accounted for 89.4% of the total length of all genomic islands. These islands harboured divergent genes related to olfactory receptors and olfaction signals that may play important roles in food selection and assortative mating in fishes. Although the genomic islands clearly indicated speciation with gene flow and rejected micro-allopatric speciation, they were too large to support the hypothesis of sympatric speciation. Theoretical and recent empirical studies suggested that continual gene flow in sympatry should give rise to many small genomic islands (as small as a few kilobases in size). Thus, the observed pattern is consistent with the extensive evidence on parapatric speciation, in which adjacent habitats facilitate divergent selection but also permit gene flow during speciation. We suggest that many, if not most, of the reported cases of sympatric speciation are likely to be micro-parapatric speciation.