Improved characterization of single-cell RNA-seq libraries with paired-end avidity sequencing.

Prevailing poly(dT)-primed 3' single-cell RNA-seq protocols generate barcoded cDNA fragments containing the reverse transcriptase priming site, which is expected to be the poly(A) tail or a genomic adenine homopolymer. Direct sequencing across this priming site was historically difficult because of DNA sequencing errors induced by the homopolymeric primer at the 'barcode' end. Here, we evaluate the capability of "avidity base chemistry" DNA sequencing from Element Biosciences to sequence through this homopolymer accurately, and the impact of the additional cDNA sequence on read alignment and precise quantification of polyadenylation site usage. We find that the Element Aviti instrument sequences through the thymine homopolymer into the subsequent cDNA sequence without detectable loss of accuracy. The resulting paired-end alignments enable direct and independent assignment of reads to polyadenylation sites, which bypasses complexities and limitations of conventional approaches but does not consistently improve read mapping rates compared to single-end alignment. We also characterize low-level artifacts and arrive at an adjusted adapter trimming and alignment workflow that significantly improves the alignment of sequence data from Element and Illumina, particularly in the context of extended read lengths. Our analyses confirm that Element avidity sequencing is an effective alternative to Illumina sequencing for standard single-cell RNA-seq, particularly for polyadenylation site analyses but do not rule out the potential for similar performance from other emerging platforms.

Ag-driven CD8+ T cell clonal expansion is a prominent feature of MASH in humans and mice.

BACKGROUND AND AIMS: Chronic liver disease due to metabolic dysfunction-associated steatohepatitis (MASH) is a rapidly increasing global epidemic. MASH progression is a consequence of the complex interplay between inflammatory insults and dysregulated hepatic immune responses. T lymphocytes have been shown to accumulate in the liver during MASH, but the cause and consequence of T cell accumulation in the liver remain unclear. Our study aimed to define the phenotype and T cell receptor diversity of T cells from human cirrhotic livers and an animal model of MASH to begin resolving their function in disease. APPROACH AND RESULTS: In these studies, we evaluated differences in T cell phenotype in the context of liver disease. Accordingly, we isolated liver resident T cell populations from humans with cirrhosis and from mice with diet-induced MASH. Using both 5' single-cell sequencing and flow cytometry, we defined the phenotype and T cell receptor repertoire of liver resident T cells during health and disease. CONCLUSIONS: MASH-induced human cirrhosis and diet-induced MASH in mice resulted in the accumulation of activated and clonally expanded T cells in the liver. The clonally expanded T cells in the liver expressed markers of chronic antigenic stimulation, including PD1, TIGIT, and TOX. Overall, this study establishes for the first time that T cells undergo Ag-dependent clonal expansion and functional differentiation during the progression of MASH. These studies could lead to the identification of antigenic targets that drive T cell activation, clonal expansion, and recruitment to the liver during MASH.

Targeting Acute Myeloid Leukemia Stem Cells through Perturbation of Mitochondrial Calcium.

Acute myeloid leukemia stem cells (LSCs) are uniquely reliant on oxidative phosphorylation (OXPHOS) for survival. Moreover, maintenance of OXPHOS is dependent on BCL-2, creating a therapeutic opportunity to target LSCs using the BCL-2 inhibitor venetoclax. Although venetoclax-based regimens have shown promising clinical activity, the emergence of drug resistance is prevalent. Thus, in the present study, we investigated how mitochondrial properties may influence venetoclax responsiveness. Our data show that utilization of mitochondrial calcium is fundamentally different between drug-responsive and nonresponsive LSCs. By comparison, venetoclax-resistant LSCs demonstrate an active metabolic (i.e., OXPHOS) status with relatively high levels of calcium. Consequently, we tested genetic and pharmacological approaches to target the mitochondrial calcium uniporter. We demonstrate that inhibition of calcium uptake reduces OXPHOS and leads to eradication of venetoclax-resistant LSCs. These findings demonstrate a central role for calcium signaling in LSCs and provide an avenue for clinical management of venetoclax resistance. Significance: We identify increased utilization of mitochondrial calcium as a distinct metabolic requirement of venetoclax-resistant LSCs and demonstrate the potential of targeting mitochondrial calcium uptake as a therapeutic strategy.

Antigen-driven CD8 + T cell clonal expansion is a prominent feature of MASH in humans and mice.

Background and Aims: Chronic liver disease due to metabolic dysfunction-associated steatohepatitis (MASH) is a rapidly increasing global epidemic. MASH progression is a consequence of the complex interplay between inflammatory insults and dysregulated hepatic immune responses. T lymphocytes have been shown to accumulate in the liver during MASH, but the cause and consequence of T cell accumulation in the liver remain unclear. Our study aimed to define the phenotype and T cell receptor diversity of T cells from human cirrhotic livers and an animal model of MASH to begin resolving their function in disease. Approach and Results: In these studies, we evaluated differences in T cell phenotype in the context of liver disease we isolated liver resident T cell populations from individuals with cirrhosis and a murine model of MASH. Using both 5' single cell sequencing and flow cytometry we defined the phenotype and T cell receptor repertoire of liver resident T cells during health and disease. Conclusions: MASH-induced cirrhosis and diet-induced MASH in mice resulted in the accumulation of activated and clonally expanded T cells in the liver. The clonally expanded T cells in the liver expressed markers of chronic antigenic stimulation, including PD1 , TIGIT and TOX . Overall, this study establishes for the first time that T cells undergo antigen-dependent clonal expansion and functional differentiation during the progression of MASH. These studies could lead to the identification of potential antigenic targets that drive T cell activation, clonal expansion, and recruitment to the liver during MASH.

An immunophenotype-coupled transcriptomic atlas of human hematopoietic progenitors.

Analysis of the human hematopoietic progenitor compartment is being transformed by single-cell multimodal approaches. Cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) enables coupled surface protein and transcriptome profiling, thereby revealing genomic programs underlying progenitor states. To perform CITE-seq systematically on primary human bone marrow cells, we used titrations with 266 CITE-seq antibodies (antibody-derived tags) and machine learning to optimize a panel of 132 antibodies. Multimodal analysis resolved >80 stem, progenitor, immune, stromal and transitional cells defined by distinctive surface markers and transcriptomes. This dataset enables flow cytometry solutions for in silico-predicted cell states and identifies dozens of cell surface markers consistently detected across donors spanning race and sex. Finally, aligning annotations from this atlas, we nominate normal marrow equivalents for acute myeloid leukemia stem cell populations that differ in clinical response. This atlas serves as an advanced digital resource for hematopoietic progenitor analyses in human health and disease.

Compromised nonsense-mediated RNA decay results in truncated RNA-binding protein production upon DUX4 expression.

Nonsense-mediated RNA decay (NMD) degrades transcripts carrying premature termination codons. NMD is thought to prevent the synthesis of toxic truncated proteins. However, whether loss of NMD results in widespread production of truncated proteins is unclear. A human genetic disease, facioscapulohumeral muscular dystrophy (FSHD), features acute inhibition of NMD upon expression of the disease-causing transcription factor, DUX4. Using a cell-based model of FSHD, we show production of truncated proteins from physiological NMD targets and find that RNA-binding proteins are enriched for aberrant truncations. The NMD isoform of one RNA-binding protein, SRSF3, is translated to produce a stable truncated protein, which is detected in FSHD patient-derived myotubes. Ectopic expression of truncated SRSF3 confers toxicity, and its downregulation is cytoprotective. Our results delineate the genome-scale impact of NMD loss. This widespread production of potentially deleterious truncated proteins has implications for FSHD biology as well as other genetic diseases where NMD is therapeutically modulated.

A Novel Type of Monocytic Leukemia Stem Cell Revealed by the Clinical Use of Venetoclax-Based Therapy.

The BCL2 inhibitor venetoclax has recently emerged as an important component of acute myeloid leukemia (AML) therapy. Notably, use of this agent has revealed a previously unrecognized form of pathogenesis characterized by monocytic disease progression. We demonstrate that this form of disease arises from a fundamentally different type of leukemia stem cell (LSC), which we designate as monocytic LSC (m-LSC), that is developmentally and clinically distinct from the more well-described primitive LSC (p-LSC). The m-LSC is distinguished by a unique immunophenotype (CD34-, CD4+, CD11b-, CD14-, CD36-), unique transcriptional state, reliance on purine metabolism, and selective sensitivity to cladribine. Critically, in some instances, m-LSC and p-LSC subtypes can co-reside in the same patient with AML and simultaneously contribute to overall tumor biology. Thus, our findings demonstrate that LSC heterogeneity has direct clinical significance and highlight the need to distinguish and target m-LSCs as a means to improve clinical outcomes with venetoclax-based regimens. SIGNIFICANCE: These studies identify and characterize a new type of human acute myeloid LSC that is responsible for monocytic disease progression in patients with AML treated with venetoclax-based regimens. Our studies describe the phenotype, molecular properties, and drug sensitivities of this unique LSC subclass. This article is featured in Selected Articles from This Issue, p. 1949.

Adrenal gene expression dynamics support hibernation in 13-lined ground squirrels.

Hibernation is a natural model of extreme physiology in a mammal. Throughout winter, small hibernators repeatedly undergo rapid, dramatic swings in body temperature, perfusion, and oxygen delivery. To gain insight into the molecular mechanisms that support homeostasis despite the numerous challenges posed by this dynamic physiology, we collected 13-lined ground squirrel adrenal glands from at least five individuals representing six key timepoints across the year using body temperature telemetry. Differentially expressed genes were identified using RNA-seq, revealing both strong seasonal and torpor-arousal cycle effects on gene expression. Two novel findings emerge from this study. First, transcripts encoding multiple genes involved in steroidogenesis decreased seasonally. Taken together with morphometric analyses, the data are consistent with preservation of mineralocorticoids but suppression of glucocorticoid and androgen output throughout winter hibernation. Second, a temporally orchestrated, serial gene expression program unfolds across the brief arousal periods. This program initiates during early rewarming with the transient activation of a set of immediate early response (IER) genes, comprised of both transcription factors and the RNA degradation proteins that assure their rapid turnover. This pulse in turn activates a cellular stress response program to restore proteostasis comprised of protein turnover, synthesis, and folding machinery. These and other data support a general model for gene expression across the torpor-arousal cycle that is facilitated in synchrony with whole body temperature shifts; induction of the immediate early response upon rewarming activates a proteostasis program followed by a restored tissue-specific gene expression profile enabling renewal, repair, and survival of the torpid state.NEW & NOTEWORTHY This pioneer study of adrenal gland gene expression dynamics in hibernating ground squirrels leverages the power of RNA-seq on multiple precisely timed samples to demonstrate: 1) steroidogenesis is seasonally reorganized to preserve aldosterone at the expense of glucocorticoids and androgens throughout winter hibernation; 2) a serial gene expression program unfolds during each short arousal whereby immediate early response genes induce the gene expression machinery that restores proteostasis and the cell-specific expression profile before torpor reentry.

Diverse cell-specific patterns of alternative polyadenylation in Drosophila.

Most genes in higher eukaryotes express isoforms with distinct 3' untranslated regions (3' UTRs), generated by alternative polyadenylation (APA). Since 3' UTRs are predominant locations of post-transcriptional regulation, APA can render such programs conditional, and can also alter protein sequences via alternative last exon (ALE) isoforms. We previously used 3'-sequencing from diverse Drosophila samples to define multiple tissue-specific APA landscapes. Here, we exploit comprehensive single nucleus RNA-sequencing data (Fly Cell Atlas) to elucidate cell-type expression of 3' UTRs across >250 adult Drosophila cell types. We reveal the cellular bases of multiple tissue-specific APA/ALE programs, such as 3' UTR lengthening in differentiated neurons and 3' UTR shortening in spermatocytes and spermatids. We trace dynamic 3' UTR patterns across cell lineages, including in the male germline, and discover new APA patterns in the intestinal stem cell lineage. Finally, we correlate expression of RNA binding proteins (RBPs), miRNAs and global levels of cleavage and polyadenylation (CPA) factors in several cell types that exhibit characteristic APA landscapes, yielding candidate regulators of transcriptome complexity. These analyses provide a comprehensive foundation for future investigations of mechanisms and biological impacts of alternative 3' isoforms across the major cell types of this widely-studied model organism.

Single-cell RNA sequencing and binary hierarchical clustering define lung interstitial macrophage heterogeneity in response to hypoxia.

Few studies have examined lung interstitial macrophage (IM) molecular phenotypes after being exposed to hypoxia in vivo at the single-cell level, even though macrophages contribute to hypoxic pulmonary hypertension (PH). We aimed to determine IM diversity and its association with hypoxia-induced PH. We hypothesized that integrating single-cell RNA sequencing (scRNAseq) and binary hierarchal clustering (BHC) could resolve IM heterogeneity under normal homeostatic conditions and changes induced by hypoxia exposure. Cx3cr1GFP/+ reporter mice were exposed to normoxic conditions (∼21% [Formula: see text]) or exposed to 1 day (D1) or 7 days (D7) of hypoxia (∼10% [Formula: see text]). We used flow cytometry to isolate Cx3cr1+ IMs and the 10X Genomics platform for scRNAseq, Cell Ranger, Seurat, ClusterMap, monocle, ingenuity pathway analysis, and Fisher's exact test (q value < 0.05) for functional investigations. n = 374 (normoxia), n = 2,526 (D1), and n = 1,211 (D7) IMs were included in the analyses. We identified three normoxia-related cell types, five hypoxia-associated cell types that emerged at D1, and three that appeared at D7. We describe the existence of a putative resident trained innate IM, which is present in normoxia, transiently depleted at D1, and recovered after 7 days of sustained hypoxia. We also define a rare putative pathogenic population associated with transcripts implicated in PH development that emerges at D7. In closing, we describe the successful integration of BHC with scRNAseq to determine IM heterogeneity and its association with PH. These results shed light on how resident-trained innate IMs become more heterogeneous but ultimately accustomed to hypoxia.

Breast Cancer Endocrine Therapy Promotes Weight Gain With Distinct Adipose Tissue Effects in Lean and Obese Female Mice.

Breast cancer survivors treated with tamoxifen and aromatase inhibitors report weight gain and have an elevated risk of type 2 diabetes, especially if they have obesity. These patient experiences are inconsistent with, preclinical studies using high doses of tamoxifen which reported acute weight loss. We investigated the impact of breast cancer endocrine therapies in a preclinical model of obesity and in a small group of breast adipose tissue samples from women taking tamoxifen to understand the clinical findings. Mature female mice were housed at thermoneutrality and fed either a low-fat/low-sucrose (LFLS) or a high-fat/high-sucrose (HFHS) diet. Consistent with the high expression of Esr1 observed in mesenchymal stem cells from adipose tissue, endocrine therapy was associated with adipose accumulation and more preadipocytes compared with estrogen-treated control mice but resulted in fewer adipocyte progenitors only in the context of HFHS. Analysis of subcutaneous adipose stromal cells revealed diet- and treatment-dependent effects of endocrine therapies on various cell types and genes, illustrating the complexity of adipose tissue estrogen receptor signaling. Breast cancer therapies supported adipocyte hypertrophy and associated with hepatic steatosis, hyperinsulinemia, and glucose intolerance, particularly in obese females. Current tamoxifen use associated with larger breast adipocyte diameter only in women with obesity. Our translational studies suggest that endocrine therapies may disrupt adipocyte progenitors and support adipocyte hypertrophy, potentially leading to ectopic lipid deposition that may be linked to a greater type 2 diabetes risk. Monitoring glucose tolerance and potential interventions that target insulin action should be considered for some women receiving life-saving endocrine therapies for breast cancer.

Liver Transcriptome Dynamics During Hibernation Are Shaped by a Shifting Balance Between Transcription and RNA Stability.

Hibernators dramatically lower metabolism to save energy while fasting for months. Prolonged fasting challenges metabolic homeostasis, yet small-bodied hibernators emerge each spring ready to resume all aspects of active life, including immediate reproduction. The liver is the body's metabolic hub, processing and detoxifying macromolecules to provide essential fuels to brain, muscle and other organs throughout the body. Here we quantify changes in liver gene expression across several distinct physiological states of hibernation in 13-lined ground squirrels, using RNA-seq to measure the steady-state transcriptome and GRO-seq to measure transcription for the first time in a hibernator. Our data capture key timepoints in both the seasonal and torpor-arousal cycles of hibernation. Strong positive correlation between transcription and the transcriptome indicates that transcriptional control dominates the known seasonal reprogramming of metabolic gene expression in liver for hibernation. During the torpor-arousal cycle, however, discordance develops between transcription and the steady-state transcriptome by at least two mechanisms: 1) although not transcribed during torpor, some transcripts are unusually stable across the torpor bout; and 2) unexpectedly, on some genes, our data suggest continuing, slow elongation with a failure to terminate transcription across the torpor bout. While the steady-state RNAs corresponding to these read through transcripts did not increase during torpor, they did increase shortly after rewarming despite their simultaneously low transcription. Both of these mechanisms would assure the immediate availability of functional transcripts upon rewarming. Integration of transcriptional, post-transcriptional and RNA stability control mechanisms, all demonstrated in these data, likely initiate a serial gene expression program across the short euthermic period that restores the tissue and prepares the animal for the next bout of torpor.

LABRAT reveals association of alternative polyadenylation with transcript localization, RNA binding protein expression, transcription speed, and cancer survival.

BACKGROUND: The sequence content of the 3' UTRs of many mRNA transcripts is regulated through alternative polyadenylation (APA). The study of this process using RNAseq data, though, has been historically challenging. RESULTS: To combat this problem, we developed LABRAT, an APA isoform quantification method. LABRAT takes advantage of newly developed transcriptome quantification techniques to accurately determine relative APA site usage and how it varies across conditions. Using LABRAT, we found consistent relationships between gene-distal APA and subcellular RNA localization in multiple cell types. We also observed connections between transcription speed and APA site choice as well as tumor-specific transcriptome-wide shifts in APA isoform abundance in hundreds of patient-derived tumor samples that were associated with patient prognosis. We investigated the effects of APA on transcript expression and found a weak overall relationship, although many individual genes showed strong correlations between relative APA isoform abundance and overall gene expression. We interrogated the roles of 191 RNA-binding proteins in the regulation of APA isoforms, finding that dozens promote broad, directional shifts in relative APA isoform abundance both in vitro and in patient-derived samples. Finally, we find that APA site shifts in the two classes of APA, tandem UTRs and alternative last exons, are strongly correlated across many contexts, suggesting that they are coregulated. CONCLUSIONS: We conclude that LABRAT has the ability to accurately quantify APA isoform ratios from RNAseq data across a variety of sample types. Further, LABRAT is able to derive biologically meaningful insights that connect APA isoform regulation to cellular and molecular phenotypes.

Quantifying alternative polyadenylation in RNAseq data with LABRAT.

Alternative polyadenylation (APA) generates transcript isoforms that differ in their 3' UTR content and may therefore be subject to different regulatory fates. Although the existence of APA has been known for decades, quantification of APA isoforms from high-throughput RNA sequencing data has been difficult. To facilitate the study of APA in large datasets, we developed an APA quantification technique called LABRAT (Lightweight Alignment-Based Reckoning of Alternative Three-prime ends). LABRAT leverages modern transcriptome quantification approaches to determine the relative abundances of APA isoforms. In this manuscript we describe how LABRAT produces its calculations, provide a step-by-step protocol for its use, and demonstrate its ability to quantify APA in single-cell RNAseq data.

The alx3 gene shapes the zebrafish neurocranium by regulating frontonasal neural crest cell differentiation timing.

During craniofacial development, different populations of cartilage- and bone-forming cells develop in precise locations in the head. Most of these cells are derived from pluripotent cranial neural crest cells and differentiate with distinct developmental timing and cellular morphologies. The mechanisms that divide neural crest cells into discrete populations are not fully understood. Here, we use single-cell RNA sequencing to transcriptomically define different populations of cranial neural crest cells. We discovered that the gene family encoding the Alx transcription factors is enriched in the frontonasal population of neural crest cells. Genetic mutant analyses indicate that alx3 functions to regulate the distinct differentiation timing and cellular morphologies among frontonasal neural crest cell subpopulations. This study furthers our understanding of how genes controlling developmental timing shape craniofacial skeletal elements.

Oxidized Low-Density Lipoprotein Drives Dysfunction of the Liver Lymphatic System.

BACKGROUND AND AIMS: As the incidence of nonalcoholic steatohepatitis (NASH) continues to rise, understanding how normal liver functions are affected during disease is required before developing novel therapeutics which could reduce morbidity and mortality. However, very little is understood about how the transport of proteins and cells from the liver by the lymphatic vasculature is affected by inflammatory mediators or during disease. METHODS: To answer these questions, we utilized a well-validated mouse model of NASH and exposure to highly oxidized low density lipoprotein (oxLDL). In addition to single cell sequencing, multiplexed immunofluorescence and metabolomic analysis of liver lymphatic endothelial cells (LEC)s we evaluated lymphatic permeability and transport both in vitro and in vivo. RESULTS: Confirming similarities between human and mouse liver lymphatic vasculature in NASH, we found that the lymphatic vasculature expands as disease progresses and results in the downregulation of genes important to lymphatic identity and function. We also demonstrate, in mice with NASH, that fluorescein isothiocyanate (FITC) dextran does not accumulate in the liver draining lymph node upon intrahepatic injection, a defect that was rescued with therapeutic administration of the lymphatic growth factor, recombinant vascular endothelial growth factor C (rVEGFC). Similarly, exposure to oxLDL reduced the amount of FITC-dextran in the portal draining lymph node and through an LEC monolayer. We provide evidence that the mechanism by which oxLDL impacts lymphatic permeability is via a reduction in Prox1 expression which decreases lymphatic specific gene expression, impedes LEC metabolism and reorganizes the highly permeable lymphatic cell-cell junctions which are a defining feature of lymphatic capillaries. CONCLUSIONS: We identify oxLDL as a major contributor to decreased lymphatic permeability in the liver, a change which is consistent with decreased protein homeostasis and increased inflammation during chronic liver disease.

Analysis of circulating breast cancer cell heterogeneity and interactions with peripheral blood mononuclear cells.

For solid tumors, extravasation of cancer cells and their survival in circulation represents a critical stage of the metastatic process that lacks complete understanding. Gaining insight into interactions between circulating tumor cells (CTCs) and other peripheral blood mononuclear cells (PBMCs) may provide valuable prognostic information. The purpose of this study was to use single-cell RNA-sequencing (scRNA-seq) of liquid biopsies from breast cancer patients to begin defining intravascular interactions. We captured CTCs from the peripheral blood of breast cancer patients using size-exclusion membranes followed by scRNA-seq of enriched CTCs and carry-over PBMCs. Transcriptome analysis identified two populations of CTCs: one enriched for transcripts indicative of estrogen responsiveness and increased proliferation and another enriched for transcripts characteristic of reduced proliferation and epithelial-mesenchymal transition (EMT). We applied interactome and pathway analysis to determine interactions between CTCs and other captured cells. Our analysis predicted for enhanced immune evasion in the CTC population with EMT characteristics. In addition, PD-1/PD-L1 pathway activation and T cell exhaustion were predicted in T cells isolated from breast cancer patients compared with normal T cells. We conclude that scRNA-seq of breast cancer CTCs generally stratifies them into two types based on their proliferative and epithelial state and differential potential to interact with PBMCs. Better understanding of CTC subtypes and their intravascular interactions may help design treatments directed against CTCs with high metastatic and immune-evasive competence.

Single-Cell RNA Sequencing of Childhood Ependymoma Reveals Neoplastic Cell Subpopulations That Impact Molecular Classification and Etiology.

Ependymoma (EPN) is a brain tumor commonly presenting in childhood that remains fatal in most children. Intra-tumoral cellular heterogeneity in bulk-tumor samples significantly confounds our understanding of EPN biology, impeding development of effective therapy. We, therefore, use single-cell RNA sequencing, histology, and deconvolution to catalog cellular heterogeneity of the major childhood EPN subgroups. Analysis of PFA subgroup EPN reveals evidence of an undifferentiated progenitor subpopulation that either differentiates into subpopulations with ependymal cell characteristics or transitions into a mesenchymal subpopulation. Histological analysis reveals that progenitor and mesenchymal subpopulations co-localize in peri-necrotic zones. In conflict with current classification paradigms, relative PFA subpopulation proportions are shown to determine bulk-tumor-assigned subgroups. We provide an interactive online resource that facilitates exploration of the EPN single-cell dataset. This atlas of EPN cellular heterogeneity increases understanding of EPN biology.

clustifyr: an R package for automated single-cell RNA sequencing cluster classification.

Assignment of cell types from single-cell RNA sequencing (scRNA-seq) data remains a time-consuming and error-prone process. Current packages for identity assignment use limited types of reference data and often have rigid data structure requirements. We developed the clustifyr R package to leverage several external data types, including gene expression profiles to assign likely cell types using data from scRNA-seq, bulk RNA-seq, microarray expression data, or signature gene lists. We benchmark various parameters of a correlation-based approach and implement gene list enrichment methods. clustifyr is a lightweight and effective cell-type assignment tool developed for compatibility with various scRNA-seq analysis workflows. clustifyr is publicly available at https://github.com/rnabioco/clustifyr.

Fibroblast subtypes define a metastatic matrisome in breast cancer.

Small primary breast cancers can show surprisingly high potential for metastasis. Clinical decision-making for tumor aggressiveness, including molecular profiling, relies primarily on analysis of the cancer cells. Here we show that this analysis is insufficient - that the stromal microenvironment of the primary tumor plays a key role in tumor cell dissemination and implantation at distant sites. We previously described 2 cancer-associated fibroblasts (CAFs) that either express (CD146+) or lack (CD146-) CD146 (official symbol MCAM, alias MUC18). We now find that when mixed with human breast cancer cells, each fibroblast subtype determines the fate of cancer cells: CD146- fibroblasts promoted increased metastasis compared with CD146+ fibroblasts. Potentially novel quantitative and qualitative proteomic analyses showed that CD146+ CAFs produced an environment rich in basement membrane proteins, while CD146- CAFs exhibited increases in fibronectin 1, lysyl oxidase, and tenascin C, all overexpressed in aggressive disease. We also show clinically that CD146- CAFs predicted for likelihood of lymph node involvement even in small primary tumors (<5 cm). Clearly small tumors enriched for CD146- CAFs require aggressive treatments.