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.     

Independent introduction of two lactase-persistence alleles into human populations reflects different history of adaptation to milk culture

The T₋₁₃₉₁₀ variant located in the enhancer element of the lactase (LCT) gene correlates perfectly with lactase persistence (LP) in Eurasian populations whereas the variant is almost nonexistent among Sub-Saharan African populations, showing high prevalence of LP. Here, we report identification of two new mutations among Saudis, also known for the high prevalence of LP. We confirmed the absence of the European T₋₁₃₉₁₀ and established two new mutations found as a compound allele: T/G₋₁₃₉₁₅ within the −13910 enhancer region and a synonymous SNP in the exon 17 of the MCM6 gene T/C₋₃₇₁₂, −3712 bp from the LCT gene. The compound allele is driven to a high prevalence among Middle East population(s). Our functional analyses in vitro showed that both SNPs of the compound allele, located 10 kb apart, are required for the enhancer effect, most probably mediated through the binding of the hepatic nuclear factor 1 α (HNF1α). High selection coefficient (s) ~0.04 for LP phenotype was found for both T₋₁₃₉₁₀ and the compound allele. The European T₋₁₃₉₁₀ and the earlier identified East African G₋₁₃₉₀₇ LP allele share the same ancestral background and most likely the same history, probably related to the same cattle domestication event. In contrast, the compound Arab allele shows a different, highly divergent ancestral haplotype, suggesting that these two major global LP alleles have arisen independently, the latter perhaps in response to camel milk consumption. These results support the convergent evolution of the LP in diverse populations, most probably reflecting different histories of adaptation to milk culture.     

Association of Parkinson Disease with Structural and Regulatory Variants in the HLA Region

Historically, association of disease with the major histocompatibility complex (HLA) genes has been tested with HLA alleles that encode antigen-binding affinity. The association with Parkinson disease (PD), however, was discovered with noncoding SNPs in a genome-wide association study (GWAS). We show here that several HLA-region SNPs that have since been associated with PD form two blocks tagged by rs3129882 (p = 9 × 10⁻¹¹) and by rs9268515 and/or rs2395163 (p = 3 × 10⁻¹¹). We investigated whether these SNP-associations were driven by HLA-alleles at adjacent loci. We imputed class I and class II HLA-alleles for 2000 PD cases and 1986 controls from the NeuroGenetics Research Consortium GWAS and sequenced a subset of 194 cases and 204 controls. We were therefore able to assess accuracy of two imputation algorithms against next-generation-sequencing while taking advantage of the larger imputed data sets for disease study. Additionally, we imputed HLA alleles for 843 cases and 856 controls from another GWAS for replication. PD risk was positively associated with the B^∗07:02_C^∗07:02_DRB5^∗01_DRB1^∗15:01_DQA1^∗01:02_DQB1^∗06:02 haplotype and negatively associated with the C^∗03:04, DRB1^∗04:04 and DQA1^∗03:01 alleles. The risk haplotype and DQA1^∗03:01 lost significance when conditioned on the SNPs, but C^∗03:04 (OR = 0.72, p = 8 × 10⁻⁶) and DRB1^∗04:04 (OR = 0.65, p = 4 × 10⁻⁵) remained significant. Similarly, rs3129882 and the closely linked rs9268515 and rs2395163 remained significant irrespective of HLA alleles. rs3129882 and rs2395163 are expression quantitative trait loci (eQTLs) for HLA-DR and HLA-DQ (9 × 10⁻⁵ ≥ P_eQTL ≥ 2 × 10⁻⁷⁹), suggesting that HLA gene expression might influence PD. Our data suggest that PD is associated with both structural and regulatory elements in HLA. Furthermore, our study demonstrates that noncoding SNPs in the HLA region can be associated with disease irrespective of HLA alleles, and that observed associations with HLA alleles can sometimes be secondary to a noncoding variant.     

Identification of Multiple Genetic Susceptibility Loci in Takayasu Arteritis

Takayasu arteritis is a rare inflammatory disease of large arteries. The etiology of Takayasu arteritis remains poorly understood, but genetic contribution to the disease pathogenesis is supported by the genetic association with HLA-B^∗52. We genotyped ∼200,000 genetic variants in two ethnically divergent Takayasu arteritis cohorts from Turkey and North America by using a custom-designed genotyping platform (Immunochip). Additional genetic variants and the classical HLA alleles were imputed and analyzed. We identified and confirmed two independent susceptibility loci within the HLA region (r² < 0.2): HLA-B/MICA (rs12524487, OR = 3.29, p = 5.57 × 10⁻¹⁶) and HLA-DQB1/HLA-DRB1 (rs113452171, OR = 2.34, p = 3.74 × 10⁻⁹; and rs189754752, OR = 2.47, p = 4.22 × 10⁻⁹). In addition, we identified and confirmed a genetic association between Takayasu arteritis and the FCGR2A/FCGR3A locus on chromosome 1 (rs10919543, OR = 1.81, p = 5.89 × 10⁻¹²). The risk allele in this locus results in increased mRNA expression of FCGR2A. We also established the genetic association between IL12B and Takayasu arteritis (rs56167332, OR = 1.54, p = 2.18 × 10⁻⁸).     

Trans-ancestral fine-mapping of MHC reveals key amino acids associated with spontaneous clearance of hepatitis C in HLA-DQβ1

Spontaneous clearance of acute hepatitis C virus (HCV) infection is associated with single nucleotide polymorphisms (SNPs) on the MHC class II. We fine-mapped the MHC region in European (n = 1,600; 594 HCV clearance/1,006 HCV persistence) and African (n = 1,869; 340 HCV clearance/1,529 HCV persistence) ancestry individuals and evaluated HCV peptide binding affinity of classical alleles. In both populations, HLA-DQβ1Leu26 (p value_Meta = 1.24 × 10⁻¹⁴) located in pocket 4 was negatively associated with HCV spontaneous clearance and HLA-DQβ1Pro55 (p value_Meta = 8.23 × 10⁻¹¹) located in the peptide binding region was positively associated, independently of HLA-DQβ1Leu26. These two amino acids are not in linkage disequilibrium (r² < 0.1) and explain the SNPs and classical allele associations represented by rs2647011, rs9274711, HLA-DQB1^∗03:01, and HLA-DRB1^∗01:01. Additionally, HCV persistence classical alleles tagged by HLA-DQβ1Leu26 had fewer HCV binding epitopes and lower predicted binding affinities compared to clearance alleles (geometric mean of combined IC₅₀ nM of persistence versus clearance; 2,321 nM versus 761.7 nM, p value = 1.35 × 10⁻³⁸). In summary, MHC class II fine-mapping revealed key amino acids in HLA-DQβ1 explaining allelic and SNP associations with HCV outcomes. This mechanistic advance in understanding of natural recovery and immunogenetics of HCV might set the stage for much needed enhancement and design of vaccine to promote spontaneous clearance of HCV infection.     

HLA-DPB1 and HLA Class I Confer Risk of and Protection from Narcolepsy

Type 1 narcolepsy, a disorder caused by a lack of hypocretin (orexin), is so strongly associated with human leukocyte antigen (HLA) class II HLA-DQA1^∗01:02-DQB1^∗06:02 (DQ0602) that very few non-DQ0602 cases have been reported. A known triggering factor for narcolepsy is pandemic 2009 influenza H1N1, suggesting autoimmunity triggered by upper-airway infections. Additional effects of other HLA-DQ alleles have been reported consistently across multiple ethnic groups. Using over 3,000 case and 10,000 control individuals of European and Chinese background, we examined the effects of other HLA loci. After careful matching of HLA-DR and HLA-DQ in case and control individuals, we found strong protective effects of HLA-DPA1^∗01:03-DPB1^∗04:02 (DP0402; odds ratio [OR] = 0.51 [0.38–0.67], p = 1.01 × 10⁻⁶) and HLA-DPA1^∗01:03-DPB1^∗04:01 (DP0401; OR = 0.61 [0.47–0.80], p = 2.07 × 10⁻⁴) and predisposing effects of HLA-DPB1^∗05:01 in Asians (OR = 1.76 [1.34–2.31], p = 4.71 × 10⁻⁰⁵). Similar effects were found by conditional analysis controlling for HLA-DR and HLA-DQ with DP0402 (OR = 0.45 [0.38–0.55] p = 8.99 × 10⁻¹⁷) and DP0501 (OR = 1.38 [1.18–1.61], p = 7.11 × 10⁻⁵). HLA-class-II-independent associations with HLA-A^∗11:01 (OR = 1.32 [1.13–1.54], p = 4.92 × 10⁻⁴), HLA-B^∗35:03 (OR = 1.96 [1.41–2.70], p = 5.14 × 10⁻⁵), and HLA-B^∗51:01 (OR = 1.49 [1.25–1.78], p = 1.09 × 10⁻⁵) were also seen across ethnic groups in the HLA class I region. These effects might reflect modulation of autoimmunity or indirect effects of HLA class I and HLA-DP alleles on response to viral infections such as that of influenza.     

Two Susceptibility Loci to Takayasu Arteritis Reveal a Synergistic Role of the IL12B and HLA-B Regions in a Japanese Population

Takayasu arteritis (TAK) is an autoimmune systemic vasculitis of unknown etiology. Although previous studies have revealed that HLA-B^∗52:01 has an effect on TAK susceptibility, no other genetic determinants have been established so far. Here, we performed genome scanning of 167 TAK cases and 663 healthy controls via Illumina Infinium Human Exome BeadChip arrays, followed by a replication study consisting of 212 TAK cases and 1,322 controls. As a result, we found that the IL12B region on chromosome 5 (rs6871626, overall p = 1.7 × 10⁻¹³, OR = 1.75, 95% CI 1.42–2.16) and the MLX region on chromosome 17 (rs665268, overall p = 5.2 × 10⁻⁷, OR = 1.50, 95% CI 1.28–1.76) as well as the HLA-B region (rs9263739, a proxy of HLA-B^∗52:01, overall p = 2.8 × 10⁻²¹, OR = 2.44, 95% CI 2.03–2.93) exhibited significant associations. A significant synergistic effect of rs6871626 and rs9263739 was found with a relative excess risk of 3.45, attributable proportion of 0.58, and synergy index of 3.24 (p ≤ 0.00028) in addition to a suggestive synergistic effect between rs665268 and rs926379 (p ≤ 0.027). We also found that rs6871626 showed a significant association with clinical manifestations of TAK, including increased risk and severity of aortic regurgitation, a representative severe complication of TAK. Detection of these susceptibility loci will provide new insights to the basic mechanisms of TAK pathogenesis. Our findings indicate that IL12B plays a fundamental role on the pathophysiology of TAK in combination with HLA-B^∗52:01 and that common autoimmune mechanisms underlie the pathology of TAK and other autoimmune disorders such as psoriasis and inflammatory bowel diseases in which IL12B is involved as a genetic predisposing factor.     

Unraveling Multiple MHC Gene Associations with Systemic Lupus Erythematosus: Model Choice Indicates a Role for HLA Alleles and Non-HLA Genes in Europeans

We have performed a meta-analysis of the major-histocompatibility-complex (MHC) region in systemic lupus erythematosus (SLE) to determine the association with both SNPs and classical human-leukocyte-antigen (HLA) alleles. More specifically, we combined results from six studies and well-known out-of-study control data sets, providing us with 3,701 independent SLE cases and 12,110 independent controls of European ancestry. This study used genotypes for 7,199 SNPs within the MHC region and for classical HLA alleles (typed and imputed). Our results from conditional analysis and model choice with the use of the Bayesian information criterion show that the best model for SLE association includes both classical loci (HLA-DRB1^∗03:01, HLA-DRB1^∗08:01, and HLA-DQA1^∗01:02) and two SNPs, rs8192591 (in class III and upstream of NOTCH4) and rs2246618 (MICB in class I). Our approach was to perform a stepwise search from multiple baseline models deduced from a priori evidence on HLA-DRB1 lupus-associated alleles, a stepwise regression on SNPs alone, and a stepwise regression on HLA alleles. With this approach, we were able to identify a model that was an overwhelmingly better fit to the data than one identified by simple stepwise regression either on SNPs alone (Bayes factor [BF] > 50) or on classical HLA alleles alone (BF > 1,000).     

A Large-Scale Genetic Analysis Reveals a Strong Contribution of the HLA Class II Region to Giant Cell Arteritis Susceptibility

We conducted a large-scale genetic analysis on giant cell arteritis (GCA), a polygenic immune-mediated vasculitis. A case-control cohort, comprising 1,651 case subjects with GCA and 15,306 unrelated control subjects from six different countries of European ancestry, was genotyped by the Immunochip array. We also imputed HLA data with a previously validated imputation method to perform a more comprehensive analysis of this genomic region. The strongest association signals were observed in the HLA region, with rs477515 representing the highest peak (p = 4.05 × 10⁻⁴⁰, OR = 1.73). A multivariate model including class II amino acids of HLA-DRβ1 and HLA-DQα1 and one class I amino acid of HLA-B explained most of the HLA association with GCA, consistent with previously reported associations of classical HLA alleles like HLA-DRB1^∗04. An omnibus test on polymorphic amino acid positions highlighted DRβ1 13 (p = 4.08 × 10⁻⁴³) and HLA-DQα1 47 (p = 4.02 × 10⁻⁴⁶), 56, and 76 (both p = 1.84 × 10⁻⁴⁵) as relevant positions for disease susceptibility. Outside the HLA region, the most significant loci included PTPN22 (rs2476601, p = 1.73 × 10⁻⁶, OR = 1.38), LRRC32 (rs10160518, p = 4.39 × 10⁻⁶, OR = 1.20), and REL (rs115674477, p = 1.10 × 10⁻⁵, OR = 1.63). Our study provides evidence of a strong contribution of HLA class I and II molecules to susceptibility to GCA. In the non-HLA region, we confirmed a key role for the functional PTPN22 rs2476601 variant and proposed other putative risk loci for GCA involved in Th1, Th17, and Treg cell function.     

Genetic Variation in the 3'-Untranslated Region of NBN Gene Is Associated with Gastric Cancer Risk in a Chinese Population

NBN plays a crucial role in carcinogenesis as a core component for both homologous recombination (HR) and non-homologous end-joining (NHEJ) DNA double-strand breaks (DSBs) repair pathways. Genetic variants in the NBN gene have been associated with multiple cancers risk, suggesting pleiotropic effect on cancer. We hypothesized that genetic variants in the NBN gene may modify the risk of gastric cancer. To test this hypothesis, we evaluated the association between four potentially functional single nucleotide polymorphisms in NBN and gastric cancer risk in a case–control study of 1,140 gastric cancer cases and 1,547 controls in a Chinese population. We found that the A allele of rs10464867 (G>A) was significantly associated with a decreased risk of gastric cancer (odds ratio [OR] = 0.81, 95% confidence interval [95% CI] = 0.71–0.94; P = 4.71×10⁻³). Furthermore, the association between A allele of rs10464867 and decreased risk of gastric cancer was more significantly in elder individuals (per-allele OR = 0.72[0.59–0.88], P = 1.07×10⁻³), and male individuals (per-allele OR = 0.73[0.62–0.87], P = 3.68×10⁻⁴). We further conducted a haplotype analysis and identified that the NBN A_rs10464867G_rs14448G_rs1063053 haplotype conferred stronger protective effect on gastric cancer (OR = 0.76[0.65–0.89], P = 6.39×10⁻⁴). In summary, these findings indicate that genetic variants at NBN gene may contribute to gastric cancer susceptibility and may further advance our understanding of NBN gene in cancer development.     

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.     

Characterization of the major histocompatibility complex locus association with Beh?et’s disease in Iran

IntroductionThe aim of this study was to characterize the association of human leukocyte antigen (HLA) B alleles and major histocompatibility complex (MHC) single nucleotide polymorphisms (SNPs) with Behçet’s disease (BD) in an Iranian dataset.MethodsThe association of three SNPs in the MHC region previously identified as the most associated in high-density genotyping studies was tested in a case–control study on 973 BD patients and 825 controls from Iran, and the association of HLA-B alleles was tested in a subset of 681 patients and 414 controls.ResultsWe found that HLA-B*51 (P = 4.11 × 10⁻⁴¹, OR [95% CI] = 4.63[3.66-5.85]) and B*15 confer risk for BD (P = 2.83 × 10⁻², OR [95% CI] = 1.75[1.08-2.84]) in Iranian, and in B*51 negative individuals, only the B*15 allele is significantly associated with BD (P = 2.51 × 10⁻³, OR [95% CI] = 2.40[1.37-4.20]). rs76546355, formerly known as rs116799036, located between HLA-B and MICA (MHC class I polypeptide-related sequence A), demonstrated the same level of association with BD as HLA-B*51 (P_adj = 1.78 × 10⁻⁴⁶, OR [95% CI] = 5.46[4.21-7.09], and P_adj = 8.34 × 10⁻⁴⁸, OR [95% CI] = 5.44[4.20-7.05], respectively) in the HLA-B allelotyped subset, while rs2848713 was less associated (P_adj = 7.14 × 10⁻³⁵, OR [95% CI] = 3.73[2.97-4.69]) and rs9260997 was not associated (P_adj = 1.00 × 10⁻¹). Additionally, we found that B*51 genotype-phenotype correlations do not survive Bonferroni correction, while carriers of the rs76546355 risk allele predominate in BD cases with genital ulcers, positive pathergy test and positive BD family history (2.31 × 10⁻⁴ ≤ P ≤ 1.59 × 10⁻³).ConclusionsWe found that the HLA-B*51 allele and the rs76546355/rs116799036 MHC SNP are independent genetic risk factors for BD in Iranian, and that positivity for the rs76546355/rs116799036 risk allele, but not for B*51, does correlate with specific demographic characteristics or clinical manifestations in BD patients.

Electronic supplementary material

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

Stability or variation? Patterns of lactase gene and its enhancer region distributions in Brazilian Amerindians

Lactase persistence (LP) is the phenotypic trait in which lactase secretion is maintained during adulthood. LP is due to mutations in the LCT enhancer region, located 14‐kb upstream of the gene. In Europeans, the ?13910*T allele is associated with LP. In Africans this allele is rare while other mutations in this same region were related to LP. The LCT is highly polymorphic in human populations, but so far Brazilian Amerindians had not been investigated for these polymorphisms or for the presence of LP mutations. We describe the genetic diversity of the LCT region and the presence of LP enhancer mutations in four native Brazilian populations (Guarani‐Kaiowá, Guarani–Ñandeva, Kaingang, and Xavante). Twelve polymorphisms were genotyped by PCR‐based methods. The ?13910*T allele varied from 0.5% in the Xavante to 7.6% in the Guarani–Ñandeva. These frequencies probably derive from European sources and they correlate with non‐native admixture proportions previously estimated for these groups. But since admixture is virtually absent in the Xavante, we suggest that the presence of the LP allele could have been determined by a de novo mutation. No other mutations in the ?14 kb enhancer region were found. The LCT was highly polymorphic in the present sample showing 15 haplotypes with a heterogeneous distribution among the four Amerindian populations. This diversity could be due to drift, as indicated by the neutrality test performed. Am J Phys Anthropol, 2012. © 2012 Wiley Periodicals, Inc.     

Association of SNPs of CD40 Gene with Multiple Sclerosis in Russians

Multiple sclerosis (MS) is a serious, incurable neurological disease. In 2009, the ANZgene studies detected the suggestive association of located upstream of CD40 gene in chromosome 20q13 (p = 1.3×10⁻⁷). Identification of the causal variant(s) in the CD40 locus leads to a better understanding of the mechanism underlying the development of autoimmune pathologies. We determined the genotypes of rs6074022, rs1883832, rs1535045, and rs11086996 in patients with MS (n = 1684) and in the control group (n = 879). Two SNPs were significantly associated with MS: rs6074022 (additive model C allele OR = 1.27, 95% CI = [1.12–1.45], p = 3×10⁻⁴) and rs1883832 (additive model T allele OR = 1.20, 95% CI = [1.05–1.38], p = 7×10⁻³). In the meta-analysis of our results and the results of four previous studies, we obtain the association p-value of 2.34×10⁻¹², which confirmed the association between MS and rs6074022 at a genome-wide significant level. Next, we demonstrated that the model including rs6074022 only sufficiently described the association. From our analysis, we can speculate that the association between rs1883832 and MS was induced by LD, whereas rs6074022 was a marker in stronger LD with the functional variant or was the functional variant itself. Our results indicated that the functional variants were located in the upstream region of the gene CD40 and were in higher LD with rs6074022 than LD with rs1883832.