The harmonizome: a collection of processed datasets gathered to serve and mine knowledge about genes and proteins

Crosstalk between degradation and bioenergetics: how autophagy and endolysosomal processes regulate energy production

Cells undergo metabolic reprogramming to adapt to changes in nutrient availability, cellular activity, and transitions in cell states. The balance between glycolysis and mitochondrial respiration is crucial for energy production, and metabolic reprogramming stipulates a shift in such balance to optimize both bioenergetic efficiency and anabolic requirements. Failure in switching bioenergetic dependence can lead to maladaptation and pathogenesis. While cellular degradation is known to recycle precursor molecules for anabolism, its potential role in regulating energy production remains less explored. The bioenergetic switch between glycolysis and mitochondrial respiration involves transcription factors and organelle homeostasis, which are both regulated by the cellular degradation pathways. A growing body of studies has demonstrated that both stem cells and differentiated cells exhibit bioenergetic switch upon perturbations of autophagic activity or endolysosomal processes. Here, we highlighted the current understanding of the interplay between degradation processes, specifically autophagy and endolysosomes, transcription factors, endolysosomal signaling, and mitochondrial homeostasis in shaping cellular bioenergetics. This review aims to summarize the relationship between degradation processes and bioenergetics, providing a foundation for future research to unveil deeper mechanistic insights into bioenergetic regulation.

Dimerization of hub protein DYNLL1 and bZIP transcription factor CREB3L1 enhances transcriptional activation of CREB3L1 target genes like arginine vasopressin

bZIP transcription factors can function as homodimers or heterodimers through interactions with their disordered coiled-coil domain. Such dimer assemblies are known to influence DNA-binding specificity and/or the recruitment of binding partners, which can cause a functional switch of a transcription factor from being an activator to a repressor. We recently identified the genomic targets of a bZIP transcription factor called CREB3L1 in rat hypothalamic supraoptic nucleus by ChIP-seq. The objective of this study was to investigate the CREB3L1 protein-to-protein interactome of which little is known. For this approach, we created and screened a rat supraoptic nucleus yeast two-hybrid prey library with the bZIP region of rat CREB3L1 as the bait. Our yeast two-hybrid approach captured five putative CREB3L1 interacting prey proteins in the supraoptic nucleus. One interactor was selected by bioinformatic analyses for more detailed investigation by co-immunoprecipitation, immunofluorescent cellular localisation, and reporter assays in vitro. Here we identify dimerisation hub protein Dynein Light Chain LC8-Type 1 as a CREB3L1 interacting protein that in vitro enhances CREB3L1 activation of target genes.

A multiomic characterization of the leukemia cell line REH using short- and long-read sequencing

The B-cell acute lymphoblastic leukemia (ALL) cell line REH, with the t(12;21) ETV6::RUNX1 translocation, is known to have a complex karyotype defined by a series of large-scale chromosomal rearrangements. Taken from a 15-yr-old at relapse, the cell line offers a practical model for the study of pediatric B-ALL. In recent years, short- and long-read DNA and RNA sequencing have emerged as a complement to karyotyping techniques in the resolution of structural variants in an oncological context. Here, we explore the integration of long-read PacBio and Oxford Nanopore whole-genome sequencing, IsoSeq RNA sequencing, and short-read Illumina sequencing to create a detailed genomic and transcriptomic characterization of the REH cell line. Whole-genome sequencing clarified the molecular traits of disrupted ALL-associated genes including CDKN2A, PAX5, BTG1, VPREB1, and TBL1XR1, as well as the glucocorticoid receptor NR3C1 Meanwhile, transcriptome sequencing identified seven fusion genes within the genomic breakpoints. Together, our extensive whole-genome investigation makes high-quality open-source data available to the leukemia genomics community.

Transcription factor 7 like 2 promotes metastasis in hepatocellular carcinoma via NEDD9-mediated activation of AKT/mTOR signaling pathway

BACKGROUND

Hepatocellular carcinoma (HCC) is one of the most common malignant tumors of the digestive system, and the exact mechanism of HCC is still unclear. Transcription factor 7 like 2 (TCF7L2) plays a pivotal role in cell proliferation and stemness maintenance. However, the exact mechanism of TCF7L2 in HCC remains unclear.

METHODS

Clinical samples and public databases were used to analyze the expression and prognosis of TCF7L2 in HCC. The function of TCF7L2 in HCC was studied in vitro and in vivo. ChIP and luciferase assays were used to explore the molecular mechanism of TCF7L2. The relationship between TCF7L2 and NEDD9 was verified in HCC clinical samples by tissue microarrays.

RESULTS

The expression of TCF7L2 was upregulated in HCC, and high expression of TCF7L2 was associated with poor prognosis of HCC patients. Overexpression of TCF7L2 promoted the metastasis of HCC in vitro and in vivo, while Knockdown of TCF7L2 showed the opposite effect. Mechanically, TCF7L2 activated neural precursor cell expressed developmentally downregulated protein 9 (NEDD9) transcription by binding to the -1522/-1509 site of the NEDD9 promoter region, thereby increasing the phosphorylation levels of AKT and mTOR. The combination of TCF7L2 and NEDD9 could distinguish the survival of HCC patients.

CONCLUSIONS

This study demonstrated that TCF7L2 promotes HCC metastasis by activating AKT/mTOR pathway in a NEDD9-dependent manner, suggesting that potential of TCF7L2 and NEDD9 as prognostic markers and therapeutic targets for HCC.

The Xenorhabdus nematophila LrhA transcriptional regulator modulates production of γ-keto-N-acyl amides with inhibitory activity against mutualistic host nematode egg hatching

Xenorhabdus nematophila is a symbiotic Gammaproteobacterium that produces diverse natural products that facilitate mutualistic and pathogenic interactions in their nematode and insect hosts, respectively. The interplay between X. nematophila secondary metabolism and symbiosis stage is tuned by various global regulators. An example of such a regulator is the LysR-type protein transcription factor LrhA, which regulates amino acid metabolism and is necessary for virulence in insects and normal nematode progeny production. Here, we utilized comparative metabolomics and molecular networking to identify small molecule factors regulated by LrhA and characterized a rare γ-ketoacid (GKA) and two new N-acyl amides, GKA-Arg (1) and GKA-Pro (2) which harbor a γ-keto acyl appendage. A lrhA null mutant produced elevated levels of compound 1 and reduced levels of compound 2 relative to wild type. N-acyl amides 1 and 2 were shown to be selective agonists for the human G-protein-coupled receptors (GPCRs) C3AR1 and CHRM2, respectively. The CHRM2 agonist 2 deleteriously affected the hatch rate and length of Steinernema nematodes. This work further highlights the utility of exploiting regulators of host-bacteria interactions for the identification of the bioactive small molecule signals that they control.

IMPORTANCE

Xenorhabdus bacteria are of interest due to their symbiotic relationship with Steinernema nematodes and their ability to produce a variety of natural bioactive compounds. Despite their importance, the regulatory hierarchy connecting specific natural products and their regulators is poorly understood. In this study, comparative metabolomic profiling was utilized to identify the secondary metabolites modulated by the X. nematophila global regulator LrhA. This analysis led to the discovery of three metabolites, including an N-acyl amide that inhibited the egg hatching rate and length of Steinernema carpocapsae nematodes. These findings support the notion that X. nematophila LrhA influences the symbiosis between X. nematophila and S. carpocapsae through N-acyl amide signaling. A deeper understanding of the regulatory hierarchy of these natural products could contribute to a better comprehension of the symbiotic relationship between X. nematophila and S. carpocapsae.

Involvement of the choroid plexus in Alzheimer's disease pathophysiology: findings from mouse and human proteomic studies

BACKGROUND

Structural and functional changes of the choroid plexus (ChP) have been reported in Alzheimer's disease (AD). Nonetheless, the role of the ChP in the pathogenesis of AD remains largely unknown. We aim to unravel the relation between ChP functioning and core AD pathogenesis using a unique proteomic approach in mice and humans.

METHODS

We used an APP knock-in mouse model, APPNL-G-F, exhibiting amyloid pathology, to study the association between AD brain pathology and protein changes in mouse ChP tissue and CSF using liquid chromatography mass spectrometry. Mouse proteomes were investigated at the age of 7 weeks (n = 5) and 40 weeks (n = 5). Results were compared with previously published human AD CSF proteomic data (n = 496) to identify key proteins and pathways associated with ChP changes in AD.

RESULTS

ChP tissue proteome was dysregulated in APPNL-G-F mice relative to wild-type mice at both 7 and 40 weeks. At both ages, ChP tissue proteomic changes were associated with epithelial cells, mitochondria, protein modification, extracellular matrix and lipids. Nonetheless, some ChP tissue proteomic changes were different across the disease trajectory; pathways related to lysosomal function, endocytosis, protein formation, actin and complement were uniquely dysregulated at 7 weeks, while pathways associated with nervous system, immune system, protein degradation and vascular system were uniquely dysregulated at 40 weeks. CSF proteomics in both mice and humans showed similar ChP-related dysregulated pathways.

CONCLUSIONS

Together, our findings support the hypothesis of ChP dysfunction in AD. These ChP changes were related to amyloid pathology. Therefore, the ChP could become a novel promising therapeutic target for AD.

Analysis of the Actions of RARγ Agonists on Growing Osteochondromas in a Mouse Model

The actions of the retinoic acid nuclear receptor gamma (RARγ) agonist, palovarotene, on pre-existing osteochondromas were investigated using a mouse multiple osteochondroma model. This approach was based on the knowledge that patients often present to the clinic after realizing the existence of osteochondroma masses, and the findings from preclinical investigations are the effects of drugs on the initial formation of osteochondromas. Systemic administration of palovarotene, with increased doses (from 1.76 to 4.0 mg/kg) over time, fully inhibited tumor growth, keeping the tumor size (0.31 ± 0.049 mm3) similar to the initial size (0.27 ± 0.031 mm3, p = 0.66) while the control group tumor grew (1.03 ± 0.23 mm3, p = 0.023 to the drug-treated group). Nanoparticle (NP)-based local delivery of the RARγ agonist also inhibited the growth of osteochondromas at an early stage (Control: 0.52 ± 0.11 mm3; NP: 0.26 ± 0.10, p = 0.008). Transcriptome analysis revealed that the osteoarthritis pathway was activated in cultured chondrocytes treated with palovarotene (Z-score = 2.29), with the upregulation of matrix catabolic genes and the downregulation of matrix anabolic genes, consistent with the histology of palovarotene-treated osteochondromas. A reporter assay performed in cultured chondrocytes demonstrated that the Stat3 pathway, but not the Stat1/2 pathway, was stimulated by RARγ agonists. The activation of Stat3 by palovarotene was confirmed using immunoblotting and immunohistochemistry. These findings suggest that palovarotene treatment is effective against pre-existing osteochondromas and that the Stat3 pathway is involved in the antitumor actions of palovarotene.

CSF proteomic profiles of neurodegeneration biomarkers in Alzheimer's disease

INTRODUCTION

We aimed to unravel the underlying pathophysiology of the neurodegeneration (N) markers neurogranin (Ng), neurofilament light (NfL), and hippocampal volume (HCV), in Alzheimer's disease (AD) using cerebrospinal fluid (CSF) proteomics.

METHODS

Individuals without dementia were classified as A+ (CSF amyloid beta [Aβ]42), T+ (CSF phosphorylated tau181), and N+ or N- based on Ng, NfL, or HCV separately. CSF proteomics were generated and compared between groups using analysis of covariance.

RESULTS

Only a few individuals were A+T+Ng-. A+T+Ng+ and A+T+NfL+ showed different proteomic profiles compared to A+T+Ng- and A+T+NfL-, respectively. Both Ng+ and NfL+ were associated with neuroplasticity, though in opposite directions. Compared to A+T+HCV-, A+T+HCV+ showed few proteomic changes, associated with oxidative stress.

DISCUSSION

Different N markers are associated with distinct neurodegenerative processes and should not be equated. N markers may differentially complement disease staging beyond amyloid and tau. Our findings suggest that Ng may not be an optimal N marker, given its low incongruency with tau pathophysiology.

HIGHLIGHTS

In Alzheimer's disease, neurogranin (Ng)+, neurofilament light (NfL)+, and hippocampal volume (HCV)+ showed differential protein expression in cerebrospinal fluid. Ng+ and NfL+ were associated with neuroplasticity, although in opposite directions. HCV+ showed few proteomic changes, related to oxidative stress. Neurodegeneration (N) markers may differentially refine disease staging beyond amyloid and tau. Ng might not be an optimal N marker, as it relates more closely to tau.

Elevating PLK1 overcomes BETi resistance in prostate cancer via triggering BRD4 phosphorylation-dependent degradation in mitosis

Bromodomain-containing protein 4 (BRD4) has emerged as a promising therapeutic target in prostate cancer (PCa). Understanding the mechanisms of BRD4 stability could enhance the clinical response to BRD4-targeted therapy. In this study, we report that BRD4 protein levels are significantly decreased during mitosis in a PLK1-dependent manner. Mechanistically, we show that BRD4 is primarily phosphorylated at T1186 by the CDK1/cyclin B complex, recruiting PLK1 to phosphorylate BRD4 at S24/S1100, which are recognized by the APC/CCdh1 complex for proteasome pathway degradation. We find that PLK1 overexpression lowers SPOP mutation-stabilized BRD4, consequently rendering PCa cells re-sensitized to BRD4 inhibitors. Intriguingly, we report that sequential treatment of docetaxel and JQ1 resulted in significant inhibition of PCa. Collectively, the results support that PLK1-phosphorylated BRD4 triggers its degradation at M phase. Sequential treatment of docetaxel and JQ1 overcomes BRD4 accumulation-associated bromodomain and extra-terminal inhibitor (BETi) resistance, which may shed light on the development of strategies to treat PCa.

A liver immune rheostat regulates CD8 T cell immunity in chronic HBV infection

Chronic hepatitis B virus (HBV) infection affects 300 million patients worldwide1,2, in whom virus-specific CD8 T cells by still ill-defined mechanisms lose their function and cannot eliminate HBV-infected hepatocytes3-7. Here we demonstrate that a liver immune rheostat renders virus-specific CD8 T cells refractory to activation and leads to their loss of effector functions. In preclinical models of persistent infection with hepatotropic viruses such as HBV, dysfunctional virus-specific CXCR6+ CD8 T cells accumulated in the liver and, as a characteristic hallmark, showed enhanced transcriptional activity of cAMP-responsive element modulator (CREM) distinct from T cell exhaustion. In patients with chronic hepatitis B, circulating and intrahepatic HBV-specific CXCR6+ CD8 T cells with enhanced CREM expression and transcriptional activity were detected at a frequency of 12-22% of HBV-specific CD8 T cells. Knocking out the inhibitory CREM/ICER isoform in T cells, however, failed to rescue T cell immunity. This indicates that CREM activity was a consequence, rather than the cause, of loss in T cell function, further supported by the observation of enhanced phosphorylation of protein kinase A (PKA) which is upstream of CREM. Indeed, we found that enhanced cAMP-PKA-signalling from increased T cell adenylyl cyclase activity augmented CREM activity and curbed T cell activation and effector function in persistent hepatic infection. Mechanistically, CD8 T cells recognizing their antigen on hepatocytes established close and extensive contact with liver sinusoidal endothelial cells, thereby enhancing adenylyl cyclase-cAMP-PKA signalling in T cells. In these hepatic CD8 T cells, which recognize their antigen on hepatocytes, phosphorylation of key signalling kinases of the T cell receptor signalling pathway was impaired, which rendered them refractory to activation. Thus, close contact with liver sinusoidal endothelial cells curbs the activation and effector function of HBV-specific CD8 T cells that target hepatocytes expressing viral antigens by means of the adenylyl cyclase-cAMP-PKA axis in an immune rheostat-like fashion.

SFINN: inferring gene regulatory network from single-cell and spatial transcriptomic data with shared factor neighborhood and integrated neural network

MOTIVATION

The rise of single-cell RNA sequencing (scRNA-seq) technology presents new opportunities for constructing detailed cell type-specific gene regulatory networks (GRNs) to study cell heterogeneity. However, challenges caused by noises, technical errors, and dropout phenomena in scRNA-seq data pose significant obstacles to GRN inference, making the design of accurate GRN inference algorithms still essential. The recent growth of both single-cell and spatial transcriptomic sequencing data enables the development of supervised deep learning methods to infer GRNs on these diverse single-cell datasets.

RESULTS

In this study, we introduce a novel deep learning framework based on shared factor neighborhood and integrated neural network (SFINN) for inferring potential interactions and causalities between transcription factors and target genes from single-cell and spatial transcriptomic data. SFINN utilizes shared factor neighborhood to construct cellular neighborhood network based on gene expression data and additionally integrates cellular network generated from spatial location information. Subsequently, the cell adjacency matrix and gene pair expression are fed into an integrated neural network framework consisting of a graph convolutional neural network and a fully-connected neural network to determine whether the genes interact. Performance evaluation in the tasks of gene interaction and causality prediction against the existing GRN reconstruction algorithms demonstrates the usability and competitiveness of SFINN across different kinds of data. SFINN can be applied to infer GRNs from conventional single-cell sequencing data and spatial transcriptomic data.

AVAILABILITY AND IMPLEMENTATION

SFINN can be accessed at GitHub: https://github.com/JGuan-lab/SFINN.

Identifying new cancer genes based on the integration of annotated gene sets via hypergraph neural networks

MOTIVATION

Identifying cancer genes remains a significant challenge in cancer genomics research. Annotated gene sets encode functional associations among multiple genes, and cancer genes have been shown to cluster in hallmark signaling pathways and biological processes. The knowledge of annotated gene sets is critical for discovering cancer genes but remains to be fully exploited.

RESULTS

Here, we present the DIsease-Specific Hypergraph neural network (DISHyper), a hypergraph-based computational method that integrates the knowledge from multiple types of annotated gene sets to predict cancer genes. First, our benchmark results demonstrate that DISHyper outperforms the existing state-of-the-art methods and highlight the advantages of employing hypergraphs for representing annotated gene sets. Second, we validate the accuracy of DISHyper-predicted cancer genes using functional validation results and multiple independent functional genomics data. Third, our model predicts 44 novel cancer genes, and subsequent analysis shows their significant associations with multiple types of cancers. Overall, our study provides a new perspective for discovering cancer genes and reveals previously undiscovered cancer genes.

AVAILABILITY AND IMPLEMENTATION

DISHyper is freely available for download at https://github.com/genemine/DISHyper.

Elucidating prognosis in cervical squamous cell carcinoma and endocervical adenocarcinoma: a novel anoikis-related gene signature model

Background

Cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC) are among the most prevalent gynecologic malignancies globally. The prognosis is abysmal once cervical cancer progresses to lymphatic metastasis. Anoikis, a specialized form of apoptosis induced by loss of cell adhesion to the extracellular matrix, plays a critical role. The prediction model based on anoikis-related genes (ARGs) expression and clinical data could greatly aid clinical decision-making. However, the relationship between ARGs and CESC remains unclear.

Methods

ARGs curated from the GeneCards and Harmonizome portals were instrumental in delineating CESC subtypes and in developing a prognostic framework for patients afflicted with this condition. We further delved into the intricacies of the immune microenvironment and pathway enrichment across the identified subtypes. Finally, our efforts culminated in the creation of an innovative nomogram that integrates ARGs. The utility of this prognostic tool was underscored by Decision Curve Analysis (DCA), which illuminate its prospective benefits in guiding clinical interventions.

Results

In our study, We discerned a set of 17 survival-pertinent, anoikis-related differentially expressed genes (DEGs) in CESC, from which nine were meticulously selected for the construction of prognostic models. The derived prognostic risk score was subsequently validated as an autonomous prognostic determinant. Through comprehensive functional analyses, we observed distinct immune profiles and drug response patterns among divergent prognostic stratifications. Further, we integrated the risk scores with the clinicopathological characteristics of CESC to develop a robust nomogram. DCA corroborated the utility of our model, demonstrating its potential to enhance patient outcomes through tailored clinical treatment strategies.

Conclusion

The predictive signature, encompassing nine pivotal genes, alongside the meticulously constructed nomogram developed in this research, furnishes clinicians with a sophisticated tool for tailoring treatment strategies to individual patients diagnosed with CESC.

Endonucleolytic RNA cleavage drives changes in gene expression during the innate immune response

Viral infection triggers several double-stranded RNA (dsRNA) sensors that lead to changes in gene expression in the cell. One of these sensors activates an endonuclease, ribonuclease L (RNase L), that cleaves single-stranded RNA. However, how the resultant widespread RNA fragmentation affects gene expression is not fully understood. Here, we show that this fragmentation induces the ribotoxic stress response via ZAKα, potentially through stalled ribosomes and/or ribosome collisions. The p38 and JNK pathways that are activated as part of this response promote outcomes that inhibit the virus, such as programmed cell death. We also show that RNase L limits the translation of stress-responsive genes. Intriguingly, we found that the activity of the generic endonuclease, RNase A, recapitulates many of the same molecular phenotypes as activated RNase L, demonstrating how widespread RNA cleavage can evoke an antiviral program.

The non-mitotic role of HMMR in regulating the localization of TPX2 and the dynamics of microtubules in neurons

A functional nervous system is built upon the proper morphogenesis of neurons to establish the intricate connection between them. The microtubule cytoskeleton is known to play various essential roles in this morphogenetic process. While many microtubule-associated proteins (MAPs) have been demonstrated to participate in neuronal morphogenesis, the function of many more remains to be determined. This study focuses on a MAP called HMMR in mice, which was originally identified as a hyaluronan binding protein and later found to possess microtubule and centrosome binding capacity. HMMR exhibits high abundance on neuronal microtubules and altering the level of HMMR significantly affects the morphology of neurons. Instead of confining to the centrosome(s) like cells in mitosis, HMMR localizes to microtubules along axons and dendrites. Furthermore, transiently expressing HMMR enhances the stability of neuronal microtubules and increases the formation frequency of growing microtubules along the neurites. HMMR regulates the microtubule localization of a non-centrosomal microtubule nucleator TPX2 along the neurite, offering an explanation for how HMMR contributes to the promotion of growing microtubules. This study sheds light on how cells utilize proteins involved in mitosis for non-mitotic functions.

Systematic Analysis of Biological Processes Reveals Gene Co-expression Modules Driving Pathway Dysregulation in Alzheimer's Disease

Alzheimer's disease (AD) manifests as a complex systems pathology with intricate interplay among various genes and biological processes. Traditional differential gene expression (DEG) analysis, while commonly employed to characterize AD-driven perturbations, does not sufficiently capture the full spectrum of underlying biological processes. Utilizing single-nucleus RNA-sequencing data from postmortem brain samples across key regions-middle temporal gyrus, superior frontal gyrus, and entorhinal cortex-we provide a comprehensive systematic analysis of disrupted processes in AD. We go beyond the DEG-centric analysis by integrating pathway activity analysis with weighted gene co-expression patterns to comprehensively map gene interconnectivity, identifying region- and cell-type-specific drivers of biological processes associated with AD. Our analysis reveals profound modular heterogeneity in neurons and glia as well as extensive AD-related functional disruptions. Co-expression networks highlighted the extended involvement of astrocytes and microglia in biological processes beyond neuroinflammation, such as calcium homeostasis, glutamate regulation, lipid metabolism, vesicle-mediated transport, and TOR signaling. We find limited representation of DEGs within dysregulated pathways across neurons and glial cells, suggesting that differential gene expression alone may not adequately represent the disease complexity. Further dissection of inferred gene modules revealed distinct dynamics of hub DEGs in neurons versus glia, suggesting that DEGs exert more impact on neurons compared to glial cells in driving modular dysregulations underlying perturbed biological processes. Interestingly, we observe an overall downregulation of astrocyte and microglia modules across all brain regions in AD, indicating a prevailing trend of functional repression in glial cells across these regions. Notable genes from the CALM and HSP90 families emerged as hub genes across neuronal modules in all brain regions, suggesting conserved roles as drivers of synaptic dysfunction in AD. Our findings demonstrate the importance of an integrated, systems-oriented approach combining pathway and network analysis to comprehensively understand the cell-type-specific roles of genes in AD-related biological processes.

Inhibiting AGTR1 reduces AML burden and protects the heart from cardiotoxicity in mouse models

Clinical treatment of acute myeloid leukemia (AML) largely relies on intensive chemotherapy. However, the application of chemotherapy is often hindered by cardiotoxicity. Patient sequence data revealed that angiotensin II receptor type 1 (AGTR1) is a shared target between AML and cardiovascular disease (CVD). We found that inhibiting AGTR1 sensitized AML to chemotherapy and protected the heart against chemotherapy-induced cardiotoxicity in a human AML cell-transplanted mouse model. These effects were regulated by the AGTR1-Notch1 axis in AML cells and cardiomyocytes from mice. In mouse cardiomyocytes, AGTR1 was hyperactivated by AML and chemotherapy. AML leukemogenesis increased the expression of the angiotensin-converting enzyme and led to increased production of angiotensin II, the ligand of AGTR1, in an MLL-AF9-driven AML mouse model. In this model, the AGTR1-Notch1 axis regulated a variety of genes involved with cell stemness and chemotherapy resistance. AML cell stemness was reduced after Agtr1a deletion in the mouse AML cell transplant model. Mechanistically, Agtr1a deletion decreased γ-secretase formation, which is required for transmembrane Notch1 cleavage and release of the Notch1 intracellular domain into the nucleus. Using multiomics, we identified AGTR1-Notch1 signaling downstream genes and found decreased binding between these gene sequences with Notch1 and chromatin enhancers, as well as increased binding with silencers. These findings describe an AML/CVD association that may be used to improve AML treatment.

Discovery of genomic loci for liver health and steatosis reveals overlap with glutathione redox genetics

Non-alcoholic fatty liver disease (NAFLD) is the most common chronic liver condition in the United States, encompassing a wide spectrum of liver pathologies including steatosis, steatohepatitis, fibrosis, and cirrhosis. Despite its high prevalence, there are no medications currently approved by the Food and Drug Administration for the treatment of NAFLD. Recent work has suggested that NAFLD has a strong genetic component and identifying causative genes will improve our understanding of the molecular mechanisms contributing to NAFLD and yield targets for future therapeutic investigations. Oxidative stress is known to play an important role in NAFLD pathogenesis, yet the underlying mechanisms accounting for disturbances in redox status are not entirely understood. To better understand the relationship between the glutathione redox system and signs of NAFLD in a genetically-diverse population, we measured liver weight, serum biomarkers aspartate aminotransferase (AST) and alanine aminotransferase (ALT), and graded liver pathology in a large cohort of Diversity Outbred mice. We compared hepatic endpoints to those of the glutathione redox system previously measured in the livers and kidneys of the same mice, and we screened for statistical and genetic associations using the R/qtl2 software. We discovered several novel genetic loci associated with markers of liver health, including loci that were associated with both liver steatosis and glutathione redox status. Candidate genes within each locus point to possible new mechanisms underlying the complex relationship between NAFLD and the glutathione redox system, which could have translational implications for future studies targeting NAFLD pathology.

Plasma Extracellular MicroRNAs Associated With Cardiovascular Disease Risk Factors in Middle-Aged and Older Adults

BACKGROUND

Extracellular microRNAs (miRNAs) are a class of noncoding RNAs that remain stable in the extracellular milieu, where they contribute to various physiological and pathological processes by facilitating intercellular signaling. Previous studies have reported associations between miRNAs and cardiovascular diseases (CVDs); however, the plasma miRNA signatures of CVD and its risk factors have not been fully elucidated at the population level.

METHODS AND RESULTS

Plasma miRNA levels were measured in 4440 FHS (Framingham Heart Study) participants. Linear regression analyses were conducted to test the cross-sectional associations of each miRNA with 8 CVD risk factors. Prospective analyses of the associations of miRNAs with new-onset obesity, hypertension, type 2 diabetes, CVD, and all-cause mortality were conducted using proportional hazards regression. Replication was carried out in 1999 RS (Rotterdam Study) participants. Pathway enrichment analyses were conducted and target genes were predicted for miRNAs associated with ≥5 risk factors in the FHS. In the FHS, 6 miRNAs (miR-193b-3p, miR-122-5p, miR-365a-3p, miR-194-5p, miR-192-5p, and miR-193a-5p) were associated with ≥5 risk factors. This miRNA signature was enriched for pathways associated with CVD and several genes annotated to these pathways were predicted targets of the identified miRNAs. Furthermore, miR-193b-3p, miR-194-5p, and miR-193a-5p were each associated with ≥2 risk factors in the RS. Prospective analysis revealed 8 miRNAs associated with all-cause mortality in the FHS.

CONCLUSIONS

These findings highlight associations between miRNAs and CVD risk factors that may provide valuable insights into the underlying pathogenesis of CVD.

Single-cell signatures identify microenvironment factors in tumors associated with patient outcomes

The cellular components of tumors and their microenvironment play pivotal roles in tumor progression, patient survival, and the response to cancer treatments. Unveiling a comprehensive cellular profile within bulk tumors via single-cell RNA sequencing (scRNA-seq) data is crucial, as it unveils intrinsic tumor cellular traits that elude identification through conventional cancer subtyping methods. Our contribution, scBeacon, is a tool that derives cell-type signatures by integrating and clustering multiple scRNA-seq datasets to extract signatures for deconvolving unrelated tumor datasets on bulk samples. Through the employment of scBeacon on the The Cancer Genome Atlas (TCGA) cohort, we find cellular and molecular attributes within specific tumor categories, many with patient outcome relevance. We developed a tumor cell-type map to visually depict the relationships among TCGA samples based on the cell-type inferences.

Progressive senescence programs induce intrinsic vulnerability to aging-related female breast cancer

Cancer incidence escalates exponentially with advancing age; however, the underlying mechanism remains unclear. In this study, we build a chronological molecular clock at single-cell transcription level with a mammary stem cell-enriched population to depict physiological aging dynamics in female mice. We find that the mammary aging process is asynchronous and progressive, initiated by an early senescence program, succeeded by an entropic late senescence program with elevated cancer associated pathways, vulnerable to cancer predisposition. The transition towards senescence program is governed by a stem cell factor Bcl11b, loss of which accelerates mammary ageing with enhanced DMBA-induced tumor formation. We have identified a drug TPCA-1 that can rejuvenate mammary cells and significantly reduce aging-related cancer incidence. Our findings establish a molecular portrait of progressive mammary cell aging and elucidate the transcriptional regulatory network bridging mammary aging and cancer predisposition, which has potential implications for the management of cancer prevalence in the aged.

Repurposing Castanea sativa Spiny Burr By-Products Extract as a Potentially Effective Anti-Inflammatory Agent for Novel Future Biotechnological Applications

The concept of a "circular bioeconomy" holds great promise for the health, cosmetic, and nutrition sectors by re-using Castanea sativa (Mill.) by-products. This sustainable resource is rich in bioactive secondary metabolites with antioxidant and anti-inflammatory properties. By transforming these by-products into high-value products for human health, we can promote sustainable economic growth and reduce the environmental impact of traditional waste disposal, adding value to previously underutilized resources. In the present study, we investigated the antioxidant capacity, phytochemical composition, and in vitro antioxidant and anti-inflammatory activity of C. sativa burr (CSB) aqueous extract. The spectrophotometric study revealed high total phenolic content (TPC) values with significant antioxidant and anti-radical properties. Using UPLC-MS/MS techniques, the phytochemical investigation identified 56 metabolites, confirming the presence of phenolic compounds in CSBs. In addition, CSBs significantly downregulated pro-inflammatory mediators in LPS-stimulated RAW 264.7 macrophage cells without significant cell toxicity. Lastly, in silico studies pinpointed three kinases from RAW 264.7 cells as binding partners with ellagic acid, the predominant compound found in our extract. These findings strongly advocate for the recycling and valorization of C. sativa by-products, challenging their conventional classification as mere "waste".

Context-Dependent Regulation of Peripheral Nerve Abundance by the PI3K Pathway in the Tumor Microenvironment of Head and Neck Squamous Cell Carcinoma

Recent studies have highlighted neurons and their associated Schwann cells (SCs) as key regulators of cancer development. However, the mode of their interaction with tumor cells or other components of the tumor microenvironment (TME) remains elusive. We established an SC-related 43-gene set as a surrogate for peripheral nerves in the TME. Head and neck squamous cell carcinoma (HNSCC) from The Cancer Genome Atlas (TCGA) were classified into low, intermediate and high SC score groups based on the expression of this gene set. Perineural invasion (PNI) and TGF-β signaling were hallmarks of SChigh tumors, whereas SClow tumors were enriched for HPV16-positive OPSCC and higher PI3K-MTOR activity. The latter activity was partially explained by a higher frequency of PTEN mutation and PIK3CA copy number gain. The inverse association between PI3K-MTOR activity and peripheral nerve abundance was context-dependent and influenced by the TP53 mutation status. An in silico drug screening approach highlighted the potential vulnerabilities of HNSCC with variable SC scores and predicted a higher sensitivity of SClow tumors to DNA topoisomerase inhibitors. In conclusion, we have established a tool for assessing peripheral nerve abundance in the TME and provided new clinical and biological insights into their regulation. This knowledge may pave the way for new therapeutic strategies and impart proof of concept in appropriate preclinical models.

System analysis based on Anoikis-related genes identifies MAPK1 as a novel therapy target for osteosarcoma with neoadjuvant chemotherapy

BACKGROUND

Osteosarcoma (OS) is the most common bone malignant tumor in children, and its prognosis is often poor. Anoikis is a unique mode of cell death.However, the effects of Anoikis in OS remain unexplored.

METHOD

Differential analysis of Anoikis-related genes was performed based on the metastatic and non-metastatic groups. Then LASSO logistic regression and SVM-RFE algorithms were applied to screen out the characteristic genes. Later, Univariate and multivariate Cox regression was conducted to identify prognostic genes and further develop the Anoikis-based risk score. In addition, correlation analysis was performed to analyze the relationship between tumor microenvironment, drug sensitivity, and prognostic models.

RESULTS

We established novel Anoikis-related subgroups and developed a prognostic model based on three Anoikis-related genes (MAPK1, MYC, and EDIL3). The survival and ROC analysis results showed that the prognostic model was reliable. Besides, the results of single-cell sequencing analysis suggested that the three prognostic genes were closely related to immune cell infiltration. Subsequently, aberrant expression of two prognostic genes was identified in osteosarcoma cells. Nilotinib can promote the apoptosis of osteosarcoma cells and down-regulate the expression of MAPK1.

CONCLUSIONS

We developed a novel Anoikis-related risk score model, which can assist clinicians in evaluating the prognosis of osteosarcoma patients in clinical practice. Analysis of the tumor immune microenvironment and chemotherapeutic drug sensitivity can provide necessary insights into subsequent mechanisms. MAPK1 may be a valuable therapeutic target for neoadjuvant chemotherapy in osteosarcoma.

Chronic TNF exposure induces glucocorticoid-like immunosuppression in the alveolar macrophages of aged mice that enhances their susceptibility to pneumonia

Chronic low-grade inflammation, particularly elevated tumor necrosis factor (TNF) levels, occurs due to advanced age and is associated with greater susceptibility to infection. One reason for this is age-dependent macrophage dysfunction (ADMD). Herein, we use the adoptive transfer of alveolar macrophages (AM) from aged mice into the airway of young mice to show that inherent age-related defects in AM were sufficient to increase the susceptibility to Streptococcus pneumoniae, a Gram-positive bacterium and the leading cause of community-acquired pneumonia. MAPK phosphorylation arrays using AM lysates from young and aged wild-type (WT) and TNF knockout (KO) mice revealed multilevel TNF-mediated suppression of kinase activity in aged mice. RNAseq analyses of AM validated the suppression of MAPK signaling as a consequence of TNF during aging. Two regulatory phosphatases that suppress MAPK signaling, Dusp1 and Ptprs, were confirmed to be upregulated with age and as a result of TNF exposure both ex vivo and in vitro. Dusp1 is known to be responsible for glucocorticoid-mediated immune suppression, and dexamethasone treatment increased Dusp1 and Ptprs expression in cells and recapitulated the ADMD phenotype. In young mice, treatment with dexamethasone increased the levels of Dusp1 and Ptprs and their susceptibility to infection. TNF-neutralizing antibody reduced Dusp1 and Ptprs levels in AM from aged mice and reduced pneumonia severity following bacterial challenge. We conclude that chronic exposure to TNF increases the expression of the glucocorticoid-associated MAPK signaling suppressors, Dusp1 and Ptprs, which inhibits AM activation and increases susceptibility to bacterial pneumonia in older adults.

Chemo-small extracellular vesicles released in cisplatin-resistance ovarian cancer cells are regulated by the lysosomal function

Chemoresistance is a common problem in ovarian cancer (OvCa) treatment, where resistant cells, in response to chemotherapy, secrete small extracellular vesicles (sEVs), known as chemo-sEVs, that transfer resistance to recipient cells. sEVs are formed as intraluminal vesicles (ILVs) within multivesicular endosomes (MVEs), whose trafficking is regulated by Ras-associated binding (RAB) GTPases that mediate sEVs secretion or lysosomal degradation. A decrease in lysosomal function can promote sEVs secretion, but the relationship between MVEs trafficking pathways and sEVs secretion in OvCa chemoresistance is unclear. Here, we show that A2780cis cisplatin (CCDP) resistant OvCa cells had an increased number of MVEs and ILVs structures, higher levels of Endosomal Sorting Complex Required for Transport (ESCRTs) machinery components, and RAB27A compared to A2780 CDDP-sensitive OvCa cells. CDDP promoted the secretion of chemo-sEVs in A2780cis cells, enriched in DNA damage response proteins. A2780cis cells exhibited poor lysosomal function with reduced levels of RAB7, essential in MVEs-Lysosomal trafficking. The silencing of RAB27A in A2780cis cells prevents the Chemo-EVs secretion, reduces its chemoresistance and restores lysosomal function and levels of RAB7, switching them into an A2780-like cellular phenotype. Enhancing lysosomal function with rapamycin reduced chemo-sEVs secretion. Our results suggest that adjusting the balance between secretory MVEs and lysosomal MVEs trafficking could be a promising strategy for overcoming CDDP chemoresistance in OvCa.

Novel Approach for Pancreas Transcriptomics Reveals the Cellular Landscape in Homeostasis and Acute Pancreatitis

BACKGROUND & AIMS

Acinar cells produce digestive enzymes that impede transcriptomic characterization of the exocrine pancreas. Thus, single-cell RNA-sequencing studies of the pancreas underrepresent acinar cells relative to histological expectations, and a robust approach to capture pancreatic cell responses in disease states is needed. We sought to innovate a method that overcomes these challenges to accelerate study of the pancreas in health and disease.

METHODS

We leverage FixNCut, a single-cell RNA-sequencing approach in which tissue is reversibly fixed with dithiobis(succinimidyl propionate) before dissociation and single-cell preparation. We apply FixNCut to an established mouse model of acute pancreatitis, validate findings using GeoMx whole transcriptome atlas profiling, and integrate our data with prior studies to compare our method in both mouse and human pancreas datasets.

RESULTS

FixNCut achieves unprecedented definition of challenging pancreatic cells, including acinar and immune populations in homeostasis and acute pancreatitis, and identifies changes in all major cell types during injury and recovery. We define the acinar transcriptome during homeostasis and acinar-to-ductal metaplasia and establish a unique gene set to measure deviation from normal acinar identity. We characterize pancreatic immune cells, and analysis of T-cell subsets reveals a polarization of the homeostatic pancreas toward type-2 immunity. We report immune responses during acute pancreatitis and recovery, including early neutrophil infiltration, expansion of dendritic cell subsets, and a substantial shift in the transcriptome of macrophages due to both resident macrophage activation and monocyte infiltration.

CONCLUSIONS

FixNCut preserves pancreatic transcriptomes to uncover novel cell states during homeostasis and following pancreatitis, establishing a broadly applicable approach and reference atlas for study of pancreas biology and disease.

A Diagnostic Model for Parkinson's Disease Based on Anoikis-Related Genes

Parkinson's disease (PD) is the second most prevalent neurodegenerative disease, and its pathological mechanisms are thought to be closely linked to apoptosis. Anoikis, a specific type of apoptosis, has recently been suggested to play a role in the progression of Parkinson's disease; however, the underlying mechanisms are not well understood. To explore the potential mechanisms involved in PD, we selected genes from the GSE28894 dataset and compared their expression in PD patients and healthy controls to identify differentially expressed genes (DEGs), and selected anoikis-related genes (ANRGs) from the DEGs. Furthermore, the least absolute shrinkage and selection operator (LASSO) regression approach and multivariate logistic regression highlighted five key genes-GSK3B, PCNA, CDC42, DAPK2, and SRC-as biomarker candidates. Subsequently, we developed a nomogram model incorporating these 5 genes along with age and sex to predict and diagnose PD. To evaluate the model's coherence, clinical applicability, and distinguishability, we utilized receiver operating characteristic (ROC) curves, the C-index, and calibration curves and validated it in both the GSE20295 dataset and our center's external clinical data. In addition, we confirmed the differential expression of the 5 model genes in human blood samples through qRT-PCR and Western blotting. Our constructed anoikis-related PD diagnostic model exhibits satisfactory predictive accuracy and offers novel insights into both diagnosis and treatment strategies for Parkinson's disease while facilitating its implementation in clinical practice.

Anoikis and SPP1 in idiopathic pulmonary fibrosis: integrating bioinformatics, cell, and animal studies to explore prognostic biomarkers and PI3K/AKT signaling regulation

OBJECTIVE

This study aims to explore the relevance of anoikis in idiopathic pulmonary fibrosis (IPF) and identify associated biomarkers and signaling pathways.

METHOD

Unsupervised consensus cluster analysis was employed to categorize IPF patients into subtypes. We utilized Weighted Gene Co-Expression Network Analysis (WGCNA) and Protein-Protein Interaction network construction to identify anoikis-related modules and key genes. A prognostic signature was developed using Lasso and multivariate Cox regression analysis. Single-cell sequencing assessed hub gene expression in various cell types, and both cell and animal experiments confirmed IPF-related pathways.

RESULTS

We identified two distinct anoikis-associated subtypes with differing prognoses. WGCNA revealed essential hub genes, with SPP1 being prominent in the anoikis-related signature. The anoikis-related signature is effective in determining the prognosis of patients with IPF. Single-cell sequencing highlighted significant differences in SPP1 expression, notably elevated in fibroblasts derived from IPF patients. In vivo and in vitro experiments demonstrated that SPP1 enhances fibrosis in mouse lung fibroblasts by regulating p27 through the PI3K/Akt pathway.

CONCLUSION

Our research demonstrates a robust prognostic signature associated with anoikis and highlights SPP1 as a pivotal regulator of the PI3K/AKT signaling pathway in pulmonary fibrosis.

Aggressive high-grade NF2 mutant meningiomas downregulate oncogenic YAP signaling via the upregulation of VGLL4 and FAT3/4

Meningiomas are the most common primary brain tumors in adults. Although generally benign, a subset of meningiomas is of higher grade, shows aggressive growth behavior and recurs even after multiple surgeries. Around half of all meningiomas harbor inactivating mutations in NF2. While benign low-grade NF2 mutant meningiomas exhibit few genetic events in addition to NF2 inactivation, aggressive high-grade NF2 mutant meningiomas frequently harbor a highly aberrant genome. We and others have previously shown that NF2 inactivation leads to YAP1 activation and that YAP1 acts as the pivotal oncogenic driver in benign NF2 mutant meningiomas. Using bulk and single-cell RNA-Seq data from a large cohort of human meningiomas, we show that aggressive NF2 mutant meningiomas harbor decreased levels YAP1 activity compared to their benign counterparts. Decreased expression levels of YAP target genes are significantly associated with an increased risk of recurrence. We then identify the increased expression of the YAP1 competitor VGLL4 as well as the YAP1 upstream regulators FAT3/4 as a potential mechanism for the downregulation of YAP activity in aggressive NF2 mutant meningiomas. High expression of these genes is significantly associated with an increased risk of recurrence. In vitro, overexpression of VGLL4 resulted in the downregulation of YAP activity in benign NF2 mutant meningioma cells, confirming the direct link between VGLL4 expression and decreased levels of YAP activity observed in aggressive NF2 mutant meningiomas. Our results shed new insight on the biology of benign and aggressive NF2 mutant meningiomas and may have important implications for the efficacy of therapies targeting oncogenic YAP1 activity in NF2 mutant meningiomas.

Inferring upstream regulatory genes of FOXP3 in human regulatory T cells from time-series transcriptomic data

The discovery of upstream regulatory genes of a gene of interest still remains challenging. Here we applied a scalable computational method to unbiasedly predict candidate regulatory genes of critical transcription factors by searching the whole genome. We illustrated our approach with a case study on the master regulator FOXP3 of human primary regulatory T cells (Tregs). While target genes of FOXP3 have been identified, its upstream regulatory machinery still remains elusive. Our methodology selected five top-ranked candidates that were tested via proof-of-concept experiments. Following knockdown, three out of five candidates showed significant effects on the mRNA expression of FOXP3 across multiple donors. This provides insights into the regulatory mechanisms modulating FOXP3 transcriptional expression in Tregs. Overall, at the genome level this represents a high level of accuracy in predicting upstream regulatory genes of key genes of interest.

Exploring cutting-edge strategies for drug repurposing in female cancers - An insight into the tools of the trade

Female cancers, which include breast and gynaecological cancers, represent a significant global health burden for women. Despite advancements in research pertinent to unearthing crucial pathological characteristics of these cancers, challenges persist in discovering potential therapeutic strategies. This is further exacerbated by economic burdens associated with de novo drug discovery and clinical intricacies such as development of drug resistance and metastasis. Drug repurposing, an innovative approach leveraging existing FDA-approved drugs for new indications, presents a promising avenue to expedite therapeutic development. Computational techniques, including virtual screening and analysis of drug-target-disease relationships, enable the identification of potential candidate drugs. Integration of diverse data types, such as omics and clinical information, enhances the precision and efficacy of drug repurposing strategies. Experimental approaches, including high-throughput screening assays, in vitro, and in vivo models, complement computational methods, facilitating the validation of repurposed drugs. This review highlights various target mining strategies based on analysis of differential gene expression, weighted gene co-expression, protein-protein interaction network, and host-pathogen interaction, among others. To unearth drug candidates, the technicalities of leveraging information from databases such as DrugBank, STITCH, LINCS, and ChEMBL, among others are discussed. Further in silico validation techniques encompassing molecular docking, pharmacophore modelling, molecular dynamic simulations, and ADMET analysis are elaborated. Overall, this review delves into the exploration of individual case studies to offer a wide perspective of the ever-evolving field of drug repurposing, emphasizing the multifaceted approaches and methodologies employed for the same to confront female cancers.

The impact of chronic pain on brain gene expression

Background

Chronic pain affects one fifth of American adults, contributing significant public health burden. Chronic pain mechanisms can be further understood through investigating brain gene expression.

Methods

We tested differentially expressed genes (DEGs) in chronic pain, migraine, lifetime fentanyl and oxymorphone use, and with chronic pain genetic risk in four brain regions (dACC, DLPFC, MeA, BLA) and imputed cell type expression data from 304 postmortem donors. We compared findings across traits and with independent transcriptomics resources, and performed gene-set enrichment.

Results

We identified two chronic pain DEGs: B4GALT and VEGFB in bulk dACC. We found over 2000 (primarily BLA microglia) chronic pain cell type DEGs. Findings were enriched for mouse microglia pain genes, and for hypoxia and immune response. Cross-trait DEG overlap was minimal.

Conclusions

Chronic pain-associated gene expression is heterogeneous across cell type, largely distinct from that in pain-related traits, and shows BLA microglia are a key cell type.

A gray box framework that optimizes a white box logical model using a black box optimizer for simulating cellular responses to perturbations

Predicting cellular responses to perturbations requires interpretable insights into molecular regulatory dynamics to perform reliable cell fate control, despite the confounding non-linearity of the underlying interactions. There is a growing interest in developing machine learning-based perturbation response prediction models to handle the non-linearity of perturbation data, but their interpretation in terms of molecular regulatory dynamics remains a challenge. Alternatively, for meaningful biological interpretation, logical network models such as Boolean networks are widely used in systems biology to represent intracellular molecular regulation. However, determining the appropriate regulatory logic of large-scale networks remains an obstacle due to the high-dimensional and discontinuous search space. To tackle these challenges, we present a scalable derivative-free optimizer trained by meta-reinforcement learning for Boolean network models. The logical network model optimized by the trained optimizer successfully predicts anti-cancer drug responses of cancer cell lines, while simultaneously providing insight into their underlying molecular regulatory mechanisms.

Bidirectional epigenetic editing reveals hierarchies in gene regulation

CRISPR perturbation methods are limited in their ability to study non-coding elements and genetic interactions. In this study, we developed a system for bidirectional epigenetic editing, called CRISPRai, in which we apply activating (CRISPRa) and repressive (CRISPRi) perturbations to two loci simultaneously in the same cell. We developed CRISPRai Perturb-seq by coupling dual perturbation gRNA detection with single-cell RNA sequencing, enabling study of pooled perturbations in a mixed single-cell population. We applied this platform to study the genetic interaction between two hematopoietic lineage transcription factors, SPI1 and GATA1, and discovered novel characteristics of their co-regulation on downstream target genes, including differences in SPI1 and GATA1 occupancy at genes that are regulated through different modes. We also studied the regulatory landscape of IL2 (interleukin-2) in Jurkat T cells, primary T cells and chimeric antigen receptor (CAR) T cells and elucidated mechanisms of enhancer-mediated IL2 gene regulation. CRISPRai facilitates investigation of context-specific genetic interactions, provides new insights into gene regulation and will enable exploration of non-coding disease-associated variants.

Phytochemical Composition, Anti-Inflammatory Property, and Anti-Atopic Effect of Chaetomorpha linum Extract

Utilizing plant-based resources, particularly their by-products, aligns with sustainability principles and circular bioeconomy, contributing to environmental preservation. The therapeutic potential of plant extracts is garnering increasing interest, and this study aimed to demonstrate promising outcomes from an extract obtained from an underutilized plant waste. Chaetomorpha linum, an invasive macroalga found in the Orbetello Lagoon, thrives in eutrophic conditions, forming persistent mats covering approximately 400 hectares since 2005. The biomass of C. linum undergoes mechanical harvesting and is treated as waste, requiring significant human efforts and economic resources-A critical concern for municipalities. Despite posing challenges to local ecosystems, the study identified C. linum as a natural source of bioactive metabolites. Phytochemical characterization revealed lipids, amino acids, and other compounds with potential anti-inflammatory activity in C. linum extract. In vitro assays with LPS-stimulated RAW 264.7 and TNF-α/IFN-γ-stimulated HaCaT cells showed the extract inhibited reactive oxygen species (ROS), nitric oxide (NO), and prostaglandin E2 (PGE2) productions, and reduced inducible nitric oxide synthase (iNOS) and cyclooxygenase-2 (COX-2) expressions via NF-κB nuclear translocation, in RAW 264.7 cells. It also reduced chemokines (TARC/CCL17, RANTES/CCL5, MCP-1/CCL2, and IL-8) and the cytokine IL-1β production in HaCaT cells, suggesting potential as a therapeutic candidate for chronic diseases like atopic dermatitis. Finally, in silico studies indicated palmitic acid as a significant contributor to the observed effect. This research not only uncovered the untapped potential of C. linum but also laid the foundation for its integration into the circular bioeconomy, promoting sustainable practices, and innovative applications across various industries.

Disentangling cell-intrinsic and extrinsic factors underlying gene expression evolution

Chimeras have played a foundational role in biology, for example by enabling the classification of developmental processes into those driven intrinsically by individual cells versus those driven extrinsically by their extracellular environment. Here, we extend this framework to decompose evolutionary divergence in gene expression and other quantitative traits into cell-intrinsic, extrinsic, and intrinsic-extrinsic interaction components. Applying this framework to reciprocal rat-mouse chimeras, we found that the majority of gene expression divergence is attributable to cell-intrinsic factors, though extrinsic factors also play an integral role. For example, a rat-like extracellular environment extrinsically up-regulates the expression of a key transcriptional regulator of the endoplasmic reticulum (ER) stress response in some but not all cell types, which in turn strongly predicts extrinsic up-regulation of its target genes and of the ER stress response pathway as a whole. This effect is also seen at the protein level, suggesting propagation through multiple regulatory levels. We also demonstrate that our framework is applicable to a cellular trait, neuronal differentiation, and estimated the intrinsic and extrinsic contributions to its divergence. Finally, we show that imprinted genes are dramatically mis-expressed in species-mismatched environments, suggesting that mismatch between rapidly evolving intrinsic and extrinsic mechanisms controlling gene imprinting may contribute to barriers to interspecies chimerism. Overall, our conceptual framework opens new avenues to investigate the mechanistic basis of evolutionary divergence in gene expression and other quantitative traits in any multicellular organism.

S100A9 interacts with a dynamic region on CD14 to activate Toll-like receptor 4

S100A9 is a Damage Associated Molecular Pattern (DAMP) that activates inflammatory pathways via Toll-like receptor 4 (TLR4). This activity plays important homeostatic roles in tissue repair, but can also contribute to inflammatory diseases. The mechanism of activation is unknown. Here, we follow up on a previous observation that the protein CD14 is an important co-receptor that enables S100A9 to activate TLR4. Using cell-based functional assays and a combination of mutations and pharmocological perturbations, we found that CD14 must be membrane bound to potentiate TLR4 activation by S100A9. Additionally, S100A9 is sensitive to inhibitors of pathways downstream of TLR4 internalization. Together, this suggests that S100A9 induces activity via CD14-dependent internalization of TLR4. We then used mutagenesis, structural modeling, and in vitro binding experiments to establish that S100A9 binds to CD14's N-terminus in a region that overlaps with, but is not identical to, the region where CD14 binds its canonical ligand, lipopolysaccharide (LPS). In molecular dynamics simulations, this region of the protein is dynamic, allowing it to reorganize to recognize both S100A9 (a soluble protein) and LPS (a small hydrophobic molecule). Our work is the first attempt at a molecular characterization of the S100A9/CD14 interaction, bringing us one step closer to unraveling the full mechanism by which S100A9 activates TLR4/MD-2.

Early life exposure to vitamin D deficiency impairs molecular mechanisms that regulate liver cholesterol biosynthesis, energy metabolism, inflammation, and detoxification

Introduction

Emerging data suggests liver disease may be initiated during development when there is high genome plasticity and the molecular pathways supporting liver function are being developed.

Methods

Here, we leveraged our Collaborative Cross mouse model of developmental vitamin D deficiency (DVD) to investigate the role of DVD in dysregulating the molecular mechanisms underlying liver disease. We defined the effects on the adult liver transcriptome and metabolome and examined the role of epigenetic dysregulation. Given that the parental origin of the genome (POG) influences response to DVD, we used our established POG model [POG1-(CC011xCC001)F1 and POG2-(CC001xCC011)F1] to identify interindividual differences.

Results

We found that DVD altered the adult liver transcriptome, primarily downregulating genes controlling liver development, response to injury/infection (detoxification & inflammation), cholesterol biosynthesis, and energy production. In concordance with these transcriptional changes, we found that DVD decreased liver cell membrane-associated lipids (including cholesterol) and pentose phosphate pathway metabolites. Each POG also exhibited distinct responses. POG1 exhibited almost 2X more differentially expressed genes (DEGs) with effects indicative of increased energy utilization. This included upregulation of lipid and amino acid metabolism genes and increased intermediate lipid and amino acid metabolites, increased energy cofactors, and decreased energy substrates. POG2 exhibited broader downregulation of cholesterol biosynthesis genes with a metabolomics profile indicative of decreased energy utilization. Although DVD primarily caused loss of liver DNA methylation for both POGs, only one epimutation was shared, and POG2 had 6.5X more differentially methylated genes. Differential methylation was detected at DEGs regulating developmental processes such as amino acid transport (POG1) and cell growth & differentiation (e.g., Wnt & cadherin signaling, POG2).

Conclusions

These findings implicate a novel role for maternal vitamin D in programming essential offspring liver functions that are dysregulated in liver disease. Importantly, impairment of these processes was not rescued by vitamin D treatment at weaning, suggesting these effects require preventative measures. Substantial differences in POG response to DVD demonstrate that the parental genomic context of exposure determines offspring susceptibility.

SAFB regulates hippocampal stem cell fate by targeting Drosha to destabilize Nfib mRNA

Neural stem cells (NSCs) are multipotent and correct fate determination is crucial to guarantee brain formation and homeostasis. How NSCs are instructed to generate neuronal or glial progeny is not well understood. Here, we addressed how murine adult hippocampal NSC fate is regulated and described how scaffold attachment factor B (SAFB) blocks oligodendrocyte production to enable neuron generation. We found that SAFB prevents NSC expression of the transcription factor nuclear factor I/B (NFIB) by binding to sequences in the Nfib mRNA and enhancing Drosha-dependent cleavage of the transcripts. We show that increasing SAFB expression prevents oligodendrocyte production by multipotent adult NSCs, and conditional deletion of Safb increases NFIB expression and oligodendrocyte formation in the adult hippocampus. Our results provide novel insights into a mechanism that controls Drosha functions for selective regulation of NSC fate by modulating the post-transcriptional destabilization of Nfib mRNA in a lineage-specific manner.

Variant-to-function analysis of the childhood obesity chr12q13 locus implicates rs7132908 as a causal variant within the 3' UTR of FAIM2

The ch12q13 locus is among the most significant childhood obesity loci identified in genome-wide association studies. This locus resides in a non-coding region within FAIM2; thus, the underlying causal variant(s) presumably influence disease susceptibility via cis-regulation. We implicated rs7132908 as a putative causal variant by leveraging our in-house 3D genomic data and public domain datasets. Using a luciferase reporter assay, we observed allele-specific cis-regulatory activity of the immediate region harboring rs7132908. We generated isogenic human embryonic stem cell lines homozygous for either rs7132908 allele to assess changes in gene expression and chromatin accessibility throughout a differentiation to hypothalamic neurons, a key cell type known to regulate feeding behavior. The rs7132908 obesity risk allele influenced expression of FAIM2 and other genes and decreased the proportion of neurons produced by differentiation. We have functionally validated rs7132908 as a causal obesity variant that temporally regulates nearby effector genes and influences neurodevelopment and survival.

Aging hallmarks of the primate ovary revealed by spatiotemporal transcriptomics

The ovary is indispensable for female reproduction, and its age-dependent functional decline is the primary cause of infertility. However, the molecular basis of ovarian aging in higher vertebrates remains poorly understood. Herein, we apply spatiotemporal transcriptomics to benchmark architecture organization as well as cellular and molecular determinants in young primate ovaries and compare these to aged primate ovaries. From a global view, somatic cells within the non-follicle region undergo more pronounced transcriptional fluctuation relative to those in the follicle region, likely constituting a hostile microenvironment that facilitates ovarian aging. Further, we uncovered that inflammation, the senescent-associated secretory phenotype, senescence, and fibrosis are the likely primary contributors to ovarian aging (PCOA). Of note, we identified spatial co-localization between a PCOA-featured spot and an unappreciated MT2 (Metallothionein 2) highly expressing spot (MT2high) characterized by high levels of inflammation, potentially serving as an aging hotspot in the primate ovary. Moreover, with advanced age, a subpopulation of MT2high accumulates, likely disseminating and amplifying the senescent signal outward. Our study establishes the first primate spatiotemporal transcriptomic atlas, advancing our understanding of mechanistic determinants underpinning primate ovarian aging and unraveling potential biomarkers and therapeutic targets for aging and age-associated human ovarian disorders.

Transcriptomic, epigenomic, and spatial metabolomic cell profiling redefines regional human kidney anatomy

A large-scale multimodal atlas that includes major kidney regions is lacking. Here, we employed simultaneous high-throughput single-cell ATAC/RNA sequencing (SHARE-seq) and spatially resolved metabolomics to profile 54 human samples from distinct kidney anatomical regions. We generated transcriptomes of 446,267 cells and chromatin accessibility profiles of 401,875 cells and developed a package to analyze 408,218 spatially resolved metabolomes. We find that the same cell type, including thin limb, thick ascending limb loop of Henle and principal cells, display distinct transcriptomic, chromatin accessibility, and metabolomic signatures, depending on anatomic location. Surveying metabolism-associated gene profiles revealed non-overlapping metabolic signatures between nephron segments and dysregulated lipid metabolism in diseased proximal tubule (PT) cells. Integrating multimodal omics with clinical data identified PLEKHA1 as a disease marker, and its in vitro knockdown increased gene expression in PT differentiation, suggesting possible pathogenic roles. This study highlights previously underrepresented cellular heterogeneity underlying the human kidney anatomy.

CHD7 regulates craniofacial cartilage development via controlling HTR2B expression

Mutations in the Chromodomain helicase DNA-binding protein 7 - coding gene (CHD7) cause CHARGE syndrome (CS). Although craniofacial and skeletal abnormalities are major features of CS patients, the role of CHD7 in bone and cartilage development remain largely unexplored. Here, using a zebrafish (Danio rerio) CS model, we show that chd7-/- larvae display abnormal craniofacial cartilage development and spinal deformities. The craniofacial and spine defects are accompanied by a marked reduction of bone mineralization. At the molecular level, we show that these phenotypes are associated with significant reduction in the expression levels of osteoblast differentiation markers. Additionally, we detected a marked depletion of collagen 2α1 in the cartilage of craniofacial regions and vertebrae, along with significantly reduced number of chondrocytes. Chondrogenesis defects are at least in part due to downregulation of htr2b, which we found to be also dysregulated in human cells derived from an individual with CHD7 mutation-positive CS. Overall, this study thus unveils an essential role for CHD7 in cartilage and bone development, with potential clinical relevance for the craniofacial defects associated with CS.

Attenuated cell cycle and DNA damage response transcriptome signatures and overrepresented cell adhesion processes imply accelerated progression in patients with lower-risk myelodysplastic neoplasms

Patients with myelodysplastic neoplasms (MDS) are classified according to the risk of acute myeloid leukemia transformation. Some lower-risk MDS patients (LR-MDS) progress rapidly despite expected good prognosis. Using diagnostic samples, we aimed to uncover the mechanisms of this accelerated progression at the transcriptome level. RNAseq was performed on CD34+ ribodepleted RNA samples from 53 LR-MDS patients without accelerated progression (stMDS) and 8 who progressed within 20 months (prMDS); 845 genes were differentially expressed (ІlogFCІ > 1, FDR < 0.01) between these groups. stMDS CD34+ cells exhibited transcriptional signatures of actively cycling, megakaryocyte/erythrocyte lineage-primed progenitors, with upregulation of cell cycle checkpoints and stress pathways, which presumably form a tumor-suppressing barrier. Conversely, cell cycle, DNA damage response (DDR) and energy metabolism-related pathways were downregulated in prMDS samples, whereas cell adhesion processes were upregulated. Also, prMDS samples showed high levels of aberrant splicing and global lncRNA expression that may contribute to the attenuation of DDR pathways. We observed overexpression of multiple oncogenes and diminished differentiation in prMDS; the expression of ZEB1 and NEK3, genes not previously associated with MDS prognosis, might serve as potential biomarkers for LR-MDS progression. Our 19-gene DDR signature showed a significant predictive power for LR-MDS progression. In validation samples (stMDS = 3, prMDS = 4), the key markers and signatures retained their significance. Collectively, accelerated progression of LR-MDS appears to be associated with transcriptome patterns of a quiescent-like cell state, reduced lineage differentiation and suppressed DDR, inherent to CD34+ cells. The attenuation of DDR-related gene-expression signature may refine risk assessment in LR-MDS patients.

The immune cell transcriptome is modulated by vitamin D3 supplementation in people with a first demyelinating event participating in a randomized placebo-controlled trial

Vitamin D deficiency is a risk factor for developing multiple sclerosis. The PrevANZ trial was conducted to determine if vitamin D3 supplementation can prevent recurrent disease activity in people with a first demyelinating event. As a sub-study of this trial, we investigated the effect of supplementation on peripheral immune cell gene expression. Participants were randomized to 1000, 5000 or 10,000 international units daily of vitamin D3 or placebo. Peripheral blood was collected at baseline and 12 weeks and sent for ribonucleic acid sequencing. Datasets from 55 participants were included. Gene expression was modulated by high dose supplementation. Antigen presentation and viral response pathways were upregulated. Oxidative phosphorylation and immune signaling pathways, including tumor necrosis factor-alpha and interleukin-17 signaling, were downregulated. Overall, vitamin D3 supplementation for 12 weeks modulated the peripheral immune cell transcriptome with induction of anti-inflammatory gene expression profiles. Our results support a dose-dependent effect of vitamin D3 supplementation on immune gene expression.

Butyrate increases methylglyoxal production through regulation of the JAK2/Stat3/Nrf2/Glo1 pathway in castration‑resistant prostate cancer cells

Cancer cells are characterized by increased glycolysis, known as the Warburg effect, which leads to increased production of cytotoxic methylglyoxal (MGO) and apoptotic cell death. Cancer cells often activate the protective nuclear factor erythroid 2‑related factor2 (Nrf2)/glyoxalase1 (Glo1) system to detoxify MGO. The effects of sodium butyrate (NaB), a product of gut microbiota, on Nrf2/Glos/MGO pathway and the underlying mechanisms in prostate cancer (PCa) cells were investigated in the present study. Treatment with NaB induced the cell death and reduced the proliferation of PCa cells (DU145 and LNCap). Moreover, the protein kinase RNA-like endoplasmic reticulum kinase/Nrf2/Glo1 pathway was greatly inhibited by NaB, thereby accumulating MGO-derived adduct hydroimidazolone (MG-H1). In response to a high amount of MGO, the expression of Nrf2 and Glo1 was attenuated, coinciding with an increased cellular death. NaB also markedly inhibited the Janus kinase 2 (JAK2)/Signal transducer and activator of transcription 3 (Stat3) pathway. Conversely, co‑treatment with Colivelin, a Stat3 activator, significantly reversed the effects of NaB on Glo1 expression, MG-H1 production, and the cell migration and viability. As expected, overexpression of Stat3 or Glo1 reduced NaB‑induced cell death. The activation of calcium/calmodulin dependent protein kinase II gamma and reactive oxygen species production also contributed to the anticancer effect of NaB. The present study, for the first time, demonstrated that NaB greatly increases MGO production through suppression of the JAK2/Stat3/Nrf2/Glo1 pathway in DU145 cells, a cell line mimicking castration‑resistant PCa (CRPC), suggesting that NaB may be a potential agent for PCa therapy.

A lipidome landscape of aging in mice

Understanding the molecular mechanisms of aging is crucial for enhancing healthy longevity. We conducted untargeted lipidomics across 13 biological samples from mice at various life stages (2, 12, 19 and 24 months) to explore the potential link between aging and lipid metabolism, considering sex (male or female) and microbiome (specific pathogen-free or germ-free) dependencies. By analyzing 2,704 molecules from 109 lipid subclasses, we characterized common and tissue-specific lipidome alterations associated with aging. For example, the levels of bis(monoacylglycero)phosphate containing polyunsaturated fatty acids increased in various organs during aging, whereas the levels of other phospholipids containing saturated and monounsaturated fatty acids decreased. In addition, we discovered age-dependent sulfonolipid accumulation, absent in germ-free mice, correlating with Alistipes abundance determined by 16S ribosomal RNA gene amplicon sequencing. In the male kidney, glycolipids such as galactosylceramides, galabiosylceramides (Gal2Cer), trihexosylceramides (Hex3Cer), and mono- and digalactosyldiacylglycerols were detected, with two lipid classes-Gal2Cer and Hex3Cer-being significantly enriched in aged mice. Integrated analysis of the kidney transcriptome revealed uridine diphosphate galactosyltransferase 8A (UGT8a), alkylglycerone phosphate synthase and fatty acyl-coenzyme A reductase 1 as potential enzymes responsible for the male-specific glycolipid biosynthesis in vivo, which would be relevant to sex dependency in kidney diseases. Inhibiting UGT8 reduced the levels of these glycolipids and the expression of inflammatory cytokines in the kidney. Our study provides a valuable resource for clarifying potential links between lipid metabolism and aging.

Distinct roles of Bdnf I and Bdnf IV transcript variant expression in hippocampal neurons

Brain-derived neurotrophic factor (Bdnf) plays a critical role in brain development, dendritic growth, synaptic plasticity, as well as learning and memory. The rodent Bdnf gene contains nine 5' non-coding exons (I-IXa), which are spliced to a common 3' coding exon (IX). Transcription of individual Bdnf variants, which all encode the same BDNF protein, is initiated at unique promoters upstream of each non-coding exon, enabling precise spatiotemporal and activity-dependent regulation of Bdnf expression. Although prior evidence suggests that Bdnf transcripts containing exon I (Bdnf I) or exon IV (Bdnf IV) are uniquely regulated by neuronal activity, the functional significance of different Bdnf transcript variants remains unclear. To investigate functional roles of activity-dependent Bdnf I and IV transcripts, we used a CRISPR activation system in which catalytically dead Cas9 fused to a transcriptional activator (VPR) is targeted to individual Bdnf promoters with single guide RNAs, resulting in transcript-specific Bdnf upregulation. Bdnf I upregulation is associated with gene expression changes linked to dendritic growth, while Bdnf IV upregulation is associated with genes that regulate protein catabolism. Upregulation of Bdnf I, but not Bdnf IV, increased mushroom spine density, volume, length, and head diameter, and also produced more complex dendritic arbors in cultured rat hippocampal neurons. In contrast, upregulation of Bdnf IV, but not Bdnf I, in the rat hippocampus attenuated contextual fear expression. Our data suggest that while Bdnf I and IV are both activity-dependent, BDNF produced from these promoters may serve unique cellular, synaptic, and behavioral functions.

MEIS2 suppresses breast cancer development by downregulating IL10

BACKGROUND

Breast cancer (BC) is the most commonly diagnosed female cancer. Homeobox protein MEIS2, a key transcription factor, is involved in the regulation of many developmental and cellular processes. However, the role of MEIS2 in the development of breast cancer is still unclear.

AIMS

We aimed to examine the role of myeloid ecotropic insertion site (MEIS2) in breast cancer and the association of MEIS2 with breast cancer clinical stages and pathological grades. We revealed the underlying mechanism by which MEIS2 affected breast cancer cell growth and tumor development.

METHODS AND RESULTS

Using human BC cell lines, clinical samples and animal xenograft model, we reveal that MEIS2 functions as a tumor suppressor in breast cancer. The expression of MEIS2 is inversely correlated with BC clinical stages and pathological grades. MEIS2 knockdown (MEIS2-KD) promotes while MEIS2 overexpression suppresses breast cancer cell proliferation and tumor development in vitro and in animal xenograft models, respectively. To determine the biological function of MEIS2, we screen the expression of a group of MEIS2 potential targeting genes in stable-established cell lines. Results show that the knockdown of MEIS2 in breast cancer cells up-regulates the IL10 expression, but MEIS2 overexpression opposed the effect on IL10 expression. Furthermore, the suppressive role of MEIS2 in breast cancer cell proliferation is associated with the IL10 expression and myeloid cells infiltration.

CONCLUSION

Our study demonstrates that the tumor suppressor of MEIS2 in breast cancer progression is partially via down regulating the expression of IL10 and promoting myeloid cells infiltration. Targeting MEIS2 would be a potentially therapeutic avenue for BC.