Computational methods for single-cell omics across modalities

The combined use of scRNA-seq and network propagation highlights key features of pan-cancer Tumor-Infiltrating T cells

Improving the selectivity and effectiveness of drugs represents a crucial issue for future therapeutic developments in immuno-oncology. Traditional bulk transcriptomics faces limitations in this context for the early phase of target discovery as resulting gene expression levels represent the average measure from multiple cell populations. Alternatively, single cell RNA sequencing can dive into unique cell populations transcriptome, facilitating the identification of specific targets. Here, we generated Tumor-Infiltrating regulatory T cells (TI-Tregs) and exhausted T cells (Tex) gene signatures from a single cell RNA-seq pan-cancer T cell atlas. To overcome noise and sparsity inherent to single cell transcriptomics, we then propagated the gene signatures by diffusion in a protein-protein interaction network using the Patrimony high-throughput computing platform. This methodology enabled the refining of signatures by rescoring genes based on their biological connectivity and shed light not only on processes characteristics of TI-Treg and Tex development and functions but also on their immunometabolic specificities. The combined use of single cell transcriptomics and network propagation may thus represent an innovative and effective methodology for the characterization of cell populations of interest and eventually the development of new therapeutic strategies in immuno-oncology.

Cut&tag: a powerful epigenetic tool for chromatin profiling

Analysis of transcription factors and chromatin modifications at the genome-wide level provides insights into gene regulatory processes, such as transcription, cell differentiation and cellular response. Chromatin immunoprecipitation is the most popular and powerful approach for mapping chromatin, and other enzyme-tethering techniques have recently become available for living cells. Among these, Cleavage Under Targets and Tagmentation (CUT&Tag) is a relatively novel chromatin profiling method that has rapidly gained popularity in the field of epigenetics since 2019. It has also been widely adapted to map chromatin modifications and TFs in different species, illustrating the association of these chromatin epitopes with various physiological and pathological processes. Scalable single-cell CUT&Tag can be combined with distinct platforms to distinguish cellular identity, epigenetic features and even spatial chromatin profiling. In addition, CUT&Tag has been developed as a strategy for joint profiling of the epigenome, transcriptome or proteome on the same sample. In this review, we will mainly consolidate the applications of CUT&Tag and its derivatives on different platforms, give a detailed explanation of the pros and cons of this technique as well as the potential development trends and applications in the future.

Addressing scalability and managing sparsity and dropout events in single-cell representation identification with ZIGACL

Despite significant advancements in single-cell representation learning, scalability and managing sparsity and dropout events continue to challenge the field as scRNA-seq datasets expand. While current computational tools struggle to maintain both efficiency and accuracy, the accurate connection of these dropout events to specific biological functions usually requires additional, complex experiments, often hampered by potential inaccuracies in cell-type annotation. To tackle these challenges, the Zero-Inflated Graph Attention Collaborative Learning (ZIGACL) method has been developed. This innovative approach combines a Zero-Inflated Negative Binomial model with a Graph Attention Network, leveraging mutual information from neighboring cells to enhance dimensionality reduction and apply dynamic adjustments to the learning process through a co-supervised deep graph clustering model. ZIGACL's integration of denoising and topological embedding significantly improves clustering accuracy and ensures similar cells are grouped closely in the latent space. Comparative analyses across nine real scRNA-seq datasets have shown that ZIGACL significantly enhances single-cell data analysis by offering superior clustering performance and improved stability in cell representations, effectively addressing scalability and managing sparsity and dropout events, thereby advancing our understanding of cellular heterogeneity.

TriTan: an efficient triple nonnegative matrix factorization method for integrative analysis of single-cell multiomics data

Single-cell multiomics have opened up tremendous opportunities for understanding gene regulatory networks underlying cell states by simultaneously profiling transcriptomes, epigenomes, and proteomes of the same cell. However, existing computational methods for integrative analysis of these high-dimensional multiomics data are either computationally expensive or limited in interpretation. These limitations pose challenges in the implementation of these methods in large-scale studies and hinder a more in-depth understanding of the underlying regulatory mechanisms. Here, we propose TriTan (Triple inTegrative fast non-negative matrix factorization), an efficient joint factorization method for single-cell multiomics data. TriTan implements a highly efficient factorization algorithm, greatly improving its computational performance. Three matrix factorization produced by TriTan helps in clustering cells, identifying signature features for each cell type, and uncovering feature associations across omics, which facilitates the identification of domains of regulatory chromatin and the prediction of cell-type-specific regulatory networks. We applied TriTan to the single-cell multiomics data obtained from different technologies and benchmarked it against the state-of-the-art methods where it shows highly competitive performance. Furthermore, we showed a range of downstream analyses conducted utilizing TriTan outputs, highlighting its capacity to facilitate interpretation in biological discovery.

Application of spatial-omics to the classification of kidney biopsy samples in transplantation

Improvement of long-term outcomes through targeted treatment is a primary concern in kidney transplant medicine. Currently, the validation of a rejection diagnosis and subsequent treatment depends on the histological assessment of allograft biopsy samples, according to the Banff classification system. However, the lack of (early) disease-specific tissue markers hinders accurate diagnosis and thus timely intervention. This challenge mainly results from an incomplete understanding of the pathophysiological processes underlying late allograft failure. Integration of large-scale multimodal approaches for investigating allograft biopsy samples might offer new insights into this pathophysiology, which are necessary for the identification of novel therapeutic targets and the development of tailored immunotherapeutic interventions. Several omics technologies - including transcriptomic, proteomic, lipidomic and metabolomic tools (and multimodal data analysis strategies) - can be applied to allograft biopsy investigation. However, despite their successful application in research settings and their potential clinical value, several barriers limit the broad implementation of many of these tools into clinical practice. Among spatial-omics technologies, mass spectrometry imaging, which is under-represented in the transplant field, has the potential to enable multi-omics investigations that might expand the insights gained with current clinical analysis technologies.

Proximogram-A multi-omics network-based framework to capture tissue heterogeneity integrating single-cell omics and spatial profiling

The increasing availability of patient-derived multimodal biological data for various diseases has opened up avenues for finding the optimal methods for jointly leveraging the information extracted in a customizable and scalable manner. Here, we propose the Proximogram, a graph-based representation that provides a joint construct for embedding independently obtained omics and spatial data. To evaluate the representation, we generated proximograms from 2 distinct biological sources, namely, multiplexed immunofluorescence images and single-cell RNA-seq data obtained from patients across two pancreatic diseases that include normal and chronic Pancreatitis (CP) and pancreatic ductal adenocarcinoma (PDAC). The generated proximograms were used as inputs to 2 distinct graph deep-learning models. The improved classification results over simpler spatial-data-based input graphs point to the increased discriminatory power obtained by integrating structural information from single-cell ligand-receptor signaling data and the spatial architecture of cells in each disease class, which can help point to markers of high diagnostic significance.

Profiling Dynamic Patterns of Single-Cell Motility

Cell motility plays an essential role in many biological processes as cells move and interact within their local microenvironments. Current methods for quantifying cell motility typically involve tracking individual cells over time, but the results are often presented as averaged values across cell populations. While informative, these ensemble approaches have limitations in assessing cellular heterogeneity and identifying generalizable patterns of single-cell behaviors, at baseline and in response to perturbations. In this study, CaMI is introduced, a computational framework designed to leverage the single-cell nature of motility data. CaMI identifies and classifies distinct spatio-temporal behaviors of individual cells, enabling robust classification of single-cell motility patterns in a large dataset (n = 74 253 cells). This framework allows quantification of spatial and temporal heterogeneities, determination of single-cell motility behaviors across various biological conditions and provides a visualization scheme for direct interpretation of dynamic cell behaviors. Importantly, CaMI reveals insights that conventional cell motility analyses may overlook, showcasing its utility in uncovering robust biological insights. Together, a multivariate framework is presented to classify emergent patterns of single-cell motility, emphasizing the critical role of cellular heterogeneity in shaping cell behaviors across populations.

Milletomics: a metabolomics centered integrated omics approach toward genetic progression

Producing alternative staple foods like millet will be essential to feeding ten billion people by 2050. The increased demand for millet is driving researchers to improve its genetic variation. Millets include protein, dietary fiber, phenolic substances, and flavonoid components. Its climate resilience makes millet an appealing crop for agronomic sustainability. Integrative omics technologies could potentially identify and develop millets with desirable phenotypes that may have high agronomic value. Millets' salinity and drought tolerance have been enhanced using transcriptomics. In foxtail, finger, and pearl millet, proteomics has discovered salt-tolerant protein, phytohormone-focused protein, and drought tolerance. Metabolomics studies have revealed that certain metabolic pathways including those involving lignin, flavonoids, phenylpropanoid, and lysophospholipids are critical for many processes, including seed germination, photosynthesis, energy metabolism, and the synthesis of bioactive chemicals necessary for drought tolerance. Metabolomics integration with other omics revealed metabolome engineering and trait-specific metabolite creation. Integrated metabolomics and ionomics are still in the development stage, but they could potentially assist in comprehending the pathway of ionomers to control nutrient levels and biofortify millet. Epigenomic analysis has shown alterations in DNA methylation patterns and chromatin structure in foxtail and pearl millets in response to abiotic stress. Whole-genome sequencing utilizing next-generation sequencing is the most proficient method for finding stress-induced phytoconstituent genes. New genome sequencing enables novel biotechnological interventions including genome-wide association, mutation-based research, and other omics approaches. Millets can breed more effectively by employing next-generation sequencing and genotyping by sequencing, which may mitigate climate change. Millet marker-assisted breeding has advanced with high-throughput markers and combined genotyping technologies.

A Cell Cycle-Aware Network for Data Integration and Label Transferring of Single-Cell RNA-Seq and ATAC-Seq

In recent years, the integration of single-cell multi-omics data has provided a more comprehensive understanding of cell functions and internal regulatory mechanisms from a non-single omics perspective, but it still suffers many challenges, such as omics-variance, sparsity, cell heterogeneity, and confounding factors. As it is known, the cell cycle is regarded as a confounder when analyzing other factors in single-cell RNA-seq data, but it is not clear how it will work on the integrated single-cell multi-omics data. Here, a cell cycle-aware network (CCAN) is developed to remove cell cycle effects from the integrated single-cell multi-omics data while keeping the cell type-specific variations. This is the first computational model to study the cell-cycle effects in the integration of single-cell multi-omics data. Validations on several benchmark datasets show the outstanding performance of CCAN in a variety of downstream analyses and applications, including removing cell cycle effects and batch effects of scRNA-seq datasets from different protocols, integrating paired and unpaired scRNA-seq and scATAC-seq data, accurately transferring cell type labels from scRNA-seq to scATAC-seq data, and characterizing the differentiation process from hematopoietic stem cells to different lineages in the integration of differentiation data.

scMoresDB: A comprehensive database of single-cell multi-omics data for human respiratory system

The human respiratory system is a complex and important system that can suffer a variety of diseases. Single-cell sequencing technologies, applied in many respiratory disease studies, have enhanced our ability in characterizing molecular and phenotypic features at a single-cell resolution. The exponentially increasing data from these studies have consequently led to difficulties in data sharing and analysis. Here, we present scMoresDB, a single-cell multi-omics database platform with extensive omics types tailored for human respiratory diseases. scMoresDB re-analyzes single-cell multi-omics datasets, providing a user-friendly interface with cross-omics search capabilities, interactive visualizations, and analytical tools for comprehensive data sharing and integrative analysis. Our example applications highlight the potential significance of BSG receptor in SARS-CoV-2 infection as well as the involvement of HHIP and TGFB2 in the development and progression of chronic obstructive pulmonary disease. scMoresDB significantly increases accessibility and utility of single-cell data relevant to human respiratory system and associated diseases.

Joint inference of clonal structure using single-cell genome and transcriptome sequencing data

Latest advancements in the high-throughput single-cell genome (scDNA) and transcriptome (scRNA) sequencing technologies enabled cell-resolved investigation of tissue clones. However, it remains challenging to cluster and couple single cells for heterogeneous scRNA and scDNA data generated from the same specimen. In this study, we present a computational framework called CCNMF, which employs a novel Coupled-Clone Non-negative Matrix Factorization technique to jointly infer clonal structure for matched scDNA and scRNA data. CCNMF couples multi-omics single cells by linking copy number and gene expression profiles through their general concordance. It successfully resolved the underlying coexisting clones with high correlations between the clonal genome and transcriptome from the same specimen. We validated that CCNMF can achieve high accuracy and robustness using both simulated benchmarks and real-world applications, including an ovarian cancer cell lines mixture, a gastric cancer cell line, and a primary gastric cancer. In summary, CCNMF provides a powerful tool for integrating multi-omics single-cell data, enabling simultaneous resolution of genomic and transcriptomic clonal architecture. This computational framework facilitates the understanding of how cellular gene expression changes in conjunction with clonal genome alternations, shedding light on the cellular genomic difference of subclones that contributes to tumor evolution.

Recent advances in plant translational genomics for crop improvement

The growing population, climate change, and limited agricultural resources put enormous pressure on agricultural systems. A plateau in crop yields is occurring and extreme weather events and urbanization threaten the livelihood of farmers. It is imperative that immediate attention is paid to addressing the increasing food demand, ensuring resilience against emerging threats, and meeting the demand for more nutritious, safer food. Under uncertain conditions, it is essential to expand genetic diversity and discover novel crop varieties or variations to develop higher and more stable yields. Genomics plays a significant role in developing abundant and nutrient-dense food crops. An alternative to traditional breeding approach, translational genomics is able to improve breeding programs in a more efficient and precise manner by translating genomic concepts into practical tools. Crop breeding based on genomics offers potential solutions to overcome the limitations of conventional breeding methods, including improved crop varieties that provide more nutritional value and are protected from biotic and abiotic stresses. Genetic markers, such as SNPs and ESTs, contribute to the discovery of QTLs controlling agronomic traits and stress tolerance. In order to meet the growing demand for food, there is a need to incorporate QTLs into breeding programs using marker-assisted selection/breeding and transgenic technologies. This chapter primarily focuses on the recent advances that are made in translational genomics for crop improvement and various omics techniques including transcriptomics, metagenomics, pangenomics, single cell omics etc. Numerous genome editing techniques including CRISPR Cas technology and their applications in crop improvement had been discussed.

Plant biotechnology research with single-cell transcriptome: recent advancements and prospects

KEY MESSAGE

Single-cell transcriptomic techniques have emerged as powerful tools in plant biology, offering high-resolution insights into gene expression at the individual cell level. This review highlights the rapid expansion of single-cell technologies in plants, their potential in understanding plant development, and their role in advancing plant biotechnology research. Single-cell techniques have emerged as powerful tools to enhance our understanding of biological systems, providing high-resolution transcriptomic analysis at the single-cell level. In plant biology, the adoption of single-cell transcriptomics has seen rapid expansion of available technologies and applications. This review article focuses on the latest advancements in the field of single-cell transcriptomic in plants and discusses the potential role of these approaches in plant development and expediting plant biotechnology research in the near future. Furthermore, inherent challenges and limitations of single-cell technology are critically examined to overcome them and enhance our knowledge and understanding.

A Cell Cycle-aware Network for Data Integration and Label Transferring of Single-cell RNA-seq and ATAC-seq

In recent years, the integration of single-cell multi-omics data has provided a more comprehensive understanding of cell functions and internal regulatory mechanisms from a non-single omics perspective, but it still suffers many challenges, such as omics-variance, sparsity, cell heterogeneity and confounding factors. As we know, cell cycle is regarded as a confounder when analyzing other factors in single-cell RNA-seq data, but it's not clear how it will work on the integrated single-cell multi-omics data. Here, we developed a Cell Cycle-Aware Network (CCAN) to remove cell cycle effects from the integrated single-cell multi-omics data while keeping the cell type-specific variations. This is the first computational model to study the cell-cycle effects in the integration of single-cell multi-omics data. Validations on several benchmark datasets show the out-standing performance of CCAN in a variety of downstream analyses and applications, including removing cell cycle effects and batch effects of scRNA-seq datasets from different protocols, integrating paired and unpaired scRNA-seq and scATAC-seq data, accurately transferring cell type labels from scRNA-seq to scATAC-seq data, and characterizing the differentiation process from hematopoietic stem cells to different lineages in the integration of differentiation data.

scDesign3 generates realistic in silico data for multimodal single-cell and spatial omics

We present a statistical simulator, scDesign3, to generate realistic single-cell and spatial omics data, including various cell states, experimental designs and feature modalities, by learning interpretable parameters from real data. Using a unified probabilistic model for single-cell and spatial omics data, scDesign3 infers biologically meaningful parameters; assesses the goodness-of-fit of inferred cell clusters, trajectories and spatial locations; and generates in silico negative and positive controls for benchmarking computational tools.

Introducing single cell stereo-sequencing technology to transform the plant transcriptome landscape

Single cell RNA-sequencing (scRNA-seq) advancements have helped detect transcriptional heterogeneities in biological samples. However, scRNA-seq cannot currently provide high-resolution spatial transcriptome information or identify subcellular organs in biological samples. These limitations have led to the development of spatially enhanced-resolution omics-sequencing (Stereo-seq), which combines spatial information with single cell transcriptomics to address the challenges of scRNA-seq alone. In this review, we discuss the advantages of Stereo-seq technology. We anticipate that the application of such an integrated approach in plant research will advance our understanding of biological process in the plant transcriptomics era. We conclude with an outlook of how such integration will enhance crop improvement.

InClust+: the deep generative framework with mask modules for multimodal data integration, imputation, and cross-modal generation

BACKGROUND

With the development of single-cell technology, many cell traits can be measured. Furthermore, the multi-omics profiling technology could jointly measure two or more traits in a single cell simultaneously. In order to process the various data accumulated rapidly, computational methods for multimodal data integration are needed.

RESULTS

Here, we present inClust+, a deep generative framework for the multi-omics. It's built on previous inClust that is specific for transcriptome data, and augmented with two mask modules designed for multimodal data processing: an input-mask module in front of the encoder and an output-mask module behind the decoder. InClust+ was first used to integrate scRNA-seq and MERFISH data from similar cell populations, and to impute MERFISH data based on scRNA-seq data. Then, inClust+ was shown to have the capability to integrate the multimodal data (e.g. tri-modal data with gene expression, chromatin accessibility and protein abundance) with batch effect. Finally, inClust+ was used to integrate an unlabeled monomodal scRNA-seq dataset and two labeled multimodal CITE-seq datasets, transfer labels from CITE-seq datasets to scRNA-seq dataset, and generate the missing modality of protein abundance in monomodal scRNA-seq data. In the above examples, the performance of inClust+ is better than or comparable to the most recent tools in the corresponding task.

CONCLUSIONS

The inClust+ is a suitable framework for handling multimodal data. Meanwhile, the successful implementation of mask in inClust+ means that it can be applied to other deep learning methods with similar encoder-decoder architecture to broaden the application scope of these models.

scMMT: a multi-use deep learning approach for cell annotation, protein prediction and embedding in single-cell RNA-seq data

Accurate cell type annotation in single-cell RNA-sequencing data is essential for advancing biological and medical research, particularly in understanding disease progression and tumor microenvironments. However, existing methods are constrained by single feature extraction approaches, lack of adaptability to immune cell types with similar molecular profiles but distinct functions and a failure to account for the impact of cell label noise on model accuracy, all of which compromise the precision of annotation. To address these challenges, we developed a supervised approach called scMMT. We proposed a novel feature extraction technique to uncover more valuable information. Additionally, we constructed a multi-task learning framework based on the GradNorm method to enhance the recognition of challenging immune cells and reduce the impact of label noise by facilitating mutual reinforcement between cell type annotation and protein prediction tasks. Furthermore, we introduced logarithmic weighting and label smoothing mechanisms to enhance the recognition ability of rare cell types and prevent model overconfidence. Through comprehensive evaluations on multiple public datasets, scMMT has demonstrated state-of-the-art performance in various aspects including cell type annotation, rare cell identification, dropout and label noise resistance, protein expression prediction and low-dimensional embedding representation.

PancrESS - a meta-analysis resource for understanding cell-type specific expression in the human pancreas

BACKGROUND

The human pancreas is composed of specialized cell types producing hormones and enzymes critical to human health. These specialized functions are the result of cell type-specific transcriptional programs which manifest in cell-specific gene expression. Understanding these programs is essential to developing therapies for pancreatic disorders. Transcription in the human pancreas has been widely studied by single-cell RNA technologies, however the diversity of protocols and analysis methods hinders their interpretability in the aggregate.

RESULTS

In this work, we perform a meta-analysis of pancreatic single-cell RNA sequencing data. We present a database for reference transcriptome abundances and cell-type specificity metrics. This database facilitates the identification and definition of marker genes within the pancreas. Additionally, we introduce a versatile tool which is freely available as an R package, and should permit integration into existing workflows. Our tool accepts count data files generated by widely-used single-cell gene expression platforms in their original format, eliminating an additional pre-formatting step. Although we designed it to calculate expression specificity of pancreas cell types, our tool is agnostic to the biological source of count data, extending its applicability to other biological systems.

CONCLUSIONS

Our findings enhance the current understanding of expression specificity within the pancreas, surpassing previous work in terms of scope and detail. Furthermore, our database and tool enable researchers to perform similar calculations in diverse biological systems, expanding the applicability of marker gene identification and facilitating comparative analyses.

Single-cell joint profiling of multiple epigenetic proteins and gene transcription

Sculpting the epigenome with a combination of histone modifications and transcription factor occupancy determines gene transcription and cell fate specification. Here, we first develop uCoTarget, utilizing a split-pool barcoding strategy for realizing ultrahigh-throughput single-cell joint profiling of multiple epigenetic proteins. Through extensive optimization for sensitivity and multimodality resolution, we demonstrate that uCoTarget enables simultaneous detection of five histone modifications (H3K27ac, H3K4me3, H3K4me1, H3K36me3, and H3K27me3) in 19,860 single cells. We applied uCoTarget to the in vitro generation of hematopoietic stem/progenitor cells (HSPCs) from human embryonic stem cells, presenting multimodal epigenomic profiles in 26,418 single cells. uCoTarget reveals establishment of pairing of HSPC enhancers (H3K27ac) and promoters (H3K4me3) and RUNX1 engagement priming for H3K27ac activation along the HSPC path. We then develop uCoTargetX, an expansion of uCoTarget to simultaneously measure transcriptome and multiple epigenome targets. Together, our methods enable generalizable, versatile multimodal profiles for reconstructing comprehensive epigenome and transcriptome landscapes and analyzing the regulatory interplay at single-cell level.

Improving Enhancer Identification with a Multi-Classifier Stacked Ensemble Model

Enhancers are DNA regions that are responsible for controlling the expression of genes. Enhancers are usually found upstream or downstream of a gene, or even inside a gene's intron region, but are normally located at a distant location from the genes they control. By integrating experimental and computational approaches, it is possible to uncover enhancers within DNA sequences, which possess regulatory properties. Experimental techniques such as ChIP-seq and ATAC-seq can identify genomic regions that are associated with transcription factors or accessible to regulatory proteins. On the other hand, computational techniques can predict enhancers based on sequence features and epigenetic modifications. In our study, we have developed a multi-classifier stacked ensemble (MCSE-enhancer) model that can accurately identify enhancers. We utilized feature descriptors from various physiochemical properties as input for our six baseline classifiers and built a stacked classifier, which outperformed previous enhancer classification techniques in terms of accuracy, specificity, sensitivity, and Mathew's correlation coefficient. Our model achieved an accuracy of 81.5%, representing a 2-3% improvement over existing models.

MANGEM: A web app for multimodal analysis of neuronal gene expression, electrophysiology, and morphology

Single-cell techniques like Patch-seq have enabled the acquisition of multimodal data from individual neuronal cells, offering systematic insights into neuronal functions. However, these data can be heterogeneous and noisy. To address this, machine learning methods have been used to align cells from different modalities onto a low-dimensional latent space, revealing multimodal cell clusters. The use of those methods can be challenging without computational expertise or suitable computing infrastructure for computationally expensive methods. To address this, we developed a cloud-based web application, MANGEM (multimodal analysis of neuronal gene expression, electrophysiology, and morphology). MANGEM provides a step-by-step accessible and user-friendly interface to machine learning alignment methods of neuronal multimodal data. It can run asynchronously for large-scale data alignment, provide users with various downstream analyses of aligned cells, and visualize the analytic results. We demonstrated the usage of MANGEM by aligning multimodal data of neuronal cells in the mouse visual cortex.

MAGNETO: Cell type marker panel generator from single-cell transcriptomic data

Single-cell RNA sequencing experiments produce data useful to identify different cell types, including uncharacterized and rare ones. This enables us to study the specific functional roles of these cells in different microenvironments and contexts. After identifying a (novel) cell type of interest, it is essential to build succinct marker panels, composed of a few genes referring to cell surface proteins and clusters of differentiation molecules, able to discriminate the desired cells from the other cell populations. In this work, we propose a fully-automatic framework called MAGNETO, which can help construct optimal marker panels starting from a single-cell gene expression matrix and a cell type identity for each cell. MAGNETO builds effective marker panels solving a tailored bi-objective optimization problem, where the first objective regards the identification of the genes able to isolate a specific cell type, while the second conflicting objective concerns the minimization of the total number of genes included in the panel. Our results on three public datasets show that MAGNETO can identify marker panels that identify the cell populations of interest better than state-of-the-art approaches. Finally, by fine-tuning MAGNETO, our results demonstrate that it is possible to obtain marker panels with different specificity levels.

Advanced Omics Techniques for Understanding Cochlear Genome, Epigenome, and Transcriptome in Health and Disease

Advanced genomics, transcriptomics, and epigenomics techniques are providing unprecedented insights into the understanding of the molecular underpinnings of the central nervous system, including the neuro-sensory cochlea of the inner ear. Here, we report for the first time a comprehensive and updated overview of the most advanced omics techniques for the study of nucleic acids and their applications in cochlear research. We describe the available in vitro and in vivo models for hearing research and the principles of genomics, transcriptomics, and epigenomics, alongside their most advanced technologies (like single-cell omics and spatial omics), which allow for the investigation of the molecular events that occur at a single-cell resolution while retaining the spatial information.

A live-cell platform to isolate phenotypically defined subpopulations for spatial multi-omic profiling

Numerous techniques have been employed to deconstruct the heterogeneity observed in normal and diseased cellular populations, including single cell RNA sequencing, in situ hybridization, and flow cytometry. While these approaches have revolutionized our understanding of heterogeneity, in isolation they cannot correlate phenotypic information within a physiologically relevant live-cell state with molecular profiles. This inability to integrate a live-cell phenotype-such as invasiveness, cell:cell interactions, and changes in spatial positioning-with multi-omic data creates a gap in understanding cellular heterogeneity. We sought to address this gap by employing lab technologies to design a detailed protocol, termed Spatiotemporal Genomic and Cellular Analysis (SaGA), for the precise imaging-based selection, isolation, and expansion of phenotypically distinct live cells. This protocol requires cells expressing a photoconvertible fluorescent protein and employs live cell confocal microscopy to photoconvert a user-defined single cell or set of cells displaying a phenotype of interest. The total population is then extracted from its microenvironment, and the optically highlighted cells are isolated using fluorescence activated cell sorting. SaGA-isolated cells can then be subjected to multi-omics analysis or cellular propagation for in vitro or in vivo studies. This protocol can be applied to a variety of conditions, creating protocol flexibility for user-specific research interests. The SaGA technique can be accomplished in one workday by non-specialists and results in a phenotypically defined cellular subpopulations for integration with multi-omics techniques. We envision this approach providing multi-dimensional datasets exploring the relationship between live cell phenotypes and multi-omic heterogeneity within normal and diseased cellular populations.

Application of single-cell multi-omics approaches in horticulture research

Cell heterogeneity shapes the morphology and function of various tissues and organs in multicellular organisms. Elucidation of the differences among cells and the mechanism of intercellular regulation is essential for an in-depth understanding of the developmental process. In recent years, the rapid development of high-throughput single-cell transcriptome sequencing technologies has influenced the study of plant developmental biology. Additionally, the accuracy and sensitivity of tools used to study the epigenome and metabolome have significantly increased, thus enabling multi-omics analysis at single-cell resolution. Here, we summarize the currently available single-cell multi-omics approaches and their recent applications in plant research, review the single-cell based studies in fruit, vegetable, and ornamental crops, and discuss the potential of such approaches in future horticulture research.

Integrated Multimodal Omics and Dietary Approaches for the Management of Neurodegeneration

Neurodegenerative diseases, such as Alzheimer's disease and Parkinson's disease, are caused by a combination of multiple events that damage neuronal function. A well-characterized biomarker of neurodegeneration is the accumulation of proteinaceous aggregates in the brain. However, the gradually worsening symptoms of neurodegenerative diseases are unlikely to be solely due to the result of a mutation in a single gene, but rather a multi-step process involving epigenetic changes. Recently, it has been suggested that a fraction of epigenetic alternations may be correlated to neurodegeneration in the brain. Unlike DNA mutations, epigenetic alterations are reversible, and therefore raise the possibilities for therapeutic intervention, including dietary modifications. Additionally, reactive oxygen species may contribute to the pathogenesis of Alzheimer's disease and Parkinson's disease through epigenetic alternation. Given that the antioxidant properties of plant-derived phytochemicals are likely to exhibit pleiotropic effects against ROS-mediated epigenetic alternation, dietary intervention may be promising for the management of neurodegeneration in these diseases. In this review, the state-of-the-art applications using single-cell multimodal omics approaches, including epigenetics, and dietary approaches for the identification of novel biomarkers and therapeutic approaches for the treatment of neurodegenerative diseases are discussed.

Dissecting and improving gene regulatory network inference using single-cell transcriptome data

Single-cell transcriptome data has been widely used to reconstruct gene regulatory networks (GRNs) controlling critical biological processes such as development and differentiation. Although a growing list of algorithms has been developed to infer GRNs using such data, achieving an inference accuracy consistently higher than random guessing has remained challenging. To address this, it is essential to delineate how the accuracy of regulatory inference is limited. Here, we systematically characterized factors limiting the accuracy of inferred GRNs and demonstrated that using pre-mRNA information can help improve regulatory inference compared to the typically used information (i.e., mature mRNA). Using kinetic modeling and simulated single-cell data sets, we showed that target genes' mature mRNA levels often fail to accurately report upstream regulatory activities because of gene-level and network-level factors, which can be improved by using pre-mRNA levels. We tested this finding on public single-cell RNA-seq data sets using intronic reads as proxies of pre-mRNA levels and can indeed achieve a higher inference accuracy compared to using exonic reads (corresponding to mature mRNAs). Using experimental data sets, we further validated findings from the simulated data sets and identified factors such as transcription factor activity dynamics influencing the accuracy of pre-mRNA-based inference. This work delineates the fundamental limitations of gene regulatory inference and helps improve GRN inference using single-cell RNA-seq data.

Prediction of distinct populations of innate lymphoid cells by transcriptional profiles

Innate lymphoid cells (ILCs) are a unique type of lymphocyte that differ from adaptive lymphocytes in that they lack antigen receptors, which primarily reside in tissues and are closely associated with fibers. Despite their plasticity and heterogeneity, identifying ILCs in peripheral blood can be difficult due to their small numbers. Accurately and rapidly identifying ILCs is critical for studying homeostasis and inflammation. To address this challenge, we collect single-cell RNA-seq data from 647 patients, including 26,087 transcripts. Background screening, Lasso analysis, and principal component analysis (PCA) are used to select features. Finally, we employ a deep neural network to classify lymphocytes. Our method achieved the highest accuracy compared to other approaches. Furthermore, we identified four genes that play a vital role in lymphocyte development. Adding these gene transcripts into model, we were able to increase the model's AUC. In summary, our study demonstrates the effectiveness of using single-cell transcriptomic analysis combined with machine learning techniques to accurately identify congenital lymphoid cells and advance our understanding of their development and function in the body.

Epstein-Barr virus evades restrictive host chromatin closure by subverting B cell activation and germinal center regulatory loci

Chromatin accessibility fundamentally governs gene expression and biological response programs that can be manipulated by pathogens. Here we capture dynamic chromatin landscapes of individual B cells during Epstein-Barr virus (EBV) infection. EBV+ cells that exhibit arrest via antiviral sensing and proliferation-linked DNA damage experience global accessibility reduction. Proliferative EBV+ cells develop expression-linked architectures and motif accessibility profiles resembling in vivo germinal center (GC) phenotypes. Remarkably, EBV elicits dark zone (DZ), light zone (LZ), and post-GC B cell chromatin features despite BCL6 downregulation. Integration of single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq), single-cell RNA sequencing (scRNA-seq), and chromatin immunoprecipitation sequencing (ChIP-seq) data enables genome-wide cis-regulatory predictions implicating EBV nuclear antigens (EBNAs) in phenotype-specific control of GC B cell activation, survival, and immune evasion. Knockouts validate bioinformatically identified regulators (MEF2C and NFE2L2) of EBV-induced GC phenotypes and EBNA-associated loci that regulate gene expression (CD274/PD-L1). These data and methods can inform high-resolution investigations of EBV-host interactions, B cell fates, and virus-mediated lymphomagenesis.

Integration of gene expression and DNA methylation data across different experiments

Integrative analysis of multi-omic datasets has proven to be extremely valuable in cancer research and precision medicine. However, obtaining multimodal data from the same samples is often difficult. Integrating multiple datasets of different omics remains a challenge, with only a few available algorithms developed to solve it. Here, we present INTEND (IntegratioN of Transcriptomic and EpigeNomic Data), a novel algorithm for integrating gene expression and DNA methylation datasets covering disjoint sets of samples. To enable integration, INTEND learns a predictive model between the two omics by training on multi-omic data measured on the same set of samples. In comprehensive testing on 11 TCGA (The Cancer Genome Atlas) cancer datasets spanning 4329 patients, INTEND achieves significantly superior results compared with four state-of-the-art integration algorithms. We also demonstrate INTEND's ability to uncover connections between DNA methylation and the regulation of gene expression in the joint analysis of two lung adenocarcinoma single-omic datasets from different sources. INTEND's data-driven approach makes it a valuable multi-omic data integration tool. The code for INTEND is available at https://github.com/Shamir-Lab/INTEND.

Homogeneous Space Construction and Projection for Single-Cell Expression Prediction Based on Deep Learning

Predicting cellular responses to perturbations is an unsolved problem in biology. Traditional approaches assume that different cell types respond similarly to perturbations. However, this assumption does not take into account the context of genome interactions in different cell types, which leads to compromised prediction quality. More recently, deep learning models used to discover gene-gene relationships can yield more accurate predictions of cellular responses. The huge difference in biological information between different cell types makes it difficult for deep learning models to encode data into a continuous low-dimensional feature space, which means that the features captured by the latent space may not be continuous. Therefore, the mapping relationship between the two conditional spaces learned by the model can only be applied where the real reference data resides, leading to the wrong mapping of the predicted target cells because they are not in the same domain as the reference data. In this paper, we propose an information-navigated variational autoencoder (INVAE), a deep neural network for cell perturbation response prediction. INVAE filters out information that is not conducive to predictive performance. For the remaining information, INVAE constructs a homogeneous space of control conditions, and finds the mapping relationship between the control condition space and the perturbation condition space. By embedding the target unit into the control space and then mapping it to the perturbation space, we can predict the perturbed state of the target unit. Comparing our proposed method with other three state-of-the-art methods on three real datasets, experimental results show that INVAE outperforms existing methods in cell state prediction after perturbation. Furthermore, we demonstrate that filtering out useless information not only improves prediction accuracy but also reveals similarities in how genes in different cell types are regulated following perturbation.

GTM-decon: guided-topic modeling of single-cell transcriptomes enables sub-cell-type and disease-subtype deconvolution of bulk transcriptomes

Cell-type composition is an important indicator of health. We present Guided Topic Model for deconvolution (GTM-decon) to automatically infer cell-type-specific gene topic distributions from single-cell RNA-seq data for deconvolving bulk transcriptomes. GTM-decon performs competitively on deconvolving simulated and real bulk data compared with the state-of-the-art methods. Moreover, as demonstrated in deconvolving disease transcriptomes, GTM-decon can infer multiple cell-type-specific gene topic distributions per cell type, which captures sub-cell-type variations. GTM-decon can also use phenotype labels from single-cell or bulk data to infer phenotype-specific gene distributions. In a nested-guided design, GTM-decon identified cell-type-specific differentially expressed genes from bulk breast cancer transcriptomes.

Single-Cell Multiomics

Single-cell RNA sequencing methods have led to improved understanding of the heterogeneity and transcriptomic states present in complex biological systems. Recently, the development of novel single-cell technologies for assaying additional modalities, specifically genomic, epigenomic, proteomic, and spatial data, allows for unprecedented insight into cellular biology. While certain technologies collect multiple measurements from the same cells simultaneously, even when modalities are separately assayed in different cells, we can apply novel computational methods to integrate these data. The application of computational integration methods to multimodal paired and unpaired data results in rich information about the identities of the cells present and the interactions between different levels of biology, such as between genetic variation and transcription. In this review, we both discuss the single-cell technologies for measuring these modalities and describe and characterize a variety of computational integration methods for combining the resulting data to leverage multimodal information toward greater biological insight.

Methods and applications for single-cell and spatial multi-omics

The joint analysis of the genome, epigenome, transcriptome, proteome and/or metabolome from single cells is transforming our understanding of cell biology in health and disease. In less than a decade, the field has seen tremendous technological revolutions that enable crucial new insights into the interplay between intracellular and intercellular molecular mechanisms that govern development, physiology and pathogenesis. In this Review, we highlight advances in the fast-developing field of single-cell and spatial multi-omics technologies (also known as multimodal omics approaches), and the computational strategies needed to integrate information across these molecular layers. We demonstrate their impact on fundamental cell biology and translational research, discuss current challenges and provide an outlook to the future.

Fatty Acid Metabolites and the Tumor Microenvironment as Potent Regulators of Cancer Stem Cell Signaling

Individual cancer cells are not equal but are organized into a cellular hierarchy in which only a rare few leukemia cells can self-renew in a manner reminiscent of the characteristic stem cell properties. The PI3K/AKT pathway functions in a variety of cancers and plays a critical role in the survival and proliferation of healthy cells under physiologic conditions. In addition, cancer stem cells might exhibit a variety of metabolic reprogramming phenotypes that cannot be completely attributed to the intrinsic heterogeneity of cancer. Given the heterogeneity of cancer stem cells, new strategies with single-cell resolution will become a powerful tool to eradicate the aggressive cell population harboring cancer stem cell phenotypes. Here, this article will provide an overview of the most important signaling pathways of cancer stem cells regarding their relevance to the tumor microenvironment and fatty acid metabolism, suggesting valuable strategies among cancer immunotherapies to inhibit the recurrence of tumors.

Multi-task learning from multimodal single-cell omics with Matilda

Multimodal single-cell omics technologies enable multiple molecular programs to be simultaneously profiled at a global scale in individual cells, creating opportunities to study biological systems at a resolution that was previously inaccessible. However, the analysis of multimodal single-cell omics data is challenging due to the lack of methods that can integrate across multiple data modalities generated from such technologies. Here, we present Matilda, a multi-task learning method for integrative analysis of multimodal single-cell omics data. By leveraging the interrelationship among tasks, Matilda learns to perform data simulation, dimension reduction, cell type classification, and feature selection in a single unified framework. We compare Matilda with other state-of-the-art methods on datasets generated from some of the most popular multimodal single-cell omics technologies. Our results demonstrate the utility of Matilda for addressing multiple key tasks on integrative multimodal single-cell omics data analysis. Matilda is implemented in Pytorch and is freely available from https://github.com/PYangLab/Matilda.

scJVAE: A novel method for integrative analysis of multimodal single-cell data

The study of cellular decision-making can be approached comprehensively using multimodal single-cell omics technology. Recent advances in multimodal single-cell technology have enabled simultaneous profiling of more than one modality from the same cell, providing more significant insights into cell characteristics. However, learning the joint representation of multimodal single-cell data is challenging due to batch effects. Here we present a novel method, scJVAE (single-cell Joint Variational AutoEncoder), for batch effect removal and joint representation of multimodal single-cell data. The scJVAE integrates and learns joint embedding of paired scRNA-seq and scATAC-seq data modalities. We evaluate and demonstrate the ability of scJVAE to remove batch effects using various datasets with paired gene expression and open chromatin. We also consider scJVAE for downstream analysis, such as lower dimensional representation, cell-type clustering, and time and memory requirement. We find scJVAE a robust and scalable method outperforming existing state-of-the-art batch effect removal and integration methods.

Single-cell transcriptomics for the assessment of cardiac disease

Cardiovascular disease is the leading cause of death globally. An advanced understanding of cardiovascular disease mechanisms is required to improve therapeutic strategies and patient risk stratification. State-of-the-art, large-scale, single-cell and single-nucleus transcriptomics facilitate the exploration of the cardiac cellular landscape at an unprecedented level, beyond its descriptive features, and can further our understanding of the mechanisms of disease and guide functional studies. In this Review, we provide an overview of the technical challenges in the experimental design of single-cell and single-nucleus transcriptomics studies, as well as a discussion of the type of inferences that can be made from the data derived from these studies. Furthermore, we describe novel findings derived from transcriptomics studies for each major cardiac cell type in both health and disease, and from development to adulthood. This Review also provides a guide to interpreting the exhaustive list of newly identified cardiac cell types and states, and highlights the consensus and discordances in annotation, indicating an urgent need for standardization. We describe advanced applications such as integration of single-cell data with spatial transcriptomics to map genes and cells on tissue and define cellular microenvironments that regulate homeostasis and disease progression. Finally, we discuss current and future translational and clinical implications of novel transcriptomics approaches, and provide an outlook of how these technologies will change the way we diagnose and treat heart disease.

Matching single cells across modalities with contrastive learning and optimal transport

Understanding the interactions between the biomolecules that govern cellular behaviors remains an emergent question in biology. Recent advances in single-cell technologies have enabled the simultaneous quantification of multiple biomolecules in the same cell, opening new avenues for understanding cellular complexity and heterogeneity. Still, the resulting multimodal single-cell datasets present unique challenges arising from the high dimensionality and multiple sources of acquisition noise. Computational methods able to match cells across different modalities offer an appealing alternative towards this goal. In this work, we propose MatchCLOT, a novel method for modality matching inspired by recent promising developments in contrastive learning and optimal transport. MatchCLOT uses contrastive learning to learn a common representation between two modalities and applies entropic optimal transport as an approximate maximum weight bipartite matching algorithm. Our model obtains state-of-the-art performance on two curated benchmarking datasets and an independent test dataset, improving the top scoring method by 26.1% while preserving the underlying biological structure of the multimodal data. Importantly, MatchCLOT offers high gains in computational time and memory that, in contrast to existing methods, allows it to scale well with the number of cells. As single-cell datasets become increasingly large, MatchCLOT offers an accurate and efficient solution to the problem of modality matching.

MANGEM: a web app for Multimodal Analysis of Neuronal Gene expression, Electrophysiology and Morphology

Single-cell techniques have enabled the acquisition of multi-modal data, particularly for neurons, to characterize cellular functions. Patch-seq, for example, combines patch-clamp recording, cell imaging, and single-cell RNA-seq to obtain electrophysiology, morphology, and gene expression data from a single neuron. While these multi-modal data offer potential insights into neuronal functions, they can be heterogeneous and noisy. To address this, machine-learning methods have been used to align cells from different modalities onto a low-dimensional latent space, revealing multi-modal cell clusters. However, the use of those methods can be challenging for biologists and neuroscientists without computational expertise and also requires suitable computing infrastructure for computationally expensive methods. To address these issues, we developed a cloud-based web application, MANGEM (Multimodal Analysis of Neuronal Gene expression, Electrophysiology, and Morphology) at https://ctc.waisman.wisc.edu/mangem. MANGEM provides a step-by-step accessible and user-friendly interface to machine-learning alignment methods of neuronal multi-modal data while enabling real-time visualization of characteristics of raw and aligned cells. It can be run asynchronously for large-scale data alignment, provides users with various downstream analyses of aligned cells and visualizes the analytic results such as identifying multi-modal cell clusters of cells and detecting correlated genes with electrophysiological and morphological features. We demonstrated the usage of MANGEM by aligning Patch-seq multimodal data of neuronal cells in the mouse visual cortex.

Multi-view clustering by CPS-merge analysis with application to multimodal single-cell data

Multi-view data can be generated from diverse sources, by different technologies, and in multiple modalities. In various fields, integrating information from multi-view data has pushed the frontier of discovery. In this paper, we develop a new approach for multi-view clustering, which overcomes the limitations of existing methods such as the need of pooling data across views, restrictions on the clustering algorithms allowed within each view, and the disregard for complementary information between views. Our new method, called CPS-merge analysis, merges clusters formed by the Cartesian product of single-view cluster labels, guided by the principle of maximizing clustering stability as evaluated by CPS analysis. In addition, we introduce measures to quantify the contribution of each view to the formation of any cluster. CPS-merge analysis can be easily incorporated into an existing clustering pipeline because it only requires single-view cluster labels instead of the original data. We can thus readily apply advanced single-view clustering algorithms. Importantly, our approach accounts for both consensus and complementary effects between different views, whereas existing ensemble methods focus on finding a consensus for multiple clustering results, implying that results from different views are variations of one clustering structure. Through experiments on single-cell datasets, we demonstrate that our approach frequently outperforms other state-of-the-art methods.

A live-cell platform to isolate phenotypically defined subpopulations for spatial multi-omic profiling

Numerous techniques have been employed to deconstruct the heterogeneity observed in normal and diseased cellular populations, including single cell RNA sequencing, in situ hybridization, and flow cytometry. While these approaches have revolutionized our understanding of heterogeneity, in isolation they cannot correlate phenotypic information within a physiologically relevant live-cell state, with molecular profiles. This inability to integrate a historical live-cell phenotype, such as invasiveness, cell:cell interactions, and changes in spatial positioning, with multi-omic data, creates a gap in understanding cellular heterogeneity. We sought to address this gap by employing lab technologies to design a detailed protocol, termed Spatiotemporal Genomics and Cellular Analysis (SaGA), for the precise imaging-based selection, isolation, and expansion of phenotypically distinct live-cells. We begin with cells stably expressing a photoconvertible fluorescent protein and employ live cell confocal microscopy to photoconvert a user-defined single cell or set of cells displaying a phenotype of interest. The total population is then extracted from its microenvironment, and the optically highlighted cells are isolated using fluorescence activated cell sorting. SaGA-isolated cells can then be subjected to multi-omics analysis or cellular propagation for in vitro or in vivo studies. This protocol can be applied to a variety of conditions, creating protocol flexibility for user-specific research interests. The SaGA technique can be accomplished in one workday by non-specialists and results in a phenotypically defined cellular subpopulation for integration with multi-omics techniques. We envision this approach providing multi-dimensional datasets exploring the relationship between live-cell phenotype and multi-omic heterogeneity within normal and diseased cellular populations.

BRAQUE: Bayesian Reduction for Amplified Quantization in UMAP Embedding

Single-cell biology has revolutionized the way we understand biological processes. In this paper, we provide a more tailored approach to clustering and analyzing spatial single-cell data coming from immunofluorescence imaging techniques. We propose Bayesian Reduction for Amplified Quantization in UMAP Embedding (BRAQUE) as an integrative novel approach, from data preprocessing to phenotype classification. BRAQUE starts with an innovative preprocessing, named Lognormal Shrinkage, which is able to enhance input fragmentation by fitting a lognormal mixture model and shrink each component towards its median, in order to help further the clustering step in finding more separated and clear clusters. Then, BRAQUE's pipeline consists of a dimensionality reduction step performed using UMAP, and a clustering performed using HDBSCAN on UMAP embedding. In the end, clusters are assigned to a cell type by experts, using effects size measures to rank markers and identify characterizing markers (Tier 1), and possibly characterize markers (Tier 2). The number of total cell types in one lymph node detectable with these technologies is unknown and difficult to predict or estimate. Therefore, with BRAQUE, we achieved a higher granularity than other similar algorithms such as PhenoGraph, following the idea that merging similar clusters is easier than splitting unclear ones into clear subclusters.

Single-cell omics: A new direction for functional genetic research in human diseases and animal models

Over the past decade, with the development of high-throughput single-cell sequencing technology, single-cell omics has been emerged as a powerful tool to understand the molecular basis of cellular mechanisms and refine our knowledge of diverse cell states. They can reveal the heterogeneity at different genetic layers and elucidate their associations by multiple omics analysis, providing a more comprehensive genetic map of biological regulatory networks. In the post-GWAS era, the molecular biological mechanisms influencing human diseases will be further elucidated by single-cell omics. This review mainly summarizes the development and trend of single-cell omics. This involves single-cell omics technologies, single-cell multi-omics technologies, multiple omics data integration methods, applications in various human organs and diseases, classic laboratory cell lines, and animal disease models. The review will reveal some perspectives for elucidating human diseases and constructing animal models.

Multipotent Embryonic Lung Progenitors: Foundational Units of In Vitro and In Vivo Lung Organogenesis

Transient, tissue-specific, embryonic progenitors are important cell populations in vertebrate development. In the course of respiratory system development, multipotent mesenchymal and epithelial progenitors drive the diversification of fates that results to the plethora of cell types that compose the airways and alveolar space of the adult lungs. Use of mouse genetic models, including lineage tracing and loss-of-function studies, has elucidated signaling pathways that guide proliferation and differentiation of embryonic lung progenitors as well as transcription factors that underlie lung progenitor identity. Furthermore, pluripotent stem cell-derived and ex vivo expanded respiratory progenitors offer novel, tractable, high-fidelity systems that allow for mechanistic studies of cell fate decisions and developmental processes. As our understanding of embryonic progenitor biology deepens, we move closer to the goal of in vitro lung organogenesis and resulting applications in developmental biology and medicine.

A unified computational framework for single-cell data integration with optimal transport

Single-cell data integration can provide a comprehensive molecular view of cells. However, how to integrate heterogeneous single-cell multi-omics as well as spatially resolved transcriptomic data remains a major challenge. Here we introduce uniPort, a unified single-cell data integration framework that combines a coupled variational autoencoder (coupled-VAE) and minibatch unbalanced optimal transport (Minibatch-UOT). It leverages both highly variable common and dataset-specific genes for integration to handle the heterogeneity across datasets, and it is scalable to large-scale datasets. uniPort jointly embeds heterogeneous single-cell multi-omics datasets into a shared latent space. It can further construct a reference atlas for gene imputation across datasets. Meanwhile, uniPort provides a flexible label transfer framework to deconvolute heterogeneous spatial transcriptomic data using an optimal transport plan, instead of embedding latent space. We demonstrate the capability of uniPort by applying it to integrate a variety of datasets, including single-cell transcriptomics, chromatin accessibility, and spatially resolved transcriptomic data.

AIscEA: unsupervised integration of single-cell gene expression and chromatin accessibility via their biological consistency

MOTIVATION

The integrative analysis of single-cell gene expression and chromatin accessibility measurements is essential for revealing gene regulation, but it is one of the key challenges in computational biology. Gene expression and chromatin accessibility are measurements from different modalities, and no common features can be directly used to guide integration. Current state-of-the-art methods lack practical solutions for finding heterogeneous clusters. However, previous methods might not generate reliable results when cluster heterogeneity exists. More importantly, current methods lack an effective way to select hyper-parameters under an unsupervised setting. Therefore, applying computational methods to integrate single-cell gene expression and chromatin accessibility measurements remains difficult.

RESULTS

We introduce AIscEA-Alignment-based Integration of single-cell gene Expression and chromatin Accessibility-a computational method that integrates single-cell gene expression and chromatin accessibility measurements using their biological consistency. AIscEA first defines a ranked similarity score to quantify the biological consistency between cell clusters across measurements. AIscEA then uses the ranked similarity score and a novel permutation test to identify cluster alignment across measurements. AIscEA further utilizes graph alignment for the aligned cell clusters to align the cells across measurements. We compared AIscEA with the competing methods on several benchmark datasets and demonstrated that AIscEA is highly robust to the choice of hyper-parameters and can better handle the cluster heterogeneity problem. Furthermore, AIscEA significantly outperforms the state-of-the-art methods when integrating real-world SNARE-seq and scMultiome-seq datasets in terms of integration accuracy.

AVAILABILITY AND IMPLEMENTATION

AIscEA is available at https://figshare.com/articles/software/AIscEA_zip/21291135 on FigShare as well as {https://github.com/elhaam/AIscEA} onGitHub.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.