The genes underlying the process of speciation

Thrushes in Love: Extensive Gene Flow, With Differential Resistance and Selection, Obscures and Reveals the Evolutionary History of a Songbird Clade

The application of high-throughput sequencing to phylogenetic analyses is allowing authors to reconstruct the true evolutionary history of species. This work can illuminate specific mechanisms underlying divergence when combined with analyses of gene flow, recombination and selection. We conducted a phylogenomic analysis of Catharus, a songbird genus with considerable potential for gene flow, variation in migratory behaviour and genomic resources. We documented discordance among trees constructed for mitochondrial, autosomal and sex (Z) chromosome partitions. Two trees were recovered on the Z. Both trees differed from the autosomes, one matched the mitochondria, and the other was unique to the Z. Gene flow with one species likely generated much of this discordance; substantial admixture between ustulatus and the remaining species was documented and linked to at least two historic events. The tree unique to the Z likely reflects the true history of Catharus; local genomic analyses recovered the same tree in autosomal regions with reduced admixture and recombination. Genes previously connected to migration were enriched in these regions suggesting transitions between migratory and non-migratory states helped generate divergence. Migratory (vs. nonmigratory) Catharus formed a monophyletic clade in a subset of genomic regions. Gene flow was elevated in some of these regions suggesting adaptive introgression may have occurred, but the dominant pattern was of balancing selection maintaining ancestral polymorphisms important for olfaction and perhaps, by extension, adaptation to temperate climates. This work illuminates the evolutionary history of an important model in speciation and demonstrates how differential resistance to gene flow can affect local genomic patterns.

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.

Meiotic Drive and Speciation

Meiotic drive is the biased transmission of alleles from heterozygotes, contrary to Mendel's laws, and reflects intragenomic conflict rather than organism-level Darwinian selection. Theory has been developed as to how centromeric properties can promote female meiotic drive and how conflict between the X and Y chromosomes in males can promote male meiotic drive. There are empirical data that fit both the centromere drive and sex chromosome drive models. Sex chromosome drive may have relevance to speciation through the buildup of Dobzhansky-Muller incompatibilities involving drive and suppressor systems, studied particularly in Drosophila. Centromere drive may promote fixation of chromosomal rearrangements involving the centromere, and those fixed rearrangements may contribute to reproductive isolation, studied particularly in the house mouse. Genome-wide tests suggest that meiotic drive promotes allele fixation with regularity, and those studying the genomics of speciation need to be aware of the potential impact of such fixations on reproductive isolation. New species can originate in many different ways (including multiple factors acting together), and a substantial body of work on meiotic drive point to it being one of the processes involved.

The biogeographic and evolutionary processes shaping population divergence in Laupala

Speciation generates biodiversity and the mechanisms involved are thought to vary across the tree of life and across environments. For example, well-studied adaptive radiations are thought to be fuelled by divergent ecological selection, but additionally are influenced heavily by biogeographic, genomic and demographic factors. Mechanisms of non-adaptive radiations, producing ecologically cryptic taxa, have been less well-studied but should likewise be influenced by these latter factors. Comparing among contexts can help pinpoint universal mechanisms and outcomes, especially if we integrate biogeographic, ecological and evolutionary processes. We investigate population divergence in the swordtail cricket Laupala cerasina, a wide-spread endemic on Hawai'i Island and one of 38 ecologically cryptic Laupala species. The nine sampled populations show striking population genetic structure at small spatio-temporal scales. The rapid differentiation among populations and species of Laupala shows that neither a specific geographical context nor ecological opportunity are pre-requisites for rapid divergence. Spatio-temporal patterns in population divergence, population size change, and gene flow are aligned with the chronosequence of the four volcanoes on which L. cerasina occurs and reveal the composite effects of geological dynamics and Quaternary climate change on population dynamics. Spatio-temporal patterns in genetic variation along the genome reveal the interplay of genetic and genomic architecture in shaping population divergence. In early phases of divergence, we find elevated differentiation in genomic regions harbouring mating song loci. In later stages of divergence, we find a signature of linked selection that interacts with recombination rate variation. Comparing our findings with recent work on complementary systems supports the conclusion that mostly universal factors influence the speciation process.

Genome sequences and population genomics reveal climatic adaptation and genomic divergence between two closely related sweetgum species

Understanding the genetic basis of population divergence and adaptation is an important goal in population genetics and evolutionary biology. However, the relative roles of demographic history, gene flow, and/or selective regime in driving genomic divergence, climatic adaptation, and speciation in non-model tree species are not yet fully understood. To address this issue, we generated whole-genome resequencing data of Liquidambar formosana and L. acalycina, which are broadly sympatric but altitudinally segregated in the Tertiary relict forests of subtropical China. We integrated genomic and environmental data to investigate the demographic history, genomic divergence, and climatic adaptation of these two sister species. We inferred a scenario of allopatric species divergence during the late Miocene, followed by secondary contact during the Holocene. We identified multiple genomic islands of elevated divergence that mainly evolved through divergence hitchhiking and recombination rate variation, likely fostered by long-term refugial isolation and recent differential introgression in low-recombination genomic regions. We also found some candidate genes with divergent selection signatures potentially involved in climatic adaptation and reproductive isolation. Our results contribute to a better understanding of how late Tertiary/Quaternary climatic change influenced speciation, genomic divergence, climatic adaptation, and introgressive hybridization in East Asia's Tertiary relict flora. In addition, they should facilitate future evolutionary, conservation genomics, and molecular breeding studies in Liquidambar, a genus of important medicinal and ornamental values.

Yosemite toad (Anaxyrus canorus) transcriptome reveals interplay between speciation genes and adaptive introgression

Genomes are heterogeneous during the early stages of speciation, with small 'islands' of DNA appearing to reflect strong adaptive differences, surrounded by vast seas of relative homogeneity. As species diverge, secondary contact zones between them can act as an interface and selectively filter through advantageous alleles of hybrid origin. Such introgression is another important adaptive process, one that allows beneficial mosaics of recombinant DNA ('rivers') to flow from one species into another. Although genomic islands of divergence appear to be associated with reproductive isolation, and genomic rivers form by adaptive introgression, it is unknown whether islands and rivers tend to be the same or different loci. We examined three replicate secondary contact zones for the Yosemite toad (Anaxyrus canorus) using two genomic data sets and a morphometric data set to answer the questions: (1) How predictably different are islands and rivers, both in terms of genomic location and gene function? (2) Are the adaptive genetic trait loci underlying tadpole growth and development reliably islands, rivers or neither? We found that island and river loci have significant overlap within a contact zone, suggesting that some loci are first islands, and later are predictably converted into rivers. However, gene ontology enrichment analysis showed strong overlap in gene function unique to all island loci, suggesting predictability in overall gene pathways for islands. Genome-wide association study outliers for tadpole development included LPIN3, a lipid metabolism gene potentially involved in climate change adaptation, that is island-like for all three contact zones, but also appears to be introgressing (as a river) across one zone. Taken together, our results suggest that adaptive divergence and introgression may be more complementary forces than currently appreciated.

Metazoa-level USCOs as markers in species delimitation and classification

Metazoa-level universal single-copy orthologs (mzl-USCOs) are universally applicable markers for DNA taxonomy in animals that can replace or supplement single-gene barcodes. Previously, mzl-USCOs from target enrichment data were shown to reliably distinguish species. Here, we tested whether USCOs are an evenly distributed, representative sample of a given metazoan genome and therefore able to cope with past hybridization events and incomplete lineage sorting. This is relevant for coalescent-based species delimitation approaches, which critically depend on the assumption that the investigated loci do not exhibit autocorrelation due to physical linkage. Based on 239 chromosome-level assembled genomes, we confirmed that mzl-USCOs are genetically unlinked for practical purposes and a representative sample of a genome in terms of reciprocal distances between USCOs on a chromosome and of distribution across chromosomes. We tested the suitability of mzl-USCOs extracted from genomes for species delimitation and phylogeny in four case studies: Anopheles mosquitos, Drosophila fruit flies, Heliconius butterflies and Darwin's finches. In almost all instances, USCOs allowed delineating species and yielded phylogenies that corresponded to those generated from whole genome data. Our phylogenetic analyses demonstrate that USCOs may complement single-gene DNA barcodes and provide more accurate taxonomic inferences. Combining USCOs from sources that used different versions of ortholog reference libraries to infer marker orthology may be challenging and, at times, impact taxonomic conclusions. However, we expect this problem to become less severe as the rapidly growing number of reference genomes provides a better representation of the number and diversity of organismal lineages.

Gene regulation and speciation in a migratory divide between songbirds

Behavioral variation abounds in nature. This variation is important for adaptation and speciation, but its molecular basis remains elusive. Here, we use a hybrid zone between two subspecies of songbirds that differ in migration - an ecologically important and taxonomically widespread behavior---to gain insight into this topic. We measure gene expression in five brain regions. Differential expression between migratory states was dominated by circadian genes in all brain regions. The remaining patterns were largely brain-region specific. For example, expression differences between the subspecies that interact with migratory state likely help maintain reproductive isolation in this system and were documented in only three brain regions. Contrary to existing work on regulatory mechanisms underlying species-specific traits, two lines of evidence suggest that trans- (vs. cis) regulatory changes underlie these patterns - no evidence for allele-specific expression in hybrids and minimal associations between genomic differentiation and expression differences. Additional work with hybrids shows expression levels were often distinct (transgressive) from parental forms. Behavioral contrasts and functional enrichment analyses allowed us to connect these patterns to mitonuclear incompatibilities and compensatory responses to stress that could exacerbate selection on hybrids and contribute to speciation.

Quantitative Analyses of Coupling in Hybrid Zones

In hybrid zones, whether barrier loci experience selection mostly independently or as a unit depends on the ratio of selection to recombination as captured by the coupling coefficient. Theory predicts a sharper transition between an uncoupled and coupled system when more loci affect hybrid fitness. However, the extent of coupling in hybrid zones has rarely been quantified. Here, we use simulations to characterize the relationship between the coupling coefficient and variance in clines across genetic loci. We then reanalyze 25 hybrid zone data sets and find that cline variances and estimated coupling coefficients form a smooth continuum from high variance and weak coupling to low variance and strong coupling. Our results are consistent with low rates of hybridization and a strong genome-wide barrier to gene flow when the coupling coefficient is much greater than 1, but also suggest that this boundary might be approached gradually and at a near constant rate over time.

A natural gene drive system confers reproductive isolation in rice

Hybrid sterility restricts the utilization of superior heterosis of indica-japonica inter-subspecific hybrids. In this study, we report the identification of RHS12, a major locus controlling male gamete sterility in indica-japonica hybrid rice. We show that RHS12 consists of two genes (iORF3/DUYAO and iORF4/JIEYAO) that confer preferential transmission of the RHS12-i type male gamete into the progeny, thereby forming a natural gene drive. DUYAO encodes a mitochondrion-targeted protein that interacts with OsCOX11 to trigger cytotoxicity and cell death, whereas JIEYAO encodes a protein that reroutes DUYAO to the autophagosome for degradation via direct physical interaction, thereby detoxifying DUYAO. Evolutionary trajectory analysis reveals that this system likely formed de novo in the AA genome Oryza clade and contributed to reproductive isolation (RI) between different lineages of rice. Our combined results provide mechanistic insights into the genetic basis of RI as well as insights for strategic designs of hybrid rice breeding.

The meaning of wild: Genetic and adaptive consequences from large-scale releases of domestic mallards

The translocation of individuals around the world is leading to rising incidences of anthropogenic hybridization, particularly between domestic and wild congeners. We apply a landscape genomics approach for thousands of mallard (Anas platyrhynchos) samples across continental and island populations to determine the result of over a century of supplementation practices. We establish that a single domestic game-farm mallard breed is the source for contemporary release programs in Eurasia and North America, as well as for established feral populations in New Zealand and Hawaii. In particular, we identify central Europe and eastern North America as epicenters of ongoing anthropogenic hybridization, and conclude that the release of game-farm mallards continues to affect the genetic integrity of wild mallards. Conversely, self-sustaining feral populations in New Zealand and Hawaii not only show strong differentiation from their original stock, but also signatures of local adaptation occurring in less than a half-century since game-farm mallard releases have ceased. We conclude that 'wild' is not singular, and that even feral populations are capable of responding to natural processes. Although considered paradoxical to biological conservation, understanding the capacity for wildness among feral and feral admixed populations in human landscapes is critical as such interactions increase in the Anthropocene.

Selection Maintains the Phenotypic Divergence of Cave and Surface Fish

AbstractGenetic divergence in the presence of gene flow has been well documented, but there is little information on the specific factors maintaining divergence. The present study investigates this in the Mexican tetra (Astyanax mexicanus), an excellent model for studying this question because surface and cave populations differ markedly in phenotype and genotype but are interfertile. Previous population studies documented significant gene flow among cave and surface populations, but they focused on analyses of neutral markers whose evolutionary dynamics likely differ from those of genes involved in cave adaptation. The present study advances our understanding of this question by focusing specifically on the genetics responsible for eye and pigmentation reduction, signature traits of cave populations. Direct observations of two cave populations over the course of 63 years verify that surface fish frequently move into the caves and even hybridize with the cave fish. Importantly, however, historical records show that surface alleles for pigmentation and eye size do not persist but are rapidly eliminated from the cave gene pool. It has been argued that regression of eyes and pigmentation was driven by drift, but the results of this study suggest that strong selection actively eliminates surface alleles from the cave populations.

Speciation and development

Understanding general principles about the origin of species remains one of the foundational challenges in evolutionary biology. The genomic divergence between groups of individuals can spawn hybrid inviability and hybrid sterility, which presents a tantalizing developmental problem. Divergent developmental programs may yield either conserved or divergent phenotypes relative to ancestral traits, both of which can be responsible for reproductive isolation during the speciation process. The genetic mechanisms of developmental evolution involve cis- and trans-acting gene regulatory change, protein-protein interactions, genetic network structures, dosage, and epigenetic regulation, all of which also have roots in population genetic and molecular evolutionary processes. Toward the goal of demystifying Darwin's "mystery of mysteries," this review integrates microevolutionary concepts of genetic change with principles of organismal development, establishing explicit links between population genetic process and developmental mechanisms in the production of macroevolutionary pattern. This integration aims to establish a more unified view of speciation that binds process and mechanism.

Genomic regions underlying the species-specific mating songs of green lacewings

Rapid species radiations provide insight into the process of speciation and diversification. The radiation of Chrysoperla carnea-group lacewings seems to be driven, at least in part, by their species-specific pre-mating vibrational duets. We associated genetic markers from across the genome with courtship song period in the offspring of a laboratory cross between Chrysoperla plorabunda and Chrysoperla adamsi, two species primarily differentiated by their mating songs. Two genomic regions were strongly associated with the song period phenotype. Large regions of chromosomes one and two were associated with song phenotype, as fewer recombination events occurred on these chromosomes relative to the other autosomes. Candidate genes were identified by functional annotation of proteins from the C. carnea reference genome. The majority of genes that are associated with vibrational courtship signals in other insects were found within QTL for lacewing song phenotype. Together these findings suggest that decreased recombination may be acting to keep together loci important to reproductive isolation between these species. Using wild-caught individuals from both species, we identified signals of genomic divergence across the genome. We identified several candidate genes both in song-associated regions and near divergence outliers including nonA, fruitless, paralytic, period, and doublesex. Together these findings bring us one step closer to identifying the genomic basis of a mating song trait critical to the maintenance of species boundaries in green lacewings.

Polygenic adaptation from standing genetic variation allows rapid ecotype formation

Adaptive ecotype formation can be the first step to speciation, but the genetic underpinnings of this process are poorly understood. Marine midges of the genus Clunio (Diptera) have recolonized Northern European shore areas after the last glaciation. In response to local tide conditions they have formed different ecotypes with respect to timing of adult emergence, oviposition behavior and larval habitat. Genomic analysis confirms the recent establishment of these ecotypes, reflected in massive haplotype sharing between ecotypes, irrespective of whether there is ongoing gene flow or geographic isolation. QTL mapping and genome screens reveal patterns of polygenic adaptation from standing genetic variation. Ecotype-associated loci prominently include circadian clock genes, as well as genes affecting sensory perception and nervous system development, hinting to a central role of these processes in ecotype formation. Our data show that adaptive ecotype formation can occur rapidly, with ongoing gene flow and largely based on a re-assortment of existing alleles.

Extensive gene flow in secondary sympatry after allopatric speciation

In the conventional view, species are separate gene pools delineated by reproductive isolation (RI). In an alternative view, species may also be delineated by a small set of 'speciation genes' without full RI, a view that has gained broad acceptance. A recent survey, however, suggested that the extensive literature on 'speciation with gene flow' is mostly (if not all) about exchanges in the early stages of speciation. There is no definitive evidence that the observed gene flow actually happened after speciation is completed. Here, we wish to know whether 'good species' (defined by the 'secondary sympatry' test) do continue to exchange genes and, importantly, under what conditions such exchanges can be observed. De novo whole-genome assembly and re-sequencing of individuals across the range of two closely related mangrove species (Rhizophora mucronata and R. stylosa) reveal the genomes to be well delineated in allopatry. They became sympatric in northeastern Australia but remain distinct species. Nevertheless, their genomes harbor ∼4000-10 000 introgression blocks averaging only about 3-4 Kb. These fine-grained introgressions indicate continual gene flow long after speciation as non-introgressable 'genomic islets,' ∼1.4 Kb in size, often harbor diverging genes of flower or gamete development. The fine-grained introgression in secondary sympatry may help settle the debate about sympatric vs. micro-allopatric speciation. In conclusion, true 'good species' may often continue to exchange genes but the opportunity for detection is highly constrained.

Meiotic drive is associated with sexual incompatibility in Neurospora

Evolution of Bateson-Dobzhansky-Muller (BDM) incompatibilities is thought to represent a key step in the formation of separate species. They are incompatible alleles that have evolved in separate populations and are exposed in hybrid offspring as hybrid sterility or lethality. In this study, we reveal a previously unconsidered mechanism promoting the formation of BDM incompatibilities, meiotic drive. Theoretical studies have evaluated the role that meiotic drive, the phenomenon whereby selfish elements bias their transmission to progeny at ratios above 50:50, plays in speciation, and have mostly concluded that drive could not result in speciation on its own. Using the model fungus Neurospora, we demonstrate that the large meiotic drive haplotypes, Sk-2 and Sk-3, contain putative sexual incompatibilities. Our experiments revealed that although crosses between Neurospora intermedia and Neurospora metzenbergii produce viable progeny at appreciable rates, when strains of N. intermedia carry Sk-2 or Sk-3 the proportion of viable progeny drops substantially. Additionally, it appears that Sk-2 and Sk-3 have accumulated different incompatibility phenotypes, consistent with their independent evolutionary history. This research illustrates how meiotic drive can contribute to reproductive isolation between populations, and thereby speciation.

The Genomic Landscapes of Desert Birds Form over Multiple Time Scales

Spatial models show that genetic differentiation between populations can be explained by factors ranging from geographic distance to environmental resistance across the landscape. However, genomes exhibit a landscape of differentiation, indicating that multiple processes may mediate divergence in different portions of the genome. We tested this idea by comparing alternative geographic predctors of differentiation in ten bird species that co-occur in Sonoran and Chihuahuan Deserts of North America. Using population-level genomic data, we described the genomic landscapes across species and modeled conditions that represented historical and contemporary mechanisms. The characteristics of genomic landscapes differed across species, influenced by varying levels of population structuring and admixture between deserts, and the best-fit models contrasted between the whole genome and partitions along the genome. Both historical and contemporary mechanisms were important in explaining genetic distance, but particularly past and current environments, suggesting that genomic evolution was modulated by climate and habitat There were also different best-ftit models across genomic partitions of the data, indicating that these regions capture different evolutionary histories. These results show that the genomic landscape of differentiation can be associated with alternative geographic factors operating on different portions of the genome, which reflect how heterogeneous patterns of genetic differentiation can evolve across species and genomes.

Widespread genomic signatures of reproductive isolation and sex-specific selection in the Eastern Yellow Robin, Eopsaltria australis

How new species evolve is one of the most fundamental questions in biology. Population divergence, which may lead to speciation, may be occurring in the Eastern Yellow Robin, a common passerine that lives along the eastern coast of Australia. This species is composed of 2 parapatric lineages that have highly divergent mitochondrial DNA; however, similar levels of divergence have not been observed in the nuclear genome. Here we re-examine the nuclear genomes of these mitolineages to test potential mechanisms underlying the discordance between nuclear and mitochondrial divergence. We find that nuclear admixture occurs in a narrow hybrid zone, although the majority of markers across the genome show evidence of reproductive isolation between populations of opposing mitolineages. There is an 8 MB section of a previously identified putative neo-sex chromosome that is highly diverged between allopatric but not parapatric populations, which may be the result of a chromosomal inversion. The neo-sex chromosomal nature of this region, as well as the geographic patterns in which it exhibits divergence, suggest it is unlikely to be contributing to reproductive isolation through mitonuclear incompatibilities as reported in earlier studies. In addition, there are sex differences in the number of markers that are differentiated between populations of opposite mitolineages, with greater differentiation occurring in females, which are heterozygous, than males. These results suggest that, despite the absence of previously observed assortative mating, mitolineages of Eastern Yellow Robin experience at least some postzygotic isolation from each other, in a pattern consistent with Haldane’s Rule.

Immunoglobulin genes, reproductive isolation and vertebrate speciation

Reproductive isolation drives the formation of new species, and many genes contribute to this through Dobzhansky–Muller incompatibilities (DMIs). These incompatibilities occur when gene divergence affects loci encoding interacting products such as receptors and their ligands. We suggest here that the nature of vertebrate immunoglobulin (IG) genes must make them prone to DMIs. The genes of these complex loci form functional genes through the process of recombination, giving rise to a repertoire of heterodimeric receptors of incredible diversity. This repertoire, within individuals and within species, must defend against pathogens but must also avoid pathogenic self‐reactivity. We suggest that this avoidance of autoimmunity is only achieved through a coordination of evolution between heavy‐ and light‐chain genes, and between these genes and the rest of the genome. Without coordinated evolution, the hybrid offspring of two diverging populations will carry a heavy burden of DMIs, resulting in a loss of fitness. Critical incompatibilities could manifest as incompatibilities between a mother and her divergent offspring. During fetal development, biochemical differences between the parents of hybrid offspring could make their offspring a target of the maternal immune system. This hypothesis was conceived in the light of recent insights into the population genetics of IG genes. This has suggested that antibody genes are probably as susceptible to evolutionary forces as other parts of the genome. Further repertoire studies in human and nonhuman species should now help determine whether antibody genes have been part of the evolutionary forces that drive the development of species.

Restricted X chromosome introgression and support for Haldane's rule in hybridizing damselflies

Contemporary hybrid zones act as natural laboratories for the investigation of species boundaries and may shed light on the little understood roles of sex chromosomes in species divergence. Sex chromosomes are considered to function as a hotspot of genetic divergence between species; indicated by less genomic introgression compared to autosomes during hybridization. Moreover, they are thought to contribute to Haldane's rule, which states that hybrids of the heterogametic sex are more likely to be inviable or sterile. To test these hypotheses, we used contemporary hybrid zones of Ischnura elegans, a damselfly species that has been expanding its range into the northern and western regions of Spain, leading to chronic hybridization with its sister species Ischnura graellsii. We analysed genome-wide SNPs in the Spanish I. elegans and I. graellsii hybrid zone and found (i) that the X chromosome shows less genomic introgression compared to autosomes, and (ii) that males are underrepresented among admixed individuals, as predicted by Haldane's rule. This is the first study in Odonata that suggests a role of the X chromosome in reproductive isolation. Moreover, our data add to the few studies on species with X0 sex determination system and contradict the hypothesis that the absence of a Y chromosome causes exceptions to Haldane's rule.

Rapid divergence of a gamete recognition gene promoted macroevolution of Eutheria

BACKGROUND

Speciation genes contribute disproportionately to species divergence, but few examples exist, especially in vertebrates. Here we test whether Zan, which encodes the sperm acrosomal protein zonadhesin that mediates species-specific adhesion to the egg's zona pellucida, is a speciation gene in placental mammals.

RESULTS

Genomic ontogeny reveals that Zan arose by repurposing of a stem vertebrate gene that was lost in multiple lineages but retained in Eutheria on acquiring a function in egg recognition. A 112-species Zan sequence phylogeny, representing 17 of 19 placental Orders, resolves all species into monophyletic groups corresponding to recognized Orders and Suborders, with <5% unsupported nodes. Three other rapidly evolving germ cell genes (Adam2, Zp2, and Prm1), a paralogous somatic cell gene (TectA), and a mitochondrial gene commonly used for phylogenetic analyses (Cytb) all yield trees with poorer resolution than the Zan tree and inferior topologies relative to a widely accepted mammalian supertree. Zan divergence by intense positive selection produces dramatic species differences in the protein's properties, with ordinal divergence rates generally reflecting species richness of placental Orders consistent with expectations for a speciation gene that acts across a wide range of taxa. Furthermore, Zan's combined phylogenetic utility and divergence exceeds those of all other genes known to have evolved in Eutheria by positive selection, including the only other mammalian speciation gene, Prdm9.

CONCLUSIONS

Species-specific egg recognition conferred by Zan's functional divergence served as a mode of prezygotic reproductive isolation that promoted the extraordinary adaptive radiation and success of Eutheria.

Genetic basis of speciation and adaptation: from loci to causative mutations

Does evolution proceed in small steps or large leaps? How repeatable is evolution? How constrained is the evolutionary process? Answering these long-standing questions in evolutionary biology is indispensable for both understanding how extant biodiversity has evolved and predicting how organisms and ecosystems will respond to changing environments in the future. Understanding the genetic basis of phenotypic diversification and speciation in natural populations is key to properly answering these questions. The leap forward in genome sequencing technologies has made it increasingly easier to not only investigate the genetic architecture but also identify the variant sites underlying adaptation and speciation in natural populations. Furthermore, recent advances in genome editing technologies are making it possible to investigate the functions of each candidate gene in organisms from natural populations. In this article, we discuss how these recent technological advances enable the analysis of causative genes and mutations and how such analysis can help answer long-standing evolutionary biology questions.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.

Allorecognition genes drive reproductive isolation in Podospora anserina

Allorecognition, the capacity to discriminate self from conspecific non-self, is a ubiquitous organismal feature typically governed by genes evolving under balancing selection. Here, we show that in the fungus Podospora anserina, allorecognition loci controlling vegetative incompatibility (het genes), define two reproductively isolated groups through pleiotropic effects on sexual compatibility. These two groups emerge from the antagonistic interactions of the unlinked loci het-r (encoding a NOD-like receptor) and het-v (encoding a methyltransferase and an MLKL/HeLo domain protein). Using a combination of genetic and ecological data, supported by simulations, we provide a concrete and molecularly defined example whereby the origin and coexistence of reproductively isolated groups in sympatry is driven by pleiotropic genes under balancing selection.

Sympatric speciation of the spiny mouse from Evolution Canyon in Israel substantiated genomically and methylomically

SignificanceWhether sympatric speciation (SS) is rare or common is still debated. Two populations of the spiny mouse, Acomys cahirinus, from Evolution Canyon I (EC I) in Israel have been depicted earlier as speciating sympatrically by molecular markers and transcriptome. Here, we investigated SS both genomically and methylomically, demonstrating that the opposite populations of spiny mice are sister taxa and split from the common ancestor around 20,000 years ago without an allopatric history. Mate choice, olfactory receptors, and speciation genes contributed to prezygotic/postzygotic reproductive isolation. The two populations showed different methylation patterns, facilitating adaptation to their local environment. They cope with abiotic and biotic stresses, due to high solar interslope radiation differences. We conclude that our new genomic and methylomic data substantiated SS.

Widespread introgression across a phylogeny of 155 Drosophila genomes

Genome-scale sequence data have invigorated the study of hybridization and introgression, particularly in animals. However, outside of a few notable cases, we lack systematic tests for introgression at a larger phylogenetic scale across entire clades. Here, we leverage 155 genome assemblies from 149 species to generate a fossil-calibrated phylogeny and conduct multilocus tests for introgression across 9 monophyletic radiations within the genus Drosophila. Using complementary phylogenomic approaches, we identify widespread introgression across the evolutionary history of Drosophila. Mapping gene-tree discordance onto the phylogeny revealed that both ancient and recent introgression has occurred across most of the 9 clades that we examined. Our results provide the first evidence of introgression occurring across the evolutionary history of Drosophila and highlight the need to continue to study the evolutionary consequences of hybridization and introgression in this genus and across the tree of life.

Comparing the genetic diversity and population structure of sister marine snails having contrasting habitat specificity

BACKGROUND

To grasp the processes of spatial genetic structuring in open and connectable marine environments is the principal study goal in molecular biological studies. Comparative seascape genetics using multiple species are a powerful approach to understand the physical geographic and oceanographic effects on genetic variation. Besides, species-specific ecological traits such as dispersal abilities and habitat specificity are important factors for spatial genetic structuring.

METHODS AND RESULTS

We focused on the sister marine snail species Tegula kusairo and T. xanthostigma around the Japanese mainland, which have contrasting habitat specificities for wave strength. Tegula kusairo only inhabits sheltered coastal environments, while T. xanthostigma is found mainly on wave-exposed rocky shores facing the open sea. We estimated their genetic diversity indices and levels of population differentiation based on mtDNA. We found that the genetic diversity of T. kusairo was lower than that of T. xanthostigma, while their level of population genetic differentiation was higher than that of T. xanthostigma. Namely, the species specific to weak wave environments had a higher level of population genetic differentiation than the species specific to strong wave action.

CONCLUSION

Ecological traits linked not only to dispersal abilities but also to habitat specificity can influence genetic variation in a pair of closely related sister species distributed in the same seascape.

Complex basis of hybrid female sterility and Haldane's rule in Heliconius butterflies: Z-linkage and epistasis

Hybrids between species are often sterile or inviable. Hybrid unfitness usually evolves first in the heterogametic sex-a pattern known as Haldane's rule. The genetics of Haldane's rule have been extensively studied in species where the male is the heterogametic (XX/XY) sex, but its basis in taxa where the female is heterogametic (ZW/ZZ), such as Lepidoptera and birds, is largely unknown. Here, we analyse a new case of female hybrid sterility between geographic subspecies of Heliconius pardalinus. The two subspecies mate freely in captivity, but female F1 hybrids in both directions of cross are sterile. Sterility is due to arrested development of oocytes after they become differentiated from nurse cells, but before yolk deposition. We backcrossed fertile male F1 hybrids to parental females and mapped quantitative trait loci (QTLs) for female sterility. We also identified genes differentially expressed in the ovary as a function of oocyte development. The Z chromosome has a major effect, similar to the 'large X effect' in Drosophila, with strong epistatic interactions between loci at either end of the Z chromosome, and between the Z chromosome and autosomal loci on chromosomes 8 and 20. By intersecting the list of genes within these QTLs with those differentially expressed in sterile and fertile hybrids, we identified three candidate genes with relevant phenotypes. This study is the first to characterize hybrid sterility using genome mapping in the Lepidoptera and shows that it is produced by multiple complex epistatic interactions often involving the sex chromosome, as predicted by the dominance theory of Haldane's rule.

P-elements strengthen reproductive isolation within the Drosophila simulans species complex

Determining mechanisms that underlie reproductive isolation (RI) is key to understanding how species boundaries are maintained in nature. Transposable elements (TEs) are ubiquitous across eukaryotic genomes. However, the role of TEs in modulating the strength of RI between species is poorly understood. Several species of Drosophila have been found to harbor P-elements (PEs), yet only D. simulans is known to be currently polymorphic for their presence in wild populations. PEs can cause RI between PE-containing (P) and PE-lacking (M) lineages of the same species. However, it is unclear whether they also contribute to the magnitude of RI between species. Here, we use the simulans species complex to assess whether differences in PE status between D. simulans and its sister species, which do not harbor PEs, contribute to multiple barriers to gene flow between species. We show that crosses involving a P D. simulans father and an M mother from a sister species exhibit lower F1 female fecundity than crosses involving an M D. simulans father and an M sister-species mother. We also find that another TE, I-element, might play a minor role in determining the frequency of dysgenesis between species. Our results suggest that the presence of PEs in a species can strengthen isolation from its sister species, providing evidence that TEs can play a role in RI.

Sex chromosomes in meiotic, hemiclonal, clonal and polyploid hybrid vertebrates: along the 'extended speciation continuum'

We review knowledge about the roles of sex chromosomes in vertebrate hybridization and speciation, exploring a gradient of divergences with increasing reproductive isolation (speciation continuum). Under early divergence, well-differentiated sex chromosomes in meiotic hybrids may cause Haldane-effects and introgress less easily than autosomes. Undifferentiated sex chromosomes are more susceptible to introgression and form multiple (or new) sex chromosome systems with hardly predictable dominance hierarchies. Under increased divergence, most vertebrates reach complete intrinsic reproductive isolation. Slightly earlier, some hybrids (linked in 'the extended speciation continuum') exhibit aberrant gametogenesis, leading towards female clonality. This facilitates the evolution of various allodiploid and allopolyploid clonal ('asexual') hybrid vertebrates, where 'asexuality' might be a form of intrinsic reproductive isolation. A comprehensive list of 'asexual' hybrid vertebrates shows that they all evolved from parents with divergences that were greater than at the intraspecific level (K2P-distances of greater than 5-22% based on mtDNA). These 'asexual' taxa inherited genetic sex determination by mostly undifferentiated sex chromosomes. Among the few known sex-determining systems in hybrid 'asexuals', female heterogamety (ZW) occurred about twice as often as male heterogamety (XY). We hypothesize that pre-/meiotic aberrations in all-female ZW-hybrids present Haldane-effects promoting their evolution. Understanding the preconditions to produce various clonal or meiotic allopolyploids appears crucial for insights into the evolution of sex, 'asexuality' and polyploidy. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Mass of genes rather than master genes underlie the genomic architecture of amphibian speciation

The genetic architecture of speciation, i.e., how intrinsic genomic incompatibilities promote reproductive isolation (RI) between diverging lineages, is one of the best-kept secrets of evolution. To directly assess whether incompatibilities arise in a limited set of large-effect speciation genes, or in a multitude of loci, we examined the geographic and genomic landscapes of introgression across the hybrid zones of 41 pairs of frog and toad lineages in the Western Palearctic region. As the divergence between lineages increases, phylogeographic transitions progressively become narrower, and larger parts of the genome resist introgression. This suggests that anuran speciation proceeds through a gradual accumulation of multiple barrier loci scattered across the genome, which ultimately deplete hybrid fitness by intrinsic postzygotic isolation, with behavioral isolation being achieved only at later stages. Moreover, these loci were disproportionately sex linked in one group (Hyla) but not in others (Rana and Bufotes), implying that large X-effects are not necessarily a rule of speciation with undifferentiated sex chromosomes. The highly polygenic nature of RI and the lack of hemizygous X/Z chromosomes could explain why the speciation clock ticks slower in amphibians compared to other vertebrates. The clock-like dynamics of speciation combined with the analytical focus on hybrid zones offer perspectives for more standardized practices of species delimitation.

Estimating hybridization rates in the wild: Easier said than done?

Hybridization has important effects on the evolutionary trajectories of natural populations but estimates of this process in the wild and at the individual-level are lacking. Justyn et al. attempted to fill this gap using the citizen science database eBird but there are limitations to this approach. Here, we outline and directly test these limitations using literature searches, case studies, and a comparison between eBird and Birds of North America (BNA), a database that documents hybridization using the scientific literature. We use a hybrid zone between Lazuli and Indigo buntings to highlight the importance of considering geographic range when estimating rates of hybridization and two literature searches to show the importance of considering cryptic hybrids (those that cannot be identified using phenotypic traits) when quantifying these rates. We also use BNA and a case study of hybrid White-faced and Glossy Ibises to show that citizen scientists are underreporting hybrids compared with experts. Justyn et al. highlighted an important gap in the literature, but their results likely represent the lower limit of hybridization between birds and a more nuanced interpretation of their results (e.g., considering extrinsic postzygotic selection) is needed.

Isolation and reconnection: Demographic history and multiple contact zones of the green odorous frog (Odorrana margaretae) around the Sichuan Basin

The green odorous frog (Odorrana margaretae) displays a circular distribution around the Sichuan Basin of western China and possesses multiple replicate hybrid zones between lineages with high levels of divergence. To gain an understanding of the speciation process, we obtained 1540 SNPs from 29 populations and 227 individuals using ddRAD sequencing. Population structure analysis revealed three groups within the species: the West, North & South, and East groups. Demographic inference showed that they were initially isolated at ~2 million years ago, and subsequent post-glacial expansion produced the current circular distribution with four secondary contact zones. Hybridization in those zones involved lineages with various levels of divergence and produced greatly different patterns of introgression. Contact zones between the East and North & South groups (E-S and E-N) had contrast admixture levels but both showed a general lack of potential barrier loci. Meanwhile, the reconnection of the West and North & South groups produced two contact zones along the rim of the Basin. The S-W zone had extensive admixture while the N-W zone had limited admixture within a narrow geographic distance. Both showed substantial barrier effects, and a large number of potential barrier loci were shared. We also detected strong coupling among these loci. The N-W hybrid zone involved two highly-diverged lineages (F_ST = 0.704) and many loci have reached fixation around the hybrid zone. This study system offers a unique opportunity to understand the dynamics of introgression in contact zones and the architecture of reproductive isolation at different stages of speciation.

Integrative assessment of intraspecific diversification in Loggerhead Shrike (Lanius ludovicianus) provides insight on the geographic pattern of phenotypic divergence and process of speciation

Integrated studies of the geographical, ecological, and historical factors that shape intraspecific phenotypic and genetic variation can help us to decipher the processes leading to geographic patterns of population divergence and speciation. We quantify and compare morphological and genetic variation in the Loggerhead Shrike (Lanius ludovicianus Linnaeus, 1766), a broadly distributed passerine in North America with both migratory and non-migratory populations that occupy a diversity of habitats and topographies. The geographic distributions and patterns of differentiation among subspecies suggest that migration has strongly impacted population divergence, including the habit of migrating itself, but also dispersal. Patterns of mitochondrial and nuclear genetic differentiation can be attributed to female-biased dispersal and to increased dispersal rates and distances in migratory populations. Weak phenotypic differentiation among migratory versus migratory and non-migratory populations suggest that migration may more strongly affect morphology than adaptation to local habitats. Our results generally support previous subspecific designations with two notable exceptions. We found little genetic differentiation between two subspecies (Lanius ludovicianus gambeli Ridgway, 1887 and Lanius ludovicianus mexicanus C.L. Brehm, 1854), but identify a new, distinct subspecies, which we refer to as Lanius ludovicianus centralis ssp. nov.

The Origin of Niches and Species in the Bacterial World

Niches are spaces for the biological units of selection, from cells to complex communities. In a broad sense, "species" are biological units of individuation. Niches do not exist without individual organisms, and every organism has a niche. We use "niche" in the Hutchinsonian sense as an abstraction of a multidimensional environmental space characterized by a variety of conditions, both biotic and abiotic, whose quantitative ranges determine the positive or negative growth rates of the microbial individual, typically a species, but also parts of the communities of species contained in this space. Microbial organisms ("species") constantly diversify, and such diversification (radiation) depends on the possibility of opening up unexploited or insufficiently exploited niches. Niche exploitation frequently implies "niche construction," as the colonized niche evolves with time, giving rise to new potential subniches, thereby influencing the selection of a series of new variants in the progeny. The evolution of niches and organisms is the result of reciprocal interacting processes that form a single unified process. Centrifugal microbial diversification expands the limits of the species' niches while a centripetal or cohesive process occurs simultaneously, mediated by horizontal gene transfers and recombinatorial events, condensing all of the information recovered during the diversifying specialization into "novel organisms" (possible future species), thereby creating a more complex niche, where the selfishness of the new organism(s) establishes a "homeostatic power" limiting the niche's variation. Once the niche's full carrying capacity has been reached, reproductive isolation occurs, as no foreign organisms can outcompete the established population/community, thereby facilitating speciation. In the case of individualization-speciation of the microbiota, its contribution to the animal' gut structure is a type of "niche construction," the result of crosstalk between the niche (host) and microorganism(s). Lastly, there is a parallelism between the hierarchy of niches and that of microbial individuals. The increasing anthropogenic effects on the biosphere (such as globalization) might reduce the diversity of niches and bacterial individuals, with the potential emergence of highly transmissible multispecialists (which are eventually deleterious) resulting from the homogenization of the microbiosphere, a possibility that should be explored and prevented.

CYP2J19 mediates carotenoid colour introgression across a natural avian hybrid zone

It has long been of interest to identify the phenotypic traits that mediate reproductive isolation between related species, and more recently, the genes that underpin them. Much work has focused on identifying genes associated with animal colour, with the candidate gene CYP2J19 identified in laboratory studies as the ketolase converting yellow dietary carotenoids to red ketocarotenoids in birds with red pigments. However, evidence that CYP2J19 explains variation between red and yellow feather coloration in wild populations of birds is lacking. Hybrid zones provide the opportunity to identify genes associated with specific traits. Here we investigate genomic regions associated with colour in red-fronted and yellow-fronted tinkerbirds across a hybrid zone in southern Africa. We sampled 85 individuals, measuring spectral reflectance of forecrown feathers and scoring colours from photographs, while testing for carotenoid presence with Raman spectroscopy. We performed a genome-wide association study to identify associations with carotenoid-based coloration, using double-digest RAD sequencing aligned to a short-read whole genome of a Pogoniulus tinkerbird. Admixture mapping using 104,933 single nucleotide polymorphisms (SNPs) identified a region of chromosome 8 that includes CYP2J19 as the only locus with more than two SNPs significantly associated with both crown hue and crown score, while Raman spectra provided evidence of ketocarotenoids in red feathers. Asymmetric backcrossing in the hybrid zone suggests that yellow-fronted females mate more often with red-fronted males than vice versa. Female red-fronted tinkerbirds mating assortatively with red-crowned males is consistent with the hypothesis that converted carotenoids are an honest signal of quality.

A single pleiotropic locus influences the rate of hybridization between two sibling species of Lygaeus bugs

The evolution of reproductive isolation lies at the heart of understanding the process of speciation. Of particular interest is the relationship between pre- and postzygotic reproductive isolation, and the genetic architecture of traits that contribute to one or both forms of reproductive isolation. The sibling species of seed bug Lygaeus equestris and L. simulans show a classic pattern of asymmetric prezygotic reproductive isolation, with female L. equestris hybridizing with male L. simulans, but with no hybridization in the reciprocal direction. We have recently described a mutant pale color form of L. simulans, that inherits as a single Mendelian locus and is pleiotropic for a number of other life history and behavioral traits. Here, we tested whether this locus also influences pre- and postzygotic reproductive isolation. Two sets of experimental crosses revealed that behavioral isolation varied with mutant versus wild-type phenotype for male L. simulans, with the pale form less successful at mating with female L. equestris. In terms of trying to assess postzygotic isolation, levels of hybrid offspring production were uniformly low across the experiments. However, we did obtain, for the first time, hybrid offspring from a pairing between a female L. simulans and a male L. equestris. In this instance, the female was of the pale mutant genotype. Together with evidence for heterozygote advantage in terms of nymph survival, we consider our results in terms of possible mechanisms of reproductive isolation between this species pair, the role of the pale mutation, and the possible genetic architectures underlying the mutation, from a single gene to a supergene.

Emerging Frontiers in the Study of Molecular Evolution

A collection of the editors of Journal of Molecular Evolution have gotten together to pose a set of key challenges and future directions for the field of molecular evolution. Topics include challenges and new directions in prebiotic chemistry and the RNA world, reconstruction of early cellular genomes and proteins, macromolecular and functional evolution, evolutionary cell biology, genome evolution, molecular evolutionary ecology, viral phylodynamics, theoretical population genomics, somatic cell molecular evolution, and directed evolution. While our effort is not meant to be exhaustive, it reflects research questions and problems in the field of molecular evolution that are exciting to our editors.

Hybridization patterns between two marine snails, Littorina fabalis and L. obtusata

Characterizing the patterns of hybridization between closely related species is crucial to understand the role of gene flow in speciation. In particular, systems comprising multiple contacts between sister species offer an outstanding opportunity to investigate how reproductive isolation varies with environmental conditions, demography and geographic contexts of divergence. The flat periwinkles, Littorina obtusata and L. fabalis (Gastropoda), are two intertidal sister species with marked ecological differences compatible with late stages of speciation. Although hybridization between the two was previously suggested, its extent across the Atlantic shores of Europe remained largely unknown. Here, we combined genetic (microsatellites and mtDNA) and morphological data (shell and male genital morphology) from multiple populations of flat periwinkles in north-western Iberia to assess the extent of current and past hybridization between L. obtusata and L. fabalis under two contrasting geographic settings of divergence (sympatry and allopatry). Hybridization signatures based on both mtDNA and microsatellites were stronger in sympatric sites, although evidence for recent extensive admixture was found in a single location. Misidentification of individuals into species based on shell morphology was higher in sympatric than in allopatric sites. However, despite hybridization, species distinctiveness based on this phenotypic trait together with male genital morphology remained relatively high. The observed variation in the extent of hybridization among locations provides a rare opportunity for future studies on the consequences of different levels of gene flow for reinforcement, thus informing about the mechanisms underlying the completion of speciation.

Assessing changes in genomic divergence following a century of human-mediated secondary contact among wild and captive-bred ducks

Along with manipulating habitat, the direct release of domesticated individuals into the wild is a practice used worldwide to augment wildlife populations. We test between possible outcomes of human-mediated secondary contact using genomic techniques at both historical and contemporary timescales for two iconic duck species. First, we sequence several thousand ddRAD-seq loci for contemporary mallards (Anas platyrhynchos) throughout North America and two domestic mallard types (i.e., known game-farm mallards and feral Khaki Campbell's). We show that North American mallards may well be becoming a hybrid swarm due to interbreeding with domesticated game-farm mallards released for hunting. Next, to attain a historical perspective, we applied a bait-capture array targeting thousands of loci in century-old (1842-1915) and contemporary (2009-2010) mallard and American black duck (Anas rubripes) specimens. We conclude that American black ducks and mallards have always been closely related, with a divergence time of ~600,000 years before present, and likely evolved through prolonged isolation followed by limited bouts of gene flow (i.e., secondary contact). They continue to maintain genetic separation, a finding that overturns decades of prior research and speculation suggesting the genetic extinction of the American black duck due to contemporary interbreeding with mallards. Thus, despite having high rates of hybridization, actual gene flow is limited between mallards and American black ducks. Conversely, our historical and contemporary data confirm that the intensive stocking of game-farm mallards during the last ~100 years has fundamentally changed the genetic integrity of North America's wild mallard population, especially in the east. It thus becomes of great interest to ask whether the iconic North American mallard is declining in the wild due to introgression of maladaptive traits from domesticated forms. Moreover, we hypothesize that differential gene flow from domestic game-farm mallards into the wild mallard population may explain the overall temporal increase in differentiation between wild black ducks and mallards, as well as the uncoupling of genetic diversity and effective population size estimates across time in our results. Finally, our findings highlight how genomic methods can recover complex population histories by capturing DNA preserved in traditional museum specimens.

Genetic divergence outpaces phenotypic evolution among threespine stickleback populations in old freshwater habitats

Species such as threespine stickleback (Gasterosteus aculeatus) that inhabit divergent selective environments and that have diversified on different time scales can be of value for understanding evolutionary processes. Here we synthesize high-resolution genotypic and phenotypic data to explore a largely unstudied distribution of threespine stickleback populations living in oceanic and freshwater habitats along coastal and inland regions of Oregon. Many inland aquatic habitats of Oregon remained unglaciated during the last ice age, meaning that some extant Oregon lake and river stickleback may have descended from freshwater populations established long before more well-studied, post-glacial freshwater populations. To address the degree of congruence between genetic and phenotypic divergence, we directly compared Oregon stickleback to much younger (post-glacial) Alaskan populations. We found phenotypic variation in Oregon stickleback to be primarily partitioned between oceanic and freshwater habitats, as has been documented in other stickleback systems. However, the main axis of genetic divergence was between coastal and inland regions regardless of habitat type. Furthermore, when comparing patterns between Oregon and Alaska we found similar levels of phenotypic divergence, but much greater genetic divergence among Oregon's populations. The Oregon stickleback system therefore appears well suited for future studies linking genotypic and phenotypic change, further extending the utility of this small fish to provide general insights into evolutionary processes.

Genome-wide association analyses reveal polygenic genomic architecture underlying divergent shell morphology in Spanish Littorina saxatilis ecotypes

Gene flow between diverging populations experiencing dissimilar ecological conditions can theoretically constrain adaptive evolution. To minimize the effect of gene flow, alleles underlying traits essential for local adaptation are predicted to be located in linked genome regions with reduced recombination. Local reduction in gene flow caused by selection is expected to produce elevated divergence in these regions. The highly divergent crab-adapted and wave-adapted ecotypes of the marine snail Littorina saxatilis present a model system to test these predictions. We used genome-wide association (GWA) analysis of geometric morphometric shell traits associated with microgeographic divergence between the two L. saxatilis ecotypes within three separate sampling sites. A total of 477 snails that had individual geometric morphometric data and individual genotypes at 4,066 single nucleotide polymorphisms (SNPs) were analyzed using GWA methods that corrected for population structure among the three sites. This approach allowed dissection of the genomic architecture of shell shape divergence between ecotypes across a wide geographic range, spanning two glacial lineages. GWA revealed 216 quantitative trait loci (QTL) with shell size or shape differences between ecotypes, with most loci explaining a small proportion of phenotypic variation. We found that QTL were evenly distributed across 17 linkage groups, and exhibited elevated interchromosomal linkage, suggesting a genome-wide response to divergent selection on shell shape between the two ecotypes. Shell shape trait-associated loci showed partial overlap with previously identified outlier loci under divergent selection between the two ecotypes, supporting the hypothesis of diversifying selection on these genomic regions. These results suggest that divergence in shell shape between the crab-adapted and wave-adapted ecotypes is produced predominantly by a polygenic genomic architecture with positive linkage disequilibrium among loci of small effect.

HD-AGPs as Speciation Genes: Positive Selection on a Proline-Rich Domain in Non-Hybridizing Species of Petunia, Solanum, and Nicotiana

Transmitting tissue-specific proteins (TTS proteins) are abundant in the extracellular matrix of Nicotiana pistils, and vital for optimal pollen tube growth and seed set. We have identified orthologs from several species in the Solanaceae, including Petunia axillaris axillaris and Petunia integrifolia. We refer to TTS proteins and their orthologs as histidine domain-arabinogalactan proteins (HD-AGPs). HD-AGPs have distinctive domains, including a small histidine-rich region and a C-terminal PAC domain. Pairwise comparisons between HD-AGPs of 15 species belonging to Petunia, Nicotiana, and Solanum show that the his-domain and PAC domain are under purifying selection. In contrast, a proline-rich domain (HV2) is conserved among cross-hybridizing species, but variant in species-pairs that are reproductively isolated by post-pollination pre-fertilization reproductive barriers. In particular, variation in a tetrapeptide motif (XKPP) is systematically correlated with the presence of an interspecific reproductive barrier. Ka/Ks ratios are not informative at the infrageneric level, but the ratios reveal a clear signature of positive selection on two hypervariable domains (HV1 and HV2) when HD-AGPs from five solanaceous genera are compared. We propose that sequence divergence in the hypervariable domains of HD-AGPs reinforces sympatric speciation in incipient species that may have first diverged as a consequence of pollinator preferences or other ecological factors.

Searching for Sympatric Speciation in the Genomic Era

Sympatric speciation illustrates how natural and sexual selection may create new species in isolation without geographic barriers. However, recent genomic reanalyses of classic examples of sympatric speciation reveal complex histories of secondary gene flow from outgroups into the radiation. In contrast, the rich theoretical literature on this process distinguishes among a diverse range of models based on simple genetic histories and different types of reproductive isolating barriers. Thus, there is a need to revisit how to connect theoretical models of sympatric speciation and their predictions to empirical case studies in the face of widespread gene flow. Here, theoretical differences among different types of sympatric speciation and speciation-with-gene-flow models are reviewed and summarized, and genomic analyses are proposed for distinguishing which models apply to case studies based on the timing and function of adaptive introgression. Investigating whether secondary gene flow contributed to reproductive isolation is necessary to test whether predictions of theory are ultimately borne out in nature.

Genetic, phenotypic and ecological differentiation suggests incipient speciation in two Charadrius plovers along the Chinese coast

Background

Speciation with gene flow is an alternative to the nascence of new taxa in strict allopatric separation. Indeed, many taxa have parapatric distributions at present. It is often unclear if these are secondary contacts, e.g. caused by past glaciation cycles or the manifestation of speciation with gene flow, which hampers our understanding of how different forces drive diversification. Here we studied genetic, phenotypic and ecological aspects of divergence in a pair of incipient shorebird species, the Kentish (Charadrius alexandrinus) and the White-faced Plovers (C. dealbatus), shorebirds with parapatric breeding ranges along the Chinese coast. We assessed divergence based on molecular markers with different modes of inheritance and quantified phenotypic and ecological divergence in aspects of morphometric, dietary and climatic niches.

Results

Our integrative analyses revealed small to moderate levels of genetic and phenotypic distinctiveness with symmetric gene flow across the contact area at the Chinese coast. The two species diverged approximately half a million years ago in dynamic isolation with secondary contact occurring due to cycling sea level changes between the Eastern and Southern China Sea in the mid-late Pleistocene. We found evidence of character displacement and ecological niche differentiation between the two species, invoking the role of selection in facilitating divergence despite gene flow.

Conclusion

These findings imply that ecology can indeed counter gene flow through divergent selection and thus contributes to incipient speciation in these plovers. Furthermore, our study highlights the importance of using integrative datasets to reveal the evolutionary history and assist the inference of mechanisms of speciation.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1449-5) contains supplementary material, which is available to authorized users.

An asymmetric allelic interaction drives allele transmission bias in interspecific rice hybrids

Hybrid sterility (HS) between Oryza sativa (Asian rice) and O. glaberrima (African rice) is mainly controlled by the S1 locus. However, our limited understanding of the HS mechanism hampers utilization of the strong interspecific heterosis. Here, we show that three closely linked genes (S1A4, S1TPR, and S1A6) in the African S1 allele (S1-g) constitute a killer-protector system that eliminates gametes carrying the Asian allele (S1-s). In Asian–African rice hybrids (S1-gS1-s), the S1TPR-S1A4-S1A6 interaction in sporophytic tissues generates an abortion signal to male and female gametes. However, S1TPR can rescue S1-g gametes, while the S1-s gametes selectively abort for lacking S1TPR. Knockout of any of the S1-g genes eliminates the HS. Evolutionary analysis suggests that S1 may have arisen from newly evolved genes, multi-step recombination, and nucleotide variations. Our findings will help to overcome the interspecific reproductive barrier and use Asian–African hybrids for increasing rice production.

Our limited understanding of the hybrid sterility (HS) mechanism in Asian–African rice hybrids hampers utilization of the interspecific heterosis for rice production. Here, the authors identify S1-mediated HS-related tripartite gamete killer-protector system, and explore their evolutionary relationship.