Correlation AnalyzeR: functional predictions from gene co-expression correlations

PGE2 induced miR365/IL-6/STAT3 signaling mediates dendritic cell dysfunction in cancer

AIM

To understand the mechanism of prostaglandin E2 (PGE2)-mediated immunosuppression in dendritic cells (DCs).

MAIN METHODS

In vivo experiments were conducted on 4T1 tumor bearing mice (TBM). In vitro experiments were performed in bone marrow-derived DCs (BMDCs), or spleen cells. Cytokines were monitored by ELISA/ELIspot. Gene expression was monitored by RT-PCR/flow cytometry.

KEY FINDINGS

In silico, in vitro, and in vivo experiments in 4T1 TBM revealed that PGE2 induced IL-6/pSTAT3 signaling through EP4 receptors in DCs, resulting in their dysfunction. These effects were reversed by EP4 antibody neutralization, EP4 antagonist, and STAT3 inhibitory peptides. PGE2 induced IL-6 was regulated by miR-365, as its mimic inhibited PGE2 induced IL-6 and the inhibitor increased lL-6 levels in DC. Bio-informatic analysis in human mammary cancers also revealed a strong compared co-relation between PGE2 and IL-6 (Correlation AnalyzeR) (R = 0.94). Mice bearing PTGS-2 KD 4T1 tumors had decreased tumor burden, PGE2, EP4, IL-6, and pSTAT3 signaling, along with improved DCs and T cell functions. Treatment of mice with a cyclooxygenase-2 (COX-2) inhibitor or EP4 antagonist decreased tumor burden, and this effect of EP4 antagonist was abrogated upon in vivo depletion of CD11c cells, indicating the crucial role of PGE2 signaling in DCs in tumor progression.

SIGNIFICANCE

In summary, our data highlights the importance of dendritic cells in mediating PGE2-mediated immunosuppression and the use of EP4 or STAT3 inhibitors or miR365 mimics can restore immunogenicity in cancer.

Pro-Oxidant Auranofin and Glutathione-Depleting Combination Unveils Synergistic Lethality in Glioblastoma Cells with Aberrant Epidermal Growth Factor Receptor Expression

Glioblastoma (GBM) is the most prevalent and advanced malignant primary brain tumor in adults. GBM frequently harbors epidermal growth factor receptor (EGFR) wild-type (EGFRwt) gene amplification and/or EGFRvIII activating mutation. EGFR-driven GBM relies on the thioredoxin (Trx) and/or glutathione (GSH) antioxidant systems to withstand the excessive production of reactive oxygen species (ROS). The impact of EGFRwt or EGFRvIII overexpression on the response to a Trx/GSH co-targeting strategy is unknown. In this study, we investigated Trx/GSH co-targeting in the context of EGFR overexpression in GBM. Auranofin is a thioredoxin reductase (TrxR) inhibitor, FDA-approved for rheumatoid arthritis. L-buthionine-sulfoximine (L-BSO) inhibits GSH synthesis by targeting the glutamate-cysteine ligase catalytic (GCLC) enzyme subunit. We analyzed the mechanisms of cytotoxicity of auranofin and the interaction between auranofin and L-BSO in U87MG, U87/EGFRwt, and U87/EGFRvIII GBM isogenic GBM cell lines. ROS-dependent effects were assessed using the antioxidant N-acetylsteine. We show that auranofin decreased TrxR1 activity and increased ROS. Auranofin decreased cell vitality and colony formation and increased protein polyubiquitination through ROS-dependent mechanisms, suggesting the role of ROS in auranofin-induced cytotoxicity in the three cell lines. ROS-dependent PARP-1 cleavage was associated with EGFRvIII downregulation in U87/EGFRvIII cells. Remarkably, the auranofin and L-BSO combination induced the significant depletion of intracellular GSH and synergistic cytotoxicity regardless of EGFR overexpression. Nevertheless, molecular mechanisms associated with cytotoxicity were modulated to a different extent among the three cell lines. U87/EGFRvIII exhibited the most prominent ROS increase, P-AKT(Ser-473), and AKT decrease along with drastic EGFRvIII downregulation. U87/EGFRwt and U87/EGFRvIII displayed lower basal intracellular GSH levels and synergistic ROS-dependent DNA damage compared to U87MG cells. Our study provides evidence for ROS-dependent synergistic cytotoxicity of auranofin and L-BSO combination in GBM in vitro. Unraveling the sensitivity of EGFR-overexpressing cells to auranofin alone, and synergistic auranofin and L-BSO combination, supports the rationale to repurpose this promising pro-oxidant treatment strategy in GBM.

Neil 1 deficiency facilitates chemoresistance through upregulation of RAD18 expression in ovarian cancer stem cells

Over the past decades, cancer stem cells (CSCs) have emerged as a critical subset of tumor cells associated with tumor recurrence and resistance to chemotherapy. Understanding the mechanisms underlying CSC-mediated chemoresistance is imperative for improving cancer therapy outcomes. This study delves into the regulatory role of NEIL1, a DNA glycosylase, in chemoresistance in ovarian CSCs. We first observed a decreased expression of NEIL1 in ovarian CSCs, suggesting its potential involvement in CSC regulation. Using pan-cancer analysis, we confirmed the diminished NEIL1 expression in ovarian tumors compared to normal tissues. Furthermore, NEIL1 downregulation correlated with an increase in stemness markers and enrichment of CSCs, highlighting its role in modulating CSC phenotype. Further mechanistic investigation revealed an inverse correlation between NEIL1 and RAD18 expression in ovarian CSCs. NEIL1 depletion led to heightened RAD18 expression, promoting chemoresistance possibly via enhancing Translesion DNA Synthesis (TLS)-mediated DNA lesion bypass. Moreover, dowregulation of NEIL1 results in reduced DNA damage accumulation and suppressed apoptosis in ovarian cancer. Overall, our findings unveil a novel mechanism involving NEIL1 and RAD18 in regulating chemoresistance in ovarian CSCs. Targeting this NEIL1-RAD18 axis may offer promising therapeutic strategies for combating chemoresistance and improving ovarian cancer treatment outcomes.

Membrane to cortex attachment determines different mechanical phenotypes in LGR5+ and LGR5- colorectal cancer cells

Colorectal cancer (CRC) tumors are composed of heterogeneous and plastic cell populations, including a pool of cancer stem cells that express LGR5. Whether these distinct cell populations display different mechanical properties, and how these properties might contribute to metastasis is poorly understood. Using CRC patient derived organoids (PDOs), we find that compared to LGR5- cells, LGR5+ cancer stem cells are stiffer, adhere better to the extracellular matrix (ECM), move slower both as single cells and clusters, display higher nuclear YAP, show a higher survival rate in response to mechanical confinement, and form larger transendothelial gaps. These differences are largely explained by the downregulation of the membrane to cortex attachment proteins Ezrin/Radixin/Moesin (ERMs) in the LGR5+ cells. By analyzing single cell RNA-sequencing (scRNA-seq) expression patterns from a patient cohort, we show that this downregulation is a robust signature of colorectal tumors. Our results show that LGR5- cells display a mechanically dynamic phenotype suitable for dissemination from the primary tumor whereas LGR5+ cells display a mechanically stable and resilient phenotype suitable for extravasation and metastatic growth.

Investigating the potential of Juglans regia phytoconstituents for the treatment of cervical cancer utilizing network biology and molecular docking approach

The fourth most frequent type of cancer in women and the leading cause of mortality for females worldwide is cervical cancer. Traditionally, medicinal plants have been utilized to treat various illnesses and ailments. The molecular docking method is used in the current study to look into the phytoconstituents of Juglans regia's possible anticancer effects on cervical cancer target proteins. This work uses the microarray dataset analysis of GSE63678 from the NCBI Gene Expression Omnibus database to find differentially expressed genes. Furthermore, protein-protein interactions of differentially expressed genes were constructed using network biology techniques. The top five hub genes (IGF1, FGF2, ESR1, MYL9, and MYH11) are then determined by computing topological parameters with Cytohubba. In addition, molecular docking research was performed on Juglans regia phytocompounds that were extracted from the IMPPAT database versus hub genes that had been identified. Utilizing molecular dynamics, simulation confirmed that prioritized docked complexes with low binding energies were stable.

Identification of myoferlin as a mitochondria-associated membranes component required for calcium signaling in PDAC cell lines

BACKGROUND

Pancreatic ductal adenocarcinoma is an aggressive cancer type with one of the lowest survival rates due to late diagnosis and the absence of effective treatments. A better understanding of PDAC biology will help researchers to discover the Achilles' heel of cancer cells. In that regard, our research team investigated the function of an emerging oncoprotein known as myoferlin. Myoferlin is overexpressed in PDAC and its silencing/targeting has been shown to affect cancer cell proliferation, migration, mitochondrial dynamics and metabolism. Nevertheless, our comprehension of myoferlin functions in cells remains limited. In this study, we aimed to understand the molecular mechanism linking myoferlin silencing to mitochondrial dynamics.

METHODS

Experiments were performed on two pancreas cancer cell lines, Panc-1 and MiaPaCa-2. Myoferlin localization on mitochondria was evaluated by immunofluorescence, proximity ligation assay, and cell fractionation. The presence of myoferlin in mitochondria-associated membranes was assessed by cell fractionation and its function in mitochondrial calcium transfer was evaluated using calcium flow experiments, proximity ligation assays, co-immunoprecipitation, and timelapse fluorescence microscopy in living cells.

RESULTS

Myoferlin localization on mitochondria was investigated. Our results suggest that myoferlin is unlikely to be located on mitochondria. Instead, we identified myoferlin as a new component of mitochondria-associated membranes. Its silencing significantly reduces the mitochondrial calcium level upon stimulation, probably through myoferlin interaction with the inositol 1,4,5-triphosphate receptors 3.

CONCLUSIONS

For the first time, myoferlin was specifically demonstrated to be located in mitochondria-associated membranes where it participates to calcium flow. We hypothesized that this function explains our previous results on mitochondrial dynamics. This study improves our comprehension of myoferlin localization and function in cancer biology.

TEM1/endosialin/CD248 promotes pathologic scarring and TGF-β activity through its receptor stability in dermal fibroblasts

BACKGROUND

Pathologic scars, including keloids and hypertrophic scars, represent a common form of exaggerated cutaneous scarring that is difficult to prevent or treat effectively. Additionally, the pathobiology of pathologic scars remains poorly understood. We aim at investigating the impact of TEM1 (also known as endosialin or CD248), which is a glycosylated type I transmembrane protein, on development of pathologic scars.

METHODS

To investigate the expression of TEM1, we utilized immunofluorescence staining, Western blotting, and single-cell RNA-sequencing (scRNA-seq) techniques. We conducted in vitro cell culture experiments and an in vivo stretch-induced scar mouse model to study the involvement of TEM1 in TGF-β-mediated responses in pathologic scars.

RESULTS

The levels of the protein TEM1 are elevated in both hypertrophic scars and keloids in comparison to normal skin. A re-analysis of scRNA-seq datasets reveals that a major profibrotic subpopulation of keloid and hypertrophic scar fibroblasts greatly expresses TEM1, with expression increasing during fibroblast activation. TEM1 promotes activation, proliferation, and ECM production in human dermal fibroblasts by enhancing TGF-β1 signaling through binding with and stabilizing TGF-β receptors. Global deletion of Tem1 markedly reduces the amount of ECM synthesis and inflammation in a scar in a mouse model of stretch-induced pathologic scarring. The intralesional administration of ontuxizumab, a humanized IgG monoclonal antibody targeting TEM1, significantly decreased both the size and collagen density of keloids.

CONCLUSIONS

Our data indicate that TEM1 plays a role in pathologic scarring, with its synergistic effect on the TGF-β signaling contributing to dermal fibroblast activation. Targeting TEM1 may represent a novel therapeutic approach in reducing the morbidity of pathologic scars.

Metabolic switch regulates lineage plasticity and induces synthetic lethality in triple-negative breast cancer

Metabolic reprogramming is key for cancer development, yet the mechanism that sustains triple-negative breast cancer (TNBC) cell growth despite deficient pyruvate kinase M2 (PKM2) and tumor glycolysis remains to be determined. Here, we find that deficiency in tumor glycolysis activates a metabolic switch from glycolysis to fatty acid β-oxidation (FAO) to fuel TNBC growth. We show that, in TNBC cells, PKM2 directly interacts with histone methyltransferase EZH2 to coordinately mediate epigenetic silencing of a carnitine transporter, SLC16A9. Inhibition of PKM2 leads to impaired EZH2 recruitment to SLC16A9, and in turn de-represses SLC16A9 expression to increase intracellular carnitine influx, programming TNBC cells to an FAO-dependent and luminal-like cell state. Together, these findings reveal a new metabolic switch that drives TNBC from a metabolically heterogeneous-lineage plastic cell state to an FAO-dependent-lineage committed cell state, where dual targeting of EZH2 and FAO induces potent synthetic lethality in TNBC.

Integrative transcriptome- and DNA methylation analysis of brain tissue from the temporal pole in suicide decedents and their controls

Suicide rates have increased steadily world-wide over the past two decades, constituting a serious public health crisis that creates a significant burden to affected families and the society as a whole. Suicidal behavior involves a multi-factorial etiology, including psychological, social and biological factors. Since the molecular neural mechanisms of suicide remain vastly uncharacterized, we examined transcriptional- and methylation profiles of postmortem brain tissue from subjects who died from suicide as well as their neurotypical healthy controls. We analyzed temporal pole tissue from 61 subjects, largely free from antidepressant and antipsychotic medication, using RNA-sequencing and DNA-methylation profiling using an array that targets over 850,000 CpG sites. Expression of NPAS4, a key regulator of inflammation and neuroprotection, was significantly downregulated in the suicide decedent group. Moreover, we identified a total of 40 differentially methylated regions in the suicide decedent group, mapping to seven genes with inflammatory function. There was a significant association between NPAS4 DNA methylation and NPAS4 expression in the control group that was absent in the suicide decedent group, confirming its dysregulation. NPAS4 expression was significantly associated with the expression of multiple inflammatory factors in the brain tissue. Overall, gene sets and pathways closely linked to inflammation were significantly upregulated, while specific pathways linked to neuronal development were suppressed in the suicide decedent group. Excitotoxicity as well as suppressed oligodendrocyte function were also implicated in the suicide decedents. In summary, we have identified central nervous system inflammatory mechanisms that may be active during suicidal behavior, along with oligodendrocyte dysfunction and altered glutamate neurotransmission. In these processes, NPAS4 might be a master regulator, warranting further studies to validate its role as a potential biomarker or therapeutic target in suicidality.

Cryptocaryone induces apoptosis in human hepatocellular carcinoma cells by inhibiting aerobic glycolysis through Akt and c-Src signaling pathways

Hepatocellular carcinoma (HCC) is the most common form of liver cancer, with the second highest mortality rate in all cancer. Energy reprogramming is one of the hallmarks of cancer, and emerging evidence showed that targeting glycolysis is a promising strategy for HCC treatment. Cryptocaryone has been shown to display promising anti-cancer activity against numerous types of cancer. Previous study also indicated that cryptocaryone induces cytotoxicity by inhibiting glucose transport in cancer cells, but the detailed mechanism still needs to be elucidated. Therefore, this study aimed to investigate the relationship between the anti-cancer effect and glycolytic metabolism of cryptocaryone in human HCC cells. In this study, we found that cryptocaryone potently induced growth inhibition by apoptotic cell death in HCC cells. Cryptocaryone also suppressed the ATP synthesis, lactate production and glycolytic capacity of HCC cells. Mechanistic investigations showed that phosphorylation of Akt and c-Src, as well as the expression of HK1 were impeded by cryptocaryone. Moreover, cryptocaryone markedly increased the expression level of transcription factor FoxO1. Importantly, clinical database analysis confirmed the negative correlation between HK1 and FoxO1. High expression levels of HK-1 were positively correlated with poorer survival in patients with HCCs. These results suggest that cryptocaryone may promote cell apoptosis by inhibiting FoxO1-mediated aerobic glycolysis through Akt and c-Src signaling cascades in human HCC cells. This is the first study to indicate that cryptocaryone exerts anti-cancer property against human HCC cells. Cryptocaryone is a potential natural product worthy of further development into a promising candidate for HCC treatment.

Glial Cell Adhesion Molecule (GlialCAM) Determines Proliferative versus Invasive Cell States in Glioblastoma

The malignant brain cancer glioblastoma (GBM) contains groups of highly invasive cells that drive tumor progression as well as recurrence after surgery and chemotherapy. The molecular mechanisms that enable these GBM cells to exit the primary mass and disperse throughout the brain remain largely unknown. Here we report using human tumor specimens and primary spheroids from male and female patients that glial cell adhesion molecule (GlialCAM), which has normal roles in brain astrocytes and is mutated in the developmental brain disorder megalencephalic leukoencephalopathy with subcortical cysts (MLC), is differentially expressed in subpopulations of GBM cells. High levels of GlialCAM promote cell-cell adhesion and a proliferative GBM cell state in the tumor core. In contrast, GBM cells with low levels of GlialCAM display diminished proliferation and enhanced invasion into the surrounding brain parenchyma. RNAi-mediated inhibition of GlialCAM expression leads to activation of proinvasive extracellular matrix adhesion and signaling pathways. Profiling GlialCAM-regulated genes combined with cross-referencing to single-cell transcriptomic datasets validates functional links among GlialCAM, Mlc1, and aquaporin-4 in the invasive cell state. Collectively, these results reveal an important adhesion and signaling axis comprised of GlialCAM and associated proteins including Mlc1 and aquaporin-4 that is critical for control of GBM cell proliferation and invasion status in the brain cancer microenvironment.SIGNIFICANCE STATEMENT Glioblastoma (GBM) contains heterogeneous populations of cells that coordinately drive proliferation and invasion. We have discovered that glial cell adhesion molecule (GlialCAM)/hepatocyte cell adhesion molecule (HepaCAM) is highly expressed in proliferative GBM cells within the tumor core. In contrast, GBM cells with low levels of GlialCAM robustly invade into surrounding brain tissue along blood vessels and white matter. Quantitative RNA sequencing identifies various GlialCAM-regulated genes with functions in cell-cell adhesion and signaling. These data reveal that GlialCAM and associated signaling partners, including Mlc1 and aquaporin-4, are key factors that determine proliferative and invasive cell states in GBM.

The Imageable Genome

Understanding human disease on a molecular level, and translating this understanding into targeted diagnostics and therapies are central tenets of molecular medicine1. Realizing this doctrine requires an efficient adaptation of molecular discoveries into the clinic. We present an approach to facilitate this process by describing the Imageable Genome, the part of the human genome whose expression can be assessed via molecular imaging. Using a deep learning-based hybrid human-AI pipeline, we bridge individual genes and their relevance in human diseases with specific molecular imaging methods. Cross-referencing the Imageable Genome with RNA-seq data from over 60,000 individuals reveals diagnostic, prognostic and predictive imageable genes for a wide variety of major human diseases. Having both the critical size and focus to be altered in its expression during the development and progression of any human disease, the Imageable Genome will generate new imaging tools that improve the understanding, diagnosis and management of human diseases.

Endoplasmic reticulum stress in pancreatic β cells induces incretin desensitization and β-cell dysfunction via ATF4-mediated PDE4D expression

Pancreatic β-cell dysfunction and eventual loss are key steps in the progression of type 2 diabetes (T2D). Endoplasmic reticulum (ER) stress responses, especially those mediated by the protein kinase RNA-like ER kinase and activating transcription factor 4 (PERK-ATF4) pathway, have been implicated in promoting these β-cell pathologies. However, the exact molecular events surrounding the role of the PERK-ATF4 pathway in β-cell dysfunction remain unknown. Here, we report our discovery that ATF4 promotes the expression of PDE4D, which disrupts β-cell function via a downregulation of cAMP signaling. We found that β-cell-specific transgenic expression of ATF4 led to early β-cell dysfunction and loss, a phenotype that resembles accelerated T2D. Expression of ATF4, rather than C/EBP homologous protein (CHOP), promoted PDE4D expression, reduced cAMP signaling, and attenuated responses to incretins and elevated glucose. Furthermore, we found that β-cells of leptin receptor-deficient diabetic (db/db) mice had elevated nuclear localization of ATF4 and PDE4D expression, accompanied by impaired β-cell function. Accordingly, pharmacological inhibition of the ATF4 pathway attenuated PDE4D expression in the islets and promoted incretin-simulated glucose tolerance and insulin secretion in db/db mice. Finally, we found that inhibiting PDE4 activity with selective pharmacological inhibitors improved β-cell function in both db/db mice and β-cell-specific ATF4 transgenic mice. In summary, our results indicate that ER stress causes β-cell failure via ATF4-mediated PDE4D production, suggesting the ATF4-PDE4D pathway could be a therapeutic target for protecting β-cell function during the progression of T2D.NEW & NOTEWORTHY Endoplasmic reticulum stress has been implied to cause multiple β-cell pathologies during the progression of type 2 diabetes (T2D). However, the precise molecular events underlying this remain unknown. Here, we discovered that elevated ATF4 activity, which was seen in T2D β cells, attenuated β-cell proliferation and impaired insulin secretion via PDE4D-mediated downregulation of cAMP signaling. Additionally, we demonstrated that pharmacological inhibition of the ATF4 pathway or PDE4D activity alleviated β-cell dysfunction, suggesting its therapeutic usefulness against T2D.

E3 ubiquitin ligase APC/CCdh1 regulates SLC35F2 protein turnover and inhibits cancer progression in HeLa cells

BACKGROUND

The solute carrier family 35 F2 (SLC35F2), belongs to membrane-bound carrier proteins that control various physiological functions and are activated in several cancers. However, the molecular mechanism regulating SLC35F2 protein turnover and its implication in cancer progression remains unexplored. Therefore, screening for E3 ligases that promote SLC35F2 protein degradation is essential during cancer progression.

METHODS

The immunoprecipitation and Duolink proximity ligation assays (PLA) were used to determine the interaction between APC/CCdh1 and SLC35F2 proteins. A CRISPR/Cas9-mediated knockdown and rescue experiment were used to validate the functional significance of APC/CCdh1 on SLC35F2 protein stabilization. The ubiquitination function of APC/CCdh1 on SLC35F2 protein was validated using in vitro ubiquitination assay and half-life analysis. The role of APC/CCdh1 regulating SLC35F2-mediated tumorigenesis was confirmed by in vitro oncogenic experiments in HeLa cells.

RESULTS

Based on the E3 ligase screen and in vitro biochemical experiments, we identified that APC/CCdh1 interacts with and reduces SLC35F2 protein level. APC/CCdh1 promotes SLC35F2 ubiquitination and decreases the half-life of SLC35F2 protein. On the other hand, the CRISPR/Cas9-mediated depletion of APC/CCdh1 increased SLC35F2 protein levels. The mRNA expression analysis revealed a negative correlation between APC/CCdh1 and SLC35F2 across a panel of cancer cell lines tested. Additionally, we demonstrated that depletion in APC/CCdh1 promotes SLC35F2-mediated cell proliferation, colony formation, migration, and invasion in HeLa cells.

CONCLUSION

Our study highlights that APC/CCdh1 is a critical regulator of SLC35F2 protein turnover and depletion of APC/CCdh1 promotes SLC35F2-mediated tumorigenesis. Thus, we envision that APC/CCdh1-SLC35F2 axis might be a therapeutic target in cancer.

CCR8 as a Therapeutic Novel Target: Omics-Integrated Comprehensive Analysis for Systematically Prioritizing Indications

Target identification is a crucial process in drug development, aiming to identify key proteins, genes, and signal pathways involved in disease progression and their relevance in potential therapeutic interventions. While C-C chemokine receptor 8 (CCR8) has been investigated as a candidate anti-cancer target, comprehensive multi-omics analyzes across various indications are limited. In this study, we conducted an extensive bioinformatics analysis integrating genomics, proteomics, and transcriptomics data to establish CCR8 as a promising anti-cancer drug target. Our approach encompassed data collection from diverse knowledge resources, gene function analysis, differential gene expression profiling, immune cell infiltration assessment, and strategic prioritization of target indications. Our findings revealed strong correlations between CCR8 and specific cancers, notably Breast Invasive Carcinoma (BRCA), Colon Adenocarcinoma (COAD), Head and Neck Squamous Cell Carcinoma (HNSC), Rectum adenocarcinoma (READ), Stomach adenocarcinoma (STAD), and Thyroid carcinoma (THCA). This research advances our understanding of CCR8 as a potential target for anti-cancer drug development, bridging the gap between molecular insights and creating opportunities for personalized treatment of solid tumors.

GM-CSF-activated STAT5A regulates macrophage functions and inflammation in atherosclerosis

Introduction

Inhibition of STAT5 was recently reported to reduce murine atherosclerosis. However, the role of STAT5 isoforms, and more in particular STAT5A in macrophages in the context of human atherosclerosis remains unknown.

Methods and results

Here, we demonstrate reciprocal expression regulation of STAT5A and STAT5B in human atherosclerotic lesions. The former was highly upregulated in ruptured over stable plaque and correlated with macrophage presence, a finding that was corroborated by the high chromosomal accessibility of STAT5A but not B gene in plaque macrophages. Phosphorylated STAT5 correlated with macrophages confirming its activation status. As macrophage STAT5 is activated by GM-CSF, we studied the effects of its silencing in GM-CSF differentiated human macrophages. STAT5A knockdown blunted the immune response, phagocytosis, cholesterol metabolism, and augmented apoptosis terms on transcriptional levels. These changes could partially be confirmed at functional level, with significant increases in apoptosis and decreases in lipid uptake and IL-6, IL-8, and TNFa cytokine secretion after STAT5A knockdown. Finally, inhibition of general and isoform A specific STAT5 significantly reduced the secretion of TNFa, IL-8 and IL-10 in ex vivo tissue slices of advanced human atherosclerotic plaques.

Discussion

In summary, we identify STAT5A as an important determinant of macrophage functions and inflammation in the context of atherosclerosis and show its promise as therapeutic target in human atherosclerotic plaque inflammation.

Laminin α4 Expression in Human Adipose Tissue Depots and Its Association with Obesity and Obesity Related Traits

Laminin α4 (LAMA4) is one of the main structural adipocyte basement membrane (BM) components that is upregulated during adipogenesis and related to obesity in mice and humans. We conducted RNA-seq-based gene expression analysis of LAMA4 in abdominal subcutaneous (SC) and visceral (VIS) adipose tissue (AT) depots across three human sub-cohorts of the Leipzig Obesity BioBank (LOBB) to explore the relationship between LAMA4 expression and obesity (N = 1479) in the context of weight loss (N = 65) and metabolic health (N = 42). We found significant associations of LAMA4 with body fat mass (p < 0.001) in VIS AT; higher expression in VIS AT compared to SC AT; and significant relation to metabolic health parameters e.g., body fat in VIS AT, waist (p = 0.009) and interleukin 6 (p = 0.002) in male VIS AT, and hemoglobin A1c (p = 0.008) in male SC AT. AT LAMA4 expression was not significantly different between subjects with or without obesity, metabolically healthy versus unhealthy, and obesity before versus after short-term weight loss. Our results support significant associations between obesity related clinical parameters and elevated LAMA4 expression in humans. Our work offers one of the first references for understanding the meaning of LAMA4 expression specifically in relation to obesity based on large-scale RNA-seq data.

A Two-Step Transcriptome Analysis of the Human Heart Reveals Broad and Disease-Responsive Expression of Ectopic Olfactory Receptors

G-protein-coupled receptors (GPCRs) are critical regulators of cardiac physiology and a key therapeutic target for the treatment of heart disease. Ectopic olfactory receptors (ORs) are GPCRs expressed in extra-nasal tissues which have recently emerged as new mediators in the metabolic control of cardiac function. The goals of this study were to profile OR gene expression in the human heart, to identify ORs dysregulated by heart failure caused by ischemic cardiomyopathy, and to provide evidence suggestive of a role for those altered ORs in the pathogenesis of heart failure. Left ventricular tissue from heart failure patients (n = 18) and non-failing heart samples (n = 4) were subjected to a two-step transcriptome analysis consisting of the quantification of 372 distinct OR transcripts on real-time PCR arrays and simultaneous determination of global cardiac gene expression by RNA sequencing. This strategy led to the identification of >160 ORs expressed in the human heart, including 38 receptors differentially regulated with heart failure. Co-expression analyses predicted the involvement of dysregulated ORs in the alteration of mitochondrial function, extracellular matrix remodeling, and inflammation. We provide this dataset as a resource for investigating roles of ORs in the human heart, with the hope that it will assist in the identification of new therapeutic targets for the treatment of heart failure.

CRISPR screens decode cancer cell pathways that trigger γδ T cell detection

γδ T cells are potent anticancer effectors with the potential to target tumours broadly, independent of patient-specific neoantigens or human leukocyte antigen background1-5. γδ T cells can sense conserved cell stress signals prevalent in transformed cells2,3, although the mechanisms behind the targeting of stressed target cells remain poorly characterized. Vγ9Vδ2 T cells-the most abundant subset of human γδ T cells4-recognize a protein complex containing butyrophilin 2A1 (BTN2A1) and BTN3A1 (refs. 6-8), a widely expressed cell surface protein that is activated by phosphoantigens abundantly produced by tumour cells. Here we combined genome-wide CRISPR screens in target cancer cells to identify pathways that regulate γδ T cell killing and BTN3A cell surface expression. The screens showed previously unappreciated multilayered regulation of BTN3A abundance on the cell surface and triggering of γδ T cells through transcription, post-translational modifications and membrane trafficking. In addition, diverse genetic perturbations and inhibitors disrupting metabolic pathways in the cancer cells, particularly ATP-producing processes, were found to alter BTN3A levels. This induction of both BTN3A and BTN2A1 during metabolic crises is dependent on AMP-activated protein kinase (AMPK). Finally, small-molecule activation of AMPK in a cell line model and in patient-derived tumour organoids led to increased expression of the BTN2A1-BTN3A complex and increased Vγ9Vδ2 T cell receptor-mediated killing. This AMPK-dependent mechanism of metabolic stress-induced ligand upregulation deepens our understanding of γδ T cell stress surveillance and suggests new avenues available to enhance γδ T cell anticancer activity.

Defective cathepsin Z affects EGFR expression and causes autosomal dominant palmoplantar keratoderma

BACKGROUND

The abnormal function of epidermal growth factor receptor (EGFR) has recently been shown to underlie various disorders of cornification.

OBJECTIVES

To delineate the genetic basis of a novel dominant form of palmoplantar keratoderma (PPK).

METHODS

Whole-exome (WES) and direct sequencing, quantitative real-time polymerase chain reaction, protein modelling, confocal immunofluorescence microscopy, immunoblotting, three-dimensional skin equivalents and an enzyme activity assay were used to delineate the genetic basis of a novel dominant form of PPK.

RESULTS

WES revealed heterozygous variants (c.274T > C and c.305C > T) in CTSZ (encoding cathepsin Z) in four individuals (belonging to three unrelated families) with focal PPK. Bioinformatics and protein modelling predicted the variants to be pathogenic. Previous studies have suggested that EGFR expression may be subject to cathepsin regulation. Immunofluorescence revealed reduced cathepsin Z expression in the upper epidermal layers and concomitant increased epidermal EGFR expression in patients harbouring CTSZ variants. Accordingly, human keratinocytes transfected with constructs expressing PPK-causing variants in CTSZ displayed reduced cathepsin Z enzymatic activity, as well as increased EGFR expression. In line with the role played by EGFR in the regulation of keratinocyte proliferation, human keratinocytes transfected with the PPK-causing variants showed significantly increased proliferation that was abolished upon exposure to erlotinib, an EGFR inhibitor. Similarly, downregulation of CTSZ resulted in increased EGFR expression and increased proliferation in human keratinocytes, suggestive of a loss-of-function effect of the pathogenic variants. Finally, three-dimensional organotypic skin equivalents grown from CTSZ-downregulated cells showed increased epidermal thickness and EGFR expression as seen in patient skin; here, too, erlotinib was found to rescue the abnormal phenotype.

CONCLUSIONS

Taken collectively, these observations attribute to cathepsin Z a hitherto unrecognized function in epidermal differentiation.

Silencing Glypican-1 enhances the antitumor effects of Pictilisib via downregulating PI3K/Akt/ERK signaling in chemo-resistant esophageal adenocarcinoma

Poorly differentiated esophageal adenocarcinoma (PDEAC) has a dismal prognosis. Glypican-1(GPC-1) is known to be upregulated in several cancer types in contrast to healthy tissues, rendering it as a biomarker. Nevertheless, the potential therapeutic targeting of GPC-1 has not been explored in PDEAC. There is accumulating evidence that GPC-1, via upregulation of PI3K/Akt/ERK signaling, plays a crucial role in the progression and chemoresistance in cancer. Pictilisib, a class I pan PI3K inhibitor, has shown promising antitumor results in clinical trials, however, has not gained widespread success due to acquired drug resistance. This study investigated the role of GPC-1 in chemo-resistant PDEAC and appraises the impact of targeted silencing of GPC-1 on the antitumor effects of Pictilisib in PDEAC cell lines. Immunohistochemistry assays in PDEAC tissue specimens demonstrated a pronounced intensity of staining with GPC-1. Upregulation of GPC-1 was found to be correlated with advanced stage and poor prognosis. In-vitro studies examined the influence of GPC-1 knockdown and Pictilisib, both as individual agents and in combination, on cytotoxicity, cell cycle distribution, apoptosis, and gene expression profiles. Silencing GPC-1 alone showed significantly reduced cell viability, migration, colony formation, epithelial-mesenchymal transition, and stemness in PDEAC cells. Significantly, knockdown of GPC-1 combined with low-dose Pictilisib led to enhancement of cytotoxicity, cell cycle arrest, and apoptosis in ESO-26 and OE-33 cells. In the xenograft mouse model, the combination of Pictilisib and GPC-1 knockdown exhibited synergy. These findings suggest that GPC-1 represents a promising target to augment chemosensitivity in esophageal adenocarcinoma.

Nicotinamide enhances osteoblast differentiation through activation of the mitochondrial antioxidant defense system

Although the normal physiological level of oxidative stress is beneficial for maintaining bone homeostasis, imbalance between reactive oxygen species (ROS) production and antioxidant defense can cause various bone diseases. The purpose of this study was to determine whether nicotinamide (NAM), an NAD⁺ precursor, can support the maintenance of bone homeostasis by regulating osteoblasts. Here, we found that NAM enhances osteoblast differentiation and mitochondrial metabolism. NAM increases the expression of antioxidant enzymes, which is due to increased FOXO3A transcriptional activity via SIRT3 activation. NAM has not only a preventive effect against weak and chronic oxidative stress but also a therapeutic effect against strong and acute exposure to H₂O₂ in osteoblast differentiation. Collectively, the results indicate that NAM increases mitochondrial biogenesis and antioxidant enzyme expression through activation of the SIRT3-FOXO3A axis, which consequently enhances osteoblast differentiation. These results suggest that NAM could be a potential preventive or therapeutic agent for bone diseases caused by ROS.

Patterns of transcription factor binding and epigenome at promoters allow interpretable predictability of multiple functions of non-coding and coding genes

Understanding the biological roles of all genes only through experimental methods is challenging. A computational approach with reliable interpretability is needed to infer the function of genes, particularly for non-coding RNAs. We have analyzed genomic features that are present across both coding and non-coding genes like transcription factor (TF) and cofactor ChIP-seq (823), histone modifications ChIP-seq (n = 621), cap analysis gene expression (CAGE) tags (n = 255), and DNase hypersensitivity profiles (n = 255) to predict ontology-based functions of genes. Our approach for gene function prediction was reliable (>90% balanced accuracy) for 486 gene-sets. PubMed abstract mining and CRISPR screens supported the inferred association of genes with biological functions, for which our method had high accuracy. Further analysis revealed that TF-binding patterns at promoters have high predictive strength for multiple functions. TF-binding patterns at the promoter add an unexplored dimension of explainable regulatory aspects of genes and their functions. Therefore, we performed a comprehensive analysis for the functional-specificity of TF-binding patterns at promoters and used them for clustering functions to reveal many latent groups of gene-sets involved in common major cellular processes. We also showed how our approach could be used to infer the functions of non-coding genes using the CRISPR screens of coding genes, which were validated using a long non-coding RNA CRISPR screen. Thus our results demonstrated the generality of our approach by using gene-sets from CRISPR screens. Overall, our approach opens an avenue for predicting the involvement of non-coding genes in various functions.

LINC01137/miR-186-5p/WWOX: a novel axis identified from WWOX-related RNA interactome in bladder cancer

Introduction: The discovery of non-coding RNA (ncRNA) dates back to the pre-genomics era, but the progress in this field is still dynamic and leverages current post-genomics solutions. WWOX is a global gene expression modulator that is scarcely investigated for its role in regulating cancer-related ncRNAs. In bladder cancer (BLCA), the link between WWOX and ncRNA remains unexplored. The description of AP-2α and AP-2γ transcription factors, known as WWOX-interacting proteins, is more commonplace regarding ncRNA but still merits investigation. Therefore, this in vitro and in silico study aimed to construct an ncRNA-containing network with WWOX/AP-2 and to investigate the most relevant observation in the context of BLCA cell lines and patients. Methods: RT-112, HT-1376, and CAL-29 cell lines were subjected to two stable lentiviral transductions. High-throughput sequencing of cellular variants (deposited in the Gene Expression Omnibus database under the GSE193659 record) enabled the investigation of WWOX/AP-2-dependent differences using various bioinformatics tools (e.g., limma-voom, FactoMineR, multiple Support Vector Machine Recursive Feature Elimination (mSVM-RFE), miRDB, Arena-Idb, ncFANs, RNAhybrid, TargetScan, Protein Annotation Through Evolutionary Relationships (PANTHER), Gene Transcription Regulation Database (GTRD), or Evaluate Cutpoints) and repositories such as The Cancer Genome Atlas (TCGA) and Cancer Cell Line Encyclopedia. The most relevant observations from cap analysis gene expression sequencing (CAGE-seq) were confirmed using real-time PCR, whereas TCGA data were validated using the GSE31684 cohort. Results: The first stage of the whole study justified focusing solely on WWOX rather than on WWOX combined with AP-2α/γ. The most relevant observation of the developed ncRNA-containing network was LINC01137, i.e., long non-coding RNAs (lncRNAs) that unraveled the core network containing UPF1, ZC3H12A, LINC01137, WWOX, and miR-186-5p, the last three being a novel lncRNA/miRNA/mRNA axis. Patients' data confirmed the LINC01137/miR-186-5p/WWOX relationship and provided a set of dependent genes (i.e., KRT18, HES1, VCP, FTH1, IFITM3, RAB34, and CLU). Together with the core network, the gene set was subjected to survival analysis for both TCGA-BLCA and GSE31684 patients, which indicated that the increased expression of WWOX or LINC01137 is favorable, similar to their combination with each other (WWOX↑ and LINC01137↑) or with MIR186 (WWOX↑/LINC01137↑ but MIR186↓). Conclusion: WWOX is implicated in the positive feedback loop with LINC01137 that sponges WWOX-targeting miR-186-5p. This novel WWOX-containing lncRNA/miRNA/mRNA axis should be further investigated to depict its relationships in a broader context, which could contribute to BLCA research and treatment.

A Study Based on Network Pharmacology Decoding the Multi-Target Mechanism of Duhuo Jisheng Decoction for the Treatment of Intervertebral Disc Degeneration

Intervertebral disc degeneration (IDD) poses a grim public health impact. Duhuo Jisheng Decoction (DJD), a traditional Chinese medicine formula, has recently received significant attention for its efficacy and safety in treating IDD. However, the pathological processes of IDD in which DJD interferes and molecular mechanism involved are poorly understood, which brings difficulties to the clinical practice of DJD for the treatment of IDD. This study systematically investigated the underlying mechanism of DJD treatment of IDD. Network pharmacology approaches were employed, integrating molecular docking and random walk with restart (RWR) algorithm, to identify key compounds and targets for DJD in the treatment of IDD. Bioinformatics approaches were used to further explore the biological insights in DJD treatment of IDD. The analysis identifies AKT1, PIK3R1, CHUK, ALB, TP53, MYC, NR3C1, IL1B, ERBB2, CAV1, CTNNB1, AR, IGF2, and ESR1 as key targets. Responses to mechanical stress, oxidative stress, cellular inflammatory responses, autophagy, and apoptosis are identified as the critical biological processes involved in DJD treatment of IDD. The regulation of DJD targets in extracellular matrix components, ion channel regulation, transcriptional regulation, synthesis and metabolic regulation of reactive oxygen products in the respiratory chain and mitochondria, fatty acid oxidation, the metabolism of Arachidonic acid, and regulation of Rho and Ras protein activation are found to be potential mechanisms in disc tissue response to mechanical stress and oxidative stress. MAPK, PI3K/AKT, and NF-κB signaling pathways are identified as vital signaling pathways for DJD to treat IDD. Quercetin and Kaempferol are assigned a central position in the treatment of IDD. This study contributes to a more comprehensive understanding of the mechanism of DJD in treating IDD. It provides a reference for applying natural products to delay the pathological process of IDD.

Tetraspanin Tspan8 restrains interferon signaling to stabilize intestinal epithelium by directing endocytosis of interferon receptor

Inflammation can impair intestinal barrier, while increased epithelial permeability can lead to inflammation. In this study, we found that the expression of Tspan8, a tetraspanin expressed specifically in epithelial cells, is downregulated in mouse model of ulcerative disease (UC) but correlated with those of cell-cell junction components, such as claudins and E-cadherin, suggesting that Tspan8 supports intestinal epithelial barrier. Tspan8 removal increases intestinal epithelial permeability and upregulates IFN-γ-Stat1 signaling. We also demonstrated that Tspan8 coalesces with lipid rafts and facilitates IFNγ-R1 localization at or near lipid rafts. As IFN-γ induces its receptor undergoing clathrin- or lipid raft-dependent endocytosis and IFN-γR endocytosis plays an important role in Jak-Stat1 signaling, our analysis on IFN-γR endocytosis revealed that Tspan8 silencing impairs lipid raft-mediated but promotes clathrin-mediated endocytosis of IFN-γR1, leading to increased Stat1 signaling. These changes in IFN-γR1 endocytosis upon Tspan8 silencing correlates with fewer lipid raft component GM1 at the cell surface and more clathrin heavy chain in the cells. Our findings indicate that Tspan8 determines the IFN-γR1 endocytosis route, to restrain Stat1 signaling, stabilize intestine epithelium, and subsequently prevent intestine from inflammation. Our finding also implies that Tspan8 is needed for proper endocytosis through lipid rafts.

CRISPR/Cas9-based genome-wide screening of the deubiquitinase subfamily identifies USP3 as a protein stabilizer of REST blocking neuronal differentiation and promotes neuroblastoma tumorigenesis

BACKGROUND

The repressor element-1 silencing transcription factor (REST), a master transcriptional repressor, is essential for maintenance, self-renewal, and differentiation in neuroblastoma. An elevated expression of REST is associated with impaired neuronal differentiation, which results in aggressive neuroblastoma formation. E3 ligases are known to regulate REST protein abundance through the 26 S proteasomal degradation pathway in neuroblastoma. However, deubiquitinating enzymes (DUBs), which counteract the function of E3 ligase-mediated REST protein degradation and their impact on neuroblastoma tumorigenesis have remained unexplored.

METHODS

We employed a CRISPR/Cas9 system to perform a genome-wide knockout of ubiquitin-specific proteases (USPs) and used western blot analysis to screen for DUBs that regulate REST protein abundance. The interaction between USP3 and REST was confirmed by immunoprecipitation and Duolink in situ proximity assays. The deubiquitinating effect of USP3 on REST protein degradation, half-life, and neuronal differentiation was validated by immunoprecipitation, in vitro deubiquitination, protein-turnover, and immunostaining assays. The correlation between USP3 and REST expression was assessed using patient neuroblastoma datasets. The USP3 gene knockout in neuroblastoma cells was performed using CRISPR/Cas9, and the clinical relevance of USP3 regulating REST-mediated neuroblastoma tumorigenesis was confirmed by in vitro and in vivo oncogenic experiments.

RESULTS

We identified a deubiquitinase USP3 that interacts with, stabilizes, and increases the half-life of REST protein by counteracting its ubiquitination in neuroblastoma. An in silico analysis showed a correlation between USP3 and REST in multiple neuroblastoma cell lines and identified USP3 as a prognostic marker for overall survival in neuroblastoma patients. Silencing of USP3 led to a decreased self-renewal capacity and promoted retinoic acid-induced differentiation in neuroblastoma. A loss of USP3 led to attenuation of REST-mediated neuroblastoma tumorigenesis in a mouse xenograft model.

CONCLUSION

The findings of this study indicate that USP3 is a critical factor that blocks neuronal differentiation, which can lead to neuroblastoma. We envision that targeting USP3 in neuroblastoma tumors might provide an effective therapeutic differentiation strategy for improved survival rates of neuroblastoma patients.

Expression of Androgen and Estrogen Receptors in the Human Lacrimal Gland

Lacrimal gland dysfunction causes dry eye disease (DED) due to decreased tear production. Aqueous-deficient DED is more prevalent in women, suggesting that sexual dimorphism of the human lacrimal gland could be a potential cause. Sex steroid hormones are a key factor in the development of sexual dimorphism. This study aimed to quantify estrogen receptor (ER) and androgen receptor (AR) expression in the human lacrimal gland and compare it between sexes. RNA was isolated from 35 human lacrimal gland tissue samples collected from 19 cornea donors. AR, ERα, and ERβ mRNA was identified in all samples, and their expression was quantified using qPCR. Immunohistochemical staining was performed on selected samples to evaluate protein expression of the receptors. ERα mRNA expression was significantly higher than the expression of AR and ERβ. No difference in sex steroid hormone (SSH) receptor mRNA expression was observed between sexes, and no correlation was observed with age. If ERα protein expression is found to be concordant with mRNA expression, it should be investigated further as a potential target for hormone therapy of DED. Further research is needed to elucidate the role of sex steroid hormone receptors in sex-related differences of lacrimal gland structure and disease.

Transcriptomic and functional analyses reveal a tumour-promoting role for the IL-36 receptor in colon cancer and crosstalk between IL-36 signalling and the IL-17/ IL-23 axis

BACKGROUND

The interleukin (IL)-36 cytokines are a sub-family of the IL-1 family which are becoming increasingly implicated in the pathogenesis of inflammatory diseases and malignancies. Initial studies of IL-36 signalling in tumorigenesis identified an immune-mediated anti-tumorigenic function for these cytokines. However, more recent studies have shown IL-36 cytokines also contribute to the pathogenesis of lung and colorectal cancer (CRC).

METHODS

The aim of this study was to investigate IL-36 expression in CRC using transcriptomic datasets and software such as several R packages, Cytoscape, GEO2R and AnalyzeR. Validation of results was completed by qRT-PCR on both cell lines and a patient cohort. Cellular proliferation was assessed by flow cytometry and resazurin reduction.

RESULTS

We demonstrate that IL-36 gene expression increases with CRC development. Decreased tumoral IL-36 receptor expression was shown to be associated with improved patient outcome. Our differential gene expression analysis revealed a novel role for the IL-36/IL-17/IL-23 axis, with these findings validated using patient-derived samples and cell lines. IL-36γ, together with either IL-17a or IL-22, was able to synergistically induce different genes involved in the IL-17/IL-23 axis in CRC cells and additively induce colon cancer cell proliferation.

CONCLUSIONS

Collectively, this data support a pro-tumorigenic role for IL-36 signalling in colon cancer, with the IL-17/IL-23 axis influential in IL-36-mediated colon tumorigenesis.

HGCA2.0: An RNA-Seq Based Webtool for Gene Coexpression Analysis in Homo sapiens

Genes with similar expression patterns in a set of diverse samples may be considered coexpressed. Human Gene Coexpression Analysis 2.0 (HGCA2.0) is a webtool which studies the global coexpression landscape of human genes. The website is based on the hierarchical clustering of 55,431 Homo sapiens genes based on a large-scale coexpression analysis of 3500 GTEx bulk RNA-Seq samples of healthy individuals, which were selected as the best representative samples of each tissue type. HGCA2.0 presents subclades of coexpressed genes to a gene of interest, and performs various built-in gene term enrichment analyses on the coexpressed genes, including gene ontologies, biological pathways, protein families, and diseases, while also being unique in revealing enriched transcription factors driving coexpression. HGCA2.0 has been successful in identifying not only genes with ubiquitous expression patterns, but also tissue-specific genes. Benchmarking showed that HGCA2.0 belongs to the top performing coexpression webtools, as shown by STRING analysis. HGCA2.0 creates working hypotheses for the discovery of gene partners or common biological processes that can be experimentally validated. It offers a simple and intuitive website design and user interface, as well as an API endpoint.

RASA2 ablation in T cells boosts antigen sensitivity and long-term function

The efficacy of adoptive T cell therapies for cancer treatment can be limited by suppressive signals from both extrinsic factors and intrinsic inhibitory checkpoints1,2. Targeted gene editing has the potential to overcome these limitations and enhance T cell therapeutic function3-10. Here we performed multiple genome-wide CRISPR knock-out screens under different immunosuppressive conditions to identify genes that can be targeted to prevent T cell dysfunction. These screens converged on RASA2, a RAS GTPase-activating protein (RasGAP) that we identify as a signalling checkpoint in human T cells, which is downregulated upon acute T cell receptor stimulation and can increase gradually with chronic antigen exposure. RASA2 ablation enhanced MAPK signalling and chimeric antigen receptor (CAR) T cell cytolytic activity in response to target antigen. Repeated tumour antigen stimulations in vitro revealed that RASA2-deficient T cells show increased activation, cytokine production and metabolic activity compared with control cells, and show a marked advantage in persistent cancer cell killing. RASA2-knockout CAR T cells had a competitive fitness advantage over control cells in the bone marrow in a mouse model of leukaemia. Ablation of RASA2 in multiple preclinical models of T cell receptor and CAR T cell therapies prolonged survival in mice xenografted with either liquid or solid tumours. Together, our findings highlight RASA2 as a promising target to enhance both persistence and effector function in T cell therapies for cancer treatment.

CRISPR/Cas9-based genome-wide screening for deubiquitinase subfamily identifies USP1 regulating MAST1-driven cisplatin-resistance in cancer cells

Background: Cisplatin is one of the frontline anticancer agents. However, development of cisplatin-resistance limits the therapeutic efficacy of cisplatin-based treatment. The expression of microtubule-associated serine/threonine kinase 1 (MAST1) is a primary factor driving cisplatin-resistance in cancers by rewiring the MEK pathway. However, the mechanisms responsible for MAST1 regulation in conferring drug resistance is unknown.

Methods: We implemented a CRISPR/Cas9-based, genome-wide, dual screening system to identify deubiquitinating enzymes (DUBs) that govern cisplatin resistance and regulate MAST1 protein level. We analyzed K48- and K63-linked polyubiquitination of MAST1 protein and mapped the interacting domain between USP1 and MAST1 by immunoprecipitation assay. The deubiquitinating effect of USP1 on MAST1 protein was validated using rescue experiments, in vitro deubiquitination assay, immunoprecipitation assays, and half-life analysis. Furthermore, USP1-knockout A549 lung cancer cells were generated to validate the deubiquitinating activity of USP1 on MAST1 abundance. The USP1-MAST1 correlation was evaluated using bioinformatics tool and in different human clinical tissues. The potential role of USP1 in regulating MAST1-mediated cisplatin resistance was confirmed using a series of in vitro and in vivo experiments. Finally, the clinical relevance of the USP1-MAST1 axis was validated by application of small-molecule inhibitors in a lung cancer xenograft model in NSG mice.

Results: The CRISPR/Cas9-based dual screening system identified USP1 as a novel deubiquitinase that interacts, stabilizes, and extends the half-life of MAST1 by preventing its K48-linked polyubiquitination. The expression analysis across human clinical tissues revealed a positive correlation between USP1 and MAST1. USP1 promotes MAST1-mediated MEK1 activation as an underlying mechanism that contributes to cisplatin-resistance in cancers. Loss of USP1 led to attenuation of MAST1-mediated cisplatin-resistance both in vitro and in vivo. The combined pharmacological inhibition of USP1 and MAST1 using small-molecule inhibitors further abrogated MAST1 level and synergistically enhanced cisplatin efficacy in a mouse xenograft model.

Conclusions: Overall, our study highlights the role of USP1 in the development of cisplatin resistance and uncovers the regulatory mechanism of MAST1-mediated cisplatin resistance in cancers. Co-treatment with USP1 and MAST1 inhibitors abrogated tumor growth and synergistically enhanced cisplatin efficacy, suggesting a novel alternative combinatorial therapeutic strategy that could further improve MAST1-based therapy in patients with cisplatin-resistant tumors.

CDC20-Mediated hnRNPU Ubiquitination Regulates Chromatin Condensation and Anti-Cancer Drug Response

Cell division cycle 20 (CDC20) functions as a critical cell cycle regulator. It plays an important role in cancer development and drug resistance. However, the molecular mechanisms by which CDC20 regulates cellular drug response remain poorly understood. Chromatin-associated CDC20 interactome in breast cancer cells was analyzed by using affinity purification coupled with mass spectrometry. hnRNPU as a CDC20 binding partner was validated by co-immunoprecipitation and immunostaining. The molecular domain, comprising amino acid residues 461-653, on hnRNPU required for its interaction with CDC20 was identified by mapping of interactions. Co-immunoprecipitation showed that CDC20-mediated hnRNPU ubiquitination promotes its interaction with the CTCF and cohesin complex. The effects of CDC20-hnRNPU on nuclear size and chromatin condensation were investigated by analyzing DAPI and H2B-mCherry staining, respectively. The role of CDC20-hnRNPU in tumor progression and drug resistance was examined by CCK-8 cell survival and clonogenic assays. Our study indicates that CDC20-mediated ubiquitination of hnRNPU modulates chromatin condensation by regulating the interaction between hnRNPU and the CTCF-cohesin complex. Dysregulation of the CDC20-hnRNPU axis contributes to tumor progression and drug resistance.

Quality-controlled R-loop meta-analysis reveals the characteristics of R-loop consensus regions

R-loops are three-stranded nucleic acid structures formed from the hybridization of RNA and DNA. While the pathological consequences of R-loops have been well-studied to date, the locations, classes, and dynamics of physiological R-loops remain poorly understood. R-loop mapping studies provide insight into R-loop dynamics, but their findings are challenging to generalize. This is due to the narrow biological scope of individual studies, the limitations of each mapping modality, and, in some cases, poor data quality. In this study, we reprocessed 810 R-loop mapping datasets from a wide array of biological conditions and mapping modalities. From this data resource, we developed an accurate R-loop data quality control method, and we reveal the extent of poor-quality data within previously published studies. We then identified a set of high-confidence R-loop mapping samples and used them to define consensus R-loop sites called 'R-loop regions' (RL regions). In the process, we identified a stark divergence between RL regions detected by S9.6 and dRNH-based mapping methods, particularly with respect to R-loop size, location, and colocalization with RNA binding factors. Taken together, this work provides a much-needed method to assess R-loop data quality and offers novel context regarding the differences between dRNH- and S9.6-based R-loop mapping approaches.

Approaches in Gene Coexpression Analysis in Eukaryotes

Simple Summary

Genes whose expression levels rise and fall similarly in a large set of samples, may be considered coexpressed. Gene coexpression analysis refers to the en masse discovery of coexpressed genes from a large variety of transcriptomic experiments. The type of biological networks that studies gene coexpression, known as Gene Coexpression Networks, consist of an undirected graph depicting genes and their coexpression relationships. Coexpressed genes are clustered in smaller subnetworks, the predominant biological roles of which can be determined through enrichment analysis. By studying well-annotated gene partners, the attribution of new roles to genes of unknown function or assumption for participation in common metabolic pathways can be achieved, through a guilt-by-association approach. In this review, we present key issues in gene coexpression analysis, as well as the most popular tools that perform it.

Abstract

Gene coexpression analysis constitutes a widely used practice for gene partner identification and gene function prediction, consisting of many intricate procedures. The analysis begins with the collection of primary transcriptomic data and their preprocessing, continues with the calculation of the similarity between genes based on their expression values in the selected sample dataset and results in the construction and visualisation of a gene coexpression network (GCN) and its evaluation using biological term enrichment analysis. As gene coexpression analysis has been studied extensively, we present most parts of the methodology in a clear manner and the reasoning behind the selection of some of the techniques. In this review, we offer a comprehensive and comprehensible account of the steps required for performing a complete gene coexpression analysis in eukaryotic organisms. We comment on the use of RNA-Seq vs. microarrays, as well as the best practices for GCN construction. Furthermore, we recount the most popular webtools and standalone applications performing gene coexpression analysis, with details on their methods, features and outputs.

Reference Module-Based Analysis of Ovarian Cancer Transcriptome Identifies Important Modules and Potential Drugs

Ovarian cancer (OVC) is often diagnosed at the advanced stage resulting in a poor overall outcome for the patient. The disease mechanisms, prognosis, and treatment require imperative elucidation. A rank-based module-centric framework was proposed to analyze the key modules related to the development, prognosis, and treatment of OVC. The ovarian cancer cell line microarray dataset GSE43765 from the Gene Expression Omnibus database was used to construct the reference modules by weighted gene correlation network analysis. Twenty-three reference modules were tested for stability and functionally annotated. Furthermore, to demonstrate the utility of reference modules, two more OVC datasets were collected, and their gene expression profiles were projected to the reference modules to generate a module-level expression. An epithelial-mesenchymal transition module was activated in OVC compared to the normal epithelium, and a pluripotency module was activated in ovarian cancer stroma compared to ovarian cancer epithelium. Seven differentially expressed modules were identified in OVC compared to the normal ovarian epithelium, with five up-regulated, and two down-regulated. One module was identified to be predictive of patient overall survival. Four modules were enriched with SNP signals. Based on differentially expressed modules and hub genes, five candidate drugs were screened. The hub genes of those modules merit further investigation. We firstly propose the reference module-based analysis of OVC. The utility of the analysis framework can be extended to transcriptome data of other kinds of diseases.