scholarly journals Exploring and analysing single cell multi-omics data with VDJView

2019 ◽  
Author(s):  
Jerome Samir ◽  
Simone Rizzetto ◽  
Money Gupta ◽  
Fabio Luciani

Abstract Background Single cell RNA sequencing provides unprecedented opportunity to simultaneously explore the transcriptomic and immune receptor diversity of T and B cells. However, there are limited tools available that simultaneously analyse large multi-omics datasets integrated with metadata such as patient and clinical information.Results We developed VDJView, which permits the simultaneous or independent analysis and visualisation of gene expression, immune receptors, and clinical metadata of both T and B cells. This tool is implemented as an easy-to-use R shiny web-application, which integrates numerous gene expression and TCR analysis tools, and accepts data from plate-based sorted or high-throughput single cell platforms. We utilised VDJView to analyse several 10X scRNA-seq datasets, including a recent dataset of 150,000 CD8+ T cells with available gene expression, TCR sequences, quantification of 15 surface proteins, and 44 antigen specificities (across viruses, cancer, and self-antigens). We performed quality control, filtering of tetramer non-specific cells, clustering, random sampling and hypothesis testing to discover antigen specific gene signatures which were associated with immune cell differentiation states and clonal expansion across the pathogen specific T cells. We also analysed 563 single cells (plate-based sorted) obtained from 11 subjects, revealing clonally expanded T and B cells across primary cancer tissues and metastatic lymph-node. These immune cells clustered with distinct gene signatures according to the breast cancer molecular subtype. VDJView has been tested in lab meetings and peer-to-peer discussions, showing effective data generation and discussion without the need to consult bioinformaticians.Conclusions VDJView enables researchers without profound bioinformatics skills to analyse immune scRNA-seq data, integrating and visualising this with clonality and metadata profiles, thus accelerating the process of hypothesis testing, data interpretation and discovery of cellular heterogeneity. VDJView is freely available at https://bitbucket.org/kirbyvisp/vdjview .

2019 ◽  
Author(s):  
Jerome Samir ◽  
Simone Rizzetto ◽  
Money Gupta ◽  
Fabio Luciani

Abstract Background Single cell RNA sequencing provides unprecedented opportunity to simultaneously explore the transcriptomic and immune receptor diversity of T and B cells. However, there are limited tools available that simultaneously analyse large multi-omics datasets integrated with metadata such as patient and clinical information.Results We developed VDJView, which permits the simultaneous or independent analysis and visualisation of gene expression, immune receptors, and clinical metadata of both T and B cells. This tool is implemented as an easy-to-use R shiny web-application, which integrates numerous gene expression and TCR analysis tools, and accepts data from plate-based sorted or high-throughput single cell platforms. We utilised VDJView to analyse several 10X scRNA-seq datasets, including a recent dataset of 150,000 CD8+ T cells with available gene expression, TCR sequences, quantification of 15 surface proteins, and 44 antigen specificities (across viruses, cancer, and self-antigens). We performed quality control, filtering of tetramer non-specific cells, clustering, random sampling and hypothesis testing to discover antigen specific gene signatures which were associated with immune cell differentiation states and clonal expansion across the pathogen specific T cells. We also analysed 563 single cells (plate-based sorted) obtained from 11 subjects, revealing clonally expanded T and B cells across primary cancer tissues and metastatic lymph-node. These immune cells clustered with distinct gene signatures according to the breast cancer molecular subtype. VDJView has been tested in lab meetings and peer-to-peer discussions, showing effective data generation and discussion without the need to consult bioinformaticians.Conclusions VDJView enables researchers without profound bioinformatics skills to analyse immune scRNA-seq data, integrating and visualising this with clonality and metadata profiles, thus accelerating the process of hypothesis testing, data interpretation and discovery of cellular heterogeneity. VDJView is freely available at https://bitbucket.org/kirbyvisp/vdjview .


2020 ◽  
Author(s):  
Jerome Samir ◽  
Simone Rizzetto ◽  
Money Gupta ◽  
Fabio Luciani

Abstract Background Single cell RNA sequencing provides unprecedented opportunity to simultaneously explore the transcriptomic and immune receptor diversity of T and B cells. However, there are limited tools available that simultaneously analyse large multi-omics datasets integrated with metadata such as patient and clinical information.Results We developed VDJView, which permits the simultaneous or independent analysis and visualisation of gene expression, immune receptors, and clinical metadata of both T and B cells. This tool is implemented as an easy-to-use R shiny web-application, which integrates numerous gene expression and TCR analysis tools, and accepts data from plate-based sorted or high-throughput single cell platforms. We utilised VDJView to analyse several 10X scRNA-seq datasets, including a recent dataset of 150,000 CD8+ T cells with available gene expression, TCR sequences, quantification of 15 surface proteins, and 44 antigen specificities (across viruses, cancer, and self-antigens). We performed quality control, filtering of tetramer non-specific cells, clustering, random sampling and hypothesis testing to discover antigen specific gene signatures which were associated with immune cell differentiation states and clonal expansion across the pathogen specific T cells. We also analysed 563 single cells (plate-based sorted) obtained from 11 subjects, revealing clonally expanded T and B cells across primary cancer tissues and metastatic lymph-node. These immune cells clustered with distinct gene signatures according to the breast cancer molecular subtype. VDJView has been tested in lab meetings and peer-to-peer discussions, showing effective data generation and discussion without the need to consult bioinformaticians.Conclusions VDJView enables researchers without profound bioinformatics skills to analyse immune scRNA-seq data, integrating and visualising this with clonality and metadata profiles, thus accelerating the process of hypothesis testing, data interpretation and discovery of cellular heterogeneity. VDJView is freely available at https://bitbucket.org/kirbyvisp/vdjview .


2019 ◽  
Author(s):  
Jerome Samir ◽  
Simone Rizzetto ◽  
Money Gupta ◽  
Fabio Luciani

AbstractBackgroundSingle cell RNA sequencing provides unprecedented opportunity to simultaneously explore the transcriptomic and immune receptor diversity of T and B cells. However, there are limited tools available that simultaneously analyse large multi-omics datasets integrated with metadata such as patient and clinical information.ResultsWe developed VDJView, which permits the simultaneous or independent analysis and visualisation of gene expression, immune receptors, and clinical metadata of both T and B cells. This tool is implemented as an easy-to-use R shiny web-application, which integrates numerous gene expression and TCR analysis tools, and accepts data from plate-based sorted or high-throughput single cell platforms. We utilised VDJView to analyse several 10X scRNA-seq datasets, including a recent dataset of 150,000 CD8+ T cells with available gene expression, TCR sequences, quantification of 15 surface proteins, and 44 antigen specificities (across viruses, cancer, and self-antigens). We performed quality control, filtering of tetramer non-specific cells, clustering, random sampling and hypothesis testing to discover antigen specific gene signatures which were associated with immune cell differentiation states and clonal expansion across the pathogen specific T cells. We also analysed 563 single cells (plate-based sorted) obtained from 11 subjects, revealing clonally expanded T and B cells across primary cancer tissues and metastatic lymph-node. These immune cells clustered with distinct gene signatures according to the breast cancer molecular subtype. VDJView has been tested in lab meetings and peer-to-peer discussions, showing effective data generation and discussion without the need to consult bioinformaticians.ConclusionsVDJView enables researchers without profound bioinformatics skills to analyse immune scRNA-seq data, integrating and visualising this with clonality and metadata profiles, thus accelerating the process of hypothesis testing, data interpretation and discovery of cellular heterogeneity. VDJView is freely available at https://bitbucket.org/kirbyvisp/vdjview.


2020 ◽  
Author(s):  
Jerome Samir ◽  
Simone Rizzetto ◽  
Money Gupta ◽  
Fabio Luciani

Abstract Background Single cell RNA sequencing provides unprecedented opportunity to simultaneously explore the transcriptomic and immune receptor diversity of T and B cells. However, there are limited tools available that simultaneously analyse large multi-omics datasets integrated with metadata such as patient and clinical information.Results We developed VDJView, which permits the simultaneous or independent analysis and visualisation of gene expression, immune receptors, and clinical metadata of both T and B cells. This tool is implemented as an easy-to-use R shiny web-application, which integrates numerous gene expression and TCR analysis tools, and accepts data from plate-based sorted or high-throughput single cell platforms. We utilised VDJView to analyse several 10X scRNA-seq datasets, including a recent dataset of 150,000 CD8+ T cells with available gene expression, TCR sequences, quantification of 15 surface proteins, and 44 antigen specificities (across viruses, cancer, and self-antigens). We performed quality control, filtering of tetramer non-specific cells, clustering, random sampling and hypothesis testing to discover antigen specific gene signatures which were associated with immune cell differentiation states and clonal expansion across the pathogen specific T cells. We also analysed 563 single cells (plate-based sorted) obtained from 11 subjects, revealing clonally expanded T and B cells across primary cancer tissues and metastatic lymph-node. These immune cells clustered with distinct gene signatures according to the breast cancer molecular subtype. VDJView has been tested in lab meetings and peer-to-peer discussions, showing effective data generation and discussion without the need to consult bioinformaticians.Conclusions VDJView enables researchers without profound bioinformatics skills to analyse immune scRNA-seq data, integrating and visualising this with clonality and metadata profiles, thus accelerating the process of hypothesis testing, data interpretation and discovery of cellular heterogeneity. VDJView is freely available at https://bitbucket.org/kirbyvisp/vdjview .


2020 ◽  
Vol 5 (51) ◽  
pp. eabb1455 ◽  
Author(s):  
Kazuko Miyazaki ◽  
Hitomi Watanabe ◽  
Genki Yoshikawa ◽  
Kenian Chen ◽  
Reiko Hidaka ◽  
...  

Cell type–specific gene expression is driven by the interplay between lineage-specific transcription factors and cis-regulatory elements to which they bind. Adaptive immunity relies on RAG-mediated assembly of T cell receptor (TCR) and immunoglobulin (Ig) genes. Although Rag1 and Rag2 expression is largely restricted to adaptive lymphoid lineage cells, it remains unclear how Rag gene expression is regulated in a cell lineage–specific manner. Here, we identified three distinct cis-regulatory elements, a T cell lineage–specific enhancer (R-TEn) and the two B cell–specific elements, R1B and R2B. By generating mice lacking either R-TEn or R1B and R2B, we demonstrate that these distinct sets of regulatory elements drive the expression of Rag genes in developing T and B cells. What these elements have in common is their ability to bind the transcription factor E2A. By generating a mouse strain that carries a mutation within the E2A binding site of R-TEn, we demonstrate that recruitment of E2A to this site is essential for orchestrating changes in chromatin conformation that drive expression of Rag genes in T cells. By mapping cis-regulatory elements and generating multiple mouse strains lacking distinct enhancer elements, we demonstrate expression of Rag genes in developing T and B cells to be driven by distinct sets of E2A-dependent cis-regulatory modules.


2021 ◽  
Author(s):  
Zi-Hang Wen ◽  
Jeremy L. Langsam ◽  
Lu Zhang ◽  
Wenjun Shen ◽  
Xin Zhou

AbstractSingle-cell RNA-seq (scRNA-seq) offers opportunities to study gene expression of tens of thousands of single cells simultaneously, to investigate cell-to-cell variation, and to reconstruct cell-type-specific gene regulatory networks. Recovering dropout events in a sparse gene expression matrix for scRNA-seq data is a long-standing matrix completion problem. We introduce Bfimpute, a Bayesian factorization imputation algorithm that reconstructs two latent gene and cell matrices to impute final gene expression matrix within each cell group, with or without the aid of cell type labels or bulk data. Bfimpute achieves better accuracy than other six publicly notable scRNA-seq imputation methods on simulated and real scRNA-seq data, as measured by several different evaluation metrics. Bfimpute can also flexibly integrate any gene or cell related information that users provide to increase the performance. Availability: Bfimpute is implemented in R and is freely available at https://github.com/maiziezhoulab/Bfimpute.


2019 ◽  
Author(s):  
Yiliang Zhang ◽  
Kexuan Liang ◽  
Molei Liu ◽  
Yue Li ◽  
Hao Ge ◽  
...  

AbstractSingle-cell RNA sequencing technologies are widely used in recent years as a powerful tool allowing the observation of gene expression at the resolution of single cells. Two of the major challenges in scRNA-seq data analysis are dropout events and batch effects. The inflation of zero(dropout rate) varies substantially across single cells. Evidence has shown that technical noise, including batch effects, explains a notable proportion of this cell-to-cell variation. To capture biological variation, it is necessary to quantify and remove technical variation. Here, we introduce SCRIBE (Single-Cell Recovery Imputation with Batch Effects), a principled framework that imputes dropout events and corrects batch effects simultaneously. We demonstrate, through real examples, that SCRIBE outperforms existing scRNA-seq data analysis tools in recovering cell-specific gene expression patterns, removing batch effects and retaining biological variation across cells. Our software is freely available online at https://github.com/YiliangTracyZhang/SCRIBE.


2021 ◽  
Vol 7 (1) ◽  
Author(s):  
Wanlin Jin ◽  
Qi Yang ◽  
Yuyao Peng ◽  
Chengkai Yan ◽  
Yi Li ◽  
...  

AbstractMyasthenia gravis (MG) is a rare autoimmune disease. Although the impact of immune cell disorder in MG has been extensively studied, little is known about the transcriptomes of individual cells. Here, we assessed the transcriptional profiles of 39,243 cells by single-cell sequencing and identified 13 major cell clusters, along with 39 subgroups of cells derived from patients with new-onset myasthenia gravis and healthy controls. We found that B cells, CD4+ T cells, and monocytes exhibited more heterogeneity in MG patients. CD4+ T cells were expanded in MG patients. We reclustered B cells and CD4+ T cells, and predict their essential regulators. Further analyses demonstrated that B cells in MG exhibited higher transcriptional activity towards plasma cell differentiation, CD4+ T cell subsets were unbalanced, and inflammatory pathways of monocytes were highly activated. Notably, we discovered a disease-relevant subgroup, CD180− B cells. Increased CD180− B cells in MG are indicative of a high IgG composition and were associated with disease activity and the anti-AChR antibody. Together, our data further the understanding of the cellular heterogeneity involved in the pathogenesis of MG and provide large cell-type-specific markers for subsequent research.


2021 ◽  
Vol 7 (1) ◽  
Author(s):  
Fengchan Li ◽  
Kunmin Yan ◽  
Lili Wu ◽  
Zhong Zheng ◽  
Yun Du ◽  
...  

AbstractDisturbed blood flow (d-flow) has been known to induce changes of the cells in the arterial wall, increasing the risk of atherosclerosis. However, the heterogeneity of the vascular cell populations under d-flow remains less understood. To generate d-flow in vivo, partial carotid artery ligation (PCL) was performed. Seven days after ligation, single-cell RNA sequencing of nine left carotid arteries (LCA) from the PCL group (10,262 cells) or control group (14,580 cells) was applied and a single-cell atlas of gene expression was constructed. The integrated analysis identified 15 distinct carotid cell clusters, including 10 d-flow-relevant subpopulations. Among endothelial cells, at least four subpopulations were identified, including Klk8hi ECs, Lrp1hi ECs, Dkk2hi ECs, and Cd36hi ECs. Analysis of GSVA and single-cell trajectories indicated that the previously undescribed Dkk2hi ECs subpopulation was mechanosensitive and potentially transformed from Klk8hi ECs under d-flow. D-flow-induced Spp1hi VSMCs subpopulation that appeared to be endowed with osteoblast differentiation, suggesting a role in arterial stiffness. Among the infiltrating cell subpopulations, Trem2hi Mφ, Birc5hi Mφ, DCs, CD4+ T cells, CXCR6+ T cells, NK cells, and granulocytes were identified under d-flow. Of note, the novel Birc5hi Mφ was identified as a potential contributor to the accumulation of macrophages in atherosclerosis. Finally, Dkk2hi ECs, and Cd36hi ECs were also found in the proatherosclerotic area of the aorta where the d-flow occurs. In conclusion, we presented a comprehensive single-cell atlas of all cells in the carotid artery under d-flow, identified previously unrecognized cell subpopulations and their gene expression signatures, and suggested their specialized functions.


2020 ◽  
Author(s):  
Lin Li ◽  
Hao Dai ◽  
Zhaoyuan Fang ◽  
Luonan Chen

AbstractThe rapid advancement of single cell technologies has shed new light on the complex mechanisms of cellular heterogeneity. However, compared with bulk RNA sequencing (RNA-seq), single-cell RNA-seq (scRNA-seq) suffers from higher noise and lower coverage, which brings new computational difficulties. Based on statistical independence, cell-specific network (CSN) is able to quantify the overall associations between genes for each cell, yet suffering from a problem of overestimation related to indirect effects. To overcome this problem, we propose the “conditional cell-specific network” (CCSN) method, which can measure the direct associations between genes by eliminating the indirect associations. CCSN can be used for cell clustering and dimension reduction on a network basis of single cells. Intuitively, each CCSN can be viewed as the transformation from less “reliable” gene expression to more “reliable” gene-gene associations in a cell. Based on CCSN, we further design network flow entropy (NFE) to estimate the differentiation potency of a single cell. A number of scRNA-seq datasets were used to demonstrate the advantages of our approach: (1) one direct association network for one cell; (2) most existing scRNA-seq methods designed for gene expression matrices are also applicable to CCSN-transformed degree matrices; (3) CCSN-based NFE helps resolving the direction of differentiation trajectories by quantifying the potency of each cell. CCSN is publicly available at http://sysbio.sibcb.ac.cn/cb/chenlab/soft/CCSN.zip.


Sign in / Sign up

Export Citation Format

Share Document