UBE2N is essential for maintenance of skin homeostasis and suppression of inflammation.

UBE2N, a Lys63-ubiquitin conjugating enzyme, plays critical roles in embryogenesis and immune system development and function. However, its roles in adult epithelial tissue homeostasis and pathogenesis are unclear. We generated conditional mouse models that deleted Ube2n in skin cells in a temporally and spatially controlled manner. We found that Ube2n-knockout (KO) in the adult skin keratinocytes induced a range of inflammatory skin defects characteristic of psoriatic and actinic keratosis. These included inflammation, epidermal and dermal thickening, parakeratosis, and increased immune cell infiltration, as well as signs of edema and blistering. Single cell transcriptomic analyses and RT-qPCR showed that Ube2n KO keratinocytes expressed elevated myeloid cell chemo-attractants such as Cxcl1 and Cxcl2 and decreased the homeostatic T lymphocyte chemo-attractant Ccl27a. Consistently, the infiltrating immune cells were predominantly myeloid-derived cells including neutrophils and M1-like macrophages that expressed high levels of inflammatory cytokines such as Il1ß and Il24. Pharmacological blockade of the IL-1 receptor associated kinases (IRAK1/4) alleviated inflammation, epidermal and dermal thickening, and immune infiltration of the Ube2n mutant skin. Together, these findings highlight a key role of keratinocyte-UBE2N in maintenance of epidermal homeostasis and skin immunity, and identify IRAK1/4 as potential therapeutic target for inflammatory skin disorders.

Spatial transcriptomics reveals segregation of tumor cell states in glioblastoma and marked immunosuppression within the perinecrotic niche.

Glioblastoma (GBM) remains an untreatable malignant tumor with poor patient outcomes, characterized by palisading necrosis and microvascular proliferation. While single-cell technology made it possible to characterize different lineage of glioma cells into neural progenitor-like (NPC-like), oligodendrocyte-progenitor-like (OPC-like), astrocyte-like (AC-like) and mesenchymal like (MES-like) states, it does not capture the spatial localization of these tumor cell states. Spatial transcriptomics empowers the study of the spatial organization of different cell types and tumor cell states and allows for the selection of regions of interest to investigate region-specific and cell-type-specific pathways. Here, we obtained paired 10x Chromium single-nuclei RNA-sequencing (snRNA-seq) and 10x Visium spatial transcriptomics data from three GBM patients to interrogate the GBM microenvironment. Integration of the snRNA-seq and spatial transcriptomics data reveals patterns of segregation of tumor cell states. For instance, OPC-like tumor and NPC-like tumor significantly segregate in two of the three samples. Our differentially expressed gene and pathway analyses uncovered significant pathways in functionally relevant niches. Specifically, perinecrotic regions were more immunosuppressive than the endogenous GBM microenvironment, and perivascular regions were more pro-inflammatory. Our gradient analysis suggests that OPC-like tumor cells tend to reside in areas closer to the tumor vasculature compared to tumor necrosis, which may reflect increased oxygen requirements for OPC-like cells. In summary, we characterized the localization of cell types and tumor cell states, the gene expression patterns, and pathways in different niches within the GBM microenvironment. Our results provide further evidence of the segregation of tumor cell states and highlight the immunosuppressive nature of the necrotic and perinecrotic niches in GBM.

An osteoarthritis pathophysiological continuum revealed by molecular biomarkers.

We aimed to identify serum biomarkers that predict knee osteoarthritis (OA) before the appearance of radiographic abnormalities in a cohort of 200 women. As few as six serum peptides, corresponding to six proteins, reached AUC 77% probability to distinguish those who developed OA from age-matched individuals who did not develop OA up to 8 years later. Prediction based on these blood biomarkers was superior to traditional prediction based on age and BMI (AUC 51%) or knee pain (AUC 57%). These results identify a prolonged molecular derangement of joint tissue before the onset of radiographic OA abnormalities consistent with an unresolved acute phase response. Among all 24 protein biomarkers predicting incident knee OA, the majority (58%) also predicted knee OA progression, revealing the existence of a pathophysiological "OA continuum" based on considerable similarity in the molecular pathophysiology of the progression to incident OA and the progression of established OA.

Gli2 Facilitates Tumor Immune Evasion and Immunotherapeutic Resistance by Coordinating Wnt Ligand and Prostaglandin Signaling.

Therapeutic resistance to immune checkpoint blockade has been commonly linked to the process of mesenchymal transformation (MT) and remains a prevalent obstacle across many cancer types. An improved mechanistic understanding for MT-mediated immune evasion promises to lead to more effective combination therapeutic regimens. Herein, we identify the Hedgehog transcription factor, Gli2, as a key node of tumor-mediated immune evasion and immunotherapy resistance during MT. Mechanistic studies reveal that Gli2 generates an immunotolerant tumor microenvironment through the upregulation of Wnt ligand production and increased prostaglandin synthesis. This pathway drives the recruitment, viability, and function of granulocytic myeloid-derived suppressor cells (PMN-MDSCs) while also impairing type I conventional dendritic cell, CD8 + T cell, and NK cell functionality. Pharmacologic EP2/EP4 prostaglandin receptor inhibition and Wnt ligand inhibition each reverses a subset of these effects, while preventing primary and adaptive resistance to anti-PD-1 immunotherapy, respectively. A transcriptional Gli2 signature correlates with resistance to anti-PD-1 immunotherapy in stage IV melanoma patients, providing a translational roadmap to direct combination immunotherapeutics in the clinic. SIGNIFICANCE: Gli2-induced EMT promotes immune evasion and immunotherapeutic resistance via coordinated prostaglandin and Wnt signaling.

Subventricular zone stem cell niche injury is associated with intestinal perforation in preterm infants and predicts future motor impairment.

Brain injury is highly associated with preterm birth. Complications of prematurity, including spontaneous or necrotizing enterocolitis (NEC)-associated intestinal perforations, are linked to lifelong neurologic impairment, yet the mechanisms are poorly understood. Early diagnosis of preterm brain injuries remains a significant challenge. Here, we identified subventricular zone echogenicity (SVE) on cranial ultrasound in preterm infants following intestinal perforations. The development of SVE was significantly associated with motor impairment at 2 years. SVE was replicated in a neonatal mouse model of intestinal perforation. Examination of the murine echogenic subventricular zone (SVZ) revealed NLRP3-inflammasome assembly in multiciliated FoxJ1+ ependymal cells and a loss of the ependymal border in this postnatal stem cell niche. These data suggest a mechanism of preterm brain injury localized to the SVZ that has not been adequately considered. Ultrasound detection of SVE may serve as an early biomarker for neurodevelopmental impairment after inflammatory disease in preterm infants.

PI3K/mTOR is a therapeutically targetable genetic dependency in diffuse intrinsic pontine glioma.

Diffuse midline glioma (DMG), including tumors diagnosed in the brainstem (diffuse intrinsic pontine glioma; DIPG), are uniformly fatal brain tumors that lack effective treatment. Analysis of CRISPR/Cas9 loss-of-function gene deletion screens identified PIK3CA and MTOR as targetable molecular dependencies across patient derived models of DIPG, highlighting the therapeutic potential of the blood-brain barrier-penetrant PI3K/Akt/mTOR inhibitor, paxalisib. At the human-equivalent maximum tolerated dose, mice treated with paxalisib experienced systemic glucose feedback and increased insulin levels commensurate with patients using PI3K inhibitors. To exploit genetic dependence and overcome resistance while maintaining compliance and therapeutic benefit, we combined paxalisib with the antihyperglycemic drug metformin. Metformin restored glucose homeostasis and decreased phosphorylation of the insulin receptor in vivo, a common mechanism of PI3K-inhibitor resistance, extending survival of orthotopic models. DIPG models treated with paxalisib increased calcium-activated PKC signaling. The brain penetrant PKC inhibitor enzastaurin, in combination with paxalisib, synergistically extended the survival of multiple orthotopic patient-derived and immunocompetent syngeneic allograft models; benefits potentiated in combination with metformin and standard-of-care radiotherapy. Therapeutic adaptation was assessed using spatial transcriptomics and ATAC-Seq, identifying changes in myelination and tumor immune microenvironment crosstalk. Collectively, this study has identified what we believe to be a clinically relevant DIPG therapeutic combinational strategy.

UBE2N is essential for maintenance of skin homeostasis and suppression of inflammation.

UBE2N, a Lys63-ubiquitin conjugating enzyme, plays critical roles in embryogenesis and immune system development and function. However, its roles in adult epithelial tissue homeostasis and pathogenesis are unclear. We generated conditional mouse models that deleted Ube2n in skin cells in a temporally and spatially controlled manner. We found that Ube2n-knockout (KO) in the adult skin keratinocytes induced a range of inflammatory skin defects characteristic of psoriatic and actinic keratosis. These included eczematous inflammation, epidermal and dermal thickening, parakeratosis, and increased immune cell infiltration, as well as signs of edema and blistering. Single cell transcriptomic analyses and RT-qPCR showed that Ube2n KO keratinocytes expressed elevated myeloid cell chemo-attractants such as Cxcl1 and Cxcl2 and decreased the homeostatic T lymphocyte chemo-attractant, Ccl27a. Consistently, the infiltrating immune cells of Ube2n-KO skin were predominantly myeloid-derived cells including neutrophils and M1-like macrophages that were highly inflammatory, as indicated by expression of Il1ß and Il24. Pharmacological blockade of the IL-1 receptor associated kinases (IRAK1/4) alleviated eczema, epidermal and dermal thickening, and immune infiltration of the Ube2n mutant skin. Together, these findings highlight a key role of keratinocyte-UBE2N in maintenance of epidermal homeostasis and skin immunity and identify IRAK1/4 as potential therapeutic target for inflammatory skin disorders.

Comprehensive characterization of pathogenic synovial fluid extracellular vesicles from knee osteoarthritis.

Synovial fluid (SF) extracellular vesicles (EVs) play a pathogenic role in osteoarthritis (OA). However, the surface markers, cell and tissue origins, and effectors of these EVs are largely unknown. We found that SF EVs contained 692 peptides that were positively associated with knee radiographic OA severity; 57.4% of these pathogenic peptides were from 46 proteins of the immune system, predominantly the innate immune system. CSPG4, BGN, NRP1, and CD109 are the major surface markers of pathogenic SF EVs. Genes encoding surface marker CSPG4 and CD109 were highly expressed by chondrocytes from damaged cartilage, while VISG4, MARCO, CD163 and NRP1 were enriched in the synovial immune cells. The frequency of CSPG4+ and VSIG4+ EV subpopulations in OA SF was high. We conclude that pathogenic SF EVs carry knee OA severity-associated proteins and specific surface markers, which could be developed as a new source of diagnostic biomarkers or therapeutic targets in OA.

Ataxia-telangiectasia mutated ( Atm ) disruption sensitizes spatially-directed H3.3K27M/TP53 diffuse midline gliomas to radiation therapy.

Diffuse midline gliomas (DMGs) are lethal brain tumors characterized by p53-inactivating mutations and oncohistone H3.3K27M mutations that rewire the cellular response to genotoxic stress, which presents therapeutic opportunities. We used RCAS/tv-a retroviruses and Cre recombinase to inactivate p53 and induce K27M in the native H3f3a allele in a lineage- and spatially-directed manner, yielding primary mouse DMGs. Genetic or pharmacologic disruption of the DNA damage response kinase Ataxia-telangiectasia mutated (ATM) enhanced the efficacy of focal brain irradiation, extending mouse survival. This finding suggests that targeting ATM will enhance the efficacy of radiation therapy for p53-mutant DMG but not p53-wildtype DMG. We used spatial in situ transcriptomics and an allelic series of primary murine DMG models with different p53 mutations to identify transactivation-independent p53 activity as a key mediator of such radiosensitivity. These studies deeply profile a genetically faithful and versatile model of a lethal brain tumor to identify resistance mechanisms for a therapeutic strategy currently in clinical trials.

Age-Specific Diagnostic Classification of ASD Using Deep Learning Approaches.

Autism Spectrum Disorder (ASD) is a highly heterogeneous condition, due to high variance in its etiology, comorbidity, pathogenesis, severity, genetics, and brain functional connectivity (FC). This makes it devoid of any robust universal biomarker. This study aims to analyze the role of age and multivariate patterns in brain FC and their accountability in diagnosing ASD by deep learning algorithms. We utilized functional magnetic resonance imaging data of three age groups (6 to 11, 11 to 18, and 6 to 18 years), available with public databases ABIDE-I and ABIDE-II, to discriminate between ASD and typically developing. The blood-oxygen-level dependent time series were extracted using the Gordon's, Harvard Oxford and Diedrichsen's atlases, over 236 regions of interest, as 236x236 sized FC matrices for each participant, with Pearson correlations. The feature sets, in the form of FC heat maps were computed with respect to each age group and were fed to a convolutional neural network, such as MobileNetV2 and DenseNet201 to build age-specific diagnostic models. The results revealed that DenseNet201 was able to adapt and extract better features from the heat maps, and hence returned better accuracy scores. The age-specific dataset, with participants of ages 6 to 11 years, performed best, followed by 11 to 18 years and 6 to 18 years, with accuracy scores of 72.19%, 71.88%, and 69.74% respectively, when tested using the DenseNet201. Our results suggest that age-specific diagnostic models are able to counter heterogeneity present in ASD, and that enables better discrimination.

Automated microarray platform for single-cell sorting and collection of lymphocytes following HIV reactivation.

A promising strategy to cure HIV-infected individuals is to use latency reversing agents (LRAs) to reactivate latent viruses, followed by host clearance of infected reservoir cells. However, reactivation of latent proviruses within infected cells is heterogeneous and often incomplete. This fact limits strategies to cure HIV which may require complete elimination of viable virus from all cellular reservoirs. For this reason, understanding the mechanism(s) of reactivation of HIV within cellular reservoirs is critical to achieve therapeutic success. Methodologies enabling temporal tracking of single cells as they reactivate followed by sorting and molecular analysis of those cells are urgently needed. To this end, microraft arrays were adapted to image T-lymphocytes expressing mCherry under the control of the HIV long terminal repeat (LTR) promoter, in response to the application of LRAs (prostratin, iBET151, and SAHA). In response to prostratin, iBET151, and SAHA, 30.5%, 11.2%, and 12.1% percentage of cells, respectively. The arrays enabled large numbers of single cells (>25,000) to be imaged over time. mCherry fluorescence quantification identified cell subpopulations with differing reactivation kinetics. Significant heterogeneity was observed at the single-cell level between different LRAs in terms of time to reactivation, rate of mCherry fluorescence increase upon reactivation, and peak fluorescence attained. In response to prostratin, subpopulations of T lymphocytes with slow and fast reactivation kinetics were identified. Single T-lymphocytes that were either fast or slow reactivators were sorted, and single-cell RNA-sequencing was performed. Different genes associated with inflammation, immune activation, and cellular and viral transcription factors were found.

20-αHydroxycholesterol, an oxysterol in human breast milk, reverses mouse neonatal white matter injury through Gli-dependent oligodendrogenesis.

White matter injuries (WMIs) are the leading cause of neurologic impairment in infants born premature. There are no treatment options available. The most common forms of WMIs in infants occur prior to the onset of normal myelination, making its pathophysiology distinctive, thus requiring a tailored approach to treatment. Neonates present a unique opportunity to repair WMIs due to a transient abundance of neural stem/progenitor cells (NSPCs) present in the germinal matrix with oligodendrogenic potential. We identified an endogenous oxysterol, 20-αHydroxycholesterol (20HC), in human maternal breast milk that induces oligodendrogenesis through a sonic hedgehog (shh), Gli-dependent mechanism. Following WMI in neonatal mice, injection of 20HC induced subventricular zone-derived oligodendrogenesis and improved myelination in the periventricular white matter, resulting in improved motor outcomes. Targeting the oligodendrogenic potential of postnatal NSPCs in neonates with WMIs may be further developed into a novel approach to mitigate this devastating complication of preterm birth.

The RNA helicase DDX39B activates FOXP3 RNA splicing to control T regulatory cell fate.

Genes associated with increased susceptibility to multiple sclerosis (MS) have been identified, but their functions are incompletely understood. One of these genes codes for the RNA helicase DExD/H-Box Polypeptide 39B (DDX39B), which shows genetic and functional epistasis with interleukin-7 receptor-α gene (IL7R) in MS-risk. Based on evolutionary and functional arguments, we postulated that DDX39B enhances immune tolerance thereby decreasing MS risk. Consistent with such a role we show that DDX39B controls the expression of many MS susceptibility genes and important immune-related genes. Among these we identified Forkhead Box P3 (FOXP3), which codes for the master transcriptional factor in CD4+/CD25+ T regulatory cells. DDX39B knockdown led to loss of immune-regulatory and gain of immune-effector expression signatures. Splicing of FOXP3 introns, which belong to a previously unrecognized type of introns with C-rich polypyrimidine tracts, was exquisitely sensitive to DDX39B levels. Given the importance of FOXP3 in autoimmunity, this work cements DDX39B as an important guardian of immune tolerance.

Sex-Biased Expression Is Associated With Chromatin State in Drosophila melanogaster and Drosophila simulans.

In Drosophila melanogaster and D. simulans head tissue, 60% of orthologous genes show evidence of sex-biased expression in at least one species. Of these, ∼39% (2,192) are conserved in direction. We hypothesize enrichment of open chromatin in the sex where we see expression bias and closed chromatin in the opposite sex. Male-biased orthologs are significantly enriched for H3K4me3 marks in males of both species (∼89% of male-biased orthologs vs. ∼76% of unbiased orthologs). Similarly, female-biased orthologs are significantly enriched for H3K4me3 marks in females of both species (∼90% of female-biased orthologs vs. ∼73% of unbiased orthologs). The sex-bias ratio in female-biased orthologs was similar in magnitude between the two species, regardless of the closed chromatin (H3K27me2me3) marks in males. However, in male-biased orthologs, the presence of H3K27me2me3 in both species significantly reduced the correlation between D. melanogaster sex-bias ratio and the D. simulans sex-bias ratio. Male-biased orthologs are enriched for evidence of positive selection in the D. melanogaster group. There are more male-biased genes than female-biased genes in both species. For orthologs with gains/losses of sex-bias between the two species, there is an excess of male-bias compared to female-bias, but there is no consistent pattern in the relationship between H3K4me3 or H3K27me2me3 chromatin marks and expression. These data suggest chromatin state is a component of the maintenance of sex-biased expression and divergence of sex-bias between species is reflected in the complexity of the chromatin status.

JAK-STAT activation contributes to cytotoxic T cell-mediated basal cell death in human chronic lung allograft dysfunction.

Chronic lung allograft dysfunction (CLAD) is the leading cause of death in lung transplant recipients. CLAD is characterized clinically by a persistent decline in pulmonary function and histologically by the development of airway-centered fibrosis known as bronchiolitis obliterans. There are no approved therapies to treat CLAD, and the mechanisms underlying its development remain poorly understood. We performed single-cell RNA-Seq and spatial transcriptomic analysis of explanted tissues from human lung recipients with CLAD, and we performed independent validation studies to identify an important role of Janus kinase-signal transducer and activator of transcription (JAK-STAT) signaling in airway epithelial cells that contributes to airway-specific alloimmune injury. Specifically, we established that activation of JAK-STAT signaling leads to upregulation of major histocompatibility complex 1 (MHC-I) in airway basal cells, an important airway epithelial progenitor population, which leads to cytotoxic T cell-mediated basal cell death. This study provides mechanistic insight into the cell-to-cell interactions driving airway-centric alloimmune injury in CLAD, suggesting a potentially novel therapeutic strategy for CLAD prevention or treatment.

Recent advances in targeting autophagy in cancer.

Autophagy is a cellular homeostasis mechanism that fuels the proliferation and survival of advanced cancers by degrading and recycling organelles and proteins. Preclinical studies have identified that within an established tumor, tumor cell autophagy and host cell autophagy conspire to support tumor growth. A growing body of evidence suggests that autophagy inhibition can augment the efficacy of chemotherapy, targeted therapy, or immunotherapy to enhance tumor shrinkage. First-generation autophagy inhibition trials in cancer using the lysosomal inhibitor hydroxychloroquine (HCQ) have produced mixed results but have guided the way for the development of more potent and specific autophagy inhibitors in clinical trials. In this review, we will discuss the role of autophagy in cancer, newly discovered molecular mechanisms of the autophagy pathway, the effects of autophagy modulation in cancer and host cells, and novel autophagy inhibitors that are entering clinical trials.

Ganglioglioma deep transcriptomics reveals primitive neuroectoderm neural precursor-like population.

Gangliogliomas are brain tumors composed of neuron-like and macroglia-like components that occur in children and young adults. Gangliogliomas are often characterized by a rare population of immature astrocyte-appearing cells expressing CD34, a marker expressed in the neuroectoderm (neural precursor cells) during embryogenesis. New insights are needed to refine tumor classification and to identify therapeutic approaches. We evaluated five gangliogliomas with single nucleus RNA-seq, cellular indexing of transcriptomes and epitopes by sequencing, and/or spatially-resolved RNA-seq. We uncovered a population of CD34+ neoplastic cells with mixed neuroectodermal, immature astrocyte, and neuronal markers. Gene regulatory network interrogation in these neuroectoderm-like cells revealed control of transcriptional programming by TCF7L2/MEIS1-PAX6 and SOX2, similar to that found during neuroectodermal/neural development. Developmental trajectory analyses place neuroectoderm-like tumor cells as precursor cells that give rise to neuron-like and macroglia-like neoplastic cells. Spatially-resolved transcriptomics revealed a neuroectoderm-like tumor cell niche with relative lack of vascular and immune cells. We used these high resolution results to deconvolute clinically-annotated transcriptomic data, confirming that CD34+ cell-associated gene programs associate with gangliogliomas compared to other glial brain tumors. Together, these deep transcriptomic approaches characterized a ganglioglioma cellular hierarchy-confirming CD34+ neuroectoderm-like tumor precursor cells, controlling transcription programs, cell signaling, and associated immune cell states. These findings may guide tumor classification, diagnosis, prognostication, and therapeutic investigations.

Automated microarray for single-cell sorting and collection of lymphocytes following HIV reactivation.

A promising strategy to cure HIV infected individuals is to use latency reversing agents (LRAs) to reactivate latent viruses, followed by host clearance of infected reservoir cells. However, reactivation of latent proviruses within infected cells is heterogeneous and often incomplete. This fact limits strategies to cure HIV which may require complete elimination of viable virus from all cellular reservoirs. For this reason, understanding the mechanism(s) of reactivation of HIV within cellular reservoirs is critical to achieve therapeutic success. Methodologies enabling temporal tracking of single cells as they reactivate followed by sorting and molecular analysis of those cells are urgently needed. To this end, microraft arrays were adapted to image T-lymphocytes expressing mCherry under the control of the HIV long terminal repeat (LTR) promoter, in response to the application of various LRAs (prostratin, iBET151, and SAHA). In response to prostratin, iBET151, and SAHA, 30.5 %, 11.2 %, and 12.1 % percentage of cells respectively, reactivated similar to that observed in other experimental systems. The arrays enabled large numbers of single cells (>25,000) to be imaged over time. mCherry fluorescence quantification identified cell subpopulations with differing reactivation kinetics. Significant heterogeneity was observed at the single cell level between different LRAs in terms of time to reactivation, rate of mCherry fluorescence increase upon reactivation, and peak fluorescence attained. In response to prostratin, subpopulations of T lymphocytes with slow and fast reactivation kinetics were identified. Single T-lymphocytes that were either fast or slow reactivators were sorted, and single-cell RNA-sequencing was performed. Different genes associated with inflammation, immune activation, and cellular and viral transcription factors were found. These results advance our conceptual understanding of HIV reactivation dynamics at the single-cell level toward a cure for HIV.