A genome-wide gene-environmental interaction study identified novel loci for the relationship between ambient air pollution exposure and depression, anxiety.

BACKGROUND: Genetic factors and environmental exposures, including air pollution, contribute to the risk of depression and anxiety. While the association between air pollution and depression and anxiety has been established in the UK Biobank, there has been limited research exploring this relationship from a genetic perspective. METHODS: Based on individual genotypic and phenotypic data from a cohort of 104,385 participants in the UK Biobank, a polygenic risk score for depression and anxiety was constructed to explore the joint effects of nitric oxide (NO), nitrogen dioxide (NO2), particulate matter (PM) with a diameter of ⩽2.5 µm (PM2.5) and 2.5-10 µm (PMcoarse) with depression and anxiety by linear and logistic regression models. Subsequently, a genome-wide gene-environmental interaction study (GWEIS) was performed using PLINK 2.0 to identify the genes interacting with air pollution for depression and anxiety. RESULTS: A substantial risk of depression and anxiety development was detected in participants exposed to the high air pollution concomitantly with high genetic risk. GWEIS identified 166, 23, 18, and 164 significant candidate loci interacting with NO, NO2, PM2.5, and PMcoarse for Patient Health Questionnaire-9 (PHQ-9) score, and detected 44, 10, 10, and 114 candidate loci associated with NO, NO2, PM2.5, and PMcoarse for General Anxiety Disorder (GAD-7) score, respectively. And some significant genes overlapped among four air pollutants, like TSN (rs184699498, PNO2 = 3.47 × 10-9; rs139212326, PPM2.5 = 1.51 × 10-8) and HSP90AB7P(rs150987455, PNO2 = 1.63 × 10-11; rs150987455, PPM2.5 = 7.64 × 10-11), which were common genes affecting PHQ-9 score for both NO2 and PM2.5. CONCLUSION: Our study identified the joint effects of air pollution with genetic susceptibility on the risk of depression and anxiety, and provided several novel candidate genes for the interaction, contributing to an understanding of the genetic architecture of depression and anxiety.

Socioeconomic inequalities, genetic susceptibility, and risks of depression and anxiety: A large-observational study.

OBJECTIVES: This study aimed to investigate the interplay between genetic susceptibility and socioeconomic disparities on psychiatric disorders. METHODS: In this study, we utilized data from the UK Biobank to analyze the Generalized Anxiety Disorder (GAD)-7 scale (N = 74,425) and the Patient Health Questionnaire (PHQ)-9 (N = 74,101), along with the Index of Multiple Deprivation (IMD). The polygenic risk score (PRS) was calculated to assess the genetic risk associated with GAD-7/PHQ-9 scores, and the individuals were classified into low, medium, and high genetic risk groups according to tertiles of PRSs related to the GAD-7/PHQ-9. Linear regression models were used to explore the relationships between GAD-7/PHQ-9 scores and IMD scores in patients with different genetic susceptibilities. RESULTS: Disadvantaged socioeconomic status was associated with the risk of anxiety and depression across all strata of genetic risk, and stronger associations were shown for individuals with greater genetic susceptibility. From low to high genetic risk, the risk of psychiatric disorders increased for the GAD-7 (ß = 0.002 to 0.032) and PHQ-9 (ß = 0.003 to 0.045) scores. In addition, strong associations of high genetic risk with anxiety (ß = 0.875) and depression (ß = 1.152) were detected in the IMD quartile 4 group compared with the least deprivation quartile group. Furthermore, income and employment were estimated to contribute strongly to anxiety (ßemployment = 7.331, ßincome = 4.492) and depression (ßemployment = 9.951, ßincome = 6.453) in the high genetic risk group. CONCLUSION: The results suggest that we should pay more attention to psychiatric disorders with high genetic susceptibility and try to improve their socioeconomic status to prevent the development of psychiatric disorders.

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.

Air pollutants and the risk of incident hepatobiliary diseases: A large prospective cohort study in the UK Biobank.

The association between air pollutants and hepatobiliary pancreatic diseases remains inconclusive. This study analyzed up to 247,091 participants of White European ancestry (aged 37 to 73 years at recruitment) from the UK Biobank, a large-scale prospective cohort with open access. An air pollution score was utilized to assess the combined effect of PM2.5, PM2.5-10, PM10, NO2, and NOX on total hepatobiliary pancreatic diseases, liver diseases, cholecyst diseases, and pancreatic diseases. Cox proportional hazard models were employed to evaluate the relationships between air pollutants and the incidence of these diseases. Restricted cubic spline regressions were used to examine the dose-response association between air pollutants and the risk of hepatobiliary pancreatic diseases. We identified 4865 cases of total hepatobiliary pancreatic diseases, over a median follow-up of 10.86 years. The air pollution scores were moderately associated with increased liver disease risk (HR = 1.009, 95 % CI: 1.004, 1.014), but not with cholecyst and pancreatic diseases. Among the individual air pollutants, PM2.5 (HR = 1.069, 95 % CI: 1.025, 1.115) and PM10 (HR = 1.036, 95 % CI: 1.011, 1.061) significantly increased liver disease risk. Males showed a higher risk of liver diseases with PM2.5 (HR = 1.075, 95 % CI: 1.015, 1.139). Additionally, individuals with overweight (HR = 1.125, 95 % CI: 1.052, 1.203), age ≥ 60 and ≤73 (HR = 1.098, 95 % CI: 1.028, 1.172), and alcohol intake ≥ 14 unit/week (HR = 1.078, 95 % CI: 1.006, 1.155) had a higher risk of developing liver diseases at high expose to PM2.5. This study suggests that prolonged exposure to ambient air pollutants may elevate the risk of liver diseases.

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.

Disadvantaged social status contributed to sleep disorders: An observational and genome-wide gene-environment interaction analysis.

BACKGROUND: Sleep is a natural and essential physiological need for individuals. Our study aimed to research the associations between accumulated social risks and sleep disorders. METHODS: In this study, we came up with a polysocial risk score (PsRS), which is a cumulative social risk index composed of 13 social determinants of health. This research includes 239,165 individuals with sleep disorders and social determinants of health data from the UK Biobank cohort. First, logistic regression models were performed to examine the associations of social determinants of health and sleep disorders, including chronotype, narcolepsy, insomnia, snoring, short and long sleep duration. Then, PsRS was calculated based on statistically significant social determinants of health for each sleep disorder. Third, a genome-wide gene-environment interaction study was conducted to explore the interactions between single-nucleotide polymorphisms and PsRS in relation to sleep disorders. RESULTS: Higher PsRS scores were associated with worse sleep status, with the adjusted odds ratio (OR) ranging from 1.10 (95% Confidence interval [CI]: 1.09-1.11) to 1.29 (95% CI: 1.27-1.30) for sleep disorders. Emotional stress (OR = 1.36, 95% CI: 1.28-1.43) and not in paid employment (OR = 2.62, 95% CI: 2.51-2.74) were found to have significant contributions for sleep disorders. Moreover, multiple single-nucleotide polymorphisms were discovered to have interactions with PsRS, such as FRAS1 (P = 2.57 × 10-14) and CACNA1A (P = 8.62 × 10-14) for narcolepsy, and ACKR3 (P = 1.24 × 10-8) for long sleep. CONCLUSIONS: Our findings suggested that cumulative social risks was associated with sleep disorders, while the interactions between genetic susceptibility and disadvantaged social status are risk factors for the development of sleep disorders.

A genome-wide association study identifies candidate genes for sleep disturbances in depressed individuals.

OBJECTIVE: This study aimed to identify candidate loci and genes related to sleep disturbances in depressed individuals and clarify the co-occurrence of sleep disturbances and depression from the genetic perspective. METHODS: The study subjects (including 58,256 self-reported depressed individuals and 6,576 participants with PHQ-9 score ≥ 10, respectively) were collected from the UK Biobank, which were determined based on the Patient Health Questionnaire (PHQ-9) and self-reported depression status, respectively. Sleep related traits included chronotype, insomnia, snoring and daytime dozing. Genome-wide association studies (GWASs) of sleep related traits in depressed individuals were conducted by PLINK 2.0 adjusting age, sex, Townsend deprivation index and 10 principal components as covariates. The CAUSALdb database was used to explore the mental traits associated with the candidate genes identified by the GWAS. RESULTS: GWAS detected 15 loci significantly associated with chronotype in the subjects with self-reported depression, such as rs12736689 at RNASEL (P = 1.00 × 10- 09), rs509476 at RGS16 (P = 1.58 × 10- 09) and rs1006751 at RFX4 (P = 1.54 × 10- 08). 9 candidate loci were identified in the subjects with PHQ-9 ≥ 10, of which 2 loci were associated with insomnia such as rs115379847 at EVC2 (P = 3.50 × 10- 08), and 7 loci were associated with daytime dozing, such as rs140876133 at SMYD3 (P = 3.88 × 10- 08) and rs139156969 at ROBO2 (P = 3.58 × 10- 08). Multiple identified genes, such as RNASEL, RGS16, RFX4 and ROBO2 were reported to be associated with chronotype, depression or cognition in previous studies. CONCLUSION: Our study identified several candidate genes related to sleep disturbances in depressed individuals, which provided new clues for understanding the biological mechanism underlying the co-occurrence of depression and sleep disorders.

Identification of cell-specific epigenetic patterns associated with chondroitin sulfate treatment response in an endemic arthritis, Kashin-Beck disease.

Aims: To assess the alterations in cell-specific DNA methylation associated with chondroitin sulphate response using peripheral blood collected from Kashin-Beck disease (KBD) patients before initiation of chondroitin sulphate treatment. Methods: Peripheral blood samples were collected from KBD patients at baseline of chondroitin sulphate treatment. Methylation profiles were generated using reduced representation bisulphite sequencing (RRBS) from peripheral blood. Differentially methylated regions (DMRs) were identified using MethylKit, while DMR-related genes were defined as those annotated to the gene body or 2.2-kilobase upstream regions of DMRs. Selected DMR-related genes were further validated by quantitative reverse transcriptase polymerase chain reaction (qRT-PCR) to assess expression levels. Tensor composition analysis was performed to identify cell-specific differential DNA methylation from bulk tissue. Results: This study revealed 21,060 hypermethylated and 44,472 hypomethylated DMRs, and 13,194 hypermethylated and 22,448 hypomethylated CpG islands for differential global methylation for chondroitin sulphate treatment response. A total of 12,666 DMR-related genes containing DMRs were identified in their promoter regions, such as CHL1 (false discovery rate (FDR) = 2.11 × 10-11), RIC8A (FDR = 7.05 × 10-4), and SOX12 (FDR = 1.43 × 10-3). Additionally, RIC8A and CHL1 were hypermethylated in responders, while SOX12 was hypomethylated in responders, all showing decreased gene expression. The patterns of cell-specific differential global methylation associated with chondroitin sulphate response were observed. Specifically, we found that DMRs located in TESPA1 and ATP11A exhibited differential DNA methylation between responders and non-responders in granulocytes, monocytes, and B cells. Conclusion: Our study identified cell-specific changes in DNA methylation associated with chondroitin sulphate response in KBD patients.

Dietary diversity score and the acceleration of biological aging: a population-based study of 88,039 participants.

OBJECTIVES: Our study aimed to investigate the association of dietary diversity score (DDS), as reflected by five dietary categories, with biological age acceleration. DESIGN: A cross-sectional study. SETTING AND PARTICIPANTS: This study included 88,039 individuals from the UK Biobank. METHODS: Biological age (BA) was assessed using Klemerae-Doubal (KDM) and PhenoAge methods. The difference between BA and chronological age represents the age acceleration (AgeAccel), termed as "KDMAccel" and "PhenoAgeAccel". AgeAccel > 0 indicates faster aging. Generalized linear regression models were performed to assess the associations of DDS with AgeAccel. Similar analyses were performed for the five dietary categories. RESULTS: After adjusting for multiple variables, DDS was inversely associated with KDMAccel (ßHigh vs Low= -0.403, 95%CI: -0.492 to -0.314, P < 0.001) and PhenoAgeAccel (ßHigh vs Low= -0.545, 95%CI: -0.641 to -0.450, P < 0.001). Each 1-point increment in the DDS was associated with a 4.4% lower risk of KDMAccel and a 5.6% lower risk of PhenoAgeAccel. The restricted cubic spline plots demonstrated a non-linear dose-response association between DDS and the risk of AgeAccel. The consumption of grains (ßKDMAccel = -0.252, ßPhenoAgeAccel = -0.197), vegetables (ßKDMAccel = -0.044, ßPhenoAgeAccel = -0.077) and fruits (ßKDMAccel = -0.179, ßPhenoAgeAccel = -0.219) was inversely associated with the two AgeAccel, while meat and protein alternatives (ßKDMAccel = 0.091, ßPhenoAgeAccel = 0.054) had a positive association (All P < 0.001). Stratified analysis revealed stronger accelerated aging effects in males, smokers, and drinkers. A strengthening trend in the association between DDS and AgeAccel as TDI quartiles increased was noted. CONCLUSIONS: This study suggested that food consumption plays a role in aging process, and adherence to a higher diversity dietary is associated with the slowing down of the aging process.

Whole-Transcriptome Sequencing of Knee Joint Cartilage from Kashin-Beck Disease and Osteoarthritis Patients.

The aim of this study was to provide a comprehensive understanding of similarities and differences in mRNAs, lncRNAs, and circRNAs within cartilage for Kashin-Beck disease (KBD) compared to osteoarthritis (OA). We conducted a comparison of the expression profiles of mRNAs, lncRNAs, and circRNAs via whole-transcriptome sequencing in eight KBD and ten OA individuals. To facilitate functional annotation-enriched analysis for differentially expressed (DE) genes, DE lncRNAs, and DE circRNAs, we employed bioinformatic analysis utilizing Gene Ontology (GO) and KEGG. Additionally, using quantitative reverse transcriptase polymerase chain reaction (qRT-PCR), we validated the expression levels of four cartilage-related genes in chondrocytes. We identified a total of 43 DE mRNAs, 1451 DE lncRNAs, and 305 DE circRNAs in KBD cartilage tissue compared to OA (q value < 0.05; |log2FC| > 1). We also performed competing endogenous RNA network analysis, which identified a total of 65 lncRNA-mRNA interactions and 4714 miRNA-circRNA interactions. In particular, we observed that circRNA12218 had binding sites for three miRNAs targeting ACAN, while circRNA12487 had binding sites for seven miRNAs targeting COL2A1. Our results add a novel set of genes and non-coding RNAs that could potentially serve as candidate diagnostic biomarkers or therapeutic targets for KBD patients.

Enhanced polygenic risk score incorporating gene-environment interaction suggests the association of major depressive disorder with cardiac and lung function.

BACKGROUND: Depression has been linked to an increased risk of cardiovascular and respiratory diseases; however, its impact on cardiac and lung function remains unclear, especially when accounting for potential gene-environment interactions. METHODS: We developed a novel polygenic and gene-environment interaction risk score (PGIRS) integrating the major genetic effect and gene-environment interaction effect of depression-associated loci. The single nucleotide polymorphisms (SNPs) demonstrating major genetic effect or environmental interaction effect were obtained from genome-wide SNP association and SNP-environment interaction analyses of depression. We then calculated the depression PGIRS for non-depressed individuals, using smoking and alcohol consumption as environmental factors. Using linear regression analysis, we assessed the associations of PGIRS and conventional polygenic risk score (PRS) with lung function (N = 42 886) and cardiac function (N = 1791) in the subjects with or without exposing to smoking and alcohol drinking. RESULTS: We detected significant associations of depression PGIRS with cardiac and lung function, contrary to conventional depression PRS. Among smokers, forced vital capacity exhibited a negative association with PGIRS (ß = -0.037, FDR = 1.00 × 10-8), contrasting with no significant association with PRS (ß = -0.002, FDR = 0.943). In drinkers, we observed a positive association between cardiac index with PGIRS (ß = 0.088, FDR = 0.010), whereas no such association was found with PRS (ß = 0.040, FDR = 0.265). Notably, in individuals who both smoked and drank, forced expiratory volume in 1-second demonstrated a negative association with PGIRS (ß = -0.042, FDR = 6.30 × 10-9), but not with PRS (ß = -0.003, FDR = 0.857). CONCLUSIONS: Our findings underscore the profound impact of depression on cardiac and lung function, highlighting the enhanced efficacy of considering gene-environment interactions in PRS-based studies.

Mendelian randomization study of the relationship between blood and urine biomarkers and schizophrenia in the UK Biobank cohort.

BACKGROUND: The identification of suitable biomarkers is of crucial clinical importance for the early diagnosis of treatment-resistant schizophrenia (TRS). This study aims to comprehensively analyze the association between TRS and blood and urine biomarkers. METHODS: Candidate TRS-related single nucleotide polymorphisms (SNPs) were obtained from a recent genome-wide association study. The UK Biobank cohort, comprising 376,807 subjects with blood and urine biomarker testing data, was used to calculate the polygenic risk score (PRS) for TRS. Pearson correlation analyses were performed to evaluate the correlation between TRS PRS and each of the biomarkers, using calculated TRS PRS as the instrumental variables. Bidirectional two-sample Mendelian randomization (MR) was used to assess potential causal associations between candidate biomarkers with TRS. RESULTS: Here we identify a significant association between TRS PRS and phosphate (r = 0.007, P = 1.96 × 10-4). Sex subgroup analyses identify seven and three candidate biomarkers associated with TRS PRS in male and female participants, respectively. For example, total protein and phosphate for males, creatinine and phosphate for females. Bidirectional two-sample MR analyses indicate that TRS is negatively associated with cholesterol (estimate = -0.363, P = 0.008). Conversely, TRS is positively associated with total protein (estimate = 0.137, P = 0.027), mean corpuscular volume (estimate = 0.032, P = 2.25 × 10-5), and mean corpuscular hemoglobin (estimate = 0.018, P = 0.007). CONCLUSIONS: Our findings provide insights into the roles of blood and urine biomarkers in the early detection and treatment of TRS.

People with schizophrenia experience periods of time during which they misperceive reality. Some people with schizophrenia do not respond well to the usual drugs that are used to relieve their symptoms. This type of schizophrenia is known as treatment-resistant schizophrenia (TRS). We looked at differences in the genes (inherited characteristics), blood and urine of a group of people in the UK with schizophrenia to see if people with TRS have particular characteristics that would enable them to be distinguished from patients with schizophrenia who tend to respond to usual treatment. We found several differences in the blood that could be used to predict which people might get TRS, including some that were specific to men or women. These discoveries are important because they can help doctors identify people who are more likely to develop TRS earlier, enabling them to avoid using treatments that might not work well for them.

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.

Long-term ambient air pollution and the risk of musculoskeletal diseases: A prospective cohort study.

Evidence of the associations of air pollution and musculoskeletal diseases is inconsistent. This study aimed to examine the associations between air pollutants and the risk of incident musculoskeletal diseases, such as degenerative joint diseases (n = 38,850) and inflammatory arthropathies (n = 20,108). An air pollution score was constructed to assess the combined effect of PM2.5, PM2.5-10, NO2, and NOX. Cox proportional hazard model was applied to assess the relationships between air pollutants and the incidence of each musculoskeletal disease. The air pollution scores exhibited the modest association with an increased risk of osteoporosis (HR = 1.006, 95% CI: 1.002-1.011). Among the individual air pollutants, PM2.5 and PM2.5-10 exhibited the most significant effect on elevated risk of musculoskeletal diseases, such as PM2.5 on osteoporosis (HR = 1.064, 95% CI: 1.020-1.110), PM2.5-10 on inflammatory arthropathies (HR = 1.059, 95% CI: 1.037-1.081). Females were found to have a higher risk of incident musculoskeletal diseases when exposed to air pollutants. Individuals with extreme BMI or lower socioeconomic status had a higher risk of developing musculoskeletal diseases. Our findings reveal that long-term exposure to ambient air pollutants may contribute to an increased risk of musculoskeletal diseases.

Involvement of Yes-Associated Protein 1 Activation in the Matrix Degradation of Human-Induced-Pluripotent-Stem-Cell-Derived Chondrocytes Induced by T-2 Toxin and Deoxynivalenol Alone and in Combination.

T-2 toxin and deoxynivalenol (DON) are two prevalent mycotoxins that cause cartilage damage in Kashin-Beck disease (KBD). Cartilage extracellular matrix (ECM) degradation in chondrocytes is a significant pathological feature of KBD. It has been shown that the Hippo pathway is involved in cartilage ECM degradation. This study aimed to examine the effect of YAP, a major regulator of the Hippo pathway, on the ECM degradation in the hiPS-derived chondrocytes (hiPS-Ch) model of KBD. The hiPS-Ch injury models were established via treatment with T-2 toxin/DON alone or in combination. We found that T-2 toxin and DON inhibited the proliferation of hiPS-Ch in a dose-dependent manner; significantly increased the levels of YAP, SOX9, and MMP13; and decreased the levels of COL2A1 and ACAN (all p values < 0.05). Immunofluorescence revealed that YAP was primarily located in the nuclei of hiPS-Ch, and its expression level increased with toxin concentrations. The inhibition of YAP resulted in the dysregulated expression of chondrogenic markers (all p values < 0.05). These findings suggest that T-2 toxin and DON may inhibit the proliferation of, and induce the ECM degradation, of hiPS-Ch mediated by YAP, providing further insight into the cellular and molecular mechanisms contributing to cartilage damage caused by toxins.

Associations of classical HLA alleles with depression and anxiety.

Immune dysregulation has been widely observed in patients with psychiatric disorders. This study aims to examine the association between HLA alleles and depression and anxiety. Using data from the UK Biobank, we performed regression analyses to assess the association of 359 HLA alleles with depression and anxiety, as determined by Patient Health Questionnaire (PHQ) score (n = 120,033), self-reported depression (n = 121,685), general anxiety disorder (GAD-7) score (n = 120,590), and self-reported anxiety (n = 108,310). Subsequently, we conducted gene environmental interaction study (GEIS) to evaluate the potential effects of interactions between HLA alleles and environmental factors on the risk of depression and anxiety. Sex stratification was implemented in all analysis. Our study identified two significant HLA alleles associated with self-reported depression, including HLA-C*07:01 (ß = -0.015, p = 5.54 × 10-5 ) and HLA-B*08:01 (ß = -0.015, p = 7.78 × 10-5 ). Additionally, we identified four significant HLA alleles associated with anxiety score, such as HLA-DRB1*07:01 (ß = 0.084, p = 9.28 × 10-5 ) and HLA-B*57:01 (ß = 0.139, p = 1.22 × 10-4 ). GEIS revealed that certain HLA alleles interacted with environmental factors to influence mental health outcomes. For instance, HLA-A*02:07 × cigarette smoking was associated with depression score (ß = 0.976, p = 1.88 × 10-6 ). Moreover, sex stratification analysis revealed significant sex-based differences in the interaction effects of certain HLA alleles with environmental factors. Our findings indicate the considerable impact of HLA alleles on the risks of depression and anxiety, providing valuable insights into the functional relevance of immune dysfunction in these conditions.

Identification of novel rare variants for anxiety: an exome-wide association study in the UK Biobank.

BACKGROUND: Rare variants are believed to play a substantial role in the genetic architecture of mental disorders, particularly in coding regions. However, limited evidence supports the impact of rare variants on anxiety. METHODS: Using whole-exome sequencing data from 200,643 participants in the UK Biobank, we investigated the contribution of rare variants to anxiety. Firstly, we computed genetic risk score (GRS) of anxiety utilizing genotype data and summary data from a genome-wide association study (GWAS) on anxiety disorder. Subsequently, we identified individuals within the lowest 50% GRS, a subgroup more likely to carry pathogenic rare variants. Within this subgroup, we classified individuals with the highest 10% 7-item Generalized Anxiety Disorder scale (GAD-7) score as cases (N = 1869), and those with the lowest 10% GAD-7 score were designated as controls (N = 1869). Finally, we conducted gene-based burden tests and single-variant association analyses to assess the relationship between rare variants and anxiety. RESULTS: Totally, 47,800 variants with MAF ≤0.01 were annotated as non-benign coding variants, consisting of 42,698 nonsynonymous SNVs, 489 nonframeshift substitution, 236 frameshift substitution, 617 stop-gain and 40 stop-loss variants. After variation aggregation, 5066 genes were included in gene-based association analysis. Totally, 11 candidate genes were detected in burden test, such as RNF123 (PBonferroni adjusted = 3.40 × 10-6), MOAP1(PBonferroni adjusted = 4.35 × 10-4), CCDC110 (PBonferroni adjusted = 5.83 × 10-4). Single-variant test detected 9 rare variants, such as rs35726701(RNF123)(PBonferroni adjusted = 3.16 × 10-10) and rs16942615(CAMTA2) (PBonferroni adjusted = 4.04 × 10-4). Notably, RNF123, CCDC110, DNAH2, and CSKMT gene were identified in both tests. CONCLUSIONS: Our study identified novel candidate genes for anxiety in protein-coding regions, revealing the contribution of rare variants to anxiety.

A multidimensional social risk atlas of depression and anxiety: An observational and genome-wide environmental interaction study.

Background: Mental disorders are largely socially determined, yet the combined impact of multidimensional social factors on the two most common mental disorders, depression and anxiety, remains unclear. Methods: We constructed a polysocial risk score (PsRS), a multidimensional social risk indicator including components from three domains: socioeconomic status, neighborhood and living environment and psychosocial factors. Supported by the UK Biobank cohort, we randomly divided 110 332 participants into the discovery cohort (60%; n = 66 200) and the replication cohort (40%; n = 44 134). We tested the associations between 13 single social factors with Patient Health Questionnaire (PHQ) score, Generalized Anxiety Disorder Scale (GAD) score and self-reported depression and anxiety. The significant social factors were used to calculate PsRS for each mental disorder by considering weights from the multivariable linear model. Generalized linear models were applied to explore the association between PsRS and depression and anxiety. Genome-wide environmental interaction study (GWEIS) was further performed to test the effect of interactions between PsRS and SNPs on the risk of mental phenotypes. Results: In the discovery cohort, PsRS was positively associated with PHQ score (ß = 0.37; 95% CI = 0.35-0.38), GAD score (ß = 0.27; 95% CI = 0.25-0.28), risk of self-reported depression (OR = 1.29; 95% CI = 1.28-1.31) and anxiety (OR = 1.19; 95% CI = 1.19-1.23). Similar results were observed in the replication cohort. Emotional stress, lack of social support and low household income were significantly associated with the development of depression and anxiety. GWEIS identified multiple candidate loci for PHQ score, such as rs149137169 (ST18) (Pdiscovery = 1.08 × 10-8, Preplication = 3.25 × 10-6) and rs3759812 (MYO9A) (Pdiscovery = 3.87 × 10-9, Preplication = 6.21 × 10-5). Additionally, seven loci were detected for GAD score, such as rs114006170 (TMPRSS11D) (Pdiscovery = 1.14 × 10-9, Preplication = 7.36 × 10-5) and rs77927903 (PIP4K2A) (Pdiscovery = 2.40 × 10-9, Preplication = 0.002). Conclusions: Our findings reveal the positive effects of multidimensional social factors on the risk of depression and anxiety. It is important to address key social disadvantage in mental health promotion and treatment.

Insight into the Causal Relationship between Gut Microbiota and Back Pain: A Two Sample Bidirectional Mendelian Randomization Study.

Observational studies have shown that alterations in gut microbiota composition are associated with low back pain. However, it remains unclear whether the association is causal. To reveal the causal association between gut microbiota and low back pain, a two-sample bidirectional Mendelian randomization (MR) analysis is performed. The inverse variance weighted regression (IVW) is performed as the principal MR analysis. MR-Egger and Weighted Median is further conducted as complementary analysis to validate the robustness of the results. Finally, a reverse MR analysis is performed to evaluate the possibility of reverse causation. The inverse variance weighted (IVW) method suggests that Peptostreptococcaceae (odds ratio [OR] 1.056, 95% confidence interval [CI] [1.015-1.098], P IVW = 0.010), and Lactobacillaceae (OR 1.070, 95% CI [1.026-1.115], P IVW = 0.003) are positively associated with back pain. The Ruminococcaceae (OR 0.923, 95% CI [0.849-0.997], P IVW = 0.033), Butyricicoccus (OR 0.920, 95% CI [0.868 - 0.972], P IVW = 0.002), and Lachnospiraceae (OR 0.948, 95% CI [0.903-0.994], P IVW = 0.022) are negatively associated with back pain. In this study, underlying causal relationships are identified among gut microbiota and low back pain. Notably, further research is needed on the biological mechanisms by which gut microbiota influences low back pain.