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

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

The shape of things to come in the study of the origin of species?

Perhaps Darwin would agree that speciation is no longer the mystery of mysteries that it used to be. It is now generally accepted that evolution by natural selection can contribute to ecological adaptation, resulting in the evolution of reproductive barriers and, hence, to the evolution of new species (Schluter & Conte 2009 ; Meyer 2011 ; Nosil 2012 ). From genes that encode silencing proteins that cause infertility in hybrid mice (Mihola et al. 2009 ), to segregation distorters linked to speciation in fruit flies (Phadnis & Orr 2009 ), or pollinator‐mediated selection on flower colour alleles driving reinforcement in Texan wildflowers (Hopkins & Rausher 2012 ), characterization of the genes that drive speciation is providing clues to the origin of species (Nosil & Schluter 2011 ). It is becoming apparent that, while recent work continues to overturn historical ideas about sympatric speciation (e.g. Barluenga et al. 2006 ), ecological circumstances strongly influence patterns of genomic divergence, and ultimately the establishment of reproductive isolation when gene flow is present (Elmer & Meyer 2011 ). Less clear, however, are the genetic mechanisms that cause speciation, particularly when ongoing gene flow is occurring. Now, in this issue, Franchini et al. ( 2014 ) employ a classic genetic mapping approach augmented with new genomic tools to elucidate the genomic architecture of ecologically divergent body shapes in a pair of sympatric crater lake cichlid fishes. From over 450 segregating SNPs in an F2 cross, 72 SNPs were linked to 11 QTL associated with external morphology measured by means of traditional and geometric morphometrics. Annotation of two highly supported QTL further pointed to genes that might contribute to ecological divergence in body shape in Midas cichlids, overall supporting the hypothesis that genomic regions of large phenotypic effect may be contributing to early‐stage divergence in Midas cichlids.     

Putative chromosomal rearrangements are associated primarily with ecotype divergence rather than geographic separation in an intertidal,poorly dispersing snail

Littorina saxatilis is becoming a model system for understanding the genomic basis of ecological speciation. The parallel formation of crab‐adapted ecotypes that exhibit partial reproductive isolation from wave‐adapted ecotypes has enabled genomic investigation of conspicuous shell traits. Recent genomic studies suggest that chromosomal rearrangements may enable ecotype divergence by reducing gene flow. However, the genomic architecture of traits that are divergent between ecotypes remains poorly understood. Here, we use 11,504 single nucleotide polymorphism (SNP) markers called using the recently released L. saxatilis genome to genotype 462 crab ecotype, wave ecotype and phenotypically intermediate Spanish L. saxatilis individuals with scored phenotypes. We used redundancy analysis to study the genetic architecture of loci associated with shell shape, shape corrected for size, shell size and shell ornamentation, and to compare levels of co‐association among different traits. We discovered 341 SNPs associated with shell traits. Loci associated with trait divergence between ecotypes were often located inside putative chromosomal rearrangements recently characterized in Swedish L. saxatilis. In contrast, we found that shell shape corrected for size varied primarily by geographic site rather than by ecotype and showed little association with these putative rearrangements. We conclude that genomic regions of elevated divergence inside putative rearrangements were associated with divergence of L. saxatilis ecotypes along steep environmental axes—consistent with models of adaptation with gene flow—but were not associated with divergence among the three geographical sites. Our findings support predictions from models indicating the importance of genomic regions of reduced recombination allowing co‐association of loci during ecological speciation with ongoing gene flow.     

Persistent postmating,prezygotic reproductive isolation between populations

Studying reproductive barriers between populations of the same species is critical to understand how speciation may proceed. Growing evidence suggests postmating, prezygotic (PMPZ) reproductive barriers play an important role in the evolution of early taxonomic divergence. However, the contribution of PMPZ isolation to speciation is typically studied between species in which barriers that maintain isolation may not be those that contributed to reduced gene flow between populations. Moreover, in internally fertilizing animals, PMPZ isolation is related to male ejaculate—female reproductive tract incompatibilities but few studies have examined how mating history of the sexes can affect the strength of PMPZ isolation and the extent to which PMPZ isolation is repeatable or restricted to particular interacting genotypes. We addressed these outstanding questions using multiple populations of Drosophila montana. We show a recurrent pattern of PMPZ isolation, with flies from one population exhibiting reproductive incompatibility in crosses with all three other populations, while those three populations were fully fertile with each other. Reproductive incompatibility is due to lack of fertilization and is asymmetrical, affecting female fitness more than males. There was no effect of male or female mating history on reproductive incompatibility, indicating that PMPZ isolation persists between populations. We found no evidence of variability in fertilization outcomes attributable to different female × male genotype interactions, and in combination with our other results, suggests that PMPZ isolation is not driven by idiosyncratic genotype × genotype interactions. Our results show PMPZ isolation as a strong, consistent barrier to gene flow early during speciation and suggest several targets of selection known to affect ejaculate‐female reproductive tract interactions within species that may cause this PMPZ isolation.     

Patterns of divergence across the geographic and genomic landscape of a butterfly hybrid zone associated with a climatic gradient

Hybrid zones are a valuable tool for studying the process of speciation and for identifying the genomic regions undergoing divergence and the ecological (extrinsic) and nonecological (intrinsic) factors involved. Here, we explored the genomic and geographic landscape of divergence in a hybrid zone between Papilio glaucus and Papilio canadensis. Using a genome scan of 28,417 ddRAD SNPs, we identified genomic regions under possible selection and examined their distribution in the context of previously identified candidate genes for ecological adaptations. We showed that differentiation was genomewide, including multiple candidate genes for ecological adaptations, particularly those involved in seasonal adaptation and host plant detoxification. The Z chromosome and four autosomes showed a disproportionate amount of differentiation, suggesting genes on these chromosomes play a potential role in reproductive isolation. Cline analyses of significantly differentiated genomic SNPs, and of species‐diagnostic genetic markers, showed a high degree of geographic coincidence (81%) and concordance (80%) and were associated with the geographic distribution of a climate‐mediated developmental threshold (length of the growing season). A relatively large proportion (1.3%) of the outliers for divergent selection were not associated with candidate genes for ecological adaptations and may reflect the presence of previously unrecognized intrinsic barriers between these species. These results suggest that exogenous (climate‐mediated) and endogenous (unknown) clines may have become coupled and act together to reinforce reproductive isolation. This approach of assessing divergence across both the genomic and geographic landscape can provide insight about the interplay between the genetic architecture of reproductive isolation and endogenous and exogenous selection.     

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

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

Understanding the genetic basis of invasiveness

Invasive species provide excellent study systems to evaluate the ecological and evolutionary processes that contribute to the colonization of novel environments. While the ecological processes that contribute to the successful establishment of invasive plants have been studied in detail, investigation of the evolutionary processes involved in successful invasions has only recently received attention. In particular, studies investigating the genomic and gene expression differences between native and introduced populations of invasive species are just beginning and are required if we are to understand how plants become invasive. In the current issue of Molecular Ecology, Hodgins et al. ( 2013 ) tackle this unresolved question, by examining gene expression differences between native and introduced populations of annual ragweed, Ambrosia artemisiifolia. The study identifies a number of potential candidate genes based on gene expression differences that may be responsible for the success of annual ragweed in its introduced range. Furthermore, genes involved in stress response are over‐represented in the differentially expressed gene set. Future experiments could use functional studies to test whether changes in gene expression at these candidate genes do in fact underlie changes in growth characteristics and reproductive output observed in this and other invasive species.     

Sexual conflict and the maintenance of genetic variation in natural populations

Understanding the factors that maintain genetic variation in natural populations is a foundational goal of evolutionary biology. To this end, population geneticists have developed a variety of models that can produce stable polymorphisms. In one of the earliest models, Owen ( 1953 ) demonstrated that differences in selection pressures acting on males and females could maintain multiple alleles of a gene at a stable equilibrium. If the selection pressures act in opposite directions in males and females, we refer to this as (inter‐) sexual conflict or sexual antagonism (Arnqvist & Rowe, 2005 ). Testing if sexual conflict maintains genetic variation in natural populations is a tremendous challenge—it requires both identifying loci that harbor sexually antagonistic alleles and determining whether those alleles are maintained as stable polymorphisms (Mank, 2017 ). Doing so genome‐wide is even harder because it is not tractable to identify sexually antagonistic alleles and test for stable polymorphisms at all loci. Dutoit et al. ( 2018 ) confront this challenge in a paper published in this issue of Molecular Ecology. Using gene expression and population genomic data from the collared flycatcher, Dutoit et al. ( 2018 ) identify associations and correlations between genomic signatures of balanced polymorphisms and sexual conflict.     

Strong reproductive isolation and narrow genomic tracts of differentiation among three woodpecker species in secondary contact

Hybrid zones allow the measurement of gene flow across the genome, producing insight into the genomic architecture of speciation. Such analysis is particularly powerful when applied to multiple pairs of hybridizing species, as patterns of genomic differentiation can then be related to age of the hybridizing species, providing a view into the build‐up of differentiation over time. We examined 33 809 single nucleotide polymorphisms (SNPs) in three hybridizing woodpecker species: Red‐breasted, Red‐naped and Yellow‐bellied sapsuckers (Sphyrapicus ruber, Sphyrapicus nuchalis and Sphyrapicus varius), two of which (ruber and nuchalis) are much more closely related than each is to the third (varius). To identify positions of SNPs on chromosomes, we developed a localization method based on comparative genomics. We found narrow clines, bimodal distributions of hybrid indices and genomic regions with decreased rates of introgression. These results suggest moderately strong reproductive isolation among species and selection against specific hybrid genotypes. We found 19 small regions of strong differentiation between species, partly shared among species pairs, but no large regions of differentiation. An association analysis revealed a single strong‐effect candidate locus associated with plumage, possibly explaining mismatch among the three species in genomic relatedness and plumage similarity. Our comparative analysis of species pairs of different age and their hybrid zones showed that moderately strong reproductive isolation can occur with little genomic differentiation, but that reproductive isolation is incomplete even with much greater genomic differentiation, implying there are long periods of time when hybridization is possible if diverging populations are in geographic contact.     

Vocal divergence is concordant with genomic evidence for strong reproductive isolation in grasshopper mice (Onychomys)

Behavioral barriers to gene flow often evolve faster than intrinsic incompatibilities and can eliminate the opportunity for hybridization between interfertile species. While acoustic signal divergence is a common driver of premating isolation in birds and insects, its contribution to speciation in mammals is less studied. Here we characterize the incidence of, and potential barriers to, hybridization among three closely related species of grasshopper mice (genus Onychomys). All three species use long‐distance acoustic signals to attract and localize mates; Onychomys arenicola and Onychomys torridus are acoustically similar and morphologically cryptic whereas Onychomys leucogaster is larger and acoustically distinct. We used genotyping‐by‐sequencing (GBS) to test for evidence of introgression in 227 mice from allopatric and sympatric localities in the western United States and northern Mexico. We conducted laboratory mating trials for all species pairs to assess reproductive compatibility, and recorded vocalizations from O. arenicola and O. torridus in sympatry and allopatry to test for evidence of acoustic character displacement. Hybridization was rare in nature and, contrary to prior evidence for O. torridus/O. arenicola hybrids, only involved O. leucogaster and O. arenicola. In contrast, laboratory crosses between O. torridus and O. arenicola produced litters whereas O. leucogaster and O. arenicola crosses did not. Call fundamental frequency in O. torridus and O. arenicola was indistinguishable in allopatry but significantly differentiated in sympatry, a pattern consistent with reproductive character displacement. These results suggest that assortative mating based on a long‐distance signal is an important isolating mechanism between O. torridus and O. arenicola and highlight the importance of behavioral barriers in determining the permeability of species boundaries.     

Geographic patterns of postzygotic isolation between two closely related widespread dung fly species (Sepsis cynipsea and Sepsis neocynipsea; Diptera: Sepsidae)

Identifying the contribution of pre‐ and postzygotic barriers to gene flow is a key goal of speciation research. The widespread dung fly species Sepsis cynipsea and Sepsis neocynipsea offer great potential for studying the speciation process over a range of opportunities for gene exchange within and across sister species (cross‐continental allopatry, continental parapatry and sympatry). We examined the role of postcopulatory isolating barriers by comparing female fecundity and egg‐to‐adult viability of F₁ and F₂ hybrids, as well as backcrosses of F₁ hybrids with the parental species, via replicated crosses of sym‐, para‐ and allopatric populations. Egg‐to‐adult viability was strongly but not totally suppressed in hybrids, and offspring production approached nil in the F₂ generation (hybrid breakdown), indicating yet unspecified intrinsic incompatibilities. Viable F₁ hybrid offspring showed almost absolute male (the heterogametic sex) sterility while females remained largely fertile, in accordance with Haldane's rule. Hybridization between the two species in European areas of sympatry (Swiss Alps) indicated only minor reinforcement based on fecundity traits. Crossing geographically isolated European and North American S. neocynipsea showed similar albeit weaker isolating barriers that are most easily explained by random genetic drift. We conclude that in this system with a biogeographic continuum of reproductive barriers, speciation is mediated primarily by genetic drift following dispersal of flies over a wide (allopatric) geographic range, with some role of natural or sexual selection in incidental or direct reinforcement of incompatibility mechanisms in areas of European sympatry. S(ubs)pecies status of continental S. neocynipsea appears warranted.     

Population genomics and speciation in yeasts

The advent of efficient whole genome sequencing and the large molecular and genetic toolbox available for studies in Saccharomyces cerevisiae and related species have allowed unprecedented analysis of big issues such as: what causes reproductive isolation and eventual speciation? The species complex encompassing S. cerevisiae and relatives consists of six species and several naturally occurring hybrids, which have nearly collinear genomes. They fit the biological species definition of within species fertility and between species sterility. There are examples of chromosome rearrangements and of genetic incompatibilities between species of the complex, which contribute to reproductive isolation but these are not universally present. In addition, simple sequence divergence has been shown to cause reproductive isolation via the action of the mismatch repair system. Although all three of these mechanisms contribute to extant reproductive isolation, which if any, drive the speciation process is still an open question. Population genomic surveys of whole genome sequences reveal introgressions and horizontal gene transfers between species, indicating that the species barriers are not complete. This gene flow between species, although infrequent, brings into question the nature of yeast species.     

Mapping the genomic architecture of ecological speciation in the wild: does linkage disequilibrium hold the key?

The hunt for the genes underlying ecological speciation has now closed in on a number of candidates, but making the link from genotype to phenotype continues to pose a significant challenge. This is partly because genetic studies in many systems remain impeded by long generation times or an inability to perform controlled crosses. Now, in this issue of Molecular Ecology, Malek et al. (2012) demonstrate the utility of a novel admixture mapping approach that can be used to identify genomic regions contributing to adaptive trait divergence between natural populations. Remarkably, they validate their approach by mapping traits associated with mate choice in a wild limnetic and benthic threespine stickleback (Gasterosteus aculeatus) species pair, finding several loci associated with male nuptial coloration and shape. While this study benefited from tried‐and‐true microsatellites in a well‐characterized species with a detailed genetic map (and genome sequence), the field is quickly moving towards the use of next‐generation sequencing, especially for nonmodel systems. The ability to characterize molecular polymorphisms for any system suggests that molecular ecologists working on virtually any species may benefit from applying Malek et al.'s approach, if naturally admixed populations are available.     

QTL‐seq for rapid identification of candidate genes for 100‐seed weight and root/total plant dry weight ratio under rainfed conditions in chickpea

Terminal drought is a major constraint to chickpea productivity. Two component traits responsible for reduction in yield under drought stress include reduction in seeds size and root length/root density. QTL‐seq approach, therefore, was used to identify candidate genomic regions for 100‐seed weight (100SDW) and total dry root weight to total plant dry weight ratio (RTR) under rainfed conditions. Genomewide SNP profiling of extreme phenotypic bulks from the ICC 4958 × ICC 1882 population identified two significant genomic regions, one on CaLG01 (1.08 Mb) and another on CaLG04 (2.7 Mb) linkage groups for 100SDW. Similarly, one significant genomic region on CaLG04 (1.10 Mb) was identified for RTR. Comprehensive analysis revealed four and five putative candidate genes associated with 100SDW and RTR, respectively. Subsequently, two genes (Ca_04364 and Ca_04607) for 100SDW and one gene (Ca_04586) for RTR were validated using CAPS/dCAPS markers. Identified candidate genomic regions and genes may be useful for molecular breeding for chickpea improvement.