Search

Recent advances in spatial transcriptomics have enabled measurements of gene expression at cell/spot resolution meanwhile retaining both the spatial information and the histopathological images of the tissues. Deciphering the spatial domains of spots in the tissues is a vital step for various downstream tasks in spatial transcriptomics analysis. Existing methods have been developed for this purpose by combining gene expression and histopathological images to conquer noises in gene expression. However, current methods only use the histopathological images to construct spot relations without updating in the training stage, or simply concatenate the information from gene expression and images into one feature vector. Here, we propose a novel method ConGI to accurately decipher spatial domains by integrating gene expression and histopathological images, where the gene expression is adapted to image information through contrastive learning. We introduce three contrastive loss functions within and between modalities to learn the common semantic representations across all modalities while avoiding their meaningless modality-private noise information. The learned representations are then used for deciphering spatial domains through a clustering method. By comprehensive tests on tumor and normal spatial transcriptomics datasets, ConGI was shown to outperform existing methods in terms of spatial domain identification. More importantly, the learned representations from our model have also been used efficiently for various downstream tasks, including trajectory inference, clustering, and visualization.

Clustering analysis is widely utilized in single-cell RNA-sequencing (scRNA-seq) data to discover cell heterogeneity and cell states. While many clustering methods have been developed for scRNA-seq analysis, most of these methods require to provide the number of clusters. However, it is not easy to know the exact number of cell types in advance, and experienced determination is not always reliable. Here, we have developed ADClust, an automatic deep embedding clustering method for scRNA-seq data, which can accurately cluster cells without requiring a predefined number of clusters. Specifically, ADClust first obtains low-dimensional representation through pre-trained autoencoder, and uses the representations to cluster cells into initial micro-clusters. The clusters are then compared in between by a statistical test, and similar micro-clusters are merged into larger clusters. According to the clustering, cell representations are updated so that each cell will be pulled toward centres of its assigned cluster and similar clusters, while cells are separated to keep distances between clusters. This is accomplished through jointly optimizing the carefully designed clustering and autoencoder loss functions. This merging process continues until convergence. ADClust was tested on eleven real scRNA-seq datasets, and shown to outperform existing methods in terms of both clustering performance and the accuracy on the number of the determined clusters. More importantly, our model provides high speed and scalability for large datasets.

Single-cell RNA sequencing (scRNA-seq) techniques provide high-resolution data on cellular heterogeneity in diverse tissues, and a critical step for the data analysis is cell type identification. Traditional methods usually cluster the cells and manually identify cell clusters through marker genes, which is time-consuming and subjective. With the launch of several large-scale single-cell projects, millions of sequenced cells have been annotated and it is promising to transfer labels from the annotated datasets to newly generated datasets. One powerful way for the transferring is to learn cell relations through the graph neural network (GNN), but traditional GNNs are difficult to process millions of cells due to the expensive costs of the message-passing procedure at each training epoch. Here, we have developed a robust and scalable GNN-based method for accurate single-cell classification (GraphCS), where the graph is constructed to connect similar cells within and between labelled and unlabeled scRNA-seq datasets for propagation of shared information. To overcome the slow information propagation of GNN at each training epoch, the diffused information is pre-calculated via the approximate Generalized PageRank algorithm, enabling sublinear complexity over cell numbers. Compared with existing methods, GraphCS demonstrates better performance on simulated, cross-platform, cross-species and cross-omics scRNA-seq datasets. More importantly, our model provides a high speed and scalability on large datasets, and can achieve superior performance for 1 million cells within 50 min.

Background The Cox proportional hazards model with neural networks is widely used to accurately predict survival outcome for choosing cancer treatment strategies. Although this method has shown outstanding performance in many tasks, it has encountered challenges when dealing with high-dimensional datasets. In this study, we point out that the Cox network has estimation bias in processing such datasets with a large number of censored samples. The estimation bias is composed of censored estimation bias and variance estimation bias, which limit the prediction performance of the model. In order to correct this bias, this paper proposes the Deep Bayesian Perturbation Cox Network (DBP), which introduces Bayesian prior knowledge about censored samples to optimize the training process of the neural network. Specifically, the model uses a sampling module called Bayesian Perturbation to approximate the prior knowledge, which can be used as a component for other Cox-based neural networks. Results The comparison between DBP and the previous model in different kinds of genomic datasets demonstrates that our model has made significant improvements over previous state-of-the-art methods. In addition, the simulation experiments are performed to illustrate how the DBP method addresses the bias caused by Cox Network. In the case study, based on the predicted risks in BRCA data from TCGA, we identify 400 differential expressed genes and 20 KEGG pathways that are associated with breast cancer prognosis, among which 65% of the top 20 genes have been proved by literature review. Conclusion Overall, these results demonstrate that our proposed method is advanced and robust in datasets with a large proportion of censored samples. Besides, it can guide to discover disease-related genes.