Complex signatures of natural selection at the Duffy blood group locus

Cis-regulatory elements and human evolution

Modification of gene regulation has long been considered an important force in human evolution, particularly through changes to cis-regulatory elements (CREs) that function in transcriptional regulation. For decades, however, the study of cis-regulatory evolution was severely limited by the available data. New data sets describing the locations of CREs and genetic variation within and between species have now made it possible to study CRE evolution much more directly on a genome-wide scale. Here, we review recent research on the evolution of CREs in humans based on large-scale genomic data sets. We consider inferences based on primate divergence, human polymorphism, and combinations of divergence and polymorphism. We then consider 'new frontiers' in this field stemming from recent research on transcriptional regulation.

Adaptive divergence with gene flow in incipient speciation of Miscanthus floridulus/sinensis complex (Poaceae)

Young incipient species provide ideal materials for untangling the process of ecological speciation in the presence of gene flow. The Miscanthus floridulus/sinensis complex exhibits diverse phenotypic and ecological differences despite recent divergence (approximately 1.59 million years ago). To elucidate the process of genetic differentiation during early stages of ecological speciation, we analyzed genomic divergence in the Miscanthus complex using 72 randomly selected genes from a newly assembled transcriptome. In this study, rampant gene flow was detected between species, estimated as M = 3.36 × 10(-9) to 1.20 × 10(-6) , resulting in contradicting phylogenies across loci. Nevertheless, beast analyses revealed the species identity and the effects of extrinsic cohesive forces that counteracted the non-stop introgression. As expected, early in speciation with gene flow, only 3-13 loci were highly diverged; two to five outliers (approximately 2.78-6.94% of the genome) were characterized by strong linkage disequilibrium, and asymmetrically distributed among ecotypes, indicating footprints of diversifying selection. In conclusion, ecological speciation of incipient species of Miscanthus probably followed the parapatric model, whereas allopatric speciation cannot be completely ruled out, especially between the geographically isolated northern and southern M. sinensis, for which no significant gene flow across oceanic barriers was detected. Divergence between local ecotypes in early-stage speciation began at a few genomic regions under the influence of natural selection and divergence hitchhiking that overcame gene flow.

Reappraisal of known malaria resistance loci in a large multicenter study

Many human genetic associations with resistance to malaria have been reported, but few have been reliably replicated. We collected data on 11,890 cases of severe malaria due to Plasmodium falciparum and 17,441 controls from 12 locations in Africa, Asia and Oceania. We tested 55 SNPs in 27 loci previously reported to associate with severe malaria. There was evidence of association at P < 1 × 10(-4) with the HBB, ABO, ATP2B4, G6PD and CD40LG loci, but previously reported associations at 22 other loci did not replicate in the multicenter analysis. The large sample size made it possible to identify authentic genetic effects that are heterogeneous across populations or phenotypes, with a striking example being the main African form of G6PD deficiency, which reduced the risk of cerebral malaria but increased the risk of severe malarial anemia. The finding that G6PD deficiency has opposing effects on different fatal complications of P. falciparum infection indicates that the evolutionary origins of this common human genetic disorder are more complex than previously supposed.

New frontiers in the study of human cultural and genetic evolution

In this review, we discuss the dynamic linkages between culture and the genetic evolution of the human species. We begin by briefly describing the framework of gene-culture coevolutionary (or dual-inheritance) models for human evolutionary change. Until recently, the literature on gene-culture coevolution was composed primarily of mathematical models and formalized theory describing the complex dynamics underlying human behavior, adaptation, and technological evolution, but had little empirical support concerning genetics. The rapid progress in the fields of molecular genetics and genomics, however, is now providing the kinds of data needed to produce rich empirical support for gene-culture coevolutionary models. We briefly outline how theoretical and methodological progress in genome sciences has provided ways for the strength of selection on genes to be evaluated, and then outline how evidence of selection on several key genes can be directly linked to human cultural practices. We then describe some exciting new directions in the empirical study of gene-culture coevolution, and conclude with a discussion of the role of gene-culture evolutionary models in the future integration of medical, biological, and social sciences.

Natural variation in abiotic stress responsive gene expression and local adaptation to climate in Arabidopsis thaliana

Gene expression varies widely in natural populations, yet the proximate and ultimate causes of this variation are poorly known. Understanding how variation in gene expression affects abiotic stress tolerance, fitness, and adaptation is central to the field of evolutionary genetics. We tested the hypothesis that genes with natural genetic variation in their expression responses to abiotic stress are likely to be involved in local adaptation to climate in Arabidopsis thaliana. Specifically, we compared genes with consistent expression responses to environmental stress (expression stress responsive, "eSR") to genes with genetically variable responses to abiotic stress (expression genotype-by-environment interaction, "eGEI"). We found that on average genes that exhibited eGEI in response to drought or cold had greater polymorphism in promoter regions and stronger associations with climate than those of eSR genes or genomic controls. We also found that transcription factor binding sites known to respond to environmental stressors, especially abscisic acid responsive elements, showed significantly higher polymorphism in drought eGEI genes in comparison to eSR genes. By contrast, eSR genes tended to exhibit relatively greater pairwise haplotype sharing, lower promoter diversity, and fewer nonsynonymous polymorphisms, suggesting purifying selection or selective sweeps. Our results indicate that cis-regulatory evolution and genetic variation in stress responsive gene expression may be important mechanisms of local adaptation to climatic selective gradients.

Molecular evidence of Plasmodium vivax mono and mixed malaria parasite infections in Duffy-negative native Cameroonians

The malaria parasite Plasmodium vivax is known to be majorly endemic to Asian and Latin American countries with no or very few reports of Africans infected with this parasite. Since the human Duffy antigens act as receptors for P. vivax to invade human RBCs and Africans are generally Duffy-negative, non-endemicity of P. vivax in Africa has been attributed to this fact. However, recent reports describing P. vivax infections in Duffy-negative Africans from West and Central parts of Africa have been surfaced including a recent report on P. vivax infection in native Cameroonians. In order to know if Cameroonians living in the southern regions are also susceptible to P. vivax infection, we collected finger-prick blood samples from 485 malarial symptomatic patients in five locations and followed PCR diagnostic assays with DNA sequencing of the 18S ribosomal RNA gene. Out of the 201 malaria positive cases detected, 193 were pure P. falciparum, six pure P. vivax and two mixed parasite infections (P. falciparum + P. vivax). The eight P. vivax infected samples (six single + two mixed) were further subjected to DNA sequencing of the P. vivax multidrug resistance 1 (pvmdr1) and the P.vivax circumsporozoite (pvcsp) genes. Alignment of the eight Cameroonian pvmdr1 sequences with the reference sequence showed high sequence similarities, reconfirming P. vivax infection in all the eight patients. DNA sequencing of the pvcsp gene indicated all the eight P. vivax to be of VK247 type. Interestingly, DNA sequencing of a part of the human Duffy gene covering the promoter region in the eight P. vivax-infected Cameroonians to identify the T-33C mutation revealed all these patients as Duffy-negative. The results provide evidence of single P. vivax as well as mixed malaria parasite infection in native Cameroonians and add knowledge to the growing evidences of P. vivax infection in Duffy-negative Africans.

Frequencies of Duffy blood group alleles in Northern Pakistani donors

OBJECTIVE

Find the allele frequencies of Duffy blood group antigens in donor population from northern Pakistan.

DESIGN

Cross sectional study.

PLACE AND DURATION OF STUDY

Armed Forces Institute of Transfusion (AFIT), Rawalpindi in year 2012.

PATIENTS AND METHODS

A total of 1000 healthy, adult blood donors were included in the study. Blood samples were collected in ethylenediamine tetra aceticacid (EDTA) tube and then tested with anti sera Fy(a) and Fy(b) by the tube method.

RESULTS

The allele frequencies of Duffy blood group antigens were calculated. The most common phenotype was Fy(a+b+) which was present in 552 (55.2%) donors followed by the Fy(a+b-) phenotype in 228 (22.8%) donors, while 178 (17.8%) were Fy(a-b+) and the least prevalent phenotype was Fy(a-b-) which was present in 42 (4.2%) of donors.

CONCLUSION

The majority of our population is heterozygous for Duffy antigens a and b.

Deciphering the genetic blueprint behind Holstein milk proteins and production

Holstein is known to provide higher milk yields than most other cattle breeds, and the dominant position of Holstein today is the result of various selection pressures. Holstein cattle have undergone intensive selection for milk production in recent decades, which has left genome-wide footprints of domestication. To further characterize the bovine genome, we performed whole-genome resequencing analysis of 10 Holstein and 11 Hanwoo cattle to identify regions containing genes as outliers in Holstein, including CSN1S1, CSN2, CSN3, and KIT whose products are likely involved in the yield and proteins of milk and their distinctive black-and-white markings. In addition, genes indicative of positive selection were associated with cardiovascular disease, which is related to simultaneous propagation of genetic defects, also known as inbreeding depression in Holstein.

Extreme population differences in the human zinc transporter ZIP4 (SLC39A4) are explained by positive selection in Sub-Saharan Africa

Extreme differences in allele frequency between West Africans and Eurasians were observed for a leucine-to-valine substitution (Leu372Val) in the human intestinal zinc uptake transporter, ZIP4, yet no further evidence was found for a selective sweep around the ZIP4 gene (SLC39A4). By interrogating allele frequencies in more than 100 diverse human populations and resequencing Neanderthal DNA, we confirmed the ancestral state of this locus and found a strong geographical gradient for the derived allele (Val372), with near fixation in West Africa. In extensive coalescent simulations, we show that the extreme differences in allele frequency, yet absence of a classical sweep signature, can be explained by the effect of a local recombination hotspot, together with directional selection favoring the Val372 allele in Sub-Saharan Africans. The possible functional effect of the Leu372Val substitution, together with two pathological mutations at the same codon (Leu372Pro and Leu372Arg) that cause acrodermatitis enteropathica (a disease phenotype characterized by extreme zinc deficiency), was investigated by transient overexpression of human ZIP4 protein in HeLa cells. Both acrodermatitis mutations cause absence of the ZIP4 transporter cell surface expression and nearly absent zinc uptake, while the Val372 variant displayed significantly reduced surface protein expression, reduced basal levels of intracellular zinc, and reduced zinc uptake in comparison with the Leu372 variant. We speculate that reduced zinc uptake by the ZIP4-derived Val372 isoform may act by starving certain pathogens of zinc, and hence may have been advantageous in Sub-Saharan Africa. Moreover, these functional results may indicate differences in zinc homeostasis among modern human populations with possible relevance for disease risk.

Signatures of natural selection on genetic variants affecting complex human traits

It has recently been hypothesized that polygenic adaptation, resulting in modest allele frequency changes at many loci, could be a major mechanism behind the adaptation of complex phenotypes in human populations. Here we leverage the large number of variants that have been identified through genome-wide association (GWA) studies to comprehensively study signatures of natural selection on genetic variants associated with complex traits. Using population differentiation based methods, such as FST and phylogenetic branch length analyses, we systematically examined nearly 1300 SNPs associated with 38 complex phenotypes. Instead of detecting selection signatures at individual variants, we aimed to identify combined evidence of natural selection by aggregating signals across many trait associated SNPs. Our results have revealed some general features of polygenic selection on complex traits associated variants. First, natural selection acting on standing variants associated with complex traits is a common phenomenon. Second, characteristics of selection for different polygenic traits vary both temporarily and geographically. Third, some studied traits (e.g. height and urate level) could have been the primary targets of selection, as indicated by the significant correlation between the effect sizes and the estimated strength of selection in the trait associated variants; however, for most traits, the allele frequency changes in trait associated variants might have been driven by the selection on other correlated phenotypes. Fourth, the changes in allele frequencies as a result of selection can be highly stochastic, such that, polygenic adaptation may accelerate differentiation in allele frequencies among populations, but generally does not produce predictable directional changes. Fifth, multiple mechanisms (pleiotropy, hitchhiking, etc) may act together to govern the changes in allele frequencies of genetic variants associated with complex traits.

NAT2 sequence polymorphisms and acetylation profiles in Indians

BACKGROUND

NAT2, a broad-spectrum drug-metabolizing gene, is of high pharmacogenetic interest. Based on seven different mutations in the NAT2 gene, an individual can either be categorized as a slow or fast acetylator.

MATERIALS & METHODS

In order to characterize acetylation profiles of Indians, where data are poorly available, we sequenced the 873 bp NAT2 coding region in 250 Indians, covering the whole of India including three tribes.

RESULTS

Altogether, 35 NAT2 alleles forming two acetylator phenotypes (distributed almost in equal proportion in India) were found; while the alleles determining slow acetylators were highly differentiated, the fast acetylator alleles were less in number but highly frequent.

CONCLUSION

Interestingly, distribution of two different acetylation phenotypes correlated well with historical dietary pattern in India. The neighbor-joining phylogenetic tree based on NAT2 gene polymorphisms in worldwide humans revealed genetic affinities among populations with similar acetylation phenotypes, which also placed Indians and Africans together in a single cluster.

Selection and adaptation in the human genome

An enduring goal of evolutionary biology is to understand how natural selection has shaped patterns of polymorphism and divergence within and between species and to map the genetic basis of adaptations. The rapid maturation of next-generation sequencing technology has generated a deluge of genomics data from nonhuman primates, extinct hominins, and diverse human populations. These emerging genome data sets have simultaneously broadened our understanding of human evolution and sharply defined existing gaps in knowledge about the mechanistic basis of evolutionary change. In this review, we summarize recent insights into how natural selection has influenced the human genome across different timescales. Although the path to a more comprehensive understanding of selection and adaptation in humans remains arduous, some general insights are beginning to emerge, such as the importance of adaptive regulatory evolution, the absence of pervasive classic selective sweeps, and the potential roles that selection from standing variation and polygenic adaptation have likely played in recent human evolutionary history.

Patterns of nucleotide and haplotype diversity at ICAM-1 across global human populations with varying levels of malaria exposure

Malaria is one of the strongest selective pressures in recent human evolution. African populations have been and continue to be at risk for malarial infections. However, few studies have re-sequenced malaria susceptibility loci across geographically and genetically diverse groups in Africa. We examined nucleotide diversity at Intercellular adhesion molecule-1 (ICAM-1), a malaria susceptibility candidate locus, in a number of human populations with a specific focus on diverse African ethnic groups. We used tests of neutrality to assess whether natural selection has impacted this locus and tested whether SNP variation at ICAM-1 is correlated with malaria endemicity. We observe differing patterns of nucleotide and haplotype variation in global populations and higher levels of diversity in Africa. Although we do not observe a deviation from neutrality based on the allele frequency distribution, we do observe several alleles at ICAM-1, including the ICAM-1 (Kilifi) allele, that are correlated with malaria endemicity. We show that the ICAM-1 (Kilifi) allele, which is common in Africa and Asia, exists on distinct haplotype backgrounds and is likely to have arisen more recently in Asia. Our results suggest that correlation analyses of allele frequencies and malaria endemicity may be useful for identifying candidate functional variants that play a role in malaria resistance and susceptibility.

Population genetic tools for dissecting innate immunity in humans

Innate immunity involves direct interactions between the host and microorganisms, both pathogenic and symbiotic, so natural selection is expected to strongly influence genes involved in these processes. Population genetics investigates the impact of past natural selection events on the genome of present-day human populations, and it complements immunological as well as clinical and epidemiological genetic studies. Recent data show that the impact of selection on the different families of innate immune receptors and their downstream signalling molecules varies considerably. This Review discusses these findings and highlights how they help to delineate the relative functional importance of innate immune pathways, which can range from being essential to being redundant.

Genomics, population genetics and evolutionary history of Plasmodium vivax

Plasmodium vivax is part of a highly diverse clade that includes several Plasmodium species found in nonhuman primates from Southeast Asia. The diversity of primate malarias in Asia is staggering; nevertheless, their origin was relatively recent in the evolution of Plasmodium. We discuss how humans acquired the lineage leading to P. vivax from a nonhuman primate determined by the complex geological processes that took place in Southeast Asia during the last few million years. We conclude that widespread population genomic investigations are needed in order to understand the demographic processes involved in the expansion of P. vivax in the human populations. India represents one of the few countries with widespread vivax malaria. Earlier studies have indicated high genetic polymorphism at antigenic loci and no evidence for geographic structuring. However, new studies using genetic markers in selectively neutral genetic regions indicate that Indian P. vivax presents complex evolutionary history but possesses features consistent with being part of the ancestral distribution range of this species. Such studies are possible due to the availability of the first P. vivax genome sequences. Next generation sequencing technologies are now paving the way for the sequencing of more P. vivax genomes that will dramatically increase our understanding of the unique biology of this species.

Pathogen-driven selection in the human genome

Infectious diseases and epidemics have always accompanied and characterized human history, representing one of the main causes of death. Even today, despite progress in sanitation and medical research, infections are estimated to account for about 15% of deaths. The hypothesis whereby infectious diseases have been acting as a powerful selective pressure was formulated long ago, but it was not until the availability of large-scale genetic data and the development of novel methods to study molecular evolution that we could assess how pervasively infectious agents have shaped human genetic diversity. Indeed, recent evidences indicated that among the diverse environmental factors that acted as selective pressures during the evolution of our species, pathogen load had the strongest influence. Beside the textbook example of the major histocompatibility complex, selection signatures left by pathogen-exerted pressure can be identified at several human loci, including genes not directly involved in immune response. In the future, high-throughput technologies and the availability of genetic data from different populations are likely to provide novel insights into the evolutionary relationships between the human host and its pathogens. Hopefully, this will help identify the genetic determinants modulating the susceptibility to infectious diseases and will translate into new treatment strategies.

Personal and population genomics of human regulatory variation

The characteristics and evolutionary forces acting on regulatory variation in humans remains elusive because of the difficulty in defining functionally important noncoding DNA. Here, we combine genome-scale maps of regulatory DNA marked by DNase I hypersensitive sites (DHSs) from 138 cell and tissue types with whole-genome sequences of 53 geographically diverse individuals in order to better delimit the patterns of regulatory variation in humans. We estimate that individuals likely harbor many more functionally important variants in regulatory DNA compared with protein-coding regions, although they are likely to have, on average, smaller effect sizes. Moreover, we demonstrate that there is significant heterogeneity in the level of functional constraint in regulatory DNA among different cell types. We also find marked variability in functional constraint among transcription factor motifs in regulatory DNA, with sequence motifs for major developmental regulators, such as HOX proteins, exhibiting levels of constraint comparable to protein-coding regions. Finally, we perform a genome-wide scan of recent positive selection and identify hundreds of novel substrates of adaptive regulatory evolution that are enriched for biologically interesting pathways such as melanogenesis and adipocytokine signaling. These data and results provide new insights into patterns of regulatory variation in individuals and populations and demonstrate that a large proportion of functionally important variation lies beyond the exome.

The impact of natural selection on health and disease: uses of the population genetics approach in humans

Investigations of the legacy of natural selection in the human genome have proved particularly informative, pinpointing functionally important regions that have participated in our genetic adaptation to the environment. Furthermore, genetic dissection of the intensity and type of selection acting on human genes can be used to predict involvement in different forms and severities of human diseases. We review here the progress made in population genetics studies toward understanding the effects of selection, in its different forms and intensities, on human genome diversity. We discuss some outstanding, robust examples of genes and biological functions subject to strong dietary, climatic and pathogen selection pressures. We also explore the possible relationship between cancer and natural selection, a topic that has been largely neglected because cancer is generally seen as a late-onset disease. Finally, we discuss how the present-day incidence of some diseases of modern societies may represent a by-product of past adaptation to other selective forces and changes in lifestyle. This perspective thus illustrates the value of adopting a population genetics approach in delineating the biological mechanisms that have played a major evolutionary role in the way humans have genetically adapted to different environments and lifestyles over time.

Duffy Blood Group System and the malaria adaptation process in humans

Malaria is an acute infectious disease caused by the protozoa of the genus Plasmodium. The antigens of the Duffy Blood Group System, in addition to incompatibilities in transfusions and hemolytic disease of the newborn, are of great interest in medicine due to their association with the invasion of red blood cells by the parasite Plasmodium vivax. For invasions to occur an interaction between the parasites and antigens of the Duffy Blood Group System is necessary. In Caucasians six antigens are produced by the Duffy locus (Fya, Fyb, F3, F4, F5 and F6). It has been observed that Fy(a-b-) individuals are resistant to Plasmodium knowlesi and P. vivax infection, because the invasion requires at least one of these antigens. The P. vivax Duffy Binding Protein (PvDBP) is functionally important in the invasion process of these parasites in Duffy / DARC positive humans. The proteins or fractions may be considered, therefore, an important and potential inoculum to be used in immunization against malaria.

Primates, Pathogens, and Evolution: A Context for Understanding Emerging Disease

The world is rife with potential pathogens. Of those that infect humans, it is estimated that roughly 20 % are of nonhuman primate origin. The same ease characterizes pathogen transmission in the other direction, from humans to nonhuman primates. This latter problem has increasingly serious ramifications for conservation efforts, as growing numbers of ecotourists and researchers serve as potential vectors of disease. Here, we present an analysis of major cross-species transmission events between human and nonhuman primates. In particular, we consider HIV and malaria as case studies in which nonhuman primate pathogens emerged and became permanent fixtures in human populations. The human practices that facilitate such events are considered, as well as the evolutionary consequences of these events. In addition, we describe human-to-nonhuman primate transmission events and discuss the potential of human pathogens to adapt to nonhuman primate hosts. The topic of emerging infections is addressed, in both human and nonhuman species, in light of changing patterns of contact and novel adaptations on the part of pathogens and hosts.

Red blood cell polymorphism and susceptibility to Plasmodium vivax

Resistance to Plasmodium vivax blood-stage infection has been widely recognised to result from absence of the Duffy (Fy) blood group from the surface of red blood cells (RBCs) in individuals of African descent. Interestingly, recent studies from different malaria-endemic regions have begun to reveal new perspectives on the association between Duffy gene polymorphism and P. vivax malaria. In Papua New Guinea and the Americas, heterozygous carriers of a Duffy-negative allele are less susceptible to P. vivax infection than Duffy-positive homozygotes. In Brazil, studies show that the Fy(a) antigen, compared to Fy(b), is associated with lower binding to the P. vivax Duffy-binding protein and reduced susceptibility to vivax malaria. Additionally, it is interesting that numerous studies have now shown that P. vivax can infect RBCs and cause clinical disease in Duffy-negative people. This suggests that the relationship between P. vivax and the Duffy antigen is more complex than customarily described. Evidence of P. vivax Duffy-independent red cell invasion indicates that the parasite must be evolving alternative red cell invasion pathways. In this chapter, we review the evidence for P. vivax Duffy-dependent and Duffy-independent red cell invasion. We also consider the influence of further host gene polymorphism associated with malaria endemicity on susceptibility to vivax malaria. The interaction between the parasite and the RBC has significant potential to influence the effectiveness of P. vivax-specific vaccines and drug treatments. Ultimately, the relationships between red cell polymorphisms and P. vivax blood-stage infection will influence our estimates on the population at risk and efforts to eliminate vivax malaria.

Natural selection mediated association of the Duffy (FY) gene polymorphisms with Plasmodium vivax malaria in India

The Duffy (Fy) antigens act as receptors for chemokines as well as for Plasmodium vivax to invade human RBCs. A recent study has correlated the occurrence of the FY*A allele of Duffy gene with decreased susceptibility to vivax malaria, but no epidemiological correlation between the distribution of FY*A allele and incidences of vivax malaria has been established so far. Furthermore, if such correlations exist, whether natural selection has mediated the association, is an important question. Since India is highly endemic to P. vivax malaria with variable eco-climatic and varying vivax malaria epidemiology across different regions, such a question could well be answered in Indians. For this, we have genotyped the FY gene at the −33^rd and the 125^th nucleotide positions in 250 Indians sampled from six different zonal plus one tribal population covering the whole of India and studied possible correlations with eco-climatic and vivax malaria incidences. No FY*O allele was found, however, both the FY*A and FY*B alleles forming FY*A/FY*A, FY*A/FY*B and FY*B/FY*B genotypes were widely distributed among Indians. Five out of seven population samples significantly deviated from the Hardy-Weinberg equilibrium expectation, and two alleles (FY*A and FY*B) and the homozygote genotype, FY*B/FY*B were clinally distributed over the population coordinates. Furthermore, vivax malaria incidences over the past five years were significantly negatively and positively associated with the frequencies of the FY*A and FY*B alleles, respectively. The Northern Indians were highly differentiated from the other zonal population samples at the FY gene, as evidenced from the reconstructed Neighbor-Joining phylogenetic tree. The results specify the role of natural selection in the distribution of FY gene polymorphism in India. Furthermore, the hypotheses on the part of the FY*A allele in conferring protection to vivax malaria could be validated following population genetic studies in a vivax malaria epidemiological setting, such as India.

Genome-wide association and population genetic analysis of C-reactive protein in African American and Hispanic American women

C-reactive protein (CRP) is a systemic inflammation marker that predicts future cardiovascular risk. CRP levels are higher in African Americans and Hispanic Americans than in European Americans, but the genetic determinants of CRP in these admixed United States minority populations are largely unknown. We performed genome-wide association studies (GWASs) of 8,280 African American (AA) and 3,548 Hispanic American (HA) postmenopausal women from the Women's Health Initiative SNP Health Association Resource. We discovered and validated a CRP-associated variant of triggering receptors expressed by myeloid cells 2 (TREM2) in chromosomal region 6p21 (p = 10(-10)). The TREM2 variant associated with higher CRP is common in Africa but rare in other ancestral populations. In AA women, the CRP region in 1q23 contained a strong admixture association signal (p = 10(-17)), which appears to be related to several independent CRP-associated alleles; the strongest of these is present only in African ancestral populations and is associated with higher CRP. Of the other genomic loci previously associated with CRP through GWASs of European populations, most loci (LEPR, IL1RN, IL6R, GCKR, NLRP3, HNF1A, HNF4A, and APOC1) showed consistent patterns of association with CRP in AA and HA women. In summary, we have identified a common TREM2 variant associated with CRP in United States minority populations. The genetic architecture underlying the CRP phenotype in AA women is complex and involves genetic variants shared across populations, as well as variants specific to populations of African descent.

Evolutionary forces shaping genomic islands of population differentiation in humans

BACKGROUND

Levels of differentiation among populations depend both on demographic and selective factors: genetic drift and local adaptation increase population differentiation, which is eroded by gene flow and balancing selection. We describe here the genomic distribution and the properties of genomic regions with unusually high and low levels of population differentiation in humans to assess the influence of selective and neutral processes on human genetic structure.

METHODS

Individual SNPs of the Human Genome Diversity Panel (HGDP) showing significantly high or low levels of population differentiation were detected under a hierarchical-island model (HIM). A Hidden Markov Model allowed us to detect genomic regions or islands of high or low population differentiation.

RESULTS

Under the HIM, only 1.5% of all SNPs are significant at the 1% level, but their genomic spatial distribution is significantly non-random. We find evidence that local adaptation shaped high-differentiation islands, as they are enriched for non-synonymous SNPs and overlap with previously identified candidate regions for positive selection. Moreover there is a negative relationship between the size of islands and recombination rate, which is stronger for islands overlapping with genes. Gene ontology analysis supports the role of diet as a major selective pressure in those highly differentiated islands. Low-differentiation islands are also enriched for non-synonymous SNPs, and contain an overly high proportion of genes belonging to the 'Oncogenesis' biological process.

CONCLUSIONS

Even though selection seems to be acting in shaping islands of high population differentiation, neutral demographic processes might have promoted the appearance of some genomic islands since i) as much as 20% of islands are in non-genic regions ii) these non-genic islands are on average two times shorter than genic islands, suggesting a more rapid erosion by recombination, and iii) most loci are strongly differentiated between Africans and non-Africans, a result consistent with known human demographic history.

Molecular evolution of a malaria resistance gene (DARC) in primates

Genes involved in host-pathogen interactions are often strongly affected by positive natural selection. The Duffy antigen, coded by the Duffy antigen receptor for chemokines (DARC) gene, serves as a receptor for Plasmodium vivax in humans and for Plasmodium knowlesi in some nonhuman primates. In the majority of sub-Saharan Africans, a nucleic acid variant in GATA-1 of the gene promoter is responsible for the nonexpression of the Duffy antigen on red blood cells and consequently resistance to invasion by P. vivax. The Duffy antigen also acts as a receptor for chemokines and is expressed in red blood cells and many other tissues of the body. Because of this dual role, we sequenced a ~3,000-bp region encompassing the entire DARC gene as well as part of its 5' and 3' flanking regions in a phylogenetic sample of primates and used statistical methods to evaluate the nature of selection pressures acting on the gene during its evolution. We analyzed both coding and regulatory regions of the DARC gene. The regulatory analysis showed accelerated rates of substitution at several sites near known motifs. Our tests of positive selection in the coding region using maximum likelihood by branch sites and maximum likelihood by codon sites did not yield statistically significant evidence for the action of positive selection. However, the maximum likelihood test in which the gene was subdivided into different structural regions showed that the known binding region for P. vivax/P. knowlesi is under very different selective pressures than the remainder of the gene. In fact, most of the gene appears to be under strong purifying selection, but this is not evident in the binding region. We suggest that the binding region is under the influence of two opposing selective pressures, positive selection possibly exerted by the parasite and purifying selection exerted by chemokines.

Signatures of positive selection in Toll-like receptor (TLR) genes in mammals

Background

Toll-like receptors (TLRs) are a major class of pattern recognition receptors (PRRs) expressed in the cell surface or membrane compartments of immune and non-immune cells. TLRs are encoded by a multigene family and represent the first line of defense against pathogens by detecting foreigner microbial molecular motifs, the pathogen-associated molecular patterns (PAMPs). TLRs are also important by triggering the adaptive immunity in vertebrates. They are characterized by the presence of leucine-rich repeats (LRRs) in the ectodomain, which are associated with the PAMPs recognition. The direct recognition of different pathogens by TLRs might result in different evolutionary adaptations important to understand the dynamics of the host-pathogen interplay. Ten mammal TLR genes, viral (TLR3, 7, 8, 9) and non-viral (TLR1-6, 10), were selected to identify signatures of positive selection that might have been imposed by interacting pathogens and to clarify if viral and non-viral TLRs might display different patterns of molecular evolution.

Results

By using Maximum Likelihood approaches, evidence of positive selection was found in all the TLRs studied. The number of positively selected codons (PSC) ranged between 2-26 codons (0.25%-2.65%) with the non-viral TLR4 as the receptor with higher percentage of positively selected codons (2.65%), followed by the viral TLR8 (2.50%). The results indicated that viral and non-viral TLRs are similarly under positive selection. Almost all TLRs have at least one PSC located in the LRR ectodomain which underlies the importance of the pathogen recognition by this region.

Conclusions

Our results are not in line with previous studies on primates and birds that identified more codons under positive selection in non-viral TLRs. This might be explained by the fact that both primates and birds are homogeneous groups probably being affected by only a restricted number of related viruses with equivalent motifs to be recognized. The analyses performed in this work encompassed a large number of species covering some of the most representative mammalian groups - Artiodactyla, Rodents, Carnivores, Lagomorphs and Primates - that are affected by different families of viruses. This might explain the role of adaptive evolution in shaping viral TLR genes.

Exploration of signals of positive selection derived from genotype-based human genome scans using re-sequencing data

We have investigated whether regions of the genome showing signs of positive selection in scans based on haplotype structure also show evidence of positive selection when sequence-based tests are applied, whether the target of selection can be localized more precisely, and whether such extra evidence can lead to increased biological insights. We used two tools: simulations under neutrality or selection, and experimental investigation of two regions identified by the HapMap2 project as putatively selected in human populations. Simulations suggested that neutral and selected regions should be readily distinguished and that it should be possible to localize the selected variant to within 40 kb at least half of the time. Re-sequencing of two ~300 kb regions (chr4:158Mb and chr10:22Mb) lacking known targets of selection in HapMap CHB individuals provided strong evidence for positive selection within each and suggested the micro-RNA gene hsa-miR-548c as the best candidate target in one region, and changes in regulation of the sperm protein gene SPAG6 in the other.

The global distribution of the Duffy blood group

Blood group variants are characteristic of population groups, and can show conspicuous geographic patterns. Interest in the global prevalence of the Duffy blood group variants is multidisciplinary, but of particular importance to malariologists due to the resistance generally conferred by the Duffy-negative phenotype against Plasmodium vivax infection. Here we collate an extensive geo-database of surveys, forming the evidence-base for a multi-locus Bayesian geostatistical model to generate global frequency maps of the common Duffy alleles to refine the global cartography of the common Duffy variants. We show that the most prevalent allele globally was FY*A, while across sub-Saharan Africa the predominant allele was the silent FY*B(ES) variant, commonly reaching fixation across stretches of the continent. The maps presented not only represent the first spatially and genetically comprehensive description of variation at this locus, but also constitute an advance towards understanding the transmission patterns of the neglected P. vivax malaria parasite.

Human TLRs and IL-1Rs in host defense: natural insights from evolutionary, epidemiological, and clinical genetics

Toll-like receptors (TLRs) and interleukin-1 receptors (IL-1Rs) have TIR intracellular domains that engage two main signaling pathways, via the TIR-containing adaptors MyD88 (which is not used by TLR3) and TRIF (which is used only by TLR3 and TLR4). Extensive studies in inbred mice in various experimental settings have attributed key roles in immunity to TLR- and IL-1R-mediated responses, but what contribution do human TLRs and IL-1Rs actually make to host defense in the natural setting? Evolutionary genetic studies have shown that human intracellular TLRs have evolved under stronger purifying selection than surface-expressed TLRs, for which the frequency of missense and nonsense alleles is high in the general population. Epidemiological genetic studies have yet to provide convincing evidence of a major contribution of common variants of human TLRs, IL-1Rs, or their adaptors to host defense. Clinical genetic studies have revealed that rare mutations affecting the TLR3-TRIF pathway underlie herpes simplex virus encephalitis, whereas mutations in the TIR-MyD88 pathway underlie pyogenic bacterial diseases in childhood. A careful reconsideration of the contributions of TLRs and IL-1Rs to host defense in natura is required.

Population genetics of malaria resistance in humans

The high mortality and widespread impact of malaria have resulted in this disease being the strongest evolutionary selective force in recent human history, and genes that confer resistance to malaria provide some of the best-known case studies of strong positive selection in modern humans. I begin by reviewing JBS Haldane's initial contribution to the potential of malaria genetic resistance in humans. Further, I discuss the population genetics aspects of many of the variants, including globin, G6PD deficiency, Duffy, ovalocytosis, ABO and human leukocyte antigen variants. Many of the variants conferring resistance to malaria are 'loss-of-function' mutants and appear to be recent polymorphisms from the last 5000-10 000 years or less. I discuss estimation of selection coefficients from case-control data and make predictions about the change for S, C and G6PD-deficiency variants. In addition, I consider the predicted joint changes when the two β-globin alleles S and C are both variable in the same population and when there is a variation for α-thalassemia and S, two unlinked, but epistatic variants. As more becomes known about genes conferring genetic resistance to malaria in humans, population genetics approaches can contribute both to investigating past selection and predicting the consequences in future generations for these variants.

Integrating evolutionary and functional approaches to infer adaptation at specific loci

Inferences about adaptation at specific loci are often exclusively based on the static analysis of DNA sequence variation. Ideally,population-genetic evidence for positive selection serves as a stepping-off point for experimental studies to elucidate the functional significance of the putatively adaptive variation. We argue that inferences about adaptation at specific loci are best achieved by integrating the indirect, retrospective insights provided by population-genetic analyses with the more direct, mechanistic insights provided by functional experiments. Integrative studies of adaptive genetic variation may sometimes be motivated by experimental insights into molecular function, which then provide the impetus to perform population genetic tests to evaluate whether the functional variation is of adaptive significance. In other cases, studies may be initiated by genome scans of DNA variation to identify candidate loci for recent adaptation. Results of such analyses can then motivate experimental efforts to test whether the identified candidate loci do in fact contribute to functional variation in some fitness-related phenotype. Functional studies can provide corroborative evidence for positive selection at particular loci, and can potentially reveal specific molecular mechanisms of adaptation.

The role played by natural selection on Mendelian traits in humans

Detecting whether and how natural selection has targeted regions of the human genome represents a complementary strategy for identifying functionally important loci and variants involved in disease resistance and adaptation to the environment. In contrast with most complex diseases or traits, the genetic architecture of most Mendelian traits is relatively well established. Most mutations associated with Mendelian disease-related traits are highly penetrant and kept at low population frequencies because of the effects of purifying selection. However, this is not always the case. Here, we review several examples of Mendelian mutations-associated with various disease conditions or other traits of anthropological interest-that have increased in frequency in the human population as a result of past positive selection. These examples clearly illustrate the value of a population genetics approach to unravel the biological mechanisms that have been central to our past and present survival against the selective pressures imposed by diseases and other environmental factors.

Nonadaptive processes in primate and human evolution

Evolutionary biology has tended to focus on adaptive evolution by positive selection as the primum mobile of evolutionary trajectories in species while underestimating the importance of nonadaptive evolutionary processes. In this review, I describe evidence that suggests that primate and human evolution has been strongly influenced by nonadaptive processes, particularly random genetic drift and mutation. This is evidenced by three fundamental effects: a relative relaxation of selective constraints (i.e., purifying selection), a relative increase in the fixation of slightly deleterious mutations, and a general reduction in the efficacy of positive selection. These effects are observed in protein-coding, regulatory regions, and in gene expression data, as well as in an augmentation of fixation of large-scale mutations, including duplicated genes, mobile genetic elements, and nuclear mitochondrial DNA. The evidence suggests a general population-level explanation such as a reduction in effective population size (N(e)). This would have tipped the balance between the evolutionary forces of natural selection and random genetic drift toward genetic drift for variants having small selective effects. After describing these proximate effects, I describe the potential consequences of these effects for primate and human evolution. For example, an increase in the fixation of slightly deleterious mutations could potentially have led to an increase in the fixation rate of compensatory mutations that act to suppress the effects of slightly deleterious substitutions. The potential consequences of compensatory evolution for the evolution of novel gene functions and in potentially confounding the detection of positively selected genes are explored. The consequences of the passive accumulation of large-scale genomic mutations by genetic drift are unclear, though evidence suggests that new gene copies as well as insertions of transposable elements into genes can potentially lead to adaptive phenotypes. Finally, because a decrease in selective constraint at the genetic level is expected to have effects at the morphological level, I review studies that compare rates of morphological change in various mammalian and island populations where N(e) is reduced. Furthermore, I discuss evidence that suggests that craniofacial morphology in the Homo lineage has shifted from an evolutionary rate constrained by purifying selection toward a neutral evolutionary rate.

Allele-specific gene expression in a wild nonhuman primate population

Natural populations hold enormous potential for evolutionary genetic studies, especially when phenotypic, genetic and environmental data are all available on the same individuals. However, untangling the genotype-phenotype relationship in natural populations remains a major challenge. Here, we describe results of an investigation of one class of phenotype, allele-specific gene expression (ASGE), in the well-studied natural population of baboons of the Amboseli basin, Kenya. ASGE measurements identify cases in which one allele of a gene is overexpressed relative to the alternative allele of the same gene, within individuals, thus providing a control for background genetic and environmental effects. Here, we characterize the incidence of ASGE in the Amboseli baboon population, focusing on the genetic and environmental contributions to ASGE in a set of eleven genes involved in immunity and defence. Within this set, we identify evidence for common ASGE in four genes. We also present examples of two relationships between cis-regulatory genetic variants and the ASGE phenotype. Finally, we identify one case in which this relationship is influenced by a novel gene-environment interaction. Specifically, the dominance rank of an individual's mother during its early life (an aspect of that individual's social environment) influences the expression of the gene CCL5 via an interaction with cis-regulatory genetic variation. These results illustrate how environmental and ecological data can be integrated into evolutionary genetic studies of functional variation in natural populations. They also highlight the potential importance of early life environmental variation in shaping the genetic architecture of complex traits in wild mammals.

Genome-wide association study SNPs in the human genome diversity project populations: does selection affect unlinked SNPs with shared trait associations?

Genome-wide association studies (GWAS) have identified more than 2,000 trait-SNP associations, and the number continues to increase. GWAS have focused on traits with potential consequences for human fitness, including many immunological, metabolic, cardiovascular, and behavioral phenotypes. Given the polygenic nature of complex traits, selection may exert its influence on them by altering allele frequencies at many associated loci, a possibility which has yet to be explored empirically. Here we use 38 different measures of allele frequency variation and 8 iHS scores to characterize over 1,300 GWAS SNPs in 53 globally distributed human populations. We apply these same techniques to evaluate SNPs grouped by trait association. We find that groups of SNPs associated with pigmentation, blood pressure, infectious disease, and autoimmune disease traits exhibit unusual allele frequency patterns and elevated iHS scores in certain geographical locations. We also find that GWAS SNPs have generally elevated scores for measures of allele frequency variation and for iHS in Eurasia and East Asia. Overall, we believe that our results provide evidence for selection on several complex traits that has caused changes in allele frequencies and/or elevated iHS scores at a number of associated loci. Since GWAS SNPs collectively exhibit elevated allele frequency measures and iHS scores, selection on complex traits may be quite widespread. Our findings are most consistent with this selection being either positive or negative, although the relative contributions of the two are difficult to discern. Our results also suggest that trait-SNP associations identified in Eurasian samples may not be present in Africa, Oceania, and the Americas, possibly due to differences in linkage disequilibrium patterns. This observation suggests that non-Eurasian and non-East Asian sample populations should be included in future GWAS.

A hierarchical Bayesian model for next-generation population genomics

The demography of populations and natural selection shape genetic variation across the genome and understanding the genomic consequences of these evolutionary processes is a fundamental aim of population genetics. We have developed a hierarchical Bayesian model to quantify genome-wide population structure and identify candidate genetic regions affected by selection. This model improves on existing methods by accounting for stochastic sampling of sequences inherent in next-generation sequencing (with pooled or indexed individual samples) and by incorporating genetic distances among haplotypes in measures of genetic differentiation. Using simulations we demonstrate that this model has a low false-positive rate for classifying neutral genetic regions as selected genes (i.e., Φ(ST) outliers), but can detect recent selective sweeps, particularly when genetic regions in multiple populations are affected by selection. Nonetheless, selection affecting just a single population was difficult to detect and resulted in a high false-negative rate under certain conditions. We applied the Bayesian model to two large sets of human population genetic data. We found evidence of widespread positive and balancing selection among worldwide human populations, including many genetic regions previously thought to be under selection. Additionally, we identified novel candidate genes for selection, several of which have been linked to human diseases. This model will facilitate the population genetic analysis of a wide range of organisms on the basis of next-generation sequence data.

Population structure and genetic differentiation associated with breeding history and selection in tomato (Solanum lycopersicum L.)

Tomato (Solanum lycopersicum L.) has undergone intensive selection during and following domestication. We investigated population structure and genetic differentiation within a collection of 70 tomato lines representing contemporary (processing and fresh-market) varieties, vintage varieties and landraces. The model-based Bayesian clustering software, STRUCTURE, was used to detect subpopulations. Six independent analyses were conducted using all marker data (173 markers) and five subsets of markers based on marker type (single-nucleotide polymorphisms, simple sequence repeats and insertion/deletions) and location (exon and intron sequences) within genes. All of these analyses consistently separated four groups predefined by market niche and age into distinct subpopulations. Furthermore, we detected at least two subpopulations within the processing varieties. These subpopulations correspond to historical patterns of breeding conducted for specific production environments. We found no subpopulation within fresh-market varieties, vintage varieties and landraces when using all marker data. High levels of admixture were shown in several varieties representing a transition in the demarcation between processing and fresh-market breeding. The genetic clustering detected by using the STRUCTURE software was confirmed by two statistics, pairwise F(st) (θ) and Nei's standard genetic distance. We also identified a total of 19 loci under positive selection between processing, fresh-market and vintage germplasm by using an F(st)-outlier method based on the deviation from the expected distribution of F(st) and heterozygosity. The markers and genome locations we identified are consistent with known patterns of selection and linkage to traits that differentiate the market classes. These results demonstrate how human selection through breeding has shaped genetic variation within cultivated tomato.

Acquired antibody responses against Plasmodium vivax infection vary with host genotype for duffy antigen receptor for chemokines (DARC)

Background

Polymorphism of the Duffy Antigen Receptor for Chemokines (DARC) is associated with susceptibility to and the severity of Plasmodium vivax malaria in humans. P. vivax uses DARC to invade erythrocytes. Individuals lacking DARC are ‘resistant’ to P. vivax erythrocytic infection. However, susceptibility to P. vivax in DARC+ individuals is reported to vary between specific DARC genotypes. We hypothesized that the natural acquisition of antibodies to P. vivax blood stages may vary with the host genotype and the level of DARC expression. Furthermore, high parasitemia has been reported to effect the acquisition of immunity against pre-erythrocytic parasites. We investigated the correlation between host DARC genotypes and the frequency and magnitude of antibodies against P. vivax erythrocytic stage antigens.

Methodology/Findings

We assessed the frequencies and magnitudes of antibody responses against P. vivax and P. falciparum sporozoite and erythrocytic antigens in Colombian donors from malaria-endemic regions. The frequency and level of naturally-acquired antibodies against the P. vivax erythrocytic antigens merozoite surface protein 1 (PvMSP1) and Duffy binding protein (PvDBP) varied with the host DARC genotypes. Donors with one negative allele (FY*B/FY*Bnull and FY*A/FY*Bnull) were more likely to have anti-PvMSP1 and anti-PvDBP antibodies than those with two positive alleles (FY*B/FY*B and FY*A/FY*B). The lower IgG3 and IgG1 components of the total IgG response may account for the decreased responses to P. vivax erythrocytic antigens with FY*A/FY*B and FY*B/FY*B genotypes. No such association was detected with P. falciparum erythrocytic antigens, which does not use DARC for erythrocyte invasion.

Conclusion/Significance

Individuals with higher DARC expression, which is associated with higher susceptibility to P. vivax infection, exhibited low frequencies and magnitudes of P. vivax blood-stage specific antibody responses. This may indicate that one of the primary mechanisms by which P. vivax evades host immunity is through DARC indirectly down-regulating humoral responses against erythrocytic invasion and development.

Genetic structure in African populations: implications for human demographic history

The continent of Africa is the source of all anatomically modern humans that dispersed across the planet during the past 100,000 years. As such, African populations are characterized by high genetic diversity and low levels of linkage disequilibrium (LD) among loci, as compared to populations from other continents. African populations also possess a number of genetic adaptations that have evolved in response to the diverse climates, diets, geographic environments, and infectious agents that characterize the African continent. Recently, Tishkoff et al. (2009) performed a genome-wide analysis of substructure based on DNA from 2432 Africans from 121 geographically diverse populations. The authors analyzed patterns of variation at 1327 nuclear microsatellite and insertion/deletion markers and identified 14 ancestral population clusters that correlate well with self-described ethnicity and shared cultural or linguistic properties. The results suggest that African populations may have maintained a large and subdivided population structure throughout much of their evolutionary history. In this chapter, we synthesize recent work documenting evidence of African population structure and discuss the implications for inferences about evolutionary history in both African populations and anatomically modern humans as a whole.

Genome-wide scans for footprints of natural selection

Detecting recent selected ‘genomic footprints’ applies directly to the discovery of disease genes and in the imputation of the formative events that molded modern population genetic structure. The imprints of historic selection/adaptation episodes left in human and animal genomes allow one to interpret modern and ancestral gene origins and modifications. Current approaches to reveal selected regions applied in genome-wide selection scans (GWSSs) fall into eight principal categories: (I) phylogenetic footprinting, (II) detecting increased rates of functional mutations, (III) evaluating divergence versus polymorphism, (IV) detecting extended segments of linkage disequilibrium, (V) evaluating local reduction in genetic variation, (VI) detecting changes in the shape of the frequency distribution (spectrum) of genetic variation, (VII) assessing differentiating between populations (F_ST), and (VIII) detecting excess or decrease in admixture contribution from one population. Here, we review and compare these approaches using available human genome-wide datasets to provide independent verification (or not) of regions found by different methods and using different populations. The lessons learned from GWSSs will be applied to identify genome signatures of historic selective pressures on genes and gene regions in other species with emerging genome sequences. This would offer considerable potential for genome annotation in functional, developmental and evolutionary contexts.

Highly punctuated patterns of population structure on the X chromosome and implications for African evolutionary history

It is well known that average levels of population structure are higher on the X chromosome compared to autosomes in humans. However, there have been surprisingly few analyses on the spatial distribution of population structure along the X chromosome. With publicly available data from the HapMap Project and Perlegen Sciences, we show a strikingly punctuated pattern of X chromosome population structure. Specifically, 87% of X-linked HapMap SNPs within the top 1% of F(ST) values cluster into five distinct loci. The largest of these regions spans 5.4 Mb and contains 66% of the most highly differentiated HapMap SNPs on the X chromosome. We demonstrate that the extreme clustering of highly differentiated SNPs on the X chromosome is not an artifact of ascertainment bias, nor is it specific to the populations genotyped in the HapMap Project. Rather, additional analyses and resequencing data suggest that these five regions have been substrates of recent and strong adaptive evolution. Finally, we discuss the implications that patterns of X-linked population structure have on the evolutionary history of African populations.

Detecting natural selection by empirical comparison to random regions of the genome

Historical episodes of natural selection can skew the frequencies of genetic variants, leaving a signature that can persist for many tens or even hundreds of thousands of years. However, formal tests for selection based on allele frequency skew require strong assumptions about demographic history and mutation, which are rarely well understood. Here, we develop an empirical approach to test for signals of selection that compares patterns of genetic variation at a candidate locus with matched random regions of the genome collected in the same way. We apply this approach to four genes that have been implicated in syndromes of impaired neurological development, comparing the pattern of variation in our re-sequencing data with a large-scale, genomic data set that provides an empirical null distribution. We confirm a previously reported signal at FOXP2, and find a novel signal of selection centered at AHI1, a gene that is involved in motor and behavior abnormalities. The locus is marked by many high frequency derived alleles in non-Africans that are of low frequency in Africans, suggesting that selection at this or a closely neighboring gene occurred in the ancestral population of non-Africans. Our study also provides a prototype for how empirical scans for ancient selection can be carried out once many genomes are sequenced.

From evolutionary genetics to human immunology: how selection shapes host defence genes

Pathogens have always been a major cause of human mortality, so they impose strong selective pressure on the human genome. Data from population genetic studies, including genome-wide scans for selection, are providing important insights into how natural selection has shaped immunity and host defence genes in specific human populations and in the human species as a whole. These findings are helping to delineate genes that are important for host defence and to increase our understanding of how past selection has had an impact on disease susceptibility in modern populations. A tighter integration between population genetic studies and immunological phenotype studies is now necessary to reveal the mechanisms that have been crucial for our past and present survival against infection.

Assessing human variation data for signatures of natural selection

In this article, we highlight some of the different types of natural selection, their effects on patterns of DNA variation, and some of the statistical tests that are commonly used to detect such effects. We also explain some of the relative strengths and weaknesses of different strategies that can be used to detect signatures of natural selection at individual loci. These strategies are illustrated by their application to empirical data from gene variants that are often associated with differences in disease susceptibility. We briefly outline some of the methods proposed to scan the genome for evidence of selection. Finally, we discuss some of the problems associated with identifying signatures of selection and with making inferences about the nature of the selective process.

The Duffy-null state is associated with a survival advantage in leukopenic HIV-infected persons of African ancestry

Persons of African ancestry, on average, have lower white blood cell (WBC) counts than those of European descent (ethnic leukopenia), but whether this impacts negatively on HIV-1 disease course remains unknown. Here, in a large natural history cohort of HIV-infected subjects, we show that, although leukopenia (< 4000 WBC/mm(3) during infection) was associated with an accelerated HIV disease course, this effect was more prominent in leukopenic subjects of European than African ancestry. The African-specific -46C/C genotype of Duffy Antigen Receptor for Chemokines (DARC) confers the malaria-resisting, Duffy-null phenotype, and we found that the recently described association of this genotype with ethnic leukopenia extends to HIV-infected African Americans (AAs). The association of Duffy-null status with HIV disease course differed according to WBC but not CD4(+) T-cell counts, such that leukopenic but not nonleukopenic HIV(+) AAs with DARC -46C/C had a survival advantage compared with all Duffy-positive subjects. This survival advantage became increasingly pronounced in those with progressively lower WBC counts. These data highlight that the interaction between DARC genotype and the cellular milieu defined by WBC counts may influence HIV disease course, and this may provide a partial explanation of why ethnic leukopenia remains benign in HIV-infected AAs, despite immunodeficiency.

Constructing genomic maps of positive selection in humans: where do we go from here?

Identifying targets of positive selection in humans has, until recently, been frustratingly slow, relying on the analysis of individual candidate genes. Genomics, however, has provided the necessary resources to systematically interrogate the entire genome for signatures of natural selection. To date, 21 genome-wide scans for recent or ongoing positive selection have been performed in humans. A key challenge is to begin synthesizing these newly constructed maps of positive selection into a coherent narrative of human evolutionary history and derive a deeper mechanistic understanding of how natural populations evolve. Here, I chronicle the recent history of the burgeoning field of human population genomics, critically assess genome-wide scans for positive selection in humans, identify important gaps in knowledge, and discuss both short- and long-term strategies for traversing the path from the low-resolution, incomplete, and error-prone maps of selection today to the ultimate goal of a detailed molecular, mechanistic, phenotypic, and population genetics characterization of adaptive alleles.

Population differentiation as an indicator of recent positive selection in humans: an empirical evaluation

We have evaluated the extent to which SNPs identified by genomewide surveys as showing unusually high levels of population differentiation in humans have experienced recent positive selection, starting from a set of 32 nonsynonymous SNPs in 27 genes highlighted by the HapMap1 project. These SNPs were genotyped again in the HapMap samples and in the Human Genome Diversity Project-Centre d'Etude du Polymorphisme Humain (HGDP-CEPH) panel of 52 populations representing worldwide diversity; extended haplotype homozygosity was investigated around all of them, and full resequence data were examined for 9 genes (5 from public sources and 4 from new data sets). For 7 of the genes, genotyping errors were responsible for an artifactual signal of high population differentiation and for 2, the population differentiation did not exceed our significance threshold. For the 18 genes with confirmed high population differentiation, 3 showed evidence of positive selection as measured by unusually extended haplotypes within a population, and 7 more did in between-population analyses. The 9 genes with resequence data included 7 with high population differentiation, and 5 showed evidence of positive selection on the haplotype carrying the nonsynonymous SNP from skewed allele frequency spectra; in addition, 2 showed evidence of positive selection on unrelated haplotypes. Thus, in humans, high population differentiation is (apart from technical artifacts) an effective way of enriching for recently selected genes, but is not an infallible pointer to recent positive selection supported by other lines of evidence.

Genome-wide comparisons of variation in linkage disequilibrium

Current genome-wide surveys of common diseases and complex traits fundamentally aim to detect indirect associations where the single nucleotide polymorphisms (SNPs) carrying the association signals are not biologically active but are in linkage disequilibrium (LD) with some unknown functional polymorphisms. Reproducing any novel discoveries from these genome-wide scans in independent studies is now a prerequisite for the putative findings to be accepted. Significant differences in patterns of LD between populations can affect the portability of phenotypic associations when the replication effort or meta-analyses are attempted in populations that are distinct from the original population in which the genome-wide study is performed. Here, we introduce a novel method for genome-wide analyses of LD variations between populations that allow the identification of candidate regions with different patterns of LD. The evidence of LD variation provided by the introduced method correlated with the degree of differences in the frequencies of the most common haplotype across the populations. Identified regions also resulted in greater variation in the success of replication attempts compared with random regions in the genome. A separate permutation strategy introduced for assessing LD variation in the absence of genome-wide data also correctly identified the expected variation in LD patterns in two well-established regions undergoing strong population-specific evolutionary pressure. Importantly, this method addresses whether a failure to reproduce a disease association in a disparate population is due to underlying differences in LD structure with an unknown functional polymorphism, which is vital in the current climate of replicating and fine-mapping established findings from genome-wide association studies.