Several Different Lactase Persistence Associated Alleles and High Diversity of the Lactase Gene in the Admixed Brazilian Population

Adult-type hypolactasia is a common phenotype caused by the lactase enzyme deficiency. The −13910 C>T polymorphism, located 14 Kb upstream of the lactase gene (LCT) in the MCM6 gene was associated with lactase persistence (LP) in Europeans. This polymorphism is rare in Africa but several other variants associated with lactase persistence were observed in Africans. The aims of this study were to identify polymorphisms in the MCM6 region associated with the lactase persistence phenotype and to determine the distribution of LCT gene haplotypes in 981 individuals from North, Northeast and South Brazil. These polymorphisms were genotyped by PCR based methods and sequencing. The −13779*C,−13910*T, −13937*A, −14010*C, −14011*T LP alleles previously described in the MCM6 gene region that acts as an enhancer for the LCT gene were identified in Brazilians. The most common LP allele was −13910*T. Its frequency was highly correlated with European ancestry in the Brazilian populations investigated. The −13910*T was higher (0.295) in southern Brazilians of European ancestry and lower (0.175) in the Northern admixed population. LCT haplotypes were derived from the 10 LCT SNPs genotyped. Overall twenty six haplotypes previously described were identified in the four Brazilian populations studied. The Multidimensional Scaling analysis showed that Belém, in the north, was closer to Amerindians. Northeastern and southern Afro-descendants were more related with Bantu-speaking South Africans whereas the Southern population with European ancestry grouped with Southern and Northern Europeans. This study shows a high variability considering the number of LCT haplotypes observed. Due to the highly admixed nature of the Brazilian populations, the diagnosis of hypolactasia in Brazil, based only in the investigation of the −13910*T allele is an oversimplification.     

Tracing Pastoralist Migrations to Southern Africa with Lactase Persistence Alleles

1. Download : Download high-res image (202KB)
2. Download : Download full-size image

     

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 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.     

Whole-Genome Sequencing Reveals Lactase Persistence Adaptation in European Dogs

Coexistence and cooperation between dogs and humans over thousands of years have supported convergent evolutionary processes in the two species. Previous studies found that Eurasian dogs evolved into a distinct geographic cluster. In this study, we used the genomes of 242 European dogs, 38 Southeast Asian indigenous (SEAI) dogs, and 41 gray wolves to identify adaptation of European dogs . We report 86 unique positively selected genes in European dogs, among which is LCT (lactase). LCT encodes lactase, which is fundamental for the digestion of lactose. We found that an A-to-G mutation (chr19:38,609,592) is almost fixed in Middle Eastern and European dogs. The results of two-dimensional site frequency spectrum (2D SFS) support that the mutation is under soft sweep . We inferred that the onset of positive selection of the mutation is shorter than 6,535 years and behind the well-developed dairy economy in central Europe. It increases the expression of LCT by reducing its binding with ZEB1, which would enhance dog’s ability to digest milk-based diets. Our study uncovers the genetic basis of convergent evolution between humans and dogs with respect to diet, emphasizing the import of the dog as a biomedical model for studying mechanisms of the digestive system.     

To the Root of the Curl: A Signature of a Recent Selective Sweep Identifies a Mutation That Defines the Cornish Rex Cat Breed

The cat (Felis silvestris catus) shows significant variation in pelage, morphological, and behavioral phenotypes amongst its over 40 domesticated breeds. The majority of the breed specific phenotypic presentations originated through artificial selection, especially on desired novel phenotypic characteristics that arose only a few hundred years ago. Variations in coat texture and color of hair often delineate breeds amongst domestic animals. Although the genetic basis of several feline coat colors and hair lengths are characterized, less is known about the genes influencing variation in coat growth and texture, especially rexoid – curly coated types. Cornish Rex is a cat breed defined by a fixed recessive curly coat trait. Genome-wide analyses for selection (d_i, Tajima’s D and nucleotide diversity) were performed in the Cornish Rex breed and in 11 phenotypically diverse breeds and two random bred populations. Approximately 63K SNPs were used in the analysis that aimed to localize the locus controlling the rexoid hair texture. A region with a strong signature of recent selective sweep was identified in the Cornish Rex breed on chromosome A1, as well as a consensus block of homozygosity that spans approximately 3 Mb. Inspection of the region for candidate genes led to the identification of the lysophosphatidic acid receptor 6 (LPAR6). A 4 bp deletion in exon 5, c.250_253_delTTTG, which induces a premature stop codon in the receptor, was identified via Sanger sequencing. The mutation is fixed in Cornish Rex, absent in all straight haired cats analyzed, and is also segregating in the German Rex breed. LPAR6 encodes a G protein-coupled receptor essential for maintaining the structural integrity of the hair shaft; and has mutations resulting in a wooly hair phenotype in humans.     

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.     

Genetic Origins of Lactase Persistence and the Spread of Pastoralism in Africa

In humans, the ability to digest lactose, the sugar in milk, declines after weaning because of decreasing levels of the enzyme lactase-phlorizin hydrolase, encoded by LCT. However, some individuals maintain high enzyme amounts and are able to digest lactose into adulthood (i.e., they have the lactase-persistence [LP] trait). It is thought that selection has played a major role in maintaining this genetically determined phenotypic trait in different human populations that practice pastoralism. To identify variants associated with the LP trait and to study its evolutionary history in Africa, we sequenced MCM6 introns 9 and 13 and ∼2 kb of the LCT promoter region in 819 individuals from 63 African populations and in 154 non-Africans from nine populations. We also genotyped four microsatellites in an ∼198 kb region in a subset of 252 individuals to reconstruct the origin and spread of LP-associated variants in Africa. Additionally, we examined the association between LP and genetic variability at candidate regulatory regions in 513 individuals from eastern Africa. Our analyses confirmed the association between the LP trait and three common variants in intron 13 (C-14010, G-13907, and G-13915). Furthermore, we identified two additional LP-associated SNPs in intron 13 and the promoter region (G-12962 and T-956, respectively). Using neutrality tests based on the allele frequency spectrum and long-range linkage disequilibrium, we detected strong signatures of recent positive selection in eastern African populations and the Fulani from central Africa. In addition, haplotype analysis supported an eastern African origin of the C-14010 LP-associated mutation in southern Africa.     

HLA Alleles Influence the Clinical Signature of Amoxicillin-Clavulanate Hepatotoxicity

Background and Aim

The genotype-phenotype interaction in drug-induced liver injury (DILI) is a subject of growing interest. Previous studies have linked amoxicillin-clavulanate (AC) hepatotoxicity susceptibility to specific HLA alleles. In this study we aimed to examine potential associations between HLA class I and II alleles and AC DILI with regards to phenotypic characteristics, severity and time to onset in Spanish AC hepatotoxicity cases.

Methods

High resolution genotyping of HLA loci A, B, C, DRB1 and DQB1 was performed in 75 AC DILI cases and 885 controls.

Results

The distributions of class I alleles A*3002 (P/Pc = 2.6E-6/5E-5, OR 6.7) and B*1801 (P/Pc = 0.008/0.22, OR 2.9) were more frequently found in hepatocellular injury cases compared to controls. In addition, the presence of the class II allele combination DRB1*1501-DQB1*0602 (P/Pc = 5.1E-4/0.014, OR 3.0) was significantly increased in cholestatic/mixed cases. The A*3002 and/or B*1801 carriers were found to be younger (54 vs 65 years, P = 0.019) and were more frequently hospitalized than the DRB1*1501-DQB1*0602 carriers. No additional alleles outside those associated with liver injury patterns were found to affect potential severity as measured by Hy’s Law criteria. The phenotype frequencies of B*1801 (P/Pc = 0.015/0.42, OR 5.2) and DRB1*0301-DQB1*0201 (P/Pc = 0.0026/0.07, OR 15) were increased in AC DILI cases with delayed onset compared to those corresponding to patients without delayed onset, while the opposite applied to DRB1*1302-DQB1*0604 (P/Pc = 0.005/0.13, OR 0.07).

Conclusions

HLA class I and II alleles influence the AC DILI signature with regards to phenotypic expression, latency presentation and severity in Spanish patients.     

Persistence of Common Alleles in Two Related Populations or Species

Mathematical studies are conducted on three problems that arise in molecular population genetics. (1) The time required for a particular allele to become extinct in a population under the effects of mutation, selection, and random genetic drift is studied. In the absence of selection, the mean extinction time of an allele with an initial frequency close to 1 is of the order of the reciprocal of the mutation rate when 4Nv << 1, where N is the effective population size and v is the mutation rate per generation. Advantageous mutations reduce the extinction time considerably, whereas deleterious mutations increase it tremendously even if the effect on fitness is very slight. (2) Mathematical formulae are derived for the distribution and the moments of extinction time of a particular allele from one or both of two related populations or species under the assumption of no selection. When 4Nv << 1, the mean extinction time is about half that for a single population, if the two populations are descended from a common original stock. (3) The expected number as well as the proportion of common neutral alleles shared by two related species at the tth generation after their separation are studied. It is shown that if 4Nv is small, the two species are expected to share a high proportion of common alleles even 4N generations after separation. In addition to the above mathematical studies, the implications of our results for the common alleles at protein loci in related Drosophila species and for the degeneration of unused characters in cave animals are discussed.     

A Selective Sweep on a Deleterious Mutation in CPT1A in Arctic Populations

Arctic populations live in an environment characterized by extreme cold and the absence of plant foods for much of the year and are likely to have undergone genetic adaptations to these environmental conditions in the time they have been living there. Genome-wide selection scans based on genotype data from native Siberians have previously highlighted a 3 Mb chromosome 11 region containing 79 protein-coding genes as the strongest candidates for positive selection in Northeast Siberians. However, it was not possible to determine which of the genes might be driving the selection signal. Here, using whole-genome high-coverage sequence data, we identified the most likely causative variant as a nonsynonymous G>A transition (rs80356779; c.1436C>T [p.Pro479Leu] on the reverse strand) in CPT1A, a key regulator of mitochondrial long-chain fatty-acid oxidation. Remarkably, the derived allele is associated with hypoketotic hypoglycemia and high infant mortality yet occurs at high frequency in Canadian and Greenland Inuits and was also found at 68% frequency in our Northeast Siberian sample. We provide evidence of one of the strongest selective sweeps reported in humans; this sweep has driven this variant to high frequency in circum-Arctic populations within the last 6–23 ka despite associated deleterious consequences, possibly as a result of the selective advantage it originally provided to either a high-fat diet or a cold environment.     

The Yule Approximation for the Site Frequency Spectrum after a Selective Sweep

In the area of evolutionary theory, a key question is which portions of the genome of a species are targets of natural selection. Genetic hitchhiking is a theoretical concept that has helped to identify various such targets in natural populations. In the presence of recombination, a severe reduction in sequence diversity is expected around a strongly beneficial allele. The site frequency spectrum is an important tool in genome scans for selection and is composed of the numbers , where is the number of single nucleotide polymorphisms (SNPs) present in from individuals. Previous work has shown that both the number of low- and high-frequency variants are elevated relative to neutral evolution when a strongly beneficial allele fixes. Here, we follow a recent investigation of genetic hitchhiking using a marked Yule process to obtain an analytical prediction of the site frequency spectrum in a panmictic population at the time of fixation of a highly beneficial mutation. We combine standard results from the neutral case with the effects of a selective sweep. As simulations show, the resulting formula produces predictions that are more accurate than previous approaches for the whole frequency spectrum. In particular, the formula correctly predicts the elevation of low- and high-frequency variants and is significantly more accurate than previously derived formulas for intermediate frequency variants.