Identification of Non-HLA Genes Associated with Celiac Disease and Country-Specific Differences in a Large,International Pediatric Cohort

Objectives

There are significant geographical differences in the prevalence and incidence of celiac disease that cannot be explained by HLA alone. More than 40 loci outside of the HLA region have been associated with celiac disease. We investigated the roles of these non-HLA genes in the development of tissue transglutaminase autoantibodies (tTGA) and celiac disease in a large international prospective cohort study.

Methods

A total of 424,788 newborns from the US and European general populations and first-degree relatives with type 1 diabetes were screened for specific HLA genotypes. Of these, 21,589 carried 1 of the 9 HLA genotypes associated with increased risk for type 1 diabetes and celiac disease; we followed 8676 of the children in a 15 y prospective follow-up study. Genotype analyses were performed on 6010 children using the Illumina ImmunoChip. Levels of tTGA were measured in serum samples using radio-ligand binding assays; diagnoses of celiac disease were made based on persistent detection of tTGA and biopsy analysis. Data were analyzed using Cox proportional hazards analyses.

Results

We found 54 single-nucleotide polymorphisms (SNPs) in 5 genes associated with celiac disease (TAGAP, IL18R1, RGS21, PLEK, and CCR9) in time to celiac disease analyses (10⁻⁴>P>5.8x10⁻⁶). The hazard ratios (HR) for the SNPs with the smallest P values in each region were 1.59, 1.45, 2.23, 2.64, and 1.40, respectively. Outside of regions previously associated with celiac disease, we identified 10 SNPs in 8 regions that could also be associated with the disease (P<10⁻⁴). A SNP near PKIA (rs117128341, P = 6.5x10⁻⁸, HR = 2.8) and a SNP near PFKFB3 (rs117139146, P<2.8x10⁻⁷, HR = 4.9) reached the genome-wide association threshold in subjects from Sweden. Analyses of time to detection of tTGA identified 29 SNPs in 2 regions previously associated with celiac disease (CTLA4, P = 1.3x10⁻⁶, HR = 0.76 and LPP, P = 2.8x10⁻⁵, HR = .80) and 6 SNPs in 5 regions not previously associated with celiac disease (P<10⁻⁴); non-HLA genes are therefore involved in development of tTGA.

Conclusions

In conclusion, using a genetic analysis of a large international cohort of children, we associated celiac disease development with 5 non-HLA regions previously associated with the disease and 8 regions not previously associated with celiac disease. We identified 5 regions associated with development of tTGA. Two loci associated with celiac disease progression reached a genome-wide association threshold in subjects from Sweden.     

Electron-transfer rates govern product distribution in electrochemically-driven P450-catalyzed dioxygen reduction

Developing electrode-driven biocatalytic systems utilizing the P450 cytochromes for selective oxidations depends not only on achieving electron transfer (ET) but also doing so at rates that favor native-like turnover. Herein we report studies that correlate rates of heme reduction with ET pathways and resulting product distributions. We utilized single-surface cysteine mutants of the heme domain of P450 from Bacillus megaterium and modified the thiols with N-(1-pyrene)-iodoacetamide, affording proteins that could bond to basal-plane graphite. Of the proteins examined, Cys mutants at position 62, 383, and 387 were able to form electroactive monolayers with similar E_1/2 values (− 335 to − 340 mV vs AgCl/Ag). Respective ET rates (k_s^o) and heme-cysteine distances for 62, 383, and 387 are 50 s^-1 and 16 ?, 0.8 s^- 1 and 25 ?, and 650 s^- 1 and 19 ?. Experiments utilizing rotated-disk electrodes were conducted to determine the products of P450-catalyzed dioxygen reduction. We found good agreement between ET rates and product distributions for the various mutants, with larger k_s^o values correlating with more electrons transferred per dioxygen during catalysis.     

Identification of CAD candidate genes in GWAS loci and their expression in vascular cells

Recent genome-wide association studies (GWAS) have identified 35 loci that significantly associate with coronary artery disease (CAD) susceptibility. The majority of the genes represented in these loci have not previously been studied in the context of atherosclerosis. To characterize the roles of these candidate genes in the vessel wall, we determined their expression levels in endothelial, smooth muscle, and macrophage cells isolated from healthy, prelesioned, and lesioned mouse aortas. We also performed expression quantitative locus (eQTL) mapping of these genes in human endothelial cells under control and proatherogenic conditions. Of the 57 genes studied, 31 were differentially expressed in one or more cell types in disease state in mice, and the expression levels of 8 were significantly associated with the CAD SNPs in human cells, 7 of which were also differentially expressed in mice. By integrating human and mouse results, we predict that PPAP2B, GALNT4, MAPKAPK5, TCTN1, SRR, SNF8, and ICAM1 play a causal role in the susceptibility to atherosclerosis through a role in the vasculature. Additionally, we highlight the genetic complexity of a subset of CAD loci through the differential expression of multiple candidate genes per locus and the involvement of genes that lie outside linkage disequilibrium blocks.