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.     

High frequency of lactose intolerance in a prehistoric hunter-gatherer population in northern Europe

Background  

Genes and culture are believed to interact, but it has been difficult to find direct evidence for the process. One candidate example that has been put forward is lactase persistence in adulthood, i.e. the ability to continue digesting the milk sugar lactose after childhood, facilitating the consumption of raw milk. This genetic trait is believed to have evolved within a short time period and to be related with the emergence of sedentary agriculture.     

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.     

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.     

A strong association of axillary osmidrosis with the wet earwax type determined by genotyping of the ABCC11 gene

Background

Lactase non-persistence is a condition where lactase activity is decreased in the intestinal wall after weaning. In European derived populations a single nucleotide polymorphism (SNP) C/T_-13910 residing 13.9 kb upstream from the lactase gene has been shown to define lactase activity, and several other single nucleotide polymorphisms (G/C_-14010 T/G_-13915, C/G_-13907 and T/C_-13913) in the same region have been identified in African and Middle East populations.

Results

The T_-13910 allele most common in European populations was present in 21.8% mixed ancestry (N = 62) individuals and it was absent in the Xhosa (N = 109) and Ghana (N = 196) subjects. Five other substitutions were also found in the region covering the previously reported variants in African and Middle East populations. These included the G/C_-14010 variant common in Kenyan and Tanzanian populations, which was present in 12.8% of Xhosa population and in 8.1% of mixed ancestry subjects. Two novel substitutions (C/T_-14091 and A/C_-14176) and one previously reported substitution G/A_-13937 (rs4988234) were less common and present only in the Xhosa population. One novel substitution G/A_-14107 was present in the Xhosa and Ghanaian populations. None of the other previously reported variants were identified.

Conclusion

Identification of the G/C_-14010 variant in the Xhosa population, further confirms their genetic relatedness to other nomadic populations members that belong to the Bantu linguistic group in Tanzania and Kenya. Further studies are needed to confirm the possible relationship of the novel substitutions to the lactase persistence trait.     

Evolution of lactase persistence: an example of human niche construction

Niche construction is the process by which organisms construct important components of their local environment in ways that introduce novel selection pressures. Lactase persistence is one of the clearest examples of niche construction in humans. Lactase is the enzyme responsible for the digestion of the milk sugar lactose and its production decreases after the weaning phase in most mammals, including most humans. Some humans, however, continue to produce lactase throughout adulthood, a trait known as lactase persistence. In European populations, a single mutation (-13910*T) explains the distribution of the phenotype, whereas several mutations are associated with it in Africa and the Middle East. Current estimates for the age of lactase persistence-associated alleles bracket those for the origins of animal domestication and the culturally transmitted practice of dairying. We report new data on the distribution of -13910*T and summarize genetic studies on the diversity of lactase persistence worldwide. We review relevant archaeological data and describe three simulation studies that have shed light on the evolution of this trait in Europe. These studies illustrate how genetic and archaeological information can be integrated to bring new insights to the origins and spread of lactase persistence. Finally, we discuss possible improvements to these models.     

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.     

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.     

Lactase persistence in Tunisia as a result of admixture with other Mediterranean populations

Background

The ability to digest lactose after weaning, namely, lactase persistence (LP), is encoded by polymorphisms in the MCM6 gene and varies widely in frequency among different human populations. Although, evolution of LP-related genetic variants was investigated in many groups of Sub-Saharan African, Middle Eastern, and European ancestry, only few studies have focused on populations from North Africa and no data are especially available from the Tunisian one. For this reason, there is an urgent need to investigate the frequency patterns at these loci in Tunisia since this adaptive trait is implicated in health.

Methods

Forty SNPs covering the LCT/MCM6 genes and including the two functional variants ??13,910 C > T and ??22,018 G > A were genotyped in 117 Tunisian individuals using the Sequenom Mass Array technology. The observed nucleotide and haplotype patterns of variation were then compared with those of several African, European, and Mediterranean human groups for which comparable data were publicly available. Admixture analysis on a 5 Mb genomic region surrounding the LCT/MCM6 loci was also performed by extracting genotypes from a previously generated genome-wide dataset in order to deepen the reconstruction of the evolutionary history of these loci.

Results

We found that lactase non-persistence (LNP)-related alleles and haplotypes were predominantly present in the examined population. A clear differentiation between Tunisian, African, and North European/North Italian samples was found, while the Tunisian population showed more genetic affinity to Central and South Italian groups.

Conclusions

Our study provided a first report of LP-associated alleles and haplotypes in the Tunisian population. We highlighted a gradient followed by LP diffusion from Europe to North Africa. Based on the rich historic background of Tunisia, we suggest that this adaptive trait was introduced in that geographic region by a relatively recent gene flow.

     

On the edge of Bantu expansions: mtDNA, Y chromosome and lactase persistence genetic variation in southwestern Angola

Background  

Current information about the expansion of Bantu-speaking peoples is hampered by the scarcity of genetic data from well identified populations from southern Africa. Here, we fill an important gap in the analysis of the western edge of the Bantu migrations by studying for the first time the patterns of Y-chromosome, mtDNA and lactase persistence genetic variation in four representative groups living around the Namib Desert in southwestern Angola (Ovimbundu, Ganguela, Nyaneka-Nkumbi and Kuvale). We assessed the differentiation between these populations and their levels of admixture with Khoe-San groups, and examined their relationship with other sub-Saharan populations. We further combined our dataset with previously published data on Y-chromosome and mtDNA variation to explore a general isolation with migration model and infer the demographic parameters underlying current genetic diversity in Bantu populations.     

Genome annotation of a 1.5 Mb region of human chromosome 6q23 encompassing a quantitative trait locus for fetal hemoglobin expression in adults

Background  

Heterocellular hereditary persistence of fetal hemoglobin (HPFH) is a common multifactorial trait characterized by a modest increase of fetal hemoglobin levels in adults. We previously localized a Quantitative Trait Locus for HPFH in an extensive Asian-Indian kindred to chromosome 6q23. As part of the strategy of positional cloning and a means towards identification of the specific genetic alteration in this family, a thorough annotation of the candidate interval based on a strategy ofin silico/ wet biology approach with comparative genomics was conducted.     

Linked genetic variants on chromosome 10 control ear morphology and body mass among dog breeds

Background

The domestic dog is a rich resource for mapping the genetic components of phenotypic variation due to its unique population history involving strong artificial selection. Genome-wide association studies have revealed a number of chromosomal regions where genetic variation associates with morphological characters that typify dog breeds. A region on chromosome 10 is among those with the highest levels of genetic differentiation between dog breeds and is associated with body mass and ear morphology, a common motif of animal domestication. We characterised variation in this region to uncover haplotype structure and identify candidate functional variants.

Results

We first identified SNPs that strongly associate with body mass and ear type by comparing sequence variation in a 3 Mb region between 19 breeds with a variety of phenotypes. We next genotyped a subset of 123 candidate SNPs in 288 samples from 46 breeds to identify the variants most highly associated with phenotype and infer haplotype structure. A cluster of SNPs that associate strongly with the drop ear phenotype is located within a narrow interval downstream of the gene MSRB3, which is involved in human hearing. These SNPs are in strong genetic linkage with another set of variants that correlate with body mass within the gene HMGA2, which affects human height. In addition we find evidence that this region has been under selection during dog domestication, and identify a cluster of SNPs within MSRB3 that are highly differentiated between dogs and wolves.

Conclusions

We characterise genetically linked variants that potentially influence ear type and body mass in dog breeds, both key traits that have been modified by selective breeding that may also be important for domestication. The finding that variants on long haplotypes have effects on more than one trait suggests that genetic linkage can be an important determinant of the phenotypic response to selection in domestic animals.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1702-2) contains supplementary material, which is available to authorized users.     

Frequency of lactase persistence genotype in a healthy Polish population

The majority of mammals are unable to digest lactose due to post-weaning deactivation of the LCT gene, which is responsible for encoding the enzyme lactase (i.e., adult-type hypolactasia). A substitution of C with T at position −13910 bp upstream of the LCT gene has been linked to the lactase persistence phenotype in European populations. We investigated the frequency of LCT-13910C>T polymorphism in 223 blood donors from central Poland. All samples were genotyped by polymerase chain reaction and direct sequencing. The LCT-13910 T allele (lactase persistence) was present in 51% of individuals sampled from the Polish population. We did not find any non-European variants associated with lactase persistence (LCT-13907C>G, LCT-13913T>C, LCT-13915T>G), or any new polymorphisms within the sequenced region. Allele frequencies obtained are in agreement with results from other countries and confirm the unique pattern of distribution of the LCT-13910C>T genotype in Europe.     

Modeling expression quantitative trait loci in data combining ethnic populations

Background  

Combining data from different ethnic populations in a study can increase efficacy of methods designed to identify expression quantitative trait loci (eQTL) compared to analyzing each population independently. In such studies, however, the genetic diversity of minor allele frequencies among populations has rarely been taken into account. Due to the fact that allele frequency diversity and population-level expression differences are present in populations, a consensus regarding the optimal statistical approach for analysis of eQTL in data combining different populations remains inconclusive.     

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.     

PCR-RFLP genotyping assay for a lactase persistence polymorphism upstream of the lactase-phlorizin hydrolase gene

The majority of the world's human population experiences a decline of lactase gene expression during maturation, so-called lactase nonpersistence. Thus, adults with lactase nonpersistence are susceptible to developing symptoms of lactose intolerance. By contrast, lactase persistence is an autosomal dominant heritable condition that results in a high level of lactase gene expression throughout adulthood and sustained lactose tolerance. Lactase persistence has recently been correlated with a single nucleotide genetic variant (a C --> T mutation) located 13,910 bases upstream from the lactase structural gene. We aimed to develop a restriction fragment length polymorphism (RFLP) method of detecting the C/T variants as a means of identifying individuals genetically inclined toward lactase persistence or nonpersistence. Genomic DNA in a 210-bp region surrounding the -13,910-bp variant site was PCR amplified with unique primers designed to avoid or mutate adjacent restriction sites. The amplified DNA was digested with a restriction enzyme, CviJI, that recognizes the base pair sequence generated by the lactase nonpersistence variant. Restriction digest gel analysis yielded DNA fragments of the expected diagnostic molecular weight sizes for individuals that were homozygote or heterozygote for the lactase persistence and nonpersistence variants. The genotypes predicted by the RFLP-based method were confirmed by DNA sequence analysis. The RFLP-based method provides a quick and noninvasive means of molecular detection of the presence or absence of the lactase persistence variant.     

The −14010*C variant associated with lactase persistence is located between an Oct-1 and HNF1α binding site and increases lactase promoter activity

In most people worldwide intestinal lactase expression declines in childhood. In many others, particularly in Europeans, lactase expression persists into adult life. The lactase persistence phenotype is in Europe associated with the −13910*T single nucleotide variant located 13,910 bp upstream the lactase gene in an enhancer region that affects lactase promoter activity. This variant falls in an Oct-1 binding site and shows greater Oct-1 binding than the ancestral variant and increases enhancer activity. Several other variants have been identified very close to the −13910 position, which are associated with lactase persistence in the Middle East and Africa. One of them, the −14010*C, is associated with lactase persistence in Africa. Here we show by deletion analysis that the −14010 position is located in a 144 bp region that reduces the enhancer activity. In transfections the −14010*C allele shows a stronger enhancer effect than the ancestral −14010*G allele. Binding sites for Oct-1 and HNF1α surrounding the −14010 position were identified by gel shift assays, which indicated that −14010*C has greater binding affinity to Oct-1 than −14010*G.     

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.