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

Abstract

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. 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. 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. 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.