The hazards of genotype imputation in chromosomal regions under selection: A case study using the Lactase gene region.

Although imputation of missing SNP results has been widely used in genetic studies, claims about the quality and usefulness of imputation have outnumbered the few studies that have questioned its limitations. But it is becoming clear that these limitations are real-for example, disease association signals can be missed in regions of LD breakdown. Here, as a case study, using the chromosomal region of the well-known lactase gene, LCT, we address the issue of imputation in the context of variants that have become frequent in a limited number of modern population groups only recently, due to selection. We study SNPs in a 500 bp region covering the enhancer of LCT, and compare imputed genotypes with directly genotyped data. We examine the haplotype pairs of all individuals with discrepant and missing genotypes. We highlight the nonrandom nature of the allelic errors and show that most incorrect imputations and missing data result from long haplotypes that are evolutionarily closely related to those carrying the derived alleles, while some relate to rare and recombinant haplotypes. We conclude that bias of incorrectly imputed and missing genotypes can decrease the accuracy of imputed results substantially.

DNA methylation defines regional identity of human intestinal epithelial organoids and undergoes dynamic changes during development.

OBJECTIVE: Human intestinal epithelial organoids (IEOs) are increasingly being recognised as a highly promising translational research tool. However, our understanding of their epigenetic molecular characteristics and behaviour in culture remains limited. DESIGN: We performed genome-wide DNA methylation and transcriptomic profiling of human IEOs derived from paediatric/adult and fetal small and large bowel as well as matching purified human gut epithelium. Furthermore, organoids were subjected to in vitro differentiation and genome editing using CRISPR/Cas9 technology. RESULTS: We discovered stable epigenetic signatures which define regional differences in gut epithelial function, including induction of segment-specific genes during cellular differentiation. Established DNA methylation profiles were independent of cellular environment since organoids retained their regional DNA methylation over prolonged culture periods. In contrast to paediatric and adult organoids, fetal gut-derived organoids showed distinct dynamic changes of DNA methylation and gene expression in culture, indicative of an in vitro maturation. By applying CRISPR/Cas9 genome editing to fetal organoids, we demonstrate that this process is partly regulated by TET1, an enzyme involved in the DNA demethylation process. Lastly, generating IEOs from a child diagnosed with gastric heterotopia revealed persistent and distinct disease-associated DNA methylation differences, highlighting the use of organoids as disease-specific research models. CONCLUSIONS: Our study demonstrates striking similarities of epigenetic signatures in mucosa-derived IEOs with matching primary epithelium. Moreover, these results suggest that intestinal stem cell-intrinsic DNA methylation patterns establish and maintain regional gut specification and are involved in early epithelial development and disease.

Adaptation to milking agropastoralism in Chilean goat herders and nutritional benefit of lactase persistence.

The genetic trait of lactase persistence (LP) evolved as an adaptation to milking pastoralism in the Old World and is a well-known example of positive natural selection in humans. However, the specific mechanisms conferring this selective advantage are unknown. To understand the relationship between milk drinking, LP, growth, reproduction, and survival, communities of the Coquimbo Region in Chile, with recent adoption of milking agropastoralism, were used as a model population. DNA samples and data on stature, reproduction, and diet were collected from 451 participants. Lactose tolerance tests were done on 41 of them. The European -13,910*T (rs4988235) was the only LP causative variant found, showing strong association (99.6%) with LP phenotype. Models of associations of inferred LP status and milk consumption, with fertility, mortality, height, and weight were adjusted with measures of ancestry and relatedness to control for population structure. Although we found no statistically significant effect of LP on fertility, a significant effect (P = 0.002) was observed of LP on body mass index (BMI) in males and of BMI on fertility (P = 0.003). These results fail to support a causal relationship between LP and fertility yet suggest the idea of a nutritional advantage of LP. Furthermore, the proportion of European ancestry around the genetic region of -13,910*T is significantly higher (P = 0.008) than the proportion of European ancestry genome-wide, providing evidence of recent positive selection since European-Amerindian admixture. This signature was absent in nonpastoralist Latin American populations, supporting the hypothesis of specific adaptation to milking agropastoralism in the Coquimbo communities.

World-wide distributions of lactase persistence alleles and the complex effects of recombination and selection.

The genetic trait of lactase persistence (LP) is associated with at least five independent functional single nucleotide variants in a regulatory region about 14 kb upstream of the lactase gene [-13910*T (rs4988235), -13907*G (rs41525747), -13915*G (rs41380347), -14009*G (rs869051967) and -14010*C (rs145946881)]. These alleles have been inferred to have spread recently and present-day frequencies have been attributed to positive selection for the ability of adult humans to digest lactose without risk of symptoms of lactose intolerance. One of the inferential approaches used to estimate the level of past selection has been to determine the extent of haplotype homozygosity (EHH) of the sequence surrounding the SNP of interest. We report here new data on the frequencies of the known LP alleles in the 'Old World' and their haplotype lineages. We examine and confirm EHH of each of the LP alleles in relation to their distinct lineages, but also show marked EHH for one of the older haplotypes that does not carry any of the five LP alleles. The region of EHH of this (B) haplotype exactly coincides with a region of suppressed recombination that is detectable in families as well as in population data, and the results show how such suppression may have exaggerated haplotype-based measures of past selection.

In Vitro Functional Analyses of Infrequent Nucleotide Variants in the Lactase Enhancer Reveal Different Molecular Routes to Increased Lactase Promoter Activity and Lactase Persistence.

The genetic trait that allows intestinal lactase to persist into adulthood in some 35% of humans worldwide operates at the level of transcription, the effect being caused by cis-acting nucleotide changes upstream of the lactase gene (LCT). A single nucleotide substitution, -13910 C>T, the first causal variant to be identified, accounts for lactase persistence over most of Europe. Located in a region shown to have enhancer function in vitro, it causes increased activity of the LCT promoter in Caco-2 cells, and altered transcription factor binding. Three other variants in close proximity, -13907 C>G, -13915 T>C and -14010 G>C, were later shown to behave in a similar manner. Here, we study four further candidate functional variants. Two, -14009 T>G and -14011 C>T, adjacent to the well-studied -14010 G>C variant, also have a clear effect on promoter activity upregulation as assessed by transfection assays, but notably are involved in different molecular interactions. The results for the two other variants (-14028 T>C, -13779 G>C) were suggestive of function, -14028*C showing a clear change in transcription factor binding, but no obvious effect in transfections, while -13779*G showed greater effect in transfections but less on transcription factor binding. Each of the four variants arose on independent haplotypic backgrounds with different geographic distribution.

Diversity of lactase persistence alleles in Ethiopia: signature of a soft selective sweep.

The persistent expression of lactase into adulthood in humans is a recent genetic adaptation that allows the consumption of milk from other mammals after weaning. In Europe, a single allele (-13910(∗)T, rs4988235) in an upstream region that acts as an enhancer to the expression of the lactase gene LCT is responsible for lactase persistence and appears to have been under strong directional selection in the last 5,000 years, evidenced by the widespread occurrence of this allele on an extended haplotype. In Africa and the Middle East, the situation is more complicated and at least three other alleles (-13907(∗)G, rs41525747; -13915(∗)G, rs41380347; -14010(∗)C, rs145946881) in the same LCT enhancer region can cause continued lactase expression. Here we examine the LCT enhancer sequence in a large lactose-tolerance-tested Ethiopian cohort of more than 350 individuals. We show that a further SNP, -14009T>G (ss 820486563), is significantly associated with lactose-digester status, and in vitro functional tests confirm that the -14009(∗)G allele also increases expression of an LCT promoter construct. The derived alleles in the LCT enhancer region are spread through several ethnic groups, and we report a greater genetic diversity in lactose digesters than in nondigesters. By examining flanking markers to control for the effects of mutation and demography, we further describe, from empirical evidence, the signature of a soft selective sweep.

Diversity of lactase persistence in African milk drinkers.

The genetic trait of lactase persistence is attributable to allelic variants in an enhancer region upstream of the lactase gene, LCT. To date, five different functional alleles, -13910*T, -13907*G, -13915*G, -14009*G and -14010*C, have been identified. The co-occurrence of several of these alleles in Ethiopian lactose digesters leads to a pattern of sequence diversity characteristic of a 'soft selective sweep'. Here we hypothesise that throughout Africa, where multiple functional alleles co-exist, the enhancer diversity will be greater in groups who are traditional milk drinkers than in non-milk drinkers, as the result of this sort of parallel selection. Samples from 23 distinct groups from 10 different countries were examined. Each group was classified 'Yes 'or 'No' for milk-drinking, and ethnicity, language spoken and geographic location were recorded. Predicted lactase persistence frequency and enhancer diversity were, as hypothesised, higher in the milk drinkers than the non-milk-drinkers, but this was almost entirely accounted for by the Afro-Asiatic language speaking peoples of east Africa. The other groups, including the 'Nilo-Saharan language speaking' milk-drinkers, show lower frequencies of LP and lower diversity, and there was a north-east to south-west decline in overall diversity. Amongst the Afro-Asiatic (Cushitic) language speaking Oromo, however, the geographic cline was not evident and the southern pastoralist Borana showed much higher LP frequency and enhancer diversity than the other groups. Together these results reflect the effects of parallel selection, the stochastic processes of the occurrence and spread of the mutations, and time depth of milk drinking tradition.

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.

Dietary environmental factors shape the immune defense against Cryptosporidium infection.

Cryptosporidium is a leading cause of diarrheal-related deaths in children, especially in resource-poor settings. It also targets the immunocompromised, chronically infecting people living with HIV and primary immunodeficiencies. There is no vaccine or effective treatment. Although it is known from human cases and animal models that CD4+ T cells play a role in curbing Cryptosporidium, the role of CD8+ T cells remains to be defined. Using a Cryptosporidium tyzzeri mouse model, we show that gut-resident CD8+ intraepithelial lymphocytes (IELs) confer resistance to parasite growth. CD8+ IELs express and depend on the ligand-dependent transcription factor aryl hydrocarbon receptor (AHR). AHR deficiency reduces CD8+ IELs, decreases their cytotoxicity, and worsens infection. Transfer of CD8+ IELs rescues severely immunodeficient mice from death following Cryptosporidium challenge. Finally, dietary supplementation of the AHR pro-ligand indole-3-carbinol in newborn mice promotes resistance to infection. Therefore, common dietary metabolites augment the host immune response to cryptosporidiosis, protecting against disease.

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

Smarcad1 mediates microbiota-induced inflammation in mouse and coordinates gene expression in the intestinal epithelium.

BACKGROUND: How intestinal epithelial cells interact with the microbiota and how this is regulated at the gene expression level are critical questions. Smarcad1 is a conserved chromatin remodeling factor with a poorly understood tissue function. As this factor is highly expressed in the stem and proliferative zones of the intestinal epithelium, we explore its role in this tissue. RESULTS: Specific deletion of Smarcad1 in the mouse intestinal epithelium leads to colitis resistance and substantial changes in gene expression, including a striking increase of expression of several genes linked to innate immunity. Absence of Smarcad1 leads to changes in chromatin accessibility and significant changes in histone H3K9me3 over many sites, including genes that are differentially regulated upon Smarcad1 deletion. We identify candidate members of the gut microbiome that elicit a Smarcad1-dependent colitis response, including members of the poorly understood TM7 phylum. CONCLUSIONS: Our study sheds light onto the role of the chromatin remodeling machinery in intestinal epithelial cells in the colitis response and shows how a highly conserved chromatin remodeling factor has a distinct role in anti-microbial defense. This work highlights the importance of the intestinal epithelium in the colitis response and the potential of microbial species as pharmacological and probiotic targets in the context of inflammatory diseases.

Transcriptome analysis identifies a robust gene expression program in the mouse intestinal epithelium on aging.

The intestinal epithelium undergoes constant regeneration driven by intestinal stem cells. How old age affects the transcriptome in this highly dynamic tissue is an important, but poorly explored question. Using transcriptomics on sorted intestinal stem cells and adult enterocytes, we identified candidate genes, which change expression on aging. Further validation of these on intestinal epithelium of multiple middle-aged versus old-aged mice highlighted the consistent up-regulation of the expression of the gene encoding chemokine receptor Ccr2, a mediator of inflammation and several disease processes. We observed also increased expression of Strc, coding for stereocilin, and dramatically decreased expression of Rps4l, coding for a ribosome subunit. Ccr2 and Rps4l are located close to the telomeric regions of chromosome 9 and 6, respectively. As only few genes were differentially expressed and we did not observe significant protein level changes of identified ageing markers, our analysis highlights the overall robustness of murine intestinal epithelium gene expression to old age.

Microbiota derived short chain fatty acids promote histone crotonylation in the colon through histone deacetylases.

The recently discovered histone post-translational modification crotonylation connects cellular metabolism to gene regulation. Its regulation and tissue-specific functions are poorly understood. We characterize histone crotonylation in intestinal epithelia and find that histone H3 crotonylation at lysine 18 is a surprisingly abundant modification in the small intestine crypt and colon, and is linked to gene regulation. We show that this modification is highly dynamic and regulated during the cell cycle. We identify class I histone deacetylases, HDAC1, HDAC2, and HDAC3, as major executors of histone decrotonylation. We show that known HDAC inhibitors, including the gut microbiota-derived butyrate, affect histone decrotonylation. Consistent with this, we find that depletion of the gut microbiota leads to a global change in histone crotonylation in the colon. Our results suggest that histone crotonylation connects chromatin to the gut microbiota, at least in part, via short-chain fatty acids and HDACs.

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.