Stem Cell Transplantation As A Dynamical System: Are Clinical Outcomes Deterministic?

Outcomes in stem cell transplantation (SCT) are modeled using probability theory. However the clinical course following SCT appears to demonstrate many characteristics of dynamical systems, especially when outcomes are considered in the context of immune reconstitution. Dynamical systems tend to evolve over time according to mathematically determined rules. Characteristically, the future states of the system are predicated on the states preceding them, and there is sensitivity to initial conditions. In SCT, the interaction between donor T cells and the recipient may be considered as such a system in which, graft source, conditioning and early immunosuppression profoundly influence immune reconstitution over time. This eventually determines clinical outcomes, either the emergence of tolerance or the development of graft versus host disease. In this paper parallels between SCT and dynamical systems are explored and a conceptual framework for developing mathematical models to understand disparate transplant outcomes is proposed.Comment: 23 pages, 4 figures. Updated version with additional data, 2 new figures and editorial revisions. New authors adde

Assaying the regulatory potential of mammalian conserved non-coding sequences in human cells

The fraction of experimentally active conserved non-coding sequences within any given cell type is low, so classical assays are unlikely to expose their potential

A spatially anchored transcriptomic atlas of the human kidney papilla identifies significant immune injury in patients with stone disease

Kidney stone disease causes significant morbidity and increases health care utilization. In this work, we decipher the cellular and molecular niche of the human renal papilla in patients with calcium oxalate (CaOx) stone disease and healthy subjects. In addition to identifying cell types important in papillary physiology, we characterize collecting duct cell subtypes and an undifferentiated epithelial cell type that was more prevalent in stone patients. Despite the focal nature of mineral deposition in nephrolithiasis, we uncover a global injury signature characterized by immune activation, oxidative stress and extracellular matrix remodeling. We also identify the association of MMP7 and MMP9 expression with stone disease and mineral deposition, respectively. MMP7 and MMP9 are significantly increased in the urine of patients with CaOx stone disease, and their levels correlate with disease activity. Our results define the spatial molecular landscape and specific pathways contributing to stone-mediated injury in the human papilla and identify associated urinary biomarkers

Dynamical System Modeling of Immune Reconstitution after Allogeneic Stem Cell Transplantation Identifies Patients at Risk for Adverse Outcomes

AbstractSystems that evolve over time and follow mathematical laws as they evolve are called dynamical systems. Lymphocyte recovery and clinical outcomes in 41 allograft recipients conditioned using antithymocyte globulin (ATG) and 4.5-Gy total body irradiation were studied to determine if immune reconstitution could be described as a dynamical system. Survival, relapse, and graft-versus-host disease (GVHD) were not significantly different in 2 cohorts of patients receiving different doses of ATG. However, donor-derived CD3+ cell reconstitution was superior in the lower ATG dose cohort, and there were fewer instances of donor lymphocyte infusion (DLI). Lymphoid recovery was plotted in each individual over time and demonstrated 1 of 3 sigmoid growth patterns: Pattern A (n = 15) had rapid growth with high lymphocyte counts, pattern B (n = 14) had slower growth with intermediate recovery, and pattern C (n = 10) had poor lymphocyte reconstitution. There was a significant association between lymphocyte recovery patterns and both the rate of change of donor-derived CD3+ at day 30 after stem cell transplantation (SCT) and clinical outcomes. GVHD was observed more frequently with pattern A, relapse and DLI more so with pattern C, with a consequent survival advantage in patients with patterns A and B. We conclude that evaluating immune reconstitution after SCT as a dynamical system may differentiate patients at risk of adverse outcomes and allow early intervention to modulate that risk

Genome-Scale Validation of Deep-Sequencing Libraries

Chromatin immunoprecipitation followed by high-throughput (HTP) sequencing (ChIP-seq) is a powerful tool to establish protein-DNA interactions genome-wide. The primary limitation of its broad application at present is the often-limited access to sequencers. Here we report a protocol, Mab-seq, that generates genome-scale quality evaluations for nucleic acid libraries intended for deep-sequencing. We show how commercially available genomic microarrays can be used to maximize the efficiency of library creation and quickly generate reliable preliminary data on a chromosomal scale in advance of deep sequencing. We also exploit this technique to compare enriched regions identified using microarrays with those identified by sequencing, demonstrating that they agree on a core set of clearly identified enriched regions, while characterizing the additional enriched regions identifiable using HTP sequencing

A Glial Signature and Wnt7 Signaling Regulate Glioma-Vascular Interactions and Tumor Microenvironment.

Gliomas comprise heterogeneous malignant glial and stromal cells. While blood vessel co-option is a potential mechanism to escape anti-angiogenic therapy, the relevance of glial phenotype in this process is unclear. We show that Olig2+ oligodendrocyte precursor-like glioma cells invade by single-cell vessel co-option and preserve the blood-brain barrier (BBB). Conversely, Olig2-negative glioma cells form dense perivascular collections and promote angiogenesis and BBB breakdown, leading to innate immune cell activation. Experimentally, Olig2 promotes Wnt7b expression, a finding that correlates in human glioma profiling. Targeted Wnt7a/7b deletion or pharmacologic Wnt inhibition blocks Olig2+ glioma single-cell vessel co-option and enhances responses to temozolomide. Finally, Olig2 and Wnt7 become upregulated after anti-VEGF treatment in preclinical models and patients. Thus, glial-encoded pathways regulate distinct glioma-vascular microenvironmental interactions

Integration of spatial and single-cell transcriptomics localizes epithelial cell–immune cross-talk in kidney injury

Single-cell sequencing studies have characterized the transcriptomic signature of cell types within the kidney. However, the spatial distribution of acute kidney injury (AKI) is regional and affects cells heterogeneously. We first optimized coordination of spatial transcriptomics and single-nuclear sequencing data sets, mapping 30 dominant cell types to a human nephrectomy. The predicted cell-type spots corresponded with the underlying histopathology. To study the implications of AKI on transcript expression, we then characterized the spatial transcriptomic signature of 2 murine AKI models: ischemia/reperfusion injury (IRI) and cecal ligation puncture (CLP). Localized regions of reduced overall expression were associated with injury pathways. Using single-cell sequencing, we deconvoluted the signature of each spatial transcriptomic spot, identifying patterns of colocalization between immune and epithelial cells. Neutrophils infiltrated the renal medulla in the ischemia model. Atf3 was identified as a chemotactic factor in S3 proximal tubules. In the CLP model, infiltrating macrophages dominated the outer cortical signature, and Mdk was identified as a corresponding chemotactic factor. The regional distribution of these immune cells was validated with multiplexed CO-Detection by indEXing (CODEX) immunofluorescence. Spatial transcriptomic sequencing complemented single-cell sequencing by uncovering mechanisms driving immune cell infiltration and detection of relevant cell subpopulations

An expansive human regulatory lexicon encoded in transcription factor footprints.

Regulatory factor binding to genomic DNA protects the underlying sequence from cleavage by DNase I, leaving nucleotide-resolution footprints. Using genomic DNase I footprinting across 41 diverse cell and tissue types, we detected 45 million transcription factor occupancy events within regulatory regions, representing differential binding to 8.4 million distinct short sequence elements. Here we show that this small genomic sequence compartment, roughly twice the size of the exome, encodes an expansive repertoire of conserved recognition sequences for DNA-binding proteins that nearly doubles the size of the human cis-regulatory lexicon. We find that genetic variants affecting allelic chromatin states are concentrated in footprints, and that these elements are preferentially sheltered from DNA methylation. High-resolution DNase I cleavage patterns mirror nucleotide-level evolutionary conservation and track the crystallographic topography of protein-DNA interfaces, indicating that transcription factor structure has been evolutionarily imprinted on the human genome sequence. We identify a stereotyped 50-base-pair footprint that precisely defines the site of transcript origination within thousands of human promoters. Finally, we describe a large collection of novel regulatory factor recognition motifs that are highly conserved in both sequence and function, and exhibit cell-selective occupancy patterns that closely parallel major regulators of development, differentiation and pluripotency