Nora Fernandez-Jimenez

A long noncoding RNA associated with susceptibility to celiac disease.

Long noncoding RNAs in inflammation Growing evidence suggests that long noncoding RNAs (lncRNAs) are important modulators of gene expression. Castellanos-Rubio et al. identified a lncRNA, lnc13, that suppresses inflammatory gene expression in macrophages (see the Perspective by Huarte). Lnc13 interacts with proteins that regulate chromatin accessibility. Stimulating macrophages with a cell wall component from bacteria decreased expression of lnc13 and increased the expression of several inflammatory genes. Decreased levels of lnc13 in intestinal tissue from individuals with celiac disease hint that lnc13 may also play a role in the pathogenesis of immune-mediated diseases. Science, this issue p. 91; see also p. 43 A specific RNA molecule can suppress inflammatory gene expression and may protect from gluten intolerance. [Also see Perspective by Huarte] 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.

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.

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.

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.

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.

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.