TCR+/BCR+ dual-expressing cells and their associated public BCR clonotype are not enriched in type 1 diabetes.

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.

Data-driven selection of analysis decisions in single-cell RNA-seq trajectory inference.

Single-cell RNA sequencing (scRNA-seq) experiments have become instrumental in developmental and differentiation studies, enabling the profiling of cells at a single or multiple time-points to uncover subtle variations in expression profiles reflecting underlying biological processes. Benchmarking studies have compared many of the computational methods used to reconstruct cellular dynamics; however, researchers still encounter challenges in their analysis due to uncertainty with respect to selecting the most appropriate methods and parameters. Even among universal data processing steps used by trajectory inference methods such as feature selection and dimension reduction, trajectory methods' performances are highly dataset-specific. To address these challenges, we developed Escort, a novel framework for evaluating a dataset's suitability for trajectory inference and quantifying trajectory properties influenced by analysis decisions. Escort evaluates the suitability of trajectory analysis and the combined effects of processing choices using trajectory-specific metrics. Escort navigates single-cell trajectory analysis through these data-driven assessments, reducing uncertainty and much of the decision burden inherent to trajectory inference analyses. Escort is implemented in an accessible R package and R/Shiny application, providing researchers with the necessary tools to make informed decisions during trajectory analysis and enabling new insights into dynamic biological processes at single-cell resolution.

Enhancing biological signals and detection rates in single-cell RNA-seq experiments with cDNA library equalization.

Considerable effort has been devoted to refining experimental protocols to reduce levels of technical variability and artifacts in single-cell RNA-sequencing data (scRNA-seq). We here present evidence that equalizing the concentration of cDNA libraries prior to pooling, a step not consistently performed in single-cell experiments, improves gene detection rates, enhances biological signals, and reduces technical artifacts in scRNA-seq data. To evaluate the effect of equalization on various protocols, we developed Scaffold, a simulation framework that models each step of an scRNA-seq experiment. Numerical experiments demonstrate that equalization reduces variation in sequencing depth and gene-specific expression variability. We then performed a set of experiments in vitro with and without the equalization step and found that equalization increases the number of genes that are detected in every cell by 17-31%, improves discovery of biologically relevant genes, and reduces nuisance signals associated with cell cycle. Further support is provided in an analysis of publicly available data.

Normalization by distributional resampling of high throughput single-cell RNA-sequencing data.

MOTIVATION: Normalization to remove technical or experimental artifacts is critical in the analysis of single-cell RNA-sequencing experiments, even those for which unique molecular identifiers are available. The majority of methods for normalizing single-cell RNA-sequencing data adjust average expression for library size (LS), allowing the variance and other properties of the gene-specific expression distribution to be non-constant in LS. This often results in reduced power and increased false discoveries in downstream analyses, a problem which is exacerbated by the high proportion of zeros present in most datasets. RESULTS: To address this, we present Dino, a normalization method based on a flexible negative-binomial mixture model of gene expression. As demonstrated in both simulated and case study datasets, by normalizing the entire gene expression distribution, Dino is robust to shallow sequencing, sample heterogeneity and varying zero proportions, leading to improved performance in downstream analyses in a number of settings. AVAILABILITY AND IMPLEMENTATION: The R package, Dino, is available on GitHub at https://github.com/JBrownBiostat/Dino. The Dino package is further archived and freely available on Zenodo at https://doi.org/10.5281/zenodo.4897558. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Methylscaper: an R/Shiny app for joint visualization of DNA methylation and nucleosome occupancy in single-molecule and single-cell data.

SUMMARY: Differential DNA methylation and chromatin accessibility are associated with disease development, particularly cancer. Methods that allow profiling of these epigenetic mechanisms in the same reaction and at the single-molecule or single-cell level continue to emerge. However, a challenge lies in jointly visualizing and analyzing the heterogeneous nature of the data and extracting regulatory insight. Here, we present methylscaper, a visualization framework for simultaneous analysis of DNA methylation and chromatin accessibility landscapes. Methylscaper implements a weighted principal component analysis that orders DNA molecules, each providing a record of the chromatin state of one epiallele, and reveals patterns of nucleosome positioning, transcription factor occupancy, and DNA methylation. We demonstrate methylscaper's utility on a long-read, single-molecule methyltransferase accessibility protocol for individual templates (MAPit-BGS) dataset and a single-cell nucleosome, methylation, and transcription sequencing (scNMT-seq) dataset. In comparison to other procedures, methylscaper is able to readily identify chromatin features that are biologically relevant to transcriptional status while scaling to larger datasets. AVAILABILITY AND IMPLEMENTATION: Methylscaper, is implemented in R (version > 4.1) and available on Bioconductor: https://bioconductor.org/packages/methylscaper/, GitHub: https://github.com/rhondabacher/methylscaper/, and Web: https://methylscaper.com. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

SCnorm: robust normalization of single-cell RNA-seq data.

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.

Single-cell RNA sequencing reveals intrinsic and extrinsic regulatory heterogeneity in yeast responding to stress.

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.

Automated minute scale RNA-seq of pluripotent stem cell differentiation reveals early divergence of human and mouse gene expression kinetics.

Pluripotent stem cells retain the developmental timing of their species of origin in vitro, an observation that suggests the existence of a cell-intrinsic developmental clock, yet the nature and machinery of the clock remain a mystery. We hypothesize that one possible component may lie in species-specific differences in the kinetics of transcriptional responses to differentiation signals. Using a liquid-handling robot, mouse and human pluripotent stem cells were exposed to identical neural differentiation conditions and sampled for RNA-sequencing at high frequency, every 4 or 10 minutes, for the first 10 hours of differentiation to test for differences in transcriptomic response rates. The majority of initial transcriptional responses occurred within a rapid window in the first minutes of differentiation for both human and mouse stem cells. Despite similarly early onsets of gene expression changes, we observed shortened and condensed gene expression patterns in mouse pluripotent stem cells compared to protracted trends in human pluripotent stem cells. Moreover, the speed at which individual genes were upregulated, as measured by the slopes of gene expression changes over time, was significantly faster in mouse compared to human cells. These results suggest that downstream transcriptomic response kinetics to signaling cues are faster in mouse versus human cells, and may offer a partial account for the vast differences in developmental rates across species.

Spatial patterns of gene expression are unveiled in the chick primitive streak by ordering single-cell transcriptomes.

During vertebrate development, progenitor cells give rise to tissues and organs through a complex choreography that commences at gastrulation. A hallmark event of gastrulation is the formation of the primitive streak, a linear assembly of cells along the anterior-posterior (AP) axis of the developing organism. To examine the primitive streak at a single-cell resolution, we measured the transcriptomes of individual chick cells from the streak or the surrounding tissue (the rest of the area pellucida) in Hamburger-Hamilton stage 4 embryos. The single-cell transcriptomes were then ordered by the statistical method Wave-Crest to deduce both the relative position along the AP axis and the prospective lineage of single cells. The ordered transcriptomes reveal intricate patterns of gene expression along the primitive streak.

Trendy: segmented regression analysis of expression dynamics in high-throughput ordered profiling experiments.

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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 .

Species-specific developmental timing is maintained by pluripotent stem cells ex utero.

How species-specific developmental timing is controlled is largely unknown. By following human embryonic stem (ES) cell and mouse epiblast stem (EpiS) cell differentiation through detailed RNA-sequencing time courses, here we show that pluripotent stem cells closely retain in vivo species-specific developmental timing in vitro. In identical neural differentiation conditions in vitro, gene expression profiles are accelerated in mouse EpiS cells compared to human ES cells with relative rates of differentiation closely reflecting the rates of progression through the Carnegie stages in utero. Dynamic Time Warping analysis identified 3389 genes that were regulated more quickly in mouse EpiS cells and identified none that were regulated more quickly in human ES cells. Interestingly, we also find that human ES cells differentiated in teratomas maintain the same rate of differentiation observed in vitro in spite of being grown in a mouse host. These results suggest the existence of a cell autonomous, species-specific developmental clock that pluripotent stem cells maintain even out of context of an intact embryo.

TRANSCRIPTOMIC DIFFERENCES IN PERIPHERAL MONOCYTE POPULATIONS IN SEPTIC PATIENTS BASED ON OUTCOME.

ABSTRACT: Postsepsis early mortality is being replaced by survivors who experience either a rapid recovery and favorable hospital discharge or the development of chronic critical illness with suboptimal outcomes. The underlying immunological response that determines these clinical trajectories remains poorly defined at the transcriptomic level. As classical and nonclassical monocytes are key leukocytes in both the innate and adaptive immune systems, we sought to delineate the transcriptomic response of these cell types. Using single-cell RNA sequencing and pathway analyses, we identified gene expression patterns between these two groups that are consistent with differences in TNF-α production based on clinical outcome. This may provide therapeutic targets for those at risk for chronic critical illness in order to improve their phenotype/endotype, morbidity, and long-term mortality.

The post-septic peripheral myeloid compartment reveals unexpected diversity in myeloid-derived suppressor cells.

Introduction: Sepsis engenders distinct host immunologic changes that include the expansion of myeloid-derived suppressor cells (MDSCs). These cells play a physiologic role in tempering acute inflammatory responses but can persist in patients who develop chronic critical illness. Methods: Cellular Indexing of Transcriptomes and Epitopes by Sequencing and transcriptomic analysis are used to describe MDSC subpopulations based on differential gene expression, RNA velocities, and biologic process clustering. Results: We identify a unique lineage and differentiation pathway for MDSCs after sepsis and describe a novel MDSC subpopulation. Additionally, we report that the heterogeneous response of the myeloid compartment of blood to sepsis is dependent on clinical outcome. Discussion: The origins and lineage of these MDSC subpopulations were previously assumed to be discrete and unidirectional; however, these cells exhibit a dynamic phenotype with considerable plasticity.

Analysis of Single-Cell RNA-seq Data.

As single-cell RNA sequencing experiments continue to advance scientific discoveries across biological disciplines, an increasing number of analysis tools and workflows for analyzing the data have been developed. In this chapter, we describe a standard workflow and elaborate on relevant data analysis tools for analyzing single-cell RNA sequencing data. We provide recommendations for the appropriate use of commonly used methods, with code examples and analysis interpretations.

Interpretable trajectory inference with single-cell Linear Adaptive Negative-binomial Expression (scLANE) testing.

The rapid proliferation of trajectory inference methods for single-cell RNA-seq data has allowed researchers to investigate complex biological processes by examining underlying gene expression dynamics. After estimating a latent cell ordering, statistical models are used to determine which genes exhibit changes in expression that are significantly associated with progression through the biological trajectory. While a few techniques for performing trajectory differential expression exist, most rely on the flexibility of generalized additive models in order to account for the inherent nonlinearity of changes in gene expression. As such, the results can be difficult to interpret, and biological conclusions often rest on subjective visual inspections of the most dynamic genes. To address this challenge, we propose scLANE testing, which is built around an interpretable generalized linear model and handles nonlinearity with basis splines chosen empirically for each gene. In addition, extensions to estimating equations and mixed models allow for reliable trajectory testing under complex experimental designs. After validating the accuracy of scLANE under several different simulation scenarios, we apply it to a set of diverse biological datasets and display its ability to provide novel biological information when used downstream of both pseudotime and RNA velocity estimation methods.

Data-driven selection of analysis decisions in single-cell RNA-seq trajectory inference.

Single-cell RNA sequencing (scRNA-seq) experiments have become instrumental in developmental and differentiation studies, enabling the profiling of cells at a single or multiple time-points to uncover subtle variations in expression profiles reflecting underlying biological processes. Benchmarking studies have compared many of the computational methods used to reconstruct cellular dynamics, however researchers still encounter challenges in their analysis due to uncertainties in selecting the most appropriate methods and parameters. Even among universal data processing steps used by trajectory inference methods such as feature selection and dimension reduction, trajectory methods' performances are highly dataset-specific. To address these challenges, we developed Escort, a framework for evaluating a dataset's suitability for trajectory inference and quantifying trajectory properties influenced by analysis decisions. Escort navigates single-cell trajectory analysis through data-driven assessments, reducing uncertainty and much of the decision burden associated with trajectory inference. Escort is implemented in an accessible R package and R/Shiny application, providing researchers with the necessary tools to make informed decisions during trajectory analysis and enabling new insights into dynamic biological processes at single-cell resolution.

Biomaterials-based nanoparticles conjugated to regulatory T cells provide a modular system for localized delivery of pharmacotherapeutic agents.

Type 1 diabetes (T1D) presents with two therapeutic challenges: the need to correct underlying autoimmunity and restore ß-cell mass. We harnessed the unique capacity of regulatory T cells (Tregs) and the T cell receptor (TCR) to direct tolerance induction along with tissue-localized delivery of therapeutic agents to restore endogenous ß-cell function. Specifically, we designed a combinatorial therapy involving biomaterials-based poly(lactic-co-glycolic acid) nanoparticles co-loaded with the Treg growth factor, IL-2, and the ß-cell regenerative agent, harmine (a tyrosine-regulated kinase 1A [DYRK1A] inhibitor), conjugated to the surface of Tregs. We observed continuous elution of IL-2 and harmine from nanoparticles for at least 7 days in vitro. When conjugated to primary human Tregs, IL-2 nanoparticles provided sufficient IL-2 receptor signaling to support STAT5 phosphorylation for sustained phenotypic stability and viability in culture. Inclusion of poly-L-lysine (PLL) during nanoparticle-cell coupling dramatically increased conjugation efficiency, providing sufficient IL-2 to support in vitro proliferation of IL-2-dependent CTLL-2 cells and primary murine Tregs. In 12-week-old female non-obese diabetic mice, adoptive transfer of IL-2/harmine nanoparticle-conjugated NOD.BDC2.5 Tregs, which express an islet antigen-specific TCR, significantly prevented diabetes demonstrating preserved in vivo viability. These data provide the preclinical basis to develop a biomaterials-optimized cellular therapy to restore immune tolerance and promote ß-cell proliferation in T1D through receptor-targeted drug delivery within pancreatic islets.

Responders to low-dose ATG induce CD4+ T cell exhaustion in type 1 diabetes.

BACKGROUNDLow-dose anti-thymocyte globulin (ATG) transiently preserves C-peptide and lowers HbA1c in individuals with recent-onset type 1 diabetes (T1D); however, the mechanisms of action and features of the response remain unclear. Here, we characterized the post hoc immunological outcomes of ATG administration and their potential use as biomarkers of metabolic response to therapy (i.e., improved preservation of endogenous insulin production).METHODSWe assessed gene and protein expression, targeted gene methylation, and cytokine concentrations in peripheral blood following treatment with ATG (n = 29), ATG plus granulocyte colony-stimulating factor (ATG/G-CSF, n = 28), or placebo (n = 31).RESULTSTreatment with low-dose ATG preserved regulatory T cells (Tregs), as measured by stable methylation of FOXP3 Treg-specific demethylation region (TSDR) and increased proportions of CD4+FOXP3+ Tregs (P < 0.001) identified by flow cytometry. While treatment effects were consistent across participants, not all maintained C-peptide. Responders exhibited a transient rise in IL-6, IP-10, and TNF-α (P < 0.05 for all) 2 weeks after treatment and a durable CD4+ exhaustion phenotype (increased PD-1+KLRG1+CD57- on CD4+ T cells [P = 0.011] and PD1+CD4+ Temra MFI [P < 0.001] at 12 weeks, following ATG and ATG/G-CSF, respectively). ATG nonresponders displayed higher proportions of senescent T cells (at baseline and after treatment) and increased methylation of EOMES (i.e., less expression of this exhaustion marker).CONCLUSIONAltogether in these exploratory analyses, Th1 inflammation-associated serum and CD4+ exhaustion transcript and cellular phenotyping profiles may be useful for identifying signatures of clinical response to ATG in T1D.TRIAL REGISTRATIONClinicalTrials.gov NCT02215200.FUNDINGThe Leona M. and Harry B. Helmsley Charitable Trust (2019PG-T1D011), the NIH (R01 DK106191 Supplement, K08 DK128628), NIH TrialNet (U01 DK085461), and the NIH NIAID (P01 AI042288).

Human immune phenotyping reveals accelerated aging in type 1 diabetes.

The proportions and phenotypes of immune cell subsets in peripheral blood undergo continual and dramatic remodeling throughout the human life span, which complicates efforts to identify disease-associated immune signatures in type 1 diabetes (T1D). We conducted cross-sectional flow cytometric immune profiling on peripheral blood from 826 individuals (stage 3 T1D, their first-degree relatives, those with ≥2 islet autoantibodies, and autoantibody-negative unaffected controls). We constructed an immune age predictive model in unaffected participants and observed accelerated immune aging in T1D. We used generalized additive models for location, shape, and scale to obtain age-corrected data for flow cytometry and complete blood count readouts, which can be visualized in our interactive portal (ImmScape); 46 parameters were significantly associated with age only, 25 with T1D only, and 23 with both age and T1D. Phenotypes associated with accelerated immunological aging in T1D included increased CXCR3+ and programmed cell death 1-positive (PD-1+) frequencies in naive and memory T cell subsets, despite reduced PD-1 expression levels on memory T cells. Phenotypes associated with T1D after age correction were predictive of T1D status. Our findings demonstrate advanced immune aging in T1D and highlight disease-associated phenotypes for biomarker monitoring and therapeutic interventions.

Data-driven assessment of dimension reduction quality for single-cell omics data.

Dimension reduction (DR) techniques have become synonymous with single-cell omics data due to their ability to generate attractive visualizations and enable analyses of high-dimensional data. In this issue of Patterns, Johnsona et al. develop a statistical approach to assist in selecting high-quality reduced representations to improve analyses and biological interpretations.