Leticia Plaza-Izurieta

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.

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 Pearsons 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.

Analysis of β-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.

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.

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.

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.

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.

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.

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 (R2) 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.

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.