Evolutionary and Functional Analysis of Celiac Risk Loci Reveals SH2B3 as a Protective Factor against Bacterial Infection

Celiac disease (CD) is an intolerance to dietary proteins of wheat, barley, and rye. CD may have substantial morbidity, yet it is quite common with a prevalence of 1%–2% in Western populations. It is not clear why the CD phenotype is so prevalent despite its negative effects on human health, especially because appropriate treatment in the form of a gluten-free diet has only been available since the 1950s, when dietary gluten was discovered to be the triggering factor. The high prevalence of CD might suggest that genes underlying this disease may have been favored by the process of natural selection. We assessed signatures of selection for ten confirmed CD-associated loci in several genome-wide data sets, comprising 8154 controls from four European populations and 195 individuals from a North African population, by studying haplotype lengths via the integrated haplotype score (iHS) method. Consistent signs of positive selection for CD-associated derived alleles were observed in three loci: IL12A, IL18RAP, and SH2B3. For the SH2B3 risk allele, we also show a difference in allele frequency distribution (F_st) between HapMap phase II populations. Functional investigation of the effect of the SH2B3 genotype in response to lipopolysaccharide and muramyl dipeptide revealed that carriers of the SH2B3 rs3184504^∗A risk allele showed stronger activation of the NOD2 recognition pathway. This suggests that SH2B3 plays a role in protection against bacteria infection, and it provides a possible explanation for the selective sweep on SH2B3, which occurred sometime between 1200 and 1700 years ago.     

Improving Phenotypic Prediction by Combining Genetic and Epigenetic Associations

We tested whether DNA-methylation profiles account for inter-individual variation in body mass index (BMI) and height and whether they predict these phenotypes over and above genetic factors. Genetic predictors were derived from published summary results from the largest genome-wide association studies on BMI (n ∼ 350,000) and height (n ∼ 250,000) to date. We derived methylation predictors by estimating probe-trait effects in discovery samples and tested them in external samples. Methylation profiles associated with BMI in older individuals from the Lothian Birth Cohorts (LBCs, n = 1,366) explained 4.9% of the variation in BMI in Dutch adults from the LifeLines DEEP study (n = 750) but did not account for any BMI variation in adolescents from the Brisbane Systems Genetic Study (BSGS, n = 403). Methylation profiles based on the Dutch sample explained 4.9% and 3.6% of the variation in BMI in the LBCs and BSGS, respectively. Methylation profiles predicted BMI independently of genetic profiles in an additive manner: 7%, 8%, and 14% of variance of BMI in the LBCs were explained by the methylation predictor, the genetic predictor, and a model containing both, respectively. The corresponding percentages for LifeLines DEEP were 5%, 9%, and 13%, respectively, suggesting that the methylation profiles represent environmental effects. The differential effects of the BMI methylation profiles by age support previous observations of age modulation of genetic contributions. In contrast, methylation profiles accounted for almost no variation in height, consistent with a mainly genetic contribution to inter-individual variation. The BMI results suggest that combining genetic and epigenetic information might have greater utility for complex-trait prediction.     

The long-term genetic stability and individual specificity of the human gut microbiome

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

     

New Gene Variants Associated with the Risk of Chronic HBV Infection

Some patients with chronic hepatitis B virus(HBV) infection failed to clear HBV, even persistently continue to produce antibodies to HBV. Here we performed a two stage genome wide association study in a cohort of Chinese patients designed to discover single nucleotide variants that associate with HBV infection and clearance of HBV. The first stage involved genome wide exome sequencing of 101 cases(HBsAg plus anti-HBs positive) compared with 102 control patients(antiHBs positive, HBsAg negative). Over 80% of individual sequences displayed 209 sequence coverage. Adapters,uncertain bases [10% or low-quality base calls([50%) were filtered and compared to the human reference genome hg19. In the second stage, 579 chronic HBV infected cases and 439 HBV clearance controls were sequenced with selected genes from the first stage. Although there were no significant associated gene variants in the first stage, two significant gene associations were discovered when the two stages were assessed in a combined analysis. One association showed rs506121-‘‘T' allele [within the dedicator of cytokinesis 8(DOCK8) gene] was higher in chronic HBV infection group than that in clearance group(P = 0.002, OR = 0.77, 95% CI [0.65, 0.91]). The second association involved rs2071676—A allele within the Carbonic anhydrase(CA9) gene that was significantly elevated in chronic HBV infection group compared to the clearance group(P = 0.0003, OR = 1.35, 95% CI [1.15, 1.58]). Upon replication these gene associations would suggest the influence of DOCK8 and CA9 as potential risk genetic factors in the persistence of HBV infection.     

Novel association in chromosome 4q27 region with rheumatoid arthritis and confirmation of type 1 diabetes point to a general risk locus for autoimmune diseases

Recently, association of celiac disease with common single-nucleotide polymorphism (SNP) variants in an extensive linkage-disequilibrium block of 480 kb containing the KIAA1109, Tenr, IL2, and IL21 genes has been demonstrated in three independent populations (rs6822844 P_combined=1.3×10^-14). The KIAA1109/Tenr/IL2/IL21 block corresponds to the Idd3 locus in the nonobese diabetic mouse model of type 1 diabetes (T1D). This block was recently found to be associated with T1D in a genomewide association study, although this finding lacks unequivocal confirmation. We therefore aimed to investigate whether the KIAA1109/Tenr/IL2/IL21 region is involved in susceptibility to multiple autoimmune diseases. We tested SNP rs6822844 for association with disease in 350 T1D-affected and 1,047 rheumatoid arthritis (RA)–affected Dutch patients and in 929 controls. We replicated the association with T1D (P=.0006; OR 0.64 [95% CI 0.50–0.83]), and revealed a similar novel association with RA (P=.0002; OR 0.72 [95% CI 0.61–0.86]). Our results replicate and extend the association found in the KIAA1109/Tenr/IL2/IL21 gene region with autoimmune diseases, implying that this locus is a general risk factor for multiple autoimmune diseases.     

Genes involved in the regulation of vascular homeostasis determine renal survival rate in patients with chronic glomerulonephritis

Chronic glomerulonephritis (CGN) is one of the most severe kidney diseases. Genes of vascular reactivity are thought to play an important role in development and progression of CGN. In this study, we analyzed association of genes of vascular homeostasis with hypertension and renal survival of CGN patients. The study sample included 238 patients with CGN and 304 healthy subjects of population control. Ten polymorphisms of ten genes of vascular homeostasis were genotyped through polymerase chain reaction (PCR), restriction fragment length polymorphism (RFLP) analysis and TaqMan assays. Association of the genotypes with renal survival was analyzed by the Kaplan–Meier estimator. Genotypes 311SC and 311SS of the PON2 gene, (− 1166)AC and (− 1166)CC of the AGTR1 gene, (+ 46)AA of the ADRB2 gene, and 198KK and 198KN of the EDN1 gene were associated with decreased rate of renal survival of the patients. Polymorphisms S311C PON2, (− 1166)A/C AGTR1, (+ 46)G/A ADRB2, and K198N EDN1 were associated with the accelerated decline in kidney function in the CGN patients.     

Fine Mapping Seronegative and Seropositive Rheumatoid Arthritis to Shared and Distinct HLA Alleles by Adjusting for the Effects of Heterogeneity

Despite progress in defining human leukocyte antigen (HLA) alleles for anti-citrullinated-protein-autoantibody-positive (ACPA⁺) rheumatoid arthritis (RA), identifying HLA alleles for ACPA-negative (ACPA⁻) RA has been challenging because of clinical heterogeneity within clinical cohorts. We imputed 8,961 classical HLA alleles, amino acids, and SNPs from Immunochip data in a discovery set of 2,406 ACPA⁻ RA case and 13,930 control individuals. We developed a statistical approach to identify and adjust for clinical heterogeneity within ACPA⁻ RA and observed independent associations for serine and leucine at position 11 in HLA-DRβ1 (p = 1.4 × 10⁻¹³, odds ratio [OR] = 1.30) and for aspartate at position 9 in HLA-B (p = 2.7 × 10⁻¹², OR = 1.39) within the peptide binding grooves. These amino acid positions induced associations at HLA-DRB1^∗03 (encoding serine at 11) and HLA-B^∗08 (encoding aspartate at 9). We validated these findings in an independent set of 427 ACPA⁻ case subjects, carefully phenotyped with a highly sensitive ACPA assay, and 1,691 control subjects (HLA-DRβ1 Ser11+Leu11: p = 5.8 × 10⁻⁴, OR = 1.28; HLA-B Asp9: p = 2.6 × 10⁻³, OR = 1.34). Although both amino acid sites drove risk of ACPA⁺ and ACPA⁻ disease, the effects of individual residues at HLA-DRβ1 position 11 were distinct (p < 2.9 × 10⁻¹⁰⁷). We also identified an association with ACPA⁺ RA at HLA-A position 77 (p = 2.7 × 10⁻⁸, OR = 0.85) in 7,279 ACPA⁺ RA case and 15,870 control subjects. These results contribute to mounting evidence that ACPA⁺ and ACPA⁻ RA are genetically distinct and potentially have separate autoantigens contributing to pathogenesis. We expect that our approach might have broad applications in analyzing clinical conditions with heterogeneity at both major histocompatibility complex (MHC) and non-MHC regions.