Quantitative assessment of single-cell RNA-sequencing methods

Patterns of Associations with Epidemiologic Factors by High-Grade Serous Ovarian Cancer Gene Expression Subtypes

BACKGROUND

Ovarian high-grade serous carcinomas (HGSC) comprise four distinct molecular subtypes based on mRNA expression patterns, with differential survival. Understanding risk factor associations is important to elucidate the etiology of HGSC. We investigated associations between different epidemiologic risk factors and HGSC molecular subtypes.

METHODS

We pooled data from 11 case-control studies with epidemiologic and tumor gene expression data from custom NanoString CodeSets developed through a collaboration within the Ovarian Tumor Tissue Analysis consortium. The PrOTYPE-validated NanoString-based 55-gene classifier was used to assign HGSC gene expression subtypes. We examined associations between epidemiologic factors and HGSC subtypes in 2,070 cases and 16,633 controls using multivariable-adjusted polytomous regression models.

RESULTS

Among the 2,070 HGSC cases, 556 (27%) were classified as C1.MES, 340 (16%) as C5.PRO, 538 (26%) as C2.IMM, and 636 (31%) as C4.DIF. The key factors, including oral contraceptive use, parity, breastfeeding, and family history of ovarian cancer, were similarly associated with all subtypes. Heterogeneity was observed for several factors. Former smoking [OR = 1.25; 95% confidence interval (CI) = 1.03, 1.51] and genital powder use (OR = 1.42; 95% CI = 1.08, 1.86) were uniquely associated with C2.IMM. History of endometriosis was associated with C5.PRO (OR = 1.46; 95% CI = 0.98, 2.16) and C4.DIF (OR = 1.27; 95% CI = 0.94, 1.71) only. Family history of breast cancer (OR = 1.44; 95% CI = 1.16, 1.78) and current smoking (OR = 1.40; 95% CI = 1.11, 1.76) were associated with C4.DIF only.

CONCLUSIONS

This study observed heterogeneous associations of epidemiologic and modifiable factors with HGSC molecular subtypes.

IMPACT

The different patterns of associations may provide key information about the etiology of the four subtypes.

Heterogeneity of the tumor immune cell microenvironment revealed by single-cell sequencing in head and neck cancer

Head and neck cancer (HNC) is the sixth most common disease in the world. The recurrence rate of patients is relatively high, and the heterogeneity of tumor immune microenvironment (TIME) cells may be an important reason for this. Single-cell sequencing (SCS) is currently the most promising and mature application in cancer research. It can identify unique genes expressed in cells and study tumor heterogeneity. According to current research, the heterogeneity of immune cells has become an important factor affecting the occurrence and development of HNC. SCSs can provide effective therapeutic targets and prognostic factors for HNC patients through analyses of gene expression levels and cell heterogeneity. Therefore, this study analyzes the basic theory of HNC and the development of SCS technology, elaborating on the application of SCS technology in HNC and its potential value in identifying HNC therapeutic targets and biomarkers.

A systematic review of single-cell RNA sequencing applications and innovations

Bulk RNA sequencing is one type of RNA sequencing technique, as well as targeted RNA sequencing and whole transcriptome sequencing. It provides valuable insights into gene expression in specific cell populations or regions. However, these methods often miss the diversity of cells within complex tissues. This restriction is overcome by single-cell RNA sequencing, which records gene expression at the single-cell level. It offers a detailed picture of the diversity of cells. It is essential to study glucose homeostasis. It offers thorough explanations of cellular variation. Networks and Governance Dynamics The use of scRNA-seq in islet cells is reviewed in this study, along with sample preparation, sequencing, and computational analysis. It highlights advances in understanding cell types. Gene activity and cell interactions. Along with the challenges and limitations of scRNA-seq, this review highlights the importance of scRNA-seq in understanding complex biological processes and diseases. It is an essential resource for future research and method development in this field, which will help to build personalized treatment.

How to view the female reproductive tract through single-cell looking glasses

Single-cell technologies have emerged as an unprecedented tool for biologists and clinicians, allowing them to assess organs and tissues at the level of individual cells. In the field of women's reproductive biology, single-cell studies have provided insights into the cellular and molecular processes that regulate reproductive and obstetrical functions in health and disease. The knowledge that these studies generate is helping clinicians to improve the understanding and diagnosis of infertility related issues or pregnancy complications and to find new avenues for their treatment. However, navigating the expansive landscape of this type of transcriptomic data analysis represents a pivotal challenge in current research. Single cell RNA sequencing involves isolating cells into droplets, reverse transcribing RNA to generate complementary DNA, with each droplet content uniquely labeled by a barcode. Upon sequencing the complementary DNAs, the barcodes enable the reassignment of sequencing reads to individual droplets, facilitating the reconstruction of the cellular landscape of the sample obtained from a tissue or organ and beyond. Researchers, equipped with the metaphorical 'single-cell glasses,' must adequately choose from a plethora of strategies to dissect and interpret cellular information. Sophisticated algorithms and the decision-making process are often underestimated, resulting in artefactual or cumbersome interpreted results. Computational biologists apply and innovate computational tools designed to process, model, and interpret expansive datasets. The ramifications of their work extend far beyond the realm of data processing; they give shape to the outcome of analyses, playing a pivotal role in drawing meaningful conclusions from the wealth of information garnered. In this review, we describe the wide variety of approaches and analytical steps available with enough detail to gain a concise picture of what a complete examination of a single-cell dataset would be. We commence with a discussion on key points in experimental design, highlighting crucial questions one should consider. Following this, we delve into the various preprocessing and quality control steps essential for any single-cell dataset. The subsequent section offers a detailed guide on constructing a single-cell atlas, exploring nuances such as differential characteristics in visualization and clustering techniques, as well as strategies for assigning identity to cell populations through gene marker annotations. Moving beyond the creation of an atlas, we explore methods for investigating pathological conditions. This involves conducting cell population comparison tests between conditions and analyzing specific cell-to-cell communications and cellular differentiation trajectories in both health and disease scenarios. This work aims to furnish a newcomer researcher and/or clinician with essential guidelines to embark on a single-cell adventure without succumbing to common pitfalls. By bridging the gap between theory and practice, it facilitates the translation of single-cell technologies into clinically relevant applications. Throughout the manuscript, practical examples of its usage in women's reproductive health studies are provided. Various sections delve into specific clinical scenarios, demonstrating how these guidelines can be instrumental in unraveling the molecular landscapes of diseases and physiological processes related to women's reproduction.

Cell cycle checkpoint factor p15Ink4b is a novel regulator of osteoclast differentiation

Osteoclasts are specialized cells essential for bone resorption, a crucial process in bone remodeling, and dysregulation of osteoclastogenesis can lead to pathological bone loss such as osteoporosis and rheumatoid arthritis. Therefore, understanding the precise mechanisms governing osteoclast differentiation is crucial for developing effective therapies for skeletal diseases. In osteoclastogenesis, as well as other differentiated cells, it is well understood that cell cycle arrest is essential for terminal differentiation and is tightly regulated by CDK inhibitors such as Cip/Kip family and Ink4 family protein. In this manuscript, we identified p15Ink4b, a member of the Ink4 family, as a novel regulator of osteoclastogenesis by comprehensive single-cell RNA sequence data reanalyzing. Furthermore, histological analysis and in vitro osteoclast differentiation assay revealed that p15Ink4b functionally regulates osteoclastogenesis. Our findings may not only provide insights into the molecular mechanisms of osteoclast differentiation but also underscore the potential of harnessing cell cycle mechanisms to develop novel therapeutic strategies for bone diseases.

Nano-scale dangers: Unravelling the impact of nanoplastics on human trophoblast invasion

Utilizing HTR-8/SVneo cells for in vitro modeling of human trophoblast invasion, we examined how different concentrations of 40 nm and 200 nm carboxylated polystyrene particles affect early-pregnancy trophoblast phenotype and function. We focused on migration and invasion, as critical processes in placental development. Our findings revealed disruptions in extravillous trophoblast mesenchymal phenotype and invasive behavior, following acute exposure to a higher concentration of the smaller sized particles. Specifically, differential uptake of the particles by trophoblast cells was observed, as well as cytotoxicity and concentration-dependent DNA damage after 72 h of exposure. In addition, a 24 h exposure to 100 μg/ml of 40 nm particles correlated with downregulated protein expression of α5 and α1 integrin subunits, N-cadherin, matrix metalloproteinase-2 and macrophage migration inhibitory factor, alongside upregulated protein expression of the epithelial marker E-cadherin. These changes likely contributed to the diminished migration of HTR-8/SVneo cells and the invasive potential of HTR-8/SVneo spheroids. Understanding these interactions is paramount for assessing the broader implications of nanoplastics on reproductive outcomes and maternal-fetal well-being and informing public health measures.

Bacterial single-cell RNA sequencing captures biofilm transcriptional heterogeneity and differential responses to immune pressure

Biofilm formation is an important mechanism of survival and persistence for many bacterial pathogens. These multicellular communities contain subpopulations of cells that display metabolic and transcriptional diversity along with recalcitrance to antibiotics and host immune defenses. Here, we present an optimized bacterial single-cell RNA sequencing method, BaSSSh-seq, to study Staphylococcus aureus diversity during biofilm growth and transcriptional adaptations following immune cell exposure. BaSSSh-seq captures extensive transcriptional heterogeneity during biofilm compared to planktonic growth. We quantify and visualize transcriptional regulatory networks across heterogeneous biofilm subpopulations and identify gene sets that are associated with a trajectory from planktonic to biofilm growth. BaSSSh-seq also detects alterations in biofilm metabolism, stress response, and virulence induced by distinct immune cell populations. This work facilitates the exploration of biofilm dynamics at single-cell resolution, unlocking the potential for identifying biofilm adaptations to environmental signals and immune pressure.

Investigation of human aging at the single-cell level

Human aging is characterized by a gradual decline in physiological functions and an increased susceptibility to various diseases. The complex mechanisms underlying human aging are still not fully elucidated. Single-cell sequencing (SCS) technologies have revolutionized aging research by providing unprecedented resolution and detailed insights into cellular diversity and dynamics. In this review, we discuss the application of various SCS technologies in human aging research, encompassing single-cell, genomics, transcriptomics, epigenomics, and proteomics. We also discuss the combination of multiple omics layers within single cells and the integration of SCS technologies with advanced methodologies like spatial transcriptomics and mass spectrometry. These approaches have been essential in identifying aging biomarkers, elucidating signaling pathways associated with aging, discovering novel aging cell subpopulations, uncovering tissue-specific aging characteristics, and investigating aging-related diseases. Furthermore, we provide an overview of aging-related databases that offer valuable resources for enhancing our understanding of the human aging process.

Shrinkage estimation of gene interaction networks in single-cell RNA sequencing data

BACKGROUND

Gene interaction networks are graphs in which nodes represent genes and edges represent functional interactions between them. These interactions can be at multiple levels, for instance, gene regulation, protein-protein interaction, or metabolic pathways. To analyse gene interaction networks at a large scale, gene co-expression network analysis is often applied on high-throughput gene expression data such as RNA sequencing data. With the advance in sequencing technology, expression of genes can be measured in individual cells. Single-cell RNA sequencing (scRNAseq) provides insights of cellular development, differentiation and characteristics at the transcriptomic level. High sparsity and high-dimensional data structures pose challenges in scRNAseq data analysis.

RESULTS

In this study, a sparse inverse covariance matrix estimation framework for scRNAseq data is developed to capture direct functional interactions between genes. Comparative analyses highlight high performance and fast computation of Stein-type shrinkage in high-dimensional data using simulated scRNAseq data. Data transformation approaches also show improvement in performance of shrinkage methods in non-Gaussian distributed data. Zero-inflated modelling of scRNAseq data based on a negative binomial distribution enhances shrinkage performance in zero-inflated data without interference on non zero-inflated count data.

CONCLUSION

The proposed framework broadens application of graphical model in scRNAseq analysis with flexibility in sparsity of count data resulting from dropout events, high performance, and fast computational time. Implementation of the framework is in a reproducible Snakemake workflow https://github.com/calathea24/ZINBGraphicalModel and R package ZINBStein https://github.com/calathea24/ZINBStein .

SpaDiT: diffusion transformer for spatial gene expression prediction using scRNA-seq

The rapid development of spatially resolved transcriptomics (SRT) technologies has provided unprecedented opportunities for exploring the structure of specific organs or tissues. However, these techniques (such as image-based SRT) can achieve single-cell resolution, but can only capture the expression levels of tens to hundreds of genes. Such spatial transcriptomics (ST) data, carrying a large number of undetected genes, have limited its application value. To address the challenge, we develop SpaDiT, a deep learning framework for spatial reconstruction and gene expression prediction using scRNA-seq data. SpaDiT employs scRNA-seq data as an a priori condition and utilizes shared genes between ST and scRNA-seq data as latent representations to construct inputs, thereby facilitating the accurate prediction of gene expression in ST data. SpaDiT enhances the accuracy of spatial gene expression predictions over a variety of spatial transcriptomics datasets. We have demonstrated the effectiveness of SpaDiT by conducting extensive experiments on both seq-based and image-based ST data. We compared SpaDiT with eight highly effective baseline methods and found that our proposed method achieved an 8%-12% improvement in performance across multiple metrics. Source code and all datasets used in this paper are available at https://github.com/wenwenmin/SpaDiT and https://zenodo.org/records/12792074.

Craters on the melanoma surface facilitate tumor-immune interactions and demonstrate pathologic response to checkpoint blockade in humans

Immunotherapy leads to cancer eradication despite the tumor's immunosuppressive environment. Here, we used extended long-term in-vivo imaging and high-resolution spatial transcriptomics of endogenous melanoma in zebrafish, and multiplex imaging of human melanoma, to identify domains that facilitate immune response during immunotherapy. We identified crater-shaped pockets at the margins of zebrafish and human melanoma, rich with beta-2 microglobulin (B2M) and antigen recognition molecules. The craters harbor the highest density of CD8+ T cells in the tumor. In zebrafish, CD8+ T cells formed prolonged interactions with melanoma cells within craters, characteristic of antigen recognition. Following immunostimulatory treatment, the craters enlarged and became the major site of activated CD8+ T cell accumulation and tumor killing that was B2M dependent. In humans, craters predicted immune response to ICB therapy, showing response better than high T cell infiltration. This marks craters as potential new diagnostic tool for immunotherapy success and targets to enhance ICB response.

The Domestication of Wild Boar Could Result in a Relaxed Selection for Maintaining Olfactory Capacity

Domesticated animals are artificially selected to exhibit desirable traits, however not all traits of domesticated animals are the result of deliberate selection. Loss of olfactory capacity in the domesticated pig (Sus scrofa domesticus) is one example. We used whole transcriptome analysis (RNA-Seq) to compare patterns of gene expression in the olfactory mucosa of the pig and two subspecies of wild boar (Sus scrofa), and investigate candidate genes that could be responsible for the loss of olfactory capacity. We identified hundreds of genes with reductions in transcript abundance in pig relative to wild boar as well as differences between the two subspecies of wild boar. These differences were detected mainly in genes involved in the formation and motility of villi, cilia and microtubules, functions associated with olfaction. In addition, differences were found in the abundances of transcripts of genes related to immune defenses, with the highest levels in continental wild boar subspecies. Overall, the loss of olfactory capacity in pigs appears to have been accompanied by reductions in the expression of candidate genes for olfaction. These changes could have resulted from unintentional selection for reduced olfactory capacity, relaxed selection for maintaining olfactory capacity, pleiotropic effects of genes under selection, or other non-selective processes. Our findings could be a cornerstone for future researches on wild boars, pigs, feral populations, and their evolutionary trajectories, aimed to provide tools to better calibrate species management as well as guidelines for breeders.

COMSE: analysis of single-cell RNA-seq data using community detection-based feature selection

BACKGROUND

Single-cell RNA sequencing enables studying cells individually, yet high gene dimensions and low cell numbers challenge analysis. And only a subset of the genes detected are involved in the biological processes underlying cell-type specific functions.

RESULT

In this study, we present COMSE, an unsupervised feature selection framework using community detection to capture informative genes from scRNA-seq data. COMSE identified homogenous cell substates with high resolution, as demonstrated by distinguishing different cell cycle stages. Evaluations based on real and simulated scRNA-seq datasets showed COMSE outperformed methods even with high dropout rates in cell clustering assignment. We also demonstrate that by identifying communities of genes associated with batch effects, COMSE parses signals reflecting biological difference from noise arising due to differences in sequencing protocols, thereby enabling integrated analysis of scRNA-seq datasets of different sources.

CONCLUSIONS

COMSE provides an efficient unsupervised framework that selects highly informative genes in scRNA-seq data improving cell sub-states identification and cell clustering. It identifies gene subsets that reveal biological and technical heterogeneity, supporting applications like batch effect correction and pathway analysis. It also provides robust results for bulk RNA-seq data analysis.

Bacterial single-cell RNA sequencing captures biofilm transcriptional heterogeneity and differential responses to immune pressure

Biofilm formation is an important mechanism of survival and persistence for many bacterial pathogens. These multicellular communities contain subpopulations of cells that display vast metabolic and transcriptional diversity along with high recalcitrance to antibiotics and host immune defenses. Investigating the complex heterogeneity within biofilm has been hindered by the lack of a sensitive and high-throughput method to assess stochastic transcriptional activity and regulation between bacterial subpopulations, which requires single-cell resolution. We have developed an optimized bacterial single-cell RNA sequencing method, BaSSSh-seq, to study Staphylococcus aureus diversity during biofilm growth and transcriptional adaptations following immune cell exposure. We validated the ability of BaSSSh-seq to capture extensive transcriptional heterogeneity during biofilm compared to planktonic growth. Application of new computational tools revealed transcriptional regulatory networks across the heterogeneous biofilm subpopulations and identification of gene sets that were associated with a trajectory from planktonic to biofilm growth. BaSSSh-seq also detected alterations in biofilm metabolism, stress response, and virulence that were tailored to distinct immune cell populations. This work provides an innovative platform to explore biofilm dynamics at single-cell resolution, unlocking the potential for identifying biofilm adaptations to environmental signals and immune pressure.

GLDADec: marker-gene guided LDA modeling for bulk gene expression deconvolution

Inferring cell type proportions from bulk transcriptome data is crucial in immunology and oncology. Here, we introduce guided LDA deconvolution (GLDADec), a bulk deconvolution method that guides topics using cell type-specific marker gene names to estimate topic distributions for each sample. Through benchmarking using blood-derived datasets, we demonstrate its high estimation performance and robustness. Moreover, we apply GLDADec to heterogeneous tissue bulk data and perform comprehensive cell type analysis in a data-driven manner. We show that GLDADec outperforms existing methods in estimation performance and evaluate its biological interpretability by examining enrichment of biological processes for topics. Finally, we apply GLDADec to The Cancer Genome Atlas tumor samples, enabling subtype stratification and survival analysis based on estimated cell type proportions, thus proving its practical utility in clinical settings. This approach, utilizing marker gene names as partial prior information, can be applied to various scenarios for bulk data deconvolution. GLDADec is available as an open-source Python package at https://github.com/mizuno-group/GLDADec.

Native N-glycome profiling of single cells and ng-level blood isolates using label-free capillary electrophoresis-mass spectrometry

The development of reliable single-cell dispensers and substantial sensitivity improvement in mass spectrometry made proteomic profiling of individual cells achievable. Yet, there are no established methods for single-cell glycome analysis due to the inability to amplify glycans and sample losses associated with sample processing and glycan labeling. In this work, we present an integrated platform coupling online in-capillary sample processing with high-sensitivity label-free capillary electrophoresis-mass spectrometry for N-glycan profiling of single mammalian cells. Direct and unbiased quantitative characterization of single-cell surface N-glycomes are demonstrated for HeLa and U87 cells, with the detection of up to 100 N-glycans per single cell. Interestingly, N-glycome alterations are unequivocally detected at the single-cell level in HeLa and U87 cells stimulated with lipopolysaccharide. The developed workflow is also applied to the profiling of ng-level amounts (5-500 ng) of blood-derived protein, extracellular vesicle, and total plasma isolates, resulting in over 170, 220, and 370 quantitated N-glycans, respectively.

Data normalization for addressing the challenges in the analysis of single-cell transcriptomic datasets

BACKGROUND

Normalization is a critical step in the analysis of single-cell RNA-sequencing (scRNA-seq) datasets. Its main goal is to make gene counts comparable within and between cells. To do so, normalization methods must account for technical and biological variability. Numerous normalization methods have been developed addressing different sources of dispersion and making specific assumptions about the count data.

MAIN BODY

The selection of a normalization method has a direct impact on downstream analysis, for example differential gene expression and cluster identification. Thus, the objective of this review is to guide the reader in making an informed decision on the most appropriate normalization method to use. To this aim, we first give an overview of the different single cell sequencing platforms and methods commonly used including isolation and library preparation protocols. Next, we discuss the inherent sources of variability of scRNA-seq datasets. We describe the categories of normalization methods and include examples of each. We also delineate imputation and batch-effect correction methods. Furthermore, we describe data-driven metrics commonly used to evaluate the performance of normalization methods. We also discuss common scRNA-seq methods and toolkits used for integrated data analysis.

CONCLUSIONS

According to the correction performed, normalization methods can be broadly classified as within and between-sample algorithms. Moreover, with respect to the mathematical model used, normalization methods can further be classified into: global scaling methods, generalized linear models, mixed methods, and machine learning-based methods. Each of these methods depict pros and cons and make different statistical assumptions. However, there is no better performing normalization method. Instead, metrics such as silhouette width, K-nearest neighbor batch-effect test, or Highly Variable Genes are recommended to assess the performance of normalization methods.

Precision and Accuracy of Single-Cell/Nuclei RNA Sequencing Data

Single-cell/nuclei RNA sequencing (sc/snRNA-Seq) is widely used for profiling cell-type gene expressions in biomedical research. An important but underappreciated issue is the quality of sc/snRNA-Seq data that would impact the reliability of downstream analyses. Here we evaluated the precision and accuracy in 18 sc/snRNA-Seq datasets. The precision was assessed on data from human brain studies with a total of 3,483,905 cells from 297 individuals, by utilizing technical replicates. The accuracy was evaluated with sample-matched scRNA-Seq and pooled-cell RNA-Seq data of cultured mononuclear phagocytes from four species. The results revealed low precision and accuracy at the single-cell level across all evaluated data. Cell number and RNA quality were highlighted as two key factors determining the expression precision, accuracy, and reproducibility of differential expression analysis in sc/snRNA-Seq. This study underscores the necessity of sequencing enough high-quality cells per cell type per individual, preferably in the hundreds, to mitigate noise in expression quantification.

Differences in Oxygen-Induced Retinopathy Susceptibility Between Two Sprague Dawley Rat Vendors: A Comparison of Retinal Transcriptomes

PURPOSE

To determine the retinal transcriptomic differences underlying the oxygen-induced retinopathy phenotypes between Sprague Dawley rat pups from two commonly used commercial vendors. This will allow us to discover genes and pathways that may be related to differences in disease severity in similarly aged premature babies and suggest possible new treatment approaches.

METHODS

We analyzed retinal vascular morphometry and transcriptomes from Sprague Dawley rat pups from Charles River Laboratories and Envigo (previously Harlan). Room air control and oxygen-induced retinopathy groups were compared. Oxygen-induced retinopathy was induced with the rat 50/10 model.

RESULTS

Pups from Charles River Laboratories developed a more severe oxygen-induced retinopathy phenotype, with 3.6-fold larger percent avascular area at P15 and twofold larger % neovascular area at P20 than pups from Envigo. Changes in retinal transcriptomes of rat pups from both vendors were substantial at baseline and in response to oxygen-induced retinopathy. Baseline differences centered on activated pathways of neuronal development in Charles River Laboratories pups. In response to oxygen-induced retinopathy, during the neovascular phase, retinas from Charles River Laboratories pups exhibited activation of pathways regulating necrosis, neuroinflammation, and interferon signaling, supporting the observed increase of neovascularization. Conversely, retinas from Envigo pups showed decreased necrosis and increased focal adhesion kinase signaling, supporting more normal vascular development. Comparing oxygen-induced retinopathy transcriptomes at P15 to those at P20, canonical pathways such as phosphate and tensin homolog, interferon, and coordinated lysosomal expression and regulation element signaling were identified, highlighting potential novel mechanistic targets for future research.

CONCLUSION

Transcriptomic profiles differ substantially between rat pup retinas from Charles River Laboratories and Envigo at baseline and in response to oxygen-induced retinopathy, providing insight into vascular morphologic differences. Comparing transcriptomes identified new pathways for further research in oxygen-induced retinopathy pathogenesis and increased scientific rigor of this model.

Modeling single cell trajectory using forward-backward stochastic differential equations

Recent advances in single-cell sequencing technology have provided opportunities for mathematical modeling of dynamic developmental processes at the single-cell level, such as inferring developmental trajectories. Optimal transport has emerged as a promising theoretical framework for this task by computing pairings between cells from different time points. However, optimal transport methods have limitations in capturing nonlinear trajectories, as they are static and can only infer linear paths between endpoints. In contrast, stochastic differential equations (SDEs) offer a dynamic and flexible approach that can model non-linear trajectories, including the shape of the path. Nevertheless, existing SDE methods often rely on numerical approximations that can lead to inaccurate inferences, deviating from true trajectories. To address this challenge, we propose a novel approach combining forward-backward stochastic differential equations (FBSDE) with a refined approximation procedure. Our FBSDE model integrates the forward and backward movements of two SDEs in time, aiming to capture the underlying dynamics of single-cell developmental trajectories. Through comprehensive benchmarking on multiple scRNA-seq datasets, we demonstrate the superior performance of FBSDE compared to other methods, highlighting its efficacy in accurately inferring developmental trajectories.

Computational frameworks integrating deep learning and statistical models in mining multimodal omics data

BACKGROUND

In health research, multimodal omics data analysis is widely used to address important clinical and biological questions. Traditional statistical methods rely on the strong assumptions of distribution. Statistical methods such as testing and differential expression are commonly used in omics analysis. Deep learning, on the other hand, is an advanced computer science technique that is powerful in mining high-dimensional omics data for prediction tasks. Recently, integrative frameworks or methods have been developed for omics studies that combine statistical models and deep learning algorithms.

METHODS AND RESULTS

The aim of these integrative frameworks is to combine the strengths of both statistical methods and deep learning algorithms to improve prediction accuracy while also providing interpretability and explainability. This review report discusses the current state-of-the-art integrative frameworks, their limitations, and potential future directions in survival and time-to-event longitudinal analysis, dimension reduction and clustering, regression and classification, feature selection, and causal and transfer learning.

Aortic Cellular Heterogeneity in Health and Disease: Novel Insights Into Aortic Diseases From Single-Cell RNA Transcriptomic Data Sets

Aortic diseases such as atherosclerosis, aortic aneurysms, and aortic stiffening are significant complications that can have significant impact on end-stage cardiovascular disease. With limited pharmacological therapeutic strategies that target the structural changes in the aorta, surgical intervention remains the only option for some patients with these diseases. Although there have been significant contributions to our understanding of the cellular architecture of the diseased aorta, particularly in the context of atherosclerosis, furthering our insight into the cellular drivers of disease is required. The major cell types of the aorta are well defined; however, the advent of single-cell RNA sequencing provides unrivaled insights into the cellular heterogeneity of each aortic cell type and the inferred biological processes associated with each cell in health and disease. This review discusses previous concepts that have now been enhanced with recent advances made by single-cell RNA sequencing with a focus on aortic cellular heterogeneity.

Recent advances in differential expression analysis for single-cell RNA-seq and spatially resolved transcriptomic studies

Differential expression (DE) analysis is a necessary step in the analysis of single-cell RNA sequencing (scRNA-seq) and spatially resolved transcriptomics (SRT) data. Unlike traditional bulk RNA-seq, DE analysis for scRNA-seq or SRT data has unique characteristics that may contribute to the difficulty of detecting DE genes. However, the plethora of DE tools that work with various assumptions makes it difficult to choose an appropriate one. Furthermore, a comprehensive review on detecting DE genes for scRNA-seq data or SRT data from multi-condition, multi-sample experimental designs is lacking. To bridge such a gap, here, we first focus on the challenges of DE detection, then highlight potential opportunities that facilitate further progress in scRNA-seq or SRT analysis, and finally provide insights and guidance in selecting appropriate DE tools or developing new computational DE methods.

Protocol to study the inheritance and propagation of non-genetically encoded states using barcode decay lineage tracing

Here, we present a protocol to perform barcode decay lineage tracing followed by single-cell transcriptome analysis (BdLT-Seq). We describe steps for BdLT-Seq experimental design, building barcoded episome reporters, performing episome transfection, and barcode retrieval. We then describe procedures for sequencing library construction while providing options for sample multiplexing and data analysis. This BdLT-Seq technique enables the assessment of clonal evolution in a directional manner while preserving isogeneity, thus allowing the comparison of non-genetic molecular features between isogenic cell lineages. For complete details on the use and execution of this protocol, please refer to Shlyakhtina et al. (2023).1.

Highly Multiplexed, Efficient, and Automated Single-Cell MicroRNA Sequencing with Digital Microfluidics

Single-cell microRNA (miRNA) sequencing has allowed for comprehensively studying the abundance and complex networks of miRNAs, which provides insights beyond single-cell heterogeneity into the dynamic regulation of cellular events. Current benchtop-based technologies for single-cell miRNA sequencing are low throughput, limited reaction efficiency, tedious manual operations, and high reagent costs. Here, a highly multiplexed, efficient, integrated, and automated sample preparation platform is introduced for single-cell miRNA sequencing based on digital microfluidics (DMF), named Hiper-seq. The platform integrates major steps and automates the iterative operations of miRNA sequencing library construction by digital control of addressable droplets on the DMF chip. Based on the design of hydrophilic micro-structures and the capability of handling droplets of DMF, multiple single cells can be selectively isolated and subject to sample processing in a highly parallel way, thus increasing the throughput and efficiency for single-cell miRNA measurement. The nanoliter reaction volume of this platform enables a much higher miRNA detection ability and lower reagent cost compared to benchtop methods. It is further applied Hiper-seq to explore miRNAs involved in the ossification of mouse skeletal stem cells after bone fracture and discovered unreported miRNAs that regulate bone repairing.

Integrating single-cell transcriptomics with cellular phenotypes: cell morphology, Ca2+ imaging and electrophysiology

I review recent technological advancements in coupling single-cell transcriptomics with cellular phenotypes including morphology, calcium signaling, and electrophysiology. Single-cell RNA sequencing (scRNAseq) has revolutionized cell type classifications by capturing the transcriptional diversity of cells. A new wave of methods to integrate scRNAseq and biophysical measurements is facilitating the linkage of transcriptomic data to cellular function, which provides physiological insight into cellular states. I briefly discuss critical factors of these phenotypical characterizations such as timescales, information content, and analytical tools. Dedicated sections focus on the integration with cell morphology, calcium imaging, and electrophysiology (patch-seq), emphasizing their complementary roles. I discuss their application in elucidating cellular states, refining cell type classifications, and uncovering functional differences in cell subtypes. To illustrate the practical applications and benefits of these methods, I highlight their use in tissues with excitable cell-types such as the brain, pancreatic islets, and the retina. The potential of combining functional phenotyping with spatial transcriptomics for a detailed mapping of cell phenotypes in situ is explored. Finally, I discuss open questions and future perspectives, emphasizing the need for a shift towards broader accessibility through increased throughput.

CellCommuNet: an atlas of cell-cell communication networks from single-cell RNA sequencing of human and mouse tissues in normal and disease states

Cell-cell communication, as a basic feature of multicellular organisms, is crucial for maintaining the biological functions and microenvironmental homeostasis of cells, organs, and whole organisms. Alterations in cell-cell communication contribute to many diseases, including cancers. Single-cell RNA sequencing (scRNA-seq) provides a powerful method for studying cell-cell communication by enabling the analysis of ligand-receptor interactions. Here, we introduce CellCommuNet (http://www.inbirg.com/cellcommunet/), a comprehensive data resource for exploring cell-cell communication networks in scRNA-seq data from human and mouse tissues in normal and disease states. CellCommuNet currently includes 376 single datasets from multiple sources, and 118 comparison datasets between disease and normal samples originating from the same study. CellCommuNet provides information on the strength of communication between cells and related signalling pathways and facilitates the exploration of differences in cell-cell communication between healthy and disease states. Users can also search for specific signalling pathways, ligand-receptor pairs, and cell types of interest. CellCommuNet provides interactive graphics illustrating cell-cell communication in different states, enabling differential analysis of communication strength between disease and control samples. This comprehensive database aims to be a valuable resource for biologists studying cell-cell communication networks.

Biomarker discovery in acetaminophen hepatotoxicity: leveraging single-cell transcriptomics and mechanistic insight

INTRODUCTION

Acetaminophen (APAP) overdose is the leading cause of drug-induced liver injury and can cause a rapid progression to acute liver failure (ALF). Therefore, the identification of prognostic biomarkers to determine which patients will require a liver transplant is critical for APAP-induced ALF.

AREAS COVERED

We begin by relating the mechanistic investigations in mouse models of APAP hepatotoxicity to the human APAP overdose pathophysiology. We draw insights from the established sequence of molecular events in mice to understand the progression of events in the APAP overdose patient. Through this mechanistic understanding, several new biomarkers, such as CXCL14, have recently been evaluated. We also explore how single-cell RNA sequencing, spatial transcriptomics, and other omics approaches have been leveraged for identifying novel biomarkers and how these approaches will continue to push the field of biomarker discovery forward.

EXPERT OPINION

Recent investigations have elucidated several new biomarkers or combination of markers such as CXCL14, a regenerative miRNA signature, a cell death miRNA signature, hepcidin, LDH, CPS1, and FABP1. While these biomarkers are promising, they all require further validation. Larger cohort studies analyzing these new biomarkers in the same patient samples, while adding these candidate biomarkers to prognostic models will further support their clinical utility.

Cell Size Contributes to Single-Cell Proteome Variation

Accurate measurements of the molecular composition of single cells will be necessary for understanding the relationship between gene expression and function in diverse cell types. One of the most important phenotypes that differs between cells is their size, which was recently shown to be an important determinant of proteome composition in populations of similarly sized cells. We, therefore, sought to test if the effects of the cell size on protein concentrations were also evident in single-cell proteomics data. Using the relative concentrations of a set of reference proteins to estimate a cell's DNA-to-cell volume ratio, we found that differences in the cell size explain a significant amount of cell-to-cell variance in two published single-cell proteome data sets.

Integrating spatial and single-cell transcriptomics data using deep generative models with SpatialScope

The rapid emergence of spatial transcriptomics (ST) technologies is revolutionizing our understanding of tissue spatial architecture and biology. Although current ST methods, whether based on next-generation sequencing (seq-based approaches) or fluorescence in situ hybridization (image-based approaches), offer valuable insights, they face limitations either in cellular resolution or transcriptome-wide profiling. To address these limitations, we present SpatialScope, a unified approach integrating scRNA-seq reference data and ST data using deep generative models. With innovation in model and algorithm designs, SpatialScope not only enhances seq-based ST data to achieve single-cell resolution, but also accurately infers transcriptome-wide expression levels for image-based ST data. We demonstrate SpatialScope's utility through simulation studies and real data analysis from both seq-based and image-based ST approaches. SpatialScope provides spatial characterization of tissue structures at transcriptome-wide single-cell resolution, facilitating downstream analysis, including detecting cellular communication through ligand-receptor interactions, localizing cellular subtypes, and identifying spatially differentially expressed genes.

Single-cell long-read sequencing in human cerebral organoids uncovers cell-type-specific and autism-associated exons

Dysregulation of alternative splicing has been repeatedly associated with neurodevelopmental disorders, but the extent of cell-type-specific splicing in human neural development remains largely uncharted. Here, single-cell long-read sequencing in induced pluripotent stem cell (iPSC)-derived cerebral organoids identifies over 31,000 uncatalogued isoforms and 4,531 cell-type-specific splicing events. Long reads uncover coordinated splicing and cell-type-specific intron retention events, which are challenging to study with short reads. Retained neuronal introns are enriched in RNA splicing regulators, showing shorter lengths, higher GC contents, and weaker 5' splice sites. We use this dataset to explore the biological processes underlying neurological disorders, focusing on autism. In comparison with prior transcriptomic data, we find that the splicing program in autistic brains is closer to the progenitor state than differentiated neurons. Furthermore, cell-type-specific exons harbor significantly more de novo mutations in autism probands than in siblings. Overall, these results highlight the importance of cell-type-specific splicing in autism and neuronal gene regulation.

In-capillary sample processing coupled to label-free capillary electrophoresis-mass spectrometry to decipher the native N-glycome of single mammalian cells and ng-level blood isolates

The development of reliable single-cell dispensers and substantial sensitivity improvement in mass spectrometry made proteomic profiling of individual cells achievable. Yet, there are no established methods for single-cell glycome analysis due to the inability to amplify glycans and sample losses associated with sample processing and glycan labeling. In this work, we developed an integrated platform coupling online in-capillary sample processing with high-sensitivity label-free capillary electrophoresis-mass spectrometry for N-glycan profiling of single mammalian cells. Direct and unbiased characterization and quantification of single-cell surface N-glycomes were demonstrated for HeLa and U87 cells, with the detection of up to 100 N-glycans per single cell. Interestingly, N-glycome alterations were unequivocally detected at the single-cell level in HeLa and U87 cells stimulated with lipopolysaccharide. The developed workflow was also applied to the profiling of ng-level amounts of blood-derived protein, extracellular vesicle, and total plasma isolates, resulting in over 170, 220, and 370 quantitated N-glycans, respectively.

scPrisma infers, filters and enhances topological signals in single-cell data using spectral template matching

Single-cell RNA sequencing has been instrumental in uncovering cellular spatiotemporal context. This task is challenging as cells simultaneously encode multiple, potentially cross-interfering, biological signals. Here we propose scPrisma, a spectral computational method that uses topological priors to decouple, enhance and filter different classes of biological processes in single-cell data, such as periodic and linear signals. We apply scPrisma to the analysis of the cell cycle in HeLa cells, circadian rhythm and spatial zonation in liver lobules, diurnal cycle in Chlamydomonas and circadian rhythm in the suprachiasmatic nucleus in the brain. scPrisma can be used to distinguish mixed cellular populations by specific characteristics such as cell type and uncover regulatory networks and cell-cell interactions specific to predefined biological signals, such as the circadian rhythm. We show scPrisma's flexibility in incorporating prior knowledge, inference of topologically informative genes and generalization to additional diverse templates and systems. scPrisma can be used as a stand-alone workflow for signal analysis and as a prior step for downstream single-cell analysis.

Single-cell RNA Sequencing Reveals Novel Cellular Factors for Response to Immunosuppressive Therapy in Aplastic Anemia

Aplastic anemia (AA) is a lethal hematological disorder; however, its pathogenesis is not fully understood. Although immunosuppressive therapy (IST) is a major treatment option for AA, one-third of patients do not respond to IST and its resistance mechanism remains elusive. To understand AA pathogenesis and IST resistance, we performed single-cell RNA sequencing (scRNA-seq) of bone marrow (BM) from healthy controls and patients with AA at diagnosis. We found that CD34+ early-stage erythroid precursor cells and PROM1+ hematopoietic stem cells were significantly depleted in AA, which suggests that the depletion of CD34+ early-stage erythroid precursor cells and PROM1+ hematopoietic stem cells might be one of the major mechanisms for AA pathogenesis related with BM-cell hypoplasia. More importantly, we observed the significant enrichment of CD8+ T cells and T cell-activating intercellular interactions in IST responders, indicating the association between the expansion and activation of T cells and the positive response of IST in AA. Taken together, our findings represent a valuable resource offering novel insights into the cellular heterogeneity in the BM of AA and reveal potential biomarkers for IST, building the foundation for future precision therapies in AA.

Non-Mass Spectrometric Targeted Single-Cell Metabolomics

Metabolic assays serve as pivotal tools in biomedical research, offering keen insights into cellular physiological and pathological states. While mass spectrometry (MS)-based metabolomics remains the gold standard for comprehensive, multiplexed analyses of cellular metabolites, innovative technologies are now emerging for the targeted, quantitative scrutiny of metabolites and metabolic pathways at the single-cell level. In this review, we elucidate an array of these advanced methodologies, spanning synthetic and surface chemistry techniques, imaging-based methods, and electrochemical approaches. We summarize the rationale, design principles, and practical applications for each method, and underscore the synergistic benefits of integrating single-cell metabolomics (scMet) with other single-cell omics technologies. Concluding, we identify prevailing challenges in the targeted scMet arena and offer a forward-looking commentary on future avenues and opportunities in this rapidly evolving field.

It's high-time to re-evaluate the value of induced-chemotherapy for reinforcing immunotherapy in colorectal cancer

Immunotherapy has made significant advances in the treatment of colorectal cancer (CRC), revolutionizing the therapeutic landscape and highlighting the indispensable role of the tumor immune microenvironment. However, some CRCs have shown poor response to immunotherapy, prompting investigation into the underlying reasons. It has been discovered that certain chemotherapeutic agents possess immune-stimulatory properties, including the induction of immunogenic cell death (ICD), the generation and processing of non-mutated neoantigens (NM-neoAgs), and the B cell follicle-driven T cell response. Based on these findings, the concept of inducing chemotherapy has been introduced, and the combination of inducing chemotherapy and immunotherapy has become a standard treatment option for certain cancers. Clinical trials have confirmed the feasibility and safety of this approach in CRC, offering a promising method for improving the efficacy of immunotherapy. Nevertheless, there are still many challenges and difficulties ahead, and further research is required to optimize its use.

Identifying Subpopulations of Cells in Single-Cell Transcriptomic Data: A Bayesian Mixture Modeling Approach to Zero Inflation of Counts

In the study of single-cell RNA-seq (scRNA-Seq) data, a key component of the analysis is to identify subpopulations of cells in the data. A variety of approaches to this have been considered, and although many machine learning-based methods have been developed, these rarely give an estimate of uncertainty in the cluster assignment. To allow for this, probabilistic models have been developed, but scRNA-Seq data exhibit a phenomenon known as dropout, whereby a large proportion of the observed read counts are zero. This poses challenges in developing probabilistic models that appropriately model the data. We develop a novel Dirichlet process mixture model that employs both a mixture at the cell level to model multiple populations of cells and a zero-inflated negative binomial mixture of counts at the transcript level. By taking a Bayesian approach, we are able to model the expression of genes within clusters, and to quantify uncertainty in cluster assignments. It is shown that this approach outperforms previous approaches that applied multinomial distributions to model scRNA-Seq counts and negative binomial models that do not take into account zero inflation. Applied to a publicly available data set of scRNA-Seq counts of multiple cell types from the mouse cortex and hippocampus, we demonstrate how our approach can be used to distinguish subpopulations of cells as clusters in the data, and to identify gene sets that are indicative of membership of a subpopulation.

Comprehensive analyses of brain cell communications based on multiple scRNA-seq and snRNA-seq datasets for revealing novel mechanism in neurodegenerative diseases

AIMS

Complex cellular communications between glial cells and neurons are critical for brain normal function and disorders, and single-cell level RNA-sequencing datasets display more advantages for analyzing cell communications. Therefore, it is necessary to systematically explore brain cell communications when considering factors such as sex and brain region.

METHODS

We extracted a total of 1,039,459 cells derived from 28 brain single-cell RNA-sequencing (scRNA-seq) or single-nucleus RNA-sequencing (snRNA-seq) datasets from the GEO database, including 12 human and 16 mouse datasets. These datasets were further divided into 71 new sub-datasets when considering disease, sex, and region conditions. In the meanwhile, we integrated four methods to evaluate ligand-receptor interaction score among six major brain cell types (microglia, neuron, astrocyte, oligodendrocyte, OPC, and endothelial cell).

RESULTS

For Alzheimer's disease (AD), disease-specific ligand-receptor pairs when compared with normal sub-datasets, such as SEMA4A-NRP1, were identified. Furthermore, we explored the sex- and region-specific cell communications and identified that WNT5A-ROR1 among microglia cells displayed close communications in male, and SPP1-ITGAV displayed close communications in the meninges region from microglia to neurons. Furthermore, based on the AD-specific cell communications, we constructed a model for AD early prediction and confirmed the predictive performance using multiple independent datasets. Finally, we developed an online platform for researchers to explore brain condition-specific cell communications.

CONCLUSION

This research provided a comprehensive study to explore brain cell communications, which could reveal novel biological mechanisms involved in normal brain function and neurodegenerative diseases such as AD.

FEED: a feature selection method based on gene expression decomposition for single cell clustering

Single-cell clustering is a critical step in biological downstream analysis. The clustering performance could be effectively improved by extracting cell-type-specific genes. The state-of-the-art feature selection methods usually calculate the importance of a single gene without considering the information contained in the gene expression distribution. Moreover, these methods ignore the intrinsic expression patterns of genes and heterogeneity within groups of different mean expression levels. In this work, we present a Feature sElection method based on gene Expression Decomposition (FEED) of scRNA-seq data, which selects informative genes to enhance clustering performance. First, the expression levels of genes are decomposed into multiple Gaussian components. Then, a novel gene correlation calculation method is proposed to measure the relationship between genes from the perspective of distribution. Finally, a permutation-based approach is proposed to determine the threshold of gene importance to obtain marker gene subsets. Compared with state-of-the-art feature selection methods, applying FEED on various scRNA-seq datasets including large datasets followed by different common clustering algorithms results in significant improvements in the accuracy of cell-type identification. The source codes for FEED are freely available at https://github.com/genemine/FEED.

Data Mining of Microarray Datasets in Translational Neuroscience

Data mining involves the computational analysis of a plethora of publicly available datasets to generate new hypotheses that can be further validated by experiments for the improved understanding of the pathogenesis of neurodegenerative diseases. Although the number of sequencing datasets is on the rise, microarray analysis conducted on diverse biological samples represent a large collection of datasets with multiple web-based programs that enable efficient and convenient data analysis. In this review, we first discuss the selection of biological samples associated with neurological disorders, and the possibility of a combination of datasets, from various types of samples, to conduct an integrated analysis in order to achieve a holistic understanding of the alterations in the examined biological system. We then summarize key approaches and studies that have made use of the data mining of microarray datasets to obtain insights into translational neuroscience applications, including biomarker discovery, therapeutic development, and the elucidation of the pathogenic mechanisms of neurodegenerative diseases. We further discuss the gap to be bridged between microarray and sequencing studies to improve the utilization and combination of different types of datasets, together with experimental validation, for more comprehensive analyses. We conclude by providing future perspectives on integrating multi-omics, to advance precision phenotyping and personalized medicine for neurodegenerative diseases.

A contamination focused approach for optimizing the single-cell RNA-seq experiment

Droplet-based single-cell RNA-seq (scRNA-seq) data are plagued by ambient contaminations caused by nucleic acid material released by dead and dying cells. This material is mixed into the buffer and is co-encapsulated with cells, leading to a lower signal-to-noise ratio. Although there exist computational methods to remove ambient contaminations post-hoc, the reliability of algorithms in generating high-quality data from low-quality sources remains uncertain. Here, we assess data quality before data filtering by a set of quantitative, contamination-based metrics that assess data quality more effectively than standard metrics. Through a series of controlled experiments, we report improvements that can minimize ambient contamination outside of tissue dissociation, via cell fixation, improved cell loading, microfluidic dilution, and nuclei versus cell preparation; many of these parameters are inaccessible on commercial platforms. We provide end-users with insights on factors that can guide their decision-making regarding optimizations that minimize ambient contamination, and metrics to assess data quality.

The ENCODE Uniform Analysis Pipelines

The Encyclopedia of DNA elements (ENCODE) project is a collaborative effort to create a comprehensive catalog of functional elements in the human genome. The current database comprises more than 19000 functional genomics experiments across more than 1000 cell lines and tissues using a wide array of experimental techniques to study the chromatin structure, regulatory and transcriptional landscape of the Homo sapiens and Mus musculus genomes. All experimental data, metadata, and associated computational analyses created by the ENCODE consortium are submitted to the Data Coordination Center (DCC) for validation, tracking, storage, and distribution to community resources and the scientific community. The ENCODE project has engineered and distributed uniform processing pipelines in order to promote data provenance and reproducibility as well as allow interoperability between genomic resources and other consortia. All data files, reference genome versions, software versions, and parameters used by the pipelines are captured and available via the ENCODE Portal. The pipeline code, developed using Docker and Workflow Description Language (WDL; https://openwdl.org/) is publicly available in GitHub, with images available on Dockerhub (https://hub.docker.com), enabling access to a diverse range of biomedical researchers. ENCODE pipelines maintained and used by the DCC can be installed to run on personal computers, local HPC clusters, or in cloud computing environments via Cromwell. Access to the pipelines and data via the cloud allows small labs the ability to use the data or software without access to institutional compute clusters. Standardization of the computational methodologies for analysis and quality control leads to comparable results from different ENCODE collections - a prerequisite for successful integrative analyses.

Inverse agonists of retinoic acid receptor/retinoid X receptor signaling as lineage-specific antitumor agents against human adenoid cystic carcinoma

BACKGROUND

Adenoid cystic carcinoma (ACC) is a lethal malignancy of exocrine glands, characterized by the coexistence within tumor tissues of 2 distinct populations of cancer cells, phenotypically similar to the myoepithelial and ductal lineages of normal salivary epithelia. The developmental relationship linking these 2 cell types, and their differential vulnerability to antitumor treatments, remains unknown.

METHODS

Using single-cell RNA sequencing, we identified cell-surface markers (CD49f, KIT) that enabled the differential purification of myoepithelial-like (CD49fhigh/KITneg) and ductal-like (CD49flow/KIT+) cells from patient-derived xenografts (PDXs) of human ACCs. Using prospective xenotransplantation experiments, we compared the tumor-initiating capacity of the 2 cell types and tested whether one could differentiate into the other. Finally, we searched for signaling pathways with differential activation between the 2 cell types and tested their role as lineage-specific therapeutic targets.

RESULTS

Myoepithelial-like cells displayed higher tumorigenicity than ductal-like cells and acted as their progenitors. Myoepithelial-like and ductal-like cells displayed differential expression of genes encoding for suppressors and activators of retinoic acid signaling, respectively. Agonists of retinoic acid receptor (RAR) or retinoid X receptor (RXR) signaling (all-trans retinoic acid, bexarotene) promoted myoepithelial-to-ductal differentiation, whereas suppression of RAR/RXR signaling with a dominant-negative RAR construct abrogated it. Inverse agonists of RAR/RXR signaling (BMS493, AGN193109) displayed selective toxicity against ductal-like cells and in vivo antitumor activity against PDX models of human ACC.

CONCLUSIONS

In human ACCs, myoepithelial-like cells act as progenitors of ductal-like cells, and myoepithelial-to-ductal differentiation is promoted by RAR/RXR signaling. Suppression of RAR/RXR signaling is lethal to ductal-like cells and represents a new therapeutic approach against human ACCs.

Robust reconstruction of single-cell RNA-seq data with iterative gene weight updates

MOTIVATION

Single-cell RNA-sequencing technologies have greatly enhanced our understanding of heterogeneous cell populations and underlying regulatory processes. However, structural (spatial or temporal) relations between cells are lost during cell dissociation. These relations are crucial for identifying associated biological processes. Many existing tissue-reconstruction algorithms use prior information about subsets of genes that are informative with respect to the structure or process to be reconstructed. When such information is not available, and in the general case when the input genes code for multiple processes, including being susceptible to noise, biological reconstruction is often computationally challenging.

RESULTS

We propose an algorithm that iteratively identifies manifold-informative genes using existing reconstruction algorithms for single-cell RNA-seq data as subroutine. We show that our algorithm improves the quality of tissue reconstruction for diverse synthetic and real scRNA-seq data, including data from the mammalian intestinal epithelium and liver lobules.

AVAILABILITY AND IMPLEMENTATION

The code and data for benchmarking are available at github.com/syq2012/iterative_weight_update_for_reconstruction.

Characterization of cell-type specific circular RNAs associated with colorectal cancer metastasis

Colorectal cancer (CRC) is the most common gastrointestinal malignancy and a leading cause of cancer deaths in the United States. More than half of CRC patients develop metastatic disease (mCRC) with an average 5-year survival rate of 13%. Circular RNAs (circRNAs) have recently emerged as important tumorigenesis regulators; however, their role in mCRC progression remains poorly characterized. Further, little is known about their cell-type specificity to elucidate their functions in the tumor microenvironment (TME). To address this, we performed total RNA sequencing (RNA-seq) on 30 matched normal, primary and metastatic samples from 14 mCRC patients. Additionally, five CRC cell lines were sequenced to construct a circRNA catalog in CRC. We detected 47 869 circRNAs, with 51% previously unannotated in CRC and 14% novel candidates when compared to existing circRNA databases. We identified 362 circRNAs differentially expressed in primary and/or metastatic tissues, termed circular RNAs associated with metastasis (CRAMS). We performed cell-type deconvolution using published single-cell RNA-seq datasets and applied a non-negative least squares statistical model to estimate cell-type specific circRNA expression. This predicted 667 circRNAs as exclusively expressed in a single cell type. Collectively, this serves as a valuable resource, TMECircDB (accessible at https://www.maherlab.com/tmecircdb-overview), for functional characterization of circRNAs in mCRC, specifically in the TME.

Transcriptome-Powered Pluripotent Stem Cell Differentiation for Regenerative Medicine

Pluripotent stem cells are endless sources for in vitro engineering human tissues for regenerative medicine. Extensive studies have demonstrated that transcription factors are the key to stem cell lineage commitment and differentiation efficacy. As the transcription factor profile varies depending on the cell type, global transcriptome analysis through RNA sequencing (RNAseq) has been a powerful tool for measuring and characterizing the success of stem cell differentiation. RNAseq has been utilized to comprehend how gene expression changes as cells differentiate and provide a guide to inducing cellular differentiation based on promoting the expression of specific genes. It has also been utilized to determine the specific cell type. This review highlights RNAseq techniques, tools for RNAseq data interpretation, RNAseq data analytic methods and their utilities, and transcriptomics-enabled human stem cell differentiation. In addition, the review outlines the potential benefits of the transcriptomics-aided discovery of intrinsic factors influencing stem cell lineage commitment, transcriptomics applied to disease physiology studies using patients' induced pluripotent stem cell (iPSC)-derived cells for regenerative medicine, and the future outlook on the technology and its implementation.

Small RNA transcriptome analysis using parallel single-cell small RNA sequencing

miRNA and other forms of small RNAs are known to regulate many biological processes. Single-cell small RNA sequencing can be used to profile small RNAs of individual cells; however, limitations of efficiency and scale prevent its widespread application. Here, we developed parallel single-cell small RNA sequencing (PSCSR-seq), which can overcome the limitations of existing methods and enable high-throughput small RNA expression profiling of individual cells. Analysis of PSCSR-seq data indicated that diverse cell types could be identified based on patterns of miRNA expression, and showed that miRNA content in nuclei is informative (for example, cell type marker miRNAs can be detected in isolated nuclei). PSCSR-seq is very sensitive: analysis of only 732 peripheral blood mononuclear cells (PBMCs) detected 774 miRNAs, whereas bulk small RNA analysis would require input RNA from approximately 10⁶ cells to detect as many miRNAs. We identified 42 miRNAs as markers for PBMC subpopulations. Moreover, we analyzed the miRNA profiles of 9,533 cells from lung cancer biopsies, and by dissecting cell subpopulations, we identified potentially diagnostic and therapeutic miRNAs for lung cancers. Our study demonstrates that PSCSR-seq is highly sensitive and reproducible, thus making it an advanced tool for miRNA analysis in cancer and life science research.

The ENCODE Uniform Analysis Pipelines

The Encyclopedia of DNA elements (ENCODE) project is a collaborative effort to create a comprehensive catalog of functional elements in the human genome. The current database comprises more than 19000 functional genomics experiments across more than 1000 cell lines and tissues using a wide array of experimental techniques to study the chromatin structure, regulatory and transcriptional landscape of the Homo sapiens and Mus musculus genomes. All experimental data, metadata, and associated computational analyses created by the ENCODE consortium are submitted to the Data Coordination Center (DCC) for validation, tracking, storage, and distribution to community resources and the scientific community. The ENCODE project has engineered and distributed uniform processing pipelines in order to promote data provenance and reproducibility as well as allow interoperability between genomic resources and other consortia. All data files, reference genome versions, software versions, and parameters used by the pipelines are captured and available via the ENCODE Portal. The pipeline code, developed using Docker and Workflow Description Language (WDL; https://openwdl.org/) is publicly available in GitHub, with images available on Dockerhub (https://hub.docker.com), enabling access to a diverse range of biomedical researchers. ENCODE pipelines maintained and used by the DCC can be installed to run on personal computers, local HPC clusters, or in cloud computing environments via Cromwell. Access to the pipelines and data via the cloud allows small labs the ability to use the data or software without access to institutional compute clusters. Standardization of the computational methodologies for analysis and quality control leads to comparable results from different ENCODE collections - a prerequisite for successful integrative analyses.

Evaluating and Validating Sunflower Reference Genes for Q-PCR Studies Under High Temperature Condition

Background

Q-PCR is the method of choice for PCR- based transcriptomics and validating microarray-based and RNA-seq results. Proper application of this technology requires proper normalization to correct as much as possible errors propagating during RNA extraction and cDNA synthesis.

Objectives

The investigation was performed to find stable reference genes in sunflower under shifting in ambient temperature.

Materials and Methods

Sequences of five well-known reference genes of Arabidopsis (Actin, Ubiquitin, Elongation factor-1, GAPDH, and SAND) and one well-known reference gene inhuman, Importin, were subjected to BLASTX against sunflower databases and the relevant genes were subjected to primer designing for q-PCR. Two sunflower inbred lines were cultivated at two dates so that anthesis occurred at nearly 30 °C and 40 °C (heat stress). The experiment was repeated for two years. Q-PCR was run on samples taken for two planting date separately at the beginning of anthesis for each genotype from leaf, taproots, receptacle base, immature and mature disc flowers and on pooled samples comprising of the tissues for each genotype, planting dates and also all tissues for both genotypes and both planting dates. Basic statistical properties of each candidate gene across all the samples were calculated. Furthermore, gene expression stability analysis was done for six candidate reference genes on Cq mean of two years using three independent algorithms, geNorm, Bestkeeper, and Refinder.

Results

Designed primers for Actin2, SAND, GAPDH, Ubiquitin, EF-1a, and Importin yielded a single peak in melting curve analysis indicating specificity of the PCR reaction. Basic statistical analysis showed that Actin2 and EF-1a had the highest and lowest expression levels across all the samples, respectively. Actin2 appeared to be the most stable reference gene across all the samples based on the three used algorithms. Pairwise variation analysis revealed that for samples taken under ambient temperature of 30 °C, Actin2, EF-1a, SAND and for those taken under ambient temperature of 40 °C, Actin2, EF-1a, Importin and SAND have to be used for normalization in q-PCR studies. Moreover, it is suggested that normalization to be based on Actin2, SAND and EF-1a for vegetative tissues and Actin2, EF-1a, SAND and Importin for reproductive tissues.

Conclusions

In the present research, proper reference genes for normalization of gene expression studies under heat stress conditions were introduced. Moreover, the presence of genotype-by- planting date interaction effects and tissue specific gene expression pattern on the behavior of the most three stable reference genes was indicated.

Femtosecond laser microdissection for isolation of regenerating C. elegans neurons for single-cell RNA sequencing

Our understanding of nerve regeneration can be enhanced by delineating its underlying molecular activities at single-neuron resolution in model organisms such as Caenorhabditis elegans. Existing cell isolation techniques cannot isolate neurons with specific regeneration phenotypes from C. elegans. We present femtosecond laser microdissection (fs-LM), a single-cell isolation method that dissects specific cells directly from living tissue by leveraging the micrometer-scale precision of fs-laser ablation. We show that fs-LM facilitates sensitive and specific gene expression profiling by single-cell RNA sequencing (scRNA-seq), while mitigating the stress-related transcriptional artifacts induced by tissue dissociation. scRNA-seq of fs-LM isolated regenerating neurons revealed transcriptional programs that are correlated with either successful or failed regeneration in wild-type and dlk-1 (0) animals, respectively. This method also allowed studying heterogeneity displayed by the same type of neuron and found gene modules with expression patterns correlated with axon regrowth rate. Our results establish fs-LM as a spatially resolved single-cell isolation method for phenotype-to-genotype mapping.