Dataset Information

Transcriptome-wide association study identified candidate genes associated with gut microbiota.

ABSTRACT:

Background

Gut microbiota is closely associated with host health and disease occurrence. Host genetic factor plays an important role in shaping gut microbial communities. The specific mechanism of host-regulated gene expression affecting gut microbiota has not been elucidated yet. Here we conducted a transcriptome-wide association study (TWAS) for gut microbiota by leveraging expression imputation from large-scale GWAS data sets.

Results

TWAS detected multiple tissue-specific candidate genes for gut microbiota, such as FUT2 for genus Bifidobacterium in transverse colon (P_PERM.ANL = 1.68 × 10^-3) and SFTPD for an unclassified genus of Proteobacteria in transverse colon (P_PERM.ANL = 5.69 × 10^-3). Fine mapping replicated 3 candidate genes in TWAS, such as HELLS for Streptococcus (PIP = 0.685) in sigmoid colon, ANO7 for Erysipelotrichaceae (PIP = 0.449) in sigmoid colon. Functional analyses detected 94 significant GO terms and 11 pathways for various taxa in total, such as GO_NUCLEOSIDE_DIPHOSPHATASE_ACTIVITY for Butyrivibrio (FDR P = 1.30 × 10^-4), KEGG_RENIN_ANGIOTENSIN_SYSTEM for Anaerostipes (FDR P = 3.16 × 10^-2). Literature search results showed 12 genes prioritized by TWAS were associated with 12 diseases. For instance, SFTPD for an unclassified genus of Proteobacteria was related to atherosclerosis, and FUT2 for Bifidobacterium was associated with Crohn's disease.

Conclusions

Our study results provided novel insights for understanding the genetic mechanism of gut microbiota, and attempted to provide clues for revealing the influence of genetic factors on gut microbiota for the occurrence and development of diseases.

SUBMITTER: Pan C

PROVIDER: S-EPMC8684646 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundOsteoarthritis (OA) is a common skeletal system disease that has been partially attributed to genetic factors. The hand is frequently affected, which seriously affects the patient's quality of life. However, the pathogenetic mechanism of hand osteoarthritis (hand OA) is still elusive.MethodsA genome-wide association study (GWAS) summary of hand OA was obtained from the UK Biobank dataset, which contains data from a total of 452,264 White British individuals, including 37,782 OA patients. The transcriptome-wide association study (TWAS) of hand OA was performed using FUnctional Summary-based ImputatiON (FUSION) with the skeletal muscle and blood as gene expression references. The significant genes identified by TWAS were further subjected to gene set enrichment analysis (GSEA) with the Database for Annotation, Visualization and Integrated Discovery (DAVID) tool. Furthermore, we compared the genes and gene sets identified by our TWAS with that of a knee OA mRNA expression profile to detect the genes and gene sets shared by TWAS and mRNA expression profiles in OA. The mRNA expression profiles of 18 normal knee cartilages and 20 OA knee cartilages were acquired from the Gene Expression Omnibus database (accession number: GSE114007).ResultsTWAS identified 177 genes with P < 0.05 for the skeletal muscle, including ANKRD44 (P = 0.0001), RIC3 (P = 0.0003), and AC005154.6 (P = 0.0004). TWAS identified 423 genes with P < 0.05 for the blood, including CRIM1 (P = 0.0002), ZNF880 (P = 0.0002), and NCKIPSD (P = 0.0003). After comparing the results of the TWAS to those of the mRNA expression profiling of OA, we identified 5 common genes, including DHRS3 (log2fold = - 1.85, P = 3.31 × 10- 9) and SKP2 (log2fold = 1.36, P = 1.62 × 10- 8). GSEA of TWAS identified 51 gene ontology (GO) terms for hand OA, for example, protein binding (P = 0.0003) and cytosol (P = 0.0020). We also detected 6 common GO terms shared by TWAS and mRNA expression profiling, including protein binding (PTWAS = 2.54 × 10- 4, PmRNA = 3.42 × 10- 8), extracellular exosome (PTWAS = 0.02, PmRNA = 1.18 × 10- 4), and cytoplasm (PTWAS = 0.0183, PmRNA = 0.0048).ConclusionIn this study, we identified 5 candidate genes and 6 GO terms related to hand OA, which may help to uncover the pathogenesis of hand OA. It should be noted that the possible difference in the gene expression profiles between hand OA and knee OA may affect our study results, which should be interpreted with caution.