Single-cell Transcriptional Profiling Research Articles

Distinct single cell transcriptional profile in CD4+ T-lymphocytes among obese children with asthma.

Introduction: Obesity is a risk factor for asthma morbidity, associated with less responsiveness to inhaled corticosteroids. CD4+ T-cells are central to the immunology of asthma and may contribute to the unique obese asthma phenotype. We sought to characterize the single cell CD4+ Transcriptional profile differences in obese children with asthma compared to normal weight children with asthma. Methods: Eight normal weight and obese participants with asthma were clinically phenotyped and matched based on asthma control. Peripheral blood (PB) CD4+ T-cells were sorted, and single cell RNA sequencing was conducted. Cell clusters were identified by canonical gene expression, differential gene expression and reactome pathway analysis was applied. The obese PB bulk transcriptomic signature from the U-BIOPRED pediatric cohort was assessed in our cohort as well. Results: Obese children with asthma have a distinct CD4+ transcriptional profile with differential gene expression. There were more activated protein tyrosine phosphate receptor type C (PTPRC)high cells and less PTPRClow in the obese children. Obese children had higher enrichment of the neutrophil degranulation, interleukin-7 (IL-7) receptor and IL-7-related janus kinase-signal transducer and activator of transcription signaling pathways. Genes previously associated with more severe asthma, IL-32, FKBP5 gene expression, IL-6 and Rho transcriptional signaling, were also enriched in obese children with asthma. Discussion: Our findings shed insight into the molecular mechanisms underpinning more severe and steroid-resistant asthma among obese children. Further investigation is needed to identify potential new therapeutic targets for this group.

Enhancer-driven Shh signaling promotes glia-to-mesenchyme transition during bone repair

Plp1-lineage Schwann cells (SCs) of peripheral nerve play a critical role in vascular remodeling and osteogenic differentiation during the early stage of bone healing, and the abnormal plasticity of SCs would jeopardize the bone regeneration. However, how Plp1-lineage cells respond to injury and initiate the vascularized osteogenesis remains incompletely understood. Here, by employing single-cell transcriptional profiling combined with lineage-specific tracing models, we uncover that Plp1-lineage cells undergoing injury-induced glia-to-MSCs transition contributed to osteogenesis and revascularization in the initial stage of bone injury. Importantly, our data demonstrated that the Sonic hedgehog (Shh) signaling was responsible for the transition process initiation, which was strongly activated by c-Jun/SIRT6/BAF170 complex-driven Shh enhancers. Collectively, these findings depict an injury-specific niche signal-mediated Plp1-lineage cells transition towards Gli1+ MSCs and may be instructive for approaches to promote bone regeneration during aging or other bone diseases.

Hypoxia is linked to acquired resistance to immune checkpoint inhibitors in lung cancer.

Despite the established use of immune checkpoint inhibitors (ICIs) to treat non-small cell lung cancer (NSCLC), only a subset of patients benefit from treatment and ∼50% of patients whose tumors respond eventually develop acquired resistance (AR). To identify novel drivers of AR, we generated murine Msh2 knock-out (KO) lung tumors that initially responded but eventually developed AR to anti-PD-1, alone or in combination with anti-CTLA-4. Resistant tumors harbored decreased infiltrating T cells and reduced cancer cell-intrinsic MHC-I and MHC-II levels, yet remained responsive to IFNγ. Resistant tumors contained extensive regions of hypoxia, and a hypoxia signature derived from single-cell transcriptional profiling of resistant cancer cells was associated with decreased progression-free survival in a cohort of NSCLC patients treated with anti-PD-1/PD-L1 therapy. Targeting hypoxic tumor regions using a hypoxia-activated pro-drug delayed AR to ICIs in murine Msh2 KO tumors. Thus, this work provides a rationale for targeting tumor metabolic features, such as hypoxia, in combination with immune checkpoint inhibition.

Bulk and single-cell transcriptome revealed the metabolic heterogeneity in human glioma.

Emerging perspectives on tumor metabolism reveal its heterogeneity, a characteristic yet to be fully explored in gliomas. To advance therapies targeting metabolic processes, it is crucial to uncover metabolic differences and identify distinct metabolic subtypes. Therefore, we aimed to develop a classification system for gliomas based on the enrichment levels of four key metabolic pathways: glutaminolysis, glycolysis, the pentose phosphate pathway, and fatty acid oxidation. Energy-related features of glioma were characterized through integrative analyses of multiple datasets, including bulk, single-cell, and spatial transcriptome profiling. The glioma energy metabolic subtypes were constructed using the R package ConsensusClusterPlus. Kaplan-Meier analysis was conducted to compare clinical outcomes between different metabolic groups. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were employed to elucidate the biological functions of genes of interest. Cell-cell communication analysis was performed at single-cell resolution using the R package CellChat and at spatial resolution using the standard stLearn pipeline. Glioma samples were stratified into two prognostic subtypes. Group 1, enriched in the glutaminolysis pathway, had better clinical outcomes. In contrast, Group 2 exhibited high activities in glycolysis, the pentose phosphate pathway, and fatty acid oxidation, correlating with decreased survival time. Group 1 samples were predominantly located in the peripheral region and had a high composition of neuron cells. Group 2, however, had increased infiltration of tumor-promoting immune cells, such as M2 macrophages, and was characterized by traits of invasion, hypoxia, and immunity. Lastly, cell-cell communications were compared across different tumor regions, and the CX3CL1/CX3CR1 ligand-receptor pair was validated using spatial transcriptomic data. Our work revealed the metabolic heterogeneity in glioma by developing a new classification system with significant prognostic and therapeutic value. Single-cell transcriptional profiles offer novel insights into tumor metabolic reprogramming, which could enhance therapies tailored to cell- or patient-specific metabolic patterns.

Diagnosing recipient- vs. donor-derived posttransplant myelodysplastic neoplasm via targeted single-cell mutational profiling.

Distinguishing donor- vs. recipient-derived myelodysplastic neoplasm (MDS) after allogeneic hematopoietic stem cell transplantation (allo-HSCT) is challenging and has direct therapeutical implications. Here, we took a translational approach that we used in addition to conventional diagnostic techniques to resolve the origin of MDS in a 38-year-old patient with acquired aplastic anemia and evolving MDS after first allo-HSCT. Specifically, we used single-cell transcriptional profiling to differentiate between donor- and recipient-derived bone marrow cells and established a strategy that additionally allows identification of cells carrying the MDS-associated U2AF1S34Y variant. The patient exhibited mixed donor chimerism combined with severely reduced erythropoiesis and dysplastic morphology within the granulocytic and megakaryocytic lineage along with the MDS-associated U2AF1S34Y mutation in the bone marrow. Single-cell transcriptional profiling together with targeted enrichment of the U2AF1S34Y-specific locus further revealed that, while the immune compartment was mainly populated by donor-derived cells, myelopoiesis was predominantly driven by the recipient. Additionally, concordant with recipient-derived MDS, we found that U2AF1S34Y-mutated cells were exclusively recipient derived with X but not Y chromosome-specific gene expression. Our study highlights the clinical potential of integrating high-resolution single-cell techniques to resolve complex cases for personalized treatment decisions. The study was funded by intramural resources of the Charité- Universitätsmedizin Berlin and the Berlin Institute of Health.

Single cell RNA-sequencing reveals no evidence for meiotic sex chromosome inactivation in the threespine stickleback fish.

As male germ cells enter meiosis, the X and Y chromosome of many species undergo a drastic repression of gene expression. In mammals, this process has been shown to be essential for fertility, and the expression of sex-linked genes can lead to meiotic arrest and cell death. This process has only been studied in a handful of organisms, which limits our understanding how conserved MSCI is across the tree of life. Teleost fish are an understudied group with many examples of independently derived sex chromosomes across closely related species. Here, we investigate whether MSCI occurs in the threespine stickleback fish, using single-cell transcriptional profiling. We found gene expression remains active throughout meiosis on the sex chromosomes, indicating MSCI does not occur. This indicates that MSCI is not a conserved feature of all taxa and is not an inevitable outcome of degenerating Y chromosomes.

IMMU-01. DE NOVO-DESIGNED BINDER-BASED CAR T CELLS ENABLE EFFECTIVE TARGETING OF GLIOBLASTOMA

Abstract Glioblastoma is a devastating disease with poor responses to all current therapies. Chimeric antigen receptor (CAR) T cell therapy has achieved impressive success against hematological malignancies but produces responses in only a limited number of glioblastoma patients. Here, we developed a novel platform for CAR constructs using de novo designed protein binders to target tumor antigens. We firstly designed a binder CAR targeting epidermal growth factor receptor (EGFR), which is highly expressed in glioblastoma, to validate the feasibility of our platform. Our EGFR binder CAR design promoted CAR T cell proliferation, cytotoxic cytokine secretion and exhaustion resistance, leading to better antitumor performance both in vitro and in vivo. Bulk and single-cell transcriptional profiling of binder CAR T cells revealed an underlying mechanism involving enhanced effector activity and reduced exhaustion responses. Mechanismly, the binder CARs exhibited higher surface expression and greater resistance to degradation triggered by tumor antigens. Moreover, we designed and engineered a novel CD276 binder CAR T cell product and verified its antitumor function, further demonstrating the general applicability of our binder-based CAR platform. Overall, our study provides an alternative avenue for the design of a CAR antigen-binding domain to potentiate CAR T cell antitumor efficacy.

MEK inhibition prevents CAR-T cell exhaustion and differentiation via downregulation of c-Fos and JunB

Clinical evidence supports the notion that T cell exhaustion and terminal differentiation pose challenges to the persistence and effectiveness of chimeric antigen receptor-T (CAR-T) cells. MEK1/2 inhibitors (MEKIs), widely used in cancer treatment due to their ability to inhibit aberrant MAPK signaling, have shown potential synergistic effects when combined with immunotherapy. However, the impact and mechanisms of MEKIs on CAR-T cells remain uncertain and controversial. To address this, we conducted a comprehensive investigation to determine whether MEKIs enhance or impair the efficacy of CAR-T cells. Our findings revealed that MEKIs attenuated CAR-T cell exhaustion and terminal differentiation induced by tonic signaling and antigen stimulation, thereby improving CAR-T cell efficacy against hematological and solid tumors. Remarkably, these effects were independent of the specific scFvs and costimulatory domains utilized in CARs. Mechanistically, analysis of bulk and single-cell transcriptional profiles demonstrates that the effect of MEK inhibition was related to diminish anabolic metabolism and downregulation of c-Fos and JunB. Additionally, the overexpression of c-Fos or JunB in CAR-T cells counteracted the effects of MEK inhibition. Furthermore, our Cut-and-Tag assay revealed that MEK inhibition downregulated the JunB-driven gene profiles associated with exhaustion, differentiation, anergy, glycolysis, and apoptosis. In summary, our research unveil the critical role of the MAPK-c-Fos-JunB axis in driving CAR-T cell exhaustion and terminal differentiation. These mechanistic insights significantly broaden the potential application of MEKIs to enhance the effectiveness of CAR-T therapy.

Targeting overexpressed antigens in glioblastoma via CAR T cells with computationally designed high-affinity protein binders.

Chimeric antigen receptor (CAR) T cells targeting receptors on tumour cells have had limited success in patients with glioblastoma. Here we report the development and therapeutic performance of CAR constructs leveraging protein binders computationally designed de novo to have high affinity for the epidermal growth factor receptor (EGFR) or the tumour-associated antigen CD276, which are overexpressed in glioblastoma. With respect to T cells with a CAR using an antibody-derived single-chain variable fragment as antigen-binding domain, the designed binders on CAR T cells promoted the proliferation of the cells, the secretion of cytotoxic cytokines and their resistance to cell exhaustion, and improved antitumour performance in vitro and in vivo. Moreover, CARs with the binders exhibited higher surface expression and greater resistance to degradation, as indicated by bulk and single-cell transcriptional profiling of the cells. The de novo design of binding domains for specific tumour antigens may potentiate the antitumour efficacy of CAR T cell therapies for other solid cancers.

Single-cell transcriptional profiling reveals cell type-specific responses to duck reovirus infection in the Bursa of Fabricius of Cairna moschata

Duck reovirus (DRV) is a universal waterfowl virus that causes significant economic losses in the duck industry. However, the role of the host innate immune response of the Bursa of Fabricius to DRV infection is largely unknown. In the present study, we constructed a single-cell resolution transcriptomic atlas of the Bursa of Fabricius of Cairna moschata after infection with HN10 (a novel DRV). Ten cell-type marker genes were used to annotate the cell type, indicating a high degree of cell heterogeneity in the Bursa of Fabricius. Most of the innate and adaptive immune system-related genes were highly expressed in T cells, B cells, neutrophils, macrophages, and DCs. In the Bursa of Fabricius, the proportions of DCs and macrophages were largely increased by HN10 infection at 14 d, suggesting that DCs and macrophages play important roles in the long-term viral response. Notably, a number of innate and adaptive immune system-related genes were highly expressed at 24 h after HN10 infection, indicating that the Bursa of Fabricius has a very strong immune function even in the early developmental stage. In the immune system, the NOD-like receptor signaling pathway and RIG-I-like receptor signaling pathway were significantly activated at the early stage of HN10 infection, while the Toll-like receptor signaling pathway was significantly activated at the late stage. Enrichment analysis suggested that different immune signaling pathways play roles in specific developmental stages. Our data provide an opportunity to reveal the immune response to DRV infection at the single-cell level.

Single-cell and bulk transcriptional profiling of mouse ovaries reveals novel genes and pathways associated with DNA damage response in oocytes

Immature oocytes enclosed in primordial follicles stored in female ovaries are under constant threat of DNA damage induced by endogenous and exogenous factors. Checkpoint kinase 2 (CHEK2) is a key mediator of the DNA damage response (DDR) in all cells. Genetic studies have shown that CHEK2 and its downstream targets, p53, and TAp63, regulate primordial follicle elimination in response to DNA damage. However, the mechanism leading to their demise is still poorly characterized. Single-cell and bulk RNA sequencing were used to determine the DDR in wild-type and Chek2-deficient ovaries. A low but oocyte-lethal dose of ionizing radiation induces ovarian DDR that is solely dependent on CHEK2. DNA damage activates multiple response pathways related to apoptosis, p53, interferon signaling, inflammation, cell adhesion, and intercellular communication. These pathways are differentially employed by different ovarian cell types, with oocytes disproportionately affected by radiation. Novel genes and pathways are induced by radiation specifically in oocytes, shedding light on their sensitivity to DNA damage, and implicating a coordinated response between oocytes and pregranulosa cells within the follicle. These findings provide a foundation for future studies on the specific mechanisms regulating oocyte survival in the context of aging, therapeutic and environmental genotoxic exposures.

VEGF-A-mediated venous endothelial cell proliferation results in neoangiogenesis during neuroinflammation.

Newly formed leaky vessels and blood-brain barrier (BBB) damage are present in demyelinating acute and chronic lesions in multiple sclerosis (MS) and experimental autoimmune encephalomyelitis (EAE). However, the endothelial cell subtypes and signaling pathways contributing to these leaky neovessels are unclear. Here, using single-cell transcriptional profiling and in vivo validation studies, we show that venous endothelial cells express neoangiogenesis gene signatures and show increased proliferation resulting in enlarged veins and higher venous coverage in acute and chronic EAE lesions in female adult mice. These changes correlate with the upregulation of vascular endothelial growth factor A (VEGF-A) signaling. We also confirmed increased expression of neoangiogenic markers in acute and chronic human MS lesions. Treatment with a VEGF-A blocking antibody diminishes the neoangiogenic transcriptomic signatures and vascular proliferation in female adult mice with EAE, but it does not restore BBB function or ameliorate EAE pathology. Our data demonstrate that venous endothelial cells contribute to neoangiogenesis in demyelinating neuroinflammatory conditions.

An Efficient Fluorescence-Activated Protoplast Sorting (FAPS) and Regeneration Protocol for Canola (Brassica napus).

Protoplast sorting and purification methods are powerful tools enabling the enrichment of cellular subpopulations for basic and applied studies in plant sciences. Fluorescence-activated protoplast sorting (FAPS) is an efficient method to isolate specific protoplast populations based on innate features (size and autofluorescence) or expression of fluorescent proteins. FAPS-based methods have recently been deployed in single-cell purification for single-cell RNA sequencing-based transcriptional profiling studies. Protoplast sorting methods integrated with the ability to culture and recover whole plants add value to functional genomics and gene editing applications. Enriching cells expressing nucleases linked to fluorescent proteins can maximize knockout or knockin editing efficiencies and minimize toxic and off-target effects. Here, we report the protocol for protoplast preparation, sterile cell sorting, culture, and downstream regeneration of plants from canola protoplasts. This protocol can be successfully applied to all totipotent protoplast methods that can regenerate into whole plants. © 2024 Wiley Periodicals LLC. Basic Protocol 1: Preparation of transfected canola protoplasts for sorting Basic Protocol 2: Fluorescence-activated protoplast sorting Basic Protocol 3: Bead culture of sorted protoplasts and recovery of plantlets.

A glia-enriched stem cell 3D model of the human brain mimics the glial-immune neurodegenerative phenotypes of multiple sclerosis

The role of central nervous system (CNS) glia in sustaining self-autonomous inflammation and driving clinical progression in multiple sclerosis (MS) is gaining scientific interest. We applied a single transcription factor (SOX10)-based protocol to accelerate oligodendrocyte differentiation from human induced pluripotent stem cell (hiPSC)-derived neural precursor cells, generating self-organizing forebrain organoids. These organoids include neurons, astrocytes, oligodendroglia, and hiPSC-derived microglia to achieve immunocompetence. Over 8weeks, organoids reproducibly generated mature CNS cell types, exhibiting single-cell transcriptional profiles similar to the adult human brain. Exposed to inflamed cerebrospinal fluid (CSF) from patients with MS, organoids properly mimic macroglia-microglia neurodegenerative phenotypes and intercellular communication seen in chronic active MS. Oligodendrocyte vulnerability emerged by day 6 post-MS-CSF exposure, with nearly 50% reduction. Temporally resolved organoid data support and expand on the role of soluble CSF mediators in sustaining downstream events leading to oligodendrocyte death and inflammatory neurodegeneration. Such findings support the implementation of this organoid model for drug screening to halt inflammatory neurodegeneration.

Single-cell RNA-sequencing reveals Eriocheir sinensis hemocyte subpopulations and their molecular responses to Spiroplasma eriocheiris infection

The hemocytes is the central part of the Eriocheir sinensis innate immune system, however, their subclasses have not been completely defined and the single-cell gene transcriptional profiles and functions are still unknown. In the present study, the single-cell mRNA sequencing (scRNA-seq) was carried out to build a single-cell gene expression map of E. sinensis hemocytes and further analysis of the interaction molecular mechanism between the crab and novel aquatic pathogens Spiroplasma eriocheiris. By scRNA-seq, we successfully identified nine cell clusters (Hem1–9) and predicted their functions according to the specific marker genes. Pseudotime inference analysis revealed that all the identified cells undergo three stages, state 1 (the initial point), state 2 and state 3 (the other two branches). After being infected by S. eriocheiris, lots of different expression genes (DEGs) were identified in each cell cluster compared with the control group (with fold change ≥2 and p-value <0.05). The bioinformatics analysis found that the Hem1 and Hem5 were the main cell clusters synthesis and release of antimicrobial peptides (AMPs) for help the host against the pathogen infection. The Hem5 also could activate the proPO system to participate in the process of host eliminate the invading S. eriocheiris. Hem7, Hem8, and Hem9 were mainly responsible for killing the invading pathogen through phagocytosis. In particular, there were also many immune-related genes identified, including β-glucan binding protein (BGBP), Clip-domain serine proteases (CLIPs), phenoloxidase 1 (PO1), anti-lipopolysaccharide factors (ALFs), syndecan, and lysozyme C, which not only played critical roles in the crab immune system but also could serve as the marker genes for identifying cell subgroup. In conclusion, the results obtained in the present study could provide pivotal insights into the crab hemocytes cell-type-specific molecular mechanisms response S. eriocheiris infection at a single-cell level.

Blockade of IL1β and PD1 with Combination Chemotherapy Reduces Systemic Myeloid Suppression in Metastatic Pancreatic Cancer with Heterogeneous Effects in the Tumor.

Innate inflammation promotes tumor development, although the role of innate inflammatory cytokines in established human tumors is unclear. Herein, we report clinical and translational results from a phase Ib trial testing whether IL1β blockade in human pancreatic cancer would alleviate myeloid immunosuppression and reveal antitumor T-cell responses to PD1 blockade. Patients with treatment-naïve advanced pancreatic ductal adenocarcinoma (n = 10) were treated with canakinumab, a high-affinity monoclonal human antiinterleukin-1β (IL1β), the PD1 blocking antibody spartalizumab, and gemcitabine/n(ab)paclitaxel. Analysis of paired peripheral blood from patients in the trial versus patients receiving multiagent chemotherapy showed a modest increase in HLA-DR+CD38+ activated CD8+ T cells and a decrease in circulating monocytic myeloid-derived suppressor cells (MDSC) by flow cytometry for patients in the trial but not in controls. Similarly, we used patient serum to differentiate monocytic MDSCs in vitro and showed that functional inhibition of T-cell proliferation was reduced when using on-treatment serum samples from patients in the trial but not when using serum from patients treated with chemotherapy alone. Within the tumor, we observed few changes in suppressive myeloid-cell populations or activated T cells as assessed by single-cell transcriptional profiling or multiplex immunofluorescence, although increases in CD8+ T cells suggest that improvements in the tumor immune microenvironment might be revealed by a larger study. Overall, the data indicate that exposure to PD1 and IL1β blockade induced a modest reactivation of peripheral CD8+ T cells and decreased circulating monocytic MDSCs; however, these changes did not lead to similarly uniform alterations in the tumor microenvironment.

Single-cell transcriptional profiling of clear cell renal cell carcinoma reveals a tumor-associated endothelial tip cell phenotype

Clear cell renal cell carcinoma (ccRCC) is the most prevalent form of renal cancer, accounting for over 75% of cases. The asymptomatic nature of the disease contributes to late-stage diagnoses and poor survival. Highly vascularized and immune infiltrated microenvironment are prominent features of ccRCC, yet the interplay between vasculature and immune cells, disease progression and response to therapy remains poorly understood. Using droplet-based single-cell RNA sequencing we profile 50,236 transcriptomes from paired tumor and healthy adjacent kidney tissues. Our analysis reveals significant heterogeneity and inter-patient variability of the tumor microenvironment. Notably, we discover a previously uncharacterized vasculature subpopulation associated with epithelial-mesenchymal transition. The cell-cell communication analysis reveals multiple modes of immunosuppressive interactions within the tumor microenvironment, including clinically relevant interactions between tumor vasculature and stromal cells with immune cells. The upregulation of the genes involved in these interactions is associated with worse survival in the TCGA KIRC cohort. Our findings demonstrate the role of tumor vasculature and stromal cell populations in shaping the ccRCC microenvironment and uncover a subpopulation of cells within the tumor vasculature that is associated with an angiogenic phenotype.

Using Human Intestinal Organoids to Understand the Small Intestine Epithelium at the Single Cell Transcriptional Level.

Single cell transcriptomics has revolutionized our understanding of the cell biology of the human body. State-of-the-art human small intestinal organoid cultures provide ex vivo model systems that bridge the gap between animal models and clinical studies. The application of single cell transcriptomics to human intestinal organoid (HIO) models is revealing previously unrecognized cell biology, biochemistry, and physiology of the GI tract. The advanced single cell transcriptomics platforms use microfluidic partitioning and barcoding to generate cDNA libraries. These barcoded cDNAs can be easily sequenced by next generation sequencing platforms and used by various visualization tools to generate maps. Here, we describe methods to culture and differentiate human small intestinal HIOs in different formats and procedures for isolating viable cells from these formats that are suitable for use in single-cell transcriptional profiling platforms. These protocols and procedures facilitate the use of small intestinal HIOs to obtain an increased understanding of the cellular response of human intestinal epithelium at the transcriptional level in the context of a variety of different environments.

Single-cell transcriptional profile of CD34+ hematopoietic progenitor cells from del(5q) myelodysplastic syndromes and impact of lenalidomide

While myelodysplastic syndromes with del(5q) (del(5q) MDS) comprises a well-defined hematological subgroup, the molecular basis underlying its origin remains unknown. Using single cell RNA-seq (scRNA-seq) on CD34+ progenitors from del(5q) MDS patients, we have identified cells harboring the deletion, characterizing the transcriptional impact of this genetic insult on disease pathogenesis and treatment response. Interestingly, both del(5q) and non-del(5q) cells present similar transcriptional lesions, indicating that all cells, and not only those harboring the deletion, may contribute to aberrant hematopoietic differentiation. However, gene regulatory network (GRN) analyses reveal a group of regulons showing aberrant activity that could trigger altered hematopoiesis exclusively in del(5q) cells, pointing to a more prominent role of these cells in disease phenotype. In del(5q) MDS patients achieving hematological response upon lenalidomide treatment, the drug reverts several transcriptional alterations in both del(5q) and non-del(5q) cells, but other lesions remain, which may be responsible for potential future relapses. Moreover, lack of hematological response is associated with the inability of lenalidomide to reverse transcriptional alterations. Collectively, this study reveals transcriptional alterations that could contribute to the pathogenesis and treatment response of del(5q) MDS.

High-throughput single-cell transcriptomics of bacteria using combinatorial barcoding.

Microbial split-pool ligation transcriptomics (microSPLiT) is a high-throughput single-cell RNA sequencing method for bacteria. With four combinatorial barcoding rounds, microSPLiT can profile transcriptional states in hundreds of thousands of Gram-negative and Gram-positive bacteria in a single experiment without specialized equipment. As bacterial samples are fixed and permeabilized before barcoding, they can be collected and stored ahead of time. During the first barcoding round, the fixed and permeabilized bacteria are distributed into a 96-well plate, where their transcripts are reverse transcribed into cDNA and labeled with the first well-specific barcode inside the cells. The cells are mixed and redistributed two more times into new 96-well plates, where the second and third barcodes are appended to the cDNA via in-cell ligation reactions. Finally, the cells are mixed and divided into aliquot sub-libraries, which can be stored until future use or prepared for sequencing with the addition of a fourth barcode. It takes 4 days to generate sequencing-ready libraries, including 1 day for collection and overnight fixation of samples. The standard plate setup enables single-cell transcriptional profiling of up to 1 million bacterial cells and up to 96 samples in a single barcoding experiment, with the possibility of expansion by adding barcoding rounds. The protocol requires experience in basic molecular biology techniques, handling of bacterial samples and preparation of DNA libraries for next-generation sequencing. It can be performed by experienced undergraduate or graduate students. Data analysis requires access to computing resources, familiarity with Unix command line and basic experience with Python or R.