Aims: Despite the interest in the association of gut microbiota with bone health, limited population-based studies of gut microbiota and bone mineral density (BMD) have been made. Our aim is to explore the possible association between gut microbiota and BMD. Methods: A total of 3,321 independent loci of gut microbiota were used to calculate the individual polygenic risk score (PRS) for 114 gut microbiota-related traits. The individual genotype data were obtained from UK Biobank cohort. Linear regressions were then conducted to evaluate the possible association of gut microbiota with L1-L4 BMD (n = 4,070), total BMD (n = 4,056), and femur total BMD (n = 4,054), respectively. PLINK 2.0 was used to detect the single-nucleotide polymorphism (SNP) × gut microbiota interaction effect on the risks of L1-L4 BMD, total BMD, and femur total BMD, respectively. Results: We detected five, three, and seven candidate gut microbiota-related traits for L1-L4 BMD, total BMD, and femur BMD, respectively, such as genus Dialister (p = 0.004) for L1-L4 BMD, and genus Eisenbergiella (p = 0.046) for total BMD. We also detected two common gut microbiota-related traits shared by L1-L4 BMD, total BMD, and femur total BMD, including genus Escherichia Shigella and genus Lactococcus. Interaction analysis of BMD detected several genes that interacted with gut microbiota, such as phospholipase D1 ( PLD1) and endomucin ( EMCN) interacting with genus Dialister in total BMD, and COL12A1 and Discs Large MAGUK Scaffold Protein 2 ( DLG2) interacting with genus Lactococcus in femur BMD. Conclusion: Our results suggest associations between gut microbiota and BMD, which will be helpful to further explore the regulation mechanism and intervention gut microbiota of BMD. Cite this article: Bone Joint Res 2021;10(11):734–741.

Previous studies have suggested that vitamin D (VD) was associated with psychiatric diseases, but efforts to elucidate the functional relevance of VD with depression and anxiety from genetic perspective have been limited. Based on the UK Biobank cohort, we first calculated polygenic risk score (PRS) for VD from genome-wide association study (GWAS) data of VD. Linear and logistic regression analysis were conducted to evaluate the associations of VD traits with depression and anxiety traits, respectively. Then, using individual genotype and phenotype data from the UK Biobank, genome-wide environment interaction studies (GWEIS) were performed to identify the potential effects of gene × VD interactions on the risks of depression and anxiety traits. In the UK Biobank cohort, we observed significant associations of blood VD level with depression and anxiety traits, as well as significant associations of VD PRS and depression and anxiety traits. GWEIS identified multiple candidate loci, such as rs114086183 (p = 4.11 × 10⁻⁸, LRRTM4) for self-reported depression status and rs149760119 (p = 3.88 × 10⁻⁸, GNB5) for self-reported anxiety status. Our study results suggested that VD was negatively associated with depression and anxiety. GWEIS identified multiple candidate genes interacting with VD, providing novel clues for understanding the biological mechanism potential associations between VD and psychiatric disorders.

The aim of this study was to investigate the possible interaction of mitochondrial dysfunction and inflammatory cytokines in the risk of anxiety and depression. We utilized the UK Biobank for the sample of this study. A mitochondria-wide association(MiWAS) and interaction analysis was performed to investigate the interaction effects of mitochondrial DNA (mtDNA)×C-reactive protein (CRP) on the risks of self-reported anxiety (N = 72,476), general anxiety disorder (GAD-7) scores (N = 80,853), self-reported depression (N = 80,778), Patient Health Questionnaire (PHQ-9) scores (N = 80,520) in total samples, females and males, respectively, adjusting for sex, age, Townsend deprivation index (TDI), education score, alcohol intake, smoking and 10 principal components. In all, 25 mtSNPs and 10 mtSNPs showed significant level of association with self-reported anxiety and GAD-7 score respectively. A total of seven significant mtDNA × CRP interactions were found for anxiety, such as m.3915G>A(MT-ND1) for self-reported anxiety in total subjects (P = 6.59 × 10⁻³), m.4561T>C(MT-ND2) (P = 3.04 × 10⁻³) for GAD-7 score in total subjects. For depression, MiWAS identified 17 significant mtSNPs for self-reported depression and 14 significant mtSNPs for PHQ-9 scores. 17 significant mtDNA associations (2 for self-reported depression and 15 for PHQ-9 score) was identified, such as m.14869G>A(MT-CYB; P = 2.22 × 10⁻³) associated with self-reported depression and m.4561T>C (MT-ND2; P value = 3.02 × 10⁻⁸) associated with PHQ-9 score in all subjects. In addition, 5 common mtDNA shared with anxiety and depression were found in MiWAS, and 4 common mtDNA variants were detected to interact with CRP for anxiety and depression, such as m.9899T>C(MT-CO3). Our study suggests the important interaction effects of mitochondrial function and CRP on the risks of anxiety and depression.

Multiple previous studies explored the associations between early life factors and the age at onset of the depression. However, they only focused on the influence of environmental or genetic factors, without considering the interactions between them. Based on previous genome-wide association study (GWAS) data, we first calculated polygenic risk score (PRS) for depression. Regression analyses were conducted to assess the interacting effects of depression PRS and 5 early life factors, including felt hated by family member (N = 40,112), physically abused by family (N = 40,464), felt loved (N = 35633), and sexually molested (N = 41,595) in childhood and maternal smoking during pregnancy (N = 38,309), on the age at onset of the depression. Genome-wide environment interaction studies (GWEIS) were then performed to identify the genes interacting with early life factors for the age at onset of the depression. In regression analyses, we observed significant interacting effects of felt loved as a child and depression PRS on the age at onset of depression in total sample (β = 0.708, P = 5.03 × 10⁻³) and males (β = 1.421, P = 7.64 × 10⁻⁴). GWEIS identified a novel candidate loci interacting with felt loved as a child at GSAP (rs2068031, P = 4.24 × 10^–8) and detected several genes with suggestive significance association, such as CMYA5 (rs7343, P = 2.03 × 10^–6) and KIRREL3 (rs535603, P = 4.84 × 10^–6) in males. Our results indicate emotional care in childhood may affect the age at onset of depression, especially in males, and GSAP plays an important role in their interaction.