Multiomic analysis reveals cellular and epigenetic plasticity in intestinal pouches of ulcerative colitis patients.

Background & Aims: Total proctocolectomy with ileal pouch anal anastomosis (IPAA) is the standard of care for patients with severe treatment resistant ulcerative colitis (UC). Despite improvements in patient outcomes, about 50% of patients will develop inflammation of the pouch within 1-2 years following surgery. Establishment of UC pouches is associated with profound histological changes of the mucosa. A detailed characterization of these changes on a cellular and molecular level is crucial for an improved understanding of pouch physiology and diseases management. Methods: We generated cell-type-resolved transcriptional and epigenetic atlases of UC pouches using scRNA-seq and scATAC-seq data from paired biopsy samples from the ileal pouch and ileal segment above the pouch (pre-pouch) of UC-IPAA patients (n=6, female=2) without symptoms. We also collected data from paired biopsies of the terminal ileum (TI) and ascending colon (AC) from healthy controls (n=6, female=3). Results: We identified novel populations of colon-like absorptive and secretory epithelial cells, constituting a significant proportion of the epithelial cell fraction in the pouch but not in matched pre-pouch samples. Pouch-specific enterocytes expressed colon-specific genes, including CEACAM5, CA2. However, in contrast to normal colonic epithelium, these cells also expressed a range of inflammatory and secretory genes, similar to previously detected gene expression signatures in IBD patients. Comparison to longitudinal bulk RNA-seq data from UC pouches demonstrated that colon-like epithelial cells are present early after pouch functionalization and independently of subsequent pouchitis. Finally, single cell chromatin accessibility revealed activation colonic transcriptional regulators, including CDX1, NFIA, and EHF. Conclusion: UC pouches are characterized by partial colonic metaplasia of the epithelium. These data constitute a resource of transcriptomic and epigenetic signatures of cell populations in the pouch and provide an anchor for understanding the underlying molecular mechanisms of pouchitis.

Ribonucleotide reductase regulatory subunit M2 drives glioblastoma TMZ resistance through modulation of dNTP production.

During therapy, adaptations driven by cellular plasticity are partly responsible for driving the inevitable recurrence of glioblastoma (GBM). To investigate plasticity-induced adaptation during standard-of-care chemotherapy temozolomide (TMZ), we performed in vivo single-cell RNA sequencing in patient-derived xenograft (PDX) tumors of GBM before, during, and after therapy. Comparing single-cell transcriptomic patterns identified distinct cellular populations present during TMZ therapy. Of interest was the increased expression of ribonucleotide reductase regulatory subunit M2 (RRM2), which we found to regulate dGTP and dCTP production vital for DNA damage response during TMZ therapy. Furthermore, multidimensional modeling of spatially resolved transcriptomic and metabolomic analysis in patients' tissues revealed strong correlations between RRM2 and dGTP. This supports our data that RRM2 regulates the demand for specific dNTPs during therapy. In addition, treatment with the RRM2 inhibitor 3-AP (Triapine) enhances the efficacy of TMZ therapy in PDX models. We present a previously unidentified understanding of chemoresistance through critical RRM2-mediated nucleotide production.

Exploring the tumor micro-environment in primary and metastatic tumors of different ovarian cancer histotypes.

Ovarian cancer is a highly heterogeneous disease consisting of at least five different histological subtypes with varying clinical features, cells of origin, molecular composition, risk factors, and treatments. While most single-cell studies have focused on High grade serous ovarian cancer, a comprehensive landscape of the constituent cell types and their interactions within the tumor microenvironment are yet to be established in the different ovarian cancer histotypes. Further characterization of tumor progression, metastasis, and various histotypes are also needed to connect molecular signatures to pathological grading for personalized diagnosis and tailored treatment. In this study, we leveraged high-resolution single-cell RNA sequencing technology to elucidate the cellular compositions on 21 solid tumor samples collected from 12 patients with six ovarian cancer histotypes and both primary (ovaries) and metastatic (omentum, rectum) sites. The diverse collection allowed us to deconstruct the histotypes and tumor site-specific expression patterns of cells in the tumor, and identify key marker genes and ligand-receptor pairs that are active in the ovarian tumor microenvironment. Our findings can be used in improving precision disease stratification and optimizing treatment options.

A molecular atlas of the human postmenopausal fallopian tube and ovary from single-cell RNA and ATAC sequencing.

As part of the Human Cell Atlas Initiative, our goal is to generate single-cell transcriptomics (single-cell RNA sequencing [scRNA-seq], 86,708 cells) and regulatory (single-cell assay on transposase accessible chromatin sequencing [scATAC-seq], 59,830 cells) profiles of the normal postmenopausal ovary and fallopian tube (FT). The FT contains 11 major cell types, and the ovary contains 6. The dominating cell type in the FT and ovary is the stromal cell, which expresses aging-associated genes. FT epithelial cells express multiple ovarian cancer risk-associated genes (CCDC170, RND3, TACC2, STK33, and ADGB) and show active communication between fimbrial epithelial cells and ovarian stromal cells. Integrated single-cell transcriptomics and chromatin accessibility data show that the regulatory landscape of the fimbriae is different from other anatomic regions. Cell types with similar gene expression in the FT display transcriptional profiles. These findings allow us to disentangle the cellular makeup of the postmenopausal FT and ovary, advancing our knowledge of gynecologic diseases in menopause.

Creating patient-specific vein models to characterize wall shear stress in hemodialysis population.

End-Stage Renal Disease (ESRD) patients require arteriovenous fistulas (AVF) that allow a mature vein to withstand hemodialysis. Unfortunately, venous thrombosis and stenosis in the cephalic vein arch after AVF placement is common and heavily influenced by hemodynamics. To better assess forces and flow behavior in the cephalic arch, we have built patient-specific millifluidic models that allow us to explore the complex interplay between patient-specific vein geometry and fluctuating hemodynamics. These 3D models were created from patient-specific intravascular ultrasound and venogram images obtained three- and twelve-months post AVF creation and fabricated into soft elastomer-based millifluidic devices. Geometric validation of fabricated phantom millifluidic device shows successful replication of original computational 3D model. Millifluidic devices were perfused with a blood-mimicking fluid containing fluorescent tracer beads under steady-state physiologic cephalic vein flow conditions (20 mL/min). Particle image velocimetry was employed to calculate wall shear stress (WSS) across the cephalic arches. Experimental WSS profile evaluation reveals that the physiologic cephalic arch model yields WSS values within physiologic range [76-760 mPa]. Moreover, upon comparing WSS profiles across all models, it is noticeable that WSS values increase as vein diameter decreases, which further supports employed experimental and analysis strategy. The presented millifluidic devices show promise for experimental WSS characterization under pathologic flow conditions to contrast from calculated physiologic hemodynamics and better understand WSS influence on thrombosis and stenosis in hemodialysis patients.

Hedgehog signaling activates a mammalian heterochronic gene regulatory network controlling differentiation timing across lineages.

Many developmental signaling pathways have been implicated in lineage-specific differentiation; however, mechanisms that explicitly control differentiation timing remain poorly defined in mammals. We report that murine Hedgehog signaling is a heterochronic pathway that determines the timing of progenitor differentiation. Hedgehog activity was necessary to prevent premature differentiation of second heart field (SHF) cardiac progenitors in mouse embryos, and the Hedgehog transcription factor GLI1 was sufficient to delay differentiation of cardiac progenitors in vitro. GLI1 directly activated a de novo progenitor-specific network in vitro, akin to that of SHF progenitors in vivo, which prevented the onset of the cardiac differentiation program. A Hedgehog signaling-dependent active-to-repressive GLI transition functioned as a differentiation timer, restricting the progenitor network to the SHF. GLI1 expression was associated with progenitor status across germ layers, and it delayed the differentiation of neural progenitors in vitro, suggesting a broad role for Hedgehog signaling as a heterochronic pathway.

Engineered exosomes targeting MYC reverse the proneural-mesenchymal transition and extend survival of glioblastoma.

Dysregulated Myc signaling is a key oncogenic pathway in glioblastoma multiforme (GBM). Yet, effective therapeutic targeting of Myc continues to be challenging. Here, we demonstrate that exosomes generated from human bone marrow mesenchymal stem cells (MSCs) engineered to encapsulate siRNAs targeting Myc (iExo-Myc) localize to orthotopic GBM tumors in mice. Treatment of late stage GBM tumors with iExo-Myc inhibits proliferation and angiogenesis, suppresses tumor growth, and extends survival. Transcriptional profiling of tumors reveals that the mesenchymal transition and estrogen receptor signaling pathways are impacted by Myc inhibition. Single nuclei RNA sequencing (snRNA-seq) shows that iExo-Myc treatment induces transcriptional repression of multiple growth factor and interleukin signaling pathways, triggering a mesenchymal to proneural transition and shifting the cellular landscape of the tumor. These data confirm that Myc is an effective anti-glioma target and that iExo-Myc offers a feasible, readily translational strategy to inhibit challenging oncogene targets for the treatment of brain tumors.

High-throughput RNA sequencing of paraformaldehyde-fixed single cells.

Single-cell transcriptomic studies that require intracellular protein staining, rare cell sorting, or inactivation of infectious pathogens are severely limited. This is because current high-throughput single-cell RNA sequencing methods are either incompatible with or necessitate laborious sample preprocessing for paraformaldehyde treatment, a common tissue and cell fixation and preservation technique. Here we present FD-seq (Fixed Droplet RNA sequencing), a high-throughput method for droplet-based RNA sequencing of paraformaldehyde-fixed, permeabilized and sorted single cells. We show that FD-seq preserves the RNA integrity and relative gene expression levels after fixation and permeabilization. Furthermore, FD-seq can detect a higher number of genes and transcripts than methanol fixation. We first apply FD-seq to analyze a rare subpopulation of cells supporting lytic reactivation of the human tumor virus KSHV, and identify TMEM119 as a potential host factor that mediates viral reactivation. Second, we find that infection with the human betacoronavirus OC43 leads to upregulation of pro-inflammatory pathways in cells that are exposed to the virus but fail to express high levels of viral genes. FD-seq thus enables integrating phenotypic with transcriptomic information in rare cell subpopulations, and preserving and inactivating pathogenic samples.

Characterizing the tumor microenvironment of metastatic ovarian cancer by single-cell transcriptomics.

Understanding the cellular composition of the tumor microenvironment and the interactions of the cells is essential to the development of successful immunotherapies in cancer. We perform single-cell RNA sequencing (scRNA-seq) of 9,885 cells isolated from the omentum in 6 patients with ovarian cancer and identify 9 major cell types, including cancer, stromal, and immune cells. Transcriptional analysis of immune cells stratifies our patient samples into 2 groups: (1) high T cell infiltration (high Tinf) and (2) low T cell infiltration (low Tinf). TOX-expressing resident memory CD8+ T (CD8+ Trm) and granulysin-expressing CD4+ T cell clusters are enriched in the high Tinf group. Concurrently, we find unique plasmablast and plasma B cell clusters, and finally, NR1H2+IRF8+ and CD274+ macrophage clusters, suggesting an anti-tumor response in the high Tinf group. Our scRNA-seq study of metastatic tumor samples provides important insights in elucidating the immune response within ovarian tumors.

Fixed single-cell RNA sequencing for understanding virus infection and host response.

Single-cell transcriptomic studies that require intracellular protein staining, rare cell sorting, or inactivation of infectious pathogens are severely limited because current high-throughput RNA sequencing methods are incompatible with paraformaldehyde treatment, a common tissue and cell fixation and preservation technique. Here we present FD-seq, a high-throughput method for droplet-based RNA sequencing of paraformaldehyde-fixed, stained and sorted single-cells. We show that FD-seq preserves the mRNA integrity and relative abundances during fixation and subsequent cell retrieval. Furthermore, FD-seq detects a higher number of genes and transcripts than methanol fixation. We applied FD-seq to investigate two important questions in Virology. First, by analyzing a rare population of cells supporting lytic reactivation of the human tumor virus KSHV, we identified TMEM119 as a host factor that mediates viral reactivation. Second, we found that upon infection with the betacoronavirus OC43, which causes the common cold and is a close relative of SARS-CoV-2, pro-inflammatory pathways are primarily upregulated in lowly-infected cells that are exposed to the virus but fail to express high levels of viral genes. FD-seq thus enables integrating phenotypic with transcriptomic information in rare cell populations, and preserving and inactivating pathogenic samples that cannot be handled under regular biosafety measures.

Hedgehog-FGF signaling axis patterns anterior mesoderm during gastrulation.

The mechanisms used by embryos to pattern tissues across their axes has fascinated developmental biologists since the founding of embryology. Here, using single-cell technology, we interrogate complex patterning defects and define a Hedgehog (Hh)-fibroblast growth factor (FGF) signaling axis required for anterior mesoderm lineage development during gastrulation. Single-cell transcriptome analysis of Hh-deficient mesoderm revealed selective deficits in anterior mesoderm populations, culminating in defects to anterior embryonic structures, including the pharyngeal arches, heart, and anterior somites. Transcriptional profiling of Hh-deficient mesoderm during gastrulation revealed disruptions to both transcriptional patterning of the mesoderm and FGF signaling for mesoderm migration. Mesoderm-specific Fgf4/Fgf8 double-mutants recapitulated anterior mesoderm defects and Hh-dependent GLI transcription factors modulated enhancers at FGF gene loci. Cellular migration defects during gastrulation induced by Hh pathway antagonism were mitigated by the addition of FGF4 protein. These findings implicate a multicomponent signaling hierarchy activated by Hh ligands from the embryonic node and executed by FGF signals in nascent mesoderm to control anterior mesoderm patterning.