The T allele of a single-nucleotide polymorphism 13.9 kb upstream of the lactase gene (LCT) (C-13.9kbT) does not predict or cause the lactase-persistence phenotype in Africans

The ability to digest the milk sugar lactose as an adult (lactase persistence) is a variable genetic trait in human populations. The lactase-persistence phenotype is found at low frequencies in the majority of populations in sub-Saharan Africa that have been tested, but, in some populations, particularly pastoral groups, it is significantly more frequent. Recently, a CT polymorphism located 13.9 kb upstream of exon 1 of the lactase gene (LCT) was shown in a Finnish population to be closely associated with the lactase-persistence phenotype (Enattah et al. 2002). We typed this polymorphism in 1,671 individuals from 20 distinct cultural groups in seven African countries. It was possible to match seven of the groups tested with groups from the literature for whom phenotypic information is available. In five of these groups, the published frequencies of lactase persistence are >/=25%. We found the T allele to be so rare that it cannot explain the frequency of the lactase-persistence phenotype throughout Africa. By use of a statistical procedure to take phenotyping and sampling errors into account, the T-allele frequency was shown to be significantly different from that predicted in five of the African groups. Only the Fulbe and Hausa from Cameroon possessed the T allele at a level consistent with phenotypic observations (as well as an Irish sample used for comparison). We conclude that the C-13.9kbT polymorphism is not a predictor of lactase persistence in sub-Saharan Africans. We also present Y-chromosome data that are consistent with previously reported evidence for a back-migration event into Cameroon, and we comment on the implications for the introgression of the -13.9kb*T allele.     

Lactase haplotype diversity in the Old World

Lactase persistence, the genetic trait in which intestinal lactase activity persists at childhood levels into adulthood, varies in frequency in different human populations, being most frequent in northern Europeans and certain African and Arabian nomadic tribes, who have a history of drinking fresh milk. Selection is likely to have played an important role in establishing these different frequencies since the development of agricultural pastoralism approximately 9,000 years ago. We have previously shown that the element responsible for the lactase persistence/nonpersistence polymorphism in humans is cis-acting to the lactase gene and that lactase persistence is associated, in Europeans, with the most common 70-kb lactase haplotype, A. We report here a study of the 11-site haplotype in 1,338 chromosomes from 11 populations that differ in lactase persistence frequency. Our data show that haplotype diversity was generated both by point mutations and recombinations. The four globally common haplotypes (A, B, C, and U) are not closely related and have different distributions; the A haplotype is at high frequencies only in northern Europeans, where lactase persistence is common; and the U haplotype is virtually absent from Indo-European populations. Much more diversity is seen in sub-Saharan Africans than in non-Africans, consistent with an "Out of Africa" model for peopling of the Old World. Analysis of recent recombinant haplotypes by allele-specific PCR, along with deduction of the root haplotype from chimpanzee sequence, allowed construction of a haplotype network that assisted in evaluation of the relative roles of drift and selection in establishing the haplotype frequencies in the different populations. We suggest that genetic drift was important in shaping the general pattern of non-African haplotype diversity, with recent directional selection in northern Europeans for the haplotype associated with lactase persistence.     

Lactose digestion and the evolutionary genetics of lactase persistence

It has been known for some 40 years that lactase production persists into adult life in some people but not in others. However, the mechanism and evolutionary significance of this variation have proved more elusive, and continue to excite the interest of investigators from different disciplines. This genetically determined trait differs in frequency worldwide and is due to cis-acting polymorphism of regulation of lactase gene expression. A single nucleotide polymorphism located 13.9 kb upstream from the lactase gene (C-13910 > T) was proposed to be the cause, and the −13910*T allele, which is widespread in Europe was found to be located on a very extended haplotype of 500 kb or more. The long region of haplotype conservation reflects a recent origin, and this, together with high frequencies, is evidence of positive selection, but also means that −13910*T might be an associated marker, rather than being causal of lactase persistence itself. Doubt about function was increased when it was shown that the original SNP did not account for lactase persistence in most African populations. However, the recent discovery that there are several other SNPs associated with lactase persistence in close proximity (within 100 bp), and that they all reside in a piece of sequence that has enhancer function in vitro, does suggest that they may each be functional, and their occurrence on different haplotype backgrounds shows that several independent mutations led to lactase persistence. Here we provide access to a database of worldwide distributions of lactase persistence and of the C-13910*T allele, as well as reviewing lactase molecular and population genetics and the role of selection in determining present day distributions of the lactase persistence phenotype. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Herders of Indian and European cattle share their predominant allele for lactase persistence

Milk consumption and lactose digestion after weaning are exclusively human traits made possible by the continued production of the enzyme lactase in adulthood. Multiple independent mutations in a 100-bp region--part of an enhancer--approximately 14-kb upstream of the LCT gene are associated with this trait in Europeans and pastoralists from Saudi Arabia and Africa. However, a single mutation of purported western Eurasian origin accounts for much of observed lactase persistence outside Africa. Given the high levels of present-day milk consumption in India, together with archaeological and genetic evidence for the independent domestication of cattle in the Indus valley roughly 7,000 years ago, we sought to determine whether lactase persistence has evolved independently in the subcontinent. Here, we present the results of the first comprehensive survey of the LCT enhancer region in south Asia. Having genotyped 2,284 DNA samples from across the Indian subcontinent, we find that the previously described west Eurasian -13910 C>T mutation accounts for nearly all the genetic variation we observed in the 400- to 700-bp LCT regulatory region that we sequenced. Geography is a significant predictor of -13910*T allele frequency, and consistent with other genomic loci, its distribution in India follows a general northwest to southeast declining pattern, although frequencies among certain neighboring populations vary substantially. We confirm that the mutation is identical by descent to the European allele and is associated with the same>1 Mb extended haplotype in both populations.     

The Origins of Lactase Persistence in Europe

Lactase persistence (LP) is common among people of European ancestry, but with the exception of some African, Middle Eastern and southern Asian groups, is rare or absent elsewhere in the world. Lactase gene haplotype conservation around a polymorphism strongly associated with LP in Europeans (−13,910 C/T) indicates that the derived allele is recent in origin and has been subject to strong positive selection. Furthermore, ancient DNA work has shown that the −13,910*T (derived) allele was very rare or absent in early Neolithic central Europeans. It is unlikely that LP would provide a selective advantage without a supply of fresh milk, and this has lead to a gene-culture coevolutionary model where lactase persistence is only favoured in cultures practicing dairying, and dairying is more favoured in lactase persistent populations. We have developed a flexible demic computer simulation model to explore the spread of lactase persistence, dairying, other subsistence practices and unlinked genetic markers in Europe and western Asia''s geographic space. Using data on −13,910*T allele frequency and farming arrival dates across Europe, and approximate Bayesian computation to estimate parameters of interest, we infer that the −13,910*T allele first underwent selection among dairying farmers around 7,500 years ago in a region between the central Balkans and central Europe, possibly in association with the dissemination of the Neolithic Linearbandkeramik culture over Central Europe. Furthermore, our results suggest that natural selection favouring a lactase persistence allele was not higher in northern latitudes through an increased requirement for dietary vitamin D. Our results provide a coherent and spatially explicit picture of the coevolution of lactase persistence and dairying in Europe.     

Lactase persistence in central Asia: phenotype, genotype, and evolution

The aim of the present study is to document the evolution of the lactase persistence trait in Central Asia, a geographical area that is thought to have been a region of long-term pastoralism. Several ethnic groups co-exist in this area: Indo-Iranian speakers who are traditionally agriculturist (Tajik) and Turkic speakers who used to be nomadic herders (Kazakh, Karakalpak, Kyrgyz, Turkmen). It was recently demonstrated that horse milking practice existed in the Botai culture of Kazakhstan as early as 5,500 BP ( Outram et al. 2009 ). However, the frequency of the lactase persistence trait and its genetic basis in Central Asian populations remain largely unknown. We propose here the first genotype-phenotype study of lactase persistence in Central Asia based on 183 individuals, as well as the estimation of the time of expansion of the lactase-persistence associated polymorphism. Our results show a remarkable genetic-phenotypic correlation, with the causal polymorphism being the same than in Europe (-13.910C>T, rs4988235). The lactase persistence trait is at low frequency in these populations: between 25% and 32% in the Kazakh population (traditionally herders), according to phenotype used, and between 11% and 30% in the Tajiko-Uzbek population (agriculturalists). The difference in lactase persistence between populations, even if small, is significant when using individuals concordant for both excretion of breath hydrogen and the lactose tolerance blood glucose test phenotypes (P = 0.018, 25% for Kazakh vs. 11% for Tajiko-Uzbeks), and the difference in frequency of the -13.910*T allele is almost significant (P = 0.06, 30% for Kazakhs vs. 19% for Tajiko-Uzbeks). Using the surrounding haplotype, we estimate a date of expansion of the T allele around 6,000-12,000 yrs ago, which is consistent with archaeological records for the emergence of agropastoralism and pastoralism in Central Asia.     

Impact of Selection and Demography on the Diffusion of Lactase Persistence

Background

The lactase enzyme allows lactose digestion in fresh milk. Its activity strongly decreases after the weaning phase in most humans, but persists at a high frequency in Europe and some nomadic populations. Two hypotheses are usually proposed to explain the particular distribution of the lactase persistence phenotype. The gene-culture coevolution hypothesis supposes a nutritional advantage of lactose digestion in pastoral populations. The calcium assimilation hypothesis suggests that carriers of the lactase persistence allele(s) (LCT*P) are favoured in high-latitude regions, where sunshine is insufficient to allow accurate vitamin-D synthesis. In this work, we test the validity of these two hypotheses on a large worldwide dataset of lactase persistence frequencies by using several complementary approaches.

Methodology

We first analyse the distribution of lactase persistence in various continents in relation to geographic variation, pastoralism levels, and the genetic patterns observed for other independent polymorphisms. Then we use computer simulations and a large database of archaeological dates for the introduction of domestication to explore the evolution of these frequencies in Europe according to different demographic scenarios and selection intensities.

Conclusions

Our results show that gene-culture coevolution is a likely hypothesis in Africa as high LCT*P frequencies are preferentially found in pastoral populations. In Europe, we show that population history played an important role in the diffusion of lactase persistence over the continent. Moreover, selection pressure on lactase persistence has been very high in the North-western part of the continent, by contrast to the South-eastern part where genetic drift alone can explain the observed frequencies. This selection pressure increasing with latitude is highly compatible with the calcium assimilation hypothesis while the gene-culture coevolution hypothesis cannot be ruled out if a positively selected lactase gene was carried at the front of the expansion wave during the Neolithic transition in Europe.     

Ancient DNA Analysis Reveals High Frequency of European Lactase Persistence Allele (T-13910) in Medieval Central Europe

Ruminant milk and dairy products are important food resources in many European, African, and Middle Eastern societies. These regions are also associated with derived genetic variants for lactase persistence. In mammals, lactase, the enzyme that hydrolyzes the milk sugar lactose, is normally down-regulated after weaning, but at least five human populations around the world have independently evolved mutations regulating the expression of the lactase-phlorizin-hydrolase gene. These mutations result in a dominant lactase persistence phenotype and continued lactase tolerance in adulthood. A single nucleotide polymorphism (SNP) at C/T-13910 is responsible for most lactase persistence in European populations, but when and where the T-13910 polymorphism originated and the evolutionary processes by which it rose to high frequency in Europe have been the subject of strong debate. A history of dairying is presumed to be a prerequisite, but archaeological evidence is lacking. In this study, DNA was extracted from the dentine of 36 individuals excavated at a medieval cemetery in Dalheim, Germany. Eighteen individuals were successfully genotyped for the C/T-13910 SNP by molecular cloning and sequencing, of which 13 (72%) exhibited a European lactase persistence genotype: 44% CT, 28% TT. Previous ancient DNA-based studies found that lactase persistence genotypes fall below detection levels in most regions of Neolithic Europe. Our research shows that by AD 1200, lactase persistence frequency had risen to over 70% in this community in western Central Europe. Given that lactase persistence genotype frequency in present-day Germany and Austria is estimated at 71–80%, our results suggest that genetic lactase persistence likely reached modern levels before the historic population declines associated with the Black Death, thus excluding plague-associated evolutionary forces in the rise of lactase persistence in this region. This new evidence sheds light on the dynamic evolutionary history of the European lactase persistence trait and its global cultural implications.     

Microsatellite variation and evolution of human lactase persistence

The levels of haplotype diversity within the lineages defined by two single-nucleotide polymorphisms (SNPs) (–13910 C/T and –22018 G/A) associated with human lactase persistence were assessed with four fast-evolving microsatellite loci in 794 chromosomes from Portugal, Italy, Fulbe from Cameroon, São Tomé and Mozambique. Age estimates based on the intraallelic microsatellite variation indicate that the –13910*T allele, which is more tightly associated with lactase persistence, originated in Eurasia before the Neolithic and after the emergence of modern humans outside Africa. We detected significant departures from neutrality for the –13910*T variant in geographically and evolutionary distant populations from southern Europe (Portuguese and Italians) and Africa (Fulbe) by using a neutrality test based on the congruence between the frequency of the allele and the levels of intraallelic variability measured by the number of mutations in adjacent microsatellites. This result supports the role of selection in the evolution of lactase persistence, ruling out possible confounding effects from recombination suppression and population history. Reevaluation of the available evidence on variation of the –13910 and –22018 loci indicates that lactase persistence probably originated from different mutations in Europe and most of Africa, even if 13910*T is not the causal allele, suggesting that selective pressure could have promoted the convergent evolution of the trait. Our study shows that a limited number of microsatellite loci may provide sufficient resolution to reconstruct key aspects of the evolutionary history of lactase persistence, providing an alternative to approaches based on large numbers of SNPs.Electronic supplementary material Supplementary material is available for this article at     

Modeling and E-M estimation of haplotype-specific relative risks from genotype data for a case-control study of unrelated individuals

The US National Cancer Institute has recently sponsored the formation of a Cohort Consortium (http://2002.cancer.gov/scpgenes.htm) to facilitate the pooling of data on very large numbers of people, concerning the effects of genes and environment on cancer incidence. One likely goal of these efforts will be generate a large population-based case-control series for which a number of candidate genes will be investigated using SNP haplotype as well as genotype analysis. The goal of this paper is to outline the issues involved in choosing a method of estimating haplotype-specific risk estimates for such data that is technically appropriate and yet attractive to epidemiologists who are already comfortable with odds ratios and logistic regression. Our interest is to develop and evaluate extensions of methods, based on haplotype imputation, that have been recently described (Schaid et al., Am J Hum Genet, 2002, and Zaykin et al., Hum Hered, 2002) as providing score tests of the null hypothesis of no effect of SNP haplotypes upon risk, which may be used for more complex tasks, such as providing confidence intervals, and tests of equivalence of haplotype-specific risks in two or more separate populations. In order to do so we (1) develop a cohort approach towards odds ratio analysis by expanding the E-M algorithm to provide maximum likelihood estimates of haplotype-specific odds ratios as well as genotype frequencies; (2) show how to correct the cohort approach, to give essentially unbiased estimates for population-based or nested case-control studies by incorporating the probability of selection as a case or control into the likelihood, based on a simplified model of case and control selection, and (3) finally, in an example data set (CYP17 and breast cancer, from the Multiethnic Cohort Study) we compare likelihood-based confidence interval estimates from the two methods with each other, and with the use of the single-imputation approach of Zaykin et al. applied under both null and alternative hypotheses. We conclude that so long as haplotypes are well predicted by SNP genotypes (we use the Rh2 criteria of Stram et al. [1]) the differences between the three methods are very small and in particular that the single imputation method may be expected to work extremely well.     

Cystic fibrosis gene variability in two southern Brazilian Amerindian populations: analysis of the deltaF508 mutation and the KM19 and XV2C haplotypes

The frequencies of the deltaF508 deletion, the most common cystic fibrosis mutation in Europeans and European-derived populations, and the XV2C and KM19 restriction fragment length polymorphisms that are tightly linked to the CFTR locus vary among populations. To determine the distribution of these extragenic markers and of the deltaF508 mutation, we analyzed 326 chromosomes of individuals from two South American Indian populations, the Guarani and the Kaingang. The allele and haplotype frequencies differed greatly between the two populations as well as among Amerindians and normal European Brazilians and European Brazilian cystic fibrosis patients. The absence of the deltaF508 mutation and the B haplotype are in agreement with the hypothesis that the deltaF508 mutation occurred after the divergence of these two populations. This finding is useful for populations containing a large Amerindian component and helps us to understand the origins of the deltaF508 deletion, the most common cystic fibrosis mutation in Europeans and European-derived populations, as well as the different incidences of cystic fibrosis in continental groups.     

Lactase Persistence and Lipid Pathway Selection in the Maasai

The Maasai are a pastoral people in Kenya and Tanzania, whose traditional diet of milk, blood and meat is rich in lactose, fat and cholesterol. In spite of this, they have low levels of blood cholesterol, and seldom suffer from gallstones or cardiac diseases. Field studies in the 1970s suggested that the Maasai have a genetic adaptation for cholesterol homeostasis. Analysis of HapMap 3 data using Fixation Index (Fst) and two metrics of haplotype diversity: the integrated Haplotype Score (iHS) and the Cross Population Extended Haplotype Homozygosity (XP-EHH), identified genomic regions and single nucleotide polymorphisms (SNPs) as strong candidates for recent selection for lactase persistence and cholesterol regulation in 143–156 founder individuals from the Maasai population in Kinyawa, Kenya (MKK). The non-synonmous SNP with the highest genome-wide Fst was the TC polymorphism at rs2241883 in Fatty Acid Binding Protein 1(FABP1), known to reduce low density lipoprotein and tri-glyceride levels in Europeans. The strongest signal identified by all three metrics was a 1.7 Mb region on Chr2q21. This region contains the genes LCT (Lactase) and MCM6 (Minichromosome Maintenance Complex Component) involved in lactase persistence, and the gene Rab3GAP1 (Rab3 GTPase-activating Protein Catalytic Subunit), which contains polymorphisms associated with total cholesterol levels in a genome-wide association study of >100,000 individuals of European ancestry. Sanger sequencing of DNA from six MKK samples showed that the GC-14010 polymorphism in the MCM6 gene, known to be associated with lactase persistence in Africans, is segregating in MKK at high frequency (∼58%). The Cytochrome P450 Family 3 Subfamily A (CYP3A) cluster of genes, involved in cholesterol metabolism, was identified by Fst and iHS as candidate loci under selection. Overall, our study identified several specific genomic regions under selection in the Maasai which contain polymorphisms in genes associated with lactase persistence and cholesterol regulation.     

Evidence of still-ongoing convergence evolution of the lactase persistence T-13910 alleles in humans

A single-nucleotide variant, C/T(-13910), located 14 kb upstream of the lactase gene (LCT), has been shown to be completely correlated with lactase persistence (LP) in northern Europeans. Here, we analyzed the background of the alleles carrying the critical variant in 1,611 DNA samples from 37 populations. Our data show that the T(-13910) variant is found on two different, highly divergent haplotype backgrounds in the global populations. The first is the most common LP haplotype (LP H98) present in all populations analyzed, whereas the others (LP H8-H12), which originate from the same ancestral allelic haplotype, are found in geographically restricted populations living west of the Urals and north of the Caucasus. The global distribution pattern of LP T(-13910) H98 supports the Caucasian origin of this allele. Age estimates based on different mathematical models show that the common LP T(-13910) H98 allele (approximately 5,000-12,000 years old) is relatively older than the other geographically restricted LP alleles (approximately 1,400-3,000 years old). Our data about global allelic haplotypes of the lactose-tolerance variant imply that the T(-13910) allele has been independently introduced more than once and that there is a still-ongoing process of convergent evolution of the LP alleles in humans.     

Linkage disequilibrium in the North American Holstein population

Linkage disequilibrium was estimated using 7119 single nucleotide polymorphism markers across the genome and 200 animals from the North American Holstein cattle population. The analysis of maternally inherited haplotypes revealed strong linkage disequilibrium ( r ²   >   0.8) in genomic regions of ∼50 kb or less. While linkage disequilibrium decays as a function of genomic distance, genomic regions within genes showed greater linkage disequilibrium and greater variation in linkage disequilibrium compared with intergenic regions. Identification of haplotype blocks could characterize the most common haplotypes. Although maximum haplotype block size was over 1 Mb, mean block size was 26–113 kb by various definitions, which was larger than that observed in humans (∼10 kb). Effective population size of the dairy cattle population was estimated from linkage disequilibrium between single nucleotide polymorphism marker pairs in various haplotype ranges. Rapid reduction of effective population size of dairy cattle was inferred from linkage disequilibrium in recent generations. This result implies a loss of genetic diversity because of the high rate of inbreeding and high selection intensity in dairy cattle. The pattern observed in this study indicated linkage disequilibrium in the current dairy cattle population could be exploited to refine mapping resolution. Changes in effective population size during past generations imply a necessity of plans to maintain polymorphism in the Holstein population.     

Detecting and Characterizing Genomic Signatures of Positive Selection in Global Populations

Natural selection is a significant force that shapes the architecture of the human genome and introduces diversity across global populations. The question of whether advantageous mutations have arisen in the human genome as a result of single or multiple mutation events remains unanswered except for the fact that there exist a handful of genes such as those that confer lactase persistence, affect skin pigmentation, or cause sickle cell anemia. We have developed a long-range-haplotype method for identifying genomic signatures of positive selection to complement existing methods, such as the integrated haplotype score (iHS) or cross-population extended haplotype homozygosity (XP-EHH), for locating signals across the entire allele frequency spectrum. Our method also locates the founder haplotypes that carry the advantageous variants and infers their corresponding population frequencies. This presents an opportunity to systematically interrogate the whole human genome whether a selection signal shared across different populations is the consequence of a single mutation process followed subsequently by gene flow between populations or of convergent evolution due to the occurrence of multiple independent mutation events either at the same variant or within the same gene. The application of our method to data from 14 populations across the world revealed that positive-selection events tend to cluster in populations of the same ancestry. Comparing the founder haplotypes for events that are present across different populations revealed that convergent evolution is a rare occurrence and that the majority of shared signals stem from the same evolutionary event.     

Genetic evidence for ongoing balanced selection at human DNA repair genes ERCC8, FANCC, and RAD51C

Using the 2.6 million single nucleotide polymorphism (SNP) genotype datasets from Perlegen Sciences and the Haplotype Map (HapMap) project (Phase I freeze), a probabilistic search for the landscape exhibited by positive Darwinian selection was conducted (Wang et al., 2006). By sorting each high frequency allele by homozygosity, we search for the expected decay of adjacent SNP linkage disequilibrium (LD) at recently selected alleles, eliminating the need for inferring haplotype. We designate this approach the LD decay (LDD) test. Cluster analysis indicates that approximately 3000 sites of recent inferred selection are present in human DNA, representing approximately 1800 genes. Prior simulation studies (Wang et al., 2006) indicate that this novel LDD test, at the Mb scale employed, effectively distinguishes selection from other causes of extensive LD, such as inversions, population bottlenecks and admixture. Based on over-representation analysis, these prior studies have shown that several predominant biological themes are common in inferred selected alleles, including genes involved with DNA metabolism and repair. Here, we show that three of these DNA repair genes, ERCC8, Fanconi Anemia Complementation Group C (FANCC), and RAD51C, exhibit genomic architectures consistent with ongoing balanced selection over the last 40,000-50,000 years.