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.