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

Reticulate evolution and rapid development of reproductive barriers upon secondary contact in a forest fungus

Reproductive barriers between sister species of the mushroom-forming fungi tend to be stronger in sympatry, leading to speculation on whether they are being reinforced by selection against hybrids. We have used population genomic analyses together with in vitro crosses of a global sample of the wood decay fungus Trichaptum abietinum to investigate reproductive barriers within this species complex and the processes that have shaped them. Our phylogeographic analyses show that T. abietinum is delimited into six major genetic groups: one in Asia, two in Europe, and three in North America. The groups present in Europe are interfertile and admixed, whereas our crosses show that the North American groups are reproductively isolated. In Asia, a more complex pattern appears, with partial intersterility between subgroups that likely originated independently and more recently than the reproductive barriers in North America. We found pre-mating barriers in T. abietinum to be moderately correlated with genomic divergence, whereas mean growth reduction of the mated hybrids showed a strong correlation with increasing genomic divergence. Genome-wide association analyses identified candidate genes with programmed cell death annotations, which are known to be involved in intersterility in distantly related fungi, although their link here remains unproven. Our demographic modeling and phylogenetic network analyses fit a scenario where reproductive barriers in Trichaptum abietinum could have been reinforced upon secondary contact between groups that diverged in allopatry during the Pleistocene glacial cycles. Our combination of experimental and genomic approaches demonstrates how T. abietinum is a tractable system for studying speciation mechanisms.

Beyond gene flow: (non)-parallelism of secondary contact in a pair of highly differentiated sibling species

Replicated secondary contact zones can provide insights into the barriers to gene flow that are important during speciation and can reveal to which degree secondary contact may result in similar evolutionary outcomes. Here, we studied two secondary contact zones between highly differentiated Alpine butterflies of the genus Erebia using whole-genome resequencing data. We assessed the genomic relationships between populations and species and found hybridization to be rare, with no to little current or historical introgression in either contact zone. There are large similarities between contact zones, consistent with an allopatric origin of interspecific differentiation, with no indications for ongoing reinforcing selection. Consistent with expected reduced effective population size, we further find that scaffolds related to the Z-chromosome show increased differentiation compared to the already high levels across the entire genome, which could also hint towards a contribution of the Z chromosome to species divergence in this system. Finally, we detected the presence of the endosymbiont Wolbachia, which can cause reproductive isolation between its hosts, in all E. cassioides, while it appears to be fully or largely absent in contact zone populations of E. tyndarus. We discuss how this rare pattern may have arisen and how it may have affected the dynamics of speciation upon secondary contact.

Genomic analysis of the rhesus macaque (Macaca mulatta) and the cynomolgus macaque (Macaca fascicularis) uncover polygenic signatures of reinforcement speciation

Speciation can involve phases of divergent adaptation in allopatry and ecological/reproductive character displacement in sympatry or parapatry. Reproductive character displacement can result as a means of preventing hybridization, a process known as reinforcement speciation. In this study, we use whole-genome sequencing (WGS) of two closely related primate species that have experienced introgression in their history, the rhesus (Macaca mulatta) and cynomolgus (M. fascicularis) macaques, to identify genes exhibiting reproductive character displacement and other patterns consistent with reinforcement speciation. Using windowed scans of various population genetic statistics to identify signatures of reinforcement, we find 184 candidate genes associated with a variety of functions, including an overrepresentation of multiple neurological functions and several genes involved in sexual development and gametogenesis. These results are consistent with a variety of genes acting in a reinforcement process between these species. We also find signatures of introgression of the Y-chromosome that confirm previous studies suggesting male-driven introgression of M. mulatta into M. fascicularis populations. This study uses WGS to find evidence of the process of reinforcement in primates that have medical and conservation relevance.

The population genetics of speciation by cascade reinforcement

Species interactions drive diverse evolutionary outcomes. Speciation by cascade reinforcement represents one example of how species interactions can contribute to the proliferation of species. This process occurs when the divergence of mating traits in response to selection against interspecific hybridization incidentally leads to reproductive isolation among populations of the same species. Here, we investigated the population genetic outcomes of cascade reinforcement in North American chorus frogs (Hylidae: Pseudacris). Specifically, we estimated the frequency of hybridization among three taxa, assessed genetic structure within the focal species, P. feriarum, and ascertained the directionality of gene flow within P. feriarum across replicated contact zones via coalescent modeling. Through field observations and preliminary experimental crosses, we assessed whether hybridization is possible under natural and laboratory conditions. We found that hybridization occurs among P. feriarum and two conspecifics at a low rate in multiple contact zones, and that gene flow within the former species is unidirectional from allopatry into sympatry with these other species in three of four contact zones studied. We found evidence of substantial genetic structuring within P. feriarum including a divergent western allopatric cluster, a behaviorally-distinct sympatric South Carolina cluster, and several genetically-overlapping clusters from the remainder of the distribution. Furthermore, we found sub-structuring between reinforced and nonreinforced populations in the two most intensely-sampled contact zones. Our literature review indicated that P. feriarum hybridizes with at least five heterospecifics at the periphery of its range providing a mechanism for further intraspecific diversification. This work strengthens the evidence for cascade reinforcement in this clade, revealing the geographic and genetic landscape upon which this process can contribute to the proliferation of species.

Regulation of Sexually Dimorphic Expression of Major Urinary Proteins

Male house mice excrete large amounts of protein in their urinary scent marks, mainly composed of Major Urinary Proteins (MUPs), and these lipocalins function as pheromones and pheromone carriers. Here, we review studies on sexually dimorphic MUP expression in house mice, including the proximate mechanisms controlling MUP gene expression and their adaptive functions. Males excrete 2 to 8 times more urinary protein than females, though there is enormous variation in gene expression across loci in both sexes. MUP expression is dynamically regulated depending upon a variety of factors. Males regulate MUP expression according to social status, whereas females do not, and males regulate expression depending upon health and condition. Male-biased MUP expression is regulated by pituitary secretion of growth hormone (GH), which binds receptors in the liver, activating the JAK2-STAT5 signaling pathway, chromatin accessibility, and MUP gene transcription. Pulsatile male GH secretion is feminized by several factors, including caloric restriction, microbiota depletion, and aging, which helps explain condition-dependent MUP expression. If MUP production has sex-specific fitness optima, then this should generate sexual antagonism over allelic expression (intra-locus sexual conflict) selectively favoring sexually dimorphic expression. MUPs influence the sexual attractiveness of male urinary odor and increased urinary protein excretion is correlated with the reproductive success of males but not females. This finding could explain the selective maintenance of sexually dimorphic MUP expression. Producing MUPs entails energetic costs, but increased excretion may reduce the net energetic costs and predation risks from male scent marking as well as prolong the release of chemical signals. MUPs may also provide physiological benefits, including regulating metabolic rate and toxin removal, which may have sex-specific effects on survival. A phylogenetic analysis on the origins of male-biased MUP gene expression in Mus musculus suggests that this sexual dimorphism evolved by increasing male MUP expression rather than reducing female expression.

Neurogenomic divergence during speciation by reinforcement of mating behaviors in chorus frogs (Pseudacris)

BACKGROUND

Species interactions can promote mating behavior divergence, particularly when these interactions are costly due to maladaptive hybridization. Selection against hybridization can indirectly cause evolution of reproductive isolation within species, a process termed cascade reinforcement. This process can drive incipient speciation by generating divergent selection pressures among populations that interact with different species assemblages. Theoretical and empirical studies indicate that divergent selection on gene expression networks has the potential to increase reproductive isolation among populations. After identifying candidate synaptic transmission genes derived from neurophysiological studies in anurans, we test for divergence of gene expression in a system undergoing cascade reinforcement, the Upland Chorus Frog (Pseudacris feriarum).

RESULTS

Our analyses identified seven candidate synaptic transmission genes that have diverged between ancestral and reinforced populations of P. feriarum, including five that encode synaptic vesicle proteins. Our gene correlation network analyses revealed four genetic modules that have diverged between these populations, two possessing a significant concentration of neurotransmission enrichment terms: one for synaptic membrane components and the other for metabolism of the neurotransmitter nitric oxide. We also ascertained that a greater number of genes have diverged in expression by geography than by sex. Moreover, we found that more genes have diverged within females as compared to males between populations. Conversely, we observed no difference in the number of differentially-expressed genes within the ancestral compared to the reinforced population between the sexes.

CONCLUSIONS

This work is consistent with the idea that divergent selection on mating behaviors via cascade reinforcement contributed to evolution of gene expression in P. feriarum. Although our study design does not allow us to fully rule out the influence of environment and demography, the fact that more genes diverged in females than males points to a role for cascade reinforcement. Our discoveries of divergent candidate genes and gene networks related to neurotransmission support the idea that neural mechanisms of acoustic mating behaviors have diverged between populations, and agree with previous neurophysiological studies in frogs. Increasing support for this hypothesis, however, will require additional experiments under common garden conditions. Our work points to the importance of future replicated and tissue-specific studies to elucidate the relative contribution of gene expression divergence to the evolution of reproductive isolation during incipient speciation.

The role of copy-number variation in the reinforcement of sexual isolation between the two European subspecies of the house mouse

Reinforcement has the potential to generate strong reproductive isolation through the evolution of barrier traits as a response to selection against maladaptive hybridization, but the genetic changes associated with this process remain largely unexplored. Building upon the increasing evidence for a role of structural variants in adaptation and speciation, we addressed the role of copy-number variation in the reinforcement of sexual isolation evidenced between the two European subspecies of the house mouse. We characterized copy-number divergence between populations of Mus musculus musculus that display assortative mate choice, and those that do not, using whole-genome resequencing data. Updating methods to detect deletions and tandem duplications (collectively: copy-number variants, CNVs) in Pool-Seq data, we developed an analytical pipeline dedicated to identifying genomic regions showing the expected pattern of copy-number displacement under a reinforcement scenario. This strategy allowed us to detect 1824 deletions and seven tandem duplications that showed extreme differences in frequency between behavioural classes across replicate comparisons. A subset of 480 deletions and four tandem duplications were specifically associated with the derived trait of assortative mate choice. These 'Choosiness-associated' CNVs occur in hundreds of genes. Consistent with our hypothesis, such genes included olfactory receptors potentially involved in the olfactory-based assortative mate choice in this system as well as one gene, Sp110, that is known to show patterns of differential expression between behavioural classes in an organ used in mate choice-the vomeronasal organ. These results demonstrate that fine-scale structural changes are common and highly variable within species, despite being under-studied, and may be important targets of reinforcing selection in this system and others. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Comparative Genomic Analysis of the Pheromone Receptor Class 1 Family (V1R) Reveals Extreme Complexity in Mouse Lemurs (Genus, Microcebus) and a Chromosomal Hotspot across Mammals

Sensory gene families are of special interest for both what they can tell us about molecular evolution and what they imply as mediators of social communication. The vomeronasal type-1 receptors (V1Rs) have often been hypothesized as playing a fundamental role in driving or maintaining species boundaries given their likely function as mediators of intraspecific mate choice, particularly in nocturnal mammals. Here, we employ a comparative genomic approach for revealing patterns of V1R evolution within primates, with a special focus on the small-bodied nocturnal mouse and dwarf lemurs of Madagascar (genera Microcebus and Cheirogaleus, respectively). By doubling the existing genomic resources for strepsirrhine primates (i.e. the lemurs and lorises), we find that the highly speciose and morphologically cryptic mouse lemurs have experienced an elaborate proliferation of V1Rs that we argue is functionally related to their capacity for rapid lineage diversification. Contrary to a previous study that found equivalent degrees of V1R diversity in diurnal and nocturnal lemurs, our study finds a strong correlation between nocturnality and V1R elaboration, with nocturnal lemurs showing elaborate V1R repertoires and diurnal lemurs showing less diverse repertoires. Recognized subfamilies among V1Rs show unique signatures of diversifying positive selection, as might be expected if they have each evolved to respond to specific stimuli. Furthermore, a detailed syntenic comparison of mouse lemurs with mouse (genus Mus) and other mammalian outgroups shows that orthologous mammalian subfamilies, predicted to be of ancient origin, tend to cluster in a densely populated region across syntenic chromosomes that we refer to as a V1R "hotspot."

The Past and Future of Experimental Speciation

Speciation is the result of evolutionary processes that generate barriers to gene flow between populations, facilitating reproductive isolation. Speciation is typically studied via theoretical models and snapshot tests in natural populations. Experimental speciation enables real-time direct tests of speciation theory and has been long touted as a critical complement to other approaches. We argue that, despite its promise to elucidate the evolution of reproductive isolation, experimental speciation has been underutilised and lags behind other contributions to speciation research. We review recent experiments and outline a framework for how experimental speciation can be implemented to address current outstanding questions that are otherwise challenging to answer. Greater uptake of this approach is necessary to rapidly advance understanding of speciation.

Genomic Signatures of Reinforcement

Reinforcement is the process by which selection against hybridization increases reproductive isolation between taxa. Much research has focused on demonstrating the existence of reinforcement, yet relatively little is known about the genetic basis of reinforcement or the evolutionary conditions under which reinforcement can occur. Inspired by reinforcement's characteristic phenotypic pattern of reproductive trait divergence in sympatry but not in allopatry, we discuss whether reinforcement also leaves a distinct genomic pattern. First, we describe three patterns of genetic variation we expect as a consequence of reinforcement. Then, we discuss a set of alternative processes and complicating factors that may make the identification of reinforcement at the genomic level difficult. Finally, we consider how genomic analyses can be leveraged to inform if and to what extent reinforcement evolved in the face of gene flow between sympatric lineages and between allopatric and sympatric populations of the same lineage. Our major goals are to understand if genome scans for particular patterns of genetic variation could identify reinforcement, isolate the genetic basis of reinforcement, or infer the conditions under which reinforcement evolved.

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.

Molecular heterogeneity in major urinary proteins of Mus musculus subspecies: potential candidates involved in speciation

When hybridisation carries a cost, natural selection is predicted to favour evolution of traits that allow assortative mating (reinforcement). Incipient speciation between the two European house mouse subspecies, Mus musculus domesticus and M.m.musculus, sharing a hybrid zone, provides an opportunity to understand evolution of assortative mating at a molecular level. Mouse urine odours allow subspecific mate discrimination, with assortative preferences evident in the hybrid zone but not in allopatry. Here we assess the potential of MUPs (major urinary proteins) as candidates for signal divergence by comparing MUP expression in urine samples from the Danish hybrid zone border (contact) and from allopatric populations. Mass spectrometric characterisation identified novel MUPs in both subspecies involving mostly new combinations of amino acid changes previously observed in M.m.domesticus. The subspecies expressed distinct MUP signatures, with most MUPs expressed by only one subspecies. Expression of at least eight MUPs showed significant subspecies divergence both in allopatry and contact zone. Another seven MUPs showed divergence in expression between the subspecies only in the contact zone, consistent with divergence by reinforcement. These proteins are candidates for the semiochemical barrier to hybridisation, providing an opportunity to characterise the nature and evolution of a putative species recognition signal.

Genetic signatures of natural selection in a model invasive ascidian

Invasive species represent promising models to study species’ responses to rapidly changing environments. Although local adaptation frequently occurs during contemporary range expansion, the associated genetic signatures at both population and genomic levels remain largely unknown. Here, we use genome-wide gene-associated microsatellites to investigate genetic signatures of natural selection in a model invasive ascidian, Ciona robusta. Population genetic analyses of 150 individuals sampled in Korea, New Zealand, South Africa and Spain showed significant genetic differentiation among populations. Based on outlier tests, we found high incidence of signatures of directional selection at 19 loci. Hitchhiking mapping analyses identified 12 directional selective sweep regions, and all selective sweep windows on chromosomes were narrow (~8.9 kb). Further analyses indentified 132 candidate genes under selection. When we compared our genetic data and six crucial environmental variables, 16 putatively selected loci showed significant correlation with these environmental variables. This suggests that the local environmental conditions have left significant signatures of selection at both population and genomic levels. Finally, we identified “plastic” genomic regions and genes that are promising regions to investigate evolutionary responses to rapid environmental change in C. robusta.

Hybridization in Plants: Old Ideas, New Techniques

Hybridization has played an important role in the evolution of many lineages. With the growing availability of genomic tools and advancements in genomic analyses, it is becoming increasingly clear that gene flow between divergent taxa can generate new phenotypic diversity, allow for adaptation to novel environments, and contribute to speciation. Hybridization can have immediate phenotypic consequences through the expression of hybrid vigor. On longer evolutionary time scales, hybridization can lead to local adaption through the introgression of novel alleles and transgressive segregation and, in some cases, result in the formation of new hybrid species. Studying both the abundance and the evolutionary consequences of hybridization has deep historical roots in plant biology. Many of the hypotheses concerning how and why hybridization contributes to biological diversity currently being investigated were first proposed tens and even hundreds of years ago. In this Update, we discuss how new advancements in genomic and genetic tools are revolutionizing our ability to document the occurrence of and investigate the outcomes of hybridization in plants.

Diversity of major urinary proteins (MUPs) in wild house mice

Major urinary proteins (MUPs) are often suggested to be highly polymorphic, and thereby provide unique chemical signatures used for individual and genetic kin recognition; however, studies on MUP variability have been lacking. We surveyed populations of wild house mice (Mus musculus musculus), and examined variation of MUP genes and proteins. We sequenced several Mup genes (9 to 11 loci) and unexpectedly found no inter-individual variation. We also found that microsatellite markers inside the MUP cluster show remarkably low levels of allelic diversity, and significantly lower than the diversity of markers flanking the cluster or other markers in the genome. We found low individual variation in the number and types of MUP proteins using a shotgun proteomic approach, even among mice with variable MUP electrophoretic profiles. We identified gel bands and spots using high-resolution mass spectrometry and discovered that gel-based methods do not separate MUP proteins, and therefore do not provide measures of MUP diversity, as generally assumed. The low diversity and high homology of Mup genes are likely maintained by purifying selection and gene conversion, and our results indicate that the type of selection on MUPs and their adaptive functions need to be re-evaluated.