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.

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.

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.

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.

MethylCal: Bayesian calibration of methylation levels.

Bisulfite amplicon sequencing has become the primary choice for single-base methylation quantification of multiple targets in parallel. The main limitation of this technology is a preferential amplification of an allele and strand in the PCR due to methylation state. This effect, known as 'PCR bias', causes inaccurate estimation of the methylation levels and calibration methods based on standard controls have been proposed to correct for it. Here, we present a Bayesian calibration tool, MethylCal, which can analyse jointly all CpGs within a CpG island (CGI) or a Differentially Methylated Region (DMR), avoiding 'one-at-a-time' CpG calibration. This enables more precise modeling of the methylation levels observed in the standard controls. It also provides accurate predictions of the methylation levels not considered in the controlled experiment, a feature that is paramount in the derivation of the corrected methylation degree. We tested the proposed method on eight independent assays (two CpG islands and six imprinting DMRs) and demonstrated its benefits, including the ability to detect outliers. We also evaluated MethylCal's calibration in two practical cases, a clinical diagnostic test on 18 patients potentially affected by Beckwith-Wiedemann syndrome, and 17 individuals with celiac disease. The calibration of the methylation levels obtained by MethylCal allows a clearer identification of patients undergoing loss or gain of methylation in borderline cases and could influence further clinical or treatment decisions.

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.

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 Diasease-associated lncRNA Named HCG14 Regulates NOD1 Expression in Intestinal Cells.

OBJECTIVE: The aim of the study is to identify additional celiac disease associated loci in the major histocompatibility complex (MHC) independent from classical HLA risk alleles (HLA-DR3-DQ2) and to characterize their potential functional impact in celiac disease pathogenesis at the intestinal level. METHODS: We performed a high-resolution single-nucleotide polymorphism (SNP) genotyping of the MHC region, comparing HLA-DR3 homozygous celiac patients and non-celiac controls carrying a single copy of the B8-DR3-DQ2 conserved extended haplotype. Expression level of potential novel risk genes was determined by RT-PCR in intestinal biopsies and in intestinal and immune cells isolated from control and celiac individuals. Small interfering RNA-driven silencing of selected genes was performed in the intestinal cell line T84. RESULTS: MHC genotyping revealed 2 associated SNPs, one located in TRIM27 gene and another in the non-coding gene HCG14. After stratification analysis, only HCG14 showed significant association independent from HLA-DR-DQ loci. Expression of HCG14 was slightly downregulated in epithelial cells isolated from duodenal biopsies of celiac patients, and eQTL analysis revealed that polymorphisms in HCG14 region were associated with decreased NOD1 expression in duodenal intestinal cells. CONCLUSIONS: We have successfully employed a conserved extended haplotype-matching strategy and identified a novel additional celiac disease risk variant in the lncRNA HCG14. This lncRNA seems to regulate the expression of NOD1 in an allele-specific manner. Further functional studies are needed to clarify the role of HCG14 in the regulation of gene expression and to determine the molecular mechanisms by which the risk variant in HCG14 contributes to celiac disease pathogenesis.

Recommendations to report and interpret HLA genetic findings in coeliac disease.

Coeliac disease (CD) is a chronic autoimmune enteropathy triggered by gluten and related prolamines in genetically predisposed individuals. Although CD is a polygenic disease, there is a strong association with genes of the human leukocyte antigen (HLA) region. Most patients present the HLA-DQ2 heterodimer, specifically the DQ2.5 isoform, which is present in around 90-96% of patients of European ancestry.

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

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

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.

Fine mapping of the celiac disease-associated LPP locus reveals a potential functional variant.

Using the Immunochip for genotyping, we identified 39 non-human leukocyte antigen (non-HLA) loci associated to celiac disease (CeD), an immune-mediated disease with a worldwide frequency of ∼1%. The most significant non-HLA signal mapped to the intronic region of 70 kb in the LPP gene. Our aim was to fine map and identify possible functional variants in the LPP locus. We performed a meta-analysis in a cohort of 25 169 individuals from six different populations previously genotyped using Immunochip. Imputation using data from the Genome of the Netherlands and 1000 Genomes projects, followed by meta-analysis, confirmed the strong association signal on the LPP locus (rs2030519, P = 1.79 × 10(-49)), without any novel associations. The conditional analysis on this top SNP-indicated association to a single common haplotype. By performing haplotype analyses in each population separately, as well as in a combined group of the four populations that reach the significant threshold after correction (P < 0.008), we narrowed down the CeD-associated region from 70 to 2.8 kb (P = 1.35 × 10(-44)). By intersecting regulatory data from the ENCODE project, we found a functional SNP, rs4686484 (P = 3.12 × 10(-49)), that maps to several B-cell enhancer elements and a highly conserved region. This SNP was also predicted to change the binding motif of the transcription factors IRF4, IRF11, Nkx2.7 and Nkx2.9, suggesting its role in transcriptional regulation. We later found significantly low levels of LPP mRNA in CeD biopsies compared with controls, thus our results suggest that rs4686484 is the functional variant in this locus, while LPP expression is decreased in CeD.

Improving coeliac disease risk prediction by testing non-HLA variants additional to HLA variants.

BACKGROUND: The majority of coeliac disease (CD) patients are not being properly diagnosed and therefore remain untreated, leading to a greater risk of developing CD-associated complications. The major genetic risk heterodimer, HLA-DQ2 and DQ8, is already used clinically to help exclude disease. However, approximately 40% of the population carry these alleles and the majority never develop CD. OBJECTIVE: We explored whether CD risk prediction can be improved by adding non-HLA-susceptible variants to common HLA testing. DESIGN: We developed an average weighted genetic risk score with 10, 26 and 57 single nucleotide polymorphisms (SNP) in 2675 cases and 2815 controls and assessed the improvement in risk prediction provided by the non-HLA SNP. Moreover, we assessed the transferability of the genetic risk model with 26 non-HLA variants to a nested case-control population (n=1709) and a prospective cohort (n=1245) and then tested how well this model predicted CD outcome for 985 independent individuals. RESULTS: Adding 57 non-HLA variants to HLA testing showed a statistically significant improvement compared to scores from models based on HLA only, HLA plus 10 SNP and HLA plus 26 SNP. With 57 non-HLA variants, the area under the receiver operator characteristic curve reached 0.854 compared to 0.823 for HLA only, and 11.1% of individuals were reclassified to a more accurate risk group. We show that the risk model with HLA plus 26 SNP is useful in independent populations. CONCLUSIONS: Predicting risk with 57 additional non-HLA variants improved the identification of potential CD patients. This demonstrates a possible role for combined HLA and non-HLA genetic testing in diagnostic work for CD.

Celiac disease in the Mediterranean area.

BACKGROUND: The World Gastroenterology Organization recommends developing national guidelines for the diagnosis of Celiac Disease (CD): hence a profile of the diagnosis of CD in each country is required. We aim to describe a cross-sectional picture of the clinical features and diagnostic facilities in 16 countries of the Mediterranean basin. Since a new ESPGHAN diagnostic protocol was recently published, our secondary aim is to estimate how many cases in the same area could be identified without a small intestinal biopsy. METHODS: By a stratified cross-sectional retrospective study design, we examined clinical, histological and laboratory data from 749 consecutive unselected CD children diagnosed by national referral centers. RESULTS: The vast majority of cases were diagnosed before the age of 10 (median: 5 years), affected by diarrhea, weight loss and food refusal, as expected. Only 59 cases (7.8%) did not suffer of major complaints. Tissue transglutaminase (tTG) assay was available, but one-third of centers reported financial constraints in the regular purchase of the assay kits. 252 cases (33.6%) showed tTG values over 10 times the local normal limit. Endomysial antibodies and HLA typing were routinely available in only half of the centers. CD was mainly diagnosed from small intestinal biopsy, available in all centers. Based on these data, only 154/749 cases (20.5%) would have qualified for a diagnosis of CD without a small intestinal biopsy, according to the new ESPGHAN protocol. CONCLUSIONS: This cross-sectional study of CD in the Mediterranean referral centers offers a puzzling picture of the capacities to deal with the emerging epidemic of CD in the area, giving a substantive support to the World Gastroenterology Organization guidelines.

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.

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.

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.

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.