Contrasting Mode of Evolution at a Coat Color Locus in Wild and Domestic Pigs

Despite having only begun ~10,000 years ago, the process of domestication has resulted in a degree of phenotypic variation within individual species normally associated with much deeper evolutionary time scales. Though many variable traits found in domestic animals are the result of relatively recent human-mediated selection, uncertainty remains as to whether the modern ubiquity of long-standing variable traits such as coat color results from selection or drift, and whether the underlying alleles were present in the wild ancestor or appeared after domestication began. Here, through an investigation of sequence diversity at the porcine melanocortin receptor 1 (MC1R) locus, we provide evidence that wild and domestic pig (Sus scrofa) haplotypes from China and Europe are the result of strikingly different selection pressures, and that coat color variation is the result of intentional selection for alleles that appeared after the advent of domestication. Asian and European wild boar (evolutionarily distinct subspecies) differed only by synonymous substitutions, demonstrating that camouflage coat color is maintained by purifying selection. In domestic pigs, however, each of nine unique mutations altered the amino acid sequence thus generating coat color diversity. Most domestic MC1R alleles differed by more than one mutation from the wild-type, implying a long history of strong positive selection for coat color variants, during which time humans have cherry-picked rare mutations that would be quickly eliminated in wild contexts. This pattern demonstrates that coat color phenotypes result from direct human selection and not via a simple relaxation of natural selective pressures.     

Chinese white Rongchang pig does not have the dominant white allele of KIT but has the dominant black allele of MC1R

The mast/stem cell growth factor receptor (KIT) and melanocortin receptor 1 (MC1R) mutations are responsible for coat color phenotypes in domestic pigs. Rongchang is a Chinese indigenous pig breed with a white coat color phenotype. To investigate the genetic variability of the KIT and MC1R genes and their possible association with the coat color phenotype in this breed, a gene duplication and splice mutation of KIT were diagnosed in a sample of 93 unrelated Rongchang animals. The results show that Rongchang pigs have a single copy of KIT without the splice mutation at the first nucleotide of intron 17, indicating that the dominant white I allele of KIT is not responsible for their white phenotype. The KIT mRNA and MC1R coding sequences were also determined in this breed. Three putative amino acid substitutions were found in the KIT gene between Rongchang and Western white pigs, their association with the Rongchang white phenotype remains unknown. For the MC1R gene, Rongchang pigs were demonstrated to have the same dominant black allele (E(D1)) as other Chinese breeds, supporting the previous conclusion that Chinese and Western pigs have independent domestication origin. We also clarified that the Rongchang white phenotype was recessive to nonwhite color phenotypes. Our results provide a good starting point for the identification of the mutations underlying the white coat color in Rongchang pigs.     

Detecting population growth, selection and inherited fertility from haplotypic data in humans

The frequency of a rare mutant allele and the level of allelic association between this allele and one or several closely linked markers are frequently measured in genetic epidemiology. Both quantities are related to the time elapsed since the appearance of the mutation in the population and the intrinsic growth rate of the mutation (which may be different from the average population growth rate). Here, we develop a method that uses these two kinds of genetic data to perform a joint estimation of the age of the mutation and the minimum growth rate that is compatible with its present frequency. In absence of demographic data, it provides a useful estimate of population growth rate. When such data are available, contrasts among estimates from several loci allow demographic processes, affecting all loci similarly, to be distinguished from selection, affecting loci differently. Testing these estimates on populations for which data are available for several disorders shows good congruence with demographic data in some cases whereas in others higher growth rates are obtained, which may be the result of selection or hidden demographic processes.     

Estimation of 2Nes from temporal allele frequency data

We develop a new method for estimating effective population sizes, Ne, and selection coefficients, s, from time-series data of allele frequencies sampled from a single diallelic locus. The method is based on calculating transition probabilities, using a numerical solution of the diffusion process, and assuming independent binomial sampling from this diffusion process at each time point. We apply the method in two example applications. First, we estimate selection coefficients acting on the CCR5-delta 32 mutation on the basis of published samples of contemporary and ancient human DNA. We show that the data are compatible with the assumption of s = 0, although moderate amounts of selection acting on this mutation cannot be excluded. In our second example, we estimate the selection coefficient acting on a mutation segregating in an experimental phage population. We show that the selection coefficient acting on this mutation is approximately 0.43.     

Bait-ER: A Bayesian method to detect targets of selection in Evolve-and-Resequence experiments

For over a decade, experimental evolution has been combined with high-throughput sequencing techniques. In so-called Evolve-and-Resequence (E&R) experiments, populations are kept in the laboratory under controlled experimental conditions where their genomes are sampled and allele frequencies monitored. However, identifying signatures of adaptation in E&R datasets is far from trivial, and it is still necessary to develop more efficient and statistically sound methods for detecting selection in genome-wide data. Here, we present Bait-ER – a fully Bayesian approach based on the Moran model of allele evolution to estimate selection coefficients from E&R experiments. The model has overlapping generations, a feature that describes several experimental designs found in the literature. We tested our method under several different demographic and experimental conditions to assess its accuracy and precision, and it performs well in most scenarios. Nevertheless, some care must be taken when analysing trajectories where drift largely dominates and starting frequencies are low. We compare our method with other available software and report that ours has generally high accuracy even for trajectories whose complexity goes beyond a classical sweep model. Furthermore, our approach avoids the computational burden of simulating an empirical null distribution, outperforming available software in terms of computational time and facilitating its use on genome-wide data. We implemented and released our method in a new open-source software package that can be accessed at https://doi.org/10.5281/zenodo.7351736 .     

Simulating genealogies of selected alleles in a population of variable size

An importance-sampling method is presented that allows the simulation of the history of a selected allele in a population of variable size. A sample path describing the number of copies of an allele that arose as a single mutant is generated by simulating backwards from the current frequency until the allele is lost. The mathematical expectation of a quantity or statistic is then estimated by taking averages over replicate simulations, weighting each replicate by the ratio of its probabilities under the Markov chains for the forward and backwards processes. This method was used to find the average age of a selected allele in an exponentially growing population. In terms of the effect on average allele age, selection in favour of an allele is not equivalent to exponential growth. To generate gene genealogies of a sample of copies of a selected allele, the neutral coalescent model is simulated for the subpopulation containing only the selected allele. From the resulting intra-allelic genealogy, it is possible to calculate the likelihood of the selection intensity as a function of the observed level of variability at marker loci closely linked to the selected allele. This method was used to estimate the intensity of selection affecting the delta 32 allele at the CCR5 locus in Europeans and a mutant at the MLH1 locus associated with colorectal cancer in the Finnish population.     

Estimates of natural selection due to protein tertiary structure inform the ancestry of biallelic loci

We consider the inference of which of two alleles is ancestral when the alleles have a single nonsynonymous difference and when natural selection acts via protein tertiary structure. Whereas the probability that an allele is ancestral under neutrality is equal to its frequency, under selection this probability depends on allele frequency and on the magnitude and direction of selection pressure. Although allele frequencies can be well estimated from intraspecific data, small fitness differences have a large evolutionary impact but can be difficult to estimate with only intraspecific data. Methods for predicting aspects of phenotype from genotype can supplement intraspecific sequence data. Recently developed statistical techniques can assess effects of phenotypes, such as protein tertiary structure on molecular evolution. While these techniques were initially designed for comparing protein-coding genes from different species, the resulting interspecific inferences can be assigned population genetic interpretations to assess the effect of selection pressure, and we use them here along with intraspecific allele frequency data to estimate the probability that an allele is ancestral. We focus on 140 nonsynonymous single nucleotide polymorphisms of humans that are in proteins with known tertiary structures. We find that our technique for employing protein tertiary structure information yields some biologically plausible results but that it does not substantially improve the inference of ancestral human allele types.     

Estimating Selection Coefficients in Spatially Structured Populations from Time Series Data of Allele Frequencies

Inferring the nature and magnitude of selection is an important problem in many biological contexts. Typically when estimating a selection coefficient for an allele, it is assumed that samples are drawn from a panmictic population and that selection acts uniformly across the population. However, these assumptions are rarely satisfied. Natural populations are almost always structured, and selective pressures are likely to act differentially. Inference about selection ought therefore to take account of structure. We do this by considering evolution in a simple lattice model of spatial population structure. We develop a hidden Markov model based maximum-likelihood approach for estimating the selection coefficient in a single population from time series data of allele frequencies. We then develop an approximate extension of this to the structured case to provide a joint estimate of migration rate and spatially varying selection coefficients. We illustrate our method using classical data sets of moth pigmentation morph frequencies, but it has wide applications in settings ranging from ecology to human evolution.     

Exclusion of melanocortin-1 receptor (mc1r) and agouti as candidates for dominant black in dogs

The domestic dog exhibits a variety of coat colors that encompass a wide range of variation among different breeds. Very little is known about the molecular biology of dog pigmentation; current understanding is based mostly on traditional breeding experiments, which in some cases have suggested genetic interactions that are different from those reported in other mammals. We have examined the molecular genetics of dominant black, a uniform coat color characteristic of black Labrador retrievers or Newfoundlands that has been proposed to be caused by either variation in the melanocortin-1 receptor gene (Mc1r) or by variation in the Agouti gene (A). We identified several coding polymorphisms within Mc1r and several simple sequence repeat polymorphisms closely linked to A, and examined their inheritance in a Labrador retriever x greyhound cross that segregates dominant black. No single Mc1r allele was found consistently in animals carrying dominant black, and neither Mc1r nor A cosegregated with dominant black. These results refine our understanding of mammalian coat color inheritance and suggest that dominant black coat color in dogs is caused by a gene not previously implicated in pigment type switching.     

Platinum coat color locus in the deer mouse

Platinum coat color in the deer mouse, Peromyscus maniculatus, is an autosomal recessive trait marking a locus, pt, distinct from silver (si), albino (c), blonde (bl), brown (b), and agouti (a). Platinum deer mice are conspicuously pale, with light ears and tail stripe. The pewter trait is allelic with and phenotypically identical to platinum, and represents an independent recurrence of this mutant. The rate of recoveries of coat color mutations from wild deer mice is consistent with available data for recurring mutation rates balanced by strong selection against the recessive phenotype.     

Classic Selective Sweeps Revealed by Massive Sequencing in Cattle

Human driven selection during domestication and subsequent breed formation has likely left detectable signatures within the genome of modern cattle. The elucidation of these signatures of selection is of interest from the perspective of evolutionary biology, and for identifying domestication-related genes that ultimately may help to further genetically improve this economically important animal. To this end, we employed a panel of more than 15 million autosomal SNPs identified from re-sequencing of 43 Fleckvieh animals. We mainly applied two somewhat complementary statistics, the integrated Haplotype Homozygosity Score (iHS) reflecting primarily ongoing selection, and the Composite of Likelihood Ratio (CLR) having the most power to detect completed selection after fixation of the advantageous allele. We find 106 candidate selection regions, many of which are harboring genes related to phenotypes relevant in domestication, such as coat coloring pattern, neurobehavioral functioning and sensory perception including KIT, MITF, MC1R, NRG4, Erbb4, TMEM132D and TAS2R16, among others. To further investigate the relationship between genes with signatures of selection and genes identified in QTL mapping studies, we use a sample of 3062 animals to perform four genome-wide association analyses using appearance traits, body size and somatic cell count. We show that regions associated with coat coloring significantly (P<0.0001) overlap with the candidate selection regions, suggesting that the selection signals we identify are associated with traits known to be affected by selection during domestication. Results also provide further evidence regarding the complexity of the genetics underlying coat coloring in cattle. This study illustrates the potential of population genetic approaches for identifying genomic regions affecting domestication-related phenotypes and further helps to identify specific regions targeted by selection during speciation, domestication and breed formation of cattle. We also show that Linkage Disequilibrium (LD) decays in cattle at a much faster rate than previously thought.     

Estimating the time since the fixation of a beneficial allele

The fixation of a beneficial allele in a population leaves a well-characterized signature in patterns of nucleotide variation at linked sites. This signature can be used to estimate the time since fixation from patterns of polymorphism in extant individuals. I introduce a method to assess the support in polymorphism data for a recent episode of directional positive selection and to estimate the time since fixation. I summarize the polymorphism data by three statistics that carry information about levels of diversity, the allele frequency spectrum, and the extent of allelic associations. Simulations are then used to obtain a sample from the posterior distribution of the time since fixation, conditional on the observed summaries. I test the performance of the approach on simulated data and apply it to the gene tb1 in maize. The data support the recent fixation of a favored allele, consistent with what is known about the importance of tb1 in the domestication process of maize.     

褐色种皮大豆与其黄色种皮衍生亲本的表型及基因型比较

大豆种皮色在从野生大豆到栽培大豆的选择过程中逐渐由黑色变成黄色,是重要的形态标记,因此,大豆种皮色相关基因的研究无论是对进化理论研究还是育种实践都具有非常重要的意义。利用褐色种皮J1265-2大豆及其衍生亲本黄色种皮大豆J1265-1为材料,通过SSR引物扩增片段,检验遗传背景的异同,同时对控制种皮的候选基因GmF3’H进行扩增和测序分析。结果表明,褐色种皮和黄色种皮材料不仅用161对SSR分子标记检测没有发现差异,其褐色种皮候选基因GmF3’H的编码区及起始密码子上游1465 bp序列也是一致的。因此,证明褐色种皮J1265-2大豆与其衍生亲本黄色种皮大豆J1265-1为近等基因系,其控制褐色种皮的基因型与已报道的基因型不同。     

Sex ratio in silver foxes: effects of domestication and the star gene

The course of changes in secondary sex ratio (proportion of males at birth) in silver foxes bred at the fur farm of this Institute was analyzed. Data collected over several years of breeding of a domesticated (experimental) population selected for amenability to domestication and of a commercial (control) were compared. A tendency to increase in secondary sex ratio was demonstrated in both populations. However, the proportion of males at birth was higher in domestic foxes. This proportion, calculated from the combined data for 1978–1993, was 0.538±0.005 and 0.511±0.007 in the selected and commercial populations, respectively. The minimal departure of the observed sex ratio from 0.5 was demonstrated for litters with five pups, which is close to the average litter size in fox populations. The proportion of males increases with both increasing and decreasing litter size. An analysis of secondary sex ratio with respect to maternal age revealed a minimal departure of sex ratio from the expected in offspring from foxes of optimal reproductive age (2–4 years). An effect of the autosomal semidominant coat color mutation star on male excess at birth was also found: secondary sex ratio was higher (0.583±0.015) in offspring of mothers heterozygous for the star mutation than from standard types of the domesticated population. The increase in secondary sex ratio in the analyzed fox populations is viewed as a correlated response to selection for domestication. The hormonal mechanisms mediating the effects of both this selection and the star mutation on sex ratio at birth are discussed.     

Estimating haplotype frequencies by combining data from large DNA pools with database information

We assume that allele frequency data have been extracted from several large DNA pools, each containing genetic material of up to hundreds of sampled individuals. Our goal is to estimate the haplotype frequencies among the sampled individuals by combining the pooled allele frequency data with prior knowledge about the set of possible haplotypes. Such prior information can be obtained, for example, from a database such as HapMap. We present a Bayesian haplotyping method for pooled DNA based on a continuous approximation of the multinomial distribution. The proposed method is applicable when the sizes of the DNA pools and/or the number of considered loci exceed the limits of several earlier methods. In the example analyses, the proposed model clearly outperforms a deterministic greedy algorithm on real data from the HapMap database. With a small number of loci, the performance of the proposed method is similar to that of an EM-algorithm, which uses a multinormal approximation for the pooled allele frequencies, but which does not utilize prior information about the haplotypes. The method has been implemented using Matlab and the code is available upon request from the authors.     

Frequency Finder: a multi-source web application for collection of public allele frequencies of SNP markers

Publicly available single nucleotide polymorphism (SNP) allele frequencies are an important resource for the selection of genetic markers that may be most useful for gene mapping and association studies. Data mining these allele frequencies through disparate public databases and Websites is time consuming and can result in inconsistent findings. We have developed a web-based software tool, Frequency Finder, to acquire SNP allele frequencies from multiple public data sources and return a summarized result to the user. Our software optimizes and automates the search of candidate markers, decreasing the amount of time it would take to extract pertinent data manually. We have included several methods to output the data, including on-screen and as a compressed text file. We show that Frequency Finder accurately retrieves available frequency data from the available sources. Using this tool, we detect significant differences between Asian, African and Caucasian populations in the allele frequency spectra of 246 097 SNPs. While limited to public databases that provide web-based access to allele frequencies, Frequency Finder provides a single, user-friendly interface for retrieving allele frequencies for large batches of SNPs from multiple data sources.     

Assigning alleles to different loci in amplifications of duplicated loci

Duplicated loci, for example those associated with major histocompatibility complex (MHC) genes, often have similar DNA sequences that can be coamplified with a pair of primers. This results in genotyping difficulties and inaccurate analyses. Here, we present a method to assign alleles to different loci in amplifications of duplicated loci. This method simultaneously considers several factors that may each affect correct allele assignment. These are the sharing of identical alleles among loci, null alleles, copy number variation, negative amplification, heterozygote excess or heterozygote deficiency, and linkage disequilibrium. The possible multilocus genotypes are extracted from the alleles for each individual and weighted to estimate the allele frequencies. The likelihood of an allele configuration is calculated and is optimized with a heuristic algorithm. Monte‐Carlo simulations and three empirical MHC data sets are used as examples to evaluate the efficacy of our method under different conditions. Our new software, mhc‐typer V1.1, is freely available at https://github.com/huangkang1987/mhc-typer .     

A missense mutation in the gene for melanocyte-stimulating hormone receptor (MCIR) is associated with the chestnut coat color in horses

The melanocyte-stimulating hormone receptor gene (MC1R) is the major candidate gene for the chestnut coat color in horses since it is assumed to be controlled by an allele at the extension locus. MC1R sequences were PCR amplified from chestnut (e/e) and non-chestnut (E/−) horses. A single-strand conformation polymorphism was found that showed a complete association to the chestnut coat color among 144 horses representing 12 breeds. Sequence analysis revealed a single missense mutation (83Ser → Phe) in the MC1R allele associated with the chestnut color. The substitution occurs in the second transmembrane region, which apparently plays a key role in the molecule since substitutions associated with coat color variants in mice and cattle as well as red hair and fair skin in humans are found in this part of the molecule. We propose that the now reported mutation is likely to be the causative mutation for the chestnut coat color. The polymorphism can be detected with a simple PCR-RFLP test, since the mutation creates a TaqI restriction site in the chestnut allele. Received: 20 May 1996 / Accepted: 31 July 1996     

Generation and mapping of SCAR and CAPS markers linked to the seed coat color gene in Brassica napus using a genome-walking technique.

The yellow seed coat trait in No. 2127-17, a resynthesized purely yellow Brassica napus line, is controlled by a single partially dominant gene, Y. A double-haploid population derived from the F1 of No. 2127-17 x 'ZY821' was used to map the seed coat color phenotype. A combination of AFLP analysis and bulked segregant analysis identified 18 AFLP markers linked to the seed coat color trait. The 18 AFLP markers were mapped to a chromosomal region of 37.0 cM with an average of 2.0 cM between adjacent markers. Two markers, AFLP-K and AFLP-H, bracketed the Y locus in an interval of 1.0 cM, such that each was 0.5 cM away from the Y locus. Two other markers, AFLP-A and AFLP-B, co-segregated with the seed color gene. For ease of use in breeding programs, these 4 most tightly linked AFLP markers were converted into reliable PCR-based markers. SCAR-K, which was derived from AFLP-K, was assigned to linkage group 9 (N9) of a B. napus reference map consisting of 150 commonly used SSR (simple sequence repeat) markers. Furthermore, 2 SSR markers (Na14-E08 and Na10-B07) linked to SCAR-K on the reference map were reversely mapped to the linkage map constructed in this study, and also showed linkage to the Y locus. These linked markers would be useful for the transfer of the dominant allele Y from No. 2127-17 to elite cultivars using a marker-assisted selection strategy and would accelerate the cloning of the seed coat color gene.     

Stable two-allele polymorphisms maintained by fluctuating fitnesses and seed banks: protecting the blues in Linanthus parryae

Motivated by data demonstrating fluctuating relative and absolute fitnesses for white- versus blue-flowered morphs of the desert annual Linanthus parryae, we present conditions under which temporally fluctuating selection and fluctuating contributions to a persistent seed bank will maintain a stable single-locus polymorphism. In L. parryae, blue flower color is determined by a single dominant allele. To disentangle the underlying diversity-maintaining mechanism from the mathematical complications associated with departures from Hardy-Weinberg genotype frequencies and dominance, we successively analyze a haploid model, a diploid model with three distinguishable genotypes, and a diploid model with complete dominance. For each model, we present conditions for the maintenance of a stable polymorphism, then use a diffusion approximation to describe the long-term fluctuations associated with these polymorphisms. Our protected polymorphism analyses show that a genotype whose arithmetic and geometric mean relative fitnesses are both less than one can persist if its relative fitness exceeds one in years that produce the most offspring. This condition is met by data from a population of L. parryae whose white morph has higher fitness (seed set) only in years of relatively heavy rain fall. The data suggest that the observed polymorphism may be explained by fluctuating selection. However, the yearly variation in flower color frequencies cannot be fully explained by our simple models, which ignore age structure and possible selection in the seed bank. We address two additional questions--one mathematical, the other biological--concerning the applicability of diffusion approximations to intense selection and the applicability of long-term predictions to datasets spanning decades for populations with long-lived seed banks.