Genomic Targets of Positive Selection in Giant Mice from Gough Island

A key challenge in understanding how organisms adapt to their environments is to identify the mutations and genes that make it possible. By comparing patterns of sequence variation to neutral predictions across genomes, the targets of positive selection can be located. We applied this logic to house mice that invaded Gough Island (GI), an unusual population that shows phenotypic and ecological hallmarks of selection. We used massively parallel short-read sequencing to survey the genomes of 14 GI mice. We computed a set of summary statistics to capture diverse aspects of variation across these genome sequences, used approximate Bayesian computation to reconstruct a null demographic model, and then applied machine learning to estimate the posterior probability of positive selection in each region of the genome. Using a conservative threshold, 1,463 5-kb windows show strong evidence for positive selection in GI mice but not in a mainland reference population of German mice. Disproportionate shares of these selection windows contain genes that harbor derived nonsynonymous mutations with large frequency differences. Over-represented gene ontologies in selection windows emphasize neurological themes. Inspection of genomic regions harboring many selection windows with high posterior probabilities pointed to genes with known effects on exploratory behavior and body size as potential targets. Some genes in these regions contain candidate adaptive variants, including missense mutations and/or putative regulatory mutations. Our results provide a genomic portrait of adaptation to island conditions and position GI mice as a powerful system for understanding the genetic component of natural selection.     

A Model-Based Approach for Identifying Signatures of Ancient Balancing Selection in Genetic Data

While much effort has focused on detecting positive and negative directional selection in the human genome, relatively little work has been devoted to balancing selection. This lack of attention is likely due to the paucity of sophisticated methods for identifying sites under balancing selection. Here we develop two composite likelihood ratio tests for detecting balancing selection. Using simulations, we show that these methods outperform competing methods under a variety of assumptions and demographic models. We apply the new methods to whole-genome human data, and find a number of previously-identified loci with strong evidence of balancing selection, including several HLA genes. Additionally, we find evidence for many novel candidates, the strongest of which is FANK1, an imprinted gene that suppresses apoptosis, is expressed during meiosis in males, and displays marginal signs of segregation distortion. We hypothesize that balancing selection acts on this locus to stabilize the segregation distortion and negative fitness effects of the distorter allele. Thus, our methods are able to reproduce many previously-hypothesized signals of balancing selection, as well as discover novel interesting candidates.     

Characterization of a divergent chromosome region in the willow warbler Phylloscopus trochilus using avian genomic resources

Genome scans have made it possible to find outlier markers thought to have been influenced by divergent selection in almost any wild population. However, the lack of genomic information in nonmodel species often makes it difficult to associate these markers with certain genes or chromosome regions. Furthermore, the extent of linkage disequilibrium (LD) in the genome will determine the density of markers required to identify the genes under selection. In this study, we investigated a chromosome region in the willow warbler Phylloscopus trochilus surrounding a single marker previously identified in a genome scan. We first located the marker in the assembled genome of another species, the zebra finch Taeniopygia guttata, and amplified surrounding sequences in Fennoscandian willow warblers. Within an investigated chromosome region of 7.3 Mb as mapped to the zebra finch genome, we observed elevated genetic differentiation between a southern and a northern population across a 2.5-Mb interval comprising numerous coding genes. Within the southern and northern populations, higher values of LD were mostly found between SNPs within the same locus, but extended across distantly situated loci when the analyses were restricted to sampling sites showing intermediate allele frequencies of southern and northern alleles. Our study shows that cross-species genome information is a useful resource to obtain candidate sequences adjacent to outlier markers in nonmodel species.     

Heterogeneous selection at specific loci in natural environments in Arabidopsis thaliana

Genetic variation for quantitative traits is often greater than that expected to be maintained by mutation in the face of purifying natural selection. One possible explanation for this observed variation is the action of heterogeneous natural selection in the wild. Here we report that selection on quantitative trait loci (QTL) for fitness traits in the model plant species Arabidopsis thaliana differs among natural ecological settings and genetic backgrounds. At one QTL, the allele that enhanced the viability of fall-germinating seedlings in North Carolina reduced the fecundity of spring-germinating seedlings in Rhode Island. Several other QTL experienced strong directional selection, but only in one site and seasonal cohort. Thus, different loci were exposed to selection in different natural environments. Selection on allelic variation also depended upon the genetic background. The allelic fitness effects of two QTL reversed direction depending on the genotype at the other locus. Moreover, alternative alleles at each of these loci caused reversals in the allelic fitness effects of a QTL closely linked to TFL1, a candidate developmental gene displaying nucleotide sequence polymorphism consistent with balancing selection. Thus, both environmental heterogeneity and epistatic selection may maintain genetic variation for fitness in wild plant species.     

Genome scans reveal candidate regions involved in the adaptation to host plant in the pea aphid complex

A major goal in evolutionary biology is to uncover the genetic basis of adaptation. Divergent selection exerted on ecological traits may result in adaptive population differentiation and reproductive isolation and affect differentially the level of genetic divergence along the genome. Genome‐wide scan of large sets of individuals from multiple populations is a powerful approach to identify loci or genomic regions under ecologically divergent selection. Here, we focused on the pea aphid, a species complex of divergent host races, to explore the organization of the genomic divergence associated with host plant adaptation and ecological speciation. We analysed 390 microsatellite markers located at variable distances from predicted genes in replicate samples of sympatric populations of the pea aphid collected on alfalfa, red clover and pea, which correspond to three common host‐adapted races reported in this species complex. Using a method that accounts for the hierarchical structure of our data set, we found a set of 11 outlier loci that show higher genetic differentiation between host races than expected under the null hypothesis of neutral evolution. Two of the outliers are close to olfactory receptor genes and three other nearby genes encoding salivary proteins. The remaining outliers are located in regions with genes of unknown functions, or which functions are unlikely to be involved in interactions with the host plant. This study reveals genetic signatures of divergent selection across the genome and provides an inventory of candidate genes responsible for plant specialization in the pea aphid, thereby setting the stage for future functional studies.     

Environmental heterogeneity and balancing selection in the acorn barnacle Semibalanus balanoides

The northern acorn barnacle Semibalans banlanoides occupies several intertidal microhabitats which vary greatly in their degree of physical stress. This environmental heterogeneity creates distinct selection regimes which can maintain genetic variation in natural populations. Despite considerable attention placed on the link between spatial variation in fitness and balancing selection at specific loci, experimental manipulations and fitness estimates for molecular polymorphisms have rarely been conducted in the wild. The aim of this transplant experiment was to manipulate the level of physical stress experienced by a cohort of barnacles in the field and then investigate the spatial variation in fitness for genotypes at three loci: two candidate allozymes and the mitochondrial DNA control region. The viability of mannose-6-phosphate isomerase (Mpi) genotypes was dependent on the level of physical stress experienced in the various treatments; alternative homozygotes were favoured in alternative high stress-low stress environments. In contrast, the fitness of genotypes at other loci was equivalent among treatments and unaffected by the manipulation. Evaluated in the light of balancing selection models, these data indicate that the presence of multiple environmental niches is sufficient to promote a stable Mpi polymorphism in barnacle populations and that allelic variation at this locus reflects the process of adaptation to the heterogeneous intertidal landscape.     

Plant genomes do a balancing act

Balancing selection is one mechanism that may explain why diversity is maintained in wild populations. However, relatively few examples of genes showing evidence of balancing selection have been identified, particularly in plants. In this issue, Reininga et al . (2009) present three Arabidopsis loci that show strong evidence of balancing selection. The loci, discovered using a genome-scanning approach, encode proteins with diverse predicted functions: starch synthesis and control of gene expression. These three genes were identified by scanning only a small fraction of the Arabidopsis genome, suggesting that balancing selection may be more prevalent than previously known.     

Identifying footprints of directional and balancing selection in marine and freshwater three-spined stickleback (Gasterosteus aculeatus) populations

Natural selection is expected to leave an imprint on the neutral polymorphisms at the adjacent genomic regions of a selected gene. While directional selection tends to reduce within-population genetic diversity and increase among-population differentiation, the reverse is expected under balancing selection. To identify targets of natural selection in the three-spined stickleback ( Gasterosteus aculeatus ) genome, 103 microsatellite and two indel markers including expressed sequence tags (EST) and quantitative trait loci (QTL)-associated loci, were genotyped in four freshwater and three marine populations. The results indicated that a high proportion of loci (14.7%) might be affected by balancing selection and a lower proportion (2.8%) by directional selection. The strongest signatures of directional selection were detected in a microsatellite locus and two indel markers located in the intronic regions of the Eda-gene coding for the number of lateral plates. Yet, other microsatellite loci previously found to be informative in QTL-mapping studies revealed no signatures of selection. Two novel microsatellite loci ( Stn12 and Stn90 ) located in chromosomes I and VIII, respectively, showed signals of directional selection and might be linked to genomic regions containing gene(s) important for adaptive divergence. Although the coverage of the total genomic content was relatively low, the predominance of balancing selection signals is in agreement with the contention that balancing, rather than directional selection is the predominant mode of selection in the wild.     

Functional and geographical differentiation of candidate balanced polymorphisms in Arabidopsis thaliana

Molecular population genetic analysis of three chromosomal regions in Arabidopsis thaliana suggested that balancing selection might operate to maintain variation at three novel candidate adaptive trait genes, including SOLUBLE STARCH SYNTHASE I (SSI) , PLASTID TRANSCRIPTIONALLY ACTIVE 7(PTAC7) , and BELL-LIKE HOMEODOMAIN 10 (BLH10). If balanced polymorphisms are indeed maintained at these loci, then we would expect to observe functional variation underlying the previously detected signatures of selection. We observe multiple replacement polymorphisms within and in the 32 amino acids just upstream of the protein–protein interacting BELL domain at the BLH10 locus. While no clear protein sequence differences are found between allele types in SSI and PTAC7, these two genes show evidence for allele-specific variation in expression levels. Geographical patterns of allelic differentiation seem consistent with population stratification in this species and a significant longitudinal cline was observed at all three candidate loci. These data support a hypothesis of balancing selection at all three candidate loci and provide a basis for more detailed functional work by identifying possible functional differences that might be selectively maintained.     

Neutral mutation as the source of genetic variation in life history traits

The mechanism underlying the maintenance of adaptive genetic variation is a long-standing question in evolutionary genetics. There are two concepts (mutation-selection balance and balancing selection) which are based on the phenotypic differences between alleles. Mutation - selection balance and balancing selection cannot properly explain the process of gene substitution, i.e. the molecular evolution of quantitative trait loci affecting fitness. I assume that such loci have non-essential functions (small effects on fitness), and that they have the potential to evolve into new functions and acquire new adaptations. Here I show that a high amount of neutral polymorphism at these loci can exist in real populations. Consistent with this, I propose a hypothesis for the maintenance of genetic variation in life history traits which can be efficient for the fixation of alleles with very small selective advantage. The hypothesis is based on neutral polymorphism at quantitative trait loci and both neutral and adaptive gene substitutions. The model of neutral - adaptive conversion (NAC) assumes that neutral alleles are not neutral indefinitely, and that in specific and very rare situations phenotypic (relative fitness) differences between them can appear. In this paper I focus on NAC due to phenotypic plasticity of neutral alleles. The important evolutionary consequence of NAC could be the increased adaptive potential of a population. Loci responsible for adaptation should be fast evolving genes with minimally discernible phenotypic effects, and the recent discovery of genes with such characteristics implicates them as suitable candidates for loci involved in adaptation.     

Genetic Variation of Microsatellite Loci in the Major Histocompatibility Complex (MHC) Region in the Southeast Asian House Mouse (<Emphasis Type="Italic">Mus musculus castaneus</Emphasis>)

Major histocompatibility complex (MHC) genes are the most polymorphic loci known for vertebrates. Here we employed five microsatellite loci closely linked to the MHC region in an attempt to study the amount of genetic variation in 19 populations of the southeast Asian house mouse (Mus musculus castaneus) in Taiwan. The overall polymorphism at the five loci was high (H_e = 0.713), and the level of polymorphism varied from locus to locus. Furthermore, in order to investigate if selection is operating on MHC genes in natural mouse populations, we compared the extent and pattern of genetic variation for the MHC-linked microsatellite loci (the MHC loci) with those for the microsatellite loci located outside the MHC region (the non-MHC loci). The number of alleles and the logarithm of variance in repeat number were significantly higher for the MHC loci than for the non-MHC loci, presumably reflecting linkage to a locus under balancing selection. Although three statistical tests used do not provide support for selection, their lack of support may be due to low statistical power of the tests, to weakness of selection, or to a profound effect of genetic drift reducing the signature of balancing selection. Our results also suggested that the populations in the central and the southwestern regions of Taiwan might be one part of a metapopulation structure.     

Polymorphism patterns in two tightly linked developmental genes,Idgf1 and Idgf3, of Drosophila melanogaster

A new developmental gene family, recently identified in D. melanogaster, has been called imaginal disc growth factors (IDGF) because the proteins promote growth of cell lineages derived from imaginal discs. These are the first genes reported that encode polypeptide factors with mitotic activity in invertebrates. Characteristics such as similar arrangement of introns and exons, small size, and different cytological localization make this family an excellent candidate for evolutionary studies. We focus on the loci Idgf1 and Idgf3, two genes that possess the most distinctive features. We examine the pattern of intra- and interspecific nucleotide variation in the sequences from 20 isogenic lines of D. melanogaster and sequences from D. simulans and D. yakuba. While MK, HKA, and Tajima's tests of neutrality fail to reject a neutral model of molecular evolution, Fu and Li's test with outgroup and McDonald's test suggest that balancing selection is modulating the evolution of the Idgf1 locus. The rate of recombination between the two loci is high enough to uncouple any linkage disequilibrium arising between Idgf1 and Idgf3, despite their close physical proximity.     

Sequence variability at three MHC loci of finless porpoises (<Emphasis Type="Italic">Neophocaena phocaenoides</Emphasis>)

Major histocompatibility complex (MHC) class II DQB and DRA genes and class I gene of finless porpoises (Neophocaena phocaenoides) were investigated by single-strand conformation polymorphism and sequence analysis. The DRA, DQB, and MHC-I loci each contained 5, 14, and 34 unique sequences, respectively, and considerable sequence variation was found at the MHC-I and DQB loci. Gene duplication was manifested as three to five distinct sequences at each of the DQB and MHC-I loci from some individuals, and these sequences at each of the two loci separately clustered into four groups (cluster A, B, C, and D) based on the phylogenetic trees. Phylogenetic reconstruction revealed a trans-species pattern of evolution. Relatively high rates of non-synonymous (dN) vs synonymous (dS) substitution in the peptide-binding region (PBR) suggested balancing selection for maintaining polymorphisms at the MHC-I and DQB loci. In contrast, one single locus with little sequence variation was detected in the DRA gene, and no non-synonymous substitutions in the PBR indicated no balancing selection on this gene.     

Genetic analysis of differentiation among breeding ponds reveals a candidate gene for local adaptation in Rana arvalis

One of the main questions in evolutionary and conservation biology is how geographical and environmental features of the landscape shape neutral and adaptive genetic variation in natural populations. The identification of genomic polymorphisms that account for adaptive variation can aid in finding candidate loci for local adaptation. Consequently, a comparison of spatial patterns in neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection at the landscape scale. Many amphibians breed in wetlands, which differ in environmental conditions and in the degree of isolation, enhancing the potential for local adaptation. We used microsatellite markers to measure genetic differentiation among 17 local populations of Rana arvalis breeding in a network of wetlands. We found that locus RC08604 deviated from neutral expectations, suggesting that it is a good candidate for directional selection. We used a genetic network analysis to show that the allele distribution in this locus is correlated with habitat characteristics, whereas this was not the case at neutral markers that displayed a different allele distribution and population network in the study area. The graph approach illustrated the genomic heterogeneity (neutral loci vs. the candidate locus for directional selection) of gene exchange and genetic divergence among populations under directional selection. Limited gene flow between wetlands was only observed at the candidate genomic region under directional selection. RC08604 is partially located inside an up‐regulated thyroid‐hormone receptor (TRβ) gene coordinating the expression of other genes during metamorphosis and appears to be linked with variation in larval life‐history traits found among R. arvalis populations. We suggest that directional selection on genes coding larval life‐history traits is strong enough to maintain the divergence in these genomic regions, reducing the effective recombination of locally adapted alleles but not in other regions of the genome. Integrating this knowledge into conservation plans at the landscape scale will improve the design of management strategies to preserve adaptive genetic diversity in wetland networks.     

A genome scan to detect candidate regions influenced by local natural selection in human populations

As human populations dispersed throughout the world, they were subjected to new selective forces, which must have led to local adaptation via natural selection and hence altered patterns of genetic variation. Yet, there are very few examples known in which such local selection has clearly influenced human genetic variation. A potential approach for detecting local selection is to screen random loci across the genome; those loci that exhibit unusually large genetic distances between human populations are then potential markers of genomic regions under local selection. We investigated this approach by genotyping 332 short tandem repeat (STR) loci in Africans and Europeans and calculating the genetic differentiation for each locus. Patterns of genetic diversity at these loci were consistent with greater variation in Africa and with local selection operating on populations as they moved out of Africa. For 11 loci exhibiting the largest genetic differences, we genotyped an additional STR locus located nearby; the genetic distances for these nearby loci were significantly larger than average. These genomic regions therefore reproducibly exhibit larger genetic distances between populations than the "average" genomic region, consistent with local selection. Our results demonstrate that genome scans are a promising means of identifying candidate regions that have been subjected to local selection.     

THE FITNESS CONSEQUENCES OF MULTIPLE-LOCUS HETEROZYGOSITY UNDER THE MULTIPLICATIVE OVERDOMINANCE AND INBREEDING DEPRESSION MODELS

There is a growing body of literature suggesting that the fitness of an individual increases with the observed number of heterozygous loci. Broad theoretical considerations indicate that under various sorts of balancing selection, this is what one should generally expect in a population of multiple-locus genotypes. To date, however, it has not been possible to distinguish between two potential explanations of the phenomenon. The first explanation is that the loci examined are themselves responsible for the fitness differences observed (or, equivalents, are very closely linked to those that do). The genetic variation in question is thought to be maintained in polymorphic equilibrium by some form of balancing selection. The second explanation assumes that the observed loci are themselves selectively irrelevant but that their heterozygosity reflects that of the total genome. Genomic heterozygosity is thought to be predictive of fitness, being an obverse measure of generalized inbreeding depression. We provide a formal derivation of an explicit relationship between fitness and multiple-locus genotype for a simple form of the first explanation, the multiplicative overdominance model. The inbreeding depression model is a degenerate special case of this more general formulation. A formal estimation and testing framework is constructed that should facilitate evaluation of the two models with empiric data on heterozygosity and fitness.     

THE OPPORTUNITY FOR BALANCING SELECTION IN EXPERIMENTAL POPULATIONS OF CAENORHABDITIS ELEGANS

The role of balancing selection in maintaining diversity during the evolution of sexual populations to novel environments is poorly understood. To address this issue, we studied the impact of two mating systems, androdioecy and dioecy, on genotype distributions during the experimental evolution of Caenorhabditis elegans. We analyzed the temporal trajectories of 334 single nucleotide polymorphisms, covering 1/3 of the genome, and found extensive allele frequency changes and little loss of heterozygosities after 100 generations. As modeled with numerical simulations, SNP differentiation was consistent with genetic drift and average fitness effects of 2%, assuming that selection acted independently at each locus. Remarkably, inbreeding by self‐fertilization was of little consequence to SNP differentiation. Modeling selection on deleterious recessive alleles suggests that the initial evolutionary dynamics can be explained by associative overdominance, but not the later stages because much lower heterozygosities would be maintained during experimental evolution. By contrast, models with selection on true overdominant loci can explain the heterozygote excess observed at all periods, particularly when negative epistasis or independent fitness effects were considered. Overall, these findings indicate that selection at single loci, including purging of recessive alleles, underlies most of the genetic differentiation accomplished during the experiment. Nonetheless, they also imply that maintenance of genetic diversity may in large part be due to balancing selection at multiple loci.     

Balancing selection in Pattern Recognition Receptor signalling pathways is associated with gene function and pleiotropy in a wild rodent

Pathogen‐mediated balancing selection is commonly considered to play an important role in the maintenance of genetic diversity, in particular in immune genes. However, the factors that may influence which immune genes are the targets of such selection are largely unknown. To address this, here we focus on Pattern Recognition Receptor (PRR) signalling pathways, which play a key role in innate immunity. We used whole‐genome resequencing data from a population of bank voles (Myodes glareolus) to test for associations between balancing selection, pleiotropy and gene function in a set of 123 PRR signalling pathway genes. To investigate the effect of gene function, we compared genes encoding (a) receptors for microbial ligands versus downstream signalling proteins, and (b) receptors recognizing components of microbial cell walls, flagella and capsids versus receptors recognizing features of microbial nucleic acids. Analyses based on the nucleotide diversity of full coding sequences showed that balancing selection primarily targeted receptor genes with a low degree of pleiotropy. Moreover, genes encoding receptors recognizing components of microbial cell walls etc. were more important targets of balancing selection than receptors recognizing nucleic acids. Tests for localized signatures of balancing selection in coding and noncoding sequences showed that such signatures were mostly located in introns, and more evenly distributed among different functional categories of PRR pathway genes. The finding that signatures of balancing selection in full coding sequences primarily occur in receptor genes, in particular those encoding receptors for components of microbial cell walls etc., is consistent with the idea that coevolution between hosts and pathogens is an important cause of balancing selection on immune genes.