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.

Genetic Susceptibility to Neurotoxicity Related to Prenatal Inorganic Arsenic Exposure in Young Spanish Children.

We explored the influence of child and maternal single nucleotide polymorphisms (SNPs) in genes related to neurological function and arsenic metabolism (i.e., ABCA1, ABCB1, PON1, CYP3A, BDNF, GSTP1, MT2A, and APOE as well as AS3MT) on the association between prenatal arsenic (As) exposure and methylation efficiency and neuropsychological development in 4-5-year-old children. Participants were 549 mother-child pairs from the INMA (Environment and Childhood) Spanish Project. We measured inorganic arsenic (iAs) and the metabolites monomethylarsonic acid (MMA) and dimethylarsinic acid (DMA) in urine samples collected during pregnancy. Neuropsychological development was assessed at the age of 4-5 years using the McCarthy Scales of Children's Abilities (MSCA). Several SNPs were determined in maternal and child DNA; AS3MT and APOE haplotypes were inferred. The median ∑As (sum of iAs, DMA, and MMA) was 7.08 µg/g creatinine. Statistically significant interactions for children's APOE haplotype were observed. Specifically, ε4-carrier children had consistently lower MSCA scores in several scales with increasing ∑As and MMA concentrations. These results provide evidence regarding the neurotoxic effects of early life exposure to As, observing that the APOE ε4 allele could make children more vulnerable to this exposure.

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.

Maternal Perfluoroalkyl Substances, Thyroid Hormones, and DIO Genes: A Spanish Cross-sectional Study.

Results of studies on perfluoroalkyl substances (PFASs) and thyroid hormones (THs) are heterogeneous, and the mechanisms underlying the action of PFASs to target THs have not been fully characterized. We examined the relation between first-trimester maternal PFAS and TH levels and the role played by polymorphisms in the iodothyronine deiodinase 1 (DIO1) and 2 (DIO2) genes in this association. Our sample comprised 919 pregnant Spanish women (recruitment = 2003-2008) with measurements of perfluorohexanesulfonic acid (PFHxS), perfluorooctanoic acid (PFOA), perfluorooctanesulfonic acid (PFOS), perfluorononanoic acid (PFNA), thyroid-stimulating hormone (TSH), total triiodothyronine (TT3), and free thyroxine (FT4), and we genotyped for single-nucleotide polymorphisms in the DIO1 (rs2235544) and DIO2 (rs12885300) genes. We performed multivariate regression analyses between PFASs and THs and included the interaction term PFAS-genotypes in the models. PFHxS was associated with an increase in TSH (% change in outcome [95% CI] per 2-fold PFAS increase = 6.09 [-0.71, 13.4]), and PFOA and PFNA were associated with a decrease in TT3 (-7.17 [-13.5, -0.39] and -6.28 [-12.3, 0.12], respectively). We found stronger associations between PFOA, PFNA, and TT3 for DIO1-CC and DIO2-CT genotypes, although interaction p-values were not significant. In conclusion, this study found evidence of an inverse association between PFOA and TT3 levels. No clear effect modification by DIO enzyme genes was observed.

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.

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.

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.

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.

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.

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.

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.

Influence of genetic polymorphisms on arsenic methylation efficiency during pregnancy: Evidence from a Spanish birth cohort.

BACKGROUND: Inorganic arsenic (iAs) is a widespread toxic metalloid. It is well-known that iAs metabolism and its toxicity are mediated by polymorphisms in AS3MT and other genes. However, studies during pregnancy are scarce. We aimed to examine the role of genetic polymorphisms in AS3MT, GSTO2, N6AMT1, MTHFR, MTR, FTCD, CBS, and FOLH1 in iAs methylation efficiency during pregnancy. METHODS: The study included 541 pregnant participants from the INMA (Environment and Childhood) Spanish cohort. Using high-performance liquid chromatography coupled to inductively coupled plasma-tandem mass, we measured arsenic (iAs and the metabolites monomethylarsonic acid (MMA) and dimethylarsinic acid (DMA)) in urine samples collected during the first trimester. iAs methylation efficiency was determined based on relative concentrations of the As metabolites in urine (%MMA, %DMA, and %iAs). Thirty-two single nucleotide polymorphisms (SNPs) in nine genes were determined in maternal DNA; AS3MT haplotypes were inferred. We assessed the association between genotypes/haplotypes and maternal As methylation efficiency using multivariate linear regression models. RESULTS: The median %MMA and %DMA were 5.3 %, and 89 %, respectively. Ancestral alleles of AS3MT SNPs (rs3740393, rs3740390, rs11191453, and rs11191454) were significantly associated with higher %MMA, %iAs, and lower %DMA. Pregnant participants with zero copies of the GGCTTCAC AS3MT haplotype presented a higher %MMA. Statistically significant associations were also found for the FOLH1 SNP rs202676 (ß 0.89 95%CI: 0.24, 1.55 for carriers of the G allele vs. the A allele). CONCLUSIONS: Our study shows that ancestral alleles in AS3MT polymorphisms were associated with lower As methylation efficiency in early pregnancy and suggests that FOLH1 also plays a role in As methylation efficiency. These results support the hypothesis that As metabolism is multigenic, being a key element for identifying susceptible populations.

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.

LncRNA ARGI Contributes to Virus-Induced Pancreatic ß Cell Inflammation Through Transcriptional Activation of IFN-Stimulated Genes.

Type 1 diabetes (T1D) is a complex autoimmune disease that develops in genetically susceptible individuals. Most T1D-associated single nucleotide polymorphisms (SNPs) are located in non-coding regions of the human genome. Interestingly, SNPs in long non-coding RNAs (lncRNAs) may result in the disruption of their secondary structure, affecting their function, and in turn, the expression of potentially pathogenic pathways. In the present work, the function of a virus-induced T1D-associated lncRNA named ARGI (Antiviral Response Gene Inducer) is characterized. Upon a viral insult, ARGI is upregulated in the nuclei of pancreatic ß cells and binds to CTCF to interact with the promoter and enhancer regions of IFNß and interferon-stimulated genes, promoting their transcriptional activation in an allele-specific manner. The presence of the T1D risk allele in ARGI induces a change in its secondary structure. Interestingly, the T1D risk genotype induces hyperactivation of type I IFN response in pancreatic ß cells, an expression signature that is present in the pancreas of T1D patients. These data shed light on the molecular mechanisms by which T1D-related SNPs in lncRNAs influence pathogenesis at the pancreatic ß cell level and opens the door for the development of therapeutic strategies based on lncRNA modulation to delay or avoid pancreatic ß cell inflammation in T1D.

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.