Seeking signatures of reinforcement at the genetic level: a hitchhiking mapping and candidate gene approach in the house mouse

Reinforcement is the process by which prezygotic isolation is strengthened as a response to selection against hybridization. Most empirical support for reinforcement comes from the observation of its possible phenotypic signature: an accentuated degree of prezygotic isolation in the hybrid zone as compared to allopatry. Here, we implemented a novel approach to this question by seeking for the signature of reinforcement at the genetic level. In the house mouse, selection against hybrids and enhanced olfactory‐based assortative mate preferences are observed in a hybrid zone between the two European subspecies Mus musculus musculus and M. m. domesticus, suggesting a possible recent reinforcement event. To test for the genetic signature of reinforcing selection and identify genes involved in sexual isolation, we adopted a hitchhiking mapping approach targeting genomic regions containing candidate genes for assortative mating in mice. We densely scanned these genomic regions in hybrid zone and allopatric samples using a large number of fast evolving microsatellite loci that allow the detection of recent selection events. We found a handful of loci showing the expected pattern of significant reduction in variability in populations close to the hybrid zone, showing assortative odour preference in mate choice experiments as compared to populations further away and displaying no such preference. These loci lie close to genes that we pinpoint as testable candidates for further investigation.     

Bayesian estimation of genomic clines

We developed a Bayesian genomic cline model to study the genetic architecture of adaptive divergence and reproductive isolation between hybridizing lineages. This model quantifies locus‐specific patterns of introgression with two cline parameters that describe the probability of locus‐specific ancestry as a function of genome‐wide admixture. ‘Outlier’ loci with extreme patterns of introgression relative to most of the genome can be identified. These loci are potentially associated with adaptive divergence or reproductive isolation. We simulated genetic data for admixed populations that included neutral introgression, as well as loci that were subject to directional, epistatic or underdominant selection, and analysed these data using the Bayesian genomic cline model. Under many demographic conditions, underdominance or directional selection had detectable and predictable effects on cline parameters, and ‘outlier’ loci were greatly enriched for genetic regions affected by selection. We also analysed previously published genetic data from two transects through a hybrid zone between Mus domesticus and M. musculus. We found considerable variation in rates of introgression across the genome and particularly low rates of introgression for two X‐linked markers. There were similarities and differences in patterns of introgression between the two transects, which likely reflects a combination of stochastic variability because of genetic drift and geographic variation in the genetic architecture of reproductive isolation. By providing a robust framework to quantify and compare patterns of introgression among genetic regions and populations, the Bayesian genomic cline model will advance our understanding of the genetics of reproductive isolation and the speciation process.     

Speciation is not necessarily easier in species with sexually monomorphic mating signals

Should we have different expectations regarding the likelihood and pace of speciation by sexual selection when considering species with sexually monomorphic mating signals? Two conditions that can facilitate rapid species divergence are Felsenstein's one‐allele mechanism and a genetic architecture that includes a genetic association between signal and preference loci. In sexually monomorphic species, the former can manifest in the form of mate choice based on phenotype matching. The latter can be promoted by selection acting upon genetic loci for divergent signals and preferences expressed simultaneously in each individual, rather than acting separately on signal loci in males and preference loci in females. Both sexes in the Chrysoperla carnea group of green lacewings (Insecta, Neuroptera, Chrysopidae) produce sexually monomorphic species‐specific mating signals. We hybridized the two species C. agilis and C. carnea to test for evidence of these speciation‐facilitating conditions. Hybrid signals were more complex than the parents and we observed a dominant influence of C. carnea. We found a dominant influence of C. agilis on preferences in the form of hybrid discrimination against C. carnea. Preferences in hybrids followed patterns predicting preference loci that determine mate choice rather than a one‐allele mechanism. The genetic association between signal and preference we detected in the segregating hybrid crosses indicates that speciation in these species with sexually monomorphic mating signals can have occurred rapidly. However, we need additional evidence to determine whether such genetic associations form more readily in sexually monomorphic species compared to dimorphic species and consequently facilitate speciation.     

The genetics of speciation: genes of small effect underlie sexual isolation in the Hawaiian cricket Laupala

Sexual behaviours often evolve rapidly and are critical for sexual isolation. We suggest that coordinated sexual signals and preferences generate stabilizing selection, favouring the accumulation of many small‐effect mutations in sexual communication traits. Rapid radiation of a sexual behaviour used in signalling, song pulse rate, has been observed in the Hawaiian cricket genus Laupala. Using marker‐assisted introgression, we isolated five known quantitative trait loci (QTL) influencing species‐level differences in pulse rate from one species, L. paranigra, into a closely related species, L. kohalensis. All five QTL were found to have a significant effect on song and appear to be largely additive in backcross introgression lines. Furthermore, all effect sizes were small in magnitude. Our data provide support for the hypothesis that stabilizing selection on sexual signals in Laupala creates genetic conditions favourable to incremental divergence during speciation, through the evolution of alleles of minor rather than major phenotypic effects.     

COMPONENTS OF REPRODUCTIVE ISOLATION BETWEEN NORTH AMERICAN PHEROMONE STRAINS OF THE EUROPEAN CORN BORER

Of 12 potential reproductive isolating barriers between closely related Z‐ and E‐pheromone strains of the European corn borer moth (Ostrinia nubilalis), seven significantly reduced gene flow but none were complete, suggesting that speciation in this lineage is a gradual process in which multiple barriers of intermediate strength accumulate. Estimation of the cumulative effect of all barriers resulted in nearly complete isolation (>99%), but geographic variation in seasonal isolation allowed as much as ～10% gene flow. With the strongest barriers arising from mate‐selection behavior or ecologically relevant traits, sexual and natural selection are the most likely evolutionary processes driving population divergence. A recent multilocus genealogical study corroborates the roles of selection and gene flow ( Dopman et al. 2005 ), because introgression is supported at all loci besides Tpi, a sex‐linked gene. Tpi reveals strains as exclusive groups, possesses signatures of selection, and is tightly linked to a QTL that contributes to seasonal isolation. With more than 98% of total cumulative isolation consisting of prezygotic barriers, Z and E strains of ECB join a growing list of taxa in which species boundaries are primarily maintained by the prevention of hybridization, possibly because premating barriers evolve during early stages of population divergence.     

Differential introgression causes genealogical discordance in host races of Acrocercops transecta (Insecta: Lepidoptera)

Recently diverged populations often exhibit incomplete reproductive isolation, with a low level of gene flow continuing between populations. Previous studies have shown that, even under a low level of gene flow, genetic divergence between populations can proceed at the loci governing local adaptation and reproductive isolation but not at other neutral loci. A leaf‐mining moth, Acrocercops transecta, consists of Juglans‐ and Lyonia‐associated host races. The two host races differ in host preferences of ovipositing females and in larval adaptation to host plants but mate readily in the laboratory, producing fertile hybrids. The Juglans and Lyonia races are often sympatric in the wild, implying that gene introgression could occur in nature between the two host races. We tested this hypothesis by combining phylogenetic analyses with coalescent simulations, focusing on mitochondrial genes (COI and ND5) and the nuclear Tpi, Per and Ldh genes located on the Z‐chromosome. The mitochondrial genes clearly distinguished the Lyonia race from the Juglnas race, whereas the Tpi, Per and Ldh genealogies did not reflect the two host races. Coalescent simulations indicated gene flow at the three Z‐linked genes in both directions, whereas there was no introgression in the mitochondrial genes. The lack of introgression in mitochondrial genes suggests that female host preference is the primary force leading to the bifurcation of maternally inherited loci. Thus, the results show that a low level of gene flow coupled with the inflexible female host preference differentiates histories of divergence between maternally and biparentally inherited genes in this host race system.     

Effects of assortative mate choice on the genomic and morphological structure of a hybrid zone between two bird subspecies

Phenotypic differentiation plays an important role in the formation and maintenance of reproductive barriers. In some cases, variation in a few key aspects of phenotype can promote and maintain divergence; hence, the identification of these traits and their associations with patterns of genomic divergence is crucial for understanding the patterns and processes of population differentiation. We studied hybridization between the alba and personata subspecies of the white wagtail (Motacilla alba), and quantified divergence and introgression of multiple morphological traits and 19,437 SNP loci on a 3,000 km transect. Our goal was to identify traits that may contribute to reproductive barriers and to assess how variation in these traits corresponds to patterns of genome‐wide divergence. Variation in only one trait—head plumage patterning—was consistent with reproductive isolation. Transitions in head plumage were steep and occurred over otherwise morphologically and genetically homogeneous populations, whereas cline centres for other traits and genomic ancestry were displaced over 100 km from the head cline. Field observational data show that social pairs mated assortatively by head plumage, suggesting that these phenotypes are maintained by divergent mating preferences. In contrast, variation in all other traits and genetic markers could be explained by neutral diffusion, although weak ecological selection cannot be ruled out. Our results emphasize that assortative mating may maintain phenotypic differences independent of other processes shaping genome‐wide variation, consistent with other recent findings that raise questions about the relative importance of mate choice, ecological selection and selectively neutral processes for divergent evolution.     

Reinforcement shapes clines in female mate discrimination in Drosophila subquinaria

Reinforcement of species boundaries may alter mate recognition in a way that also affects patterns of mate preference among conspecific populations. In the fly Drosophila subquinaria, females sympatric with the closely related species D. recens reject mating with heterospecific males as well as with conspecific males from allopatric populations. Here, we assess geographic variation in behavioral isolation within and among populations of D. subquinaria and use cline theory to understand patterns of selection on reinforced discrimination and its consequences for sexual isolation within species. We find that selection has fixed rejection of D. recens males in sympatry, while significant genetic variation in this behavior occurs within allopatric populations. In conspecific matings sexual isolation is also asymmetric and stronger in populations that are sympatric with D. recens. The clines in behavioral discrimination within and between species are similar in shape and are maintained by strong selection in the face of gene flow, and we show that some of their genetic basis may be either shared or linked. Thus, while reinforcement can drive extremely strong phenotypic divergence, the long‐term consequences for incipient speciation depend on gene flow, genetic linkage of discrimination traits, and the cost of these behaviors in allopatry.     

Physical linkage and mate preference generate linkage disequilibrium for behavioral isolation in two parapatric crickets

Behavioral isolation is a potent barrier to gene flow and a source of striking diversity in the animal kingdom. However, it remains unclear if the linkage disequilibrium (LD) between sex‐specific traits required for behavioral isolation results mostly from physical linkage between signal and preference loci or from directional mate preferences. Here, we test this in the field crickets Gryllus rubens and G. texensis. These closely related species diverged with gene flow and have strongly differentiated songs and preference functions for the mate calling song rhythm. We map quantitative trait loci for signal and preference traits (pQTL) as well as for gene expression associated with these traits (eQTL). We find strong, positive genetic covariance between song traits and between song and preference. Our results show that this is in part explained by incomplete physical linkage: although both linked pQTL and eQTL couple male and female traits, major effect loci for different traits were never on the same chromosome. We suggest that the finely tuned, highly divergent preference functions are likely an additional source of LD between male and female traits in this system. Furthermore, pleiotropy of gene expression presents an underappreciated mechanism to link sexually dimorphic phenotypes.     

The role of local ecology during hybridization at the initial stages of ecological speciation in a marine snail

Hybrid zones of ecologically divergent populations are ideal systems to study the interaction between natural selection and gene flow during the initial stages of speciation. Here, we perform an amplified fragment length polymorphism (AFLP) genome scan in parallel hybrid zones between divergent ecotypes of the marine snail Littorina saxatilis, which is considered a model case for the study of ecological speciation. Ridged‐Banded (RB) and Smooth‐Unbanded (SU) ecotypes are adapted to different shore levels and microhabitats, although they present a sympatric distribution at the mid‐shore where they meet and mate (partially assortatively). We used shell morphology, outlier and nonoutlier AFLP loci from RB, SU and hybrid specimens captured in sympatry to determine the level of phenotypic and genetic introgression. We found different levels of introgression at parallel hybrid zones and nonoutlier loci showed more gene flow with greater phenotypic introgression. These results were independent from the phylogeography of the studied populations, but not from the local ecological conditions. Genetic variation at outlier loci was highly correlated with phenotypic variation. In addition, we used the relationship between genetic and phenotypic variation to estimate the heritability of morphological traits and to identify potential Quantitative Trait Loci to be confirmed in future crosses. These results suggest that ecology (exogenous selection) plays an important role in this hybrid zone. Thus, ecologically based divergent natural selection is responsible, simultaneously, for both ecotype divergence and hybridization. On the other hand, genetic introgression occurs only at neutral loci (nonoutliers). In the future, genome‐wide studies and controlled crosses would give more valuable information about this process of speciation in the face of gene flow.     

Multiple forms of selection shape reproductive isolation in a primate hybrid zone

Speciation occurs when populations diverge and become reproductively isolated from each other. Natural selection is commonly accepted to play a large role in this process, and it has been widely assumed that reproductive isolation often results as a by‐product of divergence driven by adaptation in allopatry. When such populations come into secondary contact, reinforcement can act to strengthen reproductive isolation, but the frequency and importance of this process are still unknown. Here, we explored genomic signatures of selection in allopatry and sympatry for loci associated with reproductive isolation using a natural primate hybrid zone. By analysing reduced‐representation sequencing data, we quantified admixture and population structure across a howler monkey hybrid zone and examined the relationship between locus‐specific differentiation and introgression. We detected extensive admixture that was mostly limited to the narrow contact zone. Loci with reduced introgression into the heterospecific genomic background (the pattern expected for loci associated with reproductive isolation due to selection against hybrids) were significantly more differentiated between allopatric parental populations than loci with neutral and increased introgression, supporting the hypothesis that reproductive isolation is a by‐product of divergence in allopatry. Further, loci with reduced introgression showed greater differentiation in sympatry than in allopatry, suggesting a role for reinforcement. Thus, our results reflect multiple forms of selection that have shaped reproductive isolation in this system. We conclude that reproductive isolation may have initially been driven by divergence in allopatry, but later reinforced by divergent selection in sympatry.     

Do changes in gene expression contribute to sexual isolation and reinforcement in the house mouse?

Expression divergence, rather than sequence divergence, has been shown to be important in speciation, particularly in the early stages of divergence of traits involved in reproductive isolation. In the two European subspecies of house mice, Mus musculus musculus and Mus musculus domesticus, earlier studies have demonstrated olfactory‐based assortative mate preference in populations close to their hybrid zone. It has been suggested that this behaviour evolved following the recent secondary contact between the two taxa (~3,000 years ago) in response to selection against hybridization. To test for a role of changes in gene expression in the observed behavioural shift, we conducted a RNA sequencing experiment on mouse vomeronasal organs. Key candidate genes for pheromone‐based subspecies recognition, the vomeronasal receptors, are expressed in these organs. Overall patterns of gene expression varied significantly between samples from the two subspecies, with a large number of differentially expressed genes between the two taxa. In contrast, only ~200 genes were found repeatedly differentially expressed between populations within M. m. musculus that did or did not display assortative mate preferences (close to or more distant from the hybrid zone, respectively), with an overrepresentation of genes belonging to vomeronasal receptor family 2. These receptors are known to play a key role in recognition of chemical cues that handle information about genetic identity. Interestingly, four of five of these differentially expressed receptors belong to the same phylogenetic cluster, suggesting specialization of a group of closely related receptors in the recognition of odorant signals that may allow subspecies recognition and assortative mating.     

Host shifts and signal divergence: mating signals covary with host use in a complex of specialized plant‐feeding insects

A combination of divergent natural and sexual selection is a powerful cause of speciation. This conjunction of evolutionary forces may often occur when divergence is initiated by ecological differences between populations because local adaptation to new resources can lead to changes in sexual selection. The hypothesis that differences in resource use contribute to the evolution of reproductive isolation by altering the nature of sexual selection predicts that: (1) differences in sexual traits, such as signals and preferences, are an important source of reproductive isolation between species using different resources; (2) there are identifiable sources of selection on sexual traits that differ between species using different resources; and (3) signals vary between populations using different resources to a larger extent than between populations using the same resource at different localities. Testing these predictions requires a group of closely‐related species or populations that specialize on different resources and for which the traits involved in mate choice are known. The Enchenopa binotata species complex of treehoppers (Hemiptera: Membracidae) are host plant specialists in which speciation is associated with shifts to novel host plants. Mating in this complex is preceded by an exchange of vibrational signals transmitted through host plant stems, and the signal traits important for mate choice have been identified. In the E. binotata complex, previous work has supported the first two predictions: (1) signal differences between species are important in mate recognition and (2) host shifts can alter both the trait values favoured by sexual selection and the evolutionary response to that selection. In the present study, we tested the last prediction by conducting a large‐scale study of mating signal variation within and between the 11 species in the complex. We find that differences in host use are strongly associated with differences in signal traits important for mate recognition. This result supports the hypothesis that hosts shifts have led to speciation in this group in part through their influence on divergence in mate communication systems. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 60–72.     

Insight into the roles of selection in speciation from genomic patterns of divergence and introgression in secondary contact in venomous rattlesnakes

Investigating secondary contact of historically isolated lineages can provide insight into how selection and drift influence genomic divergence and admixture. Here, we studied the genomic landscape of divergence and introgression following secondary contact between lineages of the Western Diamondback Rattlesnake (Crotalus atrox) to determine whether genomic regions under selection in allopatry also contribute to reproductive isolation during introgression. We used thousands of nuclear loci to study genomic differentiation between two lineages that have experienced recent secondary contact following isolation, and incorporated sampling from a zone of secondary contact to identify loci that are resistant to gene flow in hybrids. Comparisons of patterns of divergence and introgression revealed a positive relationship between allelic differentiation and resistance to introgression across the genome, and greater‐than‐expected overlap between genes linked to lineage‐specific divergence and loci that resist introgression. Genes linked to putatively selected markers were related to prominent aspects of rattlesnake biology that differ between populations of Western Diamondback rattlesnakes (i.e., venom and reproductive phenotypes). We also found evidence for selection against introgression of genes that may contribute to cytonuclear incompatibility, consistent with previously observed biased patterns of nuclear and mitochondrial alleles suggestive of partial reproductive isolation due to cytonuclear incompatibilities. Our results provide a genome‐scale perspective on the relationships between divergence and introgression in secondary contact that is relevant for understanding the roles of selection in maintaining partial isolation of lineages, causing admixing lineages to not completely homogenize.     

Testing for selection on the androgen-binding protein in the Danish mouse hybrid zone

The three peripheral subspecies of the house mouse Mus musculus have fixed specific variants of the androgen-binding protein (ABP) that have been proposed to be part of a recognition mechanism that could participate in sexual isolation between the subspecies. We tested for selection on Abpa by characterizing the pattern of Abpa introgression across a transect of the hybrid zone between M. m. musculus and M. m. domesticus in Jutland. On the musculus side, the cline for Abpa resembled that of a nearly neutral allozyme more than that of strongly selected X and Y chromosome markers. However, the high central step which displaces the tail of introgression of Abpa to higher frequencies was best accounted for by linkage to a locus under strong selection against hybrids. Still, we cannot exclude that this pattern results from weak selection against Abpa in the tail of introgression, which would be compatible with an assortative choice mechanism. On the domesticus side there was little introgression close to the hybrid zone, presumably due to a geographical barrier to migration. However, substantial frequencies of musculus alleles occurred further away, suggesting mixed colonization patterns as well as flow across the hybrid zone.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 447–459.     

Insect mating signal and mate preference phenotypes covary among host plant genotypes

Sexual selection acting on small initial differences in mating signals and mate preferences can enhance signal–preference codivergence and reproductive isolation during speciation. However, the origin of initial differences in sexual traits remains unclear. We asked whether biotic environments, a source of variation in sexual traits, may provide a general solution to this problem. Specifically, we asked whether genetic variation in biotic environments provided by host plants can result in signal–preference phenotypic covariance in a host‐specific, plant‐feeding insect. We used a member of the Enchenopa binotata species complex of treehoppers (Hemiptera: Membracidae) to assess patterns of variation in male mating signals and female mate preferences induced by genetic variation in host plants. We employed a novel implementation of a quantitative genetics method, rearing field‐collected treehoppers on a sample of naturally occurring replicated host plant clone lines. We found remarkably high signal–preference covariance among host plant genotypes. Thus, genetic variation in biotic environments influences the sexual phenotypes of organisms living on those environments in a way that promotes assortative mating among environments. This consequence arises from conditions likely to be common in nature (phenotypic plasticity and variation in biotic environments). It therefore offers a general answer to how divergent sexual selection may begin.     

Divergent host plant adaptation drives the evolution of sexual isolation in the grasshopper Hesperotettix viridis (Orthoptera: Acrididae) in the absence of reinforcement

Early stages of lineage divergence in insect herbivores are often related to shifts in host plant use and divergence in mating capabilities, which may lead to sexual isolation of populations of herbivorous insects. We examined host preferences, degree of differentiation in mate choice, and divergence in cuticular morphology using near‐infrared spectroscopy in the grasshopper Hesperotettix viridis aiming to understand lineage divergence. In Kansas (USA), H. viridis is an oligophagous species feeding on Gutierrezia and Solidago host species. To identify incipient mechanisms of lineage divergence and isolation, we compared host choice, mate choice, and phenotypic divergence among natural grasshopper populations in zones of contact with populations encountering only one of the host species. A significant host‐based preference from the two host groups was detected in host‐paired feeding preference studies. No‐choice mate selection experiments revealed a preference for individuals collected from the same host species independent of geographic location, and little mating was observed between individuals collected from different host species. Female mate choice tests between males from the two host species resulted in 100% fidelity with respect to host use. Significant differentiation in colour and cuticular composition of individuals from different host plants was observed, which correlated positively with host choice and mate choice. No evidence for reinforcement in the zone of contact was detected, suggesting that divergent selection for host plant use promotes sexual isolation in this species. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 866–878.     

Reinforcement of mate preference among hybridizing Heliconius butterflies

Recent models of mate preference evolution suggest that direct selection on alleles at preference loci and correlated evolution of preference with locally adapted mating cues are more likely to drive the evolution of assortative mate preference than reinforcement. Mate preference evolution in mimetic Heliconius butterflies has been attributed to all three forms of selection, but here we show that reinforcement has been critical. By examining geographical variation in assortative mating and male mate preference among seven populations of three hybridizing Heliconius species from Costa Rica, we found pronounced character displacement of preference such that sexual isolation was enhanced in areas of interspecific contact. Of the different explanations for the evolution of assortative mate preference, only reinforcement is dependent on interspecific contact in this system. Thus, the observed pattern of reproductive character displacement of mate preference is best explained as a product of indirect selection generated by natural selection against nonmimetic hybrids.     

Divergence in coloration and ecological speciation in the Anolis marmoratus species complex

Adaptive divergence in coloration is expected to produce reproductive isolation in species that use colourful signals in mate choice and species recognition. Indeed, many adaptive radiations are characterized by differentiation in colourful signals, suggesting that divergent selection acting on coloration may be an important component of speciation. Populations in the Anolis marmoratus species complex from the Caribbean island of Guadeloupe display striking divergence in the colour and pattern of adult males that occurs over small geographic distances, suggesting strong divergent selection. Here we test the hypothesis that divergence in coloration results in reduced gene flow among populations. We quantify variation in adult male coloration across a habitat gradient between mesic and xeric habitats, use a multilocus coalescent approach to infer historical demographic parameters of divergence, and examine gene flow and population structure using microsatellite variation. We find that colour variation evolved without geographic isolation and in the face of gene flow, consistent with strong divergent selection and that both ecological and sexual selection are implicated. However, we find no significant differentiation at microsatellite loci across populations, suggesting little reproductive isolation and high levels of contemporary gene exchange. Strong divergent selection on loci affecting coloration probably maintains clinal phenotypic variation despite high gene flow at neutral loci, supporting the notion of a porous genome in which adaptive portions of the genome remain fixed whereas neutral portions are homogenized by gene flow and recombination. We discuss the impact of these findings for studies of colour evolution and ecological speciation.