Single-cell RNA-seq reveals novel regulators of human embryonic stem cell differentiation to definitive endoderm

Gene trajectory inference for single-cell data by optimal transport metrics

Single-cell RNA sequencing has been widely used to investigate cell state transitions and gene dynamics of biological processes. Current strategies to infer the sequential dynamics of genes in a process typically rely on constructing cell pseudotime through cell trajectory inference. However, the presence of concurrent gene processes in the same group of cells and technical noise can obscure the true progression of the processes studied. To address this challenge, we present GeneTrajectory, an approach that identifies trajectories of genes rather than trajectories of cells. Specifically, optimal transport distances are calculated between gene distributions across the cell-cell graph to extract gene programs and define their gene pseudotemporal order. Here we demonstrate that GeneTrajectory accurately extracts progressive gene dynamics in myeloid lineage maturation. Moreover, we show that GeneTrajectory deconvolves key gene programs underlying mouse skin hair follicle dermal condensate differentiation that could not be resolved by cell trajectory approaches. GeneTrajectory facilitates the discovery of gene programs that control the changes and activities of biological processes.

Proteomic Approach Using DIA-MS Identifies Morphogenesis-Associated Proteins during Cardiac Differentiation of Human iPS Cells

Human-induced pluripotent stem cell (hiPSC)-derived cardiomyocytes have potential applications in regenerative medicine. The quality by design (QbD) approach enables the efficiency and quality assurance in the manufacturing of hiPSC-derived products. It requires a molecular understanding of hiPSC differentiation throughout the differentiation process; however, information on cardiac differentiation remains limited. Proteins associated with the early stages of cardiac differentiation would be useful in the cardiomyocyte quality assessment. Here, we performed quantitative proteomics of hiPSC intermediate cells in the early phase of cardiac differentiation to better understand their molecular characteristics. Proteomic profiles suggested that day 5-7 cells were in the morphogenetic stage of cardiac differentiation. Trophoblast glycoprotein (TPBG) was the most up-regulated protein in the morphogenetic stage; it was previously shown to be up-regulated during differentiation into neural stem cells. Proteomics of TPBG-knockdown cells revealed that TPBG is involved in cell proliferation and is related to the cardiomyocyte yield, suggesting that it could be used as a marker in QbD development. Our approach helps us understand the molecular basis of hiPSC differentiation and could be a powerful tool in QbD-based manufacturing.

GraCEImpute: A novel graph clustering autoencoder approach for imputation of single-cell RNA-seq data

Single-cell RNA sequencing (scRNA-seq) technology establishes a unique view for elucidating cellular heterogeneity in various biological systems. Yet the scRNA-seq data is compromised by a high dropout rate due to the technological limitation, and the substantial data loss poses computational challenges on subsequent analyses. This study introduces a novel graph clustering autoencoder (GCAE)-based imputation approach (GraCEImpute) to address the challenge of missing data in scRNA-seq data. Our comprehensive evaluation demonstrates that the GraCEImpute model outperforms existing approaches in accurately imputing dropout zeros within scRNA-seq data. The proposed GraCEImpute model also demonstrates the significantly enhanced quality of downstream scRNA-seq data analyses, including clustering, differential gene expression (DEG) analysis, and cell trajectory inference. These improvements underscore the GraCEImpute model's potential to facilitate a deeper understanding of cellular processes and heterogeneity through the scRNA-seq data analyses. The source code is released at https://www.healthinformaticslab.org/supp/.

Gene regulatory network inference based on novel ensemble method

Gene regulatory networks (GRNs) contribute toward understanding the function of genes and the development of cancer or the impact of key genes on diseases. Hence, this study proposes an ensemble method based on 13 basic classification methods and a flexible neural tree (FNT) to improve GRN identification accuracy. The primary classification methods contain ridge classification, stochastic gradient descent, Gaussian process classification, Bernoulli Naive Bayes, adaptive boosting, gradient boosting decision tree, hist gradient boosting classification, eXtreme gradient boosting (XGBoost), multilayer perceptron, light gradient boosting machine, random forest, support vector machine, and k-nearest neighbor algorithm, which are regarded as the input variable set of FNT model. Additionally, a hybrid evolutionary algorithm based on a gene programming variant and particle swarm optimization is developed to search for the optimal FNT model. Experiments on three simulation datasets and three real single-cell RNA-seq datasets demonstrate that the proposed ensemble feature outperforms 13 supervised algorithms, seven unsupervised algorithms (ARACNE, CLR, GENIE3, MRNET, PCACMI, GENECI, and EPCACMI) and four single cell-specific methods (SCODE, BiRGRN, LEAP, and BiGBoost) based on the area under the receiver operating characteristic curve, area under the precision-recall curve, and F1 metrics.

CVGAE: A Self-Supervised Generative Method for Gene Regulatory Network Inference Using Single-Cell RNA Sequencing Data

Gene regulatory network (GRN) inference based on single-cell RNA sequencing data (scRNAseq) plays a crucial role in understanding the regulatory mechanisms between genes. Various computational methods have been employed for GRN inference, but their performance in terms of network accuracy and model generalization is not satisfactory, and their poor performance is caused by high-dimensional data and network sparsity. In this paper, we propose a self-supervised method for gene regulatory network inference using single-cell RNA sequencing data (CVGAE). CVGAE uses graph neural network for inductive representation learning, which merges gene expression data and observed topology into a low-dimensional vector space. The well-trained vectors will be used to calculate mathematical distance of each gene, and further predict interactions between genes. In overall framework, FastICA is implemented to relief computational complexity caused by high dimensional data, and CVGAE adopts multi-stacked GraphSAGE layers as an encoder and an improved decoder to overcome network sparsity. CVGAE is evaluated on several single cell datasets containing four related ground-truth networks, and the result shows that CVGAE achieve better performance than comparative methods. To validate learning and generalization capabilities, CVGAE is applied in few-shot environment by change the ratio of train set and test set. In condition of few-shot, CVGAE obtains comparable or superior performance.

Supervised analysis of alternative polyadenylation from single-cell and spatial transcriptomics data with spvAPA

Alternative polyadenylation (APA) is an important driver of transcriptome diversity that generates messenger RNA isoforms with distinct 3' ends. The rapid development of single-cell and spatial transcriptomic technologies opened up new opportunities for exploring APA data to discover hidden cell subpopulations invisible in conventional gene expression analysis. However, conventional gene-level analysis tools are not fully applicable to APA data, and commonly used unsupervised dimensionality reduction methods often disregard experimentally derived annotations such as cell type identities. Here, we proposed a supervised analytical framework termed spvAPA, specifically used for APA analysis from both single-cell and spatial transcriptomics data. First, an iterative imputation method based on weighted nearest neighbor was designed to recover missing APA signatures, by integrating both gene expression and APA modalities. Second, a supervised feature selection method based on sparse partial least squares discriminant analysis was devised to identify APA features distinguishing cell types or spatial morphologies. Additionally, spvAPA improves the visualization of high-dimensional data for discovering novel cell subtypes, which considers APA features and dual modalities of gene expression and APA. Evaluations across nine single-cell and spatial transcriptomics datasets demonstrate the effectiveness and applicability of spvAPA. spvAPA is available at https://github.com/BMILAB/spvAPA.

Approaches for Benchmarking Single-Cell Gene Regulatory Network Methods

Gene regulatory networks are powerful tools for modeling genetic interactions that control the expression of genes driving cell differentiation, and single-cell sequencing offers a unique opportunity to build these networks with high-resolution genomic data. There are many proposed computational methods to build these networks using single-cell data, and different approaches are used to benchmark these methods. However, a comprehensive discussion specifically focusing on benchmarking approaches is missing. In this article, we lay the GRN terminology, present an overview of common gold-standard studies and data sets, and define the performance metrics for benchmarking network construction methodologies. We also point out the advantages and limitations of different benchmarking approaches, suggest alternative ground truth data sets that can be used for benchmarking, and specify additional considerations in this context.

WENDY: Covariance dynamics based gene regulatory network inference

Determining gene regulatory network (GRN) structure is a central problem in biology, with a variety of inference methods available for different types of data. For a widely prevalent and challenging use case, namely single-cell gene expression data measured after intervention at multiple time points with unknown joint distributions, there is only one known specifically developed method, which does not fully utilize the rich information contained in this data type. We develop an inference method for the GRN in this case, netWork infErence by covariaNce DYnamics, dubbed WENDY. The core idea of WENDY is to model the dynamics of the covariance matrix, and solve this dynamics as an optimization problem to determine the regulatory relationships. To evaluate its effectiveness, we compare WENDY with other inference methods using synthetic data and experimental data. Our results demonstrate that WENDY performs well across different data sets.

Human iPSC-derived neural stem cells displaying radial glia signature exhibit long-term safety in mice

Human induced pluripotent stem cell-derived neural stem/progenitor cells (hiPSC-NSCs) hold promise for treating neurodegenerative and demyelinating disorders. However, comprehensive studies on their identity and safety remain limited. In this study, we demonstrate that hiPSC-NSCs adopt a radial glia-associated signature, sharing key epigenetic and transcriptional characteristics with human fetal neural stem cells (hfNSCs) while exhibiting divergent profiles from glioblastoma stem cells. Long-term transplantation studies in mice showed robust and stable engraftment of hiPSC-NSCs, with predominant differentiation into glial cells and no evidence of tumor formation. Additionally, we identified the Sterol Regulatory Element Binding Transcription Factor 1 (SREBF1) as a regulator of astroglial differentiation in hiPSC-NSCs. These findings provide valuable transcriptional and epigenetic reference datasets to prospectively define the maturation stage of NSCs derived from different hiPSC sources and demonstrate the long-term safety of hiPSC-NSCs, reinforcing their potential as a viable alternative to hfNSCs for clinical applications.

Spatiotemporal requirements of nuclear β-catenin define early sea urchin embryogenesis

Establishment of the 3 primordial germ layers (ectoderm, endoderm, and mesoderm) during early animal development represents an essential prerequisite for the emergence of properly patterned embryos. β-catenin is an ancient protein that is known to play essential roles in this process. However, these roles have chiefly been established through inhibition of β-catenin translation or function at the time of fertilization. Comprehensive analyses reporting the totality of functions played by nuclear β-catenin during early embryogenesis of a given animal, i.e., at different developmental stages and in different germ layers, are thus still lacking. In this study, we used an inducible, conditional knockdown system in the sea urchin to characterize all possible requirements of β-catenin for germ layer establishment and patterning. By blocking β-catenin protein production starting at 7 different time points of early development, between fertilization and 12 h post fertilization, we established a clear correlation between the position of a germ layer along the primary embryonic axis (the animal-vegetal axis) and its dependence on nuclear β-catenin activity. For example, in the vegetal hemisphere, we determined that the 3 germ layers (skeletogenic mesoderm, non-skeletogenic mesoderm, and endoderm) require distinct and highly specific durations of β-catenin production for their respective specification, with the most vegetal germ layer, the skeletogenic mesoderm, requiring the shortest duration. Likewise, for the 2 animal territories (ectoderm and anterior neuroectoderm), we established that their restriction, along the animal-vegetal axis, relies on different durations of β-catenin production and that the longest duration is required for the most animal territory, the anterior neuroectoderm. Moreover, we found that 2 of the vegetal germ layers, the non-skeletogenic mesoderm and the endoderm, further require a prolonged period of nuclear β-catenin activity after their specification to maintain their respective germ layer identities through time. Finally, we determined that restriction of the anterior neuroectoderm territory depends on at least 2 nuclear β-catenin-dependent inputs and a nuclear β-catenin-independent mechanism. Taken together, this work is the first to comprehensively define the spatiotemporal requirements of β-catenin during the early embryogenesis of a single animal, the sea urchin Paracentrotus lividus, thereby providing new experimental evidence for a better understanding of the roles played by this evolutionary conserved protein during animal development.

Comparison of the developmental competence of in vitro-produced mouse embryos cultured under 5 versus 2% O2 with in vivo-derived blastocysts

PURPOSE

The prevalence of infertility in Canada has substantially increased over 30 years, and plateaued success rates of culture systems warrant further optimization for transfer outcomes. In clinical programs, embryos commonly undergo extended culture under 5% O2 until the blastocyst stage. The aim of this study is to characterize the developmental competence and stress-related responses of embryos cultured under 5 versus 2% O2 in comparison to in vivo-derived blastocysts. We hypothesized 2% O2 compromises developmental competence through altered embryonic stress responses and induction of apoptosis-related genes relative to those cultured under 5% O2 and in vivo-derived blastocysts.

METHODS

Quantitative measures of development and relative expressions of a cohort of stress-related genes in CD1 mouse zygotes cultured to blastocysts under 5 or 2% O2 were compared to in vivo-derived embryos. Apoptotic responses were evaluated using an immunofluorescence assay for Caspase-3.

RESULTS

The mean percentage of blastocysts developed, and total cell number of embryos derived in vivo or cultured under 5% O2 was significantly higher than those cultured under 2% O2. Blastocyst expansion was greatest in embryos cultured under 5% O2. Stress response genes were significantly upregulated in embryos cultured under 2% O2, and expression of antioxidant-related genes was significantly lower in cultured versus in vivo-derived embryos. Caspase-3 immunofluorescence was significantly higher in cultured embryos versus in vivo-derived embryos.

CONCLUSION

We inferred that 5% O2 systems better approximate physiologic oxygen availability for culture of mouse embryos, warranting re-evaluation of culturing embryos under threshold or sub-physiologic oxygen concentrations during clinical IVF programs.

ALKBH5 governs human endoderm fate by regulating the DKK1/4-mediated Wnt/β-catenin activation

N6-methyladenonsine (m6A) is ubiquitously distributed in mammalian mRNA. However, the precise involvement of m6A in early development has yet to be fully elucidated. Here, we report that deletion of the m6A demethylase ALKBH5 in human embryonic stem cells (hESCs) severely impairs definitive endoderm (DE) differentiation. ALKBH5-/- hESCs fail to undergo the primitive streak (PS) intermediate transition that precedes endoderm specification. Mechanistically, we show that ALKBH5 deficiency induces m6A hypermethylation around the 3' untranslated region (3'UTR) of GATA6 transcripts and destabilizes GATA6 mRNA in a YTHDF2-dependent manner. Moreover, GATA6 binds to the promoters of critical regulatory genes involved in Wnt/β-catenin signaling transduction, including the canonical Wnt antagonist DKK1 and DKK4, which are unexpectedly repressed upon the dysregulation of GATA6 mRNA metabolism. Remarkably, DKK1 and DKK4 both exhibit a pleiotropic effect in modulating the Wnt/β-catenin cascade and guard the endogenous signaling activation underlying DE formation as potential downstream targets of the ALKBH5-GATA6 regulation. Here, we unravel a role of ALKBH5 in human endoderm formation in vitro by modulating the canonical Wnt signaling logic through the previously unrecognized functions of DKK1/4, thus capturing a more comprehensive role of m6A in early human embryogenesis.

COFFEE: consensus single cell-type specific inference for gene regulatory networks

The inference of gene regulatory networks (GRNs) is crucial to understanding the regulatory mechanisms that govern biological processes. GRNs may be represented as edges in a graph, and hence, it have been inferred computationally for scRNA-seq data. A wisdom of crowds approach to integrate edges from several GRNs to create one composite GRN has demonstrated improved performance when compared with individual algorithm implementations on bulk RNA-seq and microarray data. In an effort to extend this approach to scRNA-seq data, we present COFFEE (COnsensus single cell-type speciFic inFerence for gEnE regulatory networks), a Borda voting-based consensus algorithm that integrates information from 10 established GRN inference methods. We conclude that COFFEE has improved performance across synthetic, curated, and experimental datasets when compared with baseline methods. Additionally, we show that a modified version of COFFEE can be leveraged to improve performance on newer cell-type specific GRN inference methods. Overall, our results demonstrate that consensus-based methods with pertinent modifications continue to be valuable for GRN inference at the single cell level. While COFFEE is benchmarked on 10 algorithms, it is a flexible strategy that can incorporate any set of GRN inference algorithms according to user preference. A Python implementation of COFFEE may be found on GitHub: https://github.com/lodimk2/coffee.

A primate-specific endogenous retroviral envelope protein sequesters SFRP2 to regulate human cardiomyocyte development

Endogenous retroviruses (ERVs) occupy a significant part of the human genome, with some encoding proteins that influence the immune system or regulate cell-cell fusion in early extra-embryonic development. However, whether ERV-derived proteins regulate somatic development is unknown. Here, we report a somatic developmental function for the primate-specific ERVH48-1 (SUPYN/Suppressyn). ERVH48-1 encodes a fragment of a viral envelope that is expressed during early embryonic development. Loss of ERVH48-1 led to impaired mesoderm and cardiomyocyte commitment and diverted cells to an ectoderm-like fate. Mechanistically, ERVH48-1 is localized to sub-cellular membrane compartments through a functional N-terminal signal peptide and binds to the WNT antagonist SFRP2 to promote its polyubiquitination and degradation, thus limiting SFRP2 secretion and blocking repression of WNT/β-catenin signaling. Knockdown of SFRP2 or expression of a chimeric SFRP2 with the ERVH48-1 signal peptide rescued cardiomyocyte differentiation. This study demonstrates how ERVH48-1 modulates WNT/β-catenin signaling and cell type commitment in somatic development.

Inference of single-cell network using mutual information for scRNA-seq data analysis

BACKGROUND

With the advance in single-cell RNA sequencing (scRNA-seq) technology, deriving inherent biological system information from expression profiles at a single-cell resolution has become possible. It has been known that network modeling by estimating the associations between genes could better reveal dynamic changes in biological systems. However, accurately constructing a single-cell network (SCN) to capture the network architecture of each cell and further explore cell-to-cell heterogeneity remains challenging.

RESULTS

We introduce SINUM, a method for constructing the SIngle-cell Network Using Mutual information, which estimates mutual information between any two genes from scRNA-seq data to determine whether they are dependent or independent in a specific cell. Experiments on various scRNA-seq datasets with different cell numbers based on eight performance indexes (e.g., adjusted rand index and F-measure index) validated the accuracy and robustness of SINUM in cell type identification, superior to the state-of-the-art SCN inference method. Additionally, the SINUM SCNs exhibit high overlap with the human interactome and possess the scale-free property.

CONCLUSIONS

SINUM presents a view of biological systems at the network level to detect cell-type marker genes/gene pairs and investigate time-dependent changes in gene associations during embryo development. Codes for SINUM are freely available at https://github.com/SysMednet/SINUM .

Genome wide clustering on integrated chromatin states and Micro-C contacts reveals chromatin interaction signatures

We can now analyze 3D physical interactions of chromatin regions with chromatin conformation capture technologies, in addition to the 1D chromatin state annotations, but methods to integrate this information are lacking. We propose a method to integrate the chromatin state of interacting regions into a vector representation through the contact-weighted sum of chromatin states. Unsupervised clustering on integrated chromatin states and Micro-C contacts reveals common patterns of chromatin interaction signatures. This provides an integrated view of the complex dynamics of concurrent change occurring in chromatin state and in chromatin interaction, adding another layer of annotation beyond chromatin state or Hi-C contact separately.

Collaborative Structure-Preserved Missing Data Imputation for Single-Cell RNA-Seq Clustering

Clustering of the single-cell RNA-seq (scRNA-seq) transcriptome profiles is able to identify cell types, which is beneficial to improve the understanding of disease progression. However, in practice, the single-cell expression data often contains a significant number of missing values as a result of technical variability. Missing data is a critical challenge in scRNA-seq clustering analysis since the unknown value does not reflect the underlying true expression level and makes it difficult to discovering cell types by applying clustering algorithms directly. Various approaches have been developed to overcome missing data issue in scRNA-seq clustering. Most of them recover missing expression values by borrowing observed data from similar cells or synthesizing data via generative adversarial networks. Such that the biologically meaningful cluster structure has not been sufficiently exploited. In this work, we introduce ColImpute, a collaborative structure-preserved missing data imputation approach for the scRNA-seq clustering. Specifically, a cluster structure-preserved imputation module and a subspace clustering module, which respectively perform missing data imputation and cell subtypes identification, are integrated into a unified optimization framework to train the two networks in a collaborative manner. Consequently, the clustering module effectively contributes cluster-structure information to guide the trainning process of the missing data imputation module. Simultaneously, the cluster structure-preserved imputation module reciprocally enhances the performance of the clustering module by generating more precise recovered samples. Promising experimental results show that the proposed method is effective for both the data imputation and the cell types identification.

EPHA2 is a novel cell surface marker of OCT4-positive undifferentiated cells during the differentiation of mouse and human pluripotent stem cells

Embryonic stem cells (ESCs) and induced pluripotent stem cells (iPSCs) possess the intrinsic ability to differentiate into diverse cellular lineages, marking them as potent instruments in regenerative medicine. Nonetheless, the proclivity of these stem cells to generate teratomas post-transplantation presents a formidable obstacle to their therapeutic utility. In previous studies, we identified an array of cell surface proteins specifically expressed in the pluripotent state, as revealed through proteomic analysis. Here we focused on EPHA2, a protein found to be abundantly present on the surface of undifferentiated mouse ESCs and is diminished upon differentiation. Knock-down of Epha2 led to the spontaneous differentiation of mouse ESCs, underscoring a pivotal role of EPHA2 in maintaining an undifferentiated cell state. Further investigations revealed a strong correlation between EPHA2 and OCT4 expression during the differentiation of both mouse and human PSCs. Notably, removing EPHA2+ cells from mouse ESC-derived hepatic lineage reduced tumor formation after transplanting them into immune-deficient mice. Similarly, in human iPSCs, a larger proportion of EPHA2+ cells correlated with higher OCT4 expression, reflecting the pattern observed in mouse ESCs. Conclusively, EPHA2 emerges as a potential marker for selecting undifferentiated stem cells, providing a valuable method to decrease tumorigenesis risks after stem-cell transplantation in regenerative treatments.

Analyzing Single Cell RNA Sequencing with Topological Nonnegative Matrix Factorization

Single-cell RNA sequencing (scRNA-seq) is a relatively new technology that has stimulated enormous interest in statistics, data science, and computational biology due to the high dimensionality, complexity, and large scale associated with scRNA-seq data. Nonnegative matrix factorization (NMF) offers a unique approach due to its meta-gene interpretation of resulting low-dimensional components. However, NMF approaches suffer from the lack of multiscale analysis. This work introduces two persistent Laplacian regularized NMF methods, namely, topological NMF (TNMF) and robust topological NMF (rTNMF). By employing a total of 12 datasets, we demonstrate that the proposed TNMF and rTNMF significantly outperform all other NMF-based methods. We have also utilized TNMF and rTNMF for the visualization of popular Uniform Manifold Approximation and Projection (UMAP) and t -distributed stochastic neighbor embedding (t -SNE).

Graph-Regularized Non-Negative Matrix Factorization for Single-Cell Clustering in scRNA-Seq Data

The advent of single-cell RNA sequencing (scRNA-seq) has brought forth fresh perspectives on intricate biological processes, revealing the nuances and divergences present among distinct cells. Accurate single-cell analysis is a crucial prerequisite for in-depth investigation into the underlying mechanisms of heterogeneity. Due to various technical noises, like the impact of dropout values, scRNA-seq data remains challenging to interpret. In this work, we propose an unsupervised learning framework for scRNA-seq data analysis (aka Sc-GNNMF). Based on the non-negativity and sparsity of scRNA-seq data, we propose employing graph-regularized non-negative matrix factorization (GNNMF) algorithm for the analysis of scRNA-seq data, which involves estimating cell-cell sparse similarity and gene-gene sparse similarity through Laplacian kernels and p-nearest neighbor graphs ( p-NNG). By assuming intrinsic geometric local invariance, we use a weighted p-nearest known neighbors ( p-NKN) to optimize the scRNA-seq data. The optimized scRNA-seq data then participates in the matrix decomposition process, promoting the closeness of cells with similar types in cell-gene data space and determining a more suitable embedding space for clustering. Sc-GNNMF demonstrates superior performance compared to other methods and maintains satisfactory compatibility and robustness, as evidenced by experiments on 11 real scRNA-seq datasets. Furthermore, Sc-GNNMF yields excellent results in clustering tasks, extracting useful gene markers, and pseudo-temporal analysis.

Clustering single-cell RNA sequencing data via iterative smoothing and self-supervised discriminative embedding

Single-cell transcriptome sequencing (scRNA-seq) is a high-throughput technique used to study gene expression at the single-cell level. Clustering analysis is a commonly used method in scRNA-seq data analysis, helping researchers identify cell types and uncover interactions between cells. However, the choice of a robust similarity metric in the clustering procedure is still an open challenge due to the complex underlying structures of the data and the inherent noise in data acquisition. Here, we propose a deep clustering method for scRNA-seq data called scRISE (scRNA-seq Iterative Smoothing and self-supervised discriminative Embedding model) to resolve this challenge. The model consists of two main modules: an iterative smoothing module based on graph autoencoders designed to denoise the data and refine the pairwise similarity in turn to gradually incorporate cell structural features and enrich the data information; and a self-supervised discriminative embedding module with adaptive similarity threshold for partitioning samples into correct clusters. Our approach has shown improved quality of data representation and clustering on seventeen scRNA-seq datasets against a number of state-of-the-art deep learning clustering methods. Furthermore, utilizing the scRISE method in biological analysis against the HNSCC dataset has unveiled 62 informative genes, highlighting their potential roles as therapeutic targets and biomarkers.

SFINN: inferring gene regulatory network from single-cell and spatial transcriptomic data with shared factor neighborhood and integrated neural network

MOTIVATION

The rise of single-cell RNA sequencing (scRNA-seq) technology presents new opportunities for constructing detailed cell type-specific gene regulatory networks (GRNs) to study cell heterogeneity. However, challenges caused by noises, technical errors, and dropout phenomena in scRNA-seq data pose significant obstacles to GRN inference, making the design of accurate GRN inference algorithms still essential. The recent growth of both single-cell and spatial transcriptomic sequencing data enables the development of supervised deep learning methods to infer GRNs on these diverse single-cell datasets.

RESULTS

In this study, we introduce a novel deep learning framework based on shared factor neighborhood and integrated neural network (SFINN) for inferring potential interactions and causalities between transcription factors and target genes from single-cell and spatial transcriptomic data. SFINN utilizes shared factor neighborhood to construct cellular neighborhood network based on gene expression data and additionally integrates cellular network generated from spatial location information. Subsequently, the cell adjacency matrix and gene pair expression are fed into an integrated neural network framework consisting of a graph convolutional neural network and a fully-connected neural network to determine whether the genes interact. Performance evaluation in the tasks of gene interaction and causality prediction against the existing GRN reconstruction algorithms demonstrates the usability and competitiveness of SFINN across different kinds of data. SFINN can be applied to infer GRNs from conventional single-cell sequencing data and spatial transcriptomic data.

AVAILABILITY AND IMPLEMENTATION

SFINN can be accessed at GitHub: https://github.com/JGuan-lab/SFINN.

DeepIMAGER: Deeply Analyzing Gene Regulatory Networks from scRNA-seq Data

Understanding the dynamics of gene regulatory networks (GRNs) across diverse cell types poses a challenge yet holds immense value in unraveling the molecular mechanisms governing cellular processes. Current computational methods, which rely solely on expression changes from bulk RNA-seq and/or scRNA-seq data, often result in high rates of false positives and low precision. Here, we introduce an advanced computational tool, DeepIMAGER, for inferring cell-specific GRNs through deep learning and data integration. DeepIMAGER employs a supervised approach that transforms the co-expression patterns of gene pairs into image-like representations and leverages transcription factor (TF) binding information for model training. It is trained using comprehensive datasets that encompass scRNA-seq profiles and ChIP-seq data, capturing TF-gene pair information across various cell types. Comprehensive validations on six cell lines show DeepIMAGER exhibits superior performance in ten popular GRN inference tools and has remarkable robustness against dropout-zero events. DeepIMAGER was applied to scRNA-seq datasets of multiple myeloma (MM) and detected potential GRNs for TFs of RORC, MITF, and FOXD2 in MM dendritic cells. This technical innovation, combined with its capability to accurately decode GRNs from scRNA-seq, establishes DeepIMAGER as a valuable tool for unraveling complex regulatory networks in various cell types.

CPMI: comprehensive neighborhood-based perturbed mutual information for identifying critical states of complex biological processes

BACKGROUND

There exists a critical transition or tipping point during the complex biological process. Such critical transition is usually accompanied by the catastrophic consequences. Therefore, hunting for the tipping point or critical state is of significant importance to prevent or delay the occurrence of catastrophic consequences. However, predicting critical state based on the high-dimensional small sample data is a difficult problem, especially for single-cell expression data.

RESULTS

In this study, we propose the comprehensive neighbourhood-based perturbed mutual information (CPMI) method to detect the critical states of complex biological processes. The CPMI method takes into account the relationship between genes and neighbours, so as to reduce the noise and enhance the robustness. This method is applied to a simulated dataset and six real datasets, including an influenza dataset, two single-cell expression datasets and three bulk datasets. The method can not only successfully detect the tipping points, but also identify their dynamic network biomarkers (DNBs). In addition, the discovery of transcription factors (TFs) which can regulate DNB genes and nondifferential 'dark genes' validates the effectiveness of our method. The numerical simulation verifies that the CPMI method is robust under different noise strengths and is superior to the existing methods on identifying the critical states.

CONCLUSIONS

In conclusion, we propose a robust computational method, i.e., CPMI, which is applicable in both the bulk and single cell datasets. The CPMI method holds great potential in providing the early warning signals for complex biological processes and enabling early disease diagnosis.

scDMAE: A Generative Denoising Model Adopted Mask Strategy for scRNA-Seq Data Recovery

The advent of single-cell RNA sequencing (scRNA-seq) technology has revolutionized gene expression studies at the single-cell level. However, the presence of technical noise and data sparsity in scRNA-seq often undermines the accuracy of subsequent analyses. Existing methods for denoising and imputing scRNA-seq data often rely on stringent assumptions about data distribution, limiting the effectiveness of data recovery. In this study, we propose the scDMAE model for denoising and recovery of scRNA-seq data. First, the model fuses gene expression features and topological features to discern the primary expression patterns of genes in cells. Then, an autoencoder with a masking strategy is used to model dropout events and separate potential noise in the data. Finally, the model incorporates the original raw data to recover the true biological expression value. By conducting experiments on various types of scRNA-Seq datasets, scDMAE demonstrates superior performance compared to other comparative methods based on six distinct evaluation metrics in downstream analysis. The scDMAE method can accurately cluster similar cell populations, identify differential genes and infer cell trajectories.

Unveiling gene regulatory networks during cellular state transitions without linkage across time points

Time-stamped cross-sectional data, which lack linkage across time points, are commonly generated in single-cell transcriptional profiling. Many previous methods for inferring gene regulatory networks (GRNs) driving cell-state transitions relied on constructing single-cell temporal ordering. Introducing COSLIR (COvariance restricted Sparse LInear Regression), we presented a direct approach to reconstructing GRNs that govern cell-state transitions, utilizing only the first and second moments of samples between two consecutive time points. Simulations validated COSLIR's perfect accuracy in the oracle case and demonstrated its robust performance in real-world scenarios. When applied to single-cell RT-PCR and RNAseq datasets in developmental biology, COSLIR competed favorably with existing methods. Notably, its running time remained nearly independent of the number of cells. Therefore, COSLIR emerges as a promising addition to GRN reconstruction methods under cell-state transitions, bypassing the single-cell temporal ordering to enhance accuracy and efficiency in single-cell transcriptional profiling.

Inferring gene regulatory networks from single-cell transcriptomics based on graph embedding

MOTIVATION

Gene regulatory networks (GRNs) encode gene regulation in living organisms, and have become a critical tool to understand complex biological processes. However, due to the dynamic and complex nature of gene regulation, inferring GRNs from scRNA-seq data is still a challenging task. Existing computational methods usually focus on the close connections between genes, and ignore the global structure and distal regulatory relationships.

RESULTS

In this study, we develop a supervised deep learning framework, IGEGRNS, to infer GRNs from scRNA-seq data based on graph embedding. In the framework, contextual information of genes is captured by GraphSAGE, which aggregates gene features and neighborhood structures to generate low-dimensional embedding for genes. Then, the k most influential nodes in the whole graph are filtered through Top-k pooling. Finally, potential regulatory relationships between genes are predicted by stacking CNNs. Compared with nine competing supervised and unsupervised methods, our method achieves better performance on six time-series scRNA-seq datasets.

AVAILABILITY AND IMPLEMENTATION

Our method IGEGRNS is implemented in Python using the Pytorch machine learning library, and it is freely available at https://github.com/DHUDBlab/IGEGRNS.

Sex-biased gene expression during neural differentiation of human embryonic stem cells

Sex differences in the developing human brain are primarily attributed to hormonal influence. Recently however, genetic differences and their impact on the developing nervous system have attracted increased attention. To understand genetically driven sexual dimorphisms in neurodevelopment, we investigated genome-wide gene expression in an in vitro differentiation model of male and female human embryonic stem cell lines (hESC), independent of the effects of human sex hormones. Four male and four female-derived hESC lines were differentiated into a population of mixed neurons over 37 days. Differential gene expression and gene set enrichment analyses were conducted on bulk RNA sequencing data. While similar differentiation tendencies in all cell lines demonstrated the robustness and reproducibility of our differentiation protocol, we found sex-biased gene expression already in undifferentiated ESCs at day 0, but most profoundly after 37 days of differentiation. Male and female cell lines exhibited sex-biased expression of genes involved in neurodevelopment, suggesting that sex influences the differentiation trajectory. Interestingly, the highest contribution to sex differences was found to arise from the male transcriptome, involving both Y chromosome and autosomal genes. We propose 13 sex-biased candidate genes (10 upregulated in male cell lines and 3 in female lines) that are likely to affect neuronal development. Additionally, we confirmed gene dosage compensation of X/Y homologs escaping X chromosome inactivation through their Y homologs and identified a significant overexpression of the Y-linked demethylase UTY and KDM5D in male hESC during neuron development, confirming previous results in neural stem cells. Our results suggest that genetic sex differences affect neuronal differentiation trajectories, which could ultimately contribute to sex biases during human brain development.

Preprocessing of Single Cell RNA Sequencing Data Using Correlated Clustering and Projection

Single-cell RNA sequencing (scRNA-seq) is widely used to reveal heterogeneity in cells, which has given us insights into cell-cell communication, cell differentiation, and differential gene expression. However, analyzing scRNA-seq data is a challenge due to sparsity and the large number of genes involved. Therefore, dimensionality reduction and feature selection are important for removing spurious signals and enhancing the downstream analysis. We present Correlated Clustering and Projection (CCP), a new data-domain dimensionality reduction method, for the first time. CCP projects each cluster of similar genes into a supergene defined as the accumulated pairwise nonlinear gene-gene correlations among all cells. Using 14 benchmark data sets, we demonstrate that CCP has significant advantages over classical principal component analysis (PCA) for clustering and/or classification problems with intrinsically high dimensionality. In addition, we introduce the Residue-Similarity index (RSI) as a novel metric for clustering and classification and the R-S plot as a new visualization tool. We show that the RSI correlates with accuracy without requiring the knowledge of the true labels. The R-S plot provides a unique alternative to the uniform manifold approximation and projection (UMAP) and t-distributed stochastic neighbor embedding (t-SNE) for data with a large number of cell types.

A patterned human neural tube model using microfluidic gradients

The human nervous system is a highly complex but organized organ. The foundation of its complexity and organization is laid down during regional patterning of the neural tube, the embryonic precursor to the human nervous system. Historically, studies of neural tube patterning have relied on animal models to uncover underlying principles. Recently, models of neurodevelopment based on human pluripotent stem cells, including neural organoids1-5 and bioengineered neural tube development models6-10, have emerged. However, such models fail to recapitulate neural patterning along both rostral-caudal and dorsal-ventral axes in a three-dimensional tubular geometry, a hallmark of neural tube development. Here we report a human pluripotent stem cell-based, microfluidic neural tube-like structure, the development of which recapitulates several crucial aspects of neural patterning in brain and spinal cord regions and along rostral-caudal and dorsal-ventral axes. This structure was utilized for studying neuronal lineage development, which revealed pre-patterning of axial identities of neural crest progenitors and functional roles of neuromesodermal progenitors and the caudal gene CDX2 in spinal cord and trunk neural crest development. We further developed dorsal-ventral patterned microfluidic forebrain-like structures with spatially segregated dorsal and ventral regions and layered apicobasal cellular organizations that mimic development of the human forebrain pallium and subpallium, respectively. Together, these microfluidics-based neurodevelopment models provide three-dimensional lumenal tissue architectures with in vivo-like spatiotemporal cell differentiation and organization, which will facilitate the study of human neurodevelopment and disease.

DELVE: feature selection for preserving biological trajectories in single-cell data

Single-cell technologies can measure the expression of thousands of molecular features in individual cells undergoing dynamic biological processes. While examining cells along a computationally-ordered pseudotime trajectory can reveal how changes in gene or protein expression impact cell fate, identifying such dynamic features is challenging due to the inherent noise in single-cell data. Here, we present DELVE, an unsupervised feature selection method for identifying a representative subset of molecular features which robustly recapitulate cellular trajectories. In contrast to previous work, DELVE uses a bottom-up approach to mitigate the effects of confounding sources of variation, and instead models cell states from dynamic gene or protein modules based on core regulatory complexes. Using simulations, single-cell RNA sequencing, and iterative immunofluorescence imaging data in the context of cell cycle and cellular differentiation, we demonstrate how DELVE selects features that better define cell-types and cell-type transitions. DELVE is available as an open-source python package: https://github.com/jranek/delve .

Topological and geometric analysis of cell states in single-cell transcriptomic data

Single-cell RNA sequencing (scRNA-seq) enables dissecting cellular heterogeneity in tissues, resulting in numerous biological discoveries. Various computational methods have been devised to delineate cell types by clustering scRNA-seq data, where clusters are often annotated using prior knowledge of marker genes. In addition to identifying pure cell types, several methods have been developed to identify cells undergoing state transitions, which often rely on prior clustering results. The present computational approaches predominantly investigate the local and first-order structures of scRNA-seq data using graph representations, while scRNA-seq data frequently display complex high-dimensional structures. Here, we introduce scGeom, a tool that exploits the multiscale and multidimensional structures in scRNA-seq data by analyzing the geometry and topology through curvature and persistent homology of both cell and gene networks. We demonstrate the utility of these structural features to reflect biological properties and functions in several applications, where we show that curvatures and topological signatures of cell and gene networks can help indicate transition cells and the differentiation potential of cells. We also illustrate that structural characteristics can improve the classification of cell types.

Gastrulation-stage gene expression in Nipbl+/- mouse embryos foreshadows the development of syndromic birth defects

In animal models, Nipbl deficiency phenocopies gene expression changes and birth defects seen in Cornelia de Lange syndrome, the most common cause of which is Nipbl haploinsufficiency. Previous studies in Nipbl+/- mice suggested that heart development is abnormal as soon as cardiogenic tissue is formed. To investigate this, we performed single-cell RNA sequencing on wild-type and Nipbl+/- mouse embryos at gastrulation and early cardiac crescent stages. Nipbl+/- embryos had fewer mesoderm cells than wild-type and altered proportions of mesodermal cell subpopulations. These findings were associated with underexpression of genes implicated in driving specific mesodermal lineages. In addition, Nanog was found to be overexpressed in all germ layers, and many gene expression changes observed in Nipbl+/- embryos could be attributed to Nanog overexpression. These findings establish a link between Nipbl deficiency, Nanog overexpression, and gene expression dysregulation/lineage misallocation, which ultimately manifest as birth defects in Nipbl+/- animals and Cornelia de Lange syndrome.

DANCE: a deep learning library and benchmark platform for single-cell analysis

DANCE is the first standard, generic, and extensible benchmark platform for accessing and evaluating computational methods across the spectrum of benchmark datasets for numerous single-cell analysis tasks. Currently, DANCE supports 3 modules and 8 popular tasks with 32 state-of-art methods on 21 benchmark datasets. People can easily reproduce the results of supported algorithms across major benchmark datasets via minimal efforts, such as using only one command line. In addition, DANCE provides an ecosystem of deep learning architectures and tools for researchers to facilitate their own model development. DANCE is an open-source Python package that welcomes all kinds of contributions.

Charting cellular differentiation trajectories with Ricci flow

Complex biological processes, such as cellular differentiation, require intricate rewiring of intra-cellular signalling networks. Previous characterisations revealed a raised network entropy underlies less differentiated and malignant cell states. A connection between entropy and Ricci curvature led to applications of discrete curvatures to biological networks. However, predicting dynamic biological network rewiring remains an open problem. Here we apply Ricci curvature and Ricci flow to biological network rewiring. By investigating the relationship between network entropy and Forman-Ricci curvature, theoretically and empirically on single-cell RNA-sequencing data, we demonstrate that the two measures do not always positively correlate, as previously suggested, and provide complementary rather than interchangeable information. We next employ Ricci flow to derive network rewiring trajectories from stem cells to differentiated cells, accurately predicting true intermediate time points in gene expression time courses. In summary, we present a differential geometry toolkit for understanding dynamic network rewiring during cellular differentiation and cancer.

Generation of human excitatory forebrain neurons by cooperative binding of proneural NGN2 and homeobox factor EMX1

Generation of defined neuronal subtypes from human pluripotent stem cells remains a challenge. The proneural factor NGN2 has been shown to overcome experimental variability observed by morphogen-guided differentiation and directly converts pluripotent stem cells into neurons, but their cellular heterogeneity has not been investigated yet. Here, we found that NGN2 reproducibly produces three different kinds of excitatory neurons characterized by partial coactivation of other neurotransmitter programs. We explored two principle approaches to achieve more precise specification: prepatterning the chromatin landscape that NGN2 is exposed to and combining NGN2 with region-specific transcription factors. Unexpectedly, the chromatin context of regionalized neural progenitors only mildly altered genomic NGN2 binding and its transcriptional response and did not affect neurotransmitter specification. In contrast, coexpression of region-specific homeobox factors such as EMX1 resulted in drastic redistribution of NGN2 including recruitment to homeobox targets and resulted in glutamatergic neurons with silenced nonglutamatergic programs. These results provide the molecular basis for a blueprint for improved strategies for generating a plethora of defined neuronal subpopulations from pluripotent stem cells for therapeutic or disease-modeling purposes.

SinCWIm: An imputation method for single-cell RNA sequence dropouts using weighted alternating least squares

BACKGROUND AND OBJECTIVES

Single-cell RNA sequencing (scRNA-seq) provides a powerful tool for exploring cellular heterogeneity, discovering novel or rare cell types, distinguishing between tissue-specific cellular composition, and understanding cell differentiation during development. However, due to technological limitations, dropout events in scRNA-seq can mistakenly convert some entries in the real data to zero. This is equivalent to introducing noise into the data of cell gene expression entries. The data is contaminated, which affects the performance of downstream analyses, including clustering, cell annotation, differential gene expression analysis, and so on. Therefore, it is a crucial work to accurately determine which zeros are due to dropout events and perform imputation operations on them.

METHODS

Considering the different confidence levels of different zeros in the gene expression matrix, this paper proposes a SinCWIm method for dropout events in scRNA-seq based on weighted alternating least squares (WALS). The method utilizes Pearson correlation coefficient and hierarchical clustering to quantify the confidence of zero entries. It is then combined with WALS for matrix decomposition. And the imputation result is made close to the actual number by outlier removal and data correction operations.

RESULTS

A total of eight single-cell sequencing datasets were used for comparative experiments to demonstrate the overall superiority of SinCWIm over state-of-the-art models. SinCWIm was applied to cluster the data to obtain an adjusted RAND index evaluation, and the Usoskin, Pollen and Bladder datasets scored 94.46%, 96.48% and 76.74%, respectively. In addition, significant improvements were made in the retention of differential expression genes and visualization.

CONCLUSIONS

SinCWIm provides a valuable imputation method for handling dropout events in single-cell sequencing data. In comparison to advanced methods, SinCWIm demonstrates excellent performance in clustering, visualization and other aspects. It is applicable to various single-cell sequencing datasets.

Multiscale differential geometry learning of networks with applications to single-cell RNA sequencing data

Single-cell RNA sequencing (scRNA-seq) has emerged as a transformative technology, offering unparalleled insights into the intricate landscape of cellular diversity and gene expression dynamics. scRNA-seq analysis represents a challenging and cutting-edge frontier within the field of biological research. Differential geometry serves as a powerful mathematical tool in various applications of scientific research. In this study, we introduce, for the first time, a multiscale differential geometry (MDG) strategy for addressing the challenges encountered in scRNA-seq data analysis. We assume that intrinsic properties of cells lie on a family of low-dimensional manifolds embedded in the high-dimensional space of scRNA-seq data. Multiscale cell-cell interactive manifolds are constructed to reveal complex relationships in the cell-cell network, where curvature-based features for cells can decipher the intricate structural and biological information. We showcase the utility of our novel approach by demonstrating its effectiveness in classifying cell types. This innovative application of differential geometry in scRNA-seq analysis opens new avenues for understanding the intricacies of biological networks and holds great potential for network analysis in other fields.

TFvelo: gene regulation inspired RNA velocity estimation

RNA velocity is closely related with cell fate and is an important indicator for the prediction of cell states with elegant physical explanation derived from single-cell RNA-seq data. Most existing RNA velocity models aim to extract dynamics from the phase delay between unspliced and spliced mRNA for each individual gene. However, unspliced/spliced mRNA abundance may not provide sufficient signal for dynamic modeling, leading to poor fit in phase portraits. Motivated by the idea that RNA velocity could be driven by the transcriptional regulation, we propose TFvelo, which expands RNA velocity concept to various single-cell datasets without relying on splicing information, by introducing gene regulatory information. Our experiments on synthetic data and multiple scRNA-Seq datasets show that TFvelo can accurately fit genes dynamics on phase portraits, and effectively infer cell pseudo-time and trajectory from RNA abundance data. TFvelo opens a robust and accurate avenue for modeling RNA velocity for single cell data.

Data integration and inference of gene regulation using single-cell temporal multimodal data with scTIE

Single-cell technologies offer unprecedented opportunities to dissect gene regulatory mechanisms in context-specific ways. Although there are computational methods for extracting gene regulatory relationships from scRNA-seq and scATAC-seq data, the data integration problem, essential for accurate cell type identification, has been mostly treated as a standalone challenge. Here we present scTIE, a unified method that integrates temporal multimodal data and infers regulatory relationships predictive of cellular state changes. scTIE uses an autoencoder to embed cells from all time points into a common space by using iterative optimal transport, followed by extracting interpretable information to predict cell trajectories. Using a variety of synthetic and real temporal multimodal data sets, we show scTIE achieves effective data integration while preserving more biological signals than existing methods, particularly in the presence of batch effects and noise. Furthermore, on the exemplar multiome data set we generated from differentiating mouse embryonic stem cells over time, we show scTIE captures regulatory elements highly predictive of cell transition probabilities, providing new potentials to understand the regulatory landscape driving developmental processes.

FAVA: high-quality functional association networks inferred from scRNA-seq and proteomics data

MOTIVATION

Protein networks are commonly used for understanding how proteins interact. However, they are typically biased by data availability, favoring well-studied proteins with more interactions. To uncover functions of understudied proteins, we must use data that are not affected by this literature bias, such as single-cell RNA-seq and proteomics. Due to data sparseness and redundancy, functional association analysis becomes complex.

RESULTS

To address this, we have developed FAVA (Functional Associations using Variational Autoencoders), which compresses high-dimensional data into a low-dimensional space. FAVA infers networks from high-dimensional omics data with much higher accuracy than existing methods, across a diverse collection of real as well as simulated datasets. FAVA can process large datasets with over 0.5 million conditions and has predicted 4210 interactions between 1039 understudied proteins. Our findings showcase FAVA's capability to offer novel perspectives on protein interactions. FAVA functions within the scverse ecosystem, employing AnnData as its input source.

AVAILABILITY AND IMPLEMENTATION

Source code, documentation, and tutorials for FAVA are accessible on GitHub at https://github.com/mikelkou/fava. FAVA can also be installed and used via pip/PyPI as well as via the scverse ecosystem https://github.com/scverse/ecosystem-packages/tree/main/packages/favapy.

Hypoblast from human pluripotent stem cells regulates epiblast development

Recently, several studies using cultures of human embryos together with single-cell RNA-seq analyses have revealed differences between humans and mice, necessitating the study of human embryos1-8. Despite the importance of human embryology, ethical and legal restrictions have limited post-implantation-stage studies. Thus, recent efforts have focused on developing in vitro self-organizing models using human stem cells9-17. Here, we report genetic and non-genetic approaches to generate authentic hypoblast cells (naive hPSC-derived hypoblast-like cells (nHyCs))-known to give rise to one of the two extraembryonic tissues essential for embryonic development-from naive human pluripotent stem cells (hPSCs). Our nHyCs spontaneously assemble with naive hPSCs to form a three-dimensional bilaminar structure (bilaminoids) with a pro-amniotic-like cavity. In the presence of additional naive hPSC-derived analogues of the second extraembryonic tissue, the trophectoderm, the efficiency of bilaminoid formation increases from 20% to 40%, and the epiblast within the bilaminoids continues to develop in response to trophectoderm-secreted IL-6. Furthermore, we show that bilaminoids robustly recapitulate the patterning of the anterior-posterior axis and the formation of cells reflecting the pregastrula stage, the emergence of which can be shaped by genetically manipulating the DKK1/OTX2 hypoblast-like domain. We have therefore successfully modelled and identified the mechanisms by which the two extraembryonic tissues efficiently guide the stage-specific growth and progression of the epiblast as it establishes the post-implantation landmarks of human embryogenesis.

Bacterial phenotypic heterogeneity through the lens of single-cell RNA sequencing

Bacterial transcription is not monolithic. Microbes exist in a wide variety of cell states that help them adapt to their environment, acquire and produce essential nutrients, and engage in both competition and cooperation with their neighbors. While we typically think of bacterial adaptation as a group behavior, where all cells respond in unison, there is often a mixture of phenotypic responses within a bacterial population, where distinct cell types arise. A primary phenomenon driving these distinct cell states is transcriptional heterogeneity. Given that bacterial mRNA transcripts are extremely short-lived compared to eukaryotes, their transcriptional state is closely associated with their physiology, and thus the transcriptome of a bacterial cell acts as a snapshot of the behavior of that bacterium. Therefore, the application of single-cell transcriptomics to microbial populations will provide novel insight into cellular differentiation and bacterial ecology. In this review, we provide an overview of transcriptional heterogeneity in microbial systems, discuss the findings already provided by single-cell approaches, and plot new avenues of inquiry in transcriptional regulation, cellular biology, and mechanisms of heterogeneity that are made possible when microbial communities are analyzed at single-cell resolution.

AGImpute: imputation of scRNA-seq data based on a hybrid GAN with dropouts identification

MOTIVATION

Dropout events bring challenges in analyzing single-cell RNA sequencing data as they introduce noise and distort the true distributions of gene expression profiles. Recent studies focus on estimating dropout probability and imputing dropout events by leveraging information from similar cells or genes. However, the number of dropout events differs in different cells, due to the complex factors, such as different sequencing protocols, cell types, and batch effects. The dropout event differences are not fully considered in assessing the similarities between cells and genes, which compromises the reliability of downstream analysis.

RESULTS

This work proposes a hybrid Generative Adversarial Network with dropouts identification to impute single-cell RNA sequencing data, named AGImpute. First, the numbers of dropout events in different cells in scRNA-seq data are differentially estimated by using a dynamic threshold estimation strategy. Next, the identified dropout events are imputed by a hybrid deep learning model, combining Autoencoder with a Generative Adversarial Network. To validate the efficiency of the AGImpute, it is compared with seven state-of-the-art dropout imputation methods on two simulated datasets and seven real single-cell RNA sequencing datasets. The results show that AGImpute imputes the least number of dropout events than other methods. Moreover, AGImpute enhances the performance of downstream analysis, including clustering performance, identifying cell-specific marker genes, and inferring trajectory in the time-course dataset.

AVAILABILITY AND IMPLEMENTATION

The source code can be obtained from https://github.com/xszhu-lab/AGImpute.

A Data-Distribution and Successive Spline Points based discretization approach for evolving gene regulatory networks from scRNA-Seq time-series data using Cartesian Genetic Programming

The inference of gene regulatory networks (GRNs) is a widely addressed problem in Systems Biology. GRNs can be modeled as Boolean networks, which is the simplest approach for this task. However, Boolean models need binarized data. Several approaches have been developed for the discretization of gene expression data (GED). Also, the advance of data extraction technologies, such as single-cell RNA-Sequencing (scRNA-Seq), provides a new vision of gene expression and brings new challenges for dealing with its specificities, such as a large occurrence of zero data. This work proposes a new discretization approach for dealing with scRNA-Seq time-series data, named Distribution and Successive Spline Points Discretization (DSSPD), which considers the data distribution and a proper preprocessing step. Here, Cartesian Genetic Programming (CGP) is used to infer GRNs using the results of DSSPD. The proposal is compared with CGP with the standard data handling and five state-of-the-art algorithms on curated models and experimental data. The results show that the proposal improves the results of CGP in all tested cases and outperforms the state-of-the-art algorithms in most cases.

JLONMFSC: Clustering scRNA-seq data based on joint learning of non-negative matrix factorization and subspace clustering

The development of single cell RNA sequencing (scRNA-seq) has provided new perspectives to study biological problems at the single cell level. One of the key issues in scRNA-seq data analysis is to divide cells into several clusters for discovering the heterogeneity and diversity of cells. However, the existing scRNA-seq data are high-dimensional, sparse, and noisy, which challenges the existing single-cell clustering methods. In this study, we propose a joint learning framework (JLONMFSC) for clustering scRNA-seq data. In our method, the dimension of the original data is reduced to minimize the effect of noise. In addition, the graph regularized matrix factorization is used to learn the local features. Further, the Low-Rank Representation (LRR) subspace clustering is utilized to learn the global features. Finally, the joint learning of local features and global features is performed to obtain the results of clustering. We compare the proposed algorithm with eight state-of-the-art algorithms for clustering performance on six datasets, and the experimental results demonstrate that the JLONMFSC achieves better performance in all datasets. The code is avalable at https://github.com/lanbiolab/JLONMFSC.

From local to global gene co-expression estimation using single-cell RNA-seq data

In genomics studies, the investigation of gene relationships often brings important biological insights. Currently, the large heterogeneous datasets impose new challenges for statisticians because gene relationships are often local. They change from one sample point to another, may only exist in a subset of the sample, and can be nonlinear or even nonmonotone. Most previous dependence measures do not specifically target local dependence relationships, and the ones that do are computationally costly. In this paper, we explore a state-of-the-art network estimation technique that characterizes gene relationships at the single cell level, under the name of cell-specific gene networks. We first show that averaging the cell-specific gene relationship over a population gives a novel univariate dependence measure, the averaged Local Density Gap (aLDG), that accumulates local dependence and can detect any nonlinear, nonmonotone relationship. Together with a consistent nonparametric estimator, we establish its robustness on both the population and empirical levels. Then, we show that averaging the cell-specific gene relationship over mini-batches determined by some external structure information (eg, spatial or temporal factor) better highlights meaningful local structure change points. We explore the application of aLDG and its minibatch variant in many scenarios, including pairwise gene relationship estimation, bifurcating point detection in cell trajectory, and spatial transcriptomics structure visualization. Both simulations and real data analysis show that aLDG outperforms existing ones.

GMFGRN: a matrix factorization and graph neural network approach for gene regulatory network inference

The recent advances of single-cell RNA sequencing (scRNA-seq) have enabled reliable profiling of gene expression at the single-cell level, providing opportunities for accurate inference of gene regulatory networks (GRNs) on scRNA-seq data. Most methods for inferring GRNs suffer from the inability to eliminate transitive interactions or necessitate expensive computational resources. To address these, we present a novel method, termed GMFGRN, for accurate graph neural network (GNN)-based GRN inference from scRNA-seq data. GMFGRN employs GNN for matrix factorization and learns representative embeddings for genes. For transcription factor-gene pairs, it utilizes the learned embeddings to determine whether they interact with each other. The extensive suite of benchmarking experiments encompassing eight static scRNA-seq datasets alongside several state-of-the-art methods demonstrated mean improvements of 1.9 and 2.5% over the runner-up in area under the receiver operating characteristic curve (AUROC) and area under the precision-recall curve (AUPRC). In addition, across four time-series datasets, maximum enhancements of 2.4 and 1.3% in AUROC and AUPRC were observed in comparison to the runner-up. Moreover, GMFGRN requires significantly less training time and memory consumption, with time and memory consumed <10% compared to the second-best method. These findings underscore the substantial potential of GMFGRN in the inference of GRNs. It is publicly available at https://github.com/Lishuoyy/GMFGRN.

scGIR: deciphering cellular heterogeneity via gene ranking in single-cell weighted gene correlation networks

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for investigating cellular heterogeneity through high-throughput analysis of individual cells. Nevertheless, challenges arise from prevalent sequencing dropout events and noise effects, impacting subsequent analyses. Here, we introduce a novel algorithm, Single-cell Gene Importance Ranking (scGIR), which utilizes a single-cell gene correlation network to evaluate gene importance. The algorithm transforms single-cell sequencing data into a robust gene correlation network through statistical independence, with correlation edges weighted by gene expression levels. We then constructed a random walk model on the resulting weighted gene correlation network to rank the importance of genes. Our analysis of gene importance using PageRank algorithm across nine authentic scRNA-seq datasets indicates that scGIR can effectively surmount technical noise, enabling the identification of cell types and inference of developmental trajectories. We demonstrated that the edges of gene correlation, weighted by expression, play a critical role in enhancing the algorithm's performance. Our findings emphasize that scGIR outperforms in enhancing the clustering of cell subtypes, reverse identifying differentially expressed marker genes, and uncovering genes with potential differential importance. Overall, we proposed a promising method capable of extracting more information from single-cell RNA sequencing datasets, potentially shedding new lights on cellular processes and disease mechanisms.