The impact of natural selection on an ABCC11 SNP determining earwax type

A nonsynonymous single nucleotide polymorphism (SNP), rs17822931-G/A (538G>A; Gly180Arg), in the ABCC11 gene determines human earwax type (i.e., wet or dry) and is one of most differentiated nonsynonymous SNPs between East Asian and African populations. A recent genome-wide scan for positive selection revealed that a genomic region spanning ABCC11, LONP2, and SIAH1 genes has been subjected to a selective sweep in East Asians. Considering the potential functional significance as well as the population differentiation of SNPs located in that region, rs17822931 is the most plausible candidate polymorphism to have undergone geographically restricted positive selection. In this study, we estimated the selection intensity or selection coefficient of rs17822931-A in East Asians by analyzing two microsatellite loci flanking rs17822931 in the African (HapMap-YRI) and East Asian (HapMap-JPT and HapMap-CHB) populations. Assuming a recessive selection model, a coalescent-based simulation approach suggested that the selection coefficient of rs17822931-A had been approximately 0.01 in the East Asian population, and a simulation experiment using a pseudo-sampling variable revealed that the mutation of rs17822931-A occurred 2006 generations (95% credible interval, 1,023-3,901 generations) ago. In addition, we show that absolute latitude is significantly associated with the allele frequency of rs17822931-A in Asian, Native American, and European populations, implying that the selective advantage of rs17822931-A is related to an adaptation to a cold climate. Our results provide a striking example of how local adaptation has played a significant role in the diversification of human traits.     

Positive selection in East Asians for an EDAR allele that enhances NF-kappaB activation

Genome-wide scans for positive selection in humans provide a promising approach to establish links between genetic variants and adaptive phenotypes. From this approach, lists of hundreds of candidate genomic regions for positive selection have been assembled. These candidate regions are expected to contain variants that contribute to adaptive phenotypes, but few of these regions have been associated with phenotypic effects. Here we present evidence that a derived nonsynonymous substitution (370A) in EDAR, a gene involved in ectodermal development, was driven to high frequency in East Asia by positive selection prior to 10,000 years ago. With an in vitro transfection assay, we demonstrate that 370A enhances NF-kappaB activity. Our results suggest that 370A is a positively selected functional genetic variant that underlies an adaptive human phenotype.     

蕨类植物rpoC1内含子缺失及其分子进化速率

rpoC1基因编码RNA聚合酶β°亚基蛋白, 在转录过程中与DNA模板结合, 与β亚基形成的β-β°亚基复合体构成RNA合成的催化中心。以rpoC1基因为研究对象, 在贝叶斯因子大于20的条件下, 用HyPhy软件位点模型检测到3个正选择位点和541个负选择位点; 用PAML软件位点模型检测到10个正选择位点, 其中3个位点的后验概率超过99%。此外, 基于最大似然法构建64种蕨类植物的系统发育树, 结合HyPhy软件分析rpoC1基因的转换率、颠换率、转换率/颠换率、同义替换率、非同义替换率以及同义替换率/非同义替换率, 探讨rpoC1基因内含子丢失与分子进化速率的关系。结果表明, rpoC1基因内含子缺失对转换率、颠换率以及非同义替换率有一定影响。     

Positive Selection in East Asians for an EDAR Allele that Enhances NF-κB Activation

Genome-wide scans for positive selection in humans provide a promising approach to establish links between genetic variants and adaptive phenotypes. From this approach, lists of hundreds of candidate genomic regions for positive selection have been assembled. These candidate regions are expected to contain variants that contribute to adaptive phenotypes, but few of these regions have been associated with phenotypic effects. Here we present evidence that a derived nonsynonymous substitution (370A) in EDAR, a gene involved in ectodermal development, was driven to high frequency in East Asia by positive selection prior to 10,000 years ago. With an in vitro transfection assay, we demonstrate that 370A enhances NF-κB activity. Our results suggest that 370A is a positively selected functional genetic variant that underlies an adaptive human phenotype.     

Synonymous and nonsynonymous polymorphisms versus divergences in bacterial genomes

Comparison of the ratio of nonsynonymous to synonymous polymorphisms within species with the ratio of nonsynonymous to synonymous substitutions between species has been widely used as a supposed indicator of positive Darwinian selection, with the ratio of these 2 ratios being designated as a neutrality index (NI). Comparison of genome-wide polymorphism within 12 species of bacteria with divergence from an outgroup species showed substantial differences in NI among taxa. A low level of nonsynonymous polymorphism at a locus was the best predictor of NI < 1, rather than a high level of nonsynonymous substitution between species. Moreover, genes with NI < 1 showed a strong tendency toward the occurrence of rare nonsynonymous polymorphisms, as expected under the action of ongoing purifying selection. Thus, our results are more consistent with the hypothesis that a high relative rate of between-species nonsynonymous substitution reflects mainly the action of purifying selection within species to eliminate slightly deleterious mutations rather than positive selection between species. This conclusion is consistent with previous results highlighting an important role of slightly deleterious variants in bacterial evolution and suggests caution in the use of the McDonald-Kreitman test and related statistics as tests of positive selection.     

High polymorphism at the human melanocortin 1 receptor locus

Variation in human skin/hair pigmentation is due to varied amounts of eumelanin (brown/black melanins) and phaeomelanin (red/yellow melanins) produced by the melanocytes. The melanocortin 1 receptor (MC1R) is a regulator of eu- and phaeomelanin production in the melanocytes, and MC1R mutations causing coat color changes are known in many mammals. We have sequenced the MC1R gene in 121 individuals sampled from world populations with an emphasis on Asian populations. We found variation at five nonsynonymous sites (resulting in the variants Arg67Gln, Asp84Glu, Val92Met, Arg151Cys, and Arg163Gln), but at only one synonymous site (A942G). Interestingly, the human consensus protein sequence is observed in all 25 African individuals studied, but at lower frequencies in the other populations examined, especially in East and Southeast Asians. The Arg163Gln variant is absent in the Africans studied, almost absent in Europeans, and at a low frequency (7%) in Indians, but is at an exceptionally high frequency (70%) in East and Southeast Asians. The MC1R gene in common and pygmy chimpanzees, gorilla, orangutan, and baboon was sequenced to study the evolution of MC1R. The ancestral human MC1R sequence is identical to the human consensus protein sequence, while MC1R varies considerably among higher primates. A comparison of the rates of substitution in genes in the melanocortin receptor family indicates that MC1R has evolved the fastest. In addition, the nucleotide diversity at the MC1R locus is shown to be several times higher than the average nucleotide diversity in human populations, possibly due to diversifying selection.     

A group of V3 sequences from human immunodeficiency virus type 1 subtype E non-syncytium-inducing, CCR5-using variants are resistant to positive selection pressure

In a human immunodeficiency virus type 1 (HIV-1)-infected individual, immune-pressure-mediated positive selection operates to maintain the antigenic polymorphism on the gp120 third variable (V3) loop. Recently, we suggested on the basis of sequencing C2/V3 segments from an HIV-1 subtype E-infected family that a V3 sequence lineage group of the non-syncytium-inducing (NSI) variants (group 1) was relatively resistant to positive selection pressure (35). To better understand the relationship between the intensity of positive selection pressure and cell tropism of the virus, we determined the linkage between each V3 genotype and its function of directing coreceptor preference and MT2 cell tropism. The biological characterization of a panel of V3 recombinant viruses showed that all of the group 1 V3 sequences could confer an NSI/CCR5-using (NSI/R5) phenotype on HIV-1(LAI), whereas the group 2 V3 sequence, which was more positively charged than the group 1 sequence, dictated mainly a syncytium-inducing, CXCR4-using (SI/X4) phenotype. Phylogenetic analysis of C2/V3 sequences encoding group 1 or 2 V3 suggested that the variants carrying group 1 V3 are the ancestors of the intrafamilial infection and persisted in the family, while the variants carrying group 2 V3 evolved convergently from the group 1 V3 variants during disease progression in the individuals. Finally, a statistical test showed that the V3 sequence that could dictate an NSI/R5 phenotype had a synonymous substitution rate significantly higher than the nonsynonymous substitution rate. These data suggest that V3 sequences of the subtype E NSI/R5 variants are more resistant to positive selection pressure than those of the SI/X4 variants.     

Founder virus population related to route of virus transmission: a determinant of intrahost human immunodeficiency virus type 1 evolution?

We and others have shown that in individual human immunodeficiency virus type 1 (HIV-1) infection, the adaptive evolution of HIV-1 is influenced by host immune competence. In this study, we tested the hypothesis that in addition to selective forces operating within the host, transmission bottlenecks have an impact on HIV-1 intrahost evolution. Therefore, we studied the intrahost evolution of the V3 region of the external glycoprotein gp120 of HIV-1 during the 3- and 5-year periods following seroconversion after parenteral versus sexual (male-to-male) transmission in 41 participants of the Amsterdam prospective cohorts of homosexual men (n = 31) and intravenous drug users (IVDUs; n = 10) who were AIDS free and had comparable numbers of CD4+ cells. We observed that HIV-1 strains in homosexual men accumulated over 5 years more nonsynonymous substitutions within the V3 loop than HIV-1 strains in IVDUs as a result of lower rates of nonsynonymous evolution in both the initial 3-year period from seroconversion and the following 2-year period as well as a larger proportion of nonsynonymous back substitutions in IVDUs. The mean numbers of synonymous substitutions did not differ between the two risk groups. Since HIV-1 strains in IVDUs could be distinguished from the viruses of homosexual men based on several nucleotide substitutions of which the most conserved is a synonymous substitution at the tip of the V3 loop (GGC pattern), we studied whether the founder virus population itself has an impact on the intrahost evolution of HIV-1. The mean number of nonsynonymous substitutions accumulated over 5 years within the V3 loop was lower in 10 IVDUs infected by the HIV-1 strains with the GGC signature than in 4 IVDUs infected by HIV-1 strains lacking this pattern, while the mean numbers of synonymous substitutions were similar in the two groups.     

Analysis of Ancestral and Functionally Relevant CD5 Variants in Systemic Lupus Erythematosus Patients

Objective

CD5 plays a crucial role in autoimmunity and is a well-established genetic risk factor of developing RA. Recently, evidence of positive selection has been provided for the CD5 Pro224-Val471 haplotype in East Asian populations. The aim of the present work was to further analyze the functional relevance of non-synonymous CD5 polymorphisms conforming the ancestral and the newly derived haplotypes (Pro224-Ala471 and Pro224-Val471, respectively) as well as to investigate the potential role of CD5 on the development of SLE and/or SLE nephritis.

Methods

The CD5 SNPs rs2241002 (C/T; Pro224Leu) and rs2229177 (C/T; Ala471Val) were genotyped using TaqMan allelic discrimination assays in a total of 1,324 controls and 681 SLE patients of Spanish origin. In vitro analysis of CD3-mediated T cell proliferative and cytokine response profiles of healthy volunteers homozygous for the above mentioned CD5 haplotypes were also analyzed.

Results

T-cell proliferation and cytokine release were significantly increased showing a bias towards to a Th2 profile after CD3 cross-linking of peripheral mononuclear cells from healthy individuals homozygous for the ancestral Pro224-Ala471 (CC) haplotype, compared to the more recently derived Pro224-Val471 (CT). The same allelic combination was statistically associated with Lupus nephritis.

Conclusion

The ancestral Ala471 CD5 allele confers lymphocyte hyper-responsiveness to TCR/CD3 cross-linking and is associated with nephritis in SLE patients.     

High Trans-ethnic Replicability of GWAS Results Implies Common Causal Variants

Genome-wide association studies (GWAS) have detected many disease associations. However, the reported variants tend to explain small fractions of risk, and there are doubts about issues such as the portability of findings over different ethnic groups or the relative roles of rare versus common variants in the genetic architecture of complex disease. Studying the degree of sharing of disease-associated variants across populations can help in solving these issues. We present a comprehensive survey of GWAS replicability across 28 diseases. Most loci and SNPs discovered in Europeans for these conditions have been extensively replicated using peoples of European and East Asian ancestry, while the replication with individuals of African ancestry is much less common. We found a strong and significant correlation of Odds Ratios across Europeans and East Asians, indicating that underlying causal variants are common and shared between the two ancestries. Moreover, SNPs that failed to replicate in East Asians map into genomic regions where Linkage Disequilibrium patterns differ significantly between populations. Finally, we observed that GWAS with larger sample sizes have detected variants with weaker effects rather than with lower frequencies. Our results indicate that most GWAS results are due to common variants. In addition, the sharing of disease alleles and the high correlation in their effect sizes suggest that most of the underlying causal variants are shared between Europeans and East Asians and that they tend to map close to the associated marker SNPs.     

The Light Skin Allele of SLC24A5 in South Asians and Europeans Shares Identity by Descent

Skin pigmentation is one of the most variable phenotypic traits in humans. A non-synonymous substitution (rs1426654) in the third exon of SLC24A5 accounts for lighter skin in Europeans but not in East Asians. A previous genome-wide association study carried out in a heterogeneous sample of UK immigrants of South Asian descent suggested that this gene also contributes significantly to skin pigmentation variation among South Asians. In the present study, we have quantitatively assessed skin pigmentation for a largely homogeneous cohort of 1228 individuals from the Southern region of the Indian subcontinent. Our data confirm significant association of rs1426654 SNP with skin pigmentation, explaining about 27% of total phenotypic variation in the cohort studied. Our extensive survey of the polymorphism in 1573 individuals from 54 ethnic populations across the Indian subcontinent reveals wide presence of the derived-A allele, although the frequencies vary substantially among populations. We also show that the geospatial pattern of this allele is complex, but most importantly, reflects strong influence of language, geography and demographic history of the populations. Sequencing 11.74 kb of SLC24A5 in 95 individuals worldwide reveals that the rs1426654-A alleles in South Asian and West Eurasian populations are monophyletic and occur on the background of a common haplotype that is characterized by low genetic diversity. We date the coalescence of the light skin associated allele at 22–28 KYA. Both our sequence and genome-wide genotype data confirm that this gene has been a target for positive selection among Europeans. However, the latter also shows additional evidence of selection in populations of the Middle East, Central Asia, Pakistan and North India but not in South India.     

Evolution of cd59 gene in mammals

The CD59-coding sequences were obtained from 5 mammals by PCR and BLAST, and combined with the available sequences in GenBank, the nucleotide substitution rates of mammalian cd59 were calcu- lated. Results of synonymous and nonsynonymous substitution rates revealed that cd59 experienced negative selection in mammals overall. Four sites experiencing positive selection were found by using "site-specific" model in PAML software. These sites were distributed on the molecular surface, of which 2 sites located in the key functional domain. Furthermore, "branch-site-specific" model detected 1 positive site in cd59a and cd59b lineages which underwent accelerated evolution caused by positive selection after gene duplication in mouse.     

Germline EPHB2 receptor variants in familial colorectal cancer

Familial clustering of colorectal cancer occurs in 15-20% of cases, however recognized cancer syndromes explain only a small fraction of this disease. Thus, the genetic basis for the majority of hereditary colorectal cancer remains unknown. EPHB2 has recently been implicated as a candidate tumor suppressor gene in colorectal cancer. The aim of this study was to evaluate the contribution of EPHB2 to hereditary colorectal cancer. We screened for germline EPHB2 sequence variants in 116 population-based familial colorectal cancer cases by DNA sequencing. We then estimated the population frequencies and characterized the biological activities of the EPHB2 variants identified. Three novel nonsynonymous missense alterations were detected. Two of these variants (A438T and G787R) result in significant residue changes, while the third leads to a conservative substitution in the carboxy-terminal SAM domain (V945I). The former two variants were found once in the 116 cases, while the V945I variant was present in 2 cases. Genotyping of additional patients with colorectal cancer and control subjects revealed that A438T and G787R represent rare EPHB2 alleles. In vitro functional studies show that the G787R substitution, located in the kinase domain, causes impaired receptor kinase activity and is therefore pathogenic, whereas the A438T variant retains its receptor function and likely represents a neutral polymorphism. Tumor tissue from the G787R variant case manifested loss of heterozygosity, with loss of the wild-type allele, supporting a tumor suppressor role for EPHB2 in rare colorectal cancer cases. Rare germline EPHB2 variants may contribute to a small fraction of hereditary colorectal cancer.     

Positive Darwinian selection on crustacean hyperglycemic hormone (CHH) of the green shore crab, Carcinus maenas

The tissue-specific expression and differential function of the crustacean hyperglycemic hormone (CHH) in Carcinus maenas indicate an interesting evolutionary history. Previous studies have shown that CHH from the sinus gland X-organ (XO-type) has hyperglycemic activity, whereas the CHH from the pericardial organ (PO-type) neither shows hyperglycemic activity nor it inhibits Y-organ ecdysteroid synthesis. Here we examined the types of selective pressures operating on the variants of CHH in Carcinus maenas. Maximum likelihood-based codon substitution analyses revealed that the variants of this neuropeptide in C. maenas have been subjected to positive Darwinian selection indicating adaptive evolution and functional divergence among the CHH variants leading to two unique groups (PO and XO-type). Although the average ratio of nonsynonymous to synonymous substitution (omega) for the entire coding region is 0.5096, few codon sites showed significantly higher omega (10.95). Comparison of models that incorporate positive selection (omega > 1) with models not incorporating positive selection (omega <1) at certain codon sites failed to reject (p=0) evidence of positive Darwinian selection.     

Evolution and biological characterization of human immunodeficiency virus type 1 subtype E gp120 V3 sequences following horizontal and vertical virus transmission in a single family

It has been suggested that immune-pressure-mediated positive selection operates to maintain the antigenic polymorphism on the third variable (V3) loop of the gp120 of human immunodeficiency virus type 1 (HIV-1). Here we present evidence, on the basis of sequencing 147 independently cloned env C2/V3 segments from a single family (father, mother, and their child), that the intensity of positive selection is related to the V3 lineage. Phylogenetic analysis and amino acid comparison of env C2/V3 and gag p17/24 regions indicated that a single HIV-1 subtype E source had infected the family. The analyses of unique env C2/V3 clones revealed that two V3 lineage groups had evolved in the parents. Group 1 was maintained with low variation in all three family members regardless of the clinical state or the length of infection, whereas group 2 was only present in symptomatic individuals and was more positively charged and diverse than group 1. Only virus isolates carrying the group 2 V3 sequences infected and induced syncytia in MT2 cells, a transformed CD4(+)-T-cell line. A statistically significant excess of nonsynonymous substitutions versus synonymous substitutions was demonstrated only for the group 2 V3 region. The data suggest that HIV-1 variants, possessing the more homogeneous group 1 V3 element and exhibiting the non-syncytium-inducing phenotype, persist in infected individuals independent of clinical status and appear to be more resistant to positive selection pressure.     

Characterisation of four major histocompatibility complex class II genes of the koala (Phascolarctos cinereus)

Major histocompatibility complex (MHC) class II molecules have an integral role in the adaptive immune response, as they bind and present antigenic peptides to T helper lymphocytes. In this study of koalas, species-specific primers were designed to amplify exon 2 of the MHC class II DA and DB genes, which contain much of the peptide-binding regions of the α and β chains. A total of two DA α1 domain variants and eight DA β1 (DAB), three DB α1 and five DB β1 variants were amplified from 20 koalas from two free-living populations from South East Queensland and the Port Macquarie region in northern New South Wales. We detected greater variation in the β1 than in the α1 domains as well as evidence of positive selection in DAB. The present study provides a springboard to future investigation of the role of MHC in disease susceptibility in koalas.     

Infection-associated FUT2 (Fucosyltransferase 2) genetic variation and impact on functionality assessed by in vivo studies

The secretor (Se)/nonsecretor (se) histo-blood group variation depends on the action of the FUT2 enzyme and has major implications for human susceptibility to infections. To characterize the functionality of FUT2 variants, we assessed the correlation between saliva phenotypes and sequence variation at the FUT2 gene in sixty seven individuals from northern Portugal. While most non-secretor haplotypes were found to carry the 428G > A nonsense mutation in association with a 739G > A missense substitution, we have also identified a recombinant haplotype carrying the 739*A allele together with the efficient 428*G variant in individuals with the Se phenotype. This finding suggested, in contrast to previous results, that the 739*A allele encodes an efficient Se allele. To test this hypothesis we evaluated the in vivo enzyme activity of full coding expression constructs in transient transfection of CHO-K1 cells using FACS (fluorescence-activated cell sorting) analysis and expression of type 2 and type 3 chain H structures as read out. We detected FUT2 activity for the 739*A expression construct, demonstrating that the 739G > A substitution is indeed not inactivating. In accordance with the hypothesis that FUT2 is under long standing balancing selection, we estimated that the time depth of FUT2 global genetic variation is as old as 3 million years. Age estimates of specific variants suggest that the 428G > A mutation occurred at least 1.87 million years ago while the 739G > A substitution is about 816,000 years old. The 385A > T missense mutation underlying the non-secretor phenotype in East Asians appears to be more recent and is likely to have occurred about 256,000 years ago.     

Not so different after all: a comparison of methods for detecting amino acid sites under selection

We consider three approaches for estimating the rates of nonsynonymous and synonymous changes at each site in a sequence alignment in order to identify sites under positive or negative selection: (1) a suite of fast likelihood-based "counting methods" that employ either a single most likely ancestral reconstruction, weighting across all possible ancestral reconstructions, or sampling from ancestral reconstructions; (2) a random effects likelihood (REL) approach, which models variation in nonsynonymous and synonymous rates across sites according to a predefined distribution, with the selection pressure at an individual site inferred using an empirical Bayes approach; and (3) a fixed effects likelihood (FEL) method that directly estimates nonsynonymous and synonymous substitution rates at each site. All three methods incorporate flexible models of nucleotide substitution bias and variation in both nonsynonymous and synonymous substitution rates across sites, facilitating the comparison between the methods. We demonstrate that the results obtained using these approaches show broad agreement in levels of Type I and Type II error and in estimates of substitution rates. Counting methods are well suited for large alignments, for which there is high power to detect positive and negative selection, but appear to underestimate the substitution rate. A REL approach, which is more computationally intensive than counting methods, has higher power than counting methods to detect selection in data sets of intermediate size but may suffer from higher rates of false positives for small data sets. A FEL approach appears to capture the pattern of rate variation better than counting methods or random effects models, does not suffer from as many false positives as random effects models for data sets comprising few sequences, and can be efficiently parallelized. Our results suggest that previously reported differences between results obtained by counting methods and random effects models arise due to a combination of the conservative nature of counting-based methods, the failure of current random effects models to allow for variation in synonymous substitution rates, and the naive application of random effects models to extremely sparse data sets. We demonstrate our methods on sequence data from the human immunodeficiency virus type 1 env and pol genes and simulated alignments.     

Selection and Reduced Population Size Cannot Explain Higher Amounts of Neandertal Ancestry in East Asian than in European Human Populations

It has been hypothesized that the greater proportion of Neandertal ancestry in East Asians than in Europeans is due to the fact that purifying selection is less effective at removing weakly deleterious Neandertal alleles from East Asian populations. Using simulations of a broad range of models of selection and demography, we have shown that this hypothesis cannot account for the higher proportion of Neandertal ancestry in East Asians than in Europeans. Instead, more complex demographic scenarios, most likely involving multiple pulses of Neandertal admixture, are required to explain the data.