Modeling blood-brain barrier formation and cerebral cavernous malformations in human PSC-derived organoids.

The human blood-brain barrier (hBBB) is a highly specialized structure that regulates passage across blood and central nervous system (CNS) compartments. Despite its critical physiological role, there are no reliable in vitro models that can mimic hBBB development and function. Here, we constructed hBBB assembloids from brain and blood vessel organoids derived from human pluripotent stem cells. We validated the acquisition of blood-brain barrier (BBB)-specific molecular, cellular, transcriptomic, and functional characteristics and uncovered an extensive neuro-vascular crosstalk with a spatial pattern within hBBB assembloids. When we used patient-derived hBBB assembloids to model cerebral cavernous malformations (CCMs), we found that these assembloids recapitulated the cavernoma anatomy and BBB breakdown observed in patients. Upon comparison of phenotypes and transcriptome between patient-derived hBBB assembloids and primary human cavernoma tissues, we uncovered CCM-related molecular and cellular alterations. Taken together, we report hBBB assembloids that mimic the core properties of the hBBB and identify a potentially underlying cause of CCMs.

The adrenal stress response is an essential host response against therapy-induced lethal immune activation.

Cytokine release syndrome (CRS) is a systemic inflammatory syndrome associated with infection- or drug-induced T cell activation and can cause multiple organ failure and even death. Because current treatments are ineffective in some patients with severe CRS, we set out to identify risk factors and mechanisms behind severe CRS that might lead to preventive therapies and better clinical outcomes in patients. In mice, we found that deficiency in the adrenal stress response-with similarities to such in patients called relative adrenal insufficiency (RAI)-conferred a high risk for lethal CRS. Mice treated with CD3 antibodies were protected against lethal CRS by the production of glucocorticoids (GC) induced by the adrenal stress response in a manner dependent on the scavenger receptor B1 (SR-BI), a receptor for high-density lipoprotein (HDL). Mice with whole-body or adrenal gland-specific SR-BI deficiency exhibited impaired GC production, more severe CRS, and increased mortality in response to CD3 antibodies. Pretreatment with a low dose of GC effectively suppressed the development of CRS and rescued survival in SR-BI-deficient mice without compromising T cell function through apoptosis. Our findings suggest that RAI may be a risk factor for therapy-induced CRS and that pretreating RAI patients with GC may prevent lethal CRS.

Inhibition of CSPG receptor PTPσ promotes migration of newly born neuroblasts, axonal sprouting, and recovery from stroke.

In addition to neuroprotective strategies, neuroregenerative processes could provide targets for stroke recovery. However, the upregulation of inhibitory chondroitin sulfate proteoglycans (CSPGs) impedes innate regenerative efforts. Here, we examine the regulatory role of PTPσ (a major proteoglycan receptor) in dampening post-stroke recovery. Use of a receptor modulatory peptide (ISP) or Ptprs gene deletion leads to increased neurite outgrowth and enhanced NSCs migration upon inhibitory CSPG substrates. Post-stroke ISP treatment results in increased axonal sprouting as well as neuroblast migration deeply into the lesion scar with a transcriptional signature reflective of repair. Lastly, peptide treatment post-stroke (initiated acutely or more chronically at 7 days) results in improved behavioral recovery in both motor and cognitive functions. Therefore, we propose that CSPGs induced by stroke play a predominant role in the regulation of neural repair and that blocking CSPG signaling pathways will lead to enhanced neurorepair and functional recovery in stroke.

Multi-omics to characterize the functional relationships of R-loops with epigenetic modifications, RNAPII transcription and gene expression.

Abnormal accumulation of R-loops results in replication stress, genome instability, chromatin alterations and gene silencing. Little research has been done to characterize functional relationships among R-loops, histone marks, RNA polymerase II (RNAPII) transcription and gene regulation. We built extremely randomized trees (ETs) models to predict the genome-wide R-loops using RNAPII and multiple histone modifications chromatin immunoprecipitation (ChIP)-seq, DNase-seq, Global Run-On sequencing (GRO-seq) and R-loop profiling data. We compared the performance of ET models to multiple machine learning approaches, and the proposed ET models achieved the best and extremely robust performances. Epigenetic profiles are highly predictive of R-loops genome-widely and they are strongly associated with R-loop formation. In addition, the presence of R-loops is significantly correlated with RNAPII transcription activity, H3K4me3 and open chromatin around the transcription start site, and H3K9me1 and H3K9me3 around the transcription termination site. RNAPII pausing defects were correlated with 5'R-loops accumulation, and transcriptional termination defects and read-throughs were correlated with 3'R-loops accumulation. Furthermore, we found driver genes with 5'R-loops and RNAPII pausing defects express significantly higher and genes with 3'R-loops and read-through transcription express significantly lower than genes without R-loops. These driver genes are enriched with chromosomal instability, Hippo-Merlin signaling Dysregulation, DNA damage response and TGF-ß pathways, indicating R-loops accumulating at the 5' end of genes play oncogenic roles, whereas at the 3' end of genes play tumor-suppressive roles in tumorigenesis.

Synergistic anti-tumor efficacy of mutant isocitrate dehydrogenase 1 inhibitor SYC-435 with standard therapy in patient-derived xenograft mouse models of glioma.

Clinical outcomes in patients with WHO grade II/III astrocytoma, oligodendroglioma or secondary glioblastoma remain poor. Isocitrate dehydrogenase 1 (IDH1) is mutated in > 70% of these tumors, making it an attractive therapeutic target. To determine the efficacy of our newly developed mutant IDH1 inhibitor, SYC-435 (1-hydroxypyridin-2-one), we treated orthotopic glioma xenograft model (IC-BT142AOA) carrying R132H mutation and our newly established orthotopic patient-derived xenograft (PDX) model of recurrent anaplastic oligoastrocytoma (IC-V0914AOA) bearing R132C mutation. In addition to suppressing IDH1 mutant cell proliferation in vitro, SYC-435 (15 mg/kg, daily x 28 days) synergistically prolonged animal survival times with standard therapies (Temozolomide + fractionated radiation) mediated by reduction of H3K4/H3K9 methylation and expression of mitochondrial DNA (mtDNA)-encoded molecules. Furthermore, RNA-seq of the remnant tumors identified genes (MYO1F, CTC1 and BCL9) and pathways (base excision repair, TCA cycle II, sirtuin signaling, protein kinase A, eukaryotic initiation factor 2 and α-adrenergic signaling) as mediators of therapy resistance. Our data demonstrated the efficacy SYC-435 in targeting IDH1 mutant gliomas when combined with standard therapy and identified a novel set of genes that should be prioritized for future studies to overcome SYC-435 resistance.

Spatial Dissection of Invasive Front from Tumor Mass Enables Discovery of Novel microRNA Drivers of Glioblastoma Invasion.

Diffuse invasion is the primary cause of treatment failure of glioblastoma (GBM). Previous studies on GBM invasion have long been forced to use the resected tumor mass cells. Here, a strategy to reliably isolate matching pairs of invasive (GBMINV ) and tumor core (GBMTC ) cells from the brains of 6 highly invasive patient-derived orthotopic models is described. Direct comparison of these GBMINV and GBMTC cells reveals a significantly elevated invasion capacity in GBMINV cells, detects 23/768 miRNAs over-expressed in the GBMINV cells (miRNAINV ) and 22/768 in the GBMTC cells (miRNATC ), respectively. Silencing the top 3 miRNAsINV (miR-126, miR-369-5p, miR-487b) successfully blocks invasion of GBMINV cells in vitro and in mouse brains. Integrated analysis with mRNA expression identifies miRNAINV target genes and discovers KCNA1 as the sole common computational target gene of which 3 inhibitors significantly suppress invasion in vitro. Furthermore, in vivo treatment with 4-aminopyridine (4-AP) effectively eliminates GBM invasion and significantly prolongs animal survival times (P = 0.035). The results highlight the power of spatial dissection of functionally accurate GBMINV and GBMTC cells in identifying novel drivers of GBM invasion and provide strong rationale to support the use of biologically accurate starting materials in understanding cancer invasion and metastasis.

Tumor Microenvironment-Derived R-spondins Enhance Antitumor Immunity to Suppress Tumor Growth and Sensitize for Immune Checkpoint Blockade Therapy.

Natural killer (NK) cells and T cells are key effectors of antitumor immune responses and major targets of checkpoint inhibitors. In multiple cancer types, we find that the expression of Wnt signaling potentiator R-spondin genes (e.g., RSPO3) is associated with favorable prognosis and positively correlates with gene signatures of both NK cells and T cells. Although endothelial cells and cancer-associated fibroblasts comprise the R-spondin 3-producing cells, NK cells and T cells correspondingly express the R-spondin 3 receptor LGR6 within the tumor microenvironment (TME). Exogenous expression or intratumor injection of R-spondin 3 in tumors enhanced the infiltration and function of cytotoxic effector cells, which led to tumor regression. NK cells and CD8+ T cells independently and cooperatively contributed to R-spondin 3-induced control of distinct tumor types. The effect of R-spondin 3 was mediated in part through upregulation of MYC and ribosomal biogenesis. Importantly, R-spondin 3 expression enhanced tumor sensitivity to anti-PD-1 therapy, thereby highlighting new therapeutic avenues. SIGNIFICANCE: Our study identifies novel targets in enhancing antitumor immunity and sensitizing immune checkpoint inhibition, which provides a rationale for developing new immunotherapies against cancers. It also offers mechanistic insights on Wnt signaling-mediated modulation of anticancer immunity in the TME and implications for a putative R-spondin-LGR6 axis in regulating NK-cell biology. This article is highlighted in the In This Issue feature, p. 2945.

In situ mapping identifies distinct vascular niches for myelopoiesis.

In contrast to nearly all other tissues, the anatomy of cell differentiation in the bone marrow remains unknown. This is owing to a lack of strategies for examining myelopoiesis-the differentiation of myeloid progenitors into a large variety of innate immune cells-in situ in the bone marrow. Such strategies are required to understand differentiation and lineage-commitment decisions, and to define how spatial organizing cues inform tissue function. Here we develop approaches for imaging myelopoiesis in mice, and generate atlases showing the differentiation of granulocytes, monocytes and dendritic cells. The generation of granulocytes and dendritic cells-monocytes localizes to different blood-vessel structures known as sinusoids, and displays lineage-specific spatial and clonal architectures. Acute systemic infection with Listeria monocytogenes induces lineage-specific progenitor clusters to undergo increased self-renewal of progenitors, but the different lineages remain spatially separated. Monocyte-dendritic cell progenitors (MDPs) map with nonclassical monocytes and conventional dendritic cells; these localize to a subset of blood vessels expressing a major regulator of myelopoiesis, colony-stimulating factor 1 (CSF1, also known as M-CSF)1. Specific deletion of Csf1 in endothelium disrupts the architecture around MDPs and their localization to sinusoids. Subsequently, there are fewer MDPs and their ability to differentiate is reduced, leading to a loss of nonclassical monocytes and dendritic cells during both homeostasis and infection. These data indicate that local cues produced by distinct blood vessels are responsible for the spatial organization of definitive blood cell differentiation.

Patient-Derived Orthotopic Xenograft (PDOX) Mouse Models of Primary and Recurrent Meningioma.

BACKGROUND: Meningiomas constitute one-third of all primary brain tumors. Although typically benign, about 20% of these tumors recur despite surgery and radiation, and may ultimately prove fatal. There are currently no effective chemotherapies for meningioma. We, therefore, set out to develop patient-derived orthotopic xenograft (PDOX) mouse models of human meningioma using tumor. METHOD: Of nine patients, four had World Health Organization (WHO) grade I tumors, five had WHO grade II tumors, and in this second group two patients also had recurrent (WHO grade III) meningioma. We also classified the tumors according to our recently developed molecular classification system (Types A, B, and C, with C being the most aggressive). We transplanted all 11 surgical samples into the skull base of immunodeficient (SCID) mice. Only the primary and recurrent tumor cells from one patient-both molecular Type C, despite being WHO grades II and III, respectively-led to the formation of meningioma in the resulting mouse models. We characterized the xenografts by histopathology and RNA-seq and compared them with the original tumors. We performed an in vitro drug screen using 60 anti-cancer drugs followed by in vivo validation. RESULTS: The PDOX models established from the primary and recurrent tumors from patient K29 (K29P-PDOX and K29R-PDOX, respectively) replicated the histopathology and key gene expression profiles of the original samples. Although these xenografts could not be subtransplanted, the cryopreserved primary tumor cells were able to reliably generate PDOX tumors. Drug screening in K29P and K29R tumor cell lines revealed eight compounds that were active on both tumors, including three histone deacetylase (HDAC) inhibitors. We tested the HDAC inhibitor Panobinostat in K29R-PDOX mice, and it significantly prolonged mouse survival (p < 0.05) by inducing histone H3 acetylation and apoptosis. CONCLUSION: Meningiomas are not very amenable to PDOX modeling, for reasons that remain unclear. Yet at least some of the most malignant tumors can be modeled, and cryopreserved primary tumor cells can create large panels of tumors that can be used for preclinical drug testing.

Network-based prediction of drug-target interactions using an arbitrary-order proximity embedded deep forest.

MOTIVATION: Systematic identification of molecular targets among known drugs plays an essential role in drug repurposing and understanding of their unexpected side effects. Computational approaches for prediction of drug-target interactions (DTIs) are highly desired in comparison to traditional experimental assays. Furthermore, recent advances of multiomics technologies and systems biology approaches have generated large-scale heterogeneous, biological networks, which offer unexpected opportunities for network-based identification of new molecular targets among known drugs. RESULTS: In this study, we present a network-based computational framework, termed AOPEDF, an arbitrary-order proximity embedded deep forest approach, for prediction of DTIs. AOPEDF learns a low-dimensional vector representation of features that preserve arbitrary-order proximity from a highly integrated, heterogeneous biological network connecting drugs, targets (proteins) and diseases. In total, we construct a heterogeneous network by uniquely integrating 15 networks covering chemical, genomic, phenotypic and network profiles among drugs, proteins/targets and diseases. Then, we build a cascade deep forest classifier to infer new DTIs. Via systematic performance evaluation, AOPEDF achieves high accuracy in identifying molecular targets among known drugs on two external validation sets collected from DrugCentral [area under the receiver operating characteristic curve (AUROC) = 0.868] and ChEMBL (AUROC = 0.768) databases, outperforming several state-of-the-art methods. In a case study, we showcase that multiple molecular targets predicted by AOPEDF are associated with mechanism-of-action of substance abuse disorder for several marketed drugs (such as aripiprazole, risperidone and haloperidol). AVAILABILITY AND IMPLEMENTATION: Source code and data can be downloaded from https://github.com/ChengF-Lab/AOPEDF.

Differential regulatory network-based quantification and prioritization of key genes underlying cancer drug resistance based on time-course RNA-seq data.

Drug resistance is a major cause for the failure of cancer chemotherapy or targeted therapy. However, the molecular regulatory mechanisms controlling the dynamic evolvement of drug resistance remain poorly understood. Thus, it is important to develop methods for identifying key gene regulatory mechanisms of the resistance to specific drugs. In this study, we developed a data-driven computational framework, DryNetMC, using a differential regulatory network-based modeling and characterization strategy to quantify and prioritize key genes underlying cancer drug resistance. The DryNetMC does not only infer gene regulatory networks (GRNs) via an integrated approach, but also characterizes and quantifies dynamical network properties for measuring node importance. We used time-course RNA-seq data from glioma cells treated with dbcAMP (a cAMP activator) as a realistic case to reconstruct the GRNs for sensitive and resistant cells. Based on a novel node importance index that comprehensively quantifies network topology, network entropy and expression dynamics, the top ranked genes were verified to be predictive of the drug sensitivities of different glioma cell lines, in comparison with other existing methods. The proposed method provides a quantitative approach to gain insights into the dynamic adaptation and regulatory mechanisms of cancer drug resistance and sheds light on the design of novel biomarkers or targets for predicting or overcoming drug resistance.

Single-Cell Transcriptomics in Medulloblastoma Reveals Tumor-Initiating Progenitors and Oncogenic Cascades during Tumorigenesis and Relapse.

Progenitor heterogeneity and identities underlying tumor initiation and relapse in medulloblastomas remain elusive. Utilizing single-cell transcriptomic analysis, we demonstrated a developmental hierarchy of progenitor pools in Sonic Hedgehog (SHH) medulloblastomas, and identified OLIG2-expressing glial progenitors as transit-amplifying cells at the tumorigenic onset. Although OLIG2+ progenitors become quiescent stem-like cells in full-blown tumors, they are highly enriched in therapy-resistant and recurrent medulloblastomas. Depletion of mitotic Olig2+ progenitors or Olig2 ablation impeded tumor initiation. Genomic profiling revealed that OLIG2 modulates chromatin landscapes and activates oncogenic networks including HIPPO-YAP/TAZ and AURORA-A/MYCN pathways. Co-targeting these oncogenic pathways induced tumor growth arrest. Together, our results indicate that glial lineage-associated OLIG2+ progenitors are tumor-initiating cells during medulloblastoma tumorigenesis and relapse, suggesting OLIG2-driven oncogenic networks as potential therapeutic targets.