Molecular cartography uncovers evolutionary and microenvironmental dynamics in sporadic colorectal tumors.

Colorectal cancer exhibits dynamic cellular and genetic heterogeneity during progression from precursor lesions toward malignancy. Analysis of spatial multi-omic data from 31 human colorectal specimens enabled phylogeographic mapping of tumor evolution that revealed individualized progression trajectories and accompanying microenvironmental and clonal alterations. Phylogeographic mapping ordered genetic events, classified tumors by their evolutionary dynamics, and placed clonal regions along global pseudotemporal progression trajectories encompassing the chromosomal instability (CIN+) and hypermutated (HM) pathways. Integrated single-cell and spatial transcriptomic data revealed recurring epithelial programs and infiltrating immune states along progression pseudotime. We discovered an immune exclusion signature (IEX), consisting of extracellular matrix regulators DDR1, TGFBI, PAK4, and DPEP1, that charts with CIN+ tumor progression, is associated with reduced cytotoxic cell infiltration, and shows prognostic value in independent cohorts. This spatial multi-omic atlas provides insights into colorectal tumor-microenvironment co-evolution, serving as a resource for stratification and targeted treatments.

Differential pre-malignant programs and microenvironment chart distinct paths to malignancy in human colorectal polyps.

Colorectal cancers (CRCs) arise from precursor polyps whose cellular origins, molecular heterogeneity, and immunogenic potential may reveal diagnostic and therapeutic insights when analyzed at high resolution. We present a single-cell transcriptomic and imaging atlas of the two most common human colorectal polyps, conventional adenomas and serrated polyps, and their resulting CRC counterparts. Integrative analysis of 128 datasets from 62 participants reveals adenomas arise from WNT-driven expansion of stem cells, while serrated polyps derive from differentiated cells through gastric metaplasia. Metaplasia-associated damage is coupled to a cytotoxic immune microenvironment preceding hypermutation, driven partly by antigen-presentation differences associated with tumor cell-differentiation status. Microsatellite unstable CRCs contain distinct non-metaplastic regions where tumor cells acquire stem cell properties and cytotoxic immune cells are depleted. Our multi-omic atlas provides insights into malignant progression of colorectal polyps and their microenvironment, serving as a framework for precision surveillance and prevention of CRC.

Temporal recording of mammalian development and precancer.

Temporal ordering of cellular events offers fundamental insights into biological phenomena. Although this is traditionally achieved through continuous direct observations1,2, an alternative solution leverages irreversible genetic changes, such as naturally occurring mutations, to create indelible marks that enables retrospective temporal ordering3-5. Using a multipurpose, single-cell CRISPR platform, we developed a molecular clock approach to record the timing of cellular events and clonality in vivo, with incorporation of cell state and lineage information. Using this approach, we uncovered precise timing of tissue-specific cell expansion during mouse embryonic development, unconventional developmental relationships between cell types and new epithelial progenitor states by their unique genetic histories. Analysis of mouse adenomas, coupled to multiomic and single-cell profiling of human precancers, with clonal analysis of 418 human polyps, demonstrated the occurrence of polyclonal initiation in 15-30% of colonic precancers, showing their origins from multiple normal founders. Our study presents a multimodal framework that lays the foundation for in vivo recording, integrating synthetic or natural indelible genetic changes with single-cell analyses, to explore the origins and timing of development and tumorigenesis in mammalian systems.

Succinate Produced by Intestinal Microbes Promotes Specification of Tuft Cells to Suppress Ileal Inflammation.

BACKGROUND & AIMS: Countries endemic for parasitic infestations have a lower incidence of Crohn's disease (CD) than nonendemic countries, and there have been anecdotal reports of the beneficial effects of helminths in CD patients. Tuft cells in the small intestine sense and direct the immune response against eukaryotic parasites. We investigated the activities of tuft cells in patients with CD and mouse models of intestinal inflammation. METHODS: We used microscopy to quantify tuft cells in intestinal specimens from patients with ileal CD (n = 19), healthy individuals (n = 14), and TNFΔARE/+ mice, which develop Crohn's-like ileitis. We performed single-cell RNA sequencing, mass spectrometry, and microbiome profiling of intestinal tissues from wild-type and Atoh1-knockout mice, which have expansion of tuft cells, to study interactions between microbes and tuft cell populations. We assessed microbe dependence of tuft cell populations using microbiome depletion, organoids, and microbe transplant experiments. We used multiplex imaging and cytokine assays to assess alterations in inflammatory response following expansion of tuft cells with succinate administration in TNFΔARE/+ and anti-CD3E CD mouse models. RESULTS: Inflamed ileal tissues from patients and mice had reduced numbers of tuft cells, compared with healthy individuals or wild-type mice. Expansion of tuft cells was associated with increased expression of genes that regulate the tricarboxylic acid cycle, which resulted from microbe production of the metabolite succinate. Experiments in which we manipulated the intestinal microbiota of mice revealed the existence of an ATOH1-independent population of tuft cells that was sensitive to metabolites produced by microbes. Administration of succinate to mice expanded tuft cells and reduced intestinal inflammation in TNFΔARE/+ mice and anti-CD3E-treated mice, increased GATA3+ cells and type 2 cytokines (IL22, IL25, IL13), and decreased RORGT+ cells and type 17 cytokines (IL23) in a tuft cell-dependent manner. CONCLUSIONS: We found that tuft cell expansion reduced chronic intestinal inflammation in mice. Strategies to expand tuft cells might be developed for treatment of CD.

Dual indexed library design enables compatibility of in-Drop single-cell RNA-sequencing with exAMP chemistry sequencing platforms.

BACKGROUND: The increasing demand of single-cell RNA-sequencing (scRNA-seq) experiments, such as the number of experiments and cells queried per experiment, necessitates higher sequencing depth coupled to high data quality. New high-throughput sequencers, such as the Illumina NovaSeq 6000, enables this demand to be filled in a cost-effective manner. However, current scRNA-seq library designs present compatibility challenges with newer sequencing technologies, such as index-hopping, and their ability to generate high quality data has yet to be systematically evaluated. RESULTS: Here, we engineered a dual-indexed library structure, called TruDrop, on top of the inDrop scRNA-seq platform to solve these compatibility challenges, such that TruDrop libraries and standard Illumina libraries can be sequenced alongside each other on the NovaSeq. On scRNA-seq libraries, we implemented a previously-documented countermeasure to the well-described problem of index-hopping, demonstrated significant improvements in base-calling accuracy on the NovaSeq, and provided an example of multiplexing twenty-four scRNA-seq libraries simultaneously. We showed favorable comparisons in transcriptional diversity of TruDrop compared with prior inDrop libraries. CONCLUSIONS: Our approach enables cost-effective, high throughput generation of sequencing data with high quality, which should enable more routine use of scRNA-seq technologies.

Intestinal bile acid sequestration improves glucose control by stimulating hepatic miR-182-5p in type 2 diabetes.

Colesevelam is a bile acid sequestrant approved to treat both hyperlipidemia and type 2 diabetes, but the mechanism for its glucose-lowering effects is not fully understood. The aim of this study was to investigate the role of hepatic microRNAs (miRNAs) as regulators of metabolic disease and to investigate the link between the cholesterol and glucose-lowering effects of colesevelam. To quantify the impact of colesevelam treatment in rodent models of diabetes, metabolic studies were performed in Zucker diabetic fatty (ZDF) rats and db/db mice. Colesevelam treatments significantly decreased plasma glucose levels and increased glycolysis in the absence of changes to insulin levels in ZDF rats and db/db mice. High-throughput sequencing and real-time PCR were used to quantify hepatic miRNA and mRNA changes, and the cholesterol-sensitive miR-96/182/183 cluster was found to be significantly increased in livers from ZDF rats treated with colesevelam compared with vehicle controls. Inhibition of miR-182 in vivo attenuated colesevelam-mediated improvements to glycemic control in db/db mice. Hepatic expression of mediator complex subunit 1 (MED1), a nuclear receptor coactivator, was significantly decreased with colesevelam treatments in db/db mice, and MED1 was experimentally validated to be a direct target of miR-96/182/183 in humans and mice. In summary, these results support that colesevelam likely improves glycemic control through hepatic miR-182-5p, a mechanism that directly links cholesterol and glucose metabolism. NEW & NOTEWORTHY Colesevelam lowers systemic glucose levels in Zucker diabetic fatty rats and db/db mice and increases hepatic levels of the sterol response element binding protein 2-responsive microRNA cluster miR-96/182/183. Inhibition of miR-182 in vivo reverses the glucose-lowering effects of colesevelam in db/db mice. Mediator complex subunit 1 (MED1) is a novel, direct target of the miR-96/182/183 cluster in mice and humans.

Identification and multimodal characterization of a specialized epithelial cell type associated with Crohn's disease.

Crohn's disease (CD) is a complex chronic inflammatory disorder with both gastrointestinal and extra-intestinal manifestations associated immune dysregulation. Analyzing 202,359 cells from 170 specimens across 83 patients, we identify a distinct epithelial cell type in both terminal ileum and ascending colon (hereon as 'LND') with high expression of LCN2, NOS2, and DUOX2 and genes related to antimicrobial response and immunoregulation. LND cells, confirmed by in-situ RNA and protein imaging, are rare in non-IBD controls but expand in active CD, and actively interact with immune cells and specifically express IBD/CD susceptibility genes, suggesting a possible function in CD immunopathogenesis. Furthermore, we discover early and late LND subpopulations with different origins and developmental potential. A higher ratio of late-to-early LND cells correlates with better response to anti-TNF treatment. Our findings thus suggest a potential pathogenic role for LND cells in both Crohn's ileitis and colitis.

Temporal recording of mammalian development and precancer.

Key to understanding many biological phenomena is knowing the temporal ordering of cellular events, which often require continuous direct observations [1, 2]. An alternative solution involves the utilization of irreversible genetic changes, such as naturally occurring mutations, to create indelible markers that enables retrospective temporal ordering [3-8]. Using NSC-seq, a newly designed and validated multi-purpose single-cell CRISPR platform, we developed a molecular clock approach to record the timing of cellular events and clonality in vivo , while incorporating assigned cell state and lineage information. Using this approach, we uncovered precise timing of tissue-specific cell expansion during murine embryonic development and identified new intestinal epithelial progenitor states by their unique genetic histories. NSC-seq analysis of murine adenomas and single-cell multi-omic profiling of human precancers as part of the Human Tumor Atlas Network (HTAN), including 116 scRNA-seq datasets and clonal analysis of 418 human polyps, demonstrated the occurrence of polyancestral initiation in 15-30% of colonic precancers, revealing their origins from multiple normal founders. Thus, our multimodal framework augments existing single-cell analyses and lays the foundation for in vivo multimodal recording, enabling the tracking of lineage and temporal events during development and tumorigenesis.

A Specialized Epithelial Cell Type Regulating Mucosal Immunity and Driving Human Crohn's Disease.

Crohn's disease (CD) is a complex chronic inflammatory disorder that may affect any part of gastrointestinal tract with extra-intestinal manifestations and associated immune dysregulation. To characterize heterogeneity in CD, we profiled single-cell transcriptomics of 170 samples from 65 CD patients and 18 non-inflammatory bowel disease (IBD) controls in both the terminal ileum (TI) and ascending colon (AC). Analysis of 202,359 cells identified a novel epithelial cell type in both TI and AC, featuring high expression of LCN2, NOS2, and DUOX2, and thus is named LND. LND cells, confirmed by high-resolution in-situ RNA imaging, were rarely found in non-IBD controls, but expanded significantly in active CD. Compared to other epithelial cells, genes defining LND cells were enriched in antimicrobial response and immunoregulation. Moreover, multiplexed protein imaging demonstrated that LND cell abundance was associated with immune infiltration. Cross-talk between LND and immune cells was explored by ligand-receptor interactions and further evidenced by their spatial colocalization. LND cells showed significant enrichment of expression specificity of IBD/CD susceptibility genes, revealing its role in immunopathogenesis of CD. Investigating lineage relationships of epithelial cells detected two LND cell subpopulations with different origins and developmental potential, early and late LND. The ratio of the late to early LND cells was related to anti-TNF response. These findings emphasize the pathogenic role of the specialized LND cell type in both Crohn's ileitis and Crohn's colitis and identify novel biomarkers associated with disease activity and treatment response.

Processing single-cell RNA-seq data for dimension reduction-based analyses using open-source tools.

Single-cell RNA sequencing data require several processing procedures to arrive at interpretable results. While commercial platforms can serve as "one-stop shops" for data analysis, they relinquish the flexibility required for customized analyses and are often inflexible between experimental systems. For instance, there is no universal solution for the discrimination of informative or uninformative encapsulated cellular material; thus, pipeline flexibility takes priority. Here, we demonstrate a full data analysis pipeline, constructed modularly from open-source software, including tools that we have contributed. For complete details on the use and execution of this protocol, please refer to Petukhov et al. (2018), Heiser et al. (2020), and Heiser and Lau (2020).

Profile of Podocyte Translatome During Development of Type 2 and Type 1 Diabetic Nephropathy Using Podocyte-Specific TRAP mRNA RNA-seq.

Podocyte injury is important in development of diabetic nephropathy (DN). Although several studies have reported single-cell-based RNA sequencing (RNA-seq) of podocytes in type 1 DN (T1DN), the podocyte translating mRNA profile in type 2 DN (T2DN) has not previously been compared with that of T1DN. We analyzed the podocyte translatome in T2DN in podocin-Cre; Rosa26fsTRAP; eNOS-/-; db/db mice and compared it with that of streptozotocin-induced T1DN in podocin-Cre; Rosa26fsTRAP; eNOS-/- mice using translating ribosome affinity purification (TRAP) and RNA-seq. More than 125 genes were highly enriched in the podocyte ribosome. More podocyte TRAP genes were differentially expressed in T2DN than in T1DN. TGF-ß signaling pathway genes were upregulated, while MAPK pathway genes were downregulated only in T2DN, while ATP binding and cAMP-mediated signaling genes were downregulated only in T1DN. Genes regulating actin filament organization and apoptosis increased, while genes regulating VEGFR signaling and glomerular basement membrane components decreased in both type 1 and type 2 diabetic podocytes. A number of diabetes-induced genes not previously linked to podocyte injury were confirmed in both mouse and human DN. On the basis of differences and similarities in the podocyte translatome in T2DN and T1DN, investigators can identify factors underlying the pathophysiology of DN and novel therapeutic targets to treat diabetes-induced podocyte injury.

Dual inhibition of endothelial miR-92a-3p and miR-489-3p reduces renal injury-associated atherosclerosis.

BACKGROUND AND AIMS: Cardiovascular disease (CVD) is the leading cause of death in chronic kidney disease (CKD) patients, however, the underlying mechanisms that link CKD and CVD are not fully understood and limited treatment options exist in this high-risk population. microRNAs (miRNA) are critical regulators of gene expression for many biological processes in atherosclerosis, including endothelial dysfunction and inflammation. We hypothesized that renal injury-induced endothelial miRNAs promote atherosclerosis. Here, we demonstrate that dual inhibition of endothelial miRNAs inhibits atherosclerosis in the setting of renal injury. METHODS: Aortic endothelial miRNAs were analyzed in apolipoprotein E-deficient (Apoe-/-) mice with renal damage (5/6 nephrectomy, 5/6Nx) by real-time PCR. Endothelial miR-92a-3p and miR-489-3p were inhibited by locked-nucleic acid (LNA) miRNA inhibitors complexed to HDL. RESULTS: Renal injury significantly increased endothelial miR-92a-3p levels in Apoe-/-;5/6Nx mice. Dual inhibition of miR-92a-3p and miR-489-3p in Apoe-/-;5/6Nx with a single injection of HDL + LNA inhibitors significantly reduced atherosclerotic lesion area by 28.6% compared to HDL + LNA scramble (LNA-Scr) controls. To examine the impact of dual LNA treatment on aortic endothelial gene expression, total RNA sequencing was completed, and multiple putative target genes and pathways were identified to be significantly altered, including the STAT3 immune response pathway. Among the differentially expressed genes, Tgfb2 and Fam220a were identified as putative targets of miR-489-3p and miR-92a-3p, respectively. Both Tgfb2 and Fam220a were significantly increased in aortic endothelium after miRNA inhibition in vivo compared to HDL + LNA-Scr controls. Furthermore, Tgfb2 and Fam220a were validated with gene reporter assays as direct targets of miR-489-3p and miR-92a-3p, respectively. In human coronary artery endothelial cells, over-expression and inhibition of miR-92a-3p decreased and increased FAM220A expression, respectively. Moreover, miR-92a-3p overexpression increased STAT3 phosphorylation, likely through direct regulation of FAM220A, a negative regulator of STAT3 phosphorylation. CONCLUSIONS: These results support endothelial miRNAs as therapeutic targets and dual miRNA inhibition as viable strategy to reduce CKD-associated atherosclerosis.

Beta cell secretion of miR-375 to HDL is inversely associated with insulin secretion.

Extracellular microRNAs (miRNAs) are a new class of biomarkers for cellular phenotypes and disease, and are bioactive signals within intercellular communication networks. Previously, we reported that miRNAs are secreted from macrophage to high-density lipoproteins (HDL) and delivered to recipient cells to regulate gene expression. Despite the potential importance of HDL-miRNAs, regulation of HDL-miRNA export from cells has not been fully studied. Here, we report that pancreatic islets and beta cells abundantly export miR-375-3p to HDL and this process is inhibited by cellular mechanisms that promote insulin secretion. Small RNA sequencing and PCR approaches were used to quantify beta cell miRNA export to HDL. Strikingly, high glucose conditions were found to inhibit HDL-miR-375-3p export, which was dependent on extracellular calcium. Likewise, stimulation of cAMP was found to repress HDL-miR-375-3p export. Furthermore, we found that beta cell ATP-sensitive potassium channel (KATP) channels are required for HDL-miRNA export as chemical inhibition (tolbutamide) and global genetic knockout (Abcc8-/-) approaches inhibited HDL-miR-375-3p export. This process is not likely associated with cholesterol flux, as gain-of-function and loss-of-function studies for cholesterol transporters failed to alter HDL-miR-375-3p export. In conclusion, results support that pancreatic beta cells export miR-375-3p to HDL and this process is inversely regulated to insulin secretion.

Bioinformatic analysis of endogenous and exogenous small RNAs on lipoproteins.

To comprehensively study extracellular small RNAs (sRNA) by sequencing (sRNA-seq), we developed a novel pipeline to overcome current limitations in analysis entitled, "Tools for Integrative Genome analysis of Extracellular sRNAs (TIGER)". To demonstrate the power of this tool, sRNA-seq was performed on mouse lipoproteins, bile, urine and livers. A key advance for the TIGER pipeline is the ability to analyse both host and non-host sRNAs at genomic, parent RNA and individual fragment levels. TIGER was able to identify approximately 60% of sRNAs on lipoproteins and >85% of sRNAs in liver, bile and urine, a significant advance compared to existing software. Moreover, TIGER facilitated the comparison of lipoprotein sRNA signatures to disparate sample types at each level using hierarchical clustering, correlations, beta-dispersions, principal coordinate analysis and permutational multivariate analysis of variance. TIGER analysis was also used to quantify distinct features of exRNAs, including 5' miRNA variants, 3' miRNA non-templated additions and parent RNA positional coverage. Results suggest that the majority of sRNAs on lipoproteins are non-host sRNAs derived from bacterial sources in the microbiome and environment, specifically rRNA-derived sRNAs from Proteobacteria. Collectively, TIGER facilitated novel discoveries of lipoprotein and biofluid sRNAs and has tremendous applicability for the field of extracellular RNA.