Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for investigating cellular biology at an unprecedented resolution, enabling the characterization of cellular heterogeneity, identification of rare but significant cell types, and exploration of cell–cell communications and interactions. Its broad applications span both basic and clinical research domains. In this comprehensive review, we survey the current landscape of scRNA-seq analysis methods and tools, focusing on count modeling, cell-type annotation, data integration, including spatial transcriptomics, and the inference of cell–cell communication. We review the challenges encountered in scRNA-seq analysis, including issues of sparsity or low expression, reliability of cell annotation, and assumptions in data integration, and discuss the potential impact of suboptimal clustering and differential expression analysis tools on downstream analyses, particularly in identifying cell subpopulations. Finally, we discuss recent advancements and future directions for enhancing scRNA-seq analysis. Specifically, we highlight the development of novel tools for annotating single-cell data, integrating and interpreting multimodal datasets covering transcriptomics, epigenomics, and proteomics, and inferring cellular communication networks. By elucidating the latest progress and innovation, we provide a comprehensive overview of the rapidly advancing field of scRNA-seq analysis.

Significance: Primary template-directed amplification (PTA) is a major improvement in whole genome amplification, which is required to study intratissue cellular evolution. As presented in the manuscript, PTA produces significantly improved and reproducible genome sequencing coverage and variant detection from a single genome of a single cell. Applications of measuring genetic diversity from single cells with PTA presented in the manuscript include examining the acquisition of genetic changes during normal development and aging, measuring the consequences of specific perturbations such as genome editing, and characterizing the evolution of clonal populations during cancer formation.

Read counting and unique molecular identifier (UMI) counting are the principal gene expression quantification schemes used in single-cell RNA-sequencing (scRNA-seq) analysis. By using multiple scRNA-seq datasets, we reveal distinct distribution differences between these schemes and conclude that the negative binomial model is a good approximation for UMI counts, even in heterogeneous populations. We further propose a novel differential expression analysis algorithm based on a negative binomial model with independent dispersions in each group (NBID). Our results show that this properly controls the FDR and achieves better power for UMI counts when compared to other recently developed packages for scRNA-seq analysis.

PURPOSE: Although cure rates for childhood acute lymphoblastic leukemia (ALL) exceed 90%, ALL remains a leading cause of cancer death in children. Half of relapses arise in children initially classified with standard-risk (SR) disease. MATERIALS AND METHODS: To identify genomic determinants of relapse in children with SR ALL, we performed genome and transcriptome sequencing of diagnostic and remission samples of children with SR (n = 1,381) or high-risk B-ALL with favorable cytogenetic features (n = 115) enrolled on Children's Oncology Group trials. We used a case-control study design analyzing 439 patients who relapsed and 1,057 who remained in complete remission for at least 5 years. RESULTS: Genomic subtype was associated with relapse, which occurred in approximately 50% of cases of PAX5 -altered ALL (odds ratio [OR], 3.31 [95% CI, 2.17 to 5.03]; P = 3.18 × 10 ^–8 ). Within high-hyperdiploid ALL, gain of chromosome 10 with disomy of chromosome 7 was associated with favorable outcome (OR, 0.27 [95% CI, 0.17 to 0.42]; P = 8.02 × 10 ^–10 ; St Jude Children's Research Hospital validation cohort: OR, 0.22 [95% CI, 0.05 to 0.80]; P = .009), and disomy of chromosomes 10 and 17 with gain of chromosome 6 was associated with relapse (OR, 7.16 [95% CI, 2.63 to 21.51]; P = 2.19 × 10 ^–5 ; validation cohort: OR, 21.32 [95% CI, 3.62 to 119.30]; P = .0004). Genomic alterations were associated with relapse in a subtype-dependent manner, including alterations of INO80 in ETV6::RUNX1 ALL, IKZF1 , and CREBBP in high-hyperdiploid ALL and FHIT in BCR::ABL1 -like ALL. Genomic alterations were also associated with the presence of minimal residual disease, including NRAS and CREBBP in high-hyperdiploid ALL. CONCLUSION: Genetic subtype, patterns of aneuploidy, and secondary genomic alterations determine risk of relapse in childhood ALL. Comprehensive genomic analysis is required for optimal risk stratification.

Germline BRCA2 loss-of function variants, which can be identified through clinical genetic testing, predispose to several cancers^1,2,3,4,5. However, variants of uncertain significance limit the clinical utility of test results. Thus, there is a need for functional characterization and clinical classification of all BRCA2 variants to facilitate the clinical management of individuals with these variants. Here we analysed all possible single-nucleotide variants from exons 15 to 26 that encode the BRCA2 DNA-binding domain hotspot for pathogenic missense variants. To enable this, we used saturation genome editing CRISPR–Cas9-based knock-in endogenous targeting of human haploid HAP1 cells⁶. The assay was calibrated relative to nonsense and silent variants and was validated using pathogenic and benign standards from ClinVar and results from a homology-directed repair functional assay⁷. Variants (6,959 out of 6,960 evaluated) were assigned to seven categories of pathogenicity based on a VarCall Bayesian model⁸. Single-nucleotide variants that encode loss-of-function missense variants were associated with increased risks of breast cancer and ovarian cancer. The functional assay results were integrated into models from ClinGen, the American College of Medical Genetics and Genomics, and the Association for Molecular Pathology⁹ for clinical classification of BRCA2 variants. Using this approach, 91% were classified as pathogenic or likely pathogenic or as benign or likely benign. These classified variants can be used to improve clinical management of individuals with a BRCA2 variant.

Development of αβ and γδ T cells requires coupling of environmental signals with metabolic and redox regulation by mTORC1.