A single-cell atlas characterizes dysregulation of the bone marrow immune microenvironment associated with outcomes in multiple myeloma.

Multiple Myeloma (MM) remains incurable despite advances in treatment options. Although tumor subtypes and specific DNA abnormalities are linked to worse prognosis, the impact of immune dysfunction on disease emergence and/or treatment sensitivity remains unclear. We established a harmonized consortium to generate an Immune Atlas of MM aimed at informing disease etiology, risk stratification, and potential therapeutic strategies. We generated a transcriptome profile of 1,149,344 single cells from the bone marrow of 263 newly diagnosed patients enrolled in the CoMMpass study and characterized immune and hematopoietic cell populations. Associating cell abundances and gene expression with disease progression revealed the presence of a proinflammatory immune senescence-associated secretory phenotype in rapidly progressing patients. Furthermore, signaling analyses suggested active intercellular communication involving APRIL-BCMA, potentially promoting tumor growth and survival. Finally, we demonstrate that integrating immune cell levels with genetic information can significantly improve patient stratification.

Noncanonical TRAIL Signaling Promotes Myeloid-Derived Suppressor Cell Abundance and Tumor Growth in Cholangiocarcinoma.

BACKGROUND & AIMS: Proapoptotic tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) signaling as a cause of cancer cell death is a well-established mechanism. However, TRAIL-receptor (TRAIL-R) agonists have had very limited anticancer activity in human beings, challenging the concept of TRAIL as a potent anticancer agent. Herein, we aimed to define mechanisms by which TRAIL+ cancer cells can leverage noncanonical TRAIL signaling in myeloid-derived suppressor cells (MDSCs) promoting their abundance in murine cholangiocarcinoma (CCA). METHODS: Multiple immunocompetent syngeneic, orthotopic models of CCA were used. Single-cell RNA sequencing and cellular indexing of transcriptomes and epitopes by sequencing of CD45+ cells in murine tumors from the different CCA models was conducted. RESULTS: In multiple immunocompetent murine models of CCA, implantation of TRAIL+ murine cancer cells into Trail-r-/- mice resulted in a significant reduction in tumor volumes compared with wild-type mice. Tumor-bearing Trail-r-/- mice had a significant decrease in the abundance of MDSCs owing to attenuation of MDSC proliferation. Noncanonical TRAIL signaling with consequent nuclear factor-κB activation in MDSCs facilitated enhanced MDSC proliferation. Single-cell RNA sequencing and cellular indexing of transcriptomes and epitopes by sequencing of immune cells from murine tumors showed enrichment of a nuclear factor-κB activation signature in MDSCs. Moreover, MDSCs were resistant to TRAIL-mediated apoptosis owing to enhanced expression of cellular FLICE inhibitory protein, an inhibitor of proapoptotic TRAIL signaling. Accordingly, cellular FLICE inhibitory protein knockdown sensitized murine MDSCs to TRAIL-mediated apoptosis. Finally, cancer cell-restricted deletion of Trail significantly reduced MDSC abundance and murine tumor burden. CONCLUSIONS: Our findings highlight the therapeutic potential of targeting TRAIL+ cancer cells for treatment of a poorly immunogenic cancer.

Characterization and Optimization of Multiomic Single-Cell Epigenomic Profiling.

The snATAC + snRNA platform allows epigenomic profiling of open chromatin and gene expression with single-cell resolution. The most critical assay step is to isolate high-quality nuclei to proceed with droplet-base single nuclei isolation and barcoding. With the increasing popularity of multiomic profiling in various fields, there is a need for optimized and reliable nuclei isolation methods, mainly for human tissue samples. Herein we compared different nuclei isolation methods for cell suspensions, such as peripheral blood mononuclear cells (PBMC, n = 18) and a solid tumor type, ovarian cancer (OC, n = 18), derived from debulking surgery. Nuclei morphology and sequencing output parameters were used to evaluate the quality of preparation. Our results show that NP-40 detergent-based nuclei isolation yields better sequencing results than collagenase tissue dissociation for OC, significantly impacting cell type identification and analysis. Given the utility of applying such techniques to frozen samples, we also tested frozen preparation and digestion (n = 6). A paired comparison between frozen and fresh samples validated the quality of both specimens. Finally, we demonstrate the reproducibility of scRNA and snATAC + snRNA platform, by comparing the gene expression profiling of PBMC. Our results highlight how the choice of nuclei isolation methods is critical for obtaining quality data in multiomic assays. It also shows that the measurement of expression between scRNA and snRNA is comparable and effective for cell type identification.

Noncanonical TRAIL Signaling Promotes Myeloid-Derived Suppressor Cell Abundance and Tumor Progression in Cholangiocarcinoma.

Proapoptotic tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) signaling as a cause of cancer cell death is a well-established mechanism. However, TRAIL-receptor (TRAIL-R) agonists have had very limited anticancer activity in humans, challenging the concept of TRAIL as a potent anticancer agent. Herein, we demonstrate that TRAIL + cancer cells can leverage noncanonical TRAIL signaling in myeloid-derived suppressor cells (MDSCs) promoting their abundance in murine cholangiocarcinoma (CCA). In multiple immunocompetent syngeneic, orthotopic murine models of CCA, implantation of TRAIL + murine cancer cells into Trail-r -/- mice resulted in a significant reduction in tumor volumes compared to wild type mice. Tumor bearing Trail-r -/- mice had a significant decrease in the abundance of MDSCs due to attenuation of MDSC proliferation. Noncanonical TRAIL signaling with consequent NF-κB activation in MDSCs facilitated enhanced MDSC proliferation. Single cell RNA sequencing and cellular indexing of transcriptomes and epitopes by sequencing (CITE-Seq) of CD45 + cells in murine tumors from three distinct immunocompetent CCA models demonstrated a significant enrichment of an NF-κB activation signature in MDSCs. Moreover, MDSCs were resistant to TRAIL-mediated apoptosis due to enhanced expression of cellular FLICE inhibitory protein (cFLIP), an inhibitor of proapoptotic TRAIL signaling. Accordingly, cFLIP knockdown sensitized murine MDSCs to TRAIL-mediated apoptosis. Finally, cancer cell-restricted deletion of Trail significantly reduced MDSC abundance and murine tumor burden. In summary, our findings define a noncanonical TRAIL signal in MDSCs and highlight the therapeutic potential of targeting TRAIL + cancer cells for the treatment of a poorly immunogenic cancer.

A Protocol for the Cryopreservation of Human Intestinal Mucosal Biopsies Compatible With Single-Cell Transcriptomics and Ex Vivo Studies.

The heterogeneity of the human intestinal epithelium has hindered the understanding of the pathophysiology of distinct specialized cell types on a single-cell basis in disease states. Described here is a workflow for the cryopreservation of endoscopically obtained human intestinal mucosal biopsies, subsequent preparation of this tissue to yield highly viable fluorescence-activated cell sorting (FACS)isolated human intestinal epithelial cell (IEC) single-cell suspensions compatible with successful library preparation and deep single-cell RNA sequencing (scRNAseq). We validated this protocol in deep scRNAseq of 59,653 intestinal cells in 10 human participants. Furthermore, primary intestinal cultures were successfully generated from cryopreserved tissue, capable of surviving in short-term culture and suitable for physiological assays studying gut peptide secretion from rare hormone-producing enteroendocrine cells in humans. This study offers an accessible avenue for single-cell transcriptomics and ex vivo studies from cryopreserved intestinal mucosal biopsies. These techniques may be used in the future to dissect and define novel aberrations to the intestinal ecosystem that lead to the development and progression of disease states in humans, even in rare IEC populations.

Pediatric brain tumor cell lines exhibit miRNA-depleted, Y RNA-enriched extracellular vesicles.

BACKGROUND: Medulloblastoma (MB) and diffuse infiltrative pontine glioma (DIPG) are malignant pediatric tumors. Extracellular vesicles (EVs) and their bioactive cargoes have been implicated in tumorigenesis. Most studies have focused on adult tumors, therefore the role of EVs and the noncoding RNA (ncRNA) landscape in pediatric brain tumors is not fully characterized. The overall aim of this pilot study was to isolate EVs from MB and DIPG patient-derived cell lines and to explore the small ncRNA transcriptome. METHODS: EVs from 3 DIPG and 4 MB patient-derived cell lines were analyzed. High-throughput next generation sequencing interrogated the short non-coding RNA (ncRNA) transcriptome. Known and novel miRNAs were quantified. Differential expression analysis, in silico target prediction, and functional gene enrichment were performed. RESULTS: EV secretomes from MB and DIPG patient-derived cell lines demonstrated discrete ncRNA biotypes. Notably, miRNAs were depleted and Y RNAs were enriched in EV samples. Hierarchical cluster analysis revealed high discrimination in miRNA expression between DIPG and MB cell lines and RNA-Seq identified novel miRNAs not previously implicated in MB or DIPG pathogenesis. Known and putative target genes of dysregulated miRNAs were identified. Functional annotation analysis of the target genes for differentially expressed EV-and parental-derived miRNAs revealed significant cancer-related pathway involvement. CONCLUSIONS: This hypothesis-generating study demonstrated that pediatric brain tumor-derived cell lines secrete EVs comprised of various ncRNA cargoes. Validation of these findings in patient samples may provide new insights into the pediatric brain tumor microenvironment and identification of novel therapeutic candidates.

Common Oncogene Mutations and Novel SND1-BRAF Transcript Fusion in Lung Adenocarcinoma from Never Smokers.

Lung adenocarcinomas from never smokers account for approximately 15 to 20% of all lung cancers and these tumors often carry genetic alterations that are responsive to targeted therapy. Here we examined mutation status in 10 oncogenes among 89 lung adenocarcinomas from never smokers. We also screened for oncogene fusion transcripts in 20 of the 89 tumors by RNA-Seq. In total, 62 tumors had mutations in at least one of the 10 oncogenes, including EGFR (49 cases, 55%), K-ras (5 cases, 6%), BRAF (4 cases, 5%), PIK3CA (3 cases, 3%), and ERBB2 (4 cases, 5%). In addition to ALK fusions identified by IHC/FISH in four cases, two previously known fusions involving EZR- ROS1 and KIF5B-RET were identified by RNA-Seq as well as a third novel fusion transcript that was formed between exons 1-9 of SND1 and exons 2 to 3' end of BRAF. This in-frame fusion was observed in 3/89 tested tumors and 2/64 additional never smoker lung adenocarcinoma samples. Ectopic expression of SND1-BRAF in H1299 cells increased phosphorylation levels of MEK/ERK, cell proliferation, and spheroid formation compared to parental mock-transfected control. Jointly, our results suggest a potential role of the novel BRAF fusion in lung cancer development and therapy.

Multi-platform analysis of microRNA expression measurements in RNA from fresh frozen and FFPE tissues.

MicroRNAs play a role in regulating diverse biological processes and have considerable utility as molecular markers for diagnosis and monitoring of human disease. Several technologies are available commercially for measuring microRNA expression. However, cross-platform comparisons do not necessarily correlate well, making it difficult to determine which platform most closely represents the true microRNA expression level in a tissue. To address this issue, we have analyzed RNA derived from cell lines, as well as fresh frozen and formalin-fixed paraffin embedded tissues, using Affymetrix, Agilent, and Illumina microRNA arrays, NanoString counting, and Illumina Next Generation Sequencing. We compared the performance within- and between the different platforms, and then verified these results with those of quantitative PCR data. Our results demonstrate that the within-platform reproducibility for each method is consistently high and although the gene expression profiles from each platform show unique traits, comparison of genes that were commonly detectable showed that detection of microRNA transcripts was similar across multiple platforms.

Quantitative miRNA expression analysis using fluidigm microfluidics dynamic arrays.

BACKGROUND: MicroRNAs (miRNAs) represent a growing class of small non-coding RNAs that are important regulators of gene expression in both plants and animals. Studies have shown that miRNAs play a critical role in human cancer and they can influence the level of cell proliferation and apoptosis by modulating gene expression. Currently, methods for the detection and measurement of miRNA expression include small and moderate-throughput technologies, such as standard quantitative PCR and microarray based analysis. However, these methods have several limitations when used in large clinical studies where a high-throughput and highly quantitative technology needed for the efficient characterization of a large number of miRNA transcripts in clinical samples. Furthermore, archival formalin fixed, paraffin embedded (FFPE) samples are increasingly becoming the primary resource for gene expression studies because fresh frozen (FF) samples are often difficult to obtain and requires special storage conditions. In this study, we evaluated the miRNA expression levels in FFPE and FF samples as well as several lung cancer cell lines employing a high throughput qPCR-based microfluidic technology. The results were compared to standard qPCR and hybridization-based microarray platforms using the same samples. RESULTS: We demonstrated highly correlated Ct values between multiplex and singleplex RT reactions in standard qPCR assays for miRNA expression using total RNA from A549 (R = 0.98; p < 0.0001) and H1299 (R = 0.95; p < 0.0001) lung cancer cell lines. The Ct values generated by the microfluidic technology (Fluidigm 48.48 dynamic array systems) resulted in a left-shift toward lower Ct values compared to those observed by ABI 7900 HT (mean difference, 3.79), suggesting that the microfluidic technology exhibited a greater sensitivity. In addition, we show that as little as 10 ng total RNA can be used to reliably detect all 48 or 96 tested miRNAs using a 96-multiplexing RT reaction in both FFPE and FF samples. Finally, we compared miRNA expression measurements in both FFPE and FF samples by qPCR using the 96.96 dynamic array and Affymetrix microarrays. Fold change comparisons for comparable genes between the two platforms indicated that the overall correlation was R = 0.60. The maximum fold change detected by the Affymetrix microarray was 3.5 compared to 13 by the 96.96 dynamic array. CONCLUSION: The qPCR-array based microfluidic dynamic array platform can be used in conjunction with multiplexed RT reactions for miRNA gene expression profiling. We showed that this approach is highly reproducible and the results correlate closely with the existing singleplex qPCR platform at a throughput that is 5 to 20 times higher and a sample and reagent usage that was approximately 50-100 times lower than conventional assays. We established optimal conditions for using the Fluidigm microfluidic technology for rapid, cost effective, and customizable arrays for miRNA expression profiling and validation.