Genome-wide patterns of gene flow across a house mouse hybrid zone

Genomic networks of hybrid sterility

Hybrid dysfunction, a common feature of reproductive barriers between species, is often caused by negative epistasis between loci ("Dobzhansky-Muller incompatibilities"). The nature and complexity of hybrid incompatibilities remain poorly understood because identifying interacting loci that affect complex phenotypes is difficult. With subspecies in the early stages of speciation, an array of genetic tools, and detailed knowledge of reproductive biology, house mice (Mus musculus) provide a model system for dissecting hybrid incompatibilities. Male hybrids between M. musculus subspecies often show reduced fertility. Previous studies identified loci and several X chromosome-autosome interactions that contribute to sterility. To characterize the genetic basis of hybrid sterility in detail, we used a systems genetics approach, integrating mapping of gene expression traits with sterility phenotypes and QTL. We measured genome-wide testis expression in 305 male F2s from a cross between wild-derived inbred strains of M. musculus musculus and M. m. domesticus. We identified several thousand cis- and trans-acting QTL contributing to expression variation (eQTL). Many trans eQTL cluster into eleven 'hotspots,' seven of which co-localize with QTL for sterility phenotypes identified in the cross. The number and clustering of trans eQTL-but not cis eQTL-were substantially lower when mapping was restricted to a 'fertile' subset of mice, providing evidence that trans eQTL hotspots are related to sterility. Functional annotation of transcripts with eQTL provides insights into the biological processes disrupted by sterility loci and guides prioritization of candidate genes. Using a conditional mapping approach, we identified eQTL dependent on interactions between loci, revealing a complex system of epistasis. Our results illuminate established patterns, including the role of the X chromosome in hybrid sterility. The integrated mapping approach we employed is applicable in a broad range of organisms and we advocate for widespread adoption of a network-centered approach in speciation genetics.

X chromosome control of meiotic chromosome synapsis in mouse inter-subspecific hybrids

Hybrid sterility (HS) belongs to reproductive isolation barriers that safeguard the integrity of species in statu nascendi. Although hybrid sterility occurs almost universally among animal and plant species, most of our current knowledge comes from the classical genetic studies on Drosophila interspecific crosses or introgressions. With the house mouse subspecies Mus m. musculus and Mus m. domesticus as a model, new research tools have become available for studies of the molecular mechanisms and genetic networks underlying HS. Here we used QTL analysis and intersubspecific chromosome substitution strains to identify a 4.7 Mb critical region on Chromosome X (Chr X) harboring the Hstx2 HS locus, which causes asymmetrical spermatogenic arrest in reciprocal intersubspecific F1 hybrids. Subsequently, we mapped autosomal loci on Chrs 3, 9 and 13 that can abolish this asymmetry. Combination of immunofluorescent visualization of the proteins of synaptonemal complexes with whole-chromosome DNA FISH on pachytene spreads revealed that heterosubspecific, unlike consubspecific, homologous chromosomes are predisposed to asynapsis in F1 hybrid male and female meiosis. The asynapsis is under the trans- control of Hstx2 and Hst1/Prdm9 hybrid sterility genes in pachynemas of male but not female hybrids. The finding concurred with the fertility of intersubpecific F1 hybrid females homozygous for the Hstx2^Mmm allele and resolved the apparent conflict with the dominance theory of Haldane's rule. We propose that meiotic asynapsis in intersubspecific hybrids is a consequence of cis-acting mismatch between homologous chromosomes modulated by the trans-acting Hstx2 and Prdm9 hybrid male sterility genes.

Genomes of newly emerging species restrict their gene exchange with related taxa in order to secure integrity. Hybrid sterility is one of the reproductive isolation mechanisms restricting gene flow between closely related, sexually reproducing organisms. We showed that hybrid sterility between two closely related mouse subspecies is executed by a failure of meiotic synapsis of orthologous chromosomes in F1 hybrid males. The asynapsis of orthologous chromosomes occurred in meiosis of male and female hybrids, though only males were sterile due to trans-acting male-specific hybrid sterility genes. We located one of the two major hybrid sterility genes to a 4.7 Mb interval on Chromosome X, showed that it controls male sterility by modulating the extent of meiotic asynapsis and using the inter-subspecific chromosome substitution strains we refuted the simple interpretation of dominance theory of Haldane's rule. A new working hypothesis posits male sterility of mouse inter-subsubspecific F1 hybrids as a consequence of meiotic chromosome asynapsis caused by the cis-acting mismatch between orthologous chromosomes modulated by the trans-acting hybrid male sterility genes.

Faster-X adaptive protein evolution in house mice

The causes of the large effect of the X chromosome in reproductive isolation and speciation have long been debated. The faster-X hypothesis predicts that X-linked loci are expected to have higher rates of adaptive evolution than autosomal loci if new beneficial mutations are on average recessive. Reproductive isolation should therefore evolve faster when contributing loci are located on the X chromosome. In this study, we have analyzed genome-wide nucleotide polymorphism data from the house mouse subspecies Mus musculus castaneus and nucleotide divergence from Mus famulus and Rattus norvegicus to compare rates of adaptive evolution for autosomal and X-linked protein-coding genes. We found significantly faster adaptive evolution for X-linked loci, particularly for genes with expression in male-specific tissues, but autosomal and X-linked genes with expression in female-specific tissues evolve at similar rates. We also estimated rates of adaptive evolution for genes expressed during spermatogenesis and found that X-linked genes that escape meiotic sex chromosome inactivation (MSCI) show rapid adaptive evolution. Our results suggest that faster-X adaptive evolution is either due to net recessivity of new advantageous mutations or due to a special gene content of the X chromosome, which regulates male function and spermatogenesis. We discuss how our results help to explain the large effect of the X chromosome in speciation.

Evolution of zygotic linkage disequilibrium in a finite local population

One crucial feature of zygotic linkage disequilibrium (LD) analysis is its direct use of diploid genotyping data, irrespective of the type of mating system. Previous theories from an evolutionary perspective mainly focus on gametic LD, but the equivalent development for zygotic LD is not available. Here I study the evolution of zygotic LD and the covariances between gametic and zygotic LDs or between distinct zygotic LDs in a finite local population under constant immigration from a continent population. I derive the analytical theory under genetic hitchhiking effects or in a neutral process. Results indicate that zygotic LDs (diploid level) are more informative than gametic LD (haploid level) in indicating the effects of different evolutionary forces. Zygotic LDs may be greater than or comparable to gametic LD under the epistatic selection process, but smaller than gametic LD under the non epistatic selection process. The covariances between gametic and zygotic LDs are strongly affected by the mating system, linkage distance, and genetic drift effects, but weakly affected by seed and pollen flow and natural selection. The covariances between different zygotic LDs are generally robust to the effects of gene flow, selection, and linkage distance, but sensitive to the effects of genetic drift and mating system. Consistent patterns exist for the covariances between the zygotic LDs for the two-locus genotypes with one common genotype at one locus or without any common genotype at each locus. The results highlight that zygotic LDs can be applied to detecting natural population history.

Introgression in the Drosophila subobscura--D. Madeirensis sister species: evidence of gene flow in nuclear genes despite mitochondrial differentiation

Species hybridization, and thus the potential for gene flow, was once viewed as reproductive mistake. However, recent analysis based on large datasets and newly developed models suggest that gene exchange is not as rare as originally suspected. To investigate the history and speciation of the closely related species Drosophila subobscura, D. madeirensis, and D. guanche, we obtained polymorphism and divergence data for 26 regions throughout the genome, including the Y chromosome and mitochondrial DNA. We found that the D. subobscura X/autosome ratio of silent nucleotide diversity is significantly smaller than the 0.75 expected under neutrality. This pattern, if held genomewide, may reflect a faster accumulation of beneficial mutations on the X chromosome than on autosomes. We also detected evidence of gene flow in autosomal regions, while sex chromosomes remain distinct. This is consistent with the large X effect on hybrid male sterility seen in this system and the presence of two X chromosome inversions fixed between species. Overall, our data conform to chromosomal speciation models in which rearrangements are proposed to serve as gene flow barriers. Contrary to other observations in Drosophila, the mitochondrial genome appears resilient to gene flow in the presence of nuclear exchange.

Differential introgression in a mosaic hybrid zone reveals candidate barrier genes

Hybrid zones act as genomic sieves. Although globally advantageous alleles will spread throughout the zone and neutral alleles can be freely exchanged between species, introgression will be restricted for genes that contribute to reproductive barriers or local adaptation. Seminal fluid proteins (SFPs) are known to contribute to reproductive barriers in insects and have been proposed as candidate barrier genes in the hybridizing field crickets Gryllus pennsylvanicus and Gryllus firmus. Here, we have used 125 single nucleotide polymorphisms to characterize patterns of differential introgression and to identify genes that may contribute to prezygotic barriers between these species. Using a transcriptome scan of the male cricket accessory gland (the site of SFP synthesis), we identified genes with major allele frequency differences between the species. We then compared patterns of introgression for genes encoding SFPs with patterns for genes expressed in the same tissue that do not encode SFPs. We find no evidence that SFPs have reduced gene exchange across the cricket hybrid zone. However, a number of genes exhibit dramatically reduced introgression, and many of these genes encode proteins with functional roles consistent with known barriers.

Mate choice in Mus musculus is relative and dependent on the estrous state

Mate choice is a critical behavioral decision process with profound impact on evolution. However, the mechanistic basis of mate choice is poorly understood. In this study we focused on assortative mate choice, which is known to contribute to the reproductive isolation of the two European subspecies of house mouse, Mus musculus musculus and Mus musculus domesticus. To understand the decision process, we developed both full mating and limited-contact paradigms and tested musculus females' preference for musculus versus domesticus males, mimicking the natural musculus/domesticus contact zone. As hypothesized, when allowed to mate we found that sexually receptive musculus females exhibited a robust preference to mate with musculus males. In contrast, when non-receptive, females did not exhibit a preference and rather alternated between males in response to male mount attempts. Moreover in a no-choice condition, females mated readily with males from both subspecies. Finally, when no physical contact was allowed, and therefore male's behavior could not influence female's behavior, female's preference for its own subspecies was maintained independently of the estrous state. Together, our results suggest that the assortative preference is relative and based on a comparison of the options available rather than on an absolute preference. The results of the limited-contact experiments highlight the interplay between female's internal state and the nature of the interaction with prospective mates in the full mating conditions. With these experiments we believe we established an assortative mate preference assay that is appropriate for the investigation of its underlying substrates.

Evidence for asymmetrical divergence-gene flow of nuclear loci, but not mitochondrial loci, between seabird sister species: blue-footed (Sula nebouxii) and Peruvian (S. variegata) boobies

Understanding the process of speciation requires understanding how gene flow influences divergence. Recent analyses indicate that divergence can take place despite gene flow and that the sex chromosomes can exhibit different levels of gene flow than autosomes and mitochondrial DNA. Using an eight marker dataset including autosomal, z-linked, and mitochondrial loci we tested the hypothesis that blue-footed (Sula nebouxii) and Peruvian (S. variegata) boobies diverged from their common ancestor with gene flow, paying specific attention to the differences in gene flow estimates from nuclear and mitochondrial markers. We found no gene flow at mitochondrial markers, but found evidence from the combined autosomal and z-linked dataset that blue-footed and Peruvian boobies experienced asymmetrical gene flow during or after their initial divergence, predominantly from Peruvian boobies into blue-footed boobies. This gene exchange may have occurred either sporadically between periods of allopatry, or regularly throughout the divergence process. Our results add to growing evidence that diverging species can remain distinct but exchange genes.

Identifying loci under selection against gene flow in isolation-with-migration models

When divergence occurs in the presence of gene flow, there can arise an interesting dynamic in which selection against gene flow, at sites associated with population-specific adaptations or genetic incompatibilities, can cause net gene flow to vary across the genome. Loci linked to sites under selection may experience reduced gene flow and may experience genetic bottlenecks by the action of nearby selective sweeps. Data from histories such as these may be poorly fitted by conventional neutral model approaches to demographic inference, which treat all loci as equally subject to forces of genetic drift and gene flow. To allow for demographic inference in the face of such histories, as well as the identification of loci affected by selection, we developed an isolation-with-migration model that explicitly provides for variation among genomic regions in migration rates and/or rates of genetic drift. The method allows for loci to fall into any of multiple groups, each characterized by a different set of parameters, thus relaxing the assumption that all loci share the same demography. By grouping loci, the method can be applied to data with multiple loci and still have tractable dimensionality and statistical power. We studied the performance of the method using simulated data, and we applied the method to study the divergence of two subspecies of European rabbits (Oryctolagus cuniculus).

Detecting selection along environmental gradients: analysis of eight methods and their effectiveness for outbreeding and selfing populations

Thanks to genome-scale diversity data, present-day studies can provide a detailed view of how natural and cultivated species adapt to their environment and particularly to environmental gradients. However, due to their sensitivity, up-to-date studies might be more sensitive to undocumented demographic effects such as the pattern of migration and the reproduction regime. In this study, we provide guidelines for the use of popular or recently developed statistical methods to detect footprints of selection. We simulated 100 populations along a selective gradient and explored different migration models, sampling schemes and rates of self-fertilization. We investigated the power and robustness of eight methods to detect loci potentially under selection: three designed to detect genotype-environment correlations and five designed to detect adaptive differentiation (based on F(ST) or similar measures). We show that genotype-environment correlation methods have substantially more power to detect selection than differentiation-based methods but that they generally suffer from high rates of false positives. This effect is exacerbated whenever allele frequencies are correlated, either between populations or within populations. Our results suggest that, when the underlying genetic structure of the data is unknown, a number of robust methods are preferable. Moreover, in the simulated scenario we used, sampling many populations led to better results than sampling many individuals per population. Finally, care should be taken when using methods to identify genotype-environment correlations without correcting for allele frequency autocorrelation because of the risk of spurious signals due to allele frequency correlations between populations.

Sperm-related phenotypes implicated in both maintenance and breakdown of a natural species barrier in the house mouse

The house mouse hybrid zone (HMHZ) is a species barrier thought to be maintained by a balance between dispersal and natural selection against hybrids. While the HMHZ is characterized by frequency discontinuities for some sex chromosome markers, there is an unexpected large-scale regional introgression of a Y chromosome across the barrier, in defiance of Haldane's rule. Recent work suggests that a major force maintaining the species barrier acts through sperm traits. Here, we test whether the Y chromosome penetration of the species barrier acts through sperm traits by assessing sperm characteristics of wild-caught males directly in a field laboratory set up in a Y introgression region of the HMHZ, later calculating the hybrid index of each male using 1401 diagnostic single nucleotide polymorphisms (SNPs). We found that both sperm count (SC) and sperm velocity were significantly reduced across the natural spectrum of hybrids. However, SC was more than rescued in the presence of the invading Y. Our results imply an asymmetric advantage for Y chromosome introgression consistent with the observed large-scale introgression. We suggest that selection on sperm-related traits probably explains a large component of patterns observed in the natural hybrid zone, including the Y chromosome penetration.

Nuclear gene variation in wild brown rats

Although the brown rat (Rattus norvegicus) is widely used as a model mammal throughout biological sciences, little is known about genetic variation in wild rat populations or the relationship of commonly used inbred strains to their wild relatives. We sampled wild brown rats from the species’ presumed ancestral range in NW China and from a derived population in the UK and estimated nucleotide diversity and population subdivision, based on the sequences of 30 autosomal protein-coding loci. Neutral genetic diversity was close to 0.2% in both populations, which is about five times lower than diversity at the orthologous sites in a population of wild house mice from the species’ putative ancestral range in India. We found significant population differentiation between UK and Chinese populations, as assessed by F_st and the program STRUCTURE. Based on synonymous diversity and divergence between the brown rat and house mouse, we estimate that the recent effective population size in brown rats is approximately 130,000 (approximate 95% confidence interval 85,000-184,000), about fivefold lower than wild house mice.

Patterns of transcriptome divergence in the male accessory gland of two closely related species of field crickets

One of the central questions in evolutionary genetics is how much of the genome is involved in the early stages of divergence between populations, causing them to be reproductively isolated. In this article, we investigate genomic differentiation in a pair of closely related field crickets (Gryllus firmus and G. pennsylvanicus). These two species are the result of allopatric divergence and now interact along an extensive hybrid zone in eastern North America. Genes encoding seminal fluid proteins (SFPs) are often divergent between species, and it has been hypothesized that these proteins may play a key role in the origin and maintenance of reproductive isolation between diverging lineages. Hence, we chose to scan the accessory gland transcriptome to enable direct comparisons of differentiation for genes known to encode SFPs with differentiation in a much larger set of genes expressed in the same tissue. We have characterized differences in allele frequency between two populations for >6000 SNPs and >26,000 contigs. About 10% of all SNPs showed nearly fixed differences between the two species. Genes encoding SFPs did not have significantly elevated numbers of fixed SNPs per contig, nor did they seem to show larger differences than expected in their average allele frequencies. The distribution of allele frequency differences across the transcriptome is distinctly bimodal, but the relatively high proportion of fixed SNPs does not necessarily imply "ancient" divergence between these two lineages. Further studies of linkage disequilibrium and introgression across the hybrid zone are needed to direct our attention to those genome regions that are important for reproductive isolation.

The telocentric tandem repeat at the p-arm is not conserved in Mus musculus subspecies

Mouse chromosomes, with the exception of the Y chromosome, are telocentric. The telomere at the p-arm is separated from the centromere by the tL1 sequence and TLC tandem repeats. A previous report showed that the TLC array was also conserved in other strains of the subgenus Mus. These results suggest that the TLC arrays promote the stable evolutionary maintenance of a telocentric karyotype in the subgenus Mus. In this study, we investigated the degree of conservation of TLC arrays among a variety of wild-derived inbred strains, all of which are descendants of wild mice captured in several areas of the world. Genomic PCR analysis indicates that the sequential order of telomere-tL1 is highly conserved in all strains, whereas tL1-TLC is not. Next, Southern blot analysis of DNAs isolated from a panel of mouse subspecies showed both Mus musculus domesticus and Mus musculus castaneus subspecies possess TLC arrays. Unexpectedly, this repeat appears to be lost in almost all Mus musculus musculus and Mus musculus molossinus subspecies, which show a clear geographic divide. These results indicate that either other unknown sequences were replaced by the TLC repeat or almost all M. m. musculus and M. m. molossinus subspecies do not have any sequence between the telomere and minor satellites. Our observation suggests that the TLC array might be evolutionarily unstable and not essential for murine chromosomal conformation. This is the first example of the subspecies-specific large genome alterations in mice.

Genomic variation in cline shape across a hybrid zone

Hybrid zones are unique biological interfaces that reveal both population level and species level evolutionary processes. A genome-scale approach to assess gene flow across hybrid zones is vital, and now possible. In Mexican towhees (genus Pipilo), several morphological hybrid gradients exist. We completed a genome survey across one such gradient (9 populations, 140 birds) using mitochondrial DNA, 28 isozyme, and 377 AFLP markers. To assess variation in introgression among loci, cline parameters (i.e., width, center) for the 61 clinally varying loci were estimated and compiled into genomic distributions for tests against three empirical models spanning the range of observed cline shape. No single model accounts for observed variation in cline shape among loci. Numerous backcross individuals near the gradient center confirm a hybrid origin for these populations, contrary to a previous hypothesis based on social mimicry and character displacement. In addition, the observed variation does not bin into well-defined categories of locus types (e.g., neutral vs. highly selected). Our multi-locus analysis reveals cross-genomic variation in selective constraints on gene flow and locus-specific flexibility in the permeability of the interspecies membrane.

The language of speciation

The literature on speciation has expanded dramatically in recent years, catalyzed by the emergence of new conceptual frameworks, new theoretical approaches, and new methods for characterizing pattern and inferring process. As a consequence, the language used to describe the speciation process has become more complex. Increasing complexity may be an accurate reflection of current thinking with respect to how phenotypic differences limit gene flow, how selection results in the evolution of reproductive isolation, and genetic changes that contribute to speciation. However, increased language complexity has come at a cost; old definitions have been reconfigured and new terms have been introduced. In some instances, the introduction of new terminology has failed to recognize historical usage, leading to unnecessary ambiguity and redundancy. Although the writings of Mayr and Dobzhansky remain a reference point in the language of speciation, the last decades of the 20th century saw substantial changes in our thinking about the speciation process. During that period, the language of speciation remained relatively stable. In contrast, the first decade of the 21st century has witnessed a remarkable expansion of the language of speciation. Here, the origin and evolution of ideas about speciation are viewed through the lens of changing language use.

Genome patterns of selection and introgression of haplotypes in natural populations of the house mouse (Mus musculus)

General parameters of selection, such as the frequency and strength of positive selection in natural populations or the role of introgression, are still insufficiently understood. The house mouse (Mus musculus) is a particularly well-suited model system to approach such questions, since it has a defined history of splits into subspecies and populations and since extensive genome information is available. We have used high-density single-nucleotide polymorphism (SNP) typing arrays to assess genomic patterns of positive selection and introgression of alleles in two natural populations of each of the subspecies M. m. domesticus and M. m. musculus. Applying different statistical procedures, we find a large number of regions subject to apparent selective sweeps, indicating frequent positive selection on rare alleles or novel mutations. Genes in the regions include well-studied imprinted loci (e.g. Plagl1/Zac1), homologues of human genes involved in adaptations (e.g. alpha-amylase genes) or in genetic diseases (e.g. Huntingtin and Parkin). Haplotype matching between the two subspecies reveals a large number of haplotypes that show patterns of introgression from specific populations of the respective other subspecies, with at least 10% of the genome being affected by partial or full introgression. Using neutral simulations for comparison, we find that the size and the fraction of introgressed haplotypes are not compatible with a pure migration or incomplete lineage sorting model. Hence, it appears that introgressed haplotypes can rise in frequency due to positive selection and thus can contribute to the adaptive genomic landscape of natural populations. Our data support the notion that natural genomes are subject to complex adaptive processes, including the introgression of haplotypes from other differentiated populations or species at a larger scale than previously assumed for animals. This implies that some of the admixture found in inbred strains of mice may also have a natural origin.

Genetics and evolution of hybrid male sterility in house mice

Comparative genetic mapping provides insights into the evolution of the reproductive barriers that separate closely related species. This approach has been used to document the accumulation of reproductive incompatibilities over time, but has only been applied to a few taxa. House mice offer a powerful system to reconstruct the evolution of reproductive isolation between multiple subspecies pairs. However, studies of the primary reproductive barrier in house mice-hybrid male sterility-have been restricted to a single subspecies pair: Mus musculus musculus and Mus musculus domesticus. To provide a more complete characterization of reproductive isolation in house mice, we conducted an F(2) intercross between wild-derived inbred strains from Mus musculus castaneus and M. m. domesticus. We identified autosomal and X-linked QTL associated with a range of hybrid male sterility phenotypes, including testis weight, sperm density, and sperm morphology. The pseudoautosomal region (PAR) was strongly associated with hybrid sterility phenotypes when heterozygous. We compared QTL found in this cross with QTL identified in a previous F(2) intercross between M. m. musculus and M. m. domesticus and found three shared autosomal QTL. Most QTL were not shared, demonstrating that the genetic basis of hybrid male sterility largely differs between these closely related subspecies pairs. These results lay the groundwork for identifying genes responsible for the early stages of speciation in house mice.

Genome-wide architecture of reproductive isolation in a naturally occurring hybrid zone between Mus musculus musculus and M. m. domesticus

Studies of a hybrid zone between two house mouse subspecies (Mus musculus musculus and M. m. domesticus) along with studies using laboratory crosses reveal a large role for the X chromosome and multiple autosomal regions in reproductive isolation as a consequence of disrupted epistasis in hybrids. One limitation of previous work has been that most of the identified genomic regions have been large. The goal here is to detect and characterize precise genomic regions underlying reproductive isolation. We surveyed 1401 markers evenly spaced across the genome in 679 mice collected from two different transects. Comparisons between transects provide a means for identifying common patterns that likely reflect intrinsic incompatibilities. We used a genomic cline approach to identify patterns that correspond to epistasis. From both transects, we identified contiguous regions on the X chromosome in which markers were inferred to be involved in epistatic interactions. We then searched for autosomal regions showing the same patterns and found they constitute about 5% of autosomal markers. We discovered substantial overlap between these candidate regions underlying reproductive isolation and QTL for hybrid sterility identified in laboratory crosses. Analysis of gene content in these regions suggests a key role for several mechanisms, including the regulation of transcription, sexual conflict and sexual selection operating at both the postmating prezygotic and postzygotic stages of reproductive isolation. Taken together, these results indicate that speciation in two recently diverged (c. 0.5 Ma) house mouse subspecies is complex, involving many genes dispersed throughout the genome and associated with distinct functions.

Dissecting the genetic architecture of F1 hybrid sterility in house mice

Hybrid sterility as a postzygotic reproductive isolation mechanism has been studied for over 80 years, yet the first identifications of hybrid sterility genes in Drosophila and mouse are quite recent. To study the genetic architecture of F(1) hybrid sterility between young subspecies of house mouse Mus m. domesticus and M. m. musculus, we conducted QTL analysis of a backcross between inbred strains representing these two subspecies and probed the role of individual chromosomes in hybrid sterility using the intersubspecific chromosome substitution strains. We provide direct evidence that the asymmetry in male infertility between reciprocal crosses is conferred by the middle region of M. m. musculus Chr X, thus excluding other potential candidates such as Y, imprinted genes, and mitochondrial DNA. QTL analysis identified strong hybrid sterility loci on Chr 17 and Chr X and predicted a set of interchangeable autosomal loci, a subset of which is sufficient to activate the Dobzhansky-Muller incompatibility of the strong loci. Overall, our results indicate the oligogenic nature of F(1) hybrid sterility, which should be amenable to reconstruction by proper combination of chromosome substitution strains. Such a prefabricated model system should help to uncover the gene networks and molecular mechanisms underlying hybrid sterility.

Genomics of isolation in hybrids

Hybrid zones are common in nature and can offer critical insights into the dynamics and components of reproductive isolation. Hybrids between diverged lineages are particularly informative about the genetic architecture of reproductive isolation, because introgression in an admixed population is a direct measure of isolation. In this paper, we combine simulations and a new statistical model to determine the extent to which different genetic architectures of isolation leave different signatures on genome-level patterns of introgression. We found that reproductive isolation caused by one or several loci of large effect caused greater heterogeneity in patterns of introgression than architectures involving many loci with small fitness effects, particularly when isolating factors were closely linked. The same conditions that led to heterogeneous introgression often resulted in a reasonable correspondence between outlier loci and the genetic loci that contributed to isolation. However, demographic conditions affected both of these results, highlighting potential limitations to the study of the speciation genomics. Further progress in understanding the genomics of speciation will require large-scale empirical studies of introgression in hybrid zones and model-based analyses, as well as more comprehensive modelling of the expected levels of isolation with different demographies and genetic architectures of isolation.

Recombination rate variation and speciation: theoretical predictions and empirical results from rabbits and mice

Recently diverged taxa may continue to exchange genes. A number of models of speciation with gene flow propose that the frequency of gene exchange will be lower in genomic regions of low recombination and that these regions will therefore be more differentiated. However, several population-genetic models that focus on selection at linked sites also predict greater differentiation in regions of low recombination simply as a result of faster sorting of ancestral alleles even in the absence of gene flow. Moreover, identifying the actual amount of gene flow from patterns of genetic variation is tricky, because both ancestral polymorphism and migration lead to shared variation between recently diverged taxa. New analytic methods have been developed to help distinguish ancestral polymorphism from migration. Along with a growing number of datasets of multi-locus DNA sequence variation, these methods have spawned a renewed interest in speciation models with gene flow. Here, we review both speciation and population-genetic models that make explicit predictions about how the rate of recombination influences patterns of genetic variation within and between species. We then compare those predictions with empirical data of DNA sequence variation in rabbits and mice. We find strong support for the prediction that genomic regions experiencing low levels of recombination are more differentiated. In most cases, reduced gene flow appears to contribute to the pattern, although disentangling the relative contribution of reduced gene flow and selection at linked sites remains a challenge. We suggest fruitful areas of research that might help distinguish between different models.

No evidence for biased co-transmission of speciation islands in Anopheles gambiae

Genome-scale scans have revealed highly heterogeneous levels of divergence between closely related taxa in many systems. Generally, a small number of regions show high differentiation, with the rest of the genome showing no or only low levels of divergence. These patterns have been interpreted as evidence for ongoing speciation-with-gene-flow, with introgression homogenizing the whole genome except loci involved in reproductive isolation. However, as the number of selected loci increases, the probability of introgression at unselected loci decreases unless there is a transmission ratio distortion causing an over-representation of specific combinations of alleles. Here we examine the transmission of three 'speciation islands' that contain fixed differences between the M and S forms of the mosquito, Anopheles gambiae. We made reciprocal crosses between M and S parents and genotyped over 2000 F(2) individuals, developing a hierarchical likelihood model to identify specific genotypes that are under- or over-represented among the recombinant offspring. Though our overall results did not match the expected number of F(2) genotypes, we found no biased co-transmission among M or S alleles in the three islands. Our likelihood model did identify transmission ratio distortion at two of the three islands, but this distortion was small (approx. 3%) and in opposite directions for the two islands. We discuss how our results impinge on hypotheses of current gene flow between M and S and ongoing speciation-with-gene-flow in this system.

The strange blood: natural hybridization in primates

Hybridization between two closely related species is a natural evolutionary process that results in an admixture of previously isolated gene pools. The exchange of genes between species may accelerate adaptation and lead to the formation of new lineages. Hybridization can be regarded as one important evolutionary mechanism driving speciation processes. Although recent studies have highlighted the taxonomic breadth of natural hybridization in the primate order, information about primate hybridization is still limited compared to that about the hybridization of fish, birds, or other mammals. In primates, hybridization has occurred mainly between subspecies and species, but has also been detected between genera and even in the human lineage. Here we provide an overview of cases of natural hybridization in all major primate radiations. Our review emphasizes a phylogenetic approach. We use the data presented to discuss the impact of hybridization on taxonomy and conservation.

Strong selection against hybrids maintains a narrow contact zone between morphologically cryptic lineages in a rainforest lizard

Phenotypically cryptic lineages comprise an important yet understudied part of biodiversity; in particular, we have much to learn about how these lineages are formed and maintained. To better understand the evolutionary significance of such lineages, we studied a hybrid zone between two morphologically cryptic phylogeographic lineages in the rainforest lizard, Lampropholis coggeri. Analyzing a multilocus genetic dataset through cline inference, individual-based methods and population measures of disequilibrium and using simulations to explore our genetic results in context of theoretical expectations, we inferred the processes maintaining this hybrid zone. We find that these lineages meet in a hybrid zone that is narrow (≈400 m) relative to inferred dispersal rate. Further, the hybrid zone exhibits substantial genetic disequilibrium and sharply coincident and largely concordant clines. Based on our knowledge about the region's biogeography, the species' natural history, and our simulation results, we suggest that strong selection against hybrids structures this system. As all clines show a relatively narrow range of introgression, we posit that this hybrid zone might not yet be in equilibrium. Nonetheless, our results clearly show that phylogeographic lineages can evolve substantial reproductive isolation without concomitant morphological diversification, suggesting that such lineages can constitute a significant component of evolutionary diversity.

Introgression of mitochondrial DNA among Myodes voles: consequences for energetics?

BACKGROUND

Introgression of mitochondrial DNA (mtDNA) is among the most frequently described cases of reticulate evolution. The tendency of mtDNA to cross interspecific barriers is somewhat counter-intuitive considering the key function of enzymes that it encodes in the oxidative-phosphorylation process, which could give rise to hybrid dysfunction. How mtDNA reticulation affects the evolution of metabolic functions is, however, uncertain. Here we investigated how morpho-physiological traits vary in natural populations of a common rodent (the bank vole, Myodes glareolus) and whether this variation could be associated with mtDNA introgression. First, we confirmed that M. glareolus harbour mtDNA introgressed from M. rutilus by analyzing mtDNA (cytochrome b, 954 bp) and nuclear DNA (four markers; 2333 bp in total) sequence variation and reconstructing loci phylogenies among six natural populations in Finland. We then studied geographic variation in body size and basal metabolic rate (BMR) among the populations of M. glareolus and tested its relationship with mtDNA type.

RESULTS

Myodes glareolus and its arctic neighbour, M. rutilus, are reciprocally monophyletic at the analyzed nuclear DNA loci. In contrast, the two northernmost populations of M. glareolus have a fixed mitotype that is shared with M. rutilus, likely due to introgressive hybridization. The analyses of phenotypic traits revealed that the body mass and whole-body, but not mass corrected, BMR are significantly reduced in M. glareolus females from northern Finland that also have the introgressed mitotype. Restricting the analysis to the single population where the mitotypes coexist, the association of mtDNA type with whole-body BMR remained but those with mass corrected BMR and body mass did not. Mitochondrial sequence variation in the introgressed haplotypes is compatible with demographic growth of the populations, but may also be a result of positive selection.

CONCLUSION

Our results show that the phenotypic traits vary markedly along the north-south axis of populations of M. glareolus. This variation may be related to adaptation to local environments and coincides with the gradient of genome reticulation between M. glareolus and M. rutilus, which was assessed by mtDNA introgression. Introgression of mtDNA may have affected morpho-physiological traits but do not show strong effects on either body mass or basal metabolic rate alone. We discuss the causes and biological meaning of our results and the means to clarify these questions in future research.

The phenotypic distribution of quantitative traits in a wild mouse F1 population

The human complex diseases such as hypertension, precocious puberty, and diabetes have their own diagnostic thresholds, which are usually estimated from the epidemiological data of nature populations. In the mouse models, numerous phenotypic data of complex traits have been accumulated; however, knowledge of the phenotypic distribution of the natural mouse populations remains quite limited. In order to investigate the distribution of quantitative traits of wild mice, 170 F1 progeny aged 8-10 weeks and derived from wild mice collected from eight spots in the suburbs of Shanghai were tested for their values of anatomic, blood chemical, and blood hematological parameters. All the wild mice breeders were of Mus. m. musculus and Mus. m. castaneus maternal origin according to the single nucleotide polymorphism (SNP) markers of the mitochondrial DNA. The results showed that phenotypes in wild mice had a normal distribution with four to six times the standard deviation. For the majority of the traits, the wild outbred mice and laboratory inbred mice have significantly different ranges and mean values, whereas the wild mice did not necessarily show more phenotypic diversity than the inbred ones. Our data also showed that natural populations may have some unique phenotypes related to sugar and protein metabolism, as the mean value of wild mice differ dramatically from the inbred mice in the levels of blood glucose, BUN (blood urea nitrogen), and total blood protein. The epidemiological information of the complex traits in the nature population from our study provided valuable reference for the application of mouse models in those complex disease studies.

Higher differentiation among subspecies of the house mouse (Mus musculus) in genomic regions with low recombination

In the early stages of reproductive isolation, genomic regions of reduced recombination are expected to show greater levels of differentiation, either because gene flow between species is reduced in these regions or because the effects of selection at linked sites within species are enhanced in these regions. Here, we study the patterns of DNA sequence variation at 27 autosomal loci among populations of Mus musculus musculus, M. m. domesticus, and M. m. castaneus, three subspecies of house mice with collinear genomes. We found that some loci exhibit considerable shared variation among subspecies, while others exhibit fixed differences. We used an isolation-with-gene-flow model to estimate divergence times and effective population sizes (N(e) ) and to disentangle ancestral variation from gene flow. Estimates of divergence time indicate that all three subspecies diverged from one another within a very short period of time approximately 350,000 years ago. Overall, N(e) for each subspecies was associated with the degree of genetic differentiation: M. m. musculus had the smallest N(e) and the greatest proportion of monophyletic gene genealogies, while M. m. castaneus had the largest N(e) and the smallest proportion of monophyletic gene genealogies. M. m. domesticus and M. m. musculus were more differentiated from each other than either were from M. m. castaneus, consistent with greater reproductive isolation between M. m. domesticus and M. m. musculus. F(ST) was significantly greater at loci experiencing low recombination rates compared to loci experiencing high recombination rates in comparisons between M. m. castaneus and M. m. musculus or M. m. domesticus. These results provide evidence that genomic regions with less recombination show greater differentiation, even in the absence of chromosomal rearrangements.

Reticulate evolution: frequent introgressive hybridization among Chinese hares (genus lepus) revealed by analyses of multiple mitochondrial and nuclear DNA loci

BACKGROUND

Interspecific hybridization may lead to the introgression of genes and genomes across species barriers and contribute to a reticulate evolutionary pattern and thus taxonomic uncertainties. Since several previous studies have demonstrated that introgressive hybridization has occurred among some species within Lepus, therefore it is possible that introgressive hybridization events also occur among Chinese Lepus species and contribute to the current taxonomic confusion.

RESULTS

Data from four mtDNA genes, from 116 individuals, and one nuclear gene, from 119 individuals, provides the first evidence of frequent introgression events via historical and recent interspecific hybridizations among six Chinese Lepus species. Remarkably, the mtDNA of L. mandshuricus was completely replaced by mtDNA from L. timidus and L. sinensis. Analysis of the nuclear DNA sequence revealed a high proportion of heterozygous genotypes containing alleles from two divergent clades and that several haplotypes were shared among species, suggesting repeated and recent introgression. Furthermore, results from the present analyses suggest that Chinese hares belong to eight species.

CONCLUSION

This study provides a framework for understanding the patterns of speciation and the taxonomy of this clade. The existence of morphological intermediates and atypical mitochondrial gene genealogies resulting from frequent hybridization events likely contribute to the current taxonomic confusion of Chinese hares. The present study also demonstrated that nuclear gene sequence could offer a powerful complementary data set with mtDNA in tracing a complete evolutionary history of recently diverged species.

Adaptive introgression of anticoagulant rodent poison resistance by hybridization between old world mice

Polymorphisms in the vitamin K 2,3-epoxide reductase subcomponent 1 (vkorc1) of house mice (Mus musculus domesticus) can cause resistance to anticoagulant rodenticides such as warfarin [1-3]. Here we show that resistant house mice can also originate from selection on vkorc1 polymorphisms acquired from the Algerian mouse (M. spretus) through introgressive hybridization. We report on a polymorphic introgressed genomic region in European M. m. domesticus that stems from M. spretus, spans >10 Mb on chromosome 7, and includes the molecular target of anticoagulants vkorc1 [1-4]. We show that in the laboratory, the homozygous complete vkorc1 allele of M. spretus confers resistance when introgressed into M. m. domesticus. Consistent with selection on the introgressed allele after the introduction of rodenticides in the 1950s, we found signatures of selection in patterns of variation in M. m. domesticus. Furthermore, we detected adaptive protein evolution of vkorc1 in M. spretus (Ka/Ks = 1.54-1.93) resulting in radical amino acid substitutions that apparently cause anticoagulant tolerance in M. spretus as a pleiotropic effect. Thus, positive selection produced an adaptive, divergent, and pleiotropic vkorc1 allele in the donor species, M. spretus, which crossed a species barrier and produced an adaptive polymorphic trait in the recipient species, M. m. domesticus.

Genetic dissection of a key reproductive barrier between nascent species of house mice

Reproductive isolation between species is often caused by deleterious interactions among loci in hybrids. Finding the genes involved in these incompatibilities provides insight into the mechanisms of speciation. With recently diverged subspecies, house mice provide a powerful system for understanding the genetics of reproductive isolation early in the speciation process. Although previous studies have yielded important clues about the genetics of hybrid male sterility in house mice, they have been restricted to F1 sterility or incompatibilities involving the X chromosome. To provide a more complete characterization of this key reproductive barrier, we conducted an F2 intercross between wild-derived inbred strains from two subspecies of house mice, Mus musculus musculus and Mus musculus domesticus. We identified a suite of autosomal and X-linked QTL that underlie measures of hybrid male sterility, including testis weight, sperm density, and sperm morphology. In many cases, the autosomal loci were unique to a specific sterility trait and exhibited an effect only when homozygous, underscoring the importance of examining reproductive barriers beyond the F1 generation. We also found novel two-locus incompatibilities between the M. m. musculus X chromosome and M. m. domesticus autosomal alleles. Our results reveal a complex genetic architecture for hybrid male sterility and suggest a prominent role for reproductive barriers in advanced generations in maintaining subspecies integrity in house mice.

Applications of next generation sequencing in molecular ecology of non-model organisms

As most biologists are probably aware, technological advances in molecular biology during the last few years have opened up possibilities to rapidly generate large-scale sequencing data from non-model organisms at a reasonable cost. In an era when virtually any study organism can 'go genomic', it is worthwhile to review how this may impact molecular ecology. The first studies to put the next generation sequencing (NGS) to the test in ecologically well-characterized species without previous genome information were published in 2007 and the beginning of 2008. Since then several studies have followed in their footsteps, and a large number are undoubtedly under way. This review focuses on how NGS has been, and can be, applied to ecological, population genetic and conservation genetic studies of non-model species, in which there is no (or very limited) genomic resources. Our aim is to draw attention to the various possibilities that are opening up using the new technologies, but we also highlight some of the pitfalls and drawbacks with these methods. We will try to provide a snapshot of the current state of the art for this rapidly advancing and expanding field of research and give some likely directions for future developments.

Effects of polymorphism for locally adapted genes on rates of neutral introgression in structured populations

Adaptation to local conditions within demes balanced by migration can maintain polymorphisms for variants that reduce fitness in certain ecological contexts. Here, we address the effects of such polymorphisms on the rate of introgression of neutral marker genes, possibly genetically linked to targets of selection. Barriers to neutral gene flow are expected to increase with linkage to targets of local selection and with differences between demes in the frequencies of locally adapted alleles. This expectation is borne out under purifying and disruptive selection, regimes that promote monomorphism within demes. In contrast, overdominance within demes induces minimal barriers to neutral introgression even in the face of very large differences between demes in the frequencies of locally adapted alleles. Further, segregation distortion, a phenomenon observed in a number of interspecific hybrids, can in fact promote transmission by migrants to future generations at rates exceeding those of residents.

Measures of linkage disequilibrium among neighbouring SNPs indicate asymmetries across the house mouse hybrid zone

Theory predicts that naturally occurring hybrid zones between genetically distinct taxa can move over space and time as a result of selection and/or demographic processes, with certain types of hybrid zones being more or less likely to move. Determining whether a hybrid zone is stationary or moving has important implications for understanding evolutionary processes affecting interactions in hybrid populations. However, direct observations of hybrid zone movement are difficult to make unless the zone is moving rapidly. Here, evidence for movement in the house mouse Mus musculus domesticus × Mus musculus musculus hybrid zone is provided using measures of LD and haplotype structure among neighbouring SNP markers from across the genome. Local populations of mice across two transects in Germany and the Czech Republic were sampled, and a total of 1301 mice were genotyped at 1401 markers from the nuclear genome. Empirical measures of LD provide evidence for extinction and (re)colonization in single populations and, together with simulations, suggest hybrid zone movement because of either geography-dependent asymmetrical dispersal or selection favouring one subspecies over the other.

Subspecific origin and haplotype diversity in the laboratory mouse

Here we provide a genome-wide, high-resolution map of the phylogenetic origin of the genome of most extant laboratory mouse inbred strains. Our analysis is based on the genotypes of wild-caught mice from three subspecies of Mus musculus. We show that classical laboratory strains are derived from a few fancy mice with limited haplotype diversity. Their genomes are overwhelmingly Mus musculus domesticus in origin, and the remainder is mostly of Japanese origin. We generated genome-wide haplotype maps based on identity by descent from fancy mice and show that classical inbred strains have limited and non-randomly distributed genetic diversity. In contrast, wild-derived laboratory strains represent a broad sampling of diversity within M. musculus. Intersubspecific introgression is pervasive in these strains, and contamination by laboratory stocks has played a role in this process. The subspecific origin, haplotype diversity and identity by descent maps can be visualized using the Mouse Phylogeny Viewer (see URLs).

Reinforcement selection acting on the European house mouse hybrid zone

Behavioural isolation may lead to complete speciation when partial postzygotic isolation acts in the presence of divergent-specific mate-recognition systems. These conditions exist where Mus musculus musculus and M. m. domesticus come into contact and hybridize. We studied two mate-recognition signal systems, based on urinary and salivary proteins, across a Central European portion of the mouse hybrid zone. Introgression of the genomic regions responsible for these signals: the major urinary proteins (MUPs) and androgen binding proteins (ABPs), respectively, was compared to introgression at loci assumed to be nearly neutral and those under selection against hybridization. The preference of individuals taken from across the zone regarding these signals was measured in Y mazes, and we develop a model for the analysis of the transition of such traits under reinforcement selection. The strongest assortative preferences were found in males for urine and females for ABP. Clinal analyses confirm nearly neutral introgression of an Abp locus and two loci closely linked to the Abp gene cluster, whereas two markers flanking the Mup gene region reveal unexpected introgression. Geographic change in the preference traits matches our reinforcement selection model significantly better than standard cline models. Our study confirms that behavioural barriers are important components of reproductive isolation between the house mouse subspecies.

Speciation in the European rabbit (Oryctolagus cuniculus): islands of differentiation on the X chromosome and autosomes

Studies of gene flow between closely related taxa can provide insight into the genetic basis of speciation. To evaluate the importance of the X chromosome in reproductive isolation between subspecies of the European rabbit and to study the genomic scale over which islands of differentiation extend, we resequenced a total of 34 loci distributed along the X chromosome and chromosome 14. Previous studies based on few markers suggested that loci in centromeric regions were highly differentiated between rabbit subspecies, whereas loci in telomeric regions were less differentiated. Here, we confirmed this finding but also discovered remarkable variation in levels of differentiation among loci, with F(ST) values from nearly 0 to 1. Analyses using isolation-with-migration models suggest that this range appears to be largely explained by differential levels of gene flow among loci. The X chromosome was significantly more differentiated than the autosomes. On chromosome 14, differentiation decayed very rapidly at increasing distances from the centromere, but on the X chromosome distinct islands of differentiation encompassing several megabases were observed both at the centromeric region and along the chromosome arms. These findings support the idea that the X chromosome plays an important role in reproductive isolation between rabbit subspecies. These results also demonstrate the mosaic nature of the genome at species boundaries.

Inference of selection and stochastic effects in the house mouse hybrid zone

We explored the transition of 13 X-linked markers across two separate portions of the house mouse hybrid zone, asking whether such a comparison can distinguish the effects of selection from random factors. A heuristic search in the likelihood landscape revealed more complex likelihood profiles for data sampled in two-dimensional (2D) space relative to data sampled along a linear transect. Randomized resampling of localities analyzed for individual loci showed that deletion of sites away from the zone center can decrease cline width estimates whereas deletion of sites close to the center can significantly increase the width estimates. Deleting localities for all loci resulted in wider clines if the number of samples from the center was limited. The results suggest that, given the great variation in width estimates resulting from inclusion/exclusion of sampling sites, the geographic sampling design is important in hybrid zone studies and that our inferences should take into account measures of uncertainty such as support intervals. The comparison of the two transects indicates cline widths are narrower for loci in the central part of the X chromosome, suggesting selection is stronger in this region and genetic incompatibilities may have at least partly common architecture in the house mouse hybrid zone.