The methylome of the celiac intestinal epithelium harbours genotype-independent alterations in the HLA region.

The Human Leucocyte Antigen (HLA) locus and other DNA sequence variants identified in Genome-Wide Association (GWA) studies explain around 50% of the heritability of celiac disease (CD). However, the pathogenesis of CD could be driven by other layers of genomic information independent from sequence variation, such as DNA methylation, and it is possible that allele-specific methylation explains part of the SNP associations. Since the DNA methylation landscape is expected to be different among cell types, we analyzed the methylome of the epithelial and immune cell populations of duodenal biopsies in CD patients and controls separately. We found a cell type-specific methylation signature that includes genes mapping to the HLA region, namely TAP1 and HLA-B. We also performed Immunochip SNP genotyping of the same samples and interrogated the expression of some of the affected genes. Our analysis revealed that the epithelial methylome is characterized by the loss of CpG island (CGI) boundaries, often associated to altered gene expression, and by the increased variability of the methylation across the samples. The overlap between differentially methylated positions (DMPs) and CD-associated SNPs or variants contributing to methylation quantitative trait loci (mQTLs) is minimal. In contrast, there is a notable enrichment of mQTLs among the most significant CD-associated SNPs. Our results support the notion that DNA methylation alterations constitute a genotype-independent event and confirm its role in the HLA region (apart from the well-known, DQ allele-specific effect). Finally, we find that a fraction of the CD-associated variants could exert its phenotypic effect through DNA methylation.

Transcription Factor Binding Site Enrichment Analysis in Co-Expression Modules in Celiac Disease.

The aim of this study was to construct celiac co-expression patterns at a whole genome level and to identify transcription factors (TFs) that could drive the gliadin-related changes in coordination of gene expression observed in celiac disease (CD). Differential co-expression modules were identified in the acute and chronic responses to gliadin using expression data from a previous microarray study in duodenal biopsies. Transcription factor binding site (TFBS) and Gene Ontology (GO) annotation enrichment analyses were performed in differentially co-expressed genes (DCGs) and selection of candidate regulators was performed. Expression of candidates was measured in clinical samples and the activation of the TFs was further characterized in C2BBe1 cells upon gliadin challenge. Enrichment analyses of the DCGs identified 10 TFs and five were selected for further investigation. Expression changes related to active CD were detected in four TFs, as well as in several of their in silico predicted targets. The activation of TFs was further characterized in C2BBe1 cells upon gliadin challenge, and an increase in nuclear translocation of CAMP Responsive Element Binding Protein 1 (CREB1) and IFN regulatory factor-1 (IRF1) in response to gliadin was observed. Using transcriptome-wide co-expression analyses we are able to propose novel genes involved in CD pathogenesis that respond upon gliadin stimulation, also in non-celiac models.

Association of the IL-15 and IL-15Rα genes with celiac disease.

Celiac disease is a chronic autoimmune condition triggered by dietary gluten in genetically predisposed individuals and the treatment is a strict gluten-free diet. The major predisposing genes are HLA-DQA1 and HLA-DQB1, but these are not sufficient for disease development. One of the candidate genes worth studying is interleukin (IL)-15 gene, together with its specific receptor, IL-15Rα, as they participate in promoting lymphocyte signaling and survival, and the establishment of appropriate conditions for villous atrophy, then acting as key players in the immunopathogenesis of CD. Here we analyze IL-15 and IL-15Rα genes in samples from the Spanish Consortium for Genetics of Celiac Disease (CEGEC) collection, identifying two regulatory single-nucleotide polymorphisms (SNP) that might be associated with celiac disease: rs4956400 (p-value 0.0112, OR 1.21, 95% CI 1.04-1.40) and rs11100722 (p-value 0.0087, OR 1.24, 95% CI 1.06-1.45), both located upstream the IL15 gene. When the expression of both genes was assessed, these two SNPs were found to be correlated with IL-15 higher protein expression. Besides, rs8177655 from IL15RA was also associated to mRNA IL-15 expression in CD patients. Finally, three SNPs from IL15RA intronic regions, rs2296141, rs3136614 and rs3181148, and another from its 3'UTR region, rs2229135, could be related to the age of diagnosis of celiac disease patients.

Ancestry-based stratified analysis of Immunochip data identifies novel associations with celiac disease.

To identify candidate genes in celiac disease (CD), we reanalyzed the whole Immunochip CD cohort using a different approach that clusters individuals based on immunoancestry prior to disease association analysis, rather than by geographical origin. We detected 636 new associated SNPs (P<7.02 × 10-07) and identified 5 novel genomic regions, extended 8 others previously identified and also detected 18 isolated signals defined by one or very few significant SNPs. To test whether we could identify putative candidate genes, we performed expression analyses of several genes from the top novel region (chr2:134533564-136169524), from a previously identified locus that is now extended, and a gene marked by an isolated SNP, in duodenum biopsies of active and treated CD patients, and non-celiac controls. In the largest novel region, CCNT2 and R3HDM1 were constitutively underexpressed in disease, even after gluten removal. Moreover, several genes within this region were coexpressed in patients, but not in controls. Other novel genes like KIF21B, REL and SORD also showed altered expression in active disease. Apart from the identification of novel CD loci, these results suggest that ancestry-based stratified analysis is an efficient strategy for association studies in complex diseases.

Expression analysis in intestinal mucosa reveals complex relations among genes under the association peaks in celiac disease.

Celiac disease is a chronic immune-mediated disorder with an important genetic component. To date, there are 57 independent association signals from 39 non-HLA loci, and a total of 66 candidate genes have been proposed. We aimed to scrutinize the functional implication of 45 of those genes by analyzing their expression in the disease tissue of celiac patients (at diagnosis/treatment) compared with non-celiac controls. Moreover, we investigated the SNP genotype effect in gene expression and performed coexpression analyses. Several genes showed differential expression among disease groups, most of them related to immune response. Multiple trans-eQTLs but only four cis-eQTLs were found, and surprisingly the genotype effect seems to be stimulus dependent as it differs among groups. Coexpression levels vary from higher to lower levels in active patients at diagnosis, treated patients and non-celiac controls respectively. A subset of 18 genes tightly correlated in both groups of patients but not in controls was identified. Interestingly, this subset of genes was influenced by the genotype of three SNPs. One of the SNPs, rs1018326 on chromosome two is on top of a known lincRNA whose function is not yet described, and whose expression seems to be upregulated in active disease when comparing biopsy pairs from the same individuals. Our results strongly suggest that the effects of disease-associated SNPs go far beyond the oversimplistic idea of transcriptional control at a nearby locus. Further investigations are needed to determine how each variant disrupts fine-tuning mechanisms in the genome that eventually lead to disease.

Alteration of tight junction gene expression in celiac disease.

OBJECTIVE: The aim of the present study was to characterize the deregulation of epithelial tight junction genes and investigate its reversibility on removal of dietary gluten in small intestinal mucosa in celiac disease (CD). METHODS: The expression levels of 23 genes related to tight junctions were studied in biopsies from 16 patients with active CD and compared with biopsies from the same patients taken after 2 years on gluten-free diet (GFD) and with 16 non-CD controls. RESULTS: Nine genes showed altered expression levels in patients with active disease (CLDN2, PARD6A, ZAK, SYMPK, MYH14, and ACTB were upregulated, whereas MAGI1, TJP1, and PPP2R3A were downregulated). Alterations were reversible after 2 years on treatment, except for PPP2R3A, implicated in the negative control of cell growth and division. At the biological network level, important dysfunctions in several processes within the pathway were observed, including intestinal permeability, apicobasal polarity, and cell proliferation. CONCLUSIONS: Our work confirms the involvement of tight junction genes related to permeability, polarity, and cell proliferation in the epithelial destruction observed in CD. Coexpression patterns of several genes support the idea of a common regulatory mechanism that seems to be altered in active CD. In general, GFD normalization confirms the reversibility of the process, except for the constitutive downregulation of PPP2R3A suggestive of a genetic implication. Further studies in proteins and cells or tissues are necessary to confirm these findings.

Detection of atomic scale changes in the free volume void size of three-dimensional colorectal cancer cell culture using positron annihilation lifetime spectroscopy.

Positron annihilation lifetime spectroscopy (PALS) provides a direct measurement of the free volume void sizes in polymers and biological systems. This free volume is critical in explaining and understanding physical and mechanical properties of polymers. Moreover, PALS has been recently proposed as a potential tool in detecting cancer at early stages, probing the differences in the subnanometer scale free volume voids between cancerous/healthy skin samples of the same patient. Despite several investigations on free volume in complex cancerous tissues, no positron annihilation studies of living cancer cell cultures have been reported. We demonstrate that PALS can be applied to the study in human living 3D cell cultures. The technique is also capable to detect atomic scale changes in the size of the free volume voids due to the biological responses to TGF-ß. PALS may be developed to characterize the effect of different culture conditions in the free volume voids of cells grown in vitro.

Coregulation and modulation of NFκB-related genes in celiac disease: uncovered aspects of gut mucosal inflammation.

It is known that the NFκB route is constitutively upregulated in celiac disease (CD), an immune-mediated disorder of the gut caused by intolerance to ingested gluten. Our aim was to scrutinize the expression patterns of several of the most biologically relevant components of the NFκB route in intestinal biopsies from active and treated patients and after in vitro gliadin challenge, and to assess normalization of the expression using an inhibitor of the MALT1 paracaspase. The expression of 93 NFκB genes was measured by RT-PCR in a set of uncultured active and treated CD and control biopsies, and in cultured biopsy series challenged with gliadin, the NFκB modulator, both compounds and none. Methylation of eight genes involved in NFκB signaling was analyzed by conventional pyrosequencing. Groups were compared and Pearson's correlation matrixes were constructed to check for coexpression and co-methylation. Our results confirm the upregulation of the NFκB pathway and show that constitutively altered genes usually belong to the core of the pathway and have central roles, whereas genes overexpressed only in active CD are more peripheral. Additionally, this is the first work to detect methylation level changes in celiac intestinal mucosa. Coexpression is very common in controls, whereas gliadin challenge and especially chronic inflammation present in untreated CD result in the disruption of the regulatory equilibrium. In contrast, co-methylation occurs more often in active CD. Importantly, NFκB modulation partially restores coregulation, opening the door to future therapeutic possibilities and targets.

THEMIS and PTPRK in celiac intestinal mucosa: coexpression in disease and after in vitro gliadin challenge.

Celiac disease (CD) is an immune mediated, polygenic disorder, where HLA-DQ2/DQ8 alleles contribute around 35% to genetic risk, but several other genes are also involved. Genome-wide association studies (GWASs) and the more recent immunochip genotyping projects have fine-mapped 39 regions of genetic susceptibility to the disease, most of which harbor candidate genes that could participate in this disease process. We focused our attention to the GWAS peak on chr6: 127.99-128.38 Mb, a region including two genes, thymocyte-expressed molecule involved in selection (THEMIS) and protein tyrosine phosphatase, receptor type, kappa (PTPRK), both of which have immune-related functions. The aim of this work was to evaluate the expression levels of these two genes in duodenal mucosa of active and treated CD patients and in controls, and to determine whether SNPs (rs802734, rs55743914, rs72975916, rs10484718 and rs9491896) associated with CD have any influence on gene expression. THEMIS showed higher expression in active CD compared with treated patients and controls, whereas PTPRK showed lower expression. Our study confirmed the association of this region with CD in our population, but only the genotype of rs802734 showed some influence in the expression of THEMIS. On the other hand, we found a significant positive correlation between THEMIS and PTPRK mRNA levels in CD patients but not in controls. Our results suggest a possible role for both candidate genes in CD pathogenesis and the existence of complex, regulatory relationships that reside in the vast non-coding, functional intergenic regions of the genome. Further investigation is needed to clarify the impact of the disease-associated SNPs on gene function.

Cubic regression-based degree of correction predicts the performance of whole bisulfitome amplified DNA methylation analysis.

Epigenetic mechanisms, including DNA methylation, are important determinants in development and disease. There is a need for technologies capable of detecting small variations in methylation levels in an accurate and reproducible manner, even if only limited amounts of DNA are available (which is the case in many studies in humans). Quantitative methylation analysis of minute DNA amounts after whole bisulfitome amplification (qMAMBA) has been proposed as an alternative, but this technique has not been adequately standardized and no comparative study against conventional methods has been performed, that includes a wide range of methylation percentages and different target assays. We designed an experiment to compare the performance of qMAMBA and bisulfite-treated genomic (non-amplified) DNA pyrosequencing. Reactions were performed in duplicate for each technique in eight different target genes, using nine artificially constructed DNA samples with methylation levels ranging between 0% and 100% with intervals of 12.5%. Cubic polynomial curves were plotted from the experimental results and the real methylation values and the resulting equation was used to estimate new corrected data points. The use of the cubic regression-based correction benefits the accuracy and the power of discrimination in methylation studies. Additionally, dispersion of the new estimated data around a y = x line (R ( 2) ) served to fix a cutoff that can discriminate, with a single 9-point curve experiment, whether whole bisulfitome amplification and subsequent qMAMBA can produce accurate methylation results. Finally, even with an optimized reagent kit, DNA samples subjected to whole bisulfitome amplification enhance the preferential amplification of unmethylated alleles, and subtle changes in methylation levels cannot be detected confidently.

Angiogenesis-related gene expression analysis in celiac disease.

Celiac disease (CD) involves disturbance of the small-bowel mucosal vascular network, and transglutaminase autoantibodies (TGA) have been related to angiogenesis disturbance, a complex phenomenon probably also influenced by common genetic variants in angiogenesis-related genes. A set of genes with "angiogenesis" GO term identified in a previous expression microarray experiment (SCG2, STAB1, TGFA, ANG, ERBB2, GNA13, PML, CASP8, ECGF1, JAG1, HIF1A, TNFSF13 and TGM2) was selected for genetic and functional studies. SNPs that showed a trend for association with CD in the first GWAS were genotyped in 555 patients and 541 controls. Gene expression of all genes was quantified in 15 pairs of intestinal biopsies (diagnosis vs. GFD) and in three-dimensional HUVEC and T84 cell cultures incubated with TGA-positive and negative serum. A regulatory SNP in TNFSF13 (rs11552708) is associated with CD (p = 0.01, OR = 0.7). Expression changes in biopsies pointed to TGM2 and PML as up-regulated antiangiogenic genes and to GNA13, TGFA, ERBB2 and SCG2 as down-regulated proangiogenic factors in CD. TGA seem to enhance TGM2 expression in both cell models, but PML expression was induced only in T84 enterocytes while GNA13 and ERBB2 were repressed in HUVEC endothelial cells, with several genes showing discordant effects in each model, highlighting the complexity of gene interactions in the pathogenesis of CD. Finally, cell culture models are useful tools to help dissect complex responses observed in human explants.

Accuracy in copy number calling by qPCR and PRT: a matter of DNA.

The possible implication of copy number variation (CNV) in the genetic susceptibility to human disease needs to be assessed using robust methods that can be applied at a population scale. In this report, we analyze the performance of the two major techniques, quantitative PCR (qPCR) and paralog ratio test (PRT), and investigate the influence of input DNA amount and template integrity on the reliability of both methods. Analysis of three genes (PRELID1, SYNPO and DEFB4) in a large sample set showed that both methods are prone to false copy number assignments if sufficient attention is not paid to DNA concentration and quality. Accurate normalization of samples is essential for reproducible qPCR because it avoids the effect of differential amplification efficiencies between target and control assays, whereas PRT is generally more sensitive to template degradation due to the fact that longer amplicons are usually needed to optimize sensitivity and specificity of paralog sequence PCR. The use of normalized, high quality genomic DNA yields comparable results with both methods.

Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease.

Using variants from the 1000 Genomes Project pilot European CEU dataset and data from additional resequencing studies, we densely genotyped 183 non-HLA risk loci previously associated with immune-mediated diseases in 12,041 individuals with celiac disease (cases) and 12,228 controls. We identified 13 new celiac disease risk loci reaching genome-wide significance, bringing the number of known loci (including the HLA locus) to 40. We found multiple independent association signals at over one-third of these loci, a finding that is attributable to a combination of common, low-frequency and rare genetic variants. Compared to previously available data such as those from HapMap3, our dense genotyping in a large sample collection provided a higher resolution of the pattern of linkage disequilibrium and suggested localization of many signals to finer scale regions. In particular, 29 of the 54 fine-mapped signals seemed to be localized to single genes and, in some instances, to gene regulatory elements. Altogether, we define the complex genetic architecture of the risk regions of and refine the risk signals for celiac disease, providing the next step toward uncovering the causal mechanisms of the disease.

Upregulation of KIR3DL1 gene expression in intestinal mucosa in active celiac disease.

Killer cell immunoglobulin-like receptors (KIRs) modulate natural killer (NK) and T-cell function by human leukocyte antigen class I interaction and have been implicated in celiac disease (CD). Qualitative expression of 16 KIR genes was determined in biopsies from 22 CD patients at diagnosis and after >2 years on a gluten-free diet (GFD). Quantitative expression analysis of KIR2DL4, KIR3DL1, KIR3DL3, and KLRC2 (a marker of an NK-reprogrammed T-cell subpopulation augmented in CD) was performed in 35 additional CD biopsy pairs and 14 non-CD control biopsies. No specific KIR expression profile was observed in CD. KIR3DL1 was more frequently expressed in active CD compared with GFD (p = 0.0312) and controls (p = 0.0008), with slightly increased levels in active disease. KLRC2 was overexpressed in active (p = 0.0037) and GFD (p = 0.0469) patients compared with non-CD controls and coexpressed with KIR3DL1. Results suggest the participation of KIR3DL1 overexpression in the overall immune activation seen in CD mucosa, which could be partly explained by the NK-like T-cell subpopulation increase.

Revisiting genome wide association studies (GWAS) in coeliac disease: replication study in Spanish population and expression analysis of candidate genes.

INTRODUCTION: Recent genome wide association studies (GWAS) on coeliac disease (CD) have identified risk loci harbouring genes that fit the accepted pathogenic model and are considered aetiological candidates. METHODS: Using Taqman single nucleotide polymorphism (SNP) and expression assays, the study genotyped 11 SNPs tagging eight GWAS regions (1q31, 2q11-2q12, 3p21, 3q25-3q26, 3q28, 4q27, 6q25 and 12q24) in a Spanish cohort of 1094 CD patients and 540 controls, and performed expression analyses of candidate genes (RGS1, IL18R1/IL18RAP, CCR3, IL12A/SCHIP1, LPP, IL2/IL21-KIAA1109, TAGAP, and SH2B3) in intestinal mucosa from 29 CD children and eight controls. RESULTS: Polymorphisms in 1q31, 2q11-2q12, and 3q25 showed association in our cohort, and also 3q28 and 4q27 when combined with a previous study. Expression levels of IL12A, IL18RAP, IL21, KIAA1109, LPP, SCHIP1, and SH2B3 were affected by disease status, but the correlation between genotype and mRNA levels was observed only in IL12A, LPP, SCHIP1, and SH2B3. CONCLUSIONS: Expression differences between treated CD patients and controls along with SNP expression associations suggest a possible primary role for these four genes and their variants in pathogenesis. The lack of SNP effect in the remaining genes is probably a consequence of arbitrary candidate gene selection within association signals that are not based on functional studies.

Analysis of beta-defensin and Toll-like receptor gene copy number variation in celiac disease.

Celiac disease (CD) is an immune-mediated disorder of the gut in which innate and adaptive responses are involved. Toll-like receptor (TLR) 2 and TLR4 participate in host defense through antigen recognition, and show altered expression in CD gut mucosa. beta-defensins are inducible antimicrobial peptides, and DEFB gene copy number polymorphisms have been associated with autoimmune and inflammatory disorders. We performed copy number analysis of TLR2, TLR4, and the beta-defensin cluster (DEFB4, DEFB103 and DEFB104) by gene-specific, real-time polymerase chain reaction (PCR) in 376 CD patients and 376 controls. TLR genes did not show copy number variation, and all samples presented with two copies. beta-defensin clusters varied between 2 and 9 copies per genome, and when grouped into bins, high copy numbers (>4) were underrepresented among patients (p = 0.023; odds ratio = 0.69, 95% CI = 0.50-0.96), suggesting that increased copy numbers could protect from CD, possibly by impeding bacterial infiltration more efficiently and preserving gut epithelial integrity.