Gene Set Knowledge Discovery with Enrichr

Inhibition of K63 ubiquitination by G-Protein pathway suppressor 2 (GPS2) regulates mitochondria-associated translation

G-Protein Pathway Suppressor 2 (GPS2) is an inhibitor of non-proteolytic K63 ubiquitination mediated by the E2 ubiquitin-conjugating enzyme Ubc13. Previous studies have associated GPS2-mediated restriction of ubiquitination with the regulation of insulin signaling, inflammatory responses and mitochondria-nuclear communication across different tissues and cell types. However, a detailed understanding of the targets of GPS2/Ubc13 activity is lacking. Here, we have dissected the GPS2-regulated K63 ubiquitome in mouse embryonic fibroblasts and human breast cancer cells, unexpectedly finding an enrichment for proteins involved in RNA binding and translation on the outer mitochondrial membrane. Validation of selected targets of GPS2-mediated regulation, including the RNA-binding protein PABPC1 and translation factors RPS1, RACK1 and eIF3M, revealed a mitochondrial-specific strategy for regulating the translation of nuclear-encoded mitochondrial proteins via non-proteolytic ubiquitination. Removal of GPS2-mediated inhibition, either via genetic deletion or stress-induced nuclear translocation, promotes the import-coupled translation of selected mRNAs leading to the increased expression of an adaptive antioxidant program. In light of GPS2 role in nuclear-mitochondria communication, these findings reveal an exquisite regulatory network for modulating mitochondrial gene expression through spatially coordinated transcription and translation.

Hepatitis C virus NS5A and core protein induce fibrosis-related genes regulation on Huh7 cells through activation of LX2 cells

INTRODUCTION AND OBJECTIVES

Liver fibrosis remains a complication derived from a chronic Hepatitis C Virus (HCV) infection even when it is resolved, and no liver antifibrotic drug has been approved. Molecular mechanisms on hepatocytes and activation of hepatic stellate cells (HSCs) play a central role in liver fibrogenesis. To elucidate molecular mechanisms, it is important to analyze pathway regulation during HSC activation and HCV infection.

MATERIALS AND METHODS

We evaluate the fibrosis-associated molecular mechanisms during a co-culture of human HSCs (LX2), with human hepatocytes (Huh7) that express HCV NS5A or Core protein. We evaluated LX2 activation induced by HCV NS5A or Core expression in Huh7 cells during co-culture. We determined a fibrosis-associated gene expression profile in Huh7 that expresses NS5A or Core proteins during the co-culture with LX2.

RESULTS

We observed that NS5A induced 8.3-, 6.7- and 4-fold changes and that Core induced 6.5-, 1.8-, and 6.2-fold changes in the collagen1, TGFβ1, and timp1 gene expression, respectively, in LX2 co-cultured with transfected Huh7. In addition, NS5A induced the expression of 30 genes while Core induced 41 genes and reduced the expression of 30 genes related to fibrosis in Huh7 cells during the co-culture with LX2, compared to control. The molecular pathways enriched from the gene expression profile were involved in TGFB signaling and the organization of extracellular matrix.

CONCLUSIONS

We demonstrated that HCV NS5A and Core protein expression regulate LX2 activation. NS5A and Core-induced LX2 activation, in turn, regulates diverse fibrosis-related gene expression at different levels in Huh7, which can be further analyzed as potential antifibrotic targets during HCV infection.

A genetic and proteomic comparison of key AD biomarkers across tissues

INTRODUCTION

Plasma has been proposed as an alternative to cerebrospinal fluid (CSF) for measuring Alzheimer's disease (AD) biomarkers, but no studies have analyzed in detail which biofluid is more informative for genetics studies of AD.

METHOD

Eleven proteins associated with AD (α-synuclein, apolipoprotein E [apoE], CLU, GFAP, GRN, NfL, NRGN, SNAP-25, TREM2, VILIP-1, YKL-40) were assessed in plasma (n = 2317) and CSF (n = 3107). Both plasma and CSF genome-wide association study (GWAS) analyses were performed for each protein, followed by functional annotation. Additional characterization for each biomarker included calculation of correlations and predictive power.

RESULTS

Eighteen plasma protein quantitative train loci (pQTLs) associated with 10 proteins and 16 CSF pQTLs associated with 9 proteins were identified. Plasma and CSF shared some genetic loci, but protein levels between tissues correlated weakly. CSF protein levels better associated with AD compared to plasma.

DISCUSSION

The present results indicate that CSF is more informative than plasma for genetic studies in AD.

HIGHLIGHTS

The identification of novel protein quantitative trait loci (pQTLs) in both plasma and cerebrospinal fluid (CSF). Plasma and CSF levels of neurodegeneration-related proteins correlated weakly. CSF is more informative than plasma for genetic studies of Alzheimer's disease (AD). Neurofilament light (NfL), triggering receptor expressed on myeloid cells 2 (TREM2), and chitinase-3-like protein 1 (YKL-40) tend to show relatively strong inter-tissue associations. A novel signal in the apolipoprotein E (APOE) region was identified, which is an eQTL for APOC1.

Complement protein expression changes in various conditions of breast cancer: in-silico analyses-experimental research

Introduction

Breast cancer is the most prevalent cancer diagnosed in females worldwide. The known biomarkers are insufficient to understand the actual prognosis of breast cancer, and identifying new biomarkers is desirable and valuable data to improve the patient's survival. Many inflammatory biomarkers, such as the complement system, can be regarded as prognostic values and as potent inflammatory mediators; complement proteins have a critical role in tumorigenesis. In the current study, the authors aim to investigate complement protein expression changes, particularly complement 3 (C3), complement 7 (C7), complement factor B (CFB), and complement factor D (CFD), in various conditions of breast cancer using in-silico tools.

Methods

The intent data were extracted using webtools, including; Kaplan-Meier plotter, BcGenExMiner, UALCAN, cbioportal, GeneMania, and Enrichr. To select valid data, a P greater than 0.05 was considered.

Results

The current study clarified that 21 complement genes correlated to survival conditions. Also, down or upregulation of extracted genes and breast cancer statuses were identified. Additionally, expression level difference of complement genes in various breast cancer four stages was detected. Ultimately, co-expression genes with complement genes were extracted and networked.

Conclusion

Changes in the expression of complement proteins can strongly correlate to breast cancer's prognosis, status, and survival. Furthermore, considering the vital role of CFD and CFB complement proteins in the alternative pathway in different stages of breast cancer, CFD and CFB can be regarded as reliable prognostic values for diagnosis.

Stable inhibition of choroidal neovascularization by adeno-associated virus 2/8-vectored bispecific molecules

Neovascular age-related macular degeneration (nAMD) causes severe visual impairment. Pigment epithelium-derived factor (PEDF), soluble CD59 (sCD59), and soluble fms-like tyrosine kinase-1 (sFLT-1) are potential therapeutic agents for nAMD, which target angiogenesis and the complement system. Using the AAV2/8 vector, two bi-target gene therapy agents, AAV2/8-PEDF-P2A-sCD59 and AAV2/8-sFLT-1-P2A-sCD59, were generated, and their therapeutic efficacy was investigated in laser-induced choroidal neovascularization (CNV) and Vldlr-/- mouse models. After a single injection, AAV2/8-mediated gene expression was maintained at high levels in the retina for two months. Both AAV2/8-PEDF-P2A-sCD59 and AAV2/8-sFLT-1-P2A-sCD59 significantly reduced CNV development for an extended period without side effects and provided efficacy similar to two injections of current anti-vascular endothelial growth factor monotherapy. Mechanistically, these agents suppressed the extracellular signal-regulated kinase and nuclear factor-κB pathways, resulting in anti-angiogenic activity. This study demonstrated the safety and long-lasting effects of AAV2/8-PEDF-P2A-sCD59 and AAV2/8-sFLT-1-P2A-sCD59 in CNV treatment, providing a promising therapeutic strategy for nAMD.

The Effect of Vitamin D Deficiency on Immune-Related Hub Genes: A Network Analysis Associated With Type 1 Diabetes

Background Type 1 diabetes (T1D) is an autoimmune disorder that results in the destruction of pancreatic beta cells, causing a shortage of insulin secretion. The development of T1D is influenced by both genetic predisposition and environmental factors, such as vitamin D. This vitamin is known for its ability to regulate the immune system and has been associated with a decreased risk of T1D. However, the specific ways in which vitamin D affects immune regulation and the preservation of beta cells in T1D are not yet fully understood. Gaining a better understanding of these interactions is essential for identifying potential targets for preventing and treating T1D. Methods The analysis focused on two Gene Expression Omnibus (GEO) datasets, namely, GSE55098 and GSE50012, to detect differentially expressed genes (DEGs). Enrichr (Ma'ayan Laboratory, New York, NY) was used to perform enrichment analysis for the Gene Ontology (GO) biological process and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The Search Tool for the Retrieval of Interacting Genes 12.0 (STRING) database was used to generate a protein-protein interaction (PPI) network. The Cytoscape 3.10.1 (Cytoscape Team, San Diego, CA) was used to analyze the PPI network and discover the hub genes. Results The DEGs in both datasets were identified using the GEO2R tool, with a particular focus on genes exhibiting contrasting regulations. Enrichment analysis unveiled the participation of these oppositely regulated DEGs in processes relevant to the immune system. Cytoscape analysis of the PPI network revealed five hub genes, MNDA, LILRB2, FPR2, HCK, and FCGR2A, suggesting their potential role in the pathogenesis of T1D and the response to vitamin D. Conclusion The study elucidates the complex interaction between vitamin D metabolism and immune regulation in T1D. The identified hub genes provide important knowledge on the molecular pathways that underlie T1D and have the potential to be targeted for therapeutic intervention. This research underscores the importance of vitamin D in the immune system's modulation and its impact on T1D development.

Comprehensive molecular characterization of collecting duct carcinoma for therapeutic vulnerability

Collecting duct carcinoma (CDC) is an aggressive rare subtype of kidney cancer with unmet clinical needs. Little is known about its underlying molecular alterations and etiology, primarily due to its rarity, and lack of preclinical models. This study aims to comprehensively characterize molecular alterations in CDC and identify its therapeutic vulnerabilities. Through whole-exome and transcriptome sequencing, we identified KRAS hotspot mutations (G12A/D/V) in 3/13 (23%) of the patients, in addition to known TP53, NF2 mutations. 3/13 (23%) patients carried a mutational signature (SBS22) caused by aristolochic acid (AA) exposures, known to be more prevalent in Asia, highlighting a geologically specific disease etiology. We further discovered that cell cycle-related pathways were the most predominantly dysregulated pathways. Our drug screening with our newly established CDC preclinical models identified a CDK9 inhibitor LDC000067 that specifically inhibited CDC tumor growth and prolonged survival. Our study not only improved our understanding of oncogenic molecular alterations of Asian CDC, but also identified cell-cycle machinery as a therapeutic vulnerability, laying the foundation for clinical trials to treat patients with such aggressive cancer.

Skeletal muscle hypertrophy: cell growth is cell growth

Roberts et al. have provided an insightful counterpoint to our review article on the utility of the synergist ablation model. The purpose of this review is to provide some further dialogue regarding the strengths and weaknesses of the synergist ablation model. Specifically, we highlight that the robustness of the model overshadows surgical limitations. We also compare the transcriptomic responses to synergist ablation in mice and resistance exercise in humans to identify common pathways. We conclude that "cell growth is cell growth" and that the mechanisms available to cells to accumulate biomass and increase in size are similar across cell types and independent of the rate of growth.

Chemical profiling of Anthriscus cerefolium (L.) Hoffm., biological potential of the herbal extract, molecular modeling and KEGG pathway analysis

This study was designed to investigate chemical composition and biological activities of the Anthriscus cerefolium methanolic extract. Chemical characterization of the extracts was performed by LC-HRMS/MS analysis. Antimicrobial activities of the extract were investigated on six bacteria and eight fungi while antioxidant activity was assessed by six different assays. Anti-enzymatic activity of the methanolic extract was tested on five enzymes associated with therapy of neurodegenerative diseases and diabetes mellitus type 2. Cytotoxic properties of the extract were tested on human immortalized keratinocytes (HaCaT) and tumor cell lines (SiHa, MCF7, HepG2). Anti-inflammatory activity of the extract was assessed on bacteria mediated inflammation model using HaCaT cell line. Molecular docking studies of enzymes and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis were performed. The results showed that the obtained extract was rich in phenolic compounds (a total of seventy-two were identified), with malonyl-1,4-O-dicaffeoylquinic acid and 3,5-O-dicaffeoylquinic acid dominating in the sample. The extract expressed antimicrobial, antioxidant, anti-enzymatic, cytotoxic and anti-inflammatory properties. The identified compounds demonstrated strong binding to the acetylcholinesterase (AChE) and to a lesser extent, to the butyrylcholinesterase (BChE), glucosidase, amylase, and modestly, to tyrosinase. KEGG pathway analysis has shown that the certain phenolic compounds may be related to anti-tumor, anti-inflammatory and anti-microbial activities of the extract. The data obtained suggest that phenolic compounds of the extract and their mixtures should be considered for future research as ingredients in pharmaceutical and nutraceutical formulations.

Gene expression analysis identifies hub genes and pathways distinguishing fatal from survivor outcomes of Ebola virus disease

The Ebola virus poses a severe public health threat, yet understanding factors influencing disease outcomes remains incomplete. Our study aimed to identify critical pathways and hub genes associated with fatal and survivor Ebola disease outcomes. We analyzed differentially expressed hub genes (DEGs) between groups with fatal and survival outcomes, as well as a healthy control group. We conducted additional analysis to determine the functions and pathways associated with these DEGs. We found 13,198 DEGs in the fatal and 12,039 DEGs in the survival group compared to healthy controls, and 1873 DEGs in the acute fatal and survivor groups comparison. Upregulated DEGs in the comparison between the acute fatal and survivor groups were linked to ECM receptor interaction, complement and coagulation cascades, and PI3K-Akt signaling. Upregulated hub genes identified from the acute fatal and survivor comparison (FGB, C1QA, SERPINF2, PLAT, C9, SERPINE1, F3, VWF) were enriched in complement and coagulation cascades; the downregulated hub genes (IL1B, 1L17RE, XCL1, CXCL6, CCL4, CD8A, CD8B, CD3D) were associated with immune cell processes. Hub genes CCL2 and F2 were unique to fatal outcomes, while CXCL1, HIST1H4F, and IL1A were upregulated hub genes unique to survival outcomes compared to healthy controls. Our results demonstrate for the first time the association of EVD outcomes to specific hub genes and their associated pathways and biological processes. The identified hub genes and pathways could help better elucidate Ebola disease pathogenesis and contribute to the development of targeted interventions and personalized treatment for distinct EVD outcomes.

Construction of an immune gene expression meta signature to assess the prognostic risk of colorectal cancer patients

Despite recent advancements in colorectal cancer (CRC) treatment, particularly with the introduction of immunotherapy and checkpoint inhibitors, the efficacy of these therapies remains limited to a subset of patients. To address this challenge, our study aimed to develop a prognostic biomarker based on immune-related genes to predict better outcomes in CRC patients and aid in treatment decision-making. We comprehensively analysed immune gene expression signatures associated with CRC prognosis to construct an immune meta-signature with prognostic potential. Utilising data from The Cancer Genome Atlas (TCGA), we employed Cox regression to identify immune-related genes with prognostic significance from multiple studies. Subsequently, we compared the expression levels of immune genes, levels of immune cell infiltration, and various immune-related molecules between high-risk and low-risk patient groups. Functional analysis using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses provided insights into the biological pathways associated with the identified prognostic genes. Finally, we validated our findings using a separate CRC cohort from the Gene Expression Omnibus (GEO). Integration of the prognostic genes revealed significant disparities in survival outcomes. Differential expression analysis identified a set of immune-associated genes, which were further refined using LASSO penalisation and Cox regression. Univariate Cox regression analyses confirmed the autonomy of the gene signature as a prognostic indicator for CRC patient survival. Our risk prediction model effectively stratified CRC patients based on their prognosis, with the high-risk group showing enrichment in pro-oncogenic terms and pathways. Immune infiltration analysis revealed an augmented presence of certain immunosuppressive subsets in the high-risk group. Finally, we validated the performance of our prognostic model by applying the risk score equation to a different CRC patient dataset, confirming its prognostic potential in this new cohort. Overall, our study presents a novel immune-related gene signature with promising implications for predicting cancer progression and prognosis, thereby enabling more personalised management strategies for CRC patients.

Collagen-Heparin-FGF2-VEGF Scaffolds Induce a Regenerative Gene Expression Profile in a Fetal Sheep Wound Model

BACKGROUND

The developmental abnormality spina bifida is hallmarked by missing tissues (e.g. skin) and exposure of the spinal cord to the amniotic fluid, which can negatively impact neurological development. Surgical closure of the skin in utero limits neurological damage, but in large defects this results in scarring and contractures. Stimulating skin regeneration in utero would greatly benefit treatment outcome. Previously, we demonstrated that a porous type I collagen (COL) scaffold, functionalized with heparin (HEP), fibroblast growth factor 2 (FGF2) and vascular endothelial growth factor (VEGF) (COL-HEP/GF) improved pre- and postnatal skin regeneration in a fetal sheep full thickness wound model. In this study we uncover the early events associated with enhanced skin regeneration.

METHODS

We investigated the gene expression profiles of healing fetal skin wounds two weeks after implantation of the COL(-HEP/GF) scaffolds. Using laser dissection and microarrays, differentially expressed genes (DEG) were identified in the epidermis and dermis between untreated wounds, COL-treated wounds and wounds treated with COL-HEP/GF. Biological processes were identified using gene enrichment analysis and DEG were clustered using protein-protein-interaction networks.

RESULTS

COL-HEP/GF influences various interesting biological processes involved in wound healing. Although the changes were modest, using protein-protein-interaction networks we identified a variety of clustered genes that indicate COL-HEP/GF induces a tight but subtle control over cell signaling and extracellular matrix organization.

CONCLUSION

These data offer a novel perspective on the key processes involved in (fetal) wound healing, where a targeted and early interference during wound healing can result in long-term enhanced effects on skin regeneration.

Post-operative Crohn's Disease Recurrence and Infectious Complications: A Transcriptomic Analysis

BACKGROUND

Crohn's disease (CD) is a chronic inflammatory condition affecting the gastrointestinal tract, characterized by complications such as strictures, fistulas, and neoplasia. Despite medical advancements, a significant number of patients with Crohn's disease require surgery, and many experience post-operative complications and recurrence. Previous studies have analyzed gene expression to study recurrence and post-operative complications independently. This study aims to identify overlapping differentially expressed genes and pathways for recurrence and post-operative complications.

METHODS

A dataset including 45 patients with Crohn's disease, including gene expression from ileum and colon tissue, endoscopic recurrence, and intra-abdominal septic complications was analyzed. Gene set enrichment analysis was used to identify gene pathways associated with the outcomes. Finally, a multi-variable logistic regression model was created to assess whether gene pathways were independently associated with both outcomes.

RESULTS

In ileum tissue, several inflammatory pathways, including interferon alpha and gamma response were upregulated in patients with endoscopic recurrence and intra-abdominal septic complications. In addition, there was upregulation of the epithelial mesenchymal transition pathway. In colon tissue, metabolic processes, such as myogenesis and oxidative phosphorylation were downregulated in both outcomes. In a multivariate model, downregulation of myogenesis in colon tissue was significantly associated with both endoscopic recurrence and intra-abdominal septic complications.

CONCLUSION

These findings shed light on the underlying biology of these outcomes and suggest potential biomarkers or therapeutic targets to reduce their occurrence. Further validation and multi-institutional studies are warranted to confirm these results and improve post-operative outcomes for patients with Crohn's disease.

PIP-seq identifies novel heterogeneous lung innate lymphocyte population activation after combustion product exposure

Innate lymphoid cells (ILCs) are a heterogeneous population that play diverse roles in airway inflammation after exposure to allergens and infections. However, how ILCs respond after exposure to environmental toxins is not well understood. Here we show a novel method for studying the heterogeneity of rare lung ILC populations by magnetic enrichment for lung ILCs followed by particle-templated instant partition sequencing (PIP-seq). Using this method, we were able to identify novel group 1 and group 2 ILC subsets that exist after exposure to both fungal allergen and burn pit-related constituents (BPC) that include dioxin, aromatic hydrocarbon, and particulate matter. Toxin exposure in combination with fungal allergen induced activation of specific ILC1/NK and ILC2 populations as well as promoted neutrophilic lung inflammation. Oxidative stress pathways and downregulation of specific ribosomal protein genes (Rpl41 and Rps19) implicated in anti-inflammatory responses were present after BPC exposure. Increased IFNγ expression and other pro-neutrophilic mediator transcripts were increased in BPC-stimulated lung innate lymphoid cells. Further, the addition of BPC induced Hspa8 (encodes HSC70) and aryl hydrocarbon transcription factor activity across multiple lung ILC subsets. Overall, using an airway disease model that develops after occupational and environmental exposures, we demonstrate an effective method to better understand heterogenous ILC subset activation.

Structural and systems characterization of phosphorylation on metabolic enzymes identifies sex-specific metabolic reprogramming in obesity

Coordination of adaptive metabolism through cellular signaling networks and metabolic response is essential for balanced flow of energy and homeostasis. Post-translational modifications such as phosphorylation offer a rapid, efficient, and dynamic mechanism to regulate metabolic networks. Although numerous phosphorylation sites have been identified on metabolic enzymes, much remains unknown about their contribution to enzyme function and systemic metabolism. In this study, we stratify phosphorylation sites on metabolic enzymes based on their location with respect to functional and dimerization domains. Our analysis reveals that the majority of published phosphosites are on oxidoreductases, with particular enrichment of phosphotyrosine (pY) sites in proximity to binding domains for substrates, cofactors, active sites, or dimer interfaces. We identify phosphosites altered in obesity using a high fat diet (HFD) induced obesity model coupled to multiomics, and interrogate the functional impact of pY on hepatic metabolism. HFD induced dysregulation of redox homeostasis and reductive metabolism at the phosphoproteome and metabolome level in a sex-specific manner, which was reversed by supplementing with the antioxidant butylated hydroxyanisole (BHA). Partial least squares regression (PLSR) analysis identified pY sites that predict HFD or BHA induced changes of redox metabolites. We characterize predictive pY sites on glutathione S-transferase pi 1 (GSTP1), isocitrate dehydrogenase 1 (IDH1), and uridine monophosphate synthase (UMPS) using CRISPRi-rescue and stable isotope tracing. Our analysis revealed that sites on GSTP1 and UMPS inhibit enzyme activity while the pY site on IDH1 induces activity to promote reductive carboxylation. Overall, our approach provides insight into the convergence points where cellular signaling fine-tunes metabolism.

Modulation of blood-tumor barrier transcriptional programs improves intra-tumoral drug delivery and potentiates chemotherapy in GBM

Glioblastoma (GBM) is the most common malignant primary brain tumor. GBM has an extremely poor prognosis and new treatments are badly needed. Efficient drug delivery to GBM is a major obstacle as the blood-brain barrier (BBB) prevents passage of the majority of cancer drugs into the brain. It is also recognized that the blood-brain tumor barrier (BTB) in the growing tumor represents a challenge. The BTB is heterogeneous and poorly characterized, but similar to the BBB it can prevent therapeutics from reaching effective intra-tumoral doses, dramatically hindering their potential. Here, we identified a 12-gene signature associated with the BTB, with functions related to vasculature development, morphogenesis and cell migration. We identified CDH5 as a core molecule in this set and confirmed its over-expression in GBM vasculature using spatial transcriptomics of GBM patient specimens. We found that the indirubin-derivative, 6-bromoindirubin acetoxime (BIA), could downregulate CDH5 and other BTB signature genes, causing endothelial barrier disruption in endothelial monolayers and BBB 3D spheroids in vitro. Treatment of tumor-bearing mice with BIA enabled increased intra-tumoral accumulation of the BBB non-penetrant chemotherapeutic drug cisplatin and potentiated cisplatin-mediated DNA damage by targeting DNA repair pathways. Finally, using an injectable BIA nanoparticle formulation, PPRX-1701, we significantly improved the efficacy of cisplatin in patient-derived GBM xenograms and prolonged their survival. Overall, our work reveals potential targets at the BTB for improved chemotherapy delivery and the bifunctional properties of BIA as a BTB modulator and potentiator of chemotherapy, supporting its further development.

TFEB/3 Govern Repair Schwann Cell Generation and Function Following Peripheral Nerve Injury

TFEB and TFE3 (TFEB/3), key regulators of lysosomal biogenesis and autophagy, play diverse roles depending on cell type. This study highlights a hitherto unrecognized role of TFEB/3 crucial for peripheral nerve repair. Specifically, they promote the generation of progenitor-like repair Schwann cells after axonal injury. In Schwann cell-specific TFEB/3 double knock-out mice of either sex, the TFEB/3 loss disrupts the transcriptomic reprogramming that is essential for the formation of repair Schwann cells. Consequently, mutant mice fail to populate the injured nerve with repair Schwann cells and exhibit defects in axon regrowth, target reinnervation, and functional recovery. TFEB/3 deficiency inhibits the expression of injury-responsive repair Schwann cell genes, despite the continued expression of c-jun, a previously identified regulator of repair Schwann cell function. TFEB/3 binding motifs are enriched in the enhancer regions of injury-responsive genes, suggesting their role in repair gene activation. Autophagy-dependent myelin breakdown is not impaired despite TFEB/3 deficiency. These findings underscore a unique role of TFEB/3 in adult Schwann cells that is required for proper peripheral nerve regeneration.

Clonal evolution of the 3D chromatin landscape in patients with relapsed pediatric B-cell acute lymphoblastic leukemia

Relapsed pediatric B-cell acute lymphoblastic leukemia (B-ALL) remains one of the leading causes of cancer mortality in children. We performed Hi-C, ATAC-seq, and RNA-seq on 12 matched diagnosis/relapse pediatric leukemia specimens to uncover dynamic structural variants (SVs) and 3D chromatin rewiring that may contribute to relapse. While translocations are assumed to occur early in leukemogenesis and be maintained throughout progression, we discovered novel, dynamic translocations and confirmed several fusion transcripts, suggesting functional and therapeutic relevance. Genome-wide chromatin remodeling was observed at all organizational levels: A/B compartments, TAD interactivity, and chromatin loops, including some loci shared by 25% of patients. Shared changes were found to drive the expression of genes/pathways previously implicated in resistance as well as novel therapeutic candidates, two of which (ATXN1 and MN1) we functionally validated. Overall, these results demonstrate chromatin reorganization under the selective pressure of therapy and offer the potential for discovery of novel therapeutic interventions.

STAG2 mutations reshape the cohesin-structured spatial chromatin architecture to drive gene regulation in acute myeloid leukemia

Cohesin shapes the chromatin architecture, including enhancer-promoter interactions. Its components, especially STAG2, but not its paralog STAG1, are frequently mutated in myeloid malignancies. To elucidate the underlying mechanisms of leukemogenesis, we comprehensively characterized genetic, transcriptional, and chromatin conformational changes in acute myeloid leukemia (AML) patient samples. Specific loci displayed altered cohesin occupancy, gene expression, and local chromatin activation, which were not compensated by the remaining STAG1-cohesin. These changes could be linked to disrupted spatial chromatin looping in cohesin-mutated AMLs. Complementary depletion of STAG2 or STAG1 in primary human hematopoietic progenitors (HSPCs) revealed effects resembling STAG2-mutant AML-specific changes following STAG2 knockdown, not invoked by the depletion of STAG1. STAG2-deficient HSPCs displayed impaired differentiation capacity and maintained HSPC-like gene expression. This work establishes STAG2 as a key regulator of chromatin contacts, gene expression, and differentiation in the hematopoietic system and identifies candidate target genes that may be implicated in human leukemogenesis.

ΔNp63 bookmarks and creates an accessible epigenetic environment for TGFβ-induced cancer cell stemness and invasiveness

BACKGROUND

p63 is a transcription factor with intrinsic pioneer factor activity and pleiotropic functions. Transforming growth factor β (TGFβ) signaling via activation and cooperative action of canonical, SMAD, and non-canonical, MAP-kinase (MAPK) pathways, elicits both anti- and pro-tumorigenic properties, including cell stemness and invasiveness. TGFβ activates the ΔNp63 transcriptional program in cancer cells; however, the link between TGFβ and p63 in unmasking the epigenetic landscape during tumor progression allowing chromatin accessibility and gene transcription, is not yet reported.

METHODS

Small molecule inhibitors, including protein kinase inhibitors and RNA-silencing, provided loss of function analyses. Sphere formation assays in cancer cells, chromatin immunoprecipitation and mRNA expression assays were utilized in order to gain mechanistic evidence. Mass spectrometry analysis coupled to co-immunoprecipitation assays revealed novel p63 interactors and their involvement in p63-dependent transcription.

RESULTS

The sphere-forming capacity of breast cancer cells was enhanced upon TGFβ stimulation and significantly decreased upon ΔNp63 depletion. Activation of TGFβ signaling via p38 MAPK signaling induced ΔNp63 phosphorylation at Ser 66/68 resulting in stabilized ΔNp63 protein with enhanced DNA binding properties. TGFβ stimulation altered the ratio of H3K27ac and H3K27me3 histone modification marks, pointing towards higher H3K27ac and increased p300 acetyltransferase recruitment to chromatin. By silencing the expression of ΔNp63, the TGFβ effect on chromatin remodeling was abrogated. Inhibition of H3K27me3, revealed the important role of TGFβ as the upstream signal for guiding ΔNp63 to the TGFβ/SMAD gene loci, as well as the indispensable role of ΔNp63 in recruiting histone modifying enzymes, such as p300, to these genomic regions, regulating chromatin accessibility and gene transcription. Mechanistically, TGFβ through SMAD activation induced dissociation of ΔNp63 from NURD or NCOR/SMRT histone deacetylation complexes, while promoted the assembly of ΔNp63-p300 complexes, affecting the levels of histone acetylation and the outcome of ΔNp63-dependent transcription.

CONCLUSIONS

ΔNp63, phosphorylated and recruited by TGFβ to the TGFβ/SMAD/ΔNp63 gene loci, promotes chromatin accessibility and transcription of target genes related to stemness and cell invasion.

Deciphering Pre-existing and Induced 3D Genome Architecture Changes involved in Constricted Melanoma Migration

Metastatic cancer cells traverse constricted spaces that exert forces on their nucleus and the genomic contents within. Cancerous tumors are highly heterogeneous and not all cells within them can achieve such a feat. Here, we investigated what initial genome architecture characteristics favor the constricted migratory ability of cancer cells and which arise only after passage through multiple constrictions. We identified a cell surface protein (ITGB4) whose expression correlates with increased initial constricted migration ability in human melanoma A375 cells. Sorting out this subpopulation allowed us to identify cellular and nuclear features that pre-exist and favor migration, as well as alterations that only appear after cells have passed through constrictions. We identified specific genomic regions that experienced altered genome spatial compartment profiles only after constricted migration. Our study reveals 3D genome structure contributions to both selection and induction mechanisms of cell fate change during cancer metastasis.

Effects of Inonotus obliquus on ameliorating podocyte injury in ORG mice through TNF pathway and prediction of active compounds

Background

Podocyte injury is a common pathologic mechanism in diabetic kidney disease (DKD) and obesity-related glomerulopathy (ORG). Our previous study confirmed that Inonotus obliquus (IO) improved podocyte injury on DKD rats. The current study explored the pharmacological effects, related mechanisms and possible active components of IO on ORG mice.

Methods

Firstly, by combining ultra-high performance liquid chromatography tandem mass spectrometry analysis (UPLC-Q-TOF-MS) with network pharmacology to construct the human protein-protein interaction mechanism and enrich the pathway, which led to discover the crucial mechanism of IO against ORG. Then, ORG mice were established by high-fat diet and biochemical assays, histopathology, and Western blot were used to explore the effects of IO on obesity and podocyte injury. Finally, network pharmacology-based findings were confirmed by immunohistochemistry. The compositions of IO absorbed in mice plasma were analyzed by UPLC-Q-TOF-MS and molecular docking was used to predict the possible active compounds.

Results

The network pharmacology result suggested that IO alleviated the inflammatory response of ORG by modulating TNF signal. The 20-week in vivo experiment confirmed that IO improved glomerular hypertrophy, podocyte injury under electron microscopy, renal nephrin, synaptopodin, TNF-α and IL-6 expressions with Western blotting and immunohistochemical staining. Other indicators of ORG such as body weight, kidney weight, serum total cholesterol, liver triglyceride also improved by IO intervention. The components analysis showed that triterpenoids, including inoterpene F and trametenolic acid, might be the pharmacodynamic basis.

Conclusion

The research based on UPLC-Q-TOF-MS analysis, network pharmacology and in vivo experiment suggested that the amelioration of IO on podocyte injury in ORG mice via its modulation on TNF signal. Triterpenoids were predicated as acting components.

Whole exome sequencing identifies genetic markers of enterovirus susceptibility in East Asians

Introduction

Following acute enterovirus (EV) infection, outcomes vary based on factors like the immune response, viral cell entry receptor expression levels, tissue tropism, and genetic factors of both the host and virus. While most individuals exhibit mild, self-limited symptoms, others may suffer severe complications or prolonged infections that can lead to autoimmune disorders.

Methods

To elucidate host responses to EV infection, we performed whole exome sequencing on blood samples from both infected and uninfected individuals. Our initial focus was on genes encoding EV entry receptors-PSGL-1, SCARB2, and ANAXA2 for EV-A71, and CD155 for poliovirus-and on host genes ACBD3 and PI4KΒ, crucial for EV replication.

Results

Although no specific genetic variants directly associated with EV infection were identified, we discovered 118 variants across 116 genes enriched in East Asian populations through multi-layered variant filtering. These variants were further analyzed for their potential impacts on organs, biological processes, and molecular pathways. Phenome-wide association studies were conducted to refine our understanding of their contributions to EV infection susceptibility.

Discussion

Our findings aim to develop a predictive panel based on these 118 variants, which could help susceptible individuals during EV outbreaks, guiding targeted clinical interventions and preventative strategies.

Genomic Insights into Idiopathic Granulomatous Mastitis through Whole-Exome Sequencing: A Case Report of Eight Patients

Idiopathic granulomatous mastitis (IGM) is a rare condition characterised by chronic inflammation and granuloma formation in the breast. The aetiology of IGM is unclear. By focusing on the protein-coding regions of the genome, where most disease-related mutations often occur, whole-exome sequencing (WES) is a powerful approach for investigating rare and complex conditions, like IGM. We report WES results on paired blood and tissue samples from eight IGM patients. Samples were processed using standard genomic protocols. Somatic variants were called with two analytical pipelines: nf-core/sarek with Strelka2 and GATK4 with Mutect2. Our WES study of eight patients did not find evidence supporting a clear genetic component. The discrepancies between variant calling algorithms, along with the considerable genetic heterogeneity observed amongst the eight IGM cases, indicate that common genetic drivers are not readily identifiable. With only three genes, CHIT1, CEP170, and CTR9, recurrently altering in multiple cases, the genetic basis of IGM remains uncertain. The absence of validation for somatic variants by Sanger sequencing raises further questions about the role of genetic mutations in the disease. Other potential contributors to the disease should be explored.

Mdm2 requires Sprouty4 to regulate focal adhesion formation and metastasis independent of p53

Although the E3 ligase Mdm2 and its homologue and binding partner MdmX are the major regulators of the p53 tumor suppressor protein, it is now evident that Mdm2 and MdmX have multiple functions that do not involve p53. As one example, it is known that Mdm2 can regulate cell migration, although mechanistic insight into this function is still lacking. Here we show in cells lacking p53 expression that knockdown of Mdm2 or MdmX, as well as pharmacological inhibition of the Mdm2/MdmX complex, not only reduces cell migration and invasion, but also impairs cell spreading and focal adhesion formation. In addition, Mdm2 knockdown decreases metastasis in vivo. Interestingly, Mdm2 downregulates the expression of Sprouty4, which is required for the Mdm2 mediated effects on cell migration, focal adhesion formation and metastasis. Further, our findings indicate that Mdm2 dampening of Sprouty4 is a prerequisite for maintaining RhoA levels in the cancer cells that we have studied. Taken together we describe a molecular mechanism whereby the Mdm2/MdmX complex through Sprouty4 regulates cellular processes leading to increase metastatic capability independently of p53.

In Vitro Investigation of the Anti-Fibrotic Effects of 1-Phenyl-2-Pentanol, Identified from Moringa oleifera Lam., on Hepatic Stellate Cells

Liver fibrosis, characterized by excessive extracellular matrix deposition, is driven by activated hepatic stellate cells (HSCs). Due to the limited availability of anti-fibrotic drugs, the research on therapeutic agents continues. Here we have investigated Moringa oleifera Lam. (MO), known for its various bioactive properties, for anti-fibrotic effects. This study has focused on 1-phenyl-2-pentanol (1-PHE), a compound derived from MO leaves, and its effects on LX-2 human hepatic stellate cell activation. TGF-β1-stimulated LX-2 cells were treated with MO extract or 1-PHE, and the changes in liver fibrosis markers were assessed at both gene and protein levels. Proteomic analysis and molecular docking were employed to identify potential protein targets and signaling pathways affected by 1-PHE. Treatment with 1-PHE downregulated fibrosis markers, including collagen type I alpha 1 chain (COL1A1), collagen type IV alpha 1 chain (COL4A1), mothers against decapentaplegic homologs 2 and 3 (SMAD2/3), and matrix metalloproteinase-2 (MMP2), and reduced the secretion of matrix metalloproteinase-9 (MMP-9). Proteomic analysis data showed that 1-PHE modulates the Wnt/β-catenin pathway, providing a possible mechanism for its effects. Our results suggest that 1-PHE inhibits the TGF-β1 and Wnt/β-catenin signaling pathways and HSC activation, indicating its potential as an anti-liver-fibrosis agent.

A single-cell transcriptomic map of the developing Atoh1 lineage identifies neural fate decisions and neuronal diversity in the hindbrain

Proneural transcription factors establish molecular cascades to orchestrate neuronal diversity. One such transcription factor, Atonal homolog 1 (Atoh1), gives rise to cerebellar excitatory neurons and over 30 distinct nuclei in the brainstem critical for hearing, breathing, and balance. Although Atoh1 lineage neurons have been qualitatively described, the transcriptional programs that drive their fate decisions and the full extent of their diversity remain unknown. Here, we analyzed single-cell RNA sequencing and ATOH1 DNA binding in Atoh1 lineage neurons of the developing mouse hindbrain. This high-resolution dataset identified markers for specific brainstem nuclei and demonstrated that transcriptionally heterogeneous progenitors require ATOH1 for proper migration. Moreover, we identified a sizable population of proliferating unipolar brush cell progenitors in the mouse Atoh1 lineage, previously described in humans as the origin of one medulloblastoma subtype. Collectively, our data provide insights into the developing mouse hindbrain and markers for functional assessment of understudied neuronal populations.

Two Different Immune Profiles Are Identified in Sentinel Lymph Nodes of Early-Stage Breast Cancer

In the management of early-stage breast cancer (BC), lymph nodes (LNs) are typically characterised using the One-Step Nucleic Acid Amplification (OSNA) assay, a standard procedure for assessing subclinical metastasis in sentinel LNs (SLNs). The pivotal role of LNs in coordinating the immune response against BC is often overlooked. Our aim was to improve prognostic information provided by the OSNA assay and explore immune-related gene signatures in SLNs. The expression of an immune gene panel was analysed in SLNs from 32 patients with Luminal A early-stage BC (cT1-T2 N0). Using an unsupervised approach based on these expression values, this study identified two clusters, regardless of the SLN invasion: one evidencing an adaptive anti-tumoral immune response, characterised by an increase in naive B cells, follicular T helper cells, and activated NK cells; and another with a more undifferentiated response, with an increase in the activated-to-resting dendritic cells (DCs) ratio. Through a protein-protein interaction (PPI) network, we identified seven immunoregulatory hub genes: CD80, CD40, TNF, FCGR3A, CD163, FCGR3B, and CCR2. This study shows that, in Luminal A early-stage BC, SLNs gene expression studies enable the identification of distinct immune profiles that may influence prognosis stratification and highlight key genes that could serve as potential targets for immunotherapy.

Multiomics2Targets identifies targets from cancer cohorts profiled with transcriptomics, proteomics, and phosphoproteomics

The availability of data from profiling of cancer patients with multiomics is rapidly increasing. However, integrative analysis of such data for personalized target identification is not trivial. Multiomics2Targets is a platform that enables users to upload transcriptomics, proteomics, and phosphoproteomics data matrices collected from the same cohort of cancer patients. After uploading the data, Multiomics2Targets produces a report that resembles a research publication. The uploaded matrices are processed, analyzed, and visualized using the tools Enrichr, KEA3, ChEA3, Expression2Kinases, and TargetRanger to identify and prioritize proteins, genes, and transcripts as potential targets. Figures and tables, as well as descriptions of the methods and results, are automatically generated. Reports include an abstract, introduction, methods, results, discussion, conclusions, and references and are exportable as citable PDFs and Jupyter Notebooks. Multiomics2Targets is applied to analyze version 3 of the Clinical Proteomic Tumor Analysis Consortium (CPTAC3) pan-cancer cohort, identifying potential targets for each CPTAC3 cancer subtype. Multiomics2Targets is available from https://multiomics2targets.maayanlab.cloud/.

Transcriptional Response to Tick-Borne Flavivirus Infection in Neurons, Astrocytes and Microglia In Vivo and In Vitro

Tick-borne encephalitis virus (TBEV) is a neurotropic member of the genus Orthoflavivirus (former Flavivirus) and is of significant health concern in Europe and Asia. TBEV pathogenesis may occur directly via virus-induced damage to neurons or through immunopathology due to excessive inflammation. While primary cells isolated from the host can be used to study the immune response to TBEV, it is still unclear how well these reflect the immune response elicited in vivo. Here, we compared the transcriptional response to TBEV and the less pathogenic tick-borne flavivirus, Langat virus (LGTV), in primary monocultures of neurons, astrocytes and microglia in vitro, with the transcriptional response in vivo captured by single-nuclei RNA sequencing (snRNA-seq) of a whole mouse cortex. We detected similar transcriptional changes induced by both LGTV and TBEV infection in vitro, with the lower response to LGTV likely resulting from slower viral kinetics. Gene set enrichment analysis showed a stronger transcriptional response in vivo than in vitro for astrocytes and microglia, with a limited overlap mainly dominated by interferon signaling. Together, this adds to our understanding of neurotropic flavivirus pathogenesis and the strengths and limitations of available model systems.

The interferon/STAT1 signaling axis is a common feature of tumor-initiating cells in breast cancer

A tumor cell subpopulation of tumor-initiating cells (TIC), or "cancer stem cells", are associated with therapeutic resistance, as well as both local and distant recurrence. Enriched populations of TIC are identified by markers including aldehyde dehydrogenase (ALDH1) activity, the cell surface marker combination CD44 + /CD24 - , or fluorescent reporters for signaling pathways that regulate TIC function. We showed previously that S ignal T ransducer and A ctivator of T ranscription (STAT)-mediated transcription allows enrichment for TIC in claudin-low models of human triple-negative breast cancer using a STAT-responsive reporter. However, the molecular phenotypes of STAT TIC are not well understood, and there is no existing method to lineage-trace TIC as they undergo cell state changes. Using a new STAT-responsive lineage-tracing (LT) system in conjunction with our original reporter, we enriched for cells with enhanced mammosphere-forming potential in some, but not all, basal-like triple-negative breast cancer (TNBC) xenograft models (TNBC) indicating TIC-related and TIC-independent functions for STAT signaling. Single-cell RNA sequencing (scRNAseq) of reporter-tagged xenografts and clinical samples identified a common interferon (IFN)/STAT1-associated transcriptional state, previously linked to inflammation and macrophage differentiation, in TIC. Surprisingly, most of the genes we identified are not present in previously published TIC signatures derived using bulk RNA sequencing. Finally, we demonstrated that bone marrow stromal cell antigen 2 (BST2), is a cell surface marker of this state, and that it functionally regulates TIC frequency. These results suggest TIC may exploit the IFN/STAT1 signaling axis to promote their activity, and that targeting this pathway may help eliminate TIC.

Significance

TIC differentially express interferon response genes, which were not previously reported in bulk RNA sequencing-derived TIC signatures, highlighting the importance of coupling single-cell transcriptomics with enrichment to derive TIC signatures.

Missense Mutations in Myc Box I Influence Nucleocytoplasmic Transport to Promote Leukemogenesis

PURPOSE

Somatic missense mutations in the phosphodegron domain of the MYC gene (MYC Box I or MBI) are detected in the dominant clones of a subset of patients with acute myeloid leukemia (AML), but the mechanisms by which they contribute to AML are unknown.

EXPERIMENTAL DESIGN

To investigate the effects of MBI MYC mutations on hematopoietic cells, we employed a multi-omic approach to systematically compare the cellular and molecular consequences of expressing oncogenic doses of wild type, threonine-58 and proline-59 mutant MYC proteins in hematopoietic cells, and we developed a knockin mouse harboring the germline MBI mutation p.T58N in the Myc gene.

RESULTS

Both wild-type and MBI mutant MYC proteins promote self-renewal programs and expand highly selected subpopulations of progenitor cells in the bone marrow. Compared with their wild-type counterparts, mutant cells display decreased cell death and accelerated leukemogenesis in vivo, changes that are recapitulated in the transcriptomes of human AML-bearing MYC mutations. The mutant phenotypes feature decreased stability and translation of mRNAs encoding proapoptotic and immune-regulatory genes, increased translation of RNA binding proteins and nuclear export machinery, and distinct nucleocytoplasmic RNA profiles. MBI MYC mutant proteins also show a higher propensity to aggregate in perinuclear regions and cytoplasm. Like the overexpression model, heterozygous p.T58N knockin mice displayed similar changes in subcellular MYC localization, progenitor expansion, transcriptional signatures, and develop hematopoietic tumors.

CONCLUSIONS

This study uncovers that MBI MYC mutations alter RNA nucleocytoplasmic transport mechanisms to contribute to the development of hematopoietic malignancies.

The Loss of PPARγ Expression and Signaling Is a Key Feature of Cutaneous Actinic Disease and Squamous Cell Carcinoma: Association with Tumor Stromal Inflammation

Given the importance of peroxisome proliferator-activated receptor (PPAR)-gamma in epidermal inflammation and carcinogenesis, we analyzed the transcriptomic changes observed in epidermal PPARγ-deficient mice (Pparg-/-epi). A gene set enrichment analysis revealed a close association with epithelial malignancy, inflammatory cell chemotaxis, and cell survival. Single-cell sequencing of Pparg-/-epi mice verified changes to the stromal compartment, including increased inflammatory cell infiltrates, particularly neutrophils, and an increase in fibroblasts expressing myofibroblast marker genes. A comparison of transcriptomic data from Pparg-/-epi and publicly available human and/or mouse actinic keratoses (AKs) and cutaneous squamous cell carcinomas (SCCs) revealed a strong correlation between the datasets. Importantly, PPAR signaling was the top common inhibited canonical pathway in AKs and SCCs. Both AKs and SCCs also had significantly reduced PPARG expression and PPARγ activity z-scores. Smaller reductions in PPARA expression and PPARα activity and increased PPARD expression but reduced PPARδ activation were also observed. Reduced PPAR activity was also associated with reduced PPARα/RXRα activity, while LPS/IL1-mediated inhibition of RXR activity was significantly activated in the tumor datasets. Notably, these changes were not observed in normal sun-exposed skin relative to non-exposed skin. Finally, Ppara and Pparg were heavily expressed in sebocytes, while Ppard was highly expressed in myofibroblasts, suggesting that PPARδ has a role in myofibroblast differentiation. In conclusion, these data provide strong evidence that PPARγ and possibly PPARα represent key tumor suppressors by acting as master inhibitors of the inflammatory changes found in AKs and SCCs.

Kir6.1, a component of an ATP-sensitive potassium channel, regulates natural killer cell development

Involved in immunity and reproduction, natural killer (NK) cells offer opportunities to develop new immunotherapies to treat infections and cancer or to alleviate pregnancy complications. Most current strategies use cytokines or antibodies to enhance NK-cell function, but none use ion channel modulators, which are widely used in clinical practice to treat hypertension, diabetes, epilepsy, and other conditions. Little is known about ion channels in NK cells. We show that Kcnj8, which codes for the Kir6.1 subunit of a certain type of ATP-sensitive potassium (K ATP ) channel, is highly expressed in murine splenic and uterine NK cells compared to other K + channels previously identified in NK cells. Kcnj8 expression is highest in the most mature subset of splenic NK cells (CD27 - CD11b + ) and in NKG2A + or Ly49C/I + educated uterine NK cells. Using patch clamping, we show that a subset of NK cells expresses a current sensitive to the Kir6.1 blocker PNU-37883A. Kcnj8 does not participate in NK cell degranulation in response to tumor cells in vitro or rejection of tumor cells in vivo . Transcriptomics show that genes previously implicated in NK cell development are amongst those differentially expressed in CD27 - CD11b + NK cells deficient of Kcnj8 . Indeed, we found that mice with NK-cell specific Kcnj8 gene ablation have fewer CD11b + CD27 - and KLRG-1 + NK cells in the bone barrow and spleen. These results show that the K ATP subunit Kir6.1 has a key role in NK-cell development.

A practically efficient algorithm for identifying critical control proteins in directed probabilistic biological networks

Network controllability is unifying the traditional control theory with the structural network information rooted in many large-scale biological systems of interest, from intracellular networks in molecular biology to brain neuronal networks. In controllability approaches, the set of minimum driver nodes is not unique, and critical nodes are the most important control elements because they appear in all possible solution sets. On the other hand, a common but largely unexplored feature in network control approaches is the probabilistic failure of edges or the uncertainty in the determination of interactions between molecules. This is particularly true when directed probabilistic interactions are considered. Until now, no efficient algorithm existed to determine critical nodes in probabilistic directed networks. Here we present a probabilistic control model based on a minimum dominating set framework that integrates the probabilistic nature of directed edges between molecules and determines the critical control nodes that drive the entire network functionality. The proposed algorithm, combined with the developed mathematical tools, offers practical efficiency in determining critical control nodes in large probabilistic networks. The method is then applied to the human intracellular signal transduction network revealing that critical control nodes are associated with important biological features and perturbed sets of genes in human diseases, including SARS-CoV-2 target proteins and rare disorders. We believe that the proposed methodology can be useful to investigate multiple biological systems in which directed edges are probabilistic in nature, both in natural systems or when determined with large uncertainties in-silico.

CRISPRi-based screens in iAssembloids to elucidate neuron-glia interactions

The sheer complexity of the brain has complicated our ability to understand the cellular and molecular mechanisms underlying its function in health and disease. Genome-wide association studies have uncovered genetic variants associated with specific neurological phenotypes and diseases. In addition, single-cell transcriptomics have provided molecular descriptions of specific brain cell types and the changes they undergo during disease. Although these approaches provide a giant leap forward towards understanding how genetic variation can lead to functional changes in the brain, they do not establish molecular mechanisms. To address this need, we developed a 3D co-culture system termed iAssembloids (induced multi-lineage assembloids) that enables the rapid generation of homogenous neuron-glia spheroids. We characterize these iAssembloids with immunohistochemistry and single-cell transcriptomics and combine them with large-scale CRISPRi-based screens. In our first application, we ask how glial and neuronal cells interact to control neuronal death and survival. Our CRISPRi-based screens identified that GSK3β inhibits the protective NRF2-mediated oxidative stress response in the presence of reactive oxygen species elicited by high neuronal activity, which was not previously found in 2D monoculture neuron screens. We also apply the platform to investigate the role of APOE-ε4, a risk variant for Alzheimer's Disease, in its effect on neuronal survival. We find that APOE-ε4-expressing astrocytes may promote more neuronal activity as compared to APOE-ε3-expressing astrocytes. This platform expands the toolbox for the unbiased identification of mechanisms of cell-cell interactions in brain health and disease.

Systemic Inflammatory Proteomic Biomarkers in Atopic Dermatitis: Exploring Potential Indicators for Disease Severity

BACKGROUND

Atopic dermatitis (AD) is a chronic inflammatory cutaneous disorder, that emerges from intricate interplays among genetic predisposition, immune dysregulation, environmental factors, and compromised skin barrier. Understanding the inflammatory pathway in AD is important due to its fundamental role in the pathogenesis of AD. This study aimed to explore the diverse spectrum of proteins linked to the inflammation of AD and the relationship between systemic biomarkers and clinical severity in AD.

METHODS

We examined the blood samples from 48 patients with AD and 48 healthy controls (HCs) using the Proximity Extension Assay (Olink). Differentially expressed proteins (DEPs) were identified and Pearson correlation analysis was conducted to determine systemic proteomic biomarkers associated with severity of AD.

RESULTS

A total of 29 DEPs were significantly up-regulated and 2 DEPs were significantly down-regulated in AD compared with the HC. The MCP-4, IL-18, MCP-3, TNFRSF9, and IL-17C were the top 5 highest DEPs associated with the severity of AD.

CONCLUSION

Our study sheds light on the intricate network of inflammatory proteins in AD and their potential implications for disease severity. Our results indicate that these systemic inflammatory proteins could be valuable for assessing AD severity and enhancing our understanding of the disease's complexity and its potential management strategies.

HMGA1 orchestrates chromatin compartmentalization and sequesters genes into 3D networks coordinating senescence heterogeneity

HMGA1 is an abundant non-histone chromatin protein that has been implicated in embryonic development, cancer, and cellular senescence, but its specific role remains elusive. Here, we combine functional genomics approaches with graph theory to investigate how HMGA1 genomic deposition controls high-order chromatin networks in an oncogene-induced senescence model. While the direct role of HMGA1 in gene activation has been described previously, we find little evidence to support this. Instead, we show that the heterogeneous linear distribution of HMGA1 drives a specific 3D chromatin organization. HMGA1-dense loci form highly interactive networks, similar to, but independent of, constitutive heterochromatic loci. This, coupled with the exclusion of HMGA1-poor chromatin regions, leads to coordinated gene regulation through the repositioning of genes. In the absence of HMGA1, the whole process is largely reversed, but many regulatory interactions also emerge, amplifying the inflammatory senescence-associated secretory phenotype. Such HMGA1-mediated fine-tuning of gene expression contributes to the heterogeneous nature of senescence at the single-cell level. A similar 'buffer' effect of HMGA1 on inflammatory signalling is also detected in lung cancer cells. Our study reveals a mechanism through which HMGA1 modulates chromatin compartmentalization and gene regulation in senescence and beyond.

microRNA-1 Regulates Metabolic Flexibility in Skeletal Muscle via Pyruvate Metabolism

MicroRNA-1 (miR-1) is the most abundant miRNA in adult skeletal muscle. To determine the function of miR-1 in adult skeletal muscle, we generated an inducible, skeletal muscle-specific miR-1 knockout (KO) mouse. Integration of RNA-sequencing (RNA-seq) data from miR-1 KO muscle with Argonaute 2 enhanced crosslinking and immunoprecipitation sequencing (AGO2 eCLIP-seq) from human skeletal muscle identified miR-1 target genes involved with glycolysis and pyruvate metabolism. The loss of miR-1 in skeletal muscle induced cancer-like metabolic reprogramming, as shown by higher pyruvate kinase muscle isozyme M2 (PKM2) protein levels, which promoted glycolysis. Comprehensive bioenergetic and metabolic phenotyping combined with skeletal muscle proteomics and metabolomics further demonstrated that miR-1 KO induced metabolic inflexibility as a result of pyruvate oxidation resistance. While the genetic loss of miR-1 reduced endurance exercise performance in mice and in C. elegans, the physiological down-regulation of miR-1 expression in response to a hypertrophic stimulus in both humans and mice causes a similar metabolic reprogramming that supports muscle cell growth. Taken together, these data identify a novel post-translational mechanism of adult skeletal muscle metabolism regulation mediated by miR-1.

Circulating pancreatic enzyme levels are a causal biomarker of type 1 diabetes

Novel biomarkers of type 1 diabetes (T1D) are needed for earlier detection of disease and identifying therapeutic targets. We identified biomarkers of T1D by combining plasma cis and trans protein QTLs (pQTLs) for 2,922 proteins in the UK Biobank with a T1D genome-wide association study (GWAS) in 157k samples. T1D risk variants at over 20% of known loci colocalized with cis or trans pQTLs, and distinct sets of T1D loci colocalized with immune, pancreatic secretion, or gut-related proteins. We identified 23 proteins with evidence for a causal role in using pQTLs as genetic instruments in Mendelian Randomization which included multiple sensitivity analyses. Proteins increasing T1D risk were involved in immune processes (e.g. HLA-DRA) and, more surprisingly, T1D protective proteins were enriched in pancreatic secretions (e.g. CPA1), cholesterol metabolism (e.g. APOA1), and gut homeostasis. Genetic variants associated with plasma levels of T1D-protective pancreatic enzymes such as CPA1 were enriched in cis-regulatory elements in pancreatic exocrine and gut enteroendocrine cells, and the protective effects of CPA1 and other enzymes on T1D were consistent when using instruments specific to acinar cells. Finally, pancreatic enzymes had decreased acinar expression in T1D, including CPA1 which was altered prior to onset. Together, these results reveal causal biomarkers and highlight processes in the exocrine pancreas, immune system, and gut that modulate T1D risk.

Cytoplasmic binding partners of the Integrator endonuclease INTS11 and its paralog CPSF73 are required for their nuclear function

INTS11 and CPSF73 are metal-dependent endonucleases for Integrator and pre-mRNA 3'-end processing, respectively. Here, we show that the INTS11 binding partner BRAT1/CG7044, a factor important for neuronal fitness, stabilizes INTS11 in the cytoplasm and is required for Integrator function in the nucleus. Loss of BRAT1 in neural organoids leads to transcriptomic disruption and precocious expression of neurogenesis-driving transcription factors. The structures of the human INTS9-INTS11-BRAT1 and Drosophila dIntS11-CG7044 complexes reveal that the conserved C terminus of BRAT1/CG7044 is captured in the active site of INTS11, with a cysteine residue directly coordinating the metal ions. Inspired by these observations, we find that UBE3D is a binding partner for CPSF73, and UBE3D likely also uses a conserved cysteine residue to directly coordinate the active site metal ions. Our studies have revealed binding partners for INTS11 and CPSF73 that behave like cytoplasmic chaperones with a conserved impact on the nuclear functions of these enzymes.

mTOR activation induces endolysosomal remodeling and nonclassical secretion of IL-32 via exosomes in inflammatory reactive astrocytes

Astrocytes respond and contribute to neuroinflammation by adopting inflammatory reactive states. Although recent efforts have characterized the gene expression signatures associated with these reactive states, the cell biology underlying inflammatory reactive astrocyte phenotypes remains under-explored. Here, we used CRISPR-based screening in human iPSC-derived astrocytes to identify mTOR activation a driver of cytokine-induced endolysosomal system remodeling, manifesting as alkalinization of endolysosomal compartments, decreased autophagic flux, and increased exocytosis of certain endolysosomal cargos. Through endolysosomal proteomics, we identified and focused on one such cargo-IL-32, a disease-associated pro-inflammatory cytokine not present in rodents, whose secretion mechanism is not well understood. We found that IL-32 was partially secreted in extracellular vesicles likely to be exosomes. Furthermore, we found that IL-32 was involved in the polarization of inflammatory reactive astrocyte states and was upregulated in astrocytes in multiple sclerosis lesions. We believe that our results advance our understanding of cell biological pathways underlying inflammatory reactive astrocyte phenotypes and identify potential therapeutic targets.

Loss-of-function of the Zinc Finger Homeobox 4 (ZFHX4) gene underlies a neurodevelopmental disorder

8q21.11 microdeletions encompassing the gene encoding transcription factor ZFHX4, have previously been associated by us with a syndromic form of intellectual disability, hypotonia, decreased balance and hearing loss. Here, we report on 57 individuals, 52 probands and 5 affected family members, with protein truncating variants (n=36), (micro)deletions (n=20) or an inversion (n=1) affecting ZFHX4 with variable developmental delay and intellectual disability, distinctive facial characteristics, morphological abnormalities of the central nervous system, behavioral alterations, short stature, hypotonia, and occasionally cleft palate and anterior segment dysgenesis. The phenotypes associated with 8q21.11 microdeletions and ZFHX4 intragenic loss-of-function variants largely overlap, identifying ZFHX4 as the main driver for the microdeletion syndrome, although leukocyte-derived DNA shows a mild common methylation profile for (micro)deletions only. We identify ZFHX4 as a transcription factor that is increasingly expressed during human brain development and neuronal differentiation. Furthermore, ZFHX4 interacting factors identified via IP-MS in neural progenitor cells, suggest an important role for ZFHX4 in cellular and developmental pathways, especially during histone modifications, cytosolic transport and development. Additionally, using CUT&RUN, we observed that ZFHX4 binds with the promoter regions of genes with crucial roles in embryonic, neuron and axon development. Since loss-of-function variants in ZFHX4 are found with consistent dysmorphic facial features, we investigated whether the disruption of zfhx4 causes craniofacial abnormalities in zebrafish. First-generation (F0) zfhx4 crispant zebrafish, (mosaic) mutant for zfhx4 loss-of-function variants, have significantly shorter Meckel's cartilages and smaller ethmoid plates compared to control zebrafish. Furthermore, behavioral assays show a decreased movement frequency in the zfhx4 crispant zebrafish in comparison with control zebrafish larvae. Although further research is needed, our in vivo work suggests a role for zfhx4 in facial skeleton patterning, palatal development and behavior.

The effects of age and dysfunction on meibomian gland population dynamics

PURPOSE

While meibomian gland dysfunction (MGD) is widely recognized as a major cause of evaporative dry eye disease, little is known about normal gland differentiation and lipid synthesis or the mechanism underlying gland atrophy and abnormal lipid secretion. The purpose of this study was to use single-cell and spatial transcriptomics to probe changes in cell composition, differentiation, and gene expression associated with two murine models of MGD: age-related gland atrophy in wild-type mice and altered meibum quality in acyl-CoA wax alcohol acyltransferase 2 (Awat2) knockout (KO) mice.

METHODS

Young (6 month) and old (22 month) wild type, C57Bl/6 mice and young (3 month) and old (13 month) Awat2 KO mice were used in these studies. For single-cell analysis, the tarsal plate was dissected from the upper and lower eyelids, and single cells isolated and submitted to the UCI Genomic Core, while for the spatial analysis frozen tissue sections were shipped to Resolve Biosciences on dry ice and sections probed in duplicate using a meibomian gland specific, 100 gene Molecular Chartography panel.

RESULTS

Analysis of gene expression patterns identified the stratified expression of lipogenic genes during meibocyte differentiation, which may control the progressive synthesis of meibum lipids; an age-related decrease in meibocytes; and increased immune cell infiltration. Additionally, we detected unique immune cell populations in the Awat2 KO mouse suggesting activation of psoriasis-like, inflammatory pathways perhaps caused by ductal dilation and hyperplasia.

CONCLUSION

Together these findings support novel mechanism controlling gland function and dysfunction.

Canine hyper-sociability structural variants associated with altered three-dimensional chromatin state

Strong selection on complex traits can lead to skewed trait means and reduced trait variability in populations. An example of this phenomenon can be evidenced in allele frequency changes and skewed trait distributions driven by persistent human-directed selective pressures in domesticated species. Dog domestication is linked to several genomic variants; however, the functional impacts of these variants may not always be straightforward when found in non-coding regions of the genome. Four polymorphic transposable elements (TE) found within non-coding sites along a 5 Mb region on canine CFA6 have evolved due to directional selection associated with heightened human-directed hyper-sociability in domesticated dogs. We found that the polymorphic TE in intron 17 of the canine GTF2I gene, which was previously reported to be negatively correlated with canid human-directed hyper-sociability, is associated with altered chromatin looping and hence distinct cis-regulatory landscapes. We reported supporting evidence of an E2F1-DNA binding peak concordant with the altered loop and higher expression of GTF2I exon 18, indicative of alternative splicing. Globally, we discovered differences in pathways regulating the extra-cellular matrix with respect to TE copy number. Overall, we reported evidence suggesting an intriguing molecular convergence between the emergence of hypersocial behaviors in dogs and the same genes that, when hemizygous, produce human Williams Beuren Syndrome characterized by cranio-facial defects and heightened social behaviors. Our results additionally emphasize the often-overlooked potential role of chromatin architecture in social evolution.

Identification of Schizophrenia Susceptibility Loci in the Urban Taiwanese Population

Background and Objectives: Genomic studies have identified several SNP loci associated with schizophrenia in East Asian populations. Environmental factors, particularly urbanization, play a significant role in schizophrenia development. This study aimed to identify schizophrenia susceptibility loci and characterize their biological functions and molecular pathways in Taiwanese urban Han individuals. Materials and Methods: Participants with schizophrenia were recruited from the Taiwan Precision Medicine Initiative at Tri-Service General Hospital. Genotype-phenotype association analysis was performed, with significant variants annotated and analyzed for functional relevance. Results: A total of 137 schizophrenia patients and 26,129 controls were enrolled. Ten significant variants (p < 1 × 10-5) and 15 expressed genes were identified, including rs1010840 (SOWAHC and RGPD6), rs11083963 (TRPM4), rs11619878 (LINC00355 and LINC01052), rs117010638 (AGBL1 and MIR548AP), rs1170702 (LINC01680 and LINC01720), rs12028521 (KAZN and PRDM2), rs12859097 (DMD), rs1556812 (ATP11A), rs78144262 (LINC00977), and rs9997349 (ENPEP). These variants and associated genes are involved in immune response, blood pressure regulation, muscle function, and the cytoskeleton. Conclusions: Identified variants and associated genes suggest a potential genetic predisposition to schizophrenia in the Taiwanese urban Han population, highlighting the importance of potential comorbidities, considering population-specific genetic and environmental interactions.

Glucosylceramide synthase modulation ameliorates murine renal pathologies and promotes macrophage effector function in vitro

While significant advances have been made in understanding renal pathophysiology, less is known about the role of glycosphingolipid (GSL) metabolism in driving organ dysfunction. Here, we used a small molecule inhibitor of glucosylceramide synthase to modulate GSL levels in three mouse models of distinct renal pathologies: Alport syndrome (Col4a3 KO), polycystic kidney disease (Nek8jck), and steroid-resistant nephrotic syndrome (Nphs2 cKO). At the tissue level, we identified a core immune-enriched transcriptional signature that was shared across models and enriched in human polycystic kidney disease. Single nuclei analysis identified robust transcriptional changes across multiple kidney cell types, including epithelial and immune lineages. To further explore the role of GSL modulation in macrophage biology, we performed in vitro studies with homeostatic and inflammatory bone marrow-derived macrophages. Cumulatively, this study provides a comprehensive overview of renal dysfunction and the effect of GSL modulation on kidney-derived cells in the setting of renal dysfunction.

Investigating mechanical and inflammatory pathological mechanisms in osteoarthritis using MSC-derived osteocyte-like cells in 3D

Introduction

Changes to bone physiology play a central role in the development of osteoarthritis with the mechanosensing osteocyte releasing factors that drive disease progression. This study developed a humanised in vitro model to detect osteocyte responses to either interleukin-6, a driver of degeneration and bone remodelling in animal and human joint injury, or mechanical loading, to mimic osteoarthritis stimuli in joints.

Methods

Human MSC cells (Y201) were differentiated in 3-dimensional type I collagen gels in osteogenic media and osteocyte phenotype assessed by RTqPCR and immunostaining. Gels were subjected to a single pathophysiological load or stimulated with interleukin-6 with unloaded or unstimulated cells as controls. RNA was extracted 1-hour post-load and assessed by RNAseq. Markers of pain, bone remodelling, and inflammation were quantified by RT-qPCR and ELISA.

Results

Y201 cells embedded within 3D collagen gels assumed dendritic morphology and expressed mature osteocytes markers. Mechanical loading of the osteocyte model regulated 7564 genes (Padj p<0.05, 3026 down, 4538 up). 93% of the osteocyte transcriptome signature was expressed in the model with 38% of these genes mechanically regulated. Mechanically loaded osteocytes regulated 26% of gene ontology pathways linked to OA pain, 40% reflecting bone remodelling and 27% representing inflammation. Load regulated genes associated with osteopetrosis, osteoporosis and osteoarthritis. 42% of effector genes in a genome-wide association study meta-analysis were mechanically regulated by osteocytes with 10 genes representing potential druggable targets. Interleukin-6 stimulation of osteocytes at concentrations reported in human synovial fluids from patients with OA or following knee injury, regulated similar readouts to mechanical loading including markers of pain, bone remodelling, and inflammation.

Discussion

We have developed a reproducible model of human osteocyte like cells that express >90% of the genes in the osteocyte transcriptome signature. Mechanical loading and inflammatory stimulation regulated genes and proteins implicated in osteoarthritis symptoms of pain as well as inflammation and degeneration underlying disease progression. Nearly half of the genes classified as 'effectors' in GWAS were mechanically regulated in this model. This model will be useful in identifying new mechanisms underlying bone and joint pathologies and testing drugs targeting those mechanisms.

Histopathological-Based Analysis of Human Kidney Spatial Transcriptomics Data: Toward Precision Pathology

The application of spatial transcriptomics (ST) technologies is booming and has already yielded important insights across many different tissues and disease models. In nephrology, ST technologies have helped to decipher the cellular and molecular mechanisms at work in kidney diseases and have allowed the recent creation of spatially anchored human kidney atlases in healthy and diseased kidney tissues. During ST data analysis, the obtained computationally annotated clusters are often superimposed on a histologic image without their initial identification being based on the morphologic and spatial analyses of the tissues and lesions. In this study, we conduct a histopathologic-based analysis of ST data on a human kidney sample corresponding as closely as possible to the reality of the interpretation of a kidney biopsy sample in a health care or research context. This study shows the feasibility of a morphology-based approach to interpreting ST data, helping to improve our understanding of the lesion phenomena at work in chronic kidney disease at both the cellular and the molecular level. Finally, we show that our newly identified pathology-based clusters can be accurately projected onto other slides from nephrectomy or needle biopsy samples. They thus serve as a reference for analyzing other kidney tissues, paving the way for the future of molecular microscopy and precision pathology.

Identification of Disease-relevant, Sex-based Proteomic Differences in iPSC-derived Vascular Smooth Muscle

The prevalence of cardiovascular disease varies with sex, and the impact of intrinsic sex-based differences on vasculature is not well understood. Animal models can provide important insight into some aspects of human biology, however not all discoveries in animal systems translate well to humans. To explore the impact of chromosomal sex on proteomic phenotypes, we used iPSC-derived vascular smooth muscle cells from healthy donors of both sexes to identify sex-based proteomic differences and their possible effects on cardiovascular pathophysiology. Our analysis confirmed that differentiated cells have a proteomic profile more similar to healthy primary aortic smooth muscle than iPSCs. We also identified sex-based differences in iPSC- derived vascular smooth muscle in pathways related to ATP binding, glycogen metabolic process, and cadherin binding as well as multiple proteins relevant to cardiovascular pathophysiology and disease. Additionally, we explored the role of autosomal and sex chromosomes in protein regulation, identifying that proteins on autosomal chromosomes also show sex-based regulation that may affect the protein expression of proteins from autosomal chromosomes. This work supports the biological relevance of iPSC-derived vascular smooth muscle cells as a model for disease, and further exploration of the pathways identified here can lead to the discovery of sex-specific pharmacological targets for cardiovascular disease.

Significance

In this work, we have differentiated 4 male and 4 female iPSC lines into vascular smooth muscle cells, giving us the ability to identify statistically-significant sex-specific proteomic markers that are relevant to cardiovascular disease risk (such as PCK2, MTOR, IGFBP2, PTGR2, and SULTE1).