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.

Potential involvement of connective tissue growth factor in chondrocytes apoptosis of Kashin-Beck disease.

BACKGROUND: Kashin-Beck disease (KBD) is an endemic osteoarthropathy characterized by excessive chondrocytes apoptosis. T-2 toxin exposure has been proved to be its etiology. Connective tissue growth factor (CTGF) exerts a profound influence on cartilage growth and metabolism. We investigated the potential role of CTGF in KBD development and examined CTGF alterations under T-2 toxin stimulation. METHODS: The levels of CTGF and chondrocyte apoptosis-related markers in cartilage and primary chondrocytes from KBD and control groups were measured using qRT-PCR, Western blotting, immunohistochemistry, and immunofluorescence. We analyzed expression changes of these genes in response to T-2 toxin. Apoptosis rates of chondrocytes induced by T-2 toxin were measured by flow cytometry and TUNEL assay. The active pharmaceutical ingredient targeting CTGF was screened through Comparative Toxicogenomics Database, and molecular docking was performed using AutoDock Tools. RESULTS: The CTGF levels in KBD cartilage and chondrocytes were significantly elevated and positively associated with the levels of apoptosis-related genes. T-2 toxin exposure increased CTGF and apoptosis-related gene levels in chondrocytes, with apoptosis rates rising alongside T-2 toxin concentration. Curcumin was identified as targeting CTGF and exhibited effective binding. It could down-regulate CTGF, apoptosis-related genes, such as Cleaved caspase 3 and BAX, and also significantly reduce apoptosis rate in chondrocytes treated with T-2 toxin. CONCLUSION: CTGF plays a crucial role in the development of KBD. Curcumin has shown potential in inhibiting CTGF levels and reducing chondrocyte apoptosis, highlighting its promise as a therapeutic agent for preventing cartilage damage in KBD. Our findings provided valuable insights into the pathogenesis of KBD and could promote the development of novel therapeutic strategies for this debilitating disease.

Sex differences in the metabolic activation of and platelet response to vicagrel in mice: Androgen as a key player.

As a biological variable, sex influences the metabolism of and/or response to certain drugs. Vicagrel is being developed as an investigational new drug in China; however, it is unknown whether sex could affect its metabolic activation and platelet responsiveness. This study aimed to determine whether such differences could exist, and to elucidate the mechanisms involved. Orchiectomized (ORX) or ovariectomized (OVX) mouse models were used to investigate the effects of androgen or estrogen on the metabolic activation of and platelet response to vicagrel. Plasma vicagrel active metabolite H4 concentrations, platelet inhibition of vicagrel, and protein levels of intestinal hydrolases Aadac and Ces2 were measured, respectively. Further, p38-MAPK signaling pathway was enriched, whose role was determined using SB202190. Results showed that female mice exhibited significantly elevated systemic exposure of H4 and enhanced platelet responses to vicagrel than males, and protein expression levels of Aadac and Ces2 differed by sex. OVX mice exhibited less changes than sham mice. ORX mice exhibited increases in protein levels of intestinal hydrolases, systemic exposure of H4, and platelet inhibition of vicagrel, but dihydrotestosterone (DHT) reversed these changes in ORX mice and suppressed these changes in OVX mice. Phosphorylated p38 levels were reduced in female or ORX mice but increased in ORX mice by DHT. SB202190 reversed DHT-induced changes observed in ORX mice. We concluded that sex differences exist in metabolic activation of and platelet response to vicagrel in mice through elevation of p38 phosphorylation by androgen, suggesting sex-based vicagrel dosage adjustments for patient care.

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.

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.

Activation of Notch Signaling Pathway is involved in Extracellular Matrix Degradation in human induced pluripotent stem cells chondrocytes induced by HT-2 toxin.

Notch signaling regulates cartilage formation and homeostasis. Kashin-Beck Disease (KBD), an endemic osteochondropathy, is characterized by severe cartilage degradation. The etiology of KBD is related to the exposure of HT-2 toxin, a mycotoxin and primary metabolite of T-2 toxin. This study aims to explore the role of HT-2 toxin in the Notch signaling regulation and extracellular matrix (ECM) metabolism of hiPSCs-Chondrocytes. Immunohistochemistry and qRT-PCR were employed to investigate the expression of Notch pathway molecules in KBD articular cartilage and primary chondrocytes. hiPSCs-Chondrocytes, derived from hiPSCs, were treated with 100 ng/mL HT-2 toxin and the γ-secretase inhibitor (DAPT) for 48h, respectively. The markers related to the Notch signaling pathway and ECM were assessed using qRT-PCR and Western blot. Notch pathway dysregulation was prominent in KBD cartilage. HT-2 toxin exposure caused cytotoxicity in hiPSCs-Chondrocytes, and activated Notch signaling by increasing the mRNA and protein levels of NOTCH1 and HES1. HT-2 toxin also upregulated ECM catabolic enzymes and downregulated ECM components (COL2A1 and ACAN), indicating ECM degradation. DAPT-mediated Notch signaling inhibition suppressed the mRNA and protein level of ADAMTS5 expression while enhancing ECM component expression in hiPSCs-Chondrocytes. This study suggests that HT-2 toxin may induce ECM degradation in hiPSCs-Chondrocytes through activating Notch signaling.

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.

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.

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.

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.

Modulatory interactions of T-2 and deoxynivalenol mycotoxins on murine femoral development and osteological integrity.

In this study, we conducted a systematic assessment of the effectsof deoxynivalenol (DON) and T-2 mycotoxins (T-2) on the developmental processes and structural integrity of murine femurs, considering both the isolated and synergistic effects of these toxins. To this end, we divided 72 male mice into nine groups, each subjected to varying dosages of T-2, DON, or their combinations. Over a four-week experimental period, meticulous monitoring was undertaken regarding the mice's body weight, biochemical markers of bone formation and resorption, and the activity of relevant cells. To comprehensively evaluate alterations in bone structure, we employed biomechanical analysis, micro-computed tomography (micro-CT), and transmission electron microscopy.Our findings unveiled a significant revelation: the mice exhibited a dose-dependent decrease in body weight upon exposure to individual mycotoxins, while the combined use of these toxins manifested an atypical antagonistic effect. Furthermore, we observed variations in the levels of calcium, phosphorus, and vitamin D, as well as adjustments in the activities of osteoblasts and osteoclasts, all intricately linked to the dosage and ratio of the toxins. Alterations in biomechanical properties were also noted to correlate with the dosage and combination of toxins. Analyses via micro-CT and transmission electron microscopy further corroborated the substantial impact of toxin dosage and combinations on both cortical and trabecular bone structures.In summation, our research unequivocally demonstrates the dose- and ratio-dependent detrimental effects of DON and T-2 mycotoxins on the growth and structural integrity of murine femurs. These insights accentuate the importance of a profound understanding of the potential risks these toxins pose to bone health, offering pivotal guidance for future toxicological research and public health preventative strategies.

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.

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.

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.

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.

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.