Sex bias in celiac disease: XWAS and monocyte eQTLs in women identify TMEM187 as a functional candidate gene.

BACKGROUND: Celiac disease (CeD) is an immune-mediated disorder that develops in genetically predisposed individuals upon gluten consumption. HLA risk alleles explain 40% of the genetic component of CeD, so there have been continuing efforts to uncover non-HLA loci that can explain the remaining heritability. As in most autoimmune disorders, the prevalence of CeD is significantly higher in women. Here, we investigated the possible involvement of the X chromosome on the sex bias of CeD. METHODS: We performed a X chromosome-wide association study (XWAS) and a gene-based association study in women from the CeD Immunochip (7062 cases, 5446 controls). We also constructed a database of X chromosome cis-expression quantitative trait loci (eQTLs) in monocytes from unstimulated (n = 226) and lipopolysaccharide (LPS)-stimulated (n = 130) female donors and performed a Summary-data-based MR (SMR) analysis to integrate XWAS and eQTL information. We interrogated the expression of the potentially causal gene (TMEM187) in peripheral blood mononuclear cells (PBMCs) from celiac patients at onset, on a gluten-free diet, potential celiac patients and non-celiac controls. RESULTS: The XWAS and gene-based analyses identified 13 SNPs and 25 genes, respectively, 22 of which had not been previously associated with CeD. The X chromosome cis-eQTL analysis found 18 genes with at least one cis-eQTL in naïve female monocytes and 8 genes in LPS-stimulated female monocytes, 2 of which were common to both situations and 6 were unique to LPS stimulation. SMR identified a potentially causal association of TMEM187 expression in naïve monocytes with CeD in women, regulated by CeD-associated, eQTL-SNPs rs7350355 and rs5945386. The CeD-risk alleles were correlated with lower TMEM187 expression. These results were replicated using eQTLs from LPS-stimulated monocytes. We observed higher levels of TMEM187 expression in PBMCs from female CeD patients at onset compared to female non-celiac controls, but not in male CeD individuals. CONCLUSION: Using X chromosome genotypes and gene expression data from female monocytes, SMR has identified TMEM187 as a potentially causal candidate in CeD. Further studies are needed to understand the implication of the X chromosome in the higher prevalence of CeD in women.

Celiac disease (CeD) is an immune-related condition triggered by gluten consumption in genetically susceptible individuals. Women present higher prevalence of CeD than men, but the biological explanation of such difference has not been elucidated. In this study, we investigated whether specific genetic variations on the X chromosome were associated with CeD in each sex. Surprisingly, we found 13 genetic variants and 25 genes significantly linked to CeD in women, but not in men. Additionally, we identified genetic variants on the X chromosome associated with gene expression of monocytes, a type of immune cells that is activated in CeD after gluten intake. Integrating these data with our previous findings, we found that lower expression of a gene termed TMEM187 might be associated with a potential increase in CeD risk in women. Finally, validation experiments confirmed higher TMEM187 levels in blood cells from female CeD patients compared to non-celiac women, while no such difference was seen in males. In summary, our study suggests that the X-chromosome gene TMEM187 may play a key role in CeD development, providing insights into the higher prevalence of CeD in females.

Genome-wide association study of placental weight identifies distinct and shared genetic influences between placental and fetal growth.

A well-functioning placenta is essential for fetal and maternal health throughout pregnancy. Using placental weight as a proxy for placental growth, we report genome-wide association analyses in the fetal (n = 65,405), maternal (n = 61,228) and paternal (n = 52,392) genomes, yielding 40 independent association signals. Twenty-six signals are classified as fetal, four maternal and three fetal and maternal. A maternal parent-of-origin effect is seen near KCNQ1. Genetic correlation and colocalization analyses reveal overlap with birth weight genetics, but 12 loci are classified as predominantly or only affecting placental weight, with connections to placental development and morphology, and transport of antibodies and amino acids. Mendelian randomization analyses indicate that fetal genetically mediated higher placental weight is causally associated with preeclampsia risk and shorter gestational duration. Moreover, these analyses support the role of fetal insulin in regulating placental weight, providing a key link between fetal and placental growth.

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.

Two-Sample Mendelian Randomization detects bidirectional causality between gut microbiota and celiac disease in individuals with high genetic risk.

Background: Celiac Disease (CeD) is an autoimmune disorder triggered by gluten intake in genetically susceptible individuals. Highest risk individuals are homozygous for the Human Leucocyte Antigen (HLA) DQ2.5 haplotype or DQ2.5/DQ2.2 heterozygous. Both the HLA-DQ2-positive high genetic risk individuals and those that have developed the disease have altered intestinal microbiota, but it remains unclear whether these alterations are a cause or a consequence of CeD. Objective: To investigate a potential bidirectional causality between gut microbiota (GM) and CeD in HLA-DQ2 high genetic risk individuals. Materials and Methods: We performed a bidirectional Two-Sample Mendelian Randomization (2SMR) test using summary statistics from the largest publicly available Genome-Wide Association Study (GWAS) of GM and the summary statistics of the Immunochip CeD study of those individuals with the HLA-DQ2 high-risk haplotype. To test whether changes in GM composition were causally linked to CeD, GM data were used as exposure and CeD data as outcome; to test for reverse causation, the exposure and outcome datasets were inverted. Results: We identified several bacteria from Ruminococcaceae and Lachnospiraceae families of the Firmicutes phylum as potentially causal in both directions. In addition, our results suggest that changes in the abundance of Veillonellaceae family might be causal in the development of CeD, while alterations in Pasteurellaceae family might be a consequence of the disease itself. Conclusion: Our results suggest that the relationship between GM and HLA-DQ2 high risk individuals is highly complex and bidirectional.

Potentially causal associations between placental DNA methylation and schizophrenia and other neuropsychiatric disorders.

Increasing evidence supports the role of placenta in neurodevelopment and potentially, in the later onset of neuropsychiatric disorders. Recently, methylation quantitative trait loci (mQTL) and interaction QTL (iQTL) maps have proven useful to understand SNP-genome wide association study (GWAS) relationships, otherwise missed by conventional expression QTLs. In this context, we propose that part of the genetic predisposition to complex neuropsychiatric disorders acts through placental DNA methylation (DNAm). We constructed the first public placental cis-mQTL database including nearly eight million mQTLs calculated in 368 fetal placenta DNA samples from the INMA project, ran cell type- and gestational age-imQTL models and combined those data with the summary statistics of the largest GWAS on 10 neuropsychiatric disorders using Summary-based Mendelian Randomization (SMR) and colocalization. Finally, we evaluated the influence of the DNAm sites identified on placental gene expression in the RICHS cohort. We found that placental cis-mQTLs are highly enriched in placenta-specific active chromatin regions, and useful to map the etiology of neuropsychiatric disorders at prenatal stages. Specifically, part of the genetic burden for schizophrenia, bipolar disorder and major depressive disorder confers risk through placental DNAm. The potential causality of several of the observed associations is reinforced by secondary association signals identified in conditional analyses, regional pleiotropic methylation signals associated to the same disorder, and cell type-imQTLs, additionally associated to the expression levels of relevant immune genes in placenta. In conclusion, the genetic risk of several neuropsychiatric disorders could operate, at least in part, through DNAm and associated gene expression in placenta.

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.

Gluten-induced RNA methylation changes regulate intestinal inflammation via allele-specific XPO1 translation in epithelial cells.

OBJECTIVES: Coeliac disease (CD) is a complex autoimmune disorder that develops in genetically susceptible individuals. Dietary gluten triggers an immune response for which the only available treatment so far is a strict, lifelong gluten free diet. Human leucocyte antigen (HLA) genes and several non-HLA regions have been associated with the genetic susceptibility to CD, but their role in the pathogenesis of the disease is still essentially unknown, making it complicated to develop much needed non-dietary treatments. Here, we describe the functional involvement of a CD-associated single-nucleotide polymorphism (SNP) located in the 5'UTR of XPO1 in the inflammatory environment characteristic of the coeliac intestinal epithelium. DESIGN: The function of the CD-associated SNP was investigated using an intestinal cell line heterozygous for the SNP, N6-methyladenosine (m6A)-related knock-out and HLA-DQ2 mice, and human samples from patients with CD. RESULTS: Individuals harbouring the risk allele had higher m6A methylation in the 5'UTR of XPO1 RNA, rendering greater XPO1 protein amounts that led to downstream nuclear factor kappa B (NFkB) activity and subsequent inflammation. Furthermore, gluten exposure increased overall m6A methylation in humans as well as in in vitro and in vivo models. CONCLUSION: We identify a novel m6A-XPO1-NFkB pathway that is activated in CD patients. The findings will prompt the development of new therapeutic approaches directed at m6A proteins and XPO1, a target under evaluation for the treatment of intestinal disorders.

shinyCurves, a shiny web application to analyse multisource qPCR amplification data: a COVID-19 case study.

BACKGROUND: Quantitative, reverse transcription PCR (qRT-PCR) is currently the gold-standard for SARS-CoV-2 detection and it is also used for detection of other virus. Manual data analysis of a small number of qRT-PCR plates per day is a relatively simple task, but automated, integrative strategies are needed if a laboratory is dealing with hundreds of plates per day, as is being the case in the COVID-19 pandemic. RESULTS: Here we present shinyCurves, an online shiny-based, free software to analyze qRT-PCR amplification data from multi-plate and multi-platform formats. Our shiny application does not require any programming experience and is able to call samples Positive, Negative or Undetermined for viral infection according to a number of user-defined settings, apart from providing a complete set of melting and amplification curve plots for the visual inspection of results. CONCLUSIONS: shinyCurves is a flexible, integrative and user-friendly software that speeds-up the analysis of massive qRT-PCR data from different sources, with the possibility of automatically producing and evaluating melting and amplification curve plots.

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.

Celiac disease susceptibility: The genome and beyond.

Celiac Disease (CeD) is an immune-mediated complex disease that is triggered by the ingestion of gluten and develops in genetically susceptible individuals. It has been known for a long time that the Human Leucocyte Antigen (HLA) molecules DQ2 and DQ8 are necessary, although not sufficient, for the disease development, and therefore other susceptibility genes and (epi)genetic events must participate in CeD pathogenesis. The advances in Genomics during the last 15 years have made CeD one of the immune-related disorders with the best-characterized genetic component. In the present work, we will first review the main Genome-Wide Association Studies (GWAS) carried out in the disorder, and emphasize post-GWAS discoveries, including diverse integrative strategies, SNP prioritization approaches, and insights into the Microbiome through the host Genomics. Second, we will explore CeD-related Epigenetics and Epigenomics, mostly focusing on the emerging knowledge of the celiac methylome, and the vast but yet under-explored non-coding RNA (ncRNA) landscape. We conclude that much has been done in the field although there are still completely unvisited areas in the post-Genomics of CeD. Chromatin conformation and accessibility, and Epitranscriptomics are promising domains that need to be unveiled to complete the big picture of the celiac Genome.

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.

A Systematic Two-Sample Mendelian Randomization Analysis Identifies Shared Genetic Origin of Endometriosis and Associated Phenotypes.

Endometriosis, one of the most common gynecological disorders, is a complex disease characterized by the growth of endometrial-like tissue in extra-uterine locations and is a cause of pelvic pain and infertility. Evidence from observational studies indicate that endometriosis usually appears together with several other phenotypes. These include a list of autoimmune diseases, most of them more prevalent in women, anthropometric traits associated with leanness in the adulthood, as well as female reproductive traits, including altered hormone levels and those associated with a prolonged exposure to menstruation. However, the biological mechanisms underlying their co-morbidity remains unknown. To explore whether those phenotypes and endometriosis share a common genetic origin, we performed a systematic Two-Sample Mendelian Randomization (2SMR) analysis using public GWAS data. Our results suggest potential common genetic roots between endometriosis and female anthropometric and reproductive traits. Particularly, our data suggests that reduced weight and BMI might be mediating the genetic susceptibility to suffer endometriosis. Furthermore, data on female reproductive traits strongly suggest that genetic variants that predispose to a more frequent exposure to menstruation, through earlier age at menarche and shorter menstrual cycles, might also increase the risk to suffer from endometriosis.

Novel loci for childhood body mass index and shared heritability with adult cardiometabolic traits.

The genetic background of childhood body mass index (BMI), and the extent to which the well-known associations of childhood BMI with adult diseases are explained by shared genetic factors, are largely unknown. We performed a genome-wide association study meta-analysis of BMI in 61,111 children aged between 2 and 10 years. Twenty-five independent loci reached genome-wide significance in the combined discovery and replication analyses. Two of these, located near NEDD4L and SLC45A3, have not previously been reported in relation to either childhood or adult BMI. Positive genetic correlations of childhood BMI with birth weight and adult BMI, waist-to-hip ratio, diastolic blood pressure and type 2 diabetes were detected (Rg ranging from 0.11 to 0.76, P-values <0.002). A negative genetic correlation of childhood BMI with age at menarche was observed. Our results suggest that the biological processes underlying childhood BMI largely, but not completely, overlap with those underlying adult BMI. The well-known observational associations of BMI in childhood with cardio-metabolic diseases in adulthood may reflect partial genetic overlap, but in light of previous evidence, it is also likely that they are explained through phenotypic continuity of BMI from childhood into adulthood.

A Two-Sample Mendelian Randomization Analysis Investigates Associations Between Gut Microbiota and Celiac Disease.

Celiac disease (CeD) is a complex immune-mediated inflammatory condition triggered by the ingestion of gluten in genetically predisposed individuals. Literature suggests that alterations in gut microbiota composition and function precede the onset of CeD. Considering that microbiota is partly determined by host genetics, we speculated that the genetic makeup of CeD patients could elicit disease development through alterations in the intestinal microbiota. To evaluate potential causal relationships between gut microbiota and CeD, we performed a two-sample Mendelian randomization analysis (2SMR). Exposure data were obtained from the raw results of a previous genome-wide association study (GWAS) of gut microbiota and outcome data from summary statistics of CeD GWAS and Immunochip studies. We identified a number of putative associations between gut microbiota single nucleotide polymorphisms (SNPs) associated with CeD. Regarding bacterial composition, most of the associated SNPs were related to Firmicutes phylum, whose relative abundance has been previously reported to be altered in CeD patients. In terms of functional units, we linked a number of SNPs to several bacterial metabolic pathways that seemed to be related to CeD. Overall, this study represented the first 2SMR approach to elucidate the relationship between microbiome and CeD.

Highly Reproducible Hyperthermia Response in Water, Agar, and Cellular Environment by Discretely PEGylated Magnetite Nanoparticles.

Local heat generation from magnetic nanoparticles (MNPs) exposed to alternating magnetic fields can revolutionize cancer treatment. However, the application of MNPs as anticancer agents is limited by serious drawbacks. Foremost among these are the fast uptake and biodegradation of MNPs by cells and the unpredictable magnetic behavior of the MNPs when they accumulate within or around cells and tissues. In fact, several studies have reported that the heating power of MNPs is severely reduced in the cellular environment, probably due to a combination of increased viscosity and strong NP agglomeration. Herein, we present an optimized protocol to coat magnetite (Fe3O4) NPs larger than 20 nm (FM-NPs) with high molecular weight PEG molecules that avoid collective coatings, prevent the formation of large clusters of NPs and keep constant their high heating performance in environments with very different ionic strengths and viscosities (distilled water, physiological solutions, agar and cell culture media). The great reproducibility and reliability of the heating capacity of this FM-NP@PEG system in such different environments has been confirmed by AC magnetometry and by more conventional calorimetric measurements. The explanation of this behavior has been shown to lie in preserving as much as possible the magnetic single domain-type behavior of nearly isolated NPs. In vitro endocytosis experiments in a colon cancer-derived cell line indicate that FM-NP@PEG formulations with PEGs of higher molecular weight (20 kDa) are more resistant to endocytosis than formulations with smaller PEGs (5 kDa), showing quite large uptake mean-life (τ > 5 h) in comparison with other NP systems. The in vitro magnetic hyperthermia was performed at 21 mT and 650 kHz during 1 h in a pre-endocytosis stage and complete cell death was achieved 48 h posthyperthermia. These optimal FM-NP@PEG formulations with high resistance to endocytosis and predictable magnetic response will aid the progress and accuracy of the emerging era of theranostics.

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.

Comparison of Illumina 450K and EPIC arrays in placental DNA methylation.

Illumina HumanMethylation450 BeadChip (450K) has been commonly used to investigate DNA methylation in human tissues. Recently, it has been replaced by Illumina HumanMethylationEPIC BeadChip (EPIC) covering over 850,000 CpGs distributed genome-wide. Many consortia have now datasets coming from both arrays and aspire to analyze the two together. The placenta shows a high number of intermediate methylation levels and is often investigated for obstetric/birth outcomes, and potentially for long-term programming in offspring. We performed a systematic comparison between the two arrays using 108 duplicate placental samples from Gen3G birth cohort. We find that placenta shows a high per-sample correlation between the arrays, and higher median correlations at individual CpGs than those reported for blood. We identify 26,340 probes with absolute difference in per cent methylation >10%. We conclude that EPIC and 450K placental data can be combined, and we provide two lists of CpGs that should be excluded to avoid misleading results.

Mendelian randomization analysis of celiac GWAS reveals a blood expression signature with diagnostic potential in absence of gluten consumption.

Celiac disease (CeD) is an immune-mediated enteropathy with a strong genetic component where the main environmental trigger is dietary gluten, and currently a correct diagnosis of the disease is impossible if gluten-free diet (GFD) has already been started. We hypothesized that merging different levels of genomic information through Mendelian randomization (MR) could help discover genetic biomarkers useful for CeD diagnosis. MR was performed using public databases of expression quantitative trait loci (QTL) and methylation QTL as exposures and the largest CeD genome-wide association study conducted to date as the outcome, in order to identify potential causal genes. As a result, we identified UBE2L3, an ubiquitin ligase located in a CeD-associated region. We interrogated the expression of UBE2L3 in an independent data set of peripheral blood mononuclear cells (PBMCs) and found that its expression is altered in CeD patients on GFD when compared to non-celiac controls. The relative expression of UBE2L3 isoforms predicts CeD with 100% specificity and sensitivity and could be used as a diagnostic marker, especially in the absence of gluten consumption. This approach could be applicable to other diseases where diagnosis of asymptomatic patients can be complicated.