Search

BackgroundMany approaches have been developed to overcome technical noise in single cell RNA-sequencing (scRNAseq). As researchers dig deeper into data-looking for rare cell types, subtleties of cell states, and details of gene regulatory networks-there is a growing need for algorithms with controllable accuracy and fewer ad hoc parameters and thresholds. Impeding this goal is the fact that an appropriate null distribution for scRNAseq cannot simply be extracted from data in which ground truth about biological variation is unknown (i.e., usually).ResultsWe approach this problem analytically, assuming that scRNAseq data reflect only cell heterogeneity (what we seek to characterize), transcriptional noise (temporal fluctuations randomly distributed across cells), and sampling error (i.e., Poisson noise). We analyze scRNAseq data without normalization-a step that skews distributions, particularly for sparse data-and calculate p values associated with key statistics. We develop an improved method for selecting features for cell clustering and identifying gene-gene correlations, both positive and negative. Using simulated data, we show that this method, which we call BigSur (Basic Informatics and Gene Statistics from Unnormalized Reads), captures even weak yet significant correlation structures in scRNAseq data. Applying BigSur to data from a clonal human melanoma cell line, we identify thousands of correlations that, when clustered without supervision into gene communities, align with known cellular components and biological processes, and highlight potentially novel cell biological relationships.ConclusionsNew insights into functionally relevant gene regulatory networks can be obtained using a statistically grounded approach to the identification of gene-gene correlations.

Immune checkpoint therapies such as PD-1 blockade have vastly improved the treatment of numerous cancers, including basal cell carcinoma (BCC). However, patients afflicted with pancreatic ductal carcinoma (PDAC), one of the deadliest malignancies, overwhelmingly exhibit negative responses to checkpoint therapy. We sought to combine data analysis and machine learning to differentiate the putative mechanisms of BCC and PDAC non-response. We discover that increased MHC-I expression in malignant cells and suppression of MHC and PD-1/PD-L expression in CD8+ T cells is associated with nonresponse to treatment. Furthermore, we leverage machine learning to predict response to PD-1 blockade on a cellular level. We confirm divergent resistance mechanisms between BCC, PDAC, and melanoma and highlight the potential for rapid and affordable testing of gene expression in BCC patients to accurately predict response to checkpoint therapies. Our findings present an optimistic outlook for the use of quantitative cross-cancer analyses in characterizing immune responses and predicting immunotherapy outcomes.

The effects of adjuvants for increasing the immunogenicity of influenza vaccines are well known. However, the effect of adjuvants on increasing the breadth of cross-reactivity is less well understood. In this study we have performed a systematic screen of different toll-like receptor (TLR) agonists, with and without a squalene-in-water emulsion on the immunogenicity of a recombinant trimerized hemagglutinin (HA) vaccine in mice after single-dose administration. Antibody (Ab) cross-reactivity for other variants within and outside the immunizing subtype (homosubtypic and heterosubtypic cross-reactivity, respectively) was assessed using a protein microarray approach. Most adjuvants induced broad IgG profiles, although the response to a combination of CpG, MPLA and AddaVax (termed 'IVAX-1') appeared more quickly and reached a greater magnitude than the other formulations tested. Antigen-specific plasma cell labeling experiments show the components of IVAX-1 are synergistic. This adjuvant preferentially stimulates CD4 T cells to produce Th1>Th2 type (IgG2c>IgG1) antibodies and cytokine responses. Moreover, IVAX-1 induces identical homo- and heterosubtypic IgG and IgA cross-reactivity profiles when administered intranasally. Consistent with these observations, a single-cell transcriptomics analysis demonstrated significant increases in expression of IgG1, IgG2b and IgG2c genes of B cells in H5/IVAX-1 immunized mice relative to naïve mice, as well as significant increases in expression of the IFNγ gene of both CD4 and CD8 T cells. These data support the use of adjuvants for enhancing the breath and durability of antibody responses of influenza virus vaccines.

Characterizing genome-wide binding profiles of transcription factors (TFs) is essential for understanding biological processes. Although techniques have been developed to assess binding profiles within a population of cells, determining them at a single-cell level remains elusive. Here, we report scFAN (single-cell factor analysis network), a deep learning model that predicts genome-wide TF binding profiles in individual cells. scFAN is pretrained on genome-wide bulk assay for transposase-accessible chromatin sequencing (ATAC-seq), DNA sequence, and chromatin immunoprecipitation sequencing (ChIP-seq) data and uses single-cell ATAC-seq to predict TF binding in individual cells. We demonstrate the efficacy of scFAN by both studying sequence motifs enriched within predicted binding peaks and using predicted TFs for discovering cell types. We develop a new metric "TF activity score" to characterize each cell and show that activity scores can reliably capture cell identities. scFAN allows us to discover and study cellular identities and heterogeneity based on chromatin accessibility profiles.

The advent of immune checkpoint therapy for metastatic skin cancer has greatly improved patient survival. However, most skin cancer patients are refractory to checkpoint therapy, and furthermore, the intra-immune cell signaling driving response to checkpoint therapy remains uncharacterized. When comparing the immune transcriptome in the tumor microenvironment of melanoma and basal cell carcinoma (BCC), we found that the presence of memory B cells and macrophages negatively correlate in both cancers when stratifying patients by their response, with memory B cells more present in responders. Moreover, inhibitory immune signaling mostly decreases in melanoma responders and increases in BCC responders. We further explored the relationships between macrophages, B cells and response to checkpoint therapy by developing a stochastic differential equation model which qualitatively agrees with the data analysis. Our model predicts BCC to be more refractory to checkpoint therapy than melanoma and predicts the best qualitative ratio of memory B cells and macrophages for successful treatment.

Adjuvants play a central role in enhancing the immunogenicity of otherwise poorly immunogenic vaccine antigens. Combining adjuvants has the potential to enhance vaccine immunogenicity compared with single adjuvants, although the cellular and molecular mechanisms of combination adjuvants are not well understood. Using the influenza virus hemagglutinin H5 antigen, we define the immunological landscape of combining CpG and MPLA (TLR-9 and TLR-4 agonists, respectively) with a squalene nanoemulsion (AddaVax) using immunologic and transcriptomic profiling. Mice immunized and boosted with recombinant H5 in AddaVax, CpG+MPLA, or AddaVax plus CpG+MPLA (IVAX-1) produced comparable levels of neutralizing antibodies and were equally well protected against the H5N1 challenge. However, after challenge with H5N1 virus, H5/IVAX-1-immunized mice had 100- to 300-fold lower virus lung titers than mice receiving H5 in AddaVax or CpG+MPLA separately. Consistent with enhanced viral clearance, unsupervised expression analysis of draining lymph node cells revealed the combination adjuvant IVAX-1 significantly downregulated immune homeostasis genes, and induced higher numbers of antibody-producing plasmablasts than either AddaVax or CpG+MPLA. IVAX-1 was also more effective after single-dose administration than either AddaVax or CpG+MPLA. These data reveal a novel molecular framework for understanding the mechanisms of combination adjuvants, such as IVAX-1, and highlight their potential for the development of more effective vaccines against respiratory viruses., (© The Author(s) 2025. Published by Oxford University Press on behalf of The American Association of Immunologists.)

The high dimensionality of data in single cell transcriptomics (scRNAseq) requires investigators to choose subsets of genes (feature selection) for downstream analysis (e.g., unsupervised cell clustering). The evaluation of different approaches to feature selection is hampered by the fact that, as we show here, the performance of feature selection methods varies greatly with the task being performed. For routine cell type identification, even randomly chosen features can perform well, but for cell type differences that are subtle, both number of features and selection strategy can matter strongly. Here we present a simple feature selection method grounded in an analytical model that, without resorting to arbitrary thresholds or user-defined parameters, allows for interpretable delineation of both how many and which features to choose, facilitating identification of biologically meaningful rare cell types. We compare this method to default methods in scanpy and Seurat, as well as SCTransform, showing how greater accuracy can often be achieved with surprisingly few, well-chosen features.

Background: Many approaches have been developed to overcome technical noise in single cell RNA-sequencing (scRNAseq). As researchers dig deeper into data-looking for rare cell types, subtleties of cell states, and details of gene regulatory networks-there is a growing need for algorithms with controllable accuracy and fewer ad hoc parameters and thresholds. Impeding this goal is the fact that an appropriate null distribution for scRNAseq cannot simply be extracted from data when ground truth about biological variation is unknown (i.e., usually)., Results: We approach this problem analytically, assuming that scRNAseq data reflect only cell heterogeneity (what we seek to characterize), transcriptional noise (temporal fluctuations randomly distributed across cells), and sampling error (i.e., Poisson noise). We analyze scRNAseq data without normalization-a step that skews distributions, particularly for sparse data-and calculate p -values associated with key statistics. We develop an improved method for selecting features for cell clustering and identifying gene-gene correlations, both positive and negative. Using simulated data, we show that this method, which we call BigSur (Basic Informatics and Gene Statistics from Unnormalized Reads), captures even weak yet significant correlation structures in scRNAseq data. Applying BigSur to data from a clonal human melanoma cell line, we identify thousands of correlations that, when clustered without supervision into gene communities, align with known cellular components and biological processes, and highlight potentially novel cell biological relationships., Conclusions: New insights into functionally relevant gene regulatory networks can be obtained using a statistically grounded approach to the identification of gene-gene correlations., Competing Interests: Competing interests: The authors declare that they have no competing interests