Gut microbiota accelerates obesity in peri-/post-menopausal women via Bacteroides fragilis and acetic acid.

OBJECTIVE: Many animal experiments and epidemiological studies have shown that the gut microbiota (GM) plays an important role in the development of obesity, but the specific biological mechanism involved in the pathogenesis of disease remain unknown. We aimed to examine the relationships and functional mechanisms of GM on obesity in peri- and post-menopausal women. METHODS: We recruited 499 Chinese peri- and post-menopausal women and performed comprehensive analyses of the gut microbiome, targeted metabolomics for short-chain fatty acids in serum, and host whole-genome sequencing by various association analysis methods. RESULTS: Through constrained linear regression analysis, we found that an elevated abundance of Bacteroides fragilis (B. fragilis) was associated with obesity. We also found that serum levels of acetic acid were negatively associated with obesity, and that B. fragilis was negatively associated with serum acetic acid levels by partial Spearman correlation analysis. Mendelian randomization analysis indicated that B. fragilis increases the risk of obesity and may causally down-regulate acetic acid levels. CONCLUSIONS: We found the gut with B. fragilis may accelerate obesity, in part, by suppressing acetic acid levels. Therefore, B. fragilis and acetic acid may represent important therapeutic targets for obesity intervention in peri- and post-menopausal women.

Associations between Genetically Proxied Inhibition of Lipid-Lowering Drug Targets and Serum Micronutrients among Individuals of European Descent: A Mendelian Randomization Study.

BACKGROUND: Limited and inconclusive data exist concerning the associations between lipid-lowering drugs and serum micronutrient concentrations. METHODS: We conducted Mendelian randomization (MR) analyses to explore the associations between lipid-lowering drug targets and serum micronutrients. Single-nucleotide polymorphisms in genes encoding molecular targets of LDL cholesterol-lowering therapies were selected as instrumental variables for 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGCR; target of statins), proprotein convertase subtilisin/kexin type 9 (PCSK9; target of PCSK9 inhibitors), and Niemann-Pick C1-Like 1 (NPC1L1; target of ezetimibe). Exposure data were extracted from a published genome-wide association study (GWAS) of lipids in 188,577 European individuals, with outcome data obtained from the Integrative Epidemiology Unit (IEU) GWAS database (https://gwas.mrcieu.ac.uk). Overall, age and sex information were not calculable from the summary-level GWAS data. MR analyses were performed using the inverse-variance weighted method and MR sensitivity analysis methods. RESULTS: We found genetically proxied inhibition of HMGCR to lower iron (effect, -0.16; 95% CI: -0.27, -0.06; P value = 0.003), zinc (effect, -0.83; 95% CI: -1.36, -0.31; P value = 0.002), magnesium (effect, -0.17; 95% CI: -0.27, -0.06; P value = 0.003), potassium (effect, -0.17; 95% CI: -0.27, -0.06; P value = 0.002), genetically proxied inhibition of NPC1L1 to increase calcium (effect, 0.28; 95% CI: 0.10, 0.46; P value = 0.003), retinol (effect, 0.25; 95% CI: 0.07, 0.44; P value = 0.01), and genetically proxied inhibition of PCSK9 to increase vitamin D (effect, 0.10; 95% CI: 0.07, 0.12; P value = 1.8 × 10-19). These associations were robust in MR sensitivity analyses. However, the associations between genetically proxied inhibition of HMGCR and NPC1L1 and the micronutrients were not consistent in multiple comparisons. CONCLUSIONS: Our results provide evidence that statin use may lower serum concentrations of iron, zinc, magnesium, and potassium, PCSK9 inhibitors may increase serum vitamin D, and ezetimibe may increase serum calcium and retinol concentrations.

Improved detection of genetic loci in estimated glomerular filtration rate and type 2 diabetes using a pleiotropic cFDR method.

Genome-wide association studies (GWAS) have been shown to have the potential of explaining more of the "missing heritability" of complex human phenotypes by improving statistical approaches. Here, we applied a genetic-pleiotropy-informed conditional false discovery rate (cFDR) to capture additional polygenic effects associated with estimated glomerular filtration rate (creatinine) (eGFRcrea) and type 2 diabetes (T2D). The cFDR analysis improves the identification of pleiotropic variants by incorporating potentially shared genetic mechanisms between two related traits. The Q-Q and fold-enrichment plots were used to assess the enrichment of SNPs associated with eGFRcrea or T2D, and Manhattan plots were used for showing chromosomal locations of the significant loci detected. By applying the cFDR method, we newly identified 74 loci for eGFRcrea and 3 loci for T2D with the cFDR criterion of 0.05 compared with previous related GWAS studies. Four shared SNPs were detected to be associated with both eGFRcrea and T2D at the significant conjunction cFDR level of 0.05, and these shared SNPs had not been reported in previous studies. In addition, we used DAVID analysis to perform functional analysis of the shared SNPs' annotated genes and found their potential hidden associations with eGFRcrea and T2D. In this study, the cFDR method shows the feasibility to detect more genetic variants underlying the heritability of eGFRcrea and T2D, and the overlapping SNPs identified could be regarded as candidate loci that provide a thread of genetic mechanisms between eGFRcrea and T2D in future research.

Novel common variants associated with body mass index and coronary artery disease detected using a pleiotropic cFDR method.

Genome-wide association studies (GWAS) have been successfully applied in identifying single nucleotide polymorphisms (SNPs) associated with body mass index (BMI) and coronary heart disease (CAD). However, the SNPs to date can only explain a small percentage of the genetic variances of traits. Here, we applied a genetic pleiotropic conditional false discovery rate (cFDR) method that combines summary statistic p values from different multi-center GWAS datasets, to detect common genetic variants associated with these two traits. The enrichment of SNPs associated with BMI and CAD was assessed by conditional Q-Q plots and the common variants were identified by the cFDR method. By applying the cFDR level of 0.05, 7 variants were identified to be associated with CAD (2 variants being novel), 34 variants associated with BMI (11 variants being novel), and 3 variants associated with both BMI and CAD (2 variants being novel). The SNP rs653178 (ATXN2) is noteworthy as this variant was replicated in an independent analysis. SNP rs12411886 (CNNM2) and rs794356 (HIP1) were of note as the annotated genes may be associated with processes that are functionally important in lipid metabolism. In conclusion, the cFDR method identified novel variants associated with BMI and/or CAD by effectively incorporating different GWAS datasets.

Enteric nervous system damage caused by abnormal intestinal butyrate metabolism may lead to functional constipation.

Functional constipation (FC) is a high morbidity gastrointestinal disease for which dysfunction in the enteric nervous system is a major pathogenesis mechanism. To enhance our understanding of the involvement of intestinal microbiota and its metabolites in the pathogenesis of FC, we conducted a shotgun metagenomic sequencing analysis of gut microbiota and serum short-chain fatty acids (SCFAs) analysis in 460 Chinese women with different defecation frequencies. We observed that the abundance ofFusobacterium_varium, a butyric acid-producing bacterium, was positively correlated (P = 0.0096) with the frequency of defecation; however, the concentrations of serum butyric acid was negatively correlated (P = 3.51E-05) with defecation frequency. These results were verified in an independent cohort (6 patients with FC and 6 controls). To further study the effects of butyric acid on intestinal nerve cells, we treated mouse intestinal neurons in vitro with various concentrations of butyrate (0.1, 0.5, 1, and 2.5 mM). We found that intestinal neurons treated with 0.5 mM butyrate proliferated better than those in the other treatment groups, with significant differences in cell cycle and oxidative phosphorylation signal pathways. We suggest that the decreased butyrate production resulting from the reduced abundance of Fusobacterium in gut microbiota affects the proliferation of intestinal neurons and the energy supply of intestinal cells. However, with FC disease advancing, the consumption and excretion of butyric acid reduce, leading to its accumulation in the intestine. Moreover, the accumulation of an excessively high amount of butyric acid inhibits the proliferation of nerve cells and subsequently exacerbates the disease.

Gut microbiota impacts bone via Bacteroides vulgatus-valeric acid-related pathways.

Although the gut microbiota has been reported to influence osteoporosis risk, the individual species involved, and underlying mechanisms, remain largely unknown. We performed integrative analyses in a Chinese cohort of peri-/post-menopausal women with metagenomics/targeted metabolomics/whole-genome sequencing to identify novel microbiome-related biomarkers for bone health. Bacteroides vulgatus was found to be negatively associated with bone mineral density (BMD), which was validated in US white people. Serum valeric acid (VA), a microbiota derived metabolite, was positively associated with BMD and causally downregulated by B. vulgatus. Ovariectomized mice fed B. vulgatus demonstrated increased bone resorption and poorer bone micro-structure, while those fed VA demonstrated reduced bone resorption and better bone micro-structure. VA suppressed RELA protein production (pro-inflammatory), and enhanced IL10 mRNA expression (anti-inflammatory), leading to suppressed maturation of osteoclast-like cells and enhanced maturation of osteoblasts in vitro. The findings suggest that B. vulgatus and VA may represent promising targets for osteoporosis prevention/treatment.

A bi-directional Mendelian randomization study of the sarcopenia-related traits and osteoporosis.

Both sarcopenia and osteoporosis are common geriatric diseases causing huge socioeconomic burdens, and clinically, they often occur simultaneously. Observational studies have found a controversial correlation between sarcopenia and osteoporosis and their causal relationship is not clear. Therefore, we performed a bi-directional two-sample Mendelian randomization (MR) analysis to assess the potential causal relationship between sarcopenia-related traits (hand grip strength, lean mass, walking pace) and osteoporosis. Our analysis was performed by applying genetic variants obtained from the UK Biobank and the GEnetic Factors for OSteoporosis (GEFOS) datasets. We used inverse-variance weighted (IVW) and several sensitivity analyses to estimate and cross-validate the potential causal relationship in this study. We found that bone mineral density (BMD) was causally positively associated with left-hand grip strength (ß = 0.017, p-value = 0.001), fat-free mass (FFM; right leg FFM, ß = 0.014, p-value = 0.003; left arm FFM, ß = 0.014, p-value = 0.005), but not walking pace. Higher hand grip strength was potentially causally associated with increased LS-BMD (right-hand grip strength, ß = 0.318, p-value = 0.001; left-hand grip strength, ß = 0.358, p-value = 3.97 × 10-4). In conclusion, osteoporosis may be a risk factor for sarcopenia-related traits and muscle strength may have a site-specific effect on BMD.

Metagenomic study of the gut microbiota associated with cow milk consumption in Chinese peri-/postmenopausal women.

Cow milk consumption (CMC) and alterations of gut bacterial composition are proposed to be closely related to human health and disease. Our research aims to investigate the changes in human gut microbial composition in Chinese peri-/postmenopausal women with different CMC habits. A total of 517 subjects were recruited and questionnaires about their CMC status were collected; 394 subjects were included in the final analyses. Fecal samples were used for studying gut bacterial composition. All the subjects were divided into a control group (n = 248) and a CMC group (n = 146) according to their CMC status. Non-parametric tests and LEfSe at different taxonomic levels were used to reveal differentially abundant taxa and functional categories across different CMC groups. Relative abundance (RA) of one phylum (p_Actinobacteria), three genera (g_Bifidobacterium, g_Anaerostipes, and g_Bacteroides), and 28 species diversified significantly across groups. Specifically, taxa g_Anaerostipes (p < 0.01), g_Bacteroides (p < 0.05), s_Anaerostipes_hadrus (p < 0.01), and s_Bifidobacterium_pseudocatenulatum (p < 0.01) were positively correlated with CMC levels, but p_Actinobacteria (p < 0.01) and g_Bifidobacterium (p < 0.01) were negatively associated with CMC levels. KEGG module analysis revealed 48 gut microbiome functional modules significantly (p < 0.05) associated with CMC, including Vibrio cholerae pathogenicity signature, cholera toxins (p = 9.52e-04), and cephamycin C biosynthesis module (p = 0.0057), among others. In conclusion, CMC was associated with changes in gut microbiome patterns including beta diversity and richness of some gut microbiota. The alterations of certain bacteria including g_Anaerostipes and s_Bifidobacterium_pseudocatenulatum in the CMC group should be important for human health. This study further supports the biological value of habitual cow milk consumption.

Human gut microbiome impacts skeletal muscle mass via gut microbial synthesis of the short-chain fatty acid butyrate among healthy menopausal women.

BACKGROUND: Increasing evidence suggests that human gut microbiome plays an important role in variation of skeletal muscle mass (SMM). However, specific causal mechanistic relationship of human gut microbiome with SMM remains largely unresolved. Understanding the causal mechanistic relationship may provide a basis for novel interventions for loss of SMM. This study investigated whether human gut microbiome has a causal effect on SMM among Chinese community-dwelling healthy menopausal women. METHODS: Estimated SMM was derived from whole-body dual-energy X-ray absorptiometry. We performed integrated analyses on whole-genome sequencing, shotgun metagenomic sequencing, and serum short-chain fatty acids (SCFAs), as well as available host SMM measurements among community-dwelling healthy menopausal women (N = 482). We combined the results with summary statistics from genome-wide association analyses for human gut microbiome (N = 952) and SMM traits (N = 28 330). As a prerequisite for causality, we used a computational protocol that was proposed to measure correlations among gut metagenome, metabolome, and the host trait to investigate the relationship between human gut microbiome and SMM. Causal inference methods were applied to assess the potential causal effects of gut microbial features on SMM, through one-sample and two-sample Mendelian randomization (MR) analyses, respectively. RESULTS: In metagenomic association analyses, the increased capacity for gut microbial synthesis of the SCFA butyrate was significantly associated with serum butyrate levels [Spearman correlation coefficient (SCC) = 0.13, P = 0.02] and skeletal muscle index (SCC = 0.084, P = 0.002). Of interest was the finding that two main butyrate-producing bacterial species were both positively associated with the increased capacity for gut microbial synthesis of butyrate [Faecalibacterium prausnitzii (SCC = 0.25, P = 6.6 × 10-7 ) and Butyricimonas virosa (SCC = 0.15, P = 0.001)] and for skeletal muscle index [F. prausnitzii (SCC = 0.16, P = 6.2 × 10-4 ) and B. virosa (SCC = 0.17, P = 2.4 × 10-4 )]. One-sample MR results showed a causal effect between gut microbial synthesis of the SCFA butyrate and appendicular lean mass (ß = 0.04, 95% confidence interval 0.029 to 0.051, P = 0.003). Two-sample MR results further confirmed the causal effect between gut microbial synthesis of the SCFA butyrate and appendicular lean mass (ß = 0.06, 95% confidence interval 0 to 0.13, P = 0.06). CONCLUSIONS: Our results may help the future development of novel intervention approaches for preventing or alleviating loss of SMM.

Sarcopenia-related traits and coronary artery disease: a bi-directional Mendelian randomization study.

Previous Mendelian randomization (MR) studies have yielded a conflicting causal relationship between sarcopenia and coronary artery disease (CAD), and lack the association of CAD with sarcopenia. We performed a bi-directional MR approach to clarify the causality and causal direction between sarcopenia-related traits and CAD. In stage 1 analysis, estimates of inverse variance weighting (IVW) and several sensitivity analyses were obtained by applying genetic variants that predict sarcopenia-related traits to CAD. Conversely, we also applied genetic variants that predict CAD to sarcopenia-related traits in stage 2 analyses. IVW analysis showed that higher handgrip strength reduces risk for CAD: A 1-kilogram (kg) increase in genetically determined left handgrip strength reduced odds of CAD by 36% [odds ratio (OR) = 0.64, 95% confidence interval (CI) 0.498 - 0.821, p = 4.56E-04], and right handgrip strength reduced odds of CAD by 41.1% (OR = 0.599, 95% CI 0.476 - 0.753, p = 1.10E-05). However, genetically predicted CAD did not show any causal association with handgrip strength, and no significant causal relationship was detected between genetically instrumented body lean mass and CAD. Our results suggest that decreased muscle strength but not decreased muscle mass leads to the increased risk of CAD in sarcopenia.

Causal Effects of Genetically Predicted Cardiovascular Risk Factors on Chronic Kidney Disease: A Two-Sample Mendelian Randomization Study.

Observational studies have demonstrated that cardiovascular risk factors are associated with chronic kidney disease (CKD). However, these observational associations are potentially influenced by the residual confounding, including some unmeasured lifestyle factors and interaction risk factors. Two-sample mendelian randomization analysis was conducted in this study to evaluate whether genetically predicted cardiovascular risk factors have a causal effect on the risk of CKD. We selected genetic variants associated with cardiovascular risk factors and extracted the corresponding effect sizes from the largest GWAS summary-level dataset of CKD. Cardiovascular risk factors contain high density lipoprotein (HDL) cholesterol, low density lipoprotein (LDL) cholesterol, total cholesterol (TC), triglyceride (TG), glycated hemoglobin (HbA1c), fasting glucose, systolic blood pressure (SBP) and diastolic blood pressure (DBP). A Bonferroni corrected threshold of P = 0.006 was considered as significant, and 0.006 < P < 0.05 was considered suggestive of evidence for a potential association. Genetically predicted DBP was significantly associated with CKD [odds ratio (OR) was 1.35 (95% confidence interval (CI) (1.10, 1.65); P = 0.004)]. There was suggestive evidence for potential associations between genetically predicted higher HDL cholesterol [OR: 0.88, 95%CI (0.80, 0.98), P = 0.025] and lower adds of CKD, and between higher SBP [OR: 1.36, 95%CI (1.07, 1.73), P = 0.013] and higher adds of CKD. However, genetically predicted LDL cholesterol, TC, TG, HbA1c, and fasting glucose did not show any causal association with CKD.

Additional common variants associated with type 2 diabetes and coronary artery disease detected using a pleiotropic cFDR method.

Genome-wide association studies (GWASs) have been performed extensively in diverse populations to identify single nucleotide polymorphisms (SNPs) associated with complex diseases or traits. However, to date, the SNPs identified fail to explain a large proportion of the variance of the traits/diseases. GWASs on type 2 diabetes (T2D) and coronary artery disease (CARD) are generally performed as single-trait studies, rather than analyzing the related traits simultaneously. Despite the extensive evidence suggesting that these two phenotypes share both genetic and environmental risk factors, the shared overlapping genetic biological mechanisms between these traits remain largely unexplored. Here, we adopted a recently developed genetic pleiotropic conditional false discovery rate (cFDR) approach to discover novel loci associated with T2D and CARD by incorporating the summary statistics from existing GWASs of these two traits. Applying the cFDR level of 0.05, 33 loci were identified for T2D and 34 loci for CARD, 9 of which for both. By incorporating pleiotropic effects into a conditional analysis framework, we observed that there is significant pleiotropic enrichment between T2D and CARD. These findings may provide novel insights into the etiology of T2D and CARD, as well as the processes that may influence disease development both individually and jointly.

Genetically driven adiposity traits increase the risk of coronary artery disease independent of blood pressure, dyslipidaemia, glycaemic traits.

Adiposity has been associated with the risk of coronary artery disease (CAD) in observational studies, but their association may differ according to specific characteristics of studies. In Mendelian randomization (MR) analyses, genetic variants are used as instrumental variables (IVs) of exposures to examine causal effects to overcome confounding factors and reverse causation. We performed MR analyses for adiposity (n = 322,154) on risk of CAD (60,801 cases and 123,504 controls) based on the currently largest genome-wide association studies. The estimated associations between adiposity traits and CAD were calculated by an inverse-variance weighted method with and without excluding the IVs, which are associated with the well-known risk factors of CAD. Genetic variants are identified to be associated with the well-known risk factors of CAD by a cross-phenotype meta-analysis method. Our results furnished strong evidence for a causal role of adiposity in risk of CAD, with the odds ratios (ORs) for CAD being 1.53 (95% CI 1.36-1.72) for body mass index (BMI), 1.48 (1.20-1.96) for waist-hip ratio (WHR), and 1.34 (1.07-1.59) for WHR adjusted for BMI (WHRadjBMI), respectively. After excluding mediators-associated IVs from the MR analyses, the corresponding ORs were 1.46 (1.28-1.67) for BMI, 1.39 (1.01-1.93) for WHR, and 1.38 (1.04-1.84) for WHRadjBMI, respectively. Furthermore, our results suggested that central adiposity and general adiposity might pose a similar risk for CAD. In summary, our data supported that genetically driven adiposity traits imposed the risk of CAD independent of blood pressure, dyslipidaemia, glycaemic traits, and type 2 diabetes.

Assessing causal estimates of the association of obesity-related traits with coronary artery disease using a Mendelian randomization approach.

Obesity-related traits have been associated with coronary artery disease (CAD) in observational studies, but these associations may be biased by confounding factors and reverse causation. In this study, we specifically conducted two-sample Mendelian randomization (MR) analyses to overcome these limitations and test the associations of obesity-related traits (other than body mass index (BMI)) (n = 322,154) with CAD (22,233 cases and 64,762 controls) by using summary-level data from previous studies. The methods utilized to estimate these associations included the inverse-variance weighted method, the weighted median method and MR-Egger regression. Our results supported causal effects of BMI, hip circumference (HC), waist circumference (WC), and waist-hip ratio (WHR) on CAD. The associations of BMI-adjusted HC and WC with CAD were reversed, unlike that of WHR. In MR analyses excluding overlapping single nucleotide polymorphisms (SNPs) from obesity-related traits, the associations of these traits with CAD were preserved. The associations of BMI-adjusted HC and WC with CAD require further investigation, as collider stratification may be occurring. Additionally, central adiposity (measured by WHR) separated from general adiposity (measured by BMI) and general adiposity might pose similar risks for CAD. In clinical practice, physicians should pay attention to the potential effects of different obesity-related traits on CAD.

Novel Common Variants Associated with Obesity and Type 2 Diabetes Detected Using a cFDR Method.

Genome-wide association studies (GWASs) have been performed extensively in diverse populations to identify single nucleotide polymorphisms (SNPs) associated with complex diseases or traits. However, to date, the SNPs identified fail to explain a large proportion of the variance of the traits/diseases. GWASs on type 2 diabetes (T2D) and obesity are generally focused on individual traits independently, and genetic intercommunity (common genetic contributions or the product of over correlated phenotypic world) between them are largely unknown, despite extensive data showing that these two phenotypes share both genetic and environmental risk factors. Here, we applied a recently developed genetic pleiotropic conditional false discovery rate (cFDR) approach to discover novel loci associated with BMI and T2D by incorporating the summary statistics from existing GWASs of these two traits. Conditional Q-Q and fold enrichment plots were used to visually demonstrate the strength of pleiotropic enrichment. Adopting a cFDR nominal significance level of 0.05, 287 loci were identified for BMI and 75 loci for T2D, 23 of which for both traits. By incorporating related traits into a conditional analysis framework, we observed significant pleiotropic enrichment between obesity and T2D. These findings may provide novel insights into the etiology of obesity and T2D, individually and jointly.