Signatures of reproductive isolation in patterns of single nucleotide diversity across inbred strains of mice

Reproductive isolation is often caused by the disruption of genic interactions that evolve in geographically separate populations. Identifying the genomic regions and genes involved in these interactions, known as "Dobzhansky-Muller incompatibilities," can be challenging but is facilitated by the wealth of genetic markers now available in model systems. In recent years, the complete genome sequence and thousands of single nucleotide polymorphisms (SNPs) from laboratory mice, which are largely genetic hybrids between Mus musculus and M. domesticus, have become available. Here, we use these resources to locate genomic regions that may underlie reproductive isolation between these two species. Using genotypes from 332 SNPs that differ between wild-derived strains of M. musculus and M. domesticus, we identified several physically unlinked SNP pairs that show exceptional gametic disequilibrium across the lab strains. Conspecific alleles were associated in a disproportionate number of these cases, consistent with the action of natural selection against hybrid gene combinations. As predicted by the Dobzhansky-Muller model, this bias was differentially attributable to locus pairs for which one hybrid genotype was missing. We assembled a list of potential Dobzhansky-Muller incompatibilities from locus pairs that showed extreme associations (only three gametic types) among conspecific alleles. Two SNPs in this list map near known hybrid sterility loci on chromosome 17 and the X chromosome, allowing us to nominate partners for disrupted interactions involving these genomic regions for the first time. Together, these results indicate that patterns produced by speciation between M. musculus and M. domesticus are visible in the genomes of lab strains of mice, underscoring the potential of these genetic model organisms for addressing general questions in evolutionary biology.     

The evolutionary fate of heterogeneous gene duplications: A precarious overdominant equilibrium between environment,sublethality and complementation

Gene duplications occur at a high rate. Although most appear detrimental, some homogeneous duplications (identical gene copies) can be selected for beneficial increase in produced proteins. Heterogeneous duplications, which combine divergent alleles of a single locus, are seldom studied due to the paucity of empirical data. We investigated their role in an ongoing adaptive process at the ace‐1 locus in Culex pipiens mosquitoes. We assessed the worldwide diversity of the ace‐1 alleles (single‐copy, susceptible S and insecticide‐resistant R, and duplicated D that pair one S and one R copy), analysed their phylogeography and measured their fitness to understand their early dynamics using population genetics models. It provides a coherent and comprehensive evolutionary scenario. We show that D alleles are present in most resistant populations and display a higher diversity than R alleles (27 vs. 4). Most appear to result from independent unequal crossing‐overs between local single‐copy alleles, suggesting a recurrent process. Most duplicated alleles have a limited geographic distribution, probably resulting from their homozygous sublethality (HS phenotype). In addition, heterozygotes carrying different HS D alleles showed complementation, indicating different recessive lethal mutations. Due to mosaic insecticide control practices, balancing selection (overdominance) plays a key role in the early dynamics of heterogeneous duplicated alleles; it also favours a high local polymorphism of HS D alleles in natural populations (overdominance reinforced by complementation). Overall, our study shows that the evolutionary fate of heterogeneous duplications (and their long‐term role) depends on finely balanced selective pressures due to the environment and to their genomic structure.     

The population genetics of X-autosome synthetic lethals and steriles

Epistatic interactions are widespread, and many of these interactions involve combinations of alleles at different loci that are deleterious when present in the same individual. The average genetic environment of sex-linked genes differs from that of autosomal genes, suggesting that the population genetics of interacting X-linked and autosomal alleles may be complex. Using both analytical theory and computer simulations, we analyzed the evolutionary trajectories and mutation-selection balance conditions for X-autosome synthetic lethals and steriles. Allele frequencies follow a set of fundamental trajectories, and incompatible alleles are able to segregate at much higher frequencies than single-locus expectations. Equilibria exist, and they can involve fixation of either autosomal or X-linked alleles. The exact equilibrium depends on whether synthetic alleles are dominant or recessive and whether fitness effects are seen in males, females, or both sexes. When single-locus fitness effects and synthetic incompatibilities are both present, population dynamics depend on the dominance of alleles and historical contingency (i.e., whether X-linked or autosomal mutations occur first). Recessive synthetic lethality can result in high-frequency X-linked alleles, and dominant synthetic lethality can result in high-frequency autosomal alleles. Many X-autosome incompatibilities in natural populations may be cryptic, appearing to be single-locus effects because one locus is fixed. We also discuss the implications of these findings with respect to standing genetic variation and the origins of Haldane's rule.     

Hybrid fitness in a locally adapted parasite

The parasite (Red Queen) hypothesis for the maintenance of sexual reproduction and genetic diversity assumes that host-parasite interactions result from tight genetic specificity. Hence, hybridization between divergent parasite populations would be expected to disrupt adaptive gene combinations, leading to reduced infectivity on exposure to parental sympatric hosts, as long as gene effects are nonadditive. In contrast, hybridization would not cause reduced infectivity on allopatric hosts unless the divergent parasite populations possess alleles that are intrinsically incompatible when they are combined. In three different experiments, we compared the infectivity of locally adapted parasite (trematode) populations with that of F(1) hybrid parasites when exposed to host (snail) populations that were sympatric to one of the two parasite populations. We tested for intrinsic genetic incompatibilities in two experiments by including one host population that was allopatric to both parasite populations. As predicted, when the target host populations were sympatric to the parasite populations, the hybrids were significantly less infective than the parental average, while hybrid parasites on allopatric hosts were not, thereby ruling out intrinsic genetic incompatibilities. The results are consistent with nonadditive gene effects and tightly specific host-driven selection underlying parasite divergence, as envisioned by coevolutionary theory and the Red Queen hypothesis.     

Mapping quantitative trait loci in multiple populations of Arabidopsis thaliana identifies natural allelic variation for trichome density

The majority of biological traits are genetically complex. Mapping the quantitative trait loci (QTL) that determine these phenotypes is a powerful means for estimating many parameters of the genetic architecture for a trait and potentially identifying the genes responsible for natural variation. Typically, such experiments are conducted in a single mapping population and, therefore, have only the potential to reveal genomic regions that are polymorphic between the progenitors of the population. What remains unclear is how well the QTL identified in any one mapping experiment characterize the genetics that underlie natural variation in traits. Here we provide QTL mapping data for trichome density from four recombinant inbred mapping populations of Arabidopsis thaliana. By aligning the linkage maps for these four populations onto a common physical map, the results from each experiment were directly compared. Seven of the nine QTL identified are population specific while two were mapped in all four populations. Our results show that many lineage-specific alleles that either increase or decrease trichome density persist in natural populations and that most of this genetic variation is additive. More generally, these findings suggest that the use of multiple populations holds great promise for better understanding the genetic architecture of natural variation.     

The genomics of incompatibility factors and sex determination in hybridizing species of Cottus (Pisces)

Cottus rhenanus and Cottus perifretum have formed hybrid lineages and narrow hybrid zones that can be best explained through the action of natural selection. However, the underlying selective forces as well as their genomic targets are not well understood. This study identifies genomic regions in the parental species that cause hybrid incompatibilities and tests whether these manifest in a sex-specific manner to learn about processes that affect natural hybridization in Cottus. Interspecific F₂ crosses were analyzed for 255 markers for genetic mapping and to detect transmission distortion as a sign for genetic incompatibilities. The Cottus map consists of 24 linkage groups with a total length of 1575.4 cM. A male heterogametic (XY) sex determination region was found on different linkage groups in the two parental species. Genetic incompatibilities were incomplete, varied among individuals and populations and were not associated with the heterogametic sex. The variance between populations and individuals makes it unlikely that there are species-specific incompatibility loci that could affect the gene pool of natural hybrids in a simple and predictable way. Conserved synteny with sequenced fish genomes permits to genetically study the Cottus genome through the transfer of genomic information from the model fish species. Homology relationships of candidate genomic regions in Cottus indicate that sex determination is not based on the same genomic regions found in other fish species. This suggests a fast evolutionary turnover of the genetic basis of sex determination that, together with the small size of the heterogametic regions, may contribute to the absence of fitness effects related to the Haldane''s rule.     

Natural genetic variation in Arabidopsis: tools, traits and prospects for evolutionary ecology

BACKGROUND: The model plant Arabidopsis thaliana (Arabidopsis) shows a wide range of genetic and trait variation among wild accessions. Because of its unparalleled biological and genomic resources, the potential of Arabidopsis for molecular genetic analysis of this natural variation has increased dramatically in recent years. SCOPE: Advanced genomics has accelerated molecular phylogenetic analysis and gene identification by quantitative trait loci (QTL) mapping and/or association mapping in Arabidopsis. In particular, QTL mapping utilizing natural accessions is now becoming a major strategy of gene isolation, offering an alternative to artificial mutant lines. Furthermore, the genomic information is used by researchers to uncover the signature of natural selection acting on the genes that contribute to phenotypic variation. The evolutionary significance of such genes has been evaluated in traits such as disease resistance and flowering time. However, although molecular hallmarks of selection have been found for the genes in question, a corresponding ecological scenario of adaptive evolution has been difficult to prove. Ecological strategies, including reciprocal transplant experiments and competition experiments, and utilizing near-isogenic lines of alleles of interest will be a powerful tool to measure the relative fitness of phenotypic and/or allelic variants. CONCLUSIONS: As the plant model organism, Arabidopsis provides a wealth of molecular background information for evolutionary genetics. Because genetic diversity between and within Arabidopsis populations is much higher than anticipated, combining this background information with ecological approaches might well establish Arabidopsis as a model organism for plant evolutionary ecology.     

Temperature stress increases hybrid incompatibilities in the parasitic wasp genus Nasonia

Hybrid incompatibilities, measured as mortality and sterility, are caused by the disruption of gene interactions. They are important post-zygotic isolation barriers to species hybridization, and much effort is put into the discovery of the genes underlying these incompatibilities. In hybridization studies of the haplodiploid parasitic wasp genus Nasonia, genic incompatibilities have been shown to affect mortality and sterility. The genomic regions associated with mortality have been found to depend on the cytotype of the hybrids and thus suggest cytonuclear incompatibilities. As environmental conditions can affect gene expression and gene interaction, we here investigate the effect of developmental temperature on sterility and mortality in Nasonia hybrids. Results show that extreme temperatures strongly affect both hybrid sterility (mainly spermatogenic failure) and mortality. Molecular mapping revealed that extreme temperatures increase transmission ratio distortion of parental alleles at incompatible loci, and thus, cryptic incompatible loci surface under temperature stress that remain undiscovered under standard temperatures. Our results underline the sensitivity of hybrid incompatibilities to environmental factors and the effects of unstable epistasis.     

Invasion of gene duplication through masking for maladaptive gene flow

Gene duplication can increase an organism's ability to mask the effect of deleterious alleles present in the population, but this is typically a small effect when the source of the genetic variation is mutation. Migration can introduce orders of magnitude more deleterious alleles per generation and may therefore be an important force acting on the structure of genomes. Using formal analytical methods, we study the invasion of haplotypes containing two copies of the resident allele, assuming that a single-locus equilibrium is already established in a continent-island model of migration. Provided that the immigrant allele can be completely masked by multiple functional gene copies, a new duplication will deterministically spread so long as duplicate haplotypes are, on average, fitter than single-copy haplotypes. When fitness depends on the number of immigrant allele copies and their masking ability then the threshold for invasion depends on the rate of immigration and the rate of recombination between the gene copies. Results from several special cases, including formation of protein dimers and Dobzhansky-Muller incompatibilities, suggest that duplications can invade in a wide range of selection regimes. We hypothesize that duplication in response to gene flow may provide an explanation for the high levels of polymorphism in gene copy number observed in natural populations.     

The Pattern of Polymorphism in Arabidopsis thaliana

We resequenced 876 short fragments in a sample of 96 individuals of Arabidopsis thaliana that included stock center accessions as well as a hierarchical sample from natural populations. Although A. thaliana is a selfing weed, the pattern of polymorphism in general agrees with what is expected for a widely distributed, sexually reproducing species. Linkage disequilibrium decays rapidly, within 50 kb. Variation is shared worldwide, although population structure and isolation by distance are evident. The data fail to fit standard neutral models in several ways. There is a genome-wide excess of rare alleles, at least partially due to selection. There is too much variation between genomic regions in the level of polymorphism. The local level of polymorphism is negatively correlated with gene density and positively correlated with segmental duplications. Because the data do not fit theoretical null distributions, attempts to infer natural selection from polymorphism data will require genome-wide surveys of polymorphism in order to identify anomalous regions. Despite this, our data support the utility of A. thaliana as a model for evolutionary functional genomics.     

The pattern of polymorphism in Arabidopsis thaliana

We resequenced 876 short fragments in a sample of 96 individuals of Arabidopsis thaliana that included stock center accessions as well as a hierarchical sample from natural populations. Although A. thaliana is a selfing weed, the pattern of polymorphism in general agrees with what is expected for a widely distributed, sexually reproducing species. Linkage disequilibrium decays rapidly, within 50 kb. Variation is shared worldwide, although population structure and isolation by distance are evident. The data fail to fit standard neutral models in several ways. There is a genome-wide excess of rare alleles, at least partially due to selection. There is too much variation between genomic regions in the level of polymorphism. The local level of polymorphism is negatively correlated with gene density and positively correlated with segmental duplications. Because the data do not fit theoretical null distributions, attempts to infer natural selection from polymorphism data will require genome-wide surveys of polymorphism in order to identify anomalous regions. Despite this, our data support the utility of A. thaliana as a model for evolutionary functional genomics.     

GENE‐DOSAGE EFFECTS ON FITNESS IN RECENT ADAPTIVE DUPLICATIONS: ace‐1 IN THE MOSQUITO CULEX PIPIENS

Gene duplications have long been advocated to contribute to the evolution of new functions. The role of selection in their early spread is more controversial. Unless duplications are favored for a direct benefit of increased expression, they are likely detrimental. In this article, we investigated the case of duplications favored because they combine already functionally divergent alleles. Their gene‐dosage/fitness relations are poorly known because selection may operate on both overall expression and duplicates relative dosage. Using the well‐documented case of Culex pipiens resistance to insecticides, we compared strains with various ace‐1 allele combinations, including two duplicated alleles carrying both susceptible and resistant copies. The overall protein activity was nearly additive, but, surprisingly, fitness correlated better with the relative proportion of susceptible and resistant copies rather than any absolute measure of activity. Gene dosage is thus crucial, duplications stabilizing a “heterozygote” phenotype. It corroborates the view that these were favored because they fix a permanent heterosis, thereby solving the irreducible trade‐off between resistance and synaptic transmission. Moreover, we showed that the contrasted successes of the two duplicated alleles in natural populations depend on genetic changes unrelated to ace‐1, confirming the probable implication of recessive sublethal mutations linked to structural rearrangements in some duplications.     

Sexual conflict and the maintenance of genetic variation in natural populations

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

Rapid Parallel Adaptation to Anthropogenic Heavy Metal Pollution

The impact of human-mediated environmental change on the evolutionary trajectories of wild organisms is poorly understood. In particular, capacity of species to adapt rapidly (in hundreds of generations or less), reproducibly and predictably to extreme environmental change is unclear. Silene uniflora is predominantly a coastal species, but it has also colonized isolated, disused mines with phytotoxic, zinc-contaminated soils. To test whether rapid, parallel adaptation to anthropogenic pollution has taken place, we used reduced representation sequencing (ddRAD) to reconstruct the evolutionary history of geographically proximate mine and coastal population pairs and found largely independent colonization of mines from different coastal sites. Furthermore, our results show that parallel evolution of zinc tolerance has occurred without gene flow spreading adaptive alleles between mine populations. In genomic regions where signatures of selection were detected across multiple mine-coast pairs, we identified genes with functions linked to physiological differences between the putative ecotypes, although genetic differentiation at specific loci is only partially shared between mine populations. Our results are consistent with a complex, polygenic genetic architecture underpinning rapid adaptation. This shows that even under a scenario of strong selection and rapid adaptation, evolutionary responses to human activities (and other environmental challenges) may be idiosyncratic at the genetic level and, therefore, difficult to predict from genomic data.     

X‐AUTOSOME INCOMPATIBILITIES IN DROSOPHILA MELANOGASTER: TESTS OF HALDANE'S RULE AND GEOGRAPHIC PATTERNS WITHIN SPECIES

Substantial genetic variation exists in natural populations of Drosophila melanogaster. This segregating variation includes alleles at different loci that interact to cause lethality or sterility (synthetic incompatibilities). Fitness epistasis in natural populations has important implications for speciation and the rate of adaptive evolution. To assess the prevalence of epistatic fitness interactions, we placed naturally occurring X chromosomes into genetic backgrounds derived from different geographic locations. Considerable amounts of synthetic incompatibilities were observed between X chromosomes and autosomes: greater than 44% of all combinations were either lethal or sterile. Sex‐specific lethality and sterility were also tested to determine whether Haldane's rule holds for within‐species variation. Surprisingly, we observed an excess of female sterility in genotypes that were homozygous, but not heterozygous, for the X chromosome. The recessive nature of these incompatibilities is similar to that predicted for incompatibilities underlying Haldane's rule. Our study also found higher levels of sterility and lethality for genomes that contain chromosomes from different geographical regions. These findings are consistent with the view that genomes are coadapted gene complexes and that geography affects the likelihood of epistatic fitness interactions.     

cis-Regulatory changes in Kit ligand expression and parallel evolution of pigmentation in sticklebacks and humans

Dramatic pigmentation changes have evolved within most vertebrate groups, including fish and humans. Here we use genetic crosses in sticklebacks to investigate the parallel origin of pigmentation changes in natural populations. High-resolution mapping and expression experiments show that light gills and light ventrums map to a divergent regulatory allele of the Kit ligand (Kitlg) gene. The divergent allele reduces expression in gill and skin tissue and is shared by multiple derived freshwater populations with reduced pigmentation. In humans, Europeans and East Asians also share derived alleles at the KITLG locus. Strong signatures of selection map to regulatory regions surrounding the gene, and admixture mapping shows that the KITLG genomic region has a significant effect on human skin color. These experiments suggest that regulatory changes in Kitlg contribute to natural variation in vertebrate pigmentation, and that similar genetic mechanisms may underlie rapid evolutionary change in fish and humans.     

Duplication and Gene Conversion in the Drosophila melanogaster Genome

Using the genomic sequences of Drosophila melanogaster subgroup, the pattern of gene duplications was investigated with special attention to interlocus gene conversion. Our fine-scale analysis with careful visual inspections enabled accurate identification of a number of duplicated blocks (genomic regions). The orthologous parts of those duplicated blocks were also identified in the D. simulans and D. sechellia genomes, by which we were able to clearly classify the duplicated blocks into post- and pre-speciation blocks. We found 31 post-speciation duplicated genes, from which the rate of gene duplication (from one copy to two copies) is estimated to be 1.0×10⁻⁹ per single-copy gene per year. The role of interlocus gene conversion was observed in several respects in our data: (1) synonymous divergence between a duplicated pair is overall very low. Consequently, the gene duplication rate would be seriously overestimated by counting duplicated genes with low divergence; (2) the sizes of young duplicated blocks are generally large. We postulate that the degeneration of gene conversion around the edges could explain the shrinkage of “identifiable” duplicated regions; and (3) elevated paralogous divergence is observed around the edges in many duplicated blocks, supporting our gene conversion–degeneration model. Our analysis demonstrated that gene conversion between duplicated regions is a common and genome-wide phenomenon in the Drosophila genomes, and that its role should be especially significant in the early stages of duplicated genes. Based on a population genetic prediction, we applied a new genome-scan method to test for signatures of selection for neofunctionalization and found a strong signature in a pair of transporter genes.     

Single nucleotide polymorphism and haplotype effects associated with somatic cell score in German Holstein cattle

Background

To better understand the genetic determination of udder health, we performed a genome-wide association study (GWAS) on a population of 2354 German Holstein bulls for which daughter yield deviations (DYD) for somatic cell score (SCS) were available. For this study, we used genetic information of 44 576 informative single nucleotide polymorphisms (SNPs) and 11 725 inferred haplotype blocks.

Results

When accounting for the sub-structure of the analyzed population, 16 SNPs and 10 haplotypes in six genomic regions were significant at the Bonferroni threshold of P ≤ 1.14 × 10^-6. The size of the identified regions ranged from 0.05 to 5.62 Mb. Genomic regions on chromosomes 5, 6, 18 and 19 coincided with known QTL affecting SCS, while additional genomic regions were found on chromosomes 13 and X. Of particular interest is the region on chromosome 6 between 85 and 88 Mb, where QTL for mastitis traits and significant SNPs for SCS in different Holstein populations coincide with our results. In all identified regions, except for the region on chromosome X, significant SNPs were present in significant haplotypes. The minor alleles of identified SNPs on chromosomes 18 and 19, and the major alleles of SNPs on chromosomes 6 and X were favorable for a lower SCS. Differences in somatic cell count (SCC) between alternative SNP alleles reached 14 000 cells/mL.

Conclusions

The results support the polygenic nature of the genetic determination of SCS, confirm the importance of previously reported QTL, and provide evidence for the segregation of additional QTL for SCS in Holstein cattle. The small size of the regions identified here will facilitate the search for causal genetic variations that affect gene functions.