A single-cell atlas identifies pretreatment features of primary imatinib resistance in chronic myeloid leukemia.

Primary resistance to tyrosine kinase inhibitors (TKIs) is a significant barrier to optimal outcomes in chronic myeloid leukemia (CML), but factors contributing to response heterogeneity remain unclear. Using single-cell RNA (scRNA) sequencing, we identified 8 statistically significant features in pretreatment bone marrow, which correlated with either sensitivity (major molecular response or MMR) or extreme resistance to imatinib (eventual blast crisis [BC] transformation). Employing machine-learning, we identified leukemic stem cell (LSC) and natural killer (NK) cell gene expression profiles predicting imatinib response with >80% accuracy, including no false positives for predicting BC. A canonical erythroid-specifying (TAL1/KLF1/GATA1) regulon was a hallmark of LSCs from patients with MMR and was associated with erythroid progenitor [ERP] expansion in vivo (P < .05), and a 2- to 10-fold (6.3-fold in group A vs 1.09-fold in group C) erythroid over myeloid bias in vitro. Notably, ERPs demonstrated exquisite TKI sensitivity compared with myeloid progenitors (P < .001). These LSC features were lost with progressive resistance, and MYC- and IRF1-driven inflammatory regulons were evident in patients who progressed to transformation. Patients with MMR also exhibited a 56-fold expansion (P < .01) of a normally rare subset of hyperfunctional adaptive-like NK cells, which diminished with progressive resistance, whereas patients destined for BC accumulated inhibitory NKG2A+ NK cells favoring NK cell tolerance. Finally, we developed antibody panels to validate our scRNA-seq findings. These panels may be useful for prospective studies of primary resistance, and in assessing the contribution of predetermined vs acquired factors in TKI response heterogeneity.

Single-cell and bulk transcriptome sequencing identifies two epithelial tumor cell states and refines the consensus molecular classification of colorectal cancer.

The consensus molecular subtype (CMS) classification of colorectal cancer is based on bulk transcriptomics. The underlying epithelial cell diversity remains unclear. We analyzed 373,058 single-cell transcriptomes from 63 patients, focusing on 49,155 epithelial cells. We identified a pervasive genetic and transcriptomic dichotomy of malignant cells, based on distinct gene expression, DNA copy number and gene regulatory network. We recapitulated these subtypes in bulk transcriptomes from 3,614 patients. The two intrinsic subtypes, iCMS2 and iCMS3, refine CMS. iCMS3 comprises microsatellite unstable (MSI-H) cancers and one-third of microsatellite-stable (MSS) tumors. iCMS3 MSS cancers are transcriptomically more similar to MSI-H cancers than to other MSS cancers. CMS4 cancers had either iCMS2 or iCMS3 epithelium; the latter had the worst prognosis. We defined the intrinsic epithelial axis of colorectal cancer and propose a refined 'IMF' classification with five subtypes, combining intrinsic epithelial subtype (I), microsatellite instability status (M) and fibrosis (F).

Human Nasal Epithelial Cells Sustain Persistent SARS-CoV-2 Infection In Vitro, despite Eliciting a Prolonged Antiviral Response.

The dynamics of SARS-CoV-2 infection in COVID-19 patients are highly variable, with a subset of patients demonstrating prolonged virus shedding, which poses a significant challenge for disease management and transmission control. In this study, the long-term dynamics of SARS-CoV-2 infection were investigated using a human well-differentiated nasal epithelial cell (NEC) model of infection. NECs were observed to release SARS-CoV-2 virus onto the apical surface for up to 28 days postinfection (dpi), further corroborated by viral antigen staining. Single-cell transcriptome sequencing (sc-seq) was utilized to explore the host response from infected NECs after short-term (3-dpi) and long-term (28-dpi) infection. We identified a unique population of cells harboring high viral loads present at both 3 and 28 dpi, characterized by expression of cell stress-related genes DDIT3 and ATF3 and enriched for genes involved in tumor necrosis factor alpha (TNF-α) signaling and apoptosis. Remarkably, this sc-seq analysis revealed an antiviral gene signature within all NEC cell types even at 28 dpi. We demonstrate increased replication of basal cells, absence of widespread cell death within the epithelial monolayer, and the ability of SARS-CoV-2 to replicate despite a continuous interferon response as factors likely contributing to SARS-CoV-2 persistence. This study provides a model system for development of therapeutics aimed at improving viral clearance in immunocompromised patients and implies a crucial role for immune cells in mediating viral clearance from infected epithelia. IMPORTANCE Increasing medical attention has been drawn to the persistence of symptoms (long-COVID syndrome) or live virus shedding from subsets of COVID-19 patients weeks to months after the initial onset of symptoms. In vitro approaches to model viral or symptom persistence are needed to fully dissect the complex and likely varied mechanisms underlying these clinical observations. We show that in vitro differentiated human NECs are persistently infected with SARS-CoV-2 for up to 28 dpi. This viral replication occurred despite the presence of an antiviral gene signature across all NEC cell types even at 28 dpi. This indicates that epithelial cell intrinsic antiviral responses are insufficient for the clearance of SARS-CoV-2, implying an essential role for tissue-resident and infiltrating immune cells for eventual viral clearance from infected airway tissue in COVID-19 patients.

RCA2: a scalable supervised clustering algorithm that reduces batch effects in scRNA-seq data.

The transcriptomic diversity of cell types in the human body can be analysed in unprecedented detail using single cell (SC) technologies. Unsupervised clustering of SC transcriptomes, which is the default technique for defining cell types, is prone to group cells by technical, rather than biological, variation. Compared to de-novo (unsupervised) clustering, we demonstrate using multiple benchmarks that supervised clustering, which uses reference transcriptomes as a guide, is robust to batch effects and data quality artifacts. Here, we present RCA2, the first algorithm to combine reference projection (batch effect robustness) with graph-based clustering (scalability). In addition, RCA2 provides a user-friendly framework incorporating multiple commonly used downstream analysis modules. RCA2 also provides new reference panels for human and mouse and supports generation of custom panels. Furthermore, RCA2 facilitates cell type-specific QC, which is essential for accurate clustering of data from heterogeneous tissues. We demonstrate the advantages of RCA2 on SC data from human bone marrow, healthy PBMCs and PBMCs from COVID-19 patients. Scalable supervised clustering methods such as RCA2 will facilitate unified analysis of cohort-scale SC datasets.

Spearmint R2R3-MYB transcription factor MsMYB negatively regulates monoterpene production and suppresses the expression of geranyl diphosphate synthase large subunit (MsGPPS.LSU).

Many aromatic plants, such as spearmint, produce valuable essential oils in specialized structures called peltate glandular trichomes (PGTs). Understanding the regulatory mechanisms behind the production of these important secondary metabolites will help design new approaches to engineer them. Here, we identified a PGT-specific R2R3-MYB gene, MsMYB, from comparative RNA-Seq data of spearmint and functionally characterized it. Analysis of MsMYB-RNAi transgenic lines showed increased levels of monoterpenes, and MsMYB-overexpressing lines exhibited decreased levels of monoterpenes. These results suggest that MsMYB is a novel negative regulator of monoterpene biosynthesis. Ectopic expression of MsMYB, in sweet basil and tobacco, perturbed sesquiterpene- and diterpene-derived metabolite production. In addition, we found that MsMYB binds to cis-elements of MsGPPS.LSU and suppresses its expression. Phylogenetic analysis placed MsMYB in subgroup 7 of R2R3-MYBs whose members govern phenylpropanoid pathway and are regulated by miR858. Analysis of transgenic lines showed that MsMYB is more specific to terpene biosynthesis as it did not affect metabolites derived from phenylpropanoid pathway. Further, our results indicate that MsMYB is probably not regulated by miR858, like other members of subgroup 7.

Metabolic engineering of terpene biosynthesis in plants using a trichome-specific transcription factor MsYABBY5 from spearmint (Mentha spicata).

In many aromatic plants including spearmint (Mentha spicata), the sites of secondary metabolite production are tiny specialized structures called peltate glandular trichomes (PGT). Having high commercial values, these secondary metabolites are exploited largely as flavours, fragrances and pharmaceuticals. But, knowledge about transcription factors (TFs) that regulate secondary metabolism in PGT remains elusive. Understanding the role of TFs in secondary metabolism pathway will aid in metabolic engineering for increased yield of secondary metabolites and also the development of new production techniques for valuable metabolites. Here, we isolated and functionally characterized a novel MsYABBY5 gene that is preferentially expressed in PGT of spearmint. We generated transgenic plants in which MsYABBY5 was either overexpressed or silenced using RNA interference (RNAi). Analysis of the transgenic lines showed that the reduced expression of MsYABBY5 led to increased levels of terpenes and that overexpression decreased terpene levels. Additionally, ectopic expression of MsYABBY5 in Ocimum basilicum and Nicotiana sylvestris decreased secondary metabolite production in them, suggesting that the encoded transcription factor is probably a repressor of secondary metabolism.

A terpenoid phytoalexin plays a role in basal defense of Nicotiana benthamiana against Potato virus X.

Terpenoid phytoalexins function as defense compound against a broad spectrum of pathogens and pests in the plant kingdom. However, the role of phytoalexin in antiviral defense is still elusive. In this study, we identified the biosynthesis pathway of a sesquiterpenoid phytoalexin, capsidiol 3-acetate as an antiviral response against RNA virus Potato Virus X (PVX) in Nicotiana benthamiana. NbTPS1 and NbEAH genes were found strongly induced by PVX-infection. Enzymatic activity and genetic evidence indicated that both genes were involved in the PVX-induced biosynthesis of capsidiol 3-acetate. NbTPS1- or NbEAH-silenced plant was more susceptible to PVX. The accumulation of capsidiol 3-acetate in PVX-infected plant was partially regulated by jasmonic acid signaling receptor COI1. These findings provide an insight into a novel mechanism of how plant uses the basal arsenal machinery to mount a fight against virus attack even in susceptible species.

Over-expression of OSRIP18 increases drought and salt tolerance in transgenic rice plants.

Both drought and high salinity stresses are major abiotic factors that limit the yield of agricultural crops. Transgenic techniques have been regarded as effective ways to improve crops in their tolerance to these abiotic stresses. Functional characterization of genes is the prerequisite to identify candidates for such improvement. Here, we have investigated the biological functions of an Oryza sativa Ribosome-inactivating protein gene 18 (OSRIP18) by ectopically expressing this gene under the control of CaMV 35S promoter in the rice genome. We have generated 11 independent transgenic rice plants and all of them showed significantly increased tolerance to drought and high salinity stresses. Global gene expression changes by Microarray analysis showed that more than 100 probe sets were detected with up-regulated expression abundance while signals from only three probe sets were down-regulated after over-expression of OSRIP18. Most of them were not regulated by drought or high salinity stresses. Our data suggested that the increased tolerance to these abiotic stresses in transgenic plants might be due to up-regulation of some stress-dependent/independent genes and OSRIP18 may be potentially useful in further improving plant tolerance to various abiotic stresses by over-expression.