Single-cell RNA-sequencing (scRNA-seq) has emerged as a revolutionary tool that allows us to address scientific questions that eluded examination just a few years ago. With the advantages of ...scRNA-seq come computational challenges that are just beginning to be addressed. In this article, we highlight the computational methods available for the design and analysis of scRNA-seq experiments, their advantages and disadvantages in various settings, the open questions for which novel methods are needed, and expected future developments in this exciting area.
The normalization of RNA-seq data is essential for accurate downstream inference, but the assumptions upon which most normalization methods are based are not applicable in the single-cell setting. ...Consequently, applying existing normalization methods to single-cell RNA-seq data introduces artifacts that bias downstream analyses. To address this, we introduce SCnorm for accurate and efficient normalization of single-cell RNA-seq data.
From bacteria to humans, individual cells within isogenic populations can show significant variation in stress tolerance, but the nature of this heterogeneity is not clear. To investigate this, we ...used single-cell RNA sequencing to quantify transcript heterogeneity in single Saccharomyces cerevisiae cells treated with and without salt stress to explore population variation and identify cellular covariates that influence the stress-responsive transcriptome. Leveraging the extensive knowledge of yeast transcriptional regulation, we uncovered significant regulatory variation in individual yeast cells, both before and after stress. We also discovered that a subset of cells appears to decouple expression of ribosomal protein genes from the environmental stress response in a manner partly correlated with the cell cycle but unrelated to the yeast ultradian metabolic cycle. Live-cell imaging of cells expressing pairs of fluorescent regulators, including the transcription factor Msn2 with Dot6, Sfp1, or MAP kinase Hog1, revealed both coordinated and decoupled nucleocytoplasmic shuttling. Together with transcriptomic analysis, our results suggest that cells maintain a cellular filter against decoupled bursts of transcription factor activation but mount a stress response upon coordinated regulation, even in a subset of unstressed cells.
Regulatory T (Treg) cells are one of the major immunosuppressive cell types in cancer and a potential target for immunotherapy, but targeting tumor-infiltrating (TI) Treg cells has been challenging. ...Here, using single-cell RNA sequencing of immune cells from renal clear cell carcinoma (ccRCC) patients, we identify two distinct transcriptional fates for TI Treg cells, Fate-1 and Fate-2. The Fate-1 signature is associated with a poorer prognosis in ccRCC and several other solid cancers. CD177, a cell surface protein normally expressed on neutrophil, is specifically expressed on Fate-1 TI Treg cells in several solid cancer types, but not on other TI or peripheral Treg cells. Mechanistically, blocking CD177 reduces the suppressive activity of Treg cells in vitro, while Treg-specific deletion of Cd177 leads to decreased tumor growth and reduced TI Treg frequency in mice. Our results thus uncover a functional CD177
TI Treg population that may serve as a target for TI Treg-specific immunotherapy.
Defects in somitogenesis result in vertebral malformations at birth known as spondylocostal dysostosis (SCDO). Somites are formed with a species-specific periodicity controlled by the “segmentation ...clock,” which comprises a group of oscillatory genes in the presomitic mesoderm. Here, we report that a segmentation clock model derived from human embryonic stem cells shows many hallmarks of the mammalian segmentation clock in vivo, including a dependence on the NOTCH and WNT signaling pathways. The gene expression oscillations are highly synchronized, displaying a periodicity specific to the human clock. Introduction of a point of mutation into HES7, a specific mutation previously associated with clinical SCDO, eliminated clock gene oscillations, successfully reproducing the defects in the segmentation clock. Thus, we provide a model for studying the previously inaccessible human segmentation clock to better understand the mechanisms contributing to congenital skeletal defects.
Display omitted
•We captured a transitory somitogenesis gene-expression signature in vitro•Human ESC-derived presomitic mesoderm exhibits synchronized HES7 oscillation•WNT activation and TGF-β inhibition are required to propagate HES7 oscillation•We established a model for spondylocostal dysostosis
The segmentation clock is a molecular oscillator regulating the tempo of somite formation in a species-specific manner. Chu et al. report an embryonic-stem-cell-derived model system displaying a human-specific gene oscillation periodicity, which can shed light on human somitogenesis and model skeletal developmental disorders.
High-throughput expression profiling experiments with ordered conditions (e.g. time-course or spatial-course) are becoming more common for studying detailed differentiation processes or spatial ...patterns. Identifying dynamic changes at both the individual gene and whole transcriptome level can provide important insights about genes, pathways, and critical time points.
We present an R package, Trendy, which utilizes segmented regression models to simultaneously characterize each gene's expression pattern and summarize overall dynamic activity in ordered condition experiments. For each gene, Trendy finds the optimal segmented regression model and provides the location and direction of dynamic changes in expression. We demonstrate the utility of Trendy to provide biologically relevant results on both microarray and RNA-sequencing (RNA-seq) datasets.
Trendy is a flexible R package which characterizes gene-specific expression patterns and summarizes changes of global dynamics over ordered conditions. Trendy is freely available on Bioconductor with a full vignette at https://bioconductor.org/packages/release/bioc/html/Trendy.html .
Insulin-like growth factors (IGFs), specifically IGF1 and IGF2, promote glucose metabolism, with their availability regulated by IGF-binding proteins (IGFBPs). We hypothesized that IGF1 and IGF2 ...levels, or their bioavailability, are reduced during type 1 diabetes development. Total serum IGF1, IGF2, and IGFBP1-7 levels were measured in an age-matched, cross-sectional cohort at varying stages of progression to type 1 diabetes. IGF1 and IGF2 levels were significantly lower in autoantibody (AAb)
compared with AAb
relatives of subjects with type 1 diabetes. Most high-affinity IGFBPs were unchanged in individuals with pre-type 1 diabetes, suggesting that total IGF levels may reflect bioactivity. We also measured serum IGFs from a cohort of fasted subjects with type 1 diabetes. IGF1 levels significantly decreased with disease duration, in parallel with declining β-cell function. Additionally, plasma IGF levels were assessed in an AAb
cohort monthly for a year. IGF1 and IGF2 showed longitudinal stability in single AAb
subjects, but IGF1 levels decreased over time in subjects with multiple AAb and those who progressed to type 1 diabetes, particularly postdiagnosis. In sum, IGFs are dysregulated both before and after the clinical diagnosis of type 1 diabetes and may serve as novel biomarkers to improve disease prediction.
Ahmed and colleagues recently described a novel hybrid lymphocyte expressing both a B and T cell receptor, termed double expresser (DE) cells. DE cells in blood of type 1 diabetes (T1D) subjects were ...present at increased numbers and enriched for a public B cell clonotype. Here, we attempted to reproduce these findings. While we could identify DE cells by flow cytometry, we found no association between DE cell frequency and T1D status. We were unable to identify the reported public B cell clone, or any similar clone, in bulk B cells or sorted DE cells from T1D subjects or controls. We also did not observe increased usage of the public clone VH or DH genes in B cells or in sorted DE cells. Taken together, our findings suggest that DE cells and their alleged public clonotype are not enriched in T1D. This Matters Arising paper is in response to Ahmed et al. (2019), published in Cell. See also the response by Ahmed et al. (2021), published in this issue.
Display omitted
•DE cells not increased in blood or lymphoid tissue of T1D subjects compared to controls•Failure to replicate results of Ahmed et al. with identical FACS assay•“Public” clone-x BCR sequence not found in blood or spleen of T1D and controls•CDR3 AA sequences similar to clone-x are not enriched in B cells or sorted DE cells
There does not appear to be increased abundance of dual-expresser TCR+/BCR+ lymphocytes in type I diabetes patients, nor a public clonotype, challenging the results of a previous study.
Both severe SARS-CoV-2 infections and bacterial sepsis exhibit an immunological dyscrasia and propensity for secondary infections. The nature of the immunological dyscrasias for these differing ...etiologies and their time course remain unclear. In this study, thirty hospitalized patients with SARS-CoV-2 infection were compared with ten critically ill patients with bacterial sepsis over 21 days, as well as ten healthy control subjects. Blood was sampled between days 1 and 21 after admission for targeted plasma biomarker analysis, cellular phenotyping, and leukocyte functional analysis
enzyme-linked immunospot assay. We found that circulating inflammatory markers were significantly higher early after bacterial sepsis compared with SARS-CoV-2. Both cohorts exhibited profound immune suppression through 21 days (suppressed HLA-DR expression, reduced mononuclear cell IFN-gamma production), and expanded numbers of myeloid-derived suppressor cells (MDSCs). In addition, MDSC expansion and
production of IFN-gamma and TNF-alpha were resolving over time in bacterial sepsis, whereas in SARS-CoV-2, immunosuppression and inflammation were accelerating. Despite less severe initial physiologic derangement, SARS-CoV-2 patients had similar incidence of secondary infections (23% vs 30%) as bacterial sepsis patients. Finally, COVID patients who developed secondary bacterial infections exhibited profound immunosuppression evident by elevated sPD-L1 and depressed HLA-DR. Although both bacterial sepsis and SARS-CoV-2 are associated with inflammation and immune suppression, their immune dyscrasia temporal patterns and clinical outcomes are different. SARS-CoV-2 patients had less severe early inflammation and organ dysfunction but had persistent inflammation and immunosuppression and suffered worse clinical outcomes, especially when SARS-CoV-2 infection was followed by secondary bacterial infection.
The human T lymphocyte compartment is highly dynamic over the course of a lifetime. Of the many changes, perhaps most notable is the transition from a predominantly naïve T cell state at birth to the ...acquisition of antigen-experienced memory and effector subsets following environmental exposures. These phenotypic changes, including the induction of T cell exhaustion and senescence, have the potential to negatively impact efficacy of adoptive T cell therapies (ACT). When considering ACT with CD4
CD25
CD127
regulatory T cells (Tregs) for the induction of immune tolerance, we previously reported
expanded umbilical cord blood (CB) Tregs remained more naïve, suppressed responder T cells equivalently, and exhibited a more diverse T cell receptor (TCR) repertoire compared to expanded adult peripheral blood (APB) Tregs. Herein, we hypothesized that upon further characterization, we would observe increased lineage heterogeneity and phenotypic diversity in APB Tregs that might negatively impact lineage stability, engraftment capacity, and the potential for Tregs to home to sites of tissue inflammation following ACT. We compared the phenotypic profiles of human Tregs isolated from CB versus the more traditional source, APB. We conducted analysis of fresh and
expanded Treg subsets at both the single cell (scRNA-seq and flow cytometry) and bulk (microarray and cytokine profiling) levels. Single cell transcriptional profiles of pre-expansion APB Tregs highlighted a cluster of cells that showed increased expression of genes associated with effector and pro-inflammatory phenotypes (
,
, and NKG7) with low expression of Treg markers (
and
). CB Tregs were more diverse in TCR repertoire and homogenous in phenotype, and contained fewer effector-like cells in contrast with APB Tregs. Interestingly, expression of canonical Treg markers, such as FOXP3, TIGIT, and IKZF2, were increased in CB CD4
CD127
conventional T cells (Tconv) compared to APB Tconv, post-expansion, implying perinatal T cells may adopt a default regulatory program. Collectively, these data identify surface markers (namely CXCR3) that could be depleted to improve purity and stability of APB Tregs, and support the use of expanded CB Tregs as a potentially optimal ACT modality for the treatment of autoimmune and inflammatory diseases.