Summary-data based Mendelian randomization identifies gene expression regulatory polymorphisms associated with bovine paratuberculosis by modulation of the nuclear factor Kappa ß (NF-κß)-mediated inflammatory response.

Genome-wide association studies (GWAS) have identified host genetic variants associated with paratuberculosis (PTB) susceptibility. Most of the GWAS-identified SNPs are in non-coding regions. Connecting these non-coding variants and downstream affected genes is a challenge and, up to date, only a few functional mutations or expression quantitative loci (cis-eQTLs) associated with PTB susceptibility have been identified. In the current study, the associations between imputed whole-genome sequence genotypes and whole RNA-Sequencing data from peripheral blood (PB) and ileocecal valve (ICV) samples of Spanish Holstein cows (N = 16) were analyzed with TensorQTL. This approach allowed the identification of 88 and 37 cis-eQTLs regulating the expression levels of 90 and 37 genes in PB and ICV samples, respectively (False discorey rate, FDR ≤ 0.05). Next, we applied summary-based data Mendelian randomization (SMR) to integrate the cis-eQTL dataset with GWAS data obtained from a cohort of 813 culled cattle that were classified according to the presence or absence of PTB-associated histopathological lesions in gut tissues. After multiple testing corrections (FDR ≤ 0.05), we identified two novel cis-eQTLs affecting the expression of the early growth response factor 4 (EGR4) and the bovine neuroblastoma breakpoint family member 6-like protein isoform 2 (MGC134040) that showed pleiotropic associations with the presence of multifocal and diffuse lesions in gut tissues; P = 0.002 and P = 0.017, respectively. While EGR4 acts as a brake on T-cell proliferation and cytokine production through interaction with the nuclear factor Kappa ß (NF-κß), MGC134040 is a target gene of NF-κß. Our findings provide a better understanding of the genetic factors influencing PTB outcomes, confirm that the multifocal lesions are localized/confined lesions that have different underlying host genetics than the diffuse lesions, and highlight regulatory SNPs and regulated-gene targets to design future functional studies.

A meta-analysis of pre-pregnancy maternal body mass index and placental DNA methylation identifies 27 CpG sites with implications for mother-child health.

Higher maternal pre-pregnancy body mass index (ppBMI) is associated with increased neonatal morbidity, as well as with pregnancy complications and metabolic outcomes in offspring later in life. The placenta is a key organ in fetal development and has been proposed to act as a mediator between the mother and different health outcomes in children. The overall aim of the present work is to investigate the association of ppBMI with epigenome-wide placental DNA methylation (DNAm) in 10 studies from the PACE consortium, amounting to 2631 mother-child pairs. We identify 27 CpG sites at which we observe placental DNAm variations of up to 2.0% per 10 ppBMI-unit. The CpGs that are differentially methylated in placenta do not overlap with CpGs identified in previous studies in cord blood DNAm related to ppBMI. Many of the identified CpGs are located in open sea regions, are often close to obesity-related genes such as GPX1 and LGR4 and altogether, are enriched in cancer and oxidative stress pathways. Our findings suggest that placental DNAm could be one of the mechanisms by which maternal obesity is associated with metabolic health outcomes in newborns and children, although further studies will be needed in order to corroborate these findings.

Prenatal Maternal Smoke, DNA Methylation, and Multi-omics of Tissues and Child Health.

PURPOSE OF REVIEW: Maternal tobacco smoking during pregnancy is of public health concern, and understanding the biological mechanisms can help to promote smoking cessation campaigns. This non-systematic review focuses on the effects of maternal smoking during pregnancy on offspring's epigenome, consistent in chemical modifications of the genome that regulate gene expression. RECENT FINDINGS: Recent meta-analyses of epigenome-wide association studies have shown that maternal smoking during pregnancy is consistently associated with offspring's DNA methylation changes, both in the placenta and blood. These studies indicate that effects on blood DNA methylation can persist for years, and that the longer the duration of the exposure and the higher the dose, the larger the effects. Hence, DNA methylation scores have been developed to estimate past exposure to maternal smoking during pregnancy as biomarkers. There is robust evidence for DNA methylation alterations associated with maternal smoking during pregnancy; however, the role of sex, ethnicity, and genetic background needs further exploration. Moreover, there are no conclusive studies about exposure to low doses or during the preconception period. Similarly, studies on tissues other than the placenta and blood are scarce, and cell-type specificity within tissues needs further investigation. In addition, biological interpretation of DNA methylation findings requires multi-omics data, poorly available in epidemiological settings. Finally, although several mediation analyses link DNA methylation changes with health outcomes, they do not allow causal inference. For this, a combination of data from multiple study designs will be essential in the future to better address this topic.

Placental DNA methylation signatures of maternal smoking during pregnancy and potential impacts on fetal growth.

Maternal smoking during pregnancy (MSDP) contributes to poor birth outcomes, in part through disrupted placental functions, which may be reflected in the placental epigenome. Here we present a meta-analysis of the associations between MSDP and placental DNA methylation (DNAm) and between DNAm and birth outcomes within the Pregnancy And Childhood Epigenetics (PACE) consortium (N = 1700, 344 with MSDP). We identify 443 CpGs that are associated with MSDP, of which 142 associated with birth outcomes, 40 associated with gene expression, and 13 CpGs are associated with all three. Only two CpGs have consistent associations from a prior meta-analysis of cord blood DNAm, demonstrating substantial tissue-specific responses to MSDP. The placental MSDP-associated CpGs are enriched for environmental response genes, growth-factor signaling, and inflammation, which play important roles in placental function. We demonstrate links between placental DNAm, MSDP and poor birth outcomes, which may better inform the mechanisms through which MSDP impacts placental function and fetal growth.

Genetic Contribution of Endometriosis to the Risk of Developing Hormone-Related Cancers.

Endometriosis is a common gynecological disorder that has been associated with endometrial, breast and epithelial ovarian cancers in epidemiological studies. Since complex diseases are a result of multiple environmental and genetic factors, we hypothesized that the biological mechanism underlying their comorbidity might be explained, at least in part, by shared genetics. To assess their potential genetic relationship, we performed a two-sample mendelian randomization (2SMR) analysis on results from public genome-wide association studies (GWAS). This analysis confirmed previously reported genetic pleiotropy between endometriosis and endometrial cancer. We present robust evidence supporting a causal genetic association between endometriosis and ovarian cancer, particularly with the clear cell and endometrioid subtypes. Our study also identified genetic variants that could explain those associations, opening the door to further functional experiments. Overall, this work demonstrates the value of genomic analyses to support epidemiological data, and to identify targets of relevance in multiple disorders.

Identification of loci associated with susceptibility to bovine paratuberculosis and with the dysregulation of the MECOM, eEF1A2, and U1 spliceosomal RNA expression.

Although genome-wide association studies have identified single nucleotide polymorphisms (SNPs) associated with the susceptibility to Mycobacterium avium subsp. paratuberculosis (MAP) infection, only a few functional mutations for bovine paratuberculosis (PTB) have been characterized. Expression quantitative trait loci (eQTLs) are genetic variants typically located in gene regulatory regions that alter gene expression in an allele-specific manner. eQTLs can be considered as functional links between genomic variants, gene expression, and ultimately phenotype. In the current study, peripheral blood (PB) and ileocecal valve (ICV) gene expression was quantified by RNA-Seq from fourteen Holstein cattle with no lesions and with PTB-associated histopathological lesions in gut tissues. Genotypes were generated from the Illumina LD EuroG10K BeadChip. The associations between gene expression levels (normalized read counts) and genetic variants were analyzed by a linear regression analysis using R Matrix eQTL 2.2. This approach allowed the identification of 192 and 48 cis-eQTLs associated with the expression of 145 and 43 genes in the PB and ICV samples, respectively. To investigate potential relationships between these cis-eQTLs and MAP infection, a case-control study was performed using the genotypes for all the identified cis-eQTLs and phenotypical data (histopathology, ELISA for MAP-antibodies detection, tissue PCR, and bacteriological culture) of 986 culled cows. Our results suggested that the heterozygous genotype in the cis-eQTL-rs43744169 (T/C) was associated with the up-regulation of the MDS1 and EVI1 complex (MECOM) expression, with positive ELISA, PCR, and bacteriological culture results, and with increased risk of progression to clinical PTB. As supporting evidence, the presence of the minor allele was associated with higher MECOM levels in plasma samples from infected cows and with increased MAP survival in an ex-vivo macrophage killing assay. Moreover, the presence of the two minor alleles in the cis-eQTL-rs110345285 (C/C) was associated with the dysregulation of the eukaryotic elongation factor 1-α2 (eEF1A2) expression and with increased ELISA (OD) values. Finally, the presence of the minor allele in the cis-eQTL rs109859270 (C/T) was associated with the up-regulation of the U1 spliceosomal RNA expression and with an increased risk of progression to clinical PTB. The introduction of these novel functional variants into marker-assisted breeding programs is expected to have a relevant effect on PTB control.

DNA hypermethylation is associated with invasive phenotype of malignant melanoma.

Tumor cell invasion is one of the key processes during cancer progression, leading to life-threatening metastatic lesions in melanoma. As methylation of cancer-related genes plays a fundamental role during tumorigenesis and may lead to cellular plasticity which promotes invasion, our aim was to identify novel epigenetic markers on selected invasive melanoma cells. Using Illumina BeadChip assays and Affymetrix Human Gene 1.0 microarrays, we explored the DNA methylation landscape of selected invasive melanoma cells and examined the impact of DNA methylation on gene expression patterns. Our data revealed predominantly hypermethylated genes in the invasive cells affecting the neural crest differentiation pathway and regulation of the actin cytoskeleton. Integrative analysis of the methylation and gene expression profiles resulted in a cohort of hypermethylated genes (IL12RB2, LYPD6B, CHL1, SLC9A3, BAALC, FAM213A, SORCS1, GPR158, FBN1 and ADORA2B) with decreased expression. On the other hand, hypermethylation in the gene body of the EGFR and RBP4 genes was positively correlated with overexpression of the genes. We identified several methylation changes that can have role during melanoma progression, including hypermethylation of the promoter regions of the ARHGAP22 and NAV2 genes that are commonly altered in locally invasive primary melanomas as well as during metastasis. Interestingly, the down-regulation of the methylcytosine dioxygenase TET2 gene, which regulates DNA methylation, was associated with hypermethylated promoter region of the gene. This can probably lead to the observed global hypermethylation pattern of invasive cells and might be one of the key changes during the development of malignant melanoma cells.

Human mitochondrial DNA is extensively methylated in a non-CpG context.

Mitochondrial dysfunction plays critical roles in cancer development and related therapeutic response; however, exact molecular mechanisms remain unclear. Recently, alongside the discovery of mitochondrial-specific DNA methyltransferases, global and site-specific methylation of the mitochondrial genome has been described. Investigation of any functional consequences however remains unclear and debated due to insufficient evidence of the quantitative degree and frequency of mitochondrial DNA (mtDNA) methylation. This study uses WGBS to provide the first quantitative report of mtDNA methylation at single base pair resolution. The data show that mitochondrial genomes are extensively methylated predominantly at non-CpG sites. Importantly, these methylation patterns display notable differences between normal and cancer cells. Furthermore, knockdown of DNA methyltransferase enzymes resulted in a marked global reduction of mtDNA methylation levels, indicating these enzymes may be associated with the establishment and/or maintenance of mtDNA methylation. DNMT3B knockdown cells displayed a comparatively pronounced global reduction in mtDNA methylation with concomitant increases in gene expression, suggesting a potential functional link between methylation and gene expression. Together these results demonstrate reproducible, non-random methylation patterns of mtDNA and challenge the notion that mtDNA is lowly methylated. This study discusses key differences in methodology that suggest future investigations must allow for techniques that assess both CpG and non-CpG methylation.

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.

Genome-wide profiling of normal gastric mucosa identifies Helicobacter pylori- and cancer-associated DNA methylome changes.

The large geographic variations in the incidence of gastric cancer (GC) are likely due to differential environmental exposures, in particular to Helicobacter pylori (H. pylori) infection. We aimed to investigate the impact of H. pylori on the epigenome in normal gastric mucosa and methylation changes associated with cancer risk independent of H. pylori. A discovery set of normal gastric mucosa from GC cases (n = 42) and controls (n = 42), nested in a large case-control study and stratified by H. pylori status, were subjected to genome-wide methylation profiling. Single-nucleotide polymorphism arrays from peripheral blood leukocytes were used to conduct methylation quantitative trait loci (mQTL) analysis. A validation set of gastric mucosa samples (n = 180) was used in the replication phase. We found 1,924 differentially methylated positions (DMPs) and 438 differentially methylated regions (DMRs) associated with H. pylori infection, most of which were hypermethylated. Significant methylation alterations identified in the initial set were successfully replicated. Furthermore, the H. pylori-associated DMP/Rs showed marked stability ('epigenetic memory') after H. pylori clearance. Interestingly, we found 152 DMRs associated with cancer risk independent of the H. pylori status in normal gastric mucosa. The methylation score derived from three biomarkers was a strong predictor of GC. Finally, the mQTL analysis indicated that the H. pylori- and cancer-specific methylation signatures were minimally affected by genetic variation. The comprehensively characterized methylome changes associated with H. pylori infection and GC risk in our study might serve as potential biomarkers for early cancer progression in tumour-free gastric mucosa.

Lowly methylated region analysis identifies EBF1 as a potential epigenetic modifier in breast cancer.

Breast cancer (BC) encompasses heterogeneous pathologies with different subtypes exhibiting distinct molecular changes, including those related to DNA methylation. However, the role of these changes in mediating BC heterogeneity is poorly understood. Lowly methylated regions (LMRs), non-CpG island loci that usually contain transcription factor (TF) binding sites, have been suggested to act as regulatory elements that define cellular identity. In this study, we aimed to identify the key subtype-specific TFs that may lead to LMR generation and shape the BC methylome and transcription program. We initially used whole-genome bisulfite sequencing (WGBS) data available at The Cancer Genome Atlas (TCGA) portal to identify subtype-specific LMRs. Differentially methylated regions (DMRs) within the BC PAM50 subtype-specific LMRs were selected by comparing tumors and normal tissues in a larger TCGA cohort assessed by HumanMethylation450 BeadChip (450K) arrays and TF enrichment analyses were performed. To assess the impact of LMRs on gene expression, TCGA RNA sequencing data were downloaded and Pearson correlations between methylation levels of loci presenting subtype-specific TF motifs and expression of the nearest genes were calculated. WGBS methylome data revealed a large number of LMRs for each of the BC subtypes. Analysis of these LMRs in the 450K datasets available for a larger sample set identified 7,765, 5,657, and 19 differentially methylated positions (DMPs) between normal adjacent tissues and tumor tissues from basal, luminal, and HER2-enriched subtypes, respectively. Unsupervised clustering showed that the discriminatory power of the top DMPs was remarkably strong for basal BC. Interestingly, in this particular subtype, we found 4,409 differentially hypomethylated positions grouped into 1,185 DMRs with a strong enrichment for the early B-cell factor 1 (EBF1) motifs. The methylation levels of the DMRs containing EBF1 motifs showed a strong negative correlation with the expression of 719 nearby genes, including BTS2 and CD74, two oncogenes known to be specific for basal BC subtype and for poor outcome. This study identifies LMRs specific to the three main BC subtypes and reveals EBF1 as a potentially important regulator of BC subtype-specific methylation and gene expression program.

Unique DNA methylation signature in HPV-positive head and neck squamous cell carcinomas.

BACKGROUND: Head and neck squamous cell carcinomas (HNSCCs) represent a heterogeneous group of cancers for which human papilloma virus (HPV) infection is an emerging risk factor. Previous studies showed promoter hypermethylation in HPV(+) oropharyngeal cancers, but only few consistent target genes have been so far described, and the evidence of a functional impact on gene expression is still limited. METHODS: We performed global and stratified pooled analyses of epigenome-wide data in HNSCCs based on the Illumina HumanMethylation450 bead-array data in order to identify tissue-specific components and common viral epigenetic targets in HPV-associated tumours. RESULTS: We identified novel differentially methylated CpGs and regions associated with viral infection that are independent of the anatomic site. In particular, most hypomethylated regions were characterized by a marked loss of CpG island boundaries, which showed significant correlations with expression of neighbouring genes. Moreover, a subset of only five CpGs in a few hypomethylated regions predicted HPV status with a high level of specificity in different cohorts. Finally, this signature was a better predictor of survival compared with HPV status determined by viral gene expression by RNA sequencing in The Cancer Genome Atlas cohort. CONCLUSIONS: We identified a novel epigenetic signature of HPV infection in HNSCCs which is independent of the anatomic site, is functionally correlated with gene expression and may be leveraged for improved stratification of prognosis in HNSCCs.

A long noncoding RNA associated with susceptibility to celiac disease.

Recent studies have implicated long noncoding RNAs (lncRNAs) as regulators of many important biological processes. Here we report on the identification and characterization of a lncRNA, lnc13, that harbors a celiac disease-associated haplotype block and represses expression of certain inflammatory genes under homeostatic conditions. Lnc13 regulates gene expression by binding to hnRNPD, a member of a family of ubiquitously expressed heterogeneous nuclear ribonucleoproteins (hnRNPs). Upon stimulation, lnc13 levels are reduced, thereby allowing increased expression of the repressed genes. Lnc13 levels are significantly decreased in small intestinal biopsy samples from patients with celiac disease, which suggests that down-regulation of lnc13 may contribute to the inflammation seen in this disease. Furthermore, the lnc13 disease-associated variant binds hnRNPD less efficiently than its wild-type counterpart, thus helping to explain how these single-nucleotide polymorphisms contribute to celiac disease.

Identification of novel long non-coding RNAs deregulated in hepatocellular carcinoma using RNA-sequencing.

Functional characterization of long non-coding RNAs (lncRNAs) and their pathological relevance is still a challenging task. Abnormal expression of a few long non-coding RNAs have been found associated with hepatocellular carcinoma, with potential implications to both improve our understanding of molecular mechanism of liver carcinogenesis and to discover biomarkers for early diagnosis or therapy. However, the understanding of the global role of lncRNAs during HCC development is still in its infancy. In this study, we produced RNA-Seq data from 23 liver tissues (controls, cirrhotic and HCCs) and applied statistical and gene network analysis approaches to identify and characterize expressed lncRNAs. We detected 5,525 lncRNAs across different tissue types and identified 57 differentially expressed lncRNAs in HCC compared with adjacent non-tumour tissues using stringent criteria (FDR<0.05, Fold Change>2). Using weighted gene co-expression network analysis (WGCNA), we found that differentially expressed lncRNAs are co-expressed with genes involved in cell cycle regulation, TGF-ß signalling and liver metabolism. Furthermore, we found that more than 20% of differentially expressed lncRNAs are associated to actively transcribed enhancers and that the co-expression patterns with their closest genes change dramatically during HCC development. Our study provides the most comprehensive compendium of lncRNAs expressed in HCC, as well as in control or cirrhotic livers. Our results identified both known oncogenic lncRNAs (such as H19 and CRNDE) and novel lncRNAs involved in cell cycle deregulation and liver metabolism deficits occurring during HCC development.

Patterns of somatic uniparental disomy identify novel tumor suppressor genes in colorectal cancer.

Colorectal cancer (CRC) is characterized by specific patterns of copy number alterations (CNAs), which helped with the identification of driver oncogenes and tumor suppressor genes (TSGs). More recently, the usage of single nucleotide polymorphism arrays provided information of copy number neutral loss of heterozygosity, thus suggesting the occurrence of somatic uniparental disomy (UPD) and uniparental polysomy (UPP) events. The aim of this study is to establish an integrative profiling of recurrent UPDs/UPPs and CNAs in sporadic CRC. Our results indicate that regions showing high frequencies of UPD/UPP mostly coincide with regions typically involved in genomic losses. Among them, chromosome arms 3p, 5q, 9q, 10q, 14q, 17p, 17q, 20p, 21q and 22q preferentially showed UPDs/UPPs over genomic losses suggesting that tumor cells must maintain the disomic state of certain genes to favor cellular fitness. A meta-analysis using over 300 samples from The Cancer Genome Atlas confirmed our findings. Several regions affected by recurrent UPDs/UPPs contain well-known TSGs, as well as novel candidates such as ARID1A, DLC1, TCF7L2 and DMBT1. In addition, VCAN, FLT4, SFRP1 and GAS7 were also frequently involved in regions of UPD/UPP and displayed high levels of methylation. Finally, sequencing and fluorescence in situ hybridization analysis of the gene APC underlined that a somatic UPD event might represent the second hit to achieve biallelic inactivation of this TSG in colorectal tumors. In summary, our data define a profile of somatic UPDs/UPPs in sporadic CRC and highlights the importance of these events as a mechanism to achieve the inactivation of TSGs.

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.

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.

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.

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.