Single-cell multiomics analysis reveals cell/tissue-specific associations in bipolar disorder.

This study investigates the cellular origin and tissue heterogeneity in bipolar disorder (BD) by integrating multiomics data. Four distinct datasets were employed, including single-cell RNA sequencing (scRNA-seq) data (embryonic and fetal brain, n = 8, 1,266 cells), BD Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) data (adult brain, n = 210), BD bulk RNA-seq data (adult brain, n = 314), and BD genome-wide association study (GWAS) summary data (n = 413,466). The integration of scRNA-seq data with multiomics data relevant to BD was accomplished using the single-cell disease relevance score (scDRS) algorithm. We have identified a novel brain cell cluster named ADCY1, which exhibits distinct genetic characteristics. From a high-resolution genetic perspective, glial cells emerge as the primary cytopathology associated with BD. Specifically, astrocytes were significantly related to BD at the RNA-seq level, while microglia showed a strong association with BD across multiple panels, including the transcriptome-wide association study (TWAS), ATAC-seq, and RNA-seq. Additionally, oligodendrocyte precursor cells displayed a significant association with BD in both ATAC-seq and RNA-seq panel. Notably, our investigation of brain regions affected by BD revealed significant associations between BD and all three types of glial cells in the dorsolateral prefrontal cortex (DLPFC). Through comprehensive analyses, we identified several BD-associated genes, including CRMP1, SYT4, UCHL1, and ZBTB18. In conclusion, our findings suggest that glial cells, particularly in specific brain regions such as the DLPFC, may play a significant role in the pathogenesis of BD. The integration of multiomics data has provided valuable insights into the etiology of BD, shedding light on potential mechanisms underlying this complex psychiatric disorder.

Evaluating the associations between intelligence quotient and multi-tissue proteome from the brain, CSF and plasma.

Intelligence quotient is a vital index to evaluate the ability of an individual to think rationally, learn from experience and deal with the environment effectively. However, limited efforts have been paid to explore the potential associations of intelligence quotient traits with the tissue proteins from the brain, CSF and plasma. The information of protein quantitative trait loci was collected from a recently released genome-wide association study conducted on quantification data of proteins from the tissues including the brain, CSF and plasma. Using the individual-level genotypic data from the UK Biobank cohort, we calculated the polygenic risk scores for each protein based on the protein quantitative trait locus data sets above. Then, Pearson correlation analysis was applied to evaluate the relationships between intelligence quotient traits (including 120 330 subjects for 'fluid intelligence score' and 38 949 subjects for 'maximum digits remembered correctly') and polygenic risk scores of each protein in the brain (17 protein polygenic risk scores), CSF (116 protein polygenic risk scores) and plasma (59 protein polygenic risk scores). The Bonferroni corrected P-value threshold was P < 1.30 × 10-4 (0.05/384). Finally, Mendelian randomization analysis was conducted to test the causal relationships between 'fluid intelligence score' and pre-specific proteins from correlation analysis results. Pearson correlation analysis identified significant association signals between the protein of macrophage-stimulating protein and fluid intelligence in brain and CSF tissues (P brain = 1.21 × 10-8, P CSF = 1.10 × 10-7), as well as between B-cell lymphoma 6 protein and fluid intelligence in CSF (P CSF = 1.23 × 10-4). Other proteins showed close-to-significant associations with the trait of 'fluid intelligence score', such as plasma protease C1 inhibitor (P CSF = 4.19 × 10-4, P plasma = 6.97 × 10-4), and with the trait of 'maximum digits remembered correctly', such as tenascin (P plasma = 3.42 × 10-4). Additionally, Mendelian randomization analysis results suggested that macrophage-stimulating protein (Mendelian randomization-Egger: ß = 0.54, P = 1.64 × 10-61 in the brain; ß = 0.09, P = 1.60 × 10-12 in CSF) had causal effects on fluid intelligence score. We observed functional relevance of specific tissue proteins to intelligence quotient and identified several candidate proteins, such as macrophage-stimulating protein. This study provided a novel insight to the relationship between tissue proteins and intelligence quotient traits.

Body shape from birth to adulthood is associated with skeletal development: A Mendelian randomization study.

BACKGROUND: Observational studies have shown that childhood obesity is associated with adult bone health but yield inconsistent results. We aimed to explore the potential causal association between body shape and skeletal development. METHODS: We used two-sample Mendelian randomization (MR) to estimate causal relationships between body shape from birth to adulthood and skeletal phenotypes, with exposures including placental weight, birth weight, childhood obesity, BMI, lean mass, fat mass, waist circumference, and hip circumference. Independent genetic instruments associated with the exposures at the genome-wide significance level (P < 5 × 10-8) were selected from corresponding large-scale genome-wide association studies. The inverse-variance weighted analysis was chosen as the primary method, and complementary MR analyses included the weighted median, MR-Egger, weighted mode, and simple mode. RESULTS: The MR analysis shows strong evidence that childhood (ß = -1.29 × 10-3, P = 8.61 × 10-5) and adulthood BMI (ß = -1.28 × 10-3, P = 1.45 × 10-10) were associated with humerus length. Tibiofemoral angle was negatively associated with childhood BMI (ß = -3.60 × 10-1, P = 3.00 × 10-5) and adolescent BMI (ß = -3.62 × 10-1, P = 2.68 × 10-3). In addition, genetically predicted levels of appendicular lean mass (ß = 1.16 × 10-3, P = 1.49 × 10-13), whole body fat mass (ß = 1.66 × 10-3, P = 1.35 × 10-9), waist circumference (ß = 1.72 × 10-3, P = 6.93 × 10-8) and hip circumference (ß =1.28 × 10-3, P = 4.34 × 10-6) were all associated with tibia length. However, we found no causal association between placental weight, birth weight and bone length/width. CONCLUSIONS: This large-scale MR analysis explores changes in growth patterns in the length/width of major bone sites, highlighting the important role of childhood body shape in bone development and providing insights into factors that may drive bone maturation.

An atlas of causal association between micronutrients and osteoarthritis.

OBJECTIVE: This study examines the causal relationships between serum micronutrients and site-specific osteoarthritis (OA) using Mendelian Randomization (MR). METHODS: This study performed a two-sample MR analysis to explore causal links between 21 micronutrients and 11 OA outcomes. These outcomes encompass overall OA, seven site-specific manifestations, and three joint replacement subtypes. Sensitivity analyses using MR methods, such as the weighted median, MR-Egger, and MR-PRESSO, assessed potential horizontal pleiotropy and heterogeneity. Genome-wide association summary statistical data were utilized for both exposure and outcome data, including up to 826,690 participants with 177,517 OA cases. All data was sourced from Genome-wide association studies datasets from 2009 to 2023. RESULTS: In the analysis of associations between 21 micronutrients and 11 OA outcomes, 15 showed Bonferroni-corrected significance (P < 0.000216), without significant heterogeneity or horizontal pleiotropy. Key findings include strong links between gamma-tocopherol and spine OA (OR = 1.70), and folate with hand OA in finger joints (OR = 1.15). For joint replacements, calcium showed a notable association with a reduced likelihood of total knee replacement (TKR) (OR = 0.52) and total joint replacement (TJR) (OR = 0.56). Serum iron was significantly associated with an increased risk of total hip replacement (THR) (OR = 1.23), while folate indicated a protective effect (OR = 0.95). Various sex-specific associations were also uncovered. CONCLUSION: These findings underscore the critical role of micronutrients in osteoarthritis, providing valuable insights for preventive care and potential enhancement of treatment outcomes.

Assessing the joint effects of mitochondrial genes and physical activity on the psychiatric phenotype of subjective well-being based on the UK Biobank data.

Subjective well-being (SWB) is an important measure for mental health status. Previous research has shown that physical activity can affect an individual's well-being, yet the underlying molecular mechanism remains to be clarified. In this study, we aim to evaluate the potential interactions between mitochondrial genes and physical activity (PA) as well as their combined effects on individual well-being. SWB phenotype data in UK Biobank were enrolled for this study including nine aspects such as work/job satisfaction, health satisfaction, family relationship satisfaction, friendships satisfaction, financial situation satisfaction, ever depressed for a whole week, general happiness, general happiness with own health and belief that own life is meaningful. We made analysis for each aspects separately. Firstly, mitochondria-wide association studies (MiWAS) was conducted to assess the association of mitochondrial Single Nucleotide Polymorphisms SNP with each aspect of SWB. Then an interaction analysis of mitochondrial DNA (mtDNA) mutation and PA was performed to evaluate their joint effect on SWB status. Meanwhile, these two analysis were made for female and male group separately as well as the total samples, all under the control of possible confounding factors including gender, age, Townsend Deprivation Index (TDI), education, alcohol consumption, smoking habits, and 10 principal components. MiWAS analysis identified 45 mtSNPs associated with 9 phenotypes of SWB. For example, m.15218A > G on MT-CYB in the health satisfaction phenotype of the total subjects. Gender-specific analyses found 30 mtSNPs in females and 58 in males, involving 13 mtGenes. In mtDNA-PA interaction analysis, we also identified 10 significant mtDNA-PA interaction sets for SWB. For instance, m.13020 T > C (MT-ND5) was associated with the SWB financial situation satisfaction phenotype in all subjects (P = 0.00577). In addition, MiWAS analysis identified 12 mtGene variants associated with SWB, as MT-ND1 and MT-ND2. However, in mtDNA-PA interactions we detected 7 mtDNA affecting psychiatric disorders occurring, as in the friendships satisfaction phenotype (m.3394 T > C on MT-ND1). Our study results suggest an implication of the interaction between mitochondrial function and physical activity in the risk of psychiatric disorder development.

Assessing the interaction effects of mitochondrial DNA polymorphisms and lifestyle on heel bone mineral density.

BACKGROUND: Bone mineral density (BMD) is a major predictor of osteoporotic fractures, and previous studies have reported the effects of mitochondrial dysfunction and lifestyle on BMD, respectively. However, their interaction effects on BMD are still unclear. Therefore, we aimed to investigate the possible interaction of mitochondrial DNA (mtDNA) and common lifestyles contributing to osteoporosis. METHODS: Our analysis included 119,120 white participants (Nfemale=65,949 and Nmale=53,171) from the UK Biobank with heel BMD phenotype data. A generalized linear regression model of PLINK was performed to assess the interaction effects of mtDNA and five life environmental factors on heel BMD, including smoking, drinking, physical activity, dietary diversity score, and vitamin D. In addition, we also performed linear regression analysis for total body BMD. Finally, we assessed the potential causal relationships between mtDNA copy number (mtDNA-CN) and life environmental factors using Mendelian randomization (MR) analysis. RESULTS: Our study identified four mtDNA loci showing suggestive evidence of heel BMD, such as m.16356T>C (MT-DLOOP; P =1.50×10-3) in total samples. Multiple candidate mtDNA×lifetsyle interactions were also detected for heel BMD, such as MT-ND2×physical activity (P = 2.88×10-3) in total samples and MT-ND1×smoking (P = 8.54×10-4) in males. Notably, MT-CYB was a common candidate mtDNA loci for heel BMD to interact with five life environmental factors. Multivariable MR analysis indicated a causal effect of physical activity on heel BMD when mtDNA-CN was considered (P =1.13×10-3). CONCLUSIONS: Our study suggests the candidate interaction between mitochondria and lifestyles on heel BMD, providing novel clues for exploring the pathogenesis of osteoporosis.

Identifying novel chemical-related susceptibility genes for five psychiatric disorders through integrating genome-wide association study and tissue-specific 3'aQTL annotation datasets.

The establishment of 3'aQTLs comprehensive database provides an opportunity to help explore the functional interpretation from the genome-wide association study (GWAS) data of psychiatric disorders. In this study, we aim to search novel susceptibility genes, pathways, and related chemicals of five psychiatric disorders via GWAS and 3'aQTLs datasets. The GWAS datasets of five psychiatric disorders were collected from the open platform of Psychiatric Genomics Consortium (PGC, https://www.med.unc.edu/pgc/ ) and iPSYCH ( https://ipsych.dk/ ) (Demontis et al. in Nat Genet 51(1):63-75, 2019; Grove et al. in Nat Genet 51:431-444, 2019; Genomic Dissection of Bipolar Disorder and Schizophrenia in Cell 173: 1705-1715.e1716, 2018; Mullins et al. in Nat Genet 53: 817-829; Howard et al. in Nat Neurosci 22: 343-352, 2019). The 3'untranslated region (3'UTR) alternative polyadenylation (APA) quantitative trait loci (3'aQTLs) summary datasets of 12 brain regions were obtained from another public platform ( https://wlcb.oit.uci.edu/3aQTLatlas/ ) (Cui et al. in Nucleic Acids Res 50: D39-D45, 2022). First, we aligned the GWAS-associated SNPs of psychiatric disorders and datasets of 3'aQTLs, and then, the GWAS-associated 3'aQTLs were identified from the overlap. Second, gene ontology (GO) and pathway analysis was applied to investigate the potential biological functions of matching genes based on the methods provided by MAGMA. Finally, chemical-related gene-set analysis (GSA) was also conducted by MAGMA to explore the potential interaction of GWAS-associated 3'aQTLs and multiple chemicals in the mechanism of psychiatric disorders. A number of susceptibility genes with 3'aQTLs were found to be associated with psychiatric disorders and some of them had brain-region specificity. For schizophrenia (SCZ), HLA-A showed associated with psychiatric disorders in all 12 brain regions, such as cerebellar hemisphere (P = 1.58 × 10-36) and cortex (P = 1.58 × 10-36). GO and pathway analysis identified several associated pathways, such as Phenylpropanoid Metabolic Process (GO:0009698, P = 6.24 × 10-7 for SCZ). Chemical-related GSA detected several chemical-related gene sets associated with psychiatric disorders. For example, gene sets of Ferulic Acid (P = 6.24 × 10-7), Morin (P = 4.47 × 10-2) and Vanillic Acid (P = 6.24 × 10-7) were found to be associated with SCZ. By integrating the functional information from 3'aQTLs, we identified several susceptibility genes and associated pathways especially chemical-related gene sets for five psychiatric disorders. Our results provided new insights to understand the etiology and mechanism of psychiatric disorders.