Genome-wide associations of human gut microbiome variation and implications for causal inference analyses

David A. Hughes,Rodrigo Bacigalupe,Jun Wang,Malte C. Rühlemann,Raul Y. Tito,Gwen Falony,Marie Joossens,Sara Vieira-Silva,Liesbet Henckaerts,Leen Rymenans,Chloë Verspecht,Susan Ring,Andre Franke,Kaitlin H. Wade,Nicholas J. Timpson,Jeroen Raes

DOI: https://doi.org/10.1038/s41564-020-0743-8

IF: 28.3

2020-06-22

Nature Microbiology

Abstract:Recent population-based^{<a href="#ref-CR1">1</a>,<a href="#ref-CR2">2</a>,<a href="#ref-CR3">3</a>,<a href="/articles/s41564-020-0743-8#ref-CR4">4</a>} and clinical studies^{<a href="/articles/s41564-020-0743-8#ref-CR5">5</a>} have identified a range of factors associated with human gut microbiome variation. Murine quantitative trait loci^{<a href="/articles/s41564-020-0743-8#ref-CR6">6</a>}, human twin studies^{<a href="/articles/s41564-020-0743-8#ref-CR7">7</a>} and microbiome genome-wide association studies^{<a href="/articles/s41564-020-0743-8#ref-CR1">1</a>,<a href="/articles/s41564-020-0743-8#ref-CR3">3</a>,<a href="#ref-CR8">8</a>,<a href="#ref-CR9">9</a>,<a href="#ref-CR10">10</a>,<a href="#ref-CR11">11</a>,<a href="/articles/s41564-020-0743-8#ref-CR12">12</a>} have provided evidence for genetic contributions to microbiome composition. Despite this, there is still poor overlap in genetic association across human studies. Using appropriate taxon-specific models, along with support from independent cohorts, we show an association between human host genotype and gut microbiome variation. We also suggest that interpretation of applied analyses using genetic associations is complicated by the probable overlap between genetic contributions and heritable components of host environment. Using faecal 16S ribosomal RNA gene sequences and host genotype data from the Flemish Gut Flora Project (<i>n</i> = 2,223) and two German cohorts (FoCus, <i>n</i> = 950; PopGen, <i>n</i> = 717), we identify genetic associations involving multiple microbial traits. Two of these associations achieved a study-level threshold of <i>P</i> = 1.57 × 10⁻¹⁰; an association between <i>Ruminococcus</i> and <a href="https://www.ncbi.nlm.nih.gov/snp/?term=rs150018970">rs150018970</a> near <i>RAPGEF1</i> on chromosome 9, and between <i>Coprococcus</i> and <a href="https://www.ncbi.nlm.nih.gov/snp/?term=rs561177583">rs561177583</a> within <i>LINC01787</i> on chromosome 1. Exploratory analyses were undertaken using 11 other genome-wide associations with strong evidence for association (<i>P</i> &lt; 2.5 × 10⁻⁸) and a previously reported signal of association between <a href="https://www.ncbi.nlm.nih.gov/snp/?term=rs4988235">rs4988235</a> (<i>MCM6</i>/<i>LCT</i>) and <i>Bifidobacterium</i>. Across these 14 single-nucleotide polymorphisms there was evidence of signal overlap with other genome-wide association studies, including those for age at menarche and cardiometabolic traits. Mendelian randomization analysis was able to estimate associations between microbial traits and disease (including <i>Bifidobacterium</i> and body composition); however, in the absence of clear microbiome-driven effects, caution is needed in interpretation. Overall, this work marks a growing catalogue of genetic associations that will provide insight into the contribution of host genotype to gut microbiome. Despite this, the uncertain origin of association signals will likely complicate future work looking to dissect function or use associations for causal inference analysis.

microbiology