Intra-islet glucagon signalling regulates beta-cell connectivity, first-phase insulin secretion and glucose homeostasis

BACKGROUND

Type 2 diabetes (T2D) is characterised by the loss of first-phase insulin secretion. We studied mice with β-cell selective loss of the glucagon receptor (Gcgr fl/fl X Ins-1Cre), to investigate the role of intra-islet glucagon receptor signalling on pan-islet calcium activity and insulin secretion.

METHODS

Metabolic profiling was conducted on Gcgr β-cell-/- and littermate controls. Crossing with GCaMP6f (STOP flox) animals further allowed for β-cell specific expression of a fluorescent calcium indicator. These islets were functionally imaged in vitro and in vivo. Wild-type mice were transplanted with islets expressing GCaMP6f in β-cells into the anterior eye chamber and placed on a high fat diet. Part of the cohort received a glucagon analogue (GCG-analogue) for 40 days and the control group were fed to achieve weight matching. Calcium imaging was performed regularly during the development of hyperglycaemia and in response to GCG-analogue treatment.

RESULTS

Gcgr β-cell-/- mice exhibited higher glucose levels following intraperitoneal glucose challenge (control 12.7mmol/L ±0.6 vs. Gcgr β-cell-/- 15.4mmol/L ±0.0 at 15 min, p=0.002); fasting glycaemia was not different to controls. In vitro, Gcgr β-cell-/- islets showed profound loss of pan-islet [Ca2+]I waves in response to glucose which was only partially rescued in vivo. Diet induced obesity and hyperglycaemia also resulted in a loss of co-ordinated [Ca2+]I waves in transplanted islets. This was reversed with GCG-analogue treatment, independently of weight-loss (n=8).

CONCLUSION

These data provide novel evidence for the role of intra-islet GCGR signalling in sustaining synchronised [Ca2+]I waves and support a possible therapeutic role for glucagonergic agents to restore the insulin secretory capacity lost in T2D.

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.

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.

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.

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.

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.

Unraveling the Molecular Mechanisms of Alcohol-Mediated Skin Permeation Enhancement: Insights from Molecular Dynamics Simulations

The application of alcohols as permeation enhancers in pharmaceutical and cosmetic formulations has attracted considerable attention, owing to their skin permeation-enhancing effect. Nonetheless, the elucidation of the fundamental mechanisms underlying the skin permeation-enhancing effect remains elusive. In this study, molecular dynamics (MD) simulations were employed to investigate the effect of 1,2-propanediol (1,2-PDO), 1,2-butanediol (1,2-BDO), and ethanol (EtOH) on the stratum corneum (SC) model membrane. The results showed that the effect of alcohols on the SC model membrane displayed a concentration-dependent nature. The alcohols can interact with SC lipids and exhibit a remarkable ability to selectively extract free fatty acid (FFA) molecules from the SC model membrane and make the SC looser. Meanwhile, 1,2-BDO and EtOH can penetrate into SC lipid bilayers at higher concentrations, leading to the formation of continuous hydrophilic defects in SC. The FFA extraction and the formation of continuous hydrophilic defects induced ceramide (CER) tail chains to become more disordered and fluid and also weakened the hydrogen bonding (H-bonding) network among SC lipids. Both the FFA extraction and the continuous hydrophilic defect formation endowed alcohols with the permeation-enhancing effect. The constrained simulations revealed that the free energy barriers decreased for the permeation of the hydrophilic model molecule (COL) across the SC model membranes containing alcohols, particularly for 1,2-BDO and EtOH. The possible permeation-enhancing mechanisms of alcohols were proposed correspondingly. This work not only provided a deep understanding of the transdermal permeation-enhancing behavior of alcohols at the molecular level but also provided necessary reference information for designing effective transdermal drug delivery systems in applications.

[Epidemiological characteristics of a local cluster epidemic caused by the BA.2 evolutionary branch of Omicron variant]

Descriptive epidemiological methods were used to analyze the epidemiological characteristics of the local cluster of COVID-19 in the logistic park of Yuhang District in Hangzhou in March 2022. The cluster epidemic was detected by a case who actively visited the fever clinic. The epidemic lasted for 8 days, and a total of 58 cases (53 workers, 2 students, 1 farmer, 1 teacher and 1 unemployed) were found, including 40 males and 18 females. The age was (33.29±12.22) years. There cases were mainly in Yuhang District (48 cases, 82.77%) and Shangcheng District (7 cases, 12.07%) of Hangzhou. The real-time regeneration number peaked at 2.31 on March 10th and decreased to 0.37 on March 15th. The sequencing result of the indicated case was 100% homologous with the sequence uploaded from South Korea on March 4th, 2022.

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.

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.

WISP1 Is Involved in the Pathogenesis of Kashin-Beck Disease via the Autophagy Pathway

OBJECTIVE

Kashin-Beck disease (KBD) is a kind of endemic and chronic osteochondropathy in China. This study aims to explore the functional relevance and potential mechanism of Wnt-inducible signaling pathway protein 1 (WISP1) in the pathogenesis of KBD.

DESIGN

KBD and control cartilage specimens were collected for tissue section observation and primary chondrocyte culture. Firstly, the morphological and histopathological observations were made under a light and electron microscope. Then, the expression levels of WISP1 as well as molecular markers related to the autophagy pathway and extracellular matrix (ECM) synthesis were detected in KBD and control chondrocytes by qRT-PCR, Western blot, and immunohistochemistry. Furthermore, the lentiviral transfection technique was applied to make a WISP1 knockdown cell model based on KBD chondrocytes. In vitro intervention experiments were conducted on the C28/I2 human chondrocyte cell line using human recombinant WISP1 (rWISP1).

RESULTS

The results showed that the autolysosome appeared in the KBD chondrocytes. The expression of WISP1 was significantly higher in KBD chondrocytes. Additionally, T-2 toxin, a risk factor for KBD onset, could up-regulate the expression of WISP1 in C28/I2. The autophagy markers ATG4C and LC3II were upregulated after the low-concentration treatment of T-2 toxin and downregulated after the high-concentration treatment. After knocking down WISP1 expression in KBD chondrocytes, MAP1LC3B decreased while ATG4C and COL2A1 increased. Moreover, the rWISP1 protein treatment in C28/I2 chondrocytes could upregulate the expression of ATG4C and LC3II at the beginning and downregulate them then.

CONCLUSIONS

Our study suggested that WISP1 might play a role in the pathogenesis of KBD through autophagy.

Epigenetic analysis suggests aberrant cerebellum brain aging in old-aged adults with autism spectrum disorder and schizophrenia

The aberrant aging hypothesis of schizophrenia (SCZ) and autism spectrum disorder (ASD) has been proposed, and the DNA methylation (DNAm) clock, which is a cumulative evaluation of DNAm levels at age-related CpGs, could serve as a biological aging indicator. This study evaluated epigenetic brain aging of ASD and SCZ using Horvath's epigenetic clock, based on two public genome-wide DNA methylation datasets of post-mortem brain samples (NASD = 222; NSCZ = 142). Total subjects were further divided into subgroups by gender and age. The epigenetic age acceleration (AgeAccel) for each sample was calculated as the residual value resulting from the regression model and compared between groups. Results showed DNAm age has a strong correlation with chronological age in both datasets across multiple brain regions (P < 0.05). When divided into equally sized age groups, the AgeAccel of the cerebellum (CB) region from people over 45 years of age was greater compared to the control sample (AgeAccel of ASD vs control: 5.069 vs -6.249; P < 0.001). And a decelerated epigenetic aging process was observed in the CB region of individuals with SCZ aged 50-70 years (AgeAccel of SCZ vs control: -3.171 vs 2.418; P < 0.05). However, our results showed no significant difference in AgeAccel between ASD and control groups, and between SCZ and control groups in the total and gender-specific groups (P > 0.05). This study's results revealed some evidence for aberrant epigenetic CB brain aging in old-aged patients with ASD and SCZ, indicating a different pattern of CB aging in older adults with these two diseases. However, further studies of larger ASD and SCZ cohorts are necessary to make definitive conclusions on this observation.

Evaluating the interaction between 3'aQTL and alcohol consumption/smoking on anxiety and depression: 3'aQTL-by-environment interaction study in UK Biobank cohort

BACKGROUND

Smoking and alcohol consumption were associated with the development of depression and anxiety. 3'UTR APA quantitative trait loci (3'aQTLs) have been associated with multiple health states and conditions. Our aim is to evaluate the interactive effects of 3'aQTLs-alcohol consumption/tobacco smoking on the risk of anxiety and depression.

METHODS

The 3'aQTL data of 13 brain regions were extracted from the large-scale 3'aQTL atlas. The phenotype data (frequency of cigarette smoking and alcohol drinking, anxiety score, self-reported anxiety, depression score and self-reported depression) of 90,399-103,011 adults aged 40-69 years living in the UK and contributing to the UK Biobank during 2006-2010, were obtained from the UK Biobank cohort. The frequency of cigarette smoking and alcohol drinking of each subject were defined by the amount of smoking and alcohol drinking of self-reported, respectively. The continuous alcohol consumption/smoking terms were further categorized in tertiles. 3'aQTL-by-environmental interaction analysis was then performed to evaluate the associations of gene-smoking/alcohol consumption interactions with anxiety and depression using generalized linear model (GLM) of PLINK 2.0 with an additive mode of inheritance. Furthermore, GLM was also used to explore the relationship between alcohol consumption/smoking with hazard of anxiety/depression stratified by allele for the significant genotyped SNPs that modified the alcohol consumption/smoking-anxiety/depression association.

RESULTS

The interaction analysis identified several candidate 3'aQTLs-alcohol consumption interactions, such as rs7602638 located in PPP3R1 (β = 0.08, P = 6.50 × 10-6) for anxiety score; rs10925518 located in RYR2 (OR = 0.95, P = 3.06 × 10-5) for self-reported depression. Interestingly, we also observed that the interactions between TMOD1 (β = 0.18, P = 3.30 × 10-8 for anxiety score; β = 0.17, P = 1.42 × 10-6 for depression score), ZNF407 (β = 0.17, P = 2.11 × 10-6 for anxiety score; β = 0.15, P = 4.26 × 10-5 for depression score) and alcohol consumption was not only associated with anxiety, but related to depression. Besides, we found that relationship between alcohol consumption and hazard of anxiety/depression was significantly different for different SNPs genotypes, such as rs34505550 in TMOD1 (AA: OR = 1.03, P = 1.79 × 10-6; AG: OR = 1.00, P = 0.94; GG: OR = 1.00, P = 0.21) for self-reported anxiety.

LIMITATIONS

The identified 3'aQTLs-alcohol consumption/smoking interactions were associated with depression and anxiety, and its potential biological mechanisms need to be further revealed.

CONCLUSIONS

Our study identified important interactions between candidate 3'aQTL and alcohol consumption/smoking on depression and anxiety, and found that the 3'aQTL may modify the associations between consumption/smoking with depression and anxiety. These findings may help to further explore the pathogenesis of depression and anxiety.

Dissecting the Association between Gut Microbiota and Brain Structure Change Rate: A Two-Sample Bidirectional Mendelian Randomization Study

The connection between the gut microbiota and brain structure changes is still unclear. We conducted a Mendelian randomization (MR) study to examine the bidirectional causality between the gut microbiota (211 taxa, including 131 genera, 35 families, 20 orders, 16 classes and 9 phyla; N = 18,340 individuals) and age-independent/dependent longitudinal changes in brain structure across the lifespan (N = 15,640 individuals aged 4~99 years). We identified causal associations between the gut microbiota and age-independent/dependent longitudinal changes in brain structure, such as family Peptostreptococcaceae with age-independent longitudinal changes of cortical gray matter (GM) volume and genus Faecalibacterium with age-independent average cortical thickness and cortical GM volume. Taking age-independent longitudinal changes in brain structure across the lifespan as exposures, there were causal relationships between the surface area and genus Lachnospiraceae. Our findings may serve as fundamentals for further research on the genetic mechanisms and biological treatment of complex traits and diseases associated with the gut microbiota and the brain structure change rate.

[Tea consumption and cancer: a Mendelian randomization study]

Objective: A Mendelian randomization (MR) analysis was performed to assess the relationship between tea consumption and cancer. Methods: There were 100 639 participants with the information of gene sequencing of whole genome in the China Kadoorie Biobank. After excluding those with cancer at baseline survey, a total of 100 218 participants were included in this study. The baseline information about tea consumption were analyzed, including daily tea consumption or not, cups of daily tea consumption, and grams of daily tea consumption. We used the two-stage least square method to evaluate the associations between three tea consumption variables and incidence of cancer and some subtypes, including stomach cancer, liver and intrahepatic bile ducts cancer, colorectal cancer, tracheobronchial and lung cancer, and female breast cancer. Multivariable MR and analysis only among nondrinkers were used to control the impact of alcohol consumption. Sensitivity analyses were also performed, including inverse variance weighting, weighted median, and MR-Egger. Results: We used 54, 42, and 28 SNPs to construct non-weighted genetic risk scores as instrumental variables for daily tea consumption or not, cups of daily tea consumption, and grams of daily tea consumption, respectively. During an average of (11.4±3.0) years of follow-up, 6 886 cases of cancer were recorded. After adjusting for age, age², sex, region, array type, and the first 12 genetic principal components, there were no significant associations of three tea consumption variables with the incidence of cancer and cancer subtypes. Compared with non-daily tea drinkers, the HR (95%CI) of daily tea drinkers for cancer and some subtypes, including stomach cancer, liver and intrahepatic bile ducts cancer, colorectal cancer, tracheobronchial and lung cancer, and female breast cancer, are respectively 0.99 (0.78-1.26), 1.17 (0.58-2.36), 0.86 (0.40-1.84), 0.85 (0.42-1.73), 1.39 (0.85-2.26) and 0.63 (0.28-1.38). After controlling the impact of alcohol consumption and performing multiple sensitivity analyses, the results were similar. Conclusion: There is no causal relationship between tea consumption and risk of cancer in population in China.

Mitochondria-wide association study observed significant interactions of mitochondrial respiratory and the inflammatory in the development of anxiety and depression

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

The effects of size and surface functionalization of polystyrene nanoplastics on stratum corneum model membranes: An experimental and computational study

Nanoplastics are mainly generated from the decomposition of plastic waste and artificial production and have attracted much attention due to their wide distribution in the environment and the potential risk for humans. As the largest organ of the human body, the skin is inevitably in contact with nanoplastics. Stratum corneum is the first barrier when the skin is exposed to nanoplastics. However, little is known about the interactions between nanoplastics and stratum corneum. Here, the effects of particle size and surface functionalization (amino-modified and carboxy-modified) of polystyrene nanoplastics on the stratum corneum models were studied by Langmuir monolayer and molecular dynamics simulations. An equimolar mixture of ceramide/cholesterol/free fatty acid was used to mimic stratum corneum intercellular lipids. The Langmuir monolayer studies demonstrated that the larger size and surface functionalization of polystyrene nanoplastics significantly reduced the stability of stratum corneum lipid monolayer in a concentration-dependent fashion. Simulation results elucidated that functionalized polystyrene oligomers had a stronger interaction with lipid components of the stratum corneum model membrane. The cell experiments also indicated that functionalized polystyrene nanoplastics, especially for amino-modified polystyrene nanoplastics, had significant cytotoxicity on normal human dermal fibroblast cells. Our results provide fundamental information and the basis for a deeper understanding of the health risks of nanoplastics to humans.

Exploratory factor analysis of shared and specific genetic associations in depression and anxiety

BACKGROUND

Previous genetic studies of anxiety and depression were mostly based on independent phenotypes. This study aims to investigate the shared and specific genetic structure between anxiety and depression.

METHOD

To identify the underlying factors of Generalized Anxiety Disorder-7 (GAD-7), Patient Health Questionnaire-9 (PHQ-9), and their combined scale (joint scale), we employed exploratory factor analysis (EFA) using the eigenvalue of parallel analysis. Subsequently, we conducted a genome-wide association study (GWAS) for these factors. In addition, we utilized LD Score Regression (LDSC) to determine the genetic correlations between the identified factors and four common mental disorders, three sleep phenotypes, and other traits that have been previously linked to anxiety and depression.

RESULTS

The EFA uncovered two factors for the GAD-7 scale, namely nervousness and disturbance, two factors for the PHQ-9 scale, namely negative affect and sleep/appetite disturbance, and four factors for the joint scale, specifically nervousness, anhedonia, sleep/appetite disturbance, and fidget. We identified two genome-wide significant genomic loci, with overlap across GAD-7 factor 1 and joint scale factor 1: rs148579586 (P_GAD-7 = 1.365 × 10^-09, P_{Joint scale} = 1.434 × 1^0-09) and rs201074060 (P_GAD-7 = 3.672 × 10^-09, P_{Joint scale} = 3.824 × 10^-09). Genetic correlations in factors ranged from 0.722 to 1.000 (all p < 1.786 × 10^-3) with 27 of 28 correlations being significantly smaller than one. The genetic correlations with external phenotypes showed small variation across the eight factors.

CONCLUSION

Unidimensional structures can provide more precise scores, which can aid in identifying the shared and specific genetic associations between anxiety and depression. This is a crucial step in characterizing the genetic structure of these conditions and their co-occurrence.

Evaluating the interactive effects of dietary habits and human gut microbiome on the risks of depression and anxiety

BACKGROUND

Gut microbiome and dietary patterns have been suggested to be associated with depression/anxiety. However, limited effort has been made to explore the effects of possible interactions between diet and microbiome on the risks of depression and anxiety.

METHODS

Using the latest genome-wide association studies findings in gut microbiome and dietary habits, polygenic risk scores (PRSs) analysis of gut microbiome and dietary habits was conducted in the UK Biobank cohort. Logistic/linear regression models were applied for evaluating the associations for gut microbiome-PRS, dietary habits-PRS, and their interactions with depression/anxiety status and Patient Health Questionnaire (PHQ-9)/Generalized Anxiety Disorder-7 (GAD-7) score by R software.

RESULTS

We observed 51 common diet-gut microbiome interactions shared by both PHQ score and depression status, such as overall beef intake × genus Sporobacter [hurdle binary (HB)] (PPHQ = 7.88 × 10-4, Pdepression status = 5.86 × 10-4); carbohydrate × genus Lactococcus (HB) (PPHQ = 0.0295, Pdepression status = 0.0150). We detected 41 common diet-gut microbiome interactions shared by GAD score and anxiety status, such as sugar × genus Parasutterella (rank normal transformed) (PGAD = 5.15 × 10-3, Panxiety status = 0.0347); tablespoons of raw vegetables per day × family Coriobacteriaceae (HB) (PGAD = 6.02 × 10-4, Panxiety status = 0.0345). Some common significant interactions shared by depression and anxiety were identified, such as overall beef intake × genus Sporobacter (HB).

CONCLUSIONS

Our study results expanded our understanding of how to comprehensively consider the relationships for dietary habits-gut microbiome interactions with depression and anxiety.

[Clinicopathologic characteristics and prognostic analysis of testicular diffuse large B-cell lymphoma]

Objective: To analyze the clinicopathologic characteristics and prognosis of testicular diffuse large B-cell lymphoma (DLBCL) . Methods: A retrospective analysis was performed on 68 patients with testicular DLBCL admitted to Ruijin Hospital affiliated to Shanghai Jiao Tong University School of Medicine from October 2001 to April 2020. The gene mutation profile was evaluated by targeted sequencing (55 lymphoma-related genes) , and prognostic factors were analyzed. Results: A total of 68 patients were included, of whom 45 (66.2% ) had primary testicular DLBCL and 23 (33.8% ) had secondary testicular DLBCL. The proportion of secondary testicular DLBCL patients with Ann Arbor stage Ⅲ-Ⅳ (P<0.001) , elevated LDH (P<0.001) , ECOG score ≥ 2 points (P=0.005) , and IPI score 3-5 points (P<0.001) is higher than that of primary testicular DLBCL patients. Sixty-two (91% ) patients received rituximab in combination with cyclophosphamide, adriamycin, vincristine, and prednisone (R-CHOP) -based first-line regimen, whereas 54 cases (79% ) underwent orchiectomy prior to chemotherapy. Patients with secondary testicular DLBCL had a lower estimated 5-year progression-free survival (PFS) rate (16.5% vs 68.1% , P<0.001) and 5-year overall survival (OS) rate (63.4% vs 74.9% , P=0.008) than those with primary testicular DLBCL, and their complete remission rate (57% vs 91% , P=0.003) was also lower than that of primary testicular DLBCL. The ECOG scores of ≥2 (PFS: P=0.018; OS: P<0.001) , Ann Arbor stages Ⅲ-Ⅳ (PFS: P<0.001; OS: P=0.018) , increased LDH levels (PFS: P=0.015; OS: P=0.006) , and multiple extra-nodal involvements (PFS: P<0.001; OS: P=0.013) were poor prognostic factors in testicular DLBCL. Targeted sequencing data in 20 patients with testicular DLBCL showed that the mutation frequencies of ≥20% were PIM1 (12 cases, 60% ) , MYD88 (11 cases, 55% ) , CD79B (9 cases, 45% ) , CREBBP (5 cases, 25% ) , KMT2D (5 cases, 25% ) , ATM (4 cases, 20% ) , and BTG2 (4 cases, 20% ) . The frequency of mutations in KMT2D in patients with secondary testicular DLBCL was higher than that in patients with primary testicular DLBCL (66.7% vs 7.1% , P=0.014) and was associated with a lower 5-year PFS rate in patients with testicular DLBCL (P=0.019) . Conclusion: Patients with secondary testicular DLBCL had worse PFS and OS than those with primary testicular DLBCL. The ECOG scores of ≥2, Ann Arbor stages Ⅲ-Ⅳ, increased LDH levels, and multiple extra-nodal involvements were poor prognostic factors in testicular DLBCL. PIM1, MYD88, CD79B, CREBBP, KMT2D, ATM, and BTG2 were commonly mutated genes in testicular DLBCL, and the prognosis of patients with KMT2D mutations was poor.

Association between birth by caesarian section and anxiety, self-harm: a gene-environment interaction study using UK Biobank data

BACKGROUND

Limited efforts have been paid to explore the underlying genetic mechanisms of birth by caesarian section (CS) affecting the risks of adult anxiety and self-harm.

METHODS

Using UK Biobank cohort, the logistic regression model was first applied to evaluate the associations of adult anxiety and self-harm with birth by CS. Using birth by CS as exposure variables, genome-wide by environment interaction study (GWEIS) was then applied by PLINK2.0 to identify associated genes interacting with birth by CS for anxiety and self-harm.

RESULTS

In the observational study, significant associations were observed between birth by CS and anxiety (odds ratio (OR) = 1.24; 95% confidence interval (CI), 1.12-1.38; P = 4.86 × 10^- 5), and self-harm (OR = 1.12; 95% CI, 1.01-1.24; P = 2.90 × 10^- 2). GWEIS revealed multiple suggestive genes interacted with birth by CS for anxiety, such as DKK2 (rs13137764, P = 1.24 × 10^- 9, adjusted P = 2.68 × 10^- 7) and ATXN1 (rs62389045, P = 4.38 × 10^- 8, adjusted P = 3.55 × 10^- 6). For self-harm, significant gene-environment interactions of birth by CS on self-harm were detected, such as ALDH1A2 (rs77828167, P = 1.62 × 10^- 8; rs116899929, P = 1.92 × 10^- 8) and DAB1 (rs116124269, P = 3.20 × 10^- 8; rs191070006, P = 3.63 × 10^- 8).

CONCLUSIONS

Our results suggested that birth by CS was associated with the risk of adult anxiety and self-harm. We also discovered some genes interacted with birth by CS might influence the risk of anxiety and self-harm, which may provide novel clues for the pathogenesis of those mental disorders.

The effects of adsorbed benzo(a)pyrene on dynamic behavior of polystyrene nanoplastics through phospholipid membrane: A molecular simulation study

Nanoplastics (NPs) are mainly generated from the decomposition of plastic waste and industrial production, which have attracted much attention due to the potential risk for humans. The ability of NPs to penetrate different biological barriers has been proved, but the understanding of molecular details is very limited, especially for organic pollutant-NP combinations. Here, we investigated the uptake process of polystyrene NPs (PSNPs) combined with benzo(a)pyrene (BAP) molecules by dipalmitoylphosphatidylcholine (DPPC) bilayers by molecular dynamics (MD) simulations. The results showed that the PSNPs can adsorb and accumulate BAP molecules in water phase and then carried BAP molecules to enter DPPC bilayers. At the same time, the adsorbed BAP promoted the penetration of PSNPs into DPPC bilayers effectively by hydrophobic effect. The process of BAP-PSNP combinations penetrating into DPPC bilayers can be summarized into four steps including adhesion on the DPPC bilayer surface, uptake by the DPPC bilayer, BAP molecules detached from the PSNPs, and the PSNPs depolymerized in the bilayer interior. Furthermore, the amount of adsorbed BAP on PSNPs affected the properties of DPPC bilayers directly, especially the fluidity of DPPC bilayers that determine the physiologic function. Obviously, the combined effect of PSNPs and BAP enhanced the cytotoxicity. This work not only presented a vivid transmembrane process of BAP-PSNP combinations and revealed the nature of the effects of adsorbed benzo(a)pyrene on the dynamic behavior of polystyrene nanoplastics through phospholipid membrane, but also provide some necessary information of the potential damage for organic pollutant-nanoplastic combinations on human health at a molecular level.

Evidence of Antineutrinos from Distant Reactors Using Pure Water at SNO

The SNO+ Collaboration reports the first evidence of reactor antineutrinos in a Cherenkov detector. The nearest nuclear reactors are located 240 km away in Ontario, Canada. This analysis uses events with energies lower than in any previous analysis with a large water Cherenkov detector. Two analytical methods are used to distinguish reactor antineutrinos from background events in 190 days of data and yield consistent evidence for antineutrinos with a combined significance of 3.5σ.

[Effect of postoperative adjuvant chemotherapy on prognosis of patients with intrahepatic cholangiocarcinoma:a multicenter retrospective study]

Objectives: To examine the influence of adjuvant chemotherapy after radical resection on the survival of patients with intrahepatic cholangiocarcinoma(ICC) and to identify patients who may benefit from it. Methods: The clinical and pathological data of 654 patients with ICC diagnosed by postoperative pathology from December 2011 to December 2017 at 13 hospitals in China were collected retrospectively. According to the inclusion and exclusion criteria,455 patients were included in this study,including 69 patients (15.2%) who received adjuvant chemotherapy and 386 patients (84.8%) who did not receive adjuvant chemotherapy. There were 278 males and 177 females,with age of 59 (16) years (M(IQR))(range:23 to 88 years). Propensity score matching (PSM) method was used to balance the difference between adjuvant chemotherapy group and non-adjuvant chemotherapy group. Kaplan-Meier method was used to plot the survival curve,the Log-rank test was used to compare the difference of overall survival(OS) and recurrence free survival(RFS)between the two groups. Univariate analysis was used to determine prognostic factors for OS. Multivariate Cox proportional hazards models were then performed for prognostic factors with P<0.10 to identify potential independent risk factors. The study population were stratified by included study variables and the AJCC staging system,and a subgroup analysis was performed using the Kaplan-Meier method to explore the potential benefit subgroup population of adjuvant chemotherapy. Results: After 1∶1 PSM matching,69 patients were obtained in each group. There was no significant difference in baseline data between the two groups (all P>0.05). After PSM,Cox multivariate analysis showed that lymph node metastasis (HR=3.06,95%CI:1.52 to 6.16,P=0.039),width of resection margin (HR=0.56,95%CI:0.32 to 0.99,P=0.044) and adjuvant chemotherapy (HR=0.51,95%CI:0.29 to 0.91,P=0.022) were independent prognostic factors for OS. Kaplan-Meier analysis showed that the median OS time of adjuvant chemotherapy group was significantly longer than that of non-adjuvant chemotherapy group (P<0.05). There was no significant difference in RFS time between the adjuvant chemotherapy group and the non-adjuvant chemotherapy group (P>0.05). Subgroup analysis showed that,the OS of female patients,without HBV infection,carcinoembryonic antigen<9.6 μg/L,CA19-9≥200 U/ml,intraoperative bleeding<400 ml,tumor diameter>5 cm,microvascular invasion negative,without lymph node metastasis,and AJCC stage Ⅲ patients could benefit from adjuvant chemotherapy (all P<0.05). Conclusion: Adjuvant chemotherapy can prolong the OS of patients with ICC after radical resection,and patients with tumor diameter>5 cm,without lymph node metastasis,AJCC stage Ⅲ,and microvascular invasion negative are more likely to benefit from adjuvant chemotherapy.

Integrating genome-wide association study with regulatory SNP annotations identified novel candidate genes for osteoporosis

Osteoporosis (OP) is a metabolic bone disease, characterized by a decrease in bone mineral density (BMD). However, the research of regulatory variants has been limited for BMD. In this study, we aimed to explore novel regulatory genetic variants associated with BMD. We conducted an integrative analysis of BMD genome-wide association study (GWAS) and regulatory single nucleotide polymorphism (rSNP) annotation information. Firstly, the discovery GWAS dataset and replication GWAS dataset were integrated with rSNP annotation database to obtain BMD associated SNP regulatory elements and SNP regulatory element-target gene (E-G) pairs, respectively. Then, the common genes were further subjected to HumanNet v2 to explore the biological effects. Through discovery and replication integrative analysis for BMD GWAS and rSNP annotation database, we identified 36 common BMD-associated genes for BMD irrespective of regulatory elements, such as FAM3C (p_{discovery GWAS} = 1.21 × 10^-25, p_{replication GWAS} = 1.80 × 10^-12), CCDC170 (p_{discovery GWAS} = 1.23 × 10^-11, p_{replication GWAS} = 3.22 × 10^-9), and SOX6 (p_{discovery GWAS} = 4.41 × 10^-15, p_{replication GWAS} = 6.57 × 10^-14). Then, for the 36 common target genes, multiple gene ontology (GO) terms were detected for BMD such as positive regulation of cartilage development (p = 9.27 × 10^-3) and positive regulation of chondrocyte differentiation (p = 9.27 × 10^-3). We explored the potential roles of rSNP in the genetic mechanisms of BMD and identified multiple candidate genes. Our study results support the implication of regulatory genetic variants in the development of OP.

Assessing the interaction effects of brain structure longitudinal changes and life environmental factors on depression and anxiety

Disrupted brain structures and several life environmental factors have been shown to influence depression and anxiety, but their interactions with anxiety and depression remain elusive. Genome-wide association study datasets of 15 brain structure longitudinal changes (N = 15,640) were obtained from the published study. Genotype and phenotype-related data of depression, anxiety, and life environmental factors (including smoking, alcohol drinking, coffee intake, maternal smoking, physical activity, vitamin D, insomnia, sleep duration, and family satisfaction) were collected from UK Biobank. We calculated the polygenic risk scores (PRS) of 15 brain structure changes and then conducted linear regression analyses to explore the interactions of brain structure changes and life environmental factors on depression and anxiety using 15 brain structure change-related PRS, life environmental factors and interactions of them as instrumental variables, and depression score or anxiety score as outcomes. Sex stratification in all analyses was performed to reveal sex-specific differences in the interactions. We found 14 shared interactions related to both depression and anxiety in total sample, such as alcohol drinking × cerebellum white matter 3 (WM; beta = -.003, p = .018 for depression; beta = -003, p = .008 for anxiety) and maternal smoking × nucleus accumbens 2 (beta = .088, p = .002 for depression; beta = .070, p = .008 for anxiety). We also observed sex-specific differences in the interactions, for instance, alcohol drinking × cerebellum WM 3 was negatively associated with depression and anxiety in males (beta = -.004, p = .020 for depression; beta = -.005, p = .002 for anxiety). Our study results reveal the important interactions between brain structure changes and several life environmental factors on depression and anxiety, which may help to explore the pathogenesis of depression and anxiety.

Treatment Resistance in Schizophrenia Is Associated with Attention Deficit/Hyperactivity Disorder and Gut Microbiota: A Genetic Correlation and Mendelian Randomization Study

INTRODUCTION

Observational studies highlight associations of common diseases with individual schizophrenia symptoms. However, it is unclear whether these diseases are associated with individual treatment-resistant schizophrenia (TRS). We aimed to explore the genetic associations between common immune diseases, metabolic diseases, psychiatric disorders, gut microbiota and TRS.

METHODS

Genome-wide association study (GWAS) summary data of European participants (n = ∼456,327) included TRS, 11 psychiatric disorders, 23 immune and metabolic diseases, body mass index, height, and 211 gut microbiota. In this genetic correlation and two-sample Mendelian randomization (MR) study, linkage disequilibrium score (LDSC) regression was applied to infer genetic correlation estimates. Two-sample MR tested potential causal associations of genetic variants associated with common immune diseases, metabolic diseases, psychiatric disorders, and gut microbiota with TRS.

RESULTS

LDSC revealed candidate associations between attention deficit/hyperactivity disorder (ADHD), schizophrenia, intestinal infectious diseases, obesity and TRS (genetic correlation range, 0.230-0.702; p < 0.05). Two-sample MR analyses suggested that ADHD was positively associated with TRS (estimate [SE] = 0.204 [0.073], p = 0.005), a finding that remained stable across statistical models. Besides, schizophrenia and genus Barnesiella levels were causally associated with TRS but not consistent across MR approaches.

CONCLUSION

This study reports genetic correlations between ADHD, schizophrenia, intestinal infectious diseases, obesity and TRS. The study also found that genus Barnesiella was associated with TRS. These findings may have clinical implications, highlighting the possible strategy for TRS prevention.

Assessing the effect of interaction between lifestyle and longitudinal changes in brain structure on sleep phenotypes

Longitudinal changes in brain structure and lifestyle can affect sleep phenotypes. However, the influence of the interaction between longitudinal changes in brain structure and lifestyle on sleep phenotypes remains unclear. Genome-wide association study dataset of longitudinal changes in brain structure was obtained from published study. Phenotypic data of lifestyles and sleep phenotypes were obtained from UK Biobank cohort. Using genotype data from UK Biobank, we calculated polygenetic risk scores of longitudinal changes in brain structure phenotypes. Linear/logistic regression analysis was conducted to evaluate interactions between longitudinal changes in brain structure and lifestyles on sleep duration, chronotype, insomnia, snoring and daytime dozing. Multiple lifestyle × longitudinal changes in brain structure interactions were detected for 5 sleep phenotypes, such as physical activity×caudate_age2 for daytime dozing (OR = 1.0389, P = 8.84 × 10-3) in total samples, coffee intake×cerebellar white matter volume_age2 for daytime dozing (OR = 0.9652, P = 1.13 × 10-4) in females. Besides, we found 4 overlapping interactions in different sleep phenotypes. We conducted sex stratification analysis and identified one overlapping interaction between female and male. Our results support the moderate effects of interaction between lifestyle and longitudinal changes in brain structure on sleep phenotypes, and deepen our understanding of the pathogenesis of sleep disorders.

[Clinical characteristics and prognosis of primary and secondary diffuse large B-cell lymphoma of the pancreas]

Objective: To analyze the clinical characteristics and prognosis of primary and secondary pancreatic diffuse large B-cell lymphoma (DLBCL) . Methods: Clinical data of patients with pancreatic DLBCL admitted at Shanghai Rui Jin Hospital affiliated with Shanghai Jiao Tong University School of Medicine from April 2003 to June 2020 were analyzed. Gene mutation profiles were evaluated by targeted sequencing (55 lymphoma-related genes). Univariate and multivariate Cox regression models were used to evaluate the prognostic factors of overall survival (OS) and progression-free survival (PFS) . Results: Overall, 80 patients were included; 12 patients had primary pancreatic DLBCL (PPDLBCL), and 68 patients had secondary pancreatic DLBCL (SPDLBCL). Compared with those with PPDLBCL, patients with SPDLBCL had a higher number of affected extranodal sites (P<0.001) and had higher IPI scores (P=0.013). There was no significant difference in the OS (P=0.120) and PFS (P=0.067) between the two groups. Multivariate analysis indicated that IPI intermediate-high/high risk (P=0.025) and double expressor (DE) (P=0.017) were independent adverse prognostic factors of OS in patients with pancreatic DLBCL. IPI intermediate-high/high risk (P=0.021) was an independent adverse prognostic factor of PFS in patients with pancreatic DLBCL. Targeted sequencing of 29 patients showed that the mutation frequency of PIM1, SGK1, BTG2, FAS, MYC, and MYD88 in patients with pancreatic DLBCL were all >20%. PIM1 (P=0.006 for OS, P=0.032 for PFS) and MYD88 (P=0.001 for OS, P=0.017 for PFS) mutations were associated with poor OS and PFS in patients with SPDLBCL. Conclusion: There was no significant difference in the OS and PFS between patients with PPDLBCL and those with SPDLBCL. IPI intermediate-high/high risk and DE were adverse prognostic factors of pancreatic DLBCL. PIM1, SGK1, BTG2, FAS, MYC, and MYD88 were common mutations in pancreatic DLBCL. PIM1 and MYD88 mutations indicated worse prognosis.

Evaluating the Genetic Effects of Gut Microbiota on the Development of Neuroticism and General Happiness: A Polygenic Score Analysis and Interaction Study Using UK Biobank Data

Limited efforts have been invested in exploring the interaction effects between genetic factors and gut microbiota on neuroticism and general happiness. The polygenic risk scores (PRS) of gut microbiota were calculated from individual-level genotype data of the UK Biobank cohort. Linear regression models were then used to assess the associations between individual PRS of gut microbiota and mental traits and interaction analysis was performed by PLINK2.0. KOBAS-i was used to conduct gene ontology (GO) enrichment analysis of the identified genes. We observed suggestive significant associations between neuroticism and PRS for the genus Bifidobacterium (rank-normal transformation, RNT) (beta = -1.10, P = 4.16 × 10-3) and the genus Desulfovibrio (RNT) (beta = 0.54, P = 7.46 × 10-3). PRS for the genus Bifidobacterium (hurdle binary, HB) (beta = 1.99, P = 5.24 × 10-3) and the genus Clostridium (RNT) (beta = 1.26, P = 9.27 × 10-3) were found to be suggestive positively associated with general happiness. Interaction analysis identified several significant genes that interacted with gut microbiota, such as RORA (rs575949009, beta = -45.00, P = 1.82 × 10-9) for neuroticism and ASTN2 (rs36005728, beta = 19.15, P = 3.37 × 10-8) for general happiness. Our study results support the genetic effects of gut microbiota on the development of neuroticism and general happiness.

Evaluating the role of anxiety on the association between irritable bowel syndrome and brain volumes: a mediation analysis in the UK Biobank cohort

There is a strong link between irritable bowel syndrome and brain volumes, yet, to date, research examining the mediators of this association has been little. Based on the phenotypic data of 15 248 participants from the UK Biobank, a two-stage mediation analysis was performed to assess the association among brain volumes, anxiety, and irritable bowel syndrome. In the first stage, we identified the candidate mediating role of anxiety for irritable bowel syndrome associated with brain volumes using regression models. Then, we quantified the magnitude of the mediation effects by evaluating the average causal-mediated effect and proportion of mediation through performing mediation analyses in the R package in the second stage. In the first stage, we identified the partly mediating role of anxiety in the association between irritable bowel syndrome and the volume of thalamus (P _left = 1.16 × 10^-4, P _right = 2.41 × 10^-4), and grey matter (P _left = 3.22 × 10^-2, P _right = 1.18 × 10^-2) in the VIIIa cerebellum. In the second stage, we observed that the proportion of the total effect of irritable bowel syndrome on volume of thalamus mediated by anxiety was 14.3% for the left region (β _{Average causal-mediated effect} = -0.008, P _{Average causal-mediated effect} = 0.004) and 14.6% for the right region (β _{Average causal-mediated effect} = -0.007, P _{Average causal-mediated effect} = 0.006). Anxiety mediated 30.8% for the left region (β _{Average causal-mediated effect} = -0.013, P _{Average causal-mediated effect} = 0.002) and 21.6% for the right region (β _{Average causal-mediated effect} = -0.010, P _{Average causal-mediated effect} x= 0.018) of the total effect of irritable bowel syndrome on the volume of grey matter in the VIIIa cerebellum. Our study revealed the indirect mediating role of anxiety in the association between irritable bowel syndrome and brain volumes, promoting our understanding of the functional mechanisms of irritable bowel syndrome and its related psychosocial factors.

Prognostic Value of Four Preimplantation Malnutrition Estimation Tools in Predicting Heart Failure Hospitalization of the Older Diabetic Patients with Right Ventricular Pacing

OBJECTIVES

The prognostic value of preimplantation nutritional status is not yet known for older diabetic patients that received right ventricular pacing (RVP). The study aimed to investigate the clinical value of the four malnutrition screening tools for the prediction of heart failure hospitalization (HFH) in older diabetic patients that received RVP.

DESIGN

Retrospective observational cohort study.

SETTING AND PARTICIPANTS

This study was conducted between January 2017 and January 2018 at the Fuwai Hospital, Beijing, China, and included older (age ≥ 65 years) diabetic patients that received RVP for the first time Measurements: The Prognostic Nutritional Index (PNI), Geriatric Nutritional Risk Index (GNRI), Naples Prognostic Score (NPS), and the Controlling Nutritional Status (CONUT) score were used to estimate the preimplantation nutritional status of the patients. Univariate and multivariate Cox proportional hazard regression analyses were performed to investigate the association between preimplantation malnutrition and HFH.

RESULTS

Overall, 231 older diabetic patients receiving RVP were included. The median follow-up period after RVP was 53 months. HFH was reported for 19.9% of the included patients. Our results showed preimplantation malnutrition for 18.2%, 15.2%, 86.6% and 66.2% of the included patients based on the PNI, GNRI, NPS, and CONUT score, respectively. The cumulative rate of HFH during follow-up period was significantly higher for patients in the preimplantation malnutrition group based on the PNI (log-rank = 13.0, P = 0.001), GNRI (log-rank = 8.5, P = 0.01), and NPS (log-rank = 15.7, P < 0.001) compared to the normal nutrition group, but was not statistically significant for those in the preimplantation malnutrition group based on the CONUT score (log-rank = 2.7, P = 0.3). As continuous variables, all the nutritional indices showed significant correlation with HFH (all P < 0.05). However, multivariate analysis showed that only GNRI was independently associated with HFH (HR = 0.97, 95% CI: 0.937-0.997, P = 0.032). As categorical variables, PNI, GNRI, and NPS showed significant correlation with HFH. After adjustment of confounding factors, moderate-to-severe degree of malnutrition was an independent predictor of HFH based on the PNI (HR = 4.66, 95% CI: 1.03-21.00, P = 0.045) and GNRI (HR = 3.02, 95% CI: 1.02-9.00, P = 0.047).

CONCLUSION

Preimplantation malnutrition was highly prevalent in older diabetic patients that received RVP. The malnutrition prediction tools, PNI and GNRI, showed significant prognostic value in accurately predicting HFH in older diabetic patients with RVP.

Assessing the joint effects of mitochondrial function and human behavior on the risks of anxiety and depression

BACKGROUND

Psychiatric disorders have great health hazards and the exact pathogeny remains elusive now. We aim to explore the potential interaction effects of mitochondrial function and human behavior on the risks of anxiety and depression.

METHODS

The genome-wide association study (GWAS) data of mitochondrial function (N = 383,476-982,072) were obtained from published studies. Individual level genotype and phenotype data of anxiety, depression and behavioral factors (including drinking, smoking and physical activity) were all from the UK Biobank (N = 84,805-85,164). We first calculated the polygenic risk scores (PRS) of mitochondrial function as the instrumental variables, and then constructed linear regression analyses to systematically explore the potential interaction effects of mitochondrial function and human behavior on anxiety and depression.

RESULTS

In total samples, we observed mitochondrial heteroplasmy (MtHz) vs. Drinking (P_GAD-7 = 6.49 × 10^-3; P_PHQ-9 = 1.89 × 10^-3) was positively associated with both anxiety and depression. In males, MtHz vs. Drinking (P_Male = 3.46 × 10^-5) was positively correlated with depression. In females, blood mitochondrial DNA copy number (mtDNA-CN) vs. Drinking (P_Female = 8.63 × 10^-3) was negatively related to anxiety. Furthermore, we identified additional 6 suggestive interaction effects (P < 0.05) for anxiety and depression.

LIMITATIONS

Considering all subjects were from UK Biobank, it should be careful to extrapolate our findings to other populations with different genetic background.

CONCLUSIONS

Our results suggest the significant impacts of mitochondrial function and human behavior interactions on the development of anxiety and depression, providing new clues for clarifying the pathogenesis of anxiety and depression.

Association between gut microbiota and longevity: a genetic correlation and mendelian randomization study

BACKGROUND

Longevity is one of the most complex phenotypes, and its genetic basis remains unclear. This study aimed to explore the genetic correlation and potential causal association between gut microbiota and longevity.

RESULTS

Linkage disequilibrium score (LDSC) regression analysis and a bi-directional two-sample Mendelian Randomization (MR) analysis were performed to analyze gut microbiota and longevity-related traits. LDSC analysis detected four candidate genetic correlations, including Veillonella (genetic correlation = 0.5578, P = 4.67 × 10^- 2) and Roseburia (genetic correlation = 0.4491, P = 2.67 × 10^- 2) for longevity, Collinsella (genetic correlation = 0.3144, P = 4.07 × 10^- 2) for parental lifespan and Sporobacter (genetic correlation = 0.2092, P = 3.53 × 10^- 2) for healthspan. Further MR analysis observed suggestive causation between Collinsella and parental longevity (father's age at death) (weighted median: b = 1.79 × 10^- 3, P = 3.52 × 10^- 2). Reverse MR analysis also detected several causal effects of longevity-related traits on gut microbiota, such as longevity and Sporobacter (IVW: b = 7.02 × 10^- 1, P = 4.21 × 10^- 25). Statistical insignificance of the heterogeneity test and pleiotropy test supported the validity of the MR study.

CONCLUSION

Our study found evidence that gut microbiota is causally associated with longevity, or vice versa, providing novel clues for understanding the roles of gut microbiota in aging development.

The genetic structure of pain in depression patients: A genome-wide association study and proteome-wide association study

BACKGROUND

Comparing with the general population, the pain in depression patients has more complex biological mechanism. We aim to explore the etiological mechanism of pain in depression patients from the perspective of genetics.

METHODS

Utilizing the UK Biobank samples with self-reported depression status or PHQ score ≥10, we conducted genome-wide association studies (GWAS) of seven pain traits (N = 1,133-58,349). Additionally, we used FUSION pipeline to perform proteome-wide association study (PWAS) and transcriptome-wide association study (TWAS) by integrating GWAS summary data with two different proteome reference weights (ROS/MAP and Banner) and Rnaseq gene expression reference weights, respectively.

RESULTS

GWAS identified 3 significant genes associated with different pain traits in depression patients, including TRIOBP (P_GWAS = 4.48 × 10^-8) for stomach or abdominal pain, SLC9A9(P_GWAS = 2.77 × 10^-8) for multisite chronic pain (MCP) and ADGRF1 (P_GWAS = 1.51 × 10^-8) for neck or shoulder pain. In addition, PWAS and TWAS analysis also identified multiple candidate genes associated with different pain traits in depression patients, such as TPRG1L (P_PWAS-Banner = 3.38 × 10^-2) and SIRPA (P_PWAS-Banner = 3.65 × 10^-2) for MCP, etc. Notably, when comparing the results of PWAS and TWAS analysis, we found overlapping candidate genes in these pain traits, such as GSTM3 (P_PWAS-_Banner = 3.38 × 10^-2, P_TWAS = 6.92 × 10^-3) in the stomach or abdominal pain phenotype, ATG7 (P_PWAS-_Rosmap = 3.15 × 10^-2, P_TWAS = 2.98 × 10^-2) in the MCP, etc. CONCLUSIONS: We identified multiple novel candidate genes for pain traits in depression patients from different perspectives of genetics, which provided novel clues for understanding the genetic mechanisms underlying the pain in depression patients.

Integrated analysis of proteome-wide and transcriptome-wide association study identified novel genes and chemicals for vertigo

Vertigo is a leading symptom of various peripheral and central vestibular disorders. Although genome-wide association studies (GWAS) have identified multiple risk variants for vertigo, how these risk variants contribute to the risk of vertigo remains unknown. Discovery proteome-wide association study (PWAS) was first performed by integrating the protein quantitative trait loci from dorsolateral prefrontal cortex (DLPFC) in Banner Sun Health Research Institute dataset (N = 152) and GWAS summary of vertigo (N = 942,613), followed by replication PWAS using the protein quantitative trait loci from the DLPFC in Religious Orders Study or the Rush Memory and Aging Project dataset (N = 376). Transcriptome-wide association studies (TWASs) were then performed by integrating the same GWAS datasets of vertigo (N = 942,613) with mRNA expression reference from human fetal brain, and DLPFC. Chemical-related gene set enrichment analysis (GSEA) and Gene ontology/Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were finally conducted to further reveal the pathogenesis of vertigo. Permutation-based empirical P values were calculated in PWAS, TWAS and GSEA. By integrating the GWAS of vertigo and two independent brain proteomes from human DLPFC, 3 genes were identified to genetically regulate protein abundance levels in vertigo, and were not previously implicated by GWAS, including MTERFD2 (PBanner = 0.045, PROSMAP = 0.031), MGST1 (PBanner = 0.014, PROSMAP = 0.018), and RAB3B (PBanner = 0.045, PROSMAP = 0.035). Compared with TWAS results, we identified overlapping genes RAB3B (PTWAS = 0.017) and MTERFD2 (PTWAS = 0.003) that showed significant associations with vertigo at both proteome-wide and transcriptome-wide levels. Chemical-related GSEA identified multiple chemicals that might be associated with vertigo, such as nickel (P = 0.007), glycidamide (P = 0.005), and proanthocyanidins (P = 0.015). Our study provides novel clues for understanding the biological mechanism of vertigo, and highlights several possible risks and therapeutic chemicals for vertigo.

[Characteristics and regulatory mechanisms of lipid metabolism remodeling after malignant transformation of glioma-associated macrophages]

Objective: To observe the lipid metabolism characteristics of tumor-associated macrophages (TAM) after malignant transformation in the glioma micro-environment, and analyze the biological phenotype changes and regulatory mechanisms after inhibiting the lipid metabolism remodeling. Methods: Twelve male Balb/c mice of 6-8 weeks were used in the study. Macrophages (Mφ) were derived from mouse bone marrow, and malignantly transformed macrophages (tMφ₁ and tMφ₂) were cloned from the model of glioma stem cell (GSC) through interaction with Mφ in vivo and in vitro. Intracellular lipid droplet formation and cellular cholesterol content were measured respectively in Mφ, tMφ₁ and tMφ₂. qRT-PCR was performed to detect the genes expression level related with lipid metabolism, including sterol regulatory element binding protein (SREBP), fatty acid synthase (FASN), and 3-hydroxy-3-methylglutarate monoacyl coenzyme A reductase (HMG-CoA). Simvastatin (SIM) was used to analyze the proliferation, immigration and invasiveness ability in tMφ₁ and tMφ₂ after inhibition of the lipid metabolism. Differential expression profiles of miRNAs after SIM treatment were constructed in t-Mφ₁ and bio-informatics analysis was screened and verified for miR449a and its target gene sorting micro-tubule connectin 17 (SNX17) associated with lipid metabolism remodeling. The effect on SNX17 by up-regulated miR-449a were analyzed by qRT-PCR and Western blot, meanwhile, the biological phenotype and cholesterol content were observed after up-regulation of miR449a. Low-density lipoprotein receptor (LDLR) protein levels after SNX17 knockdown and intracellular cholesterol content after LDLR knockdown were detected respectively. Results: The numbers of intracellular lipid droplet formation in tMφ₁ and tMφ₂ were more than that in Mφ (P<0.001). Likewise, the relative contents of cholesterol (3.89±0.68 and 3.56±0.53), SREBP (4.78±0.60 and 2.84±0.41), FASN (4.65±0.70 and 3.01±0.45), and HMG-CoA (5.74±0.55 and 2.97±0.34) were significantly higher in tMφ₁ and tMφ₂ than those of Mφ (1.01 wel, 1.02 wel and 0.99 wel, respectively) (all P<0.001). The proliferation rates of tMφ₁ and tMφ₂ decreased from (47.06±5.88) % and (45.29±5.64)% to (23.53±4.70)% and (18.74±5.76)%, respectively after treatment with SIM (both P<0.05). The numbers of migrated cells decreased from 1 025±138 and 350±47 to 205±63 and 99±25, respectively (both P<0.001). And the numbers of invasiveness cells decreased from 919±45 and 527±34 to 220±23 and 114±21, respectively (both P<0.001). While the relative intracellular cholesterol content decreased to 0.52±0.08 and 0.58±0.07 (both P<0.05), respectively. MiR-449a was screened from tMφ₁ by SIM, and the target gene was analyzed and verified to be SNX17. SNX17 expression was down-regulated, and the proliferation rate, the number of migration and invasiveness was significantly decreased after miR-449a over-expression (all P<0.05). Low-density lipoprotein receptor (LDLR) expression was down-regulated after knock-down of SNX17, while the cholesterol content was decreased after knock-down of LDLR in tMφ₁ and tMφ₂ (all P<0.05). Conclusions: Malignantly transformed TAMs undergo lipid metabolism remodeling characterized with enhanced lipid metabolism. MiR-449a regulates the LDLR by targeting SNX17, thereby affecting the lipid metabolism of malignantly transformed macrophages, and subsequently inhibiting its proliferation, migration, and invasion ability. Precise intervention with miR-449a/SNX17/LDLR axis could provide an experimental basis for reversing its tumor-promoting micro-environment remodeled by GSC through metabolic intervention.

Evaluating the role of rare genetic variation in sleep duration

OBJECTIVES

To explore the roles of rare and high-impact variants in sleep duration.

DESIGN

Based on the recently released UK Biobank 200k exome dataset, an exome-wide association study was conducted to detect rare variants (minor allele frequency <0.01) contributing to sleep duration. Variant annotations were performed by the software tool ANNOVAR. Gene-based burden tests of sleep duration were conducted by the SKAT R-package. After quality control, 137,047 subjects were included in this study. CAUSALdb database was used to explore the related mental traits of identified genes.

RESULTS

We detected 730,572 variants with MAF < 1%, including 3873 frameshift variants, 3977 nonframeshift variants, 449,632 nonsynonymous variants, 1293 startloss variants, 10,254 stopgain variants, 413 stoploss variants, 261,130 synonymous variants, and 3102 variants are annotated as unknown. The burden test of exonic variants detected two exome-wide significant associations for sleep duration including TMIE at 3p21.31 (P_{Bonferroni adjusted} = 0.015) and ZIC2 at 13q32.3 (P_{Bonferroni adjusted} = 0.047). There are only nonsynonymous contained in TMIE; as for ZIC2, we detected 2 annotations of variants: nonsynonymous (P_{Bonferroni adjusted} =2.04 × 10^-4) and nonframeshift (P_{Bonferroni adjusted} =0.85). TMIE and ZIC2 were reported to be associated with several mental traits, such as chronotype, depression, and brain natriuretic peptide in published study.

CONCLUSION

This study reported 2 novel candidate genes for a sleep duration, supporting the roles of rare genetic variants in the regulation of sleep duration.

Advances in physiological mechanisms of selenium to improve heavy metal stress tolerance in plants

Selenium (Se) is a metalloid mineral nutrient for human and animal health. Plants are the main foodstuff source of the Se intake of humans. For plants, the addition of an appropriate amount of Se could promotes growth and development, and improves the tolerance to environmental stress, especially stress from some of heavy metals (HM) stress, such as cadmium (Cd) and mercury (Hg). This paper mainly reviews and summarizes the physiological mechanism of Se in enhancing HM stress tolerance in plants. The antagonistic effect of Se on HM is a comprehensive effect that includes many physiological mechanisms. Se can promote the removal of excessive reactive oxygen species and reduce the oxidative damage of plant cells under HM elements stress. Se participates in the regulation of the transportation and distribution of HM ions in plants, and alleviates the damage caused by of HM stress. Moreover, Se combine with HM elements to form Se-HM complexes and promote the production of phytochelatins (PCs), thereby reducing the accumulation of HM ions in plants. Overall, Se plays an important role in plant response to HM stress, but current studies mainly focus on physiological mechanism, and further in-depth study on the molecular mechanism is essential to confirm the participation of Se in plant response to environmental stress. This review helps to comprehensively understand the physiological mechanism of Se in plant tolerance against to HM stress of plants, and provides important theoretical support for the practical application of Se in environmental remediation and agricultural development.

New evaluation of the thermodynamics stability for bcc-Fe

The thermodynamic properties for bcc-Fe were predicted by combination of the first-principles calculations, the quasiharmonic approximation, the CALPHAD method and the Weiss molecular field theory. The hybrid method considers the effects of the lattice vibration, electron, intrinsic magnetism and external magnetic fields on the thermodynamic properties at finite temperature. Combined with experimental data, the calculated heat capacity without external magnetic fields was used to verify the validity of the hybrid method. Close to the Fermi level the high electronic density of states leads to a significant electronic contribution to free energy. Near the Curie temperature lattice vibrations dominant the Gibbs free energy. The order of the other three excitation contributions to Gibbs free energy from high to low is: intrinsic magnetism > electron > external magnetic fields. The investigation suggests that all the excitation contributions to Gibbs free energy are not negligible which provides a correct direction for tuning the thermodynamic properties for Fe-based alloy.

Index of multiple deprivation contributed to common psychiatric disorders: A systematic review and comprehensive analysis

BACKGROUND

Limited studies have been conducted to explore the interaction effects of social environmental and genetic factors on the risks of common psychiatric disorders.

METHODS

56,613-106,695 individuals were collected from the UK Biobank cohort. Logistic or liner regression models were first used to evaluate the associations of index of multiple deprivation (IMD) with bipolar disorder (BD), depression and anxiety in UK Biobank cohort. Then, for the significant IMD associated with BD, depression and anxiety, genome-wide gene-environment interaction study (GWEIS) was performed by PLINK 2.0.

RESULT

Totally, the higher levels of IMD were significantly associated with higher risks of BD, depression and anxiety. For BD, GWEIS identified multiple significant SNPs interacting with IMD, such as rs75182167 for income and rs111841503 for education. For depression and anxiety, GWEIS found significant SNPs interacting with income and education, such as rs147013419 for income and rs142366753 for education.

CONCLUSION

Social environmental deprivations contributed to the risks of psychiatric disorders. Besides, we reported multiple candidate genetic loci interacting with IMD, providing novel insights into the biological mechanism.