1
|
Moriel N, Memet E, Nitzan M. Optimal sequencing budget allocation for trajectory reconstruction of single cells. Bioinformatics 2024; 40:i446-i452. [PMID: 38940162 PMCID: PMC11211845 DOI: 10.1093/bioinformatics/btae258] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/29/2024] Open
Abstract
BACKGROUND Charting cellular trajectories over gene expression is key to understanding dynamic cellular processes and their underlying mechanisms. While advances in single-cell RNA-sequencing technologies and computational methods have pushed forward the recovery of such trajectories, trajectory inference remains a challenge due to the noisy, sparse, and high-dimensional nature of single-cell data. This challenge can be alleviated by increasing either the number of cells sampled along the trajectory (breadth) or the sequencing depth, i.e. the number of reads captured per cell (depth). Generally, these two factors are coupled due to an inherent breadth-depth tradeoff that arises when the sequencing budget is constrained due to financial or technical limitations. RESULTS Here we study the optimal allocation of a fixed sequencing budget to optimize the recovery of trajectory attributes. Empirical results reveal that reconstruction accuracy of internal cell structure in expression space scales with the logarithm of either the breadth or depth of sequencing. We additionally observe a power law relationship between the optimal number of sampled cells and the corresponding sequencing budget. For linear trajectories, non-monotonicity in trajectory reconstruction across the breadth-depth tradeoff can impact downstream inference, such as expression pattern analysis along the trajectory. We demonstrate these results for five single-cell RNA-sequencing datasets encompassing differentiation of embryonic stem cells, pancreatic beta cells, hepatoblast and multipotent hematopoietic cells, as well as induced reprogramming of embryonic fibroblasts into neurons. By addressing the challenges of single-cell data, our study offers insights into maximizing the efficiency of cellular trajectory analysis through strategic allocation of sequencing resources.
Collapse
Affiliation(s)
- Noa Moriel
- School of Computer Science and Engineering, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| | - Edvin Memet
- Department of Physics, Harvard University, Cambridge, MA 02138, United States
| | - Mor Nitzan
- School of Computer Science and Engineering, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
- Racah Institute of Physics, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
- Faculty of Medicine, The Hebrew University of Jerusalem, Jerusalem 9112102, Israel
| |
Collapse
|
2
|
Kumari P, Kaur M, Dindhoria K, Ashford B, Amarasinghe SL, Thind AS. Advances in long-read single-cell transcriptomics. Hum Genet 2024:10.1007/s00439-024-02678-x. [PMID: 38787419 DOI: 10.1007/s00439-024-02678-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2023] [Accepted: 05/07/2024] [Indexed: 05/25/2024]
Abstract
Long-read single-cell transcriptomics (scRNA-Seq) is revolutionizing the way we profile heterogeneity in disease. Traditional short-read scRNA-Seq methods are limited in their ability to provide complete transcript coverage, resolve isoforms, and identify novel transcripts. The scRNA-Seq protocols developed for long-read sequencing platforms overcome these limitations by enabling the characterization of full-length transcripts. Long-read scRNA-Seq techniques initially suffered from comparatively poor accuracy compared to short read scRNA-Seq. However, with improvements in accuracy, accessibility, and cost efficiency, long-reads are gaining popularity in the field of scRNA-Seq. This review details the advances in long-read scRNA-Seq, with an emphasis on library preparation protocols and downstream bioinformatics analysis tools.
Collapse
Affiliation(s)
- Pallawi Kumari
- Institute of Microbial Technology, Council of Scientific and Industrial Research, Chandigarh, India
| | - Manmeet Kaur
- Institute of Microbial Technology, Council of Scientific and Industrial Research, Chandigarh, India
| | - Kiran Dindhoria
- Institute of Microbial Technology, Council of Scientific and Industrial Research, Chandigarh, India
| | - Bruce Ashford
- Illawarra Shoalhaven Local Health District (ISLHD), NSW Health, Wollongong, NSW, Australia
| | - Shanika L Amarasinghe
- Monash Biomedical Discovery Institute, Monash University, Clayton, VIC, 3800, Australia
- Walter and Eliza Hall Institute of Medical Research, 1G, Royal Parade, Parkville, VIC, 3025, Australia
| | - Amarinder Singh Thind
- Illawarra Shoalhaven Local Health District (ISLHD), NSW Health, Wollongong, NSW, Australia.
- The School of Chemistry and Molecular Bioscience (SCMB), University of Wollongong, Loftus St, Wollongong, NSW, 2500, Australia.
| |
Collapse
|
3
|
Joe S, Park JL, Kim J, Kim S, Park JH, Yeo MK, Lee D, Yang JO, Kim SY. Comparison of structural variant callers for massive whole-genome sequence data. BMC Genomics 2024; 25:318. [PMID: 38549092 PMCID: PMC10976732 DOI: 10.1186/s12864-024-10239-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2023] [Accepted: 03/18/2024] [Indexed: 04/01/2024] Open
Abstract
BACKGROUND Detecting structural variations (SVs) at the population level using next-generation sequencing (NGS) requires substantial computational resources and processing time. Here, we compared the performances of 11 SV callers: Delly, Manta, GridSS, Wham, Sniffles, Lumpy, SvABA, Canvas, CNVnator, MELT, and INSurVeyor. These SV callers have been recently published and have been widely employed for processing massive whole-genome sequencing datasets. We evaluated the accuracy, sequence depth, running time, and memory usage of the SV callers. RESULTS Notably, several callers exhibited better calling performance for deletions than for duplications, inversions, and insertions. Among the SV callers, Manta identified deletion SVs with better performance and efficient computing resources, and both Manta and MELT demonstrated relatively good precision regarding calling insertions. We confirmed that the copy number variation callers, Canvas and CNVnator, exhibited better performance in identifying long duplications as they employ the read-depth approach. Finally, we also verified the genotypes inferred from each SV caller using a phased long-read assembly dataset, and Manta showed the highest concordance in terms of the deletions and insertions. CONCLUSIONS Our findings provide a comprehensive understanding of the accuracy and computational efficiency of SV callers, thereby facilitating integrative analysis of SV profiles in diverse large-scale genomic datasets.
Collapse
Grants
- NRF-2020M3E5D708517212, 2020M3A9I6A0103605713 Ministry of Science and ICT, South Korea
- NRF-2020M3E5D708517212, 2020M3A9I6A0103605713 Ministry of Science and ICT, South Korea
- NRF-2020M3E5D708517212, 2020M3A9I6A0103605713 Ministry of Science and ICT, South Korea
- NRF-2020M3E5D708517212, 2020M3A9I6A0103605713 Ministry of Science and ICT, South Korea
- NRF-2020M3E5D708517212, 2020M3A9I6A0103605713 Ministry of Science and ICT, South Korea
- NRF-2020M3E5D708517212, 2020M3A9I6A0103605713 Ministry of Science and ICT, South Korea
- NRF-2020M3E5D708517212, 2020M3A9I6A0103605713 Ministry of Science and ICT, South Korea
- NRF-2020M3E5D708517212, 2020M3A9I6A0103605713 Ministry of Science and ICT, South Korea
- NTIS-1711170620 KRIBB Research Initiative Program
- NTIS-1711170620 KRIBB Research Initiative Program
- NTIS-1711170620 KRIBB Research Initiative Program
- NTIS-1711170620 KRIBB Research Initiative Program
- NTIS-1711170620 KRIBB Research Initiative Program
- NTIS-1711170620 KRIBB Research Initiative Program
- NTIS-1711170620 KRIBB Research Initiative Program
Collapse
Affiliation(s)
- Soobok Joe
- Korea Bioinformation Center (KOBIC), Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea
| | - Jong-Lyul Park
- Aging Convergence Research Center, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea
- Department of Functional Genomics, University of Science and Technology (UST), 34113, Daejeon, Republic of Korea
| | - Jun Kim
- Department of Convergent Bioscience and Informatics, College of Bioscience and Biotechnology, Chungnam National University, Daejeon, 34134, Republic of Korea
| | - Sangok Kim
- Korea Bioinformation Center (KOBIC), Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea
| | - Ji-Hwan Park
- Korea Bioinformation Center (KOBIC), Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea
- Department of Bioscience, University of Science and Technology (UST), Daejeon, 34113, Republic of Korea
| | - Min-Kyung Yeo
- Department of Pathology, Chungnam National University School of Medicine, Daejeon, 35015, Republic of Korea
| | - Dongyoon Lee
- Korea Bioinformation Center (KOBIC), Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea
| | - Jin Ok Yang
- Korea Bioinformation Center (KOBIC), Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea.
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea.
| | - Seon-Young Kim
- Korea Bioinformation Center (KOBIC), Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, 34141, Republic of Korea.
- Department of Bioscience, University of Science and Technology (UST), Daejeon, 34113, Republic of Korea.
| |
Collapse
|
4
|
Calia GP, Chen X, Zuckerman B, Weinberger LS. Comparative analysis between single-cell RNA-seq and single-molecule RNA FISH indicates that the pyrimidine nucleobase idoxuridine (IdU) globally amplifies transcriptional noise. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.03.14.532632. [PMID: 36993609 PMCID: PMC10055090 DOI: 10.1101/2023.03.14.532632] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 04/29/2023]
Abstract
Stochastic fluctuations (noise) in transcription generate substantial cell-to-cell variability, but the physiological roles of noise have remained difficult to determine in the absence of generalized noise-modulation approaches. Previous single-cell RNA-sequencing (scRNA-seq) suggested that the pyrimidine-base analog (5'-iodo-2'-deoxyuridine, IdU) could generally amplify noise without substantially altering mean-expression levels but scRNA-seq technical drawbacks potentially obscured the penetrance of IdU-induced transcriptional noise amplification. Here we quantify global-vs.-partial penetrance of IdU-induced noise amplification by assessing scRNA-seq data using numerous normalization algorithms and directly quantifying noise using single-molecule RNA FISH (smFISH) for a panel of genes from across the transcriptome. Alternate scRNA-seq analyses indicate IdU-induced noise amplification for ~90% of genes, and smFISH data verified noise amplification for ~90% of tested genes. Collectively, this analysis indicates which scRNA-seq algorithms are appropriate for quantifying noise and argues that IdU is a globally penetrant noise-enhancer molecule that could enable investigations of the physiological impacts of transcriptional noise.
Collapse
Affiliation(s)
- Giuliana P. Calia
- Gladstone|UCSF Center for Cell Circuitry, University of California, San Francisco, CA 94158
| | - Xinyue Chen
- Gladstone|UCSF Center for Cell Circuitry, University of California, San Francisco, CA 94158
| | - Binyamin Zuckerman
- Gladstone|UCSF Center for Cell Circuitry, University of California, San Francisco, CA 94158
| | - Leor S. Weinberger
- Gladstone|UCSF Center for Cell Circuitry, University of California, San Francisco, CA 94158
- Department of Biochemistry and Biophysics, University of California, San Francisco, CA 94158
- Department of Pharmaceutical Chemistry, University of California, San Francisco, CA 94158
| |
Collapse
|
5
|
Masood F, Chang JH, Akbar A, Song A, Hu WY, Azar DT, Rosenblatt MI. Therapeutic Strategies for Restoring Perturbed Corneal Epithelial Homeostasis in Limbal Stem Cell Deficiency: Current Trends and Future Directions. Cells 2022; 11:3247. [PMID: 36291115 PMCID: PMC9600167 DOI: 10.3390/cells11203247] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2022] [Revised: 10/11/2022] [Accepted: 10/12/2022] [Indexed: 02/03/2023] Open
Abstract
Limbal stem cells constitute an important cell population required for regeneration of the corneal epithelium. If insults to limbal stem cells or their niche are sufficiently severe, a disease known as limbal stem cell deficiency occurs. In the absence of functioning limbal stem cells, vision-compromising conjunctivalization of the corneal epithelium occurs, leading to opacification, inflammation, neovascularization, and chronic scarring. Limbal stem cell transplantation is the standard treatment for unilateral cases of limbal stem cell deficiency, but bilateral cases require allogeneic transplantation. Herein we review the current therapeutic utilization of limbal stem cells. We also describe several limbal stem cell markers that impact their phenotype and function and discuss the possibility of modulating limbal stem cells and other sources of stem cells to facilitate the development of novel therapeutic interventions. We finally consider several hurdles for widespread adoption of these proposed methodologies and discuss how they can be overcome to realize vision-restoring interventions.
Collapse
Affiliation(s)
- Faisal Masood
- Department of Ophthalmology and Visual Sciences, Illinois Eye and Ear Infirmary, College of Medicine, University of Illinois at Chicago, Chicago, IL 60612, USA
| | - Jin-Hong Chang
- Department of Ophthalmology and Visual Sciences, Illinois Eye and Ear Infirmary, College of Medicine, University of Illinois at Chicago, Chicago, IL 60612, USA
| | - Anosh Akbar
- Department of Ophthalmology and Visual Sciences, Illinois Eye and Ear Infirmary, College of Medicine, University of Illinois at Chicago, Chicago, IL 60612, USA
| | - Amy Song
- Department of Ophthalmology and Visual Sciences, Illinois Eye and Ear Infirmary, College of Medicine, University of Illinois at Chicago, Chicago, IL 60612, USA
| | - Wen-Yang Hu
- Department of Urology, University of Illinois at Chicago, Chicago, IL 60612, USA
| | - Dimitri T. Azar
- Department of Ophthalmology and Visual Sciences, Illinois Eye and Ear Infirmary, College of Medicine, University of Illinois at Chicago, Chicago, IL 60612, USA
| | - Mark I. Rosenblatt
- Department of Ophthalmology and Visual Sciences, Illinois Eye and Ear Infirmary, College of Medicine, University of Illinois at Chicago, Chicago, IL 60612, USA
| |
Collapse
|
6
|
Qi J, Sheng Q, Zhou Y, Hua J, Xiao S, Jin S. scMTD: a statistical multidimensional imputation method for single-cell RNA-seq data leveraging transcriptome dynamic information. Cell Biosci 2022; 12:142. [PMID: 36056412 PMCID: PMC9440561 DOI: 10.1186/s13578-022-00886-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2022] [Accepted: 08/17/2022] [Indexed: 11/17/2022] Open
Abstract
Background Single-cell RNA sequencing (scRNA-seq) provides a powerful tool to capture transcriptomes at single-cell resolution. However, dropout events distort the gene expression levels and underlying biological signals, misleading the downstream analysis of scRNA-seq data. Results We develop a statistical model-based multidimensional imputation algorithm, scMTD, that identifies local cell neighbors and specific gene co-expression networks based on the pseudo-time of cells, leveraging information on cell-level, gene-level, and transcriptome dynamic to recover scRNA-seq data. Compared with the state-of-the-art imputation methods through several real-data-based analytical experiments, scMTD effectively recovers biological signals of transcriptomes and consistently outperforms the other algorithms in improving FISH validation, trajectory inference, differential expression analysis, clustering analysis, and identification of cell types. Conclusions scMTD maintains the gene expression characteristics, enhances the clustering of cell subpopulations, assists the study of gene expression dynamics, contributes to the discovery of rare cell types, and applies to both UMI-based and non-UMI-based data. Overall, scMTD’s reliability, applicability, and scalability make it a promising imputation approach for scRNA-seq data. Supplementary Information The online version contains supplementary material available at 10.1186/s13578-022-00886-4.
Collapse
|
7
|
A human adipose tissue cell-type transcriptome atlas. Cell Rep 2022; 40:111046. [PMID: 35830816 DOI: 10.1016/j.celrep.2022.111046] [Citation(s) in RCA: 30] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2021] [Revised: 04/29/2022] [Accepted: 06/13/2022] [Indexed: 12/19/2022] Open
Abstract
The importance of defining cell-type-specific genes is well acknowledged. Technological advances facilitate high-resolution sequencing of single cells, but practical challenges remain. Adipose tissue is composed primarily of adipocytes, large buoyant cells requiring extensive, artefact-generating processing for separation and analysis. Thus, adipocyte data are frequently absent from single-cell RNA sequencing (scRNA-seq) datasets, despite being the primary functional cell type. Here, we decipher cell-type-enriched transcriptomes from unfractionated human adipose tissue RNA-seq data. We profile all major constituent cell types, using 527 visceral adipose tissue (VAT) or 646 subcutaneous adipose tissue (SAT) samples, identifying over 2,300 cell-type-enriched transcripts. Sex-subset analysis uncovers a panel of male-only cell-type-enriched genes. By resolving expression profiles of genes differentially expressed between SAT and VAT, we identify mesothelial cells as the primary driver of this variation. This study provides an accessible method to profile cell-type-enriched transcriptomes using bulk RNA-seq, generating a roadmap for adipose tissue biology.
Collapse
|
8
|
Patel S, Ranadive I, Buch P, Khaire K, Balakrishnan S. De Novo Transcriptome Sequencing and Analysis of Differential Gene Expression among Various Stages of Tail Regeneration in Hemidactylus flaviviridis. J Dev Biol 2022; 10:jdb10020024. [PMID: 35735915 PMCID: PMC9225231 DOI: 10.3390/jdb10020024] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2022] [Revised: 06/08/2022] [Accepted: 06/10/2022] [Indexed: 11/16/2022] Open
Abstract
Across the animal kingdom, lizards are the only amniotes capable of regenerating their lost tail through epimorphosis. Of the many reptiles, the northern house gecko, Hemidactylus flaviviridis, is an excellent model system that is used for understanding the mechanism of epimorphic regeneration. A stage-specific transcriptome profile was generated in the current study following an autotomized tail with the HiSeq2500 platform. The reads obtained from de novo sequencing were filtered and high-quality reads were considered for gene ontology (GO) annotation and pathway analysis. Millions of reads were recorded for each stage upon de novo assembly. Up and down-regulated transcripts were categorized for early blastema (EBL), blastema (BL) and differentiation (DF) stages compared to the normal tail (NT) by differential gene expression analysis. The transcripts from developmentally significant pathways such as FGF, Wnt, Shh and TGF-β/BMP were present during tail regeneration. Additionally, differential expression of transcripts was recorded from biological processes, namely inflammation, cell proliferation, apoptosis and cell migration. Overall, the study reveals the stage-wise transcriptome analysis in conjunction with cellular processes as well as molecular signaling pathways during lizard tail regeneration. The knowledge obtained from the data can be extrapolated to configure regenerative responses in other amniotes, including humans, upon loss of a complex organ.
Collapse
|
9
|
Kim H, Kamm RD, Vunjak-Novakovic G, Wu JC. Progress in multicellular human cardiac organoids for clinical applications. Cell Stem Cell 2022; 29:503-514. [PMID: 35395186 PMCID: PMC9352318 DOI: 10.1016/j.stem.2022.03.012] [Citation(s) in RCA: 42] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Advances in self-organizing cardiac organoids to recapitulate human cardiogenesis have provided a powerful tool for unveiling human cardiac development, studying cardiovascular diseases, testing drugs, and transplantation. Here, we highlight the recent remarkable progress on multicellular cardiac organoids and review the current status for their practical applications. We then introduce key readouts and tools for assessing cardiac organoids for clinical applications, address major challenges, and provide suggestions for each assessment method. Lastly, we discuss the current limitations of cardiac organoids as miniature models of the human heart and suggest a direction for moving forward toward building the mini-heart from cardiac organoids.
Collapse
Affiliation(s)
- Hyeonyu Kim
- Stanford Cardiovascular Institute, Stanford, CA 94305, USA
| | - Roger D Kamm
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA; Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Gordana Vunjak-Novakovic
- Department of Biomedical Engineering, Columbia University, New York, NY 10027, USA; Department of Medicine, Columbia University, New York, NY 10032, USA
| | - Joseph C Wu
- Stanford Cardiovascular Institute, Stanford, CA 94305, USA; Department of Medicine, Stanford University, Stanford, CA 94305, USA.
| |
Collapse
|
10
|
A novel method for single-cell data imputation using subspace regression. Sci Rep 2022; 12:2697. [PMID: 35177662 PMCID: PMC8854597 DOI: 10.1038/s41598-022-06500-4] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2021] [Accepted: 01/27/2022] [Indexed: 12/13/2022] Open
Abstract
Recent advances in biochemistry and single-cell RNA sequencing (scRNA-seq) have allowed us to monitor the biological systems at the single-cell resolution. However, the low capture of mRNA material within individual cells often leads to inaccurate quantification of genetic material. Consequently, a significant amount of expression values are reported as missing, which are often referred to as dropouts. To overcome this challenge, we develop a novel imputation method, named single-cell Imputation via Subspace Regression (scISR), that can reliably recover the dropout values of scRNA-seq data. The scISR method first uses a hypothesis-testing technique to identify zero-valued entries that are most likely affected by dropout events and then estimates the dropout values using a subspace regression model. Our comprehensive evaluation using 25 publicly available scRNA-seq datasets and various simulation scenarios against five state-of-the-art methods demonstrates that scISR is better than other imputation methods in recovering scRNA-seq expression profiles via imputation. scISR consistently improves the quality of cluster analysis regardless of dropout rates, normalization techniques, and quantification schemes. The source code of scISR can be found on GitHub at https://github.com/duct317/scISR.
Collapse
|
11
|
Singh R, Zhang F, Li Q. Assessing reproducibility of high-throughput experiments in the case of missing data. Stat Med 2022; 41:1884-1899. [PMID: 35178743 PMCID: PMC9039958 DOI: 10.1002/sim.9334] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2021] [Revised: 01/08/2022] [Accepted: 01/14/2022] [Indexed: 11/11/2022]
Abstract
High-throughput experiments are an essential part of modern biological and biomedical research. The outcomes of high-throughput biological experiments often have a lot of missing observations due to signals below detection levels. For example, most single-cell RNA-seq (scRNA-seq) protocols experience high levels of dropout due to the small amount of starting material, leading to a majority of reported expression levels being zero. Though missing data contain information about reproducibility, they are often excluded in the reproducibility assessment, potentially generating misleading assessments. In this article, we develop a regression model to assess how the reproducibility of high-throughput experiments is affected by the choices of operational factors (eg, platform or sequencing depth) when a large number of measurements are missing. Using a latent variable approach, we extend correspondence curve regression, a recently proposed method for assessing the effects of operational factors to reproducibility, to incorporate missing values. Using simulations, we show that our method is more accurate in detecting differences in reproducibility than existing measures of reproducibility. We illustrate the usefulness of our method using a single-cell RNA-seq dataset collected on HCT116 cells. We compare the reproducibility of different library preparation platforms and study the effect of sequencing depth on reproducibility, thereby determining the cost-effective sequencing depth that is required to achieve sufficient reproducibility.
Collapse
Affiliation(s)
- Roopali Singh
- Department of Statistics, Pennsylvania State University, University Park, Pennsylvania, USA
| | - Feipeng Zhang
- School of Economics and Finance, Xi'an Jiaotong University, Xi'an, China
| | - Qunhua Li
- Department of Statistics, Pennsylvania State University, University Park, Pennsylvania, USA
| |
Collapse
|
12
|
Mapping the expression of an ANK3 isoform associated with bipolar disorder in the human brain. Transl Psychiatry 2022; 12:45. [PMID: 35091539 PMCID: PMC8799726 DOI: 10.1038/s41398-022-01784-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/03/2021] [Revised: 12/09/2021] [Accepted: 01/04/2022] [Indexed: 11/08/2022] Open
Abstract
The gene ankyrin-3 (ANK3) has been consistently associated with bipolar disorder (BD) in several genome-wide association studies (GWAS). The exact molecular mechanisms underlying this genetic association remain unknown. The discovery of a loss-of-function variant (rs41283526*G) in an alternatively spliced exon (ENSE00001786716) with a protective effect, suggested that elevated expression of this particular isoform could be a risk factor for developing the disorder. We developed a novel approach for measuring the expression level of all splice forms at a challenging genetic locus using a combination of droplet digital PCR and high-throughput sequencing of indexed PCR amplicons. The combined method was performed on a large collection of 568 postmortem brain samples of BD and SCZ cases and controls. We also studied the expression of the splice forms in a child-development cohort of 41 healthy males. We found that our approach can quantify the splice forms in brain samples, although with less precision than ddPCR. We detected highly significant differences in expression of splice forms and transcription start sites between brain regions, notably with higher expression of the BD-associated isoform in the corpus callosum compared to frontal tissue (mean fold change = 1.80, p < 1e-4). Although the patients in our sample expressed the BD-associated splice form at a similar level to controls, adolescents in our child-development cohort had a clearly higher expression level than younger children (mean fold change = 1.97, p = 5e-3). These results suggest that this ANK3 splice form may play a role in the myelin maturation of the human brain.
Collapse
|
13
|
Zheng J, Takahashi S, Zhou Y, Cheatham MA. Prestin and electromotility may serve multiple roles in cochlear outer hair cells. Hear Res 2021; 423:108428. [PMID: 34987016 DOI: 10.1016/j.heares.2021.108428] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/11/2021] [Revised: 12/16/2021] [Accepted: 12/23/2021] [Indexed: 11/04/2022]
Abstract
Outer hair cells (OHCs) are innervated by both medial olivocochlear (MOC) efferents and type II afferents, which also innervate supporting cells to form a local neural network. It has also been demonstrated that prestin provides the molecular basis for OHC somatic electromotility, amplifying movements within the organ of Corti. Although not anticipated, early-onset OHC loss was found in two prestin transgenic mouse models that either lack prestin protein or lack electromotility. To uncover the molecular pathways that evoke OHC death, we profiled the coding transcriptome of OHCs from wildtype (WT), prestin-knockout (KO), and 499-knockin (KI) mice using single-cell RNA sequencing (scRNA-seq). scRNA-Seq transcriptomics and pathway analyses did not reveal common pathways associated with OHC loss observed in prestin-KO and 499-KI mice. Clustering enrichment analysis showed that increased gene expression in OHCs from prestin-KO mice was associated with lipid metabolic processes and cell death pathways. These mRNA profiles likely contribute to the OHC loss observed in prestin-KO mice and support the notion that prestin is also a structural protein, important for the normal plasma membrane compartmentalization that is essential to establish MOC efferent synapses. In contrast, the mRNA profile of OHCs from 499-KI mice did not provide a rational explanation of the early-onset OHC loss in this mutant. OHCs from 499-KI mice have normal plasma membrane compartmentalization and normal OHC-MOC contacts. However, 499 prestin lacks electromotility and appears to change the local neural network around OHCs, as more synaptic markers were found near neighboring supporting cells when compared to WT and prestin-KO mice. Thus, OHCs in prestin-KOs (no prestin protein, no electromotility) and 499-KIs (prestin protein present, no electromotility) may influence local neuronal networks in different ways. Collectively, our data suggest that prestin and its motile properties are important for OHC survival and the maintenance of local afferent/efferent circuits, as well as for its role in cochlear amplification.
Collapse
Affiliation(s)
- Jing Zheng
- Departments of Otolaryngology, Feinberg School of Medicine, Northwestern University, Chicago, IL; Communication Sciences and Disorders, School of Communication; The Knowles Hearing Center, Northwestern University, Evanston, IL.
| | - Satoe Takahashi
- Departments of Otolaryngology, Feinberg School of Medicine, Northwestern University, Chicago, IL
| | - Yingjie Zhou
- Communication Sciences and Disorders, School of Communication
| | - Mary Ann Cheatham
- Communication Sciences and Disorders, School of Communication; The Knowles Hearing Center, Northwestern University, Evanston, IL
| |
Collapse
|
14
|
Li H, Huang Q, Liu Y, Garmire LX. Single cell transcriptome research in human placenta. Reproduction 2021; 160:R155-R167. [PMID: 33112783 PMCID: PMC7707799 DOI: 10.1530/rep-20-0231] [Citation(s) in RCA: 41] [Impact Index Per Article: 13.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2020] [Accepted: 09/22/2020] [Indexed: 12/30/2022]
Abstract
Human placenta is a complex and heterogeneous organ interfacing between the mother and the fetus that supports fetal development. Alterations to placental structural components are associated with various pregnancy complications. To reveal the heterogeneity among various placenta cell types in normal and diseased placentas, as well as elucidate molecular interactions within a population of placental cells, a new genomics technology called single cell RNA-seq (or scRNA-seq) has been employed in the last couple of years. Here we review the principles of scRNA-seq technology, and summarize the recent human placenta studies at scRNA-seq level across gestational ages as well as in pregnancy complications, such as preterm birth and preeclampsia. We list the computational analysis platforms and resources available for the public use. Lastly, we discuss the future areas of interest for placenta single cell studies, as well as the data analytics needed to accomplish them.
Collapse
Affiliation(s)
- Hui Li
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan, USA
| | - Qianhui Huang
- Department of Biostatistics, School of Public Health, University of Michigan, Ann Arbor, Michigan, USA
| | - Yu Liu
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan, USA
| | - Lana X Garmire
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan, USA
| |
Collapse
|
15
|
High-throughput and single-cell T cell receptor sequencing technologies. Nat Methods 2021; 18:881-892. [PMID: 34282327 PMCID: PMC9345561 DOI: 10.1038/s41592-021-01201-8] [Citation(s) in RCA: 123] [Impact Index Per Article: 41.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2020] [Accepted: 06/07/2021] [Indexed: 02/06/2023]
Abstract
T cells express T cell receptors (TCRs) composed of somatically recombined TCRα and TCRβ chains, which mediate recognition of major histocompatibility complex (MHC)-antigen complexes and drive the antigen-specific adaptive immune response to pathogens and cancer. The TCR repertoire in each individual is highly diverse, which allows for recognition of a wide array of foreign antigens, but also presents a challenge in analyzing this response using conventional methods. Recent studies have developed high-throughput sequencing technologies to identify TCR sequences, analyze their antigen specificities using experimental and computational tools, and pair TCRs with transcriptional and epigenetic cell state phenotypes in single cells. In this Review, we highlight these technological advances and describe how they have been applied to discover fundamental insights into T cell-mediated immunity.
Collapse
|
16
|
Rocque B, Barbetta A, Singh P, Goldbeck C, Helou DG, Loh YHE, Ung N, Lee J, Akbari O, Emamaullee J. Creation of a Single Cell RNASeq Meta-Atlas to Define Human Liver Immune Homeostasis. Front Immunol 2021; 12:679521. [PMID: 34335581 PMCID: PMC8322955 DOI: 10.3389/fimmu.2021.679521] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2021] [Accepted: 06/28/2021] [Indexed: 01/16/2023] Open
Abstract
The liver is unique in both its ability to maintain immune homeostasis and in its potential for immune tolerance following solid organ transplantation. Single-cell RNA sequencing (scRNA seq) is a powerful approach to generate highly dimensional transcriptome data to understand cellular phenotypes. However, when scRNA data is produced by different groups, with different data models, different standards, and samples processed in different ways, it can be challenging to draw meaningful conclusions from the aggregated data. The goal of this study was to establish a method to combine ‘human liver’ scRNA seq datasets by 1) characterizing the heterogeneity between studies and 2) using the meta-atlas to define the dominant phenotypes across immune cell subpopulations in healthy human liver. Publicly available scRNA seq data generated from liver samples obtained from a combined total of 17 patients and ~32,000 cells were analyzed. Liver-specific immune cells (CD45+) were extracted from each dataset, and immune cell subpopulations (myeloid cells, NK and T cells, plasma cells, and B cells) were examined using dimensionality reduction (UMAP), differential gene expression, and ingenuity pathway analysis. All datasets co-clustered, but cell proportions differed between studies. Gene expression correlation demonstrated similarity across all studies, and canonical pathways that differed between datasets were related to cell stress and oxidative phosphorylation rather than immune-related function. Next, a meta-atlas was generated via data integration and compared against PBMC data to define gene signatures for each hepatic immune subpopulation. This analysis defined key features of hepatic immune homeostasis, with decreased expression across immunologic pathways and enhancement of pathways involved with cell death. This method for meta-analysis of scRNA seq data provides a novel approach to broadly define the features of human liver immune homeostasis. Specific pathways and cellular phenotypes described in this human liver immune meta-atlas provide a critical reference point for further study of immune mediated disease processes within the liver.
Collapse
Affiliation(s)
- Brittany Rocque
- Division of Abdominal Organ Transplantation and Hepatobiliary Surgery, Department of Surgery, Keck School of Medicine, University of Southern California, Los Angeles, CA, United States.,Eli and Edythe Broad Center for Regenerative Medicine and Stem Cell Research, Keck School of Medicine, University of Southern California, Los Angeles, CA, United States
| | - Arianna Barbetta
- Division of Abdominal Organ Transplantation and Hepatobiliary Surgery, Department of Surgery, Keck School of Medicine, University of Southern California, Los Angeles, CA, United States
| | - Pranay Singh
- Division of Abdominal Organ Transplantation and Hepatobiliary Surgery, Department of Surgery, Keck School of Medicine, University of Southern California, Los Angeles, CA, United States
| | - Cameron Goldbeck
- Division of Abdominal Organ Transplantation and Hepatobiliary Surgery, Department of Surgery, Keck School of Medicine, University of Southern California, Los Angeles, CA, United States
| | - Doumet Georges Helou
- Department of Molecular Microbiology and Immunology, Keck School of Medicine, University of Southern California, Los Angeles, CA, United States
| | - Yong-Hwee Eddie Loh
- Norris Medical Library, University of Southern California, Los Angeles, CA, United States
| | - Nolan Ung
- Lawrence J. Ellison Institute for Transformative Medicine, University of Southern California, Los Angeles, CA, United States
| | - Jerry Lee
- Lawrence J. Ellison Institute for Transformative Medicine, University of Southern California, Los Angeles, CA, United States.,Department of Medicine, Keck School of Medicine, University of Southern California, Los Angeles, CA, United States.,Department of Chemical Engineering and Materials Sciences, University of Southern California, Los Angeles, CA, United States
| | - Omid Akbari
- Department of Molecular Microbiology and Immunology, Keck School of Medicine, University of Southern California, Los Angeles, CA, United States
| | - Juliet Emamaullee
- Division of Abdominal Organ Transplantation and Hepatobiliary Surgery, Department of Surgery, Keck School of Medicine, University of Southern California, Los Angeles, CA, United States.,Department of Molecular Microbiology and Immunology, Keck School of Medicine, University of Southern California, Los Angeles, CA, United States
| |
Collapse
|
17
|
Qi J, Zhou Y, Zhao Z, Jin S. SDImpute: A statistical block imputation method based on cell-level and gene-level information for dropouts in single-cell RNA-seq data. PLoS Comput Biol 2021; 17:e1009118. [PMID: 34138847 PMCID: PMC8266063 DOI: 10.1371/journal.pcbi.1009118] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2020] [Revised: 07/08/2021] [Accepted: 05/27/2021] [Indexed: 11/19/2022] Open
Abstract
The single-cell RNA sequencing (scRNA-seq) technologies obtain gene expression at single-cell resolution and provide a tool for exploring cell heterogeneity and cell types. As the low amount of extracted mRNA copies per cell, scRNA-seq data exhibit a large number of dropouts, which hinders the downstream analysis of the scRNA-seq data. We propose a statistical method, SDImpute (Single-cell RNA-seq Dropout Imputation), to implement block imputation for dropout events in scRNA-seq data. SDImpute automatically identifies the dropout events based on the gene expression levels and the variations of gene expression across similar cells and similar genes, and it implements block imputation for dropouts by utilizing gene expression unaffected by dropouts from similar cells. In the experiments, the results of the simulated datasets and real datasets suggest that SDImpute is an effective tool to recover the data and preserve the heterogeneity of gene expression across cells. Compared with the state-of-the-art imputation methods, SDImpute improves the accuracy of the downstream analysis including clustering, visualization, and differential expression analysis. Single-cell RNA sequencing (scRNA-seq) allows researchers to analyze gene expression of thousands of single cells simultaneously. However, the low amount of extracted mRNA leads to a large number of dropout events, which introduce computational challenges and hinder downstream analysis of data. To address this problem, we developed SDImpute, a novel statistical method to recover the scRNA-seq data based on cell-level and gene-level information in this manuscript. The goal of our algorithm is to denoise the scRNA-seq data while maintaining the biological nature of gene expression. Combining SDImpute with the downstream analysis tools, we considered the matched bulk expression data and known cell labels of the scRNA-seq data as criteria to design experiments to validate the performance of our method in both simulated and real datasets. Moreover, we offer an R package with detailed instructions and an example input dataset. We hope that SDImpute will be beneficial to researchers to identify mechanisms underlying some biological processes by analysis of the scRNA-seq data.
Collapse
Affiliation(s)
- Jing Qi
- School of Mathematics, Harbin Institute of Technology, Harbin, P.R, China
| | - Yang Zhou
- School of Mathematics, Harbin Institute of Technology, Harbin, P.R, China
| | - Zicen Zhao
- School of Mathematics, Harbin Institute of Technology, Harbin, P.R, China
| | - Shuilin Jin
- School of Mathematics, Harbin Institute of Technology, Harbin, P.R, China
- * E-mail:
| |
Collapse
|
18
|
Nguyen TM, Jeevan JJ, Xu N, Chen JY. Polar Gini Curve: A Technique to Discover Gene Expression Spatial Patterns from Single-cell RNA-seq Data. GENOMICS, PROTEOMICS & BIOINFORMATICS 2021; 19:493-503. [PMID: 34958962 PMCID: PMC8864247 DOI: 10.1016/j.gpb.2020.09.006] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/05/2020] [Revised: 07/09/2020] [Accepted: 10/29/2020] [Indexed: 12/13/2022]
Abstract
In this work, we describe the development of Polar Gini Curve, a method for characterizing cluster markers by analyzing single-cell RNA sequencing (scRNA-seq) data. Polar Gini Curve combines the gene expression and the 2D coordinates ("spatial") information to detect patterns of uniformity in any clustered cells from scRNA-seq data. We demonstrate that Polar Gini Curve can help users characterize the shape and density distribution of cells in a particular cluster, which can be generated during routine scRNA-seq data analysis. To quantify the extent to which a gene is uniformly distributed in a cell cluster space, we combine two polar Gini curves (PGCs)-one drawn upon the cell-points expressing the gene (the "foreground curve") and the other drawn upon all cell-points in the cluster (the "background curve"). We show that genes with highly dissimilar foreground and background curves tend not to uniformly distributed in the cell cluster-thus having spatially divergent gene expression patterns within the cluster. Genes with similar foreground and background curves tend to uniformly distributed in the cell cluster-thus having uniform gene expression patterns within the cluster. Such quantitative attributes of PGCs can be applied to sensitively discover biomarkers across clusters from scRNA-seq data. We demonstrate the performance of the Polar Gini Curve framework in several simulation case studies. Using this framework to analyze a real-world neonatal mouse heart cell dataset, the detected biomarkers may characterize novel subtypes of cardiac muscle cells. The source code and data for Polar Gini Curve could be found at http://discovery.informatics.uab.edu/PGC/ or https://figshare.com/projects/Polar_Gini_Curve/76749.
Collapse
Affiliation(s)
- Thanh Minh Nguyen
- Informatics Institute, the University of Alabama at Birmingham, Birmingham, AL 35294, USA
| | - Jacob John Jeevan
- Informatics Institute, the University of Alabama at Birmingham, Birmingham, AL 35294, USA
| | - Nuo Xu
- Collat School of Business, the University of Alabama at Birmingham, Birmingham, AL 35294, USA
| | - Jake Y Chen
- Informatics Institute, the University of Alabama at Birmingham, Birmingham, AL 35294, USA.
| |
Collapse
|
19
|
Louie RH, Luciani F. Recent advances in single-cell multimodal analysis to study immune cells. Immunol Cell Biol 2021; 99:157-167. [PMID: 33314406 DOI: 10.1111/imcb.12432] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2020] [Revised: 12/09/2020] [Accepted: 12/09/2020] [Indexed: 12/11/2022]
Abstract
Recent advances in single-cell technologies have enabled the profiling of the genome, epigenome, transcriptome and proteome, along with temporal and spatial information of individual cells. These technologies have provided unique opportunities to understand mechanisms underpinning the immune system, such as characterizations of the molecular cell state, how the cell state evolves along its lineage and the impact of spatial location on cell state. In this review, we discuss how these mechanisms have been studied through recent advances in single-cell multimodal technologies.
Collapse
Affiliation(s)
- Raymond Hy Louie
- School of Medical Sciences, The Kirby Institute, University of New South Wales (UNSW), Sydney, NSW, Australia
| | - Fabio Luciani
- School of Medical Sciences, The Kirby Institute, University of New South Wales (UNSW), Sydney, NSW, Australia
| |
Collapse
|
20
|
Qi J, Zhou Y, Hua J, Zhang L, Bian J, Liu B, Zhao Z, Jin S. The scRNA-seq Expression Profiling of the Receptor ACE2 and the Cellular Protease TMPRSS2 Reveals Human Organs Susceptible to SARS-CoV-2 Infection. INTERNATIONAL JOURNAL OF ENVIRONMENTAL RESEARCH AND PUBLIC HEALTH 2021; 18:E284. [PMID: 33401657 PMCID: PMC7794913 DOI: 10.3390/ijerph18010284] [Citation(s) in RCA: 63] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/29/2020] [Revised: 12/21/2020] [Accepted: 12/22/2020] [Indexed: 01/08/2023]
Abstract
COVID-19 patients always develop multiple organ dysfunction syndromes other than lungs, suggesting the novel virus SARS-CoV-2 also invades other organs. Therefore, studying the viral susceptibility of other organs is important for a deeper understanding of viral pathogenesis. Angiotensin-converting enzyme II (ACE2) is the receptor protein of SARS-CoV-2, and TMPRSS2 promotes virus proliferation and transmission. We investigated the ACE2 and TMPRSS2 expression levels of cell types from 31 organs to evaluate the risk of viral infection using single-cell RNA sequencing (scRNA-seq) data. For the first time, we found that the gall bladder and fallopian tube are vulnerable to SARS-CoV-2 infection. Besides, the nose, heart, small intestine, large intestine, esophagus, brain, testis, and kidney are also identified to be high-risk organs with high expression levels of ACE2 and TMPRSS2. Moreover, the susceptible organs are grouped into three risk levels based on the ACE2 and TMPRSS2 expression. As a result, the respiratory system, digestive system, and urinary system are at the top-risk level for SARS-CoV-2 infection. This study provides evidence for SARS-CoV-2 infection in the human nervous system, digestive system, reproductive system, respiratory system, circulatory system, and urinary system using scRNA-seq data, which helps in the clinical diagnosis and treatment of patients.
Collapse
Affiliation(s)
| | | | | | | | | | | | | | - Shuilin Jin
- School of Mathematics, Harbin Institute of Technology, Harbin 150001, China; (J.Q.); (Y.Z.); (J.H.); (L.Z.); (J.B.); (B.L.); (Z.Z.)
| |
Collapse
|
21
|
Peng Y, Qiao H. The Application of Single-Cell RNA Sequencing in Mammalian Meiosis Studies. Front Cell Dev Biol 2021; 9:673642. [PMID: 34485276 PMCID: PMC8416306 DOI: 10.3389/fcell.2021.673642] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2021] [Accepted: 07/05/2021] [Indexed: 12/14/2022] Open
Abstract
Meiosis is a cellular division process that produces gametes for sexual reproduction. Disruption of complex events throughout meiosis, such as synapsis and homologous recombination, can lead to infertility and aneuploidy. To reveal the molecular mechanisms of these events, transcriptome studies of specific substages must be conducted. However, conventional methods, such as bulk RNA-seq and RT-qPCR, are not able to detect the transcriptional variations effectively and precisely, especially for identifying cell types and stages with subtle differences. In recent years, mammalian meiotic transcriptomes have been intensively studied at the single-cell level by using single-cell RNA-seq (scRNA-seq) approaches, especially through two widely used platforms, Smart-seq2 and Drop-seq. The scRNA-seq protocols along with their downstream analysis enable researchers to accurately identify cell heterogeneities and investigate meiotic transcriptomes at a higher resolution. In this review, we compared bulk RNA-seq and scRNA-seq to show the advantages of the scRNA-seq in meiosis studies; meanwhile, we also pointed out the challenges and limitations of the scRNA-seq. We listed recent findings from mammalian meiosis (male and female) studies where scRNA-seq applied. Next, we summarized the scRNA-seq analysis methods and the meiotic marker genes from spermatocytes and oocytes. Specifically, we emphasized the different features of the two scRNA-seq protocols (Smart-seq2 and Drop-seq) in the context of meiosis studies and discussed their strengths and weaknesses in terms of different research purposes. Finally, we discussed the future applications of scRNA-seq in the meiosis field.
Collapse
Affiliation(s)
- Yiheng Peng
- Department of Comparative Biosciences, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - Huanyu Qiao
- Department of Comparative Biosciences, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| |
Collapse
|
22
|
Ding J, Smith SL, Orozco G, Barton A, Eyre S, Martin P. Characterisation of CD4+ T-cell subtypes using single cell RNA sequencing and the impact of cell number and sequencing depth. Sci Rep 2020; 10:19825. [PMID: 33188258 PMCID: PMC7666207 DOI: 10.1038/s41598-020-76972-9] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2020] [Accepted: 10/14/2020] [Indexed: 12/31/2022] Open
Abstract
CD4+ T-cells represent a heterogeneous collection of specialised sub-types and are a key cell type in the pathogenesis of many diseases due to their role in the adaptive immune system. By investigating CD4+ T-cells at the single cell level, using RNA sequencing (scRNA-seq), there is the potential to identify specific cell states driving disease or treatment response. However, the impact of sequencing depth and cell numbers, two important factors in scRNA-seq, has not been determined for a complex cell population such as CD4+ T-cells. We therefore generated a high depth, high cell number dataset to determine the effect of reduced sequencing depth and cell number on the ability to accurately identify CD4+ T-cell subtypes. Furthermore, we investigated T-cell signatures under resting and stimulated conditions to assess cluster specific effects of stimulation. We found that firstly, cell number has a much more profound effect than sequencing depth on the ability to classify cells; secondly, this effect is greater when cells are unstimulated and finally, resting and stimulated samples can be combined to leverage additional power whilst still allowing differences between samples to be observed. While based on one individual, these results could inform future scRNA-seq studies to ensure the most efficient experimental design.
Collapse
Affiliation(s)
- James Ding
- Centre for Genetics and Genomics Versus Arthritis, Centre for Musculoskeletal Research, Manchester Academic Health Science Centre, The University of Manchester, AV Hill Building, Oxford Road, Manchester, M13 9PT, UK
| | - Samantha L Smith
- Centre for Genetics and Genomics Versus Arthritis, Centre for Musculoskeletal Research, Manchester Academic Health Science Centre, The University of Manchester, AV Hill Building, Oxford Road, Manchester, M13 9PT, UK
| | - Gisela Orozco
- Centre for Genetics and Genomics Versus Arthritis, Centre for Musculoskeletal Research, Manchester Academic Health Science Centre, The University of Manchester, AV Hill Building, Oxford Road, Manchester, M13 9PT, UK
| | - Anne Barton
- Centre for Genetics and Genomics Versus Arthritis, Centre for Musculoskeletal Research, Manchester Academic Health Science Centre, The University of Manchester, AV Hill Building, Oxford Road, Manchester, M13 9PT, UK
- NIHR Manchester Biomedical Research Centre, Manchester University NHS Foundation Trust, Manchester Academic Health Science Centre, Manchester, UK
| | - Steve Eyre
- Centre for Genetics and Genomics Versus Arthritis, Centre for Musculoskeletal Research, Manchester Academic Health Science Centre, The University of Manchester, AV Hill Building, Oxford Road, Manchester, M13 9PT, UK
| | - Paul Martin
- Centre for Genetics and Genomics Versus Arthritis, Centre for Musculoskeletal Research, Manchester Academic Health Science Centre, The University of Manchester, AV Hill Building, Oxford Road, Manchester, M13 9PT, UK.
- The Lydia Becker Institute of Immunology and Inflammation, Faculty of Biology, Medicine and Health, University of Manchester, Manchester, UK.
| |
Collapse
|
23
|
Mandric I, Schwarz T, Majumdar A, Hou K, Briscoe L, Perez R, Subramaniam M, Hafemeister C, Satija R, Ye CJ, Pasaniuc B, Halperin E. Optimized design of single-cell RNA sequencing experiments for cell-type-specific eQTL analysis. Nat Commun 2020; 11:5504. [PMID: 33127880 PMCID: PMC7599215 DOI: 10.1038/s41467-020-19365-w] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2019] [Accepted: 10/06/2020] [Indexed: 02/07/2023] Open
Abstract
Single-cell RNA-sequencing (scRNA-Seq) is a compelling approach to directly and simultaneously measure cellular composition and state, which can otherwise only be estimated by applying deconvolution methods to bulk RNA-Seq estimates. However, it has not yet become a widely used tool in population-scale analyses, due to its prohibitively high cost. Here we show that given the same budget, the statistical power of cell-type-specific expression quantitative trait loci (eQTL) mapping can be increased through low-coverage per-cell sequencing of more samples rather than high-coverage sequencing of fewer samples. We use simulations starting from one of the largest available real single-cell RNA-Seq data from 120 individuals to also show that multiple experimental designs with different numbers of samples, cells per sample and reads per cell could have similar statistical power, and choosing an appropriate design can yield large cost savings especially when multiplexed workflows are considered. Finally, we provide a practical approach on selecting cost-effective designs for maximizing cell-type-specific eQTL power which is available in the form of a web tool.
Collapse
Affiliation(s)
- Igor Mandric
- Department of Computer Science, University of California Los Angeles, 404 Westwood Plaza, Los Angeles, CA, 90095, USA
| | - Tommer Schwarz
- Bioinformatics Interdepartmental Program, University of California Los Angeles, 611 Charles E. Young Drive East, Los Angeles, CA, 90095, USA
| | - Arunabha Majumdar
- Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, 10833 Le Conte Ave, Los Angeles, CA, 90095, USA
| | - Kangcheng Hou
- Bioinformatics Interdepartmental Program, University of California Los Angeles, 611 Charles E. Young Drive East, Los Angeles, CA, 90095, USA
| | - Leah Briscoe
- Bioinformatics Interdepartmental Program, University of California Los Angeles, 611 Charles E. Young Drive East, Los Angeles, CA, 90095, USA
| | - Richard Perez
- Institute for Human Genetics, University of California San Francisco, 513 Parnassus Avenue, San Francisco, CA, 94143, USA
- Bakar Computational Health Sciences Institute, University of California San Francisco, 550 16th Street, San Francisco, CA, 94158, USA
- Division of Rheumatology, Department of Medicine, University of California San Francisco, 513 Parnassus Avenue, San Francisco, CA, 94143, USA
| | - Meena Subramaniam
- Institute for Human Genetics, University of California San Francisco, 513 Parnassus Avenue, San Francisco, CA, 94143, USA
- Bakar Computational Health Sciences Institute, University of California San Francisco, 550 16th Street, San Francisco, CA, 94158, USA
- Division of Rheumatology, Department of Medicine, University of California San Francisco, 513 Parnassus Avenue, San Francisco, CA, 94143, USA
- Bioinformatics Program, University of California San Francisco, 513 Parnassus Ave, San Francisco, CA, 94143, USA
| | | | - Rahul Satija
- New York Genome Center, 101 Avenue of the Americas, New York, NY, 10013, USA
- Center for Genomics and Systems Biology, New York University, 12 Waverly Place, New York, NY, 10003, USA
| | - Chun Jimmie Ye
- Institute for Human Genetics, University of California San Francisco, 513 Parnassus Avenue, San Francisco, CA, 94143, USA
- Bakar Computational Health Sciences Institute, University of California San Francisco, 550 16th Street, San Francisco, CA, 94158, USA
- Division of Rheumatology, Department of Medicine, University of California San Francisco, 513 Parnassus Avenue, San Francisco, CA, 94143, USA
- Bioinformatics Program, University of California San Francisco, 513 Parnassus Ave, San Francisco, CA, 94143, USA
| | - Bogdan Pasaniuc
- Bioinformatics Interdepartmental Program, University of California Los Angeles, 611 Charles E. Young Drive East, Los Angeles, CA, 90095, USA.
- Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, 10833 Le Conte Ave, Los Angeles, CA, 90095, USA.
- Department of Human Genetics, David Geffen School of Medicine, University of California Los Angeles, 6506 Gonda Center, Los Angeles, CA, 90095, USA.
- Department of Computational Medicine, David Geffen School of Medicine, University of California Los Angeles, Room 5303 Life Sciences, Los Angeles, CA, 90095, USA.
| | - Eran Halperin
- Department of Computer Science, University of California Los Angeles, 404 Westwood Plaza, Los Angeles, CA, 90095, USA
- Department of Human Genetics, David Geffen School of Medicine, University of California Los Angeles, 6506 Gonda Center, Los Angeles, CA, 90095, USA
- Department of Computational Medicine, David Geffen School of Medicine, University of California Los Angeles, Room 5303 Life Sciences, Los Angeles, CA, 90095, USA
- Department of Anesthesiology and Perioperative Medicine, David Geffen School of Medicine, University of California Los Angeles, 757 Westwood Plaza, Los Angeles, CA, 90095, USA
- Institute of Precision Health, University of California, Los Angeles, CA, USA
| |
Collapse
|
24
|
Gupta R, Schrooders Y, Hauser D, van Herwijnen M, Albrecht W, Ter Braak B, Brecklinghaus T, Castell JV, Elenschneider L, Escher S, Guye P, Hengstler JG, Ghallab A, Hansen T, Leist M, Maclennan R, Moritz W, Tolosa L, Tricot T, Verfaillie C, Walker P, van de Water B, Kleinjans J, Caiment F. Comparing in vitro human liver models to in vivo human liver using RNA-Seq. Arch Toxicol 2020; 95:573-589. [PMID: 33106934 PMCID: PMC7870774 DOI: 10.1007/s00204-020-02937-6] [Citation(s) in RCA: 174] [Impact Index Per Article: 43.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2020] [Accepted: 10/12/2020] [Indexed: 01/29/2023]
Abstract
The liver plays an important role in xenobiotic metabolism and represents a primary target for toxic substances. Many different in vitro cell models have been developed in the past decades. In this study, we used RNA-sequencing (RNA-Seq) to analyze the following human in vitro liver cell models in comparison to human liver tissue: cancer-derived cell lines (HepG2, HepaRG 3D), induced pluripotent stem cell-derived hepatocyte-like cells (iPSC-HLCs), cancerous human liver-derived assays (hPCLiS, human precision cut liver slices), non-cancerous human liver-derived assays (PHH, primary human hepatocytes) and 3D liver microtissues. First, using CellNet, we analyzed whether these liver in vitro cell models were indeed classified as liver, based on their baseline expression profile and gene regulatory networks (GRN). More comprehensive analyses using non-differentially expressed genes (non-DEGs) and differential transcript usage (DTU) were applied to assess the coverage for important liver pathways. Through different analyses, we noticed that 3D liver microtissues exhibited a high similarity with in vivo liver, in terms of CellNet (C/T score: 0.98), non-DEGs (10,363) and pathway coverage (highest for 19 out of 20 liver specific pathways shown) at the beginning of the incubation period (0 h) followed by a decrease during long-term incubation for 168 and 336 h. PHH also showed a high degree of similarity with human liver tissue and allowed stable conditions for a short-term cultivation period of 24 h. Using the same metrics, HepG2 cells illustrated the lowest similarity (C/T: 0.51, non-DEGs: 5623, and pathways coverage: least for 7 out of 20) with human liver tissue. The HepG2 are widely used in hepatotoxicity studies, however, due to their lower similarity, they should be used with caution. HepaRG models, iPSC-HLCs, and hPCLiS ranged clearly behind microtissues and PHH but showed higher similarity to human liver tissue than HepG2 cells. In conclusion, this study offers a resource of RNA-Seq data of several biological replicates of human liver cell models in vitro compared to human liver tissue.
Collapse
Affiliation(s)
- Rajinder Gupta
- Department of Toxicogenomics, School of Oncology and Developmental Biology (GROW), Maastricht University, Maastricht, The Netherlands
| | - Yannick Schrooders
- Department of Toxicogenomics, School of Oncology and Developmental Biology (GROW), Maastricht University, Maastricht, The Netherlands
| | - Duncan Hauser
- Department of Toxicogenomics, School of Oncology and Developmental Biology (GROW), Maastricht University, Maastricht, The Netherlands
| | - Marcel van Herwijnen
- Department of Toxicogenomics, School of Oncology and Developmental Biology (GROW), Maastricht University, Maastricht, The Netherlands
| | - Wiebke Albrecht
- Leibniz Research Centre for Working Environment and Human Factors, Technical University of Dortmund (IfADo), Dortmund, Germany
| | - Bas Ter Braak
- Division of Drug Discovery and Safety, Leiden Academic Centre for Drug Research, Leiden University, PO Box 9503, 2300 RA, Leiden, The Netherlands
| | - Tim Brecklinghaus
- Leibniz Research Centre for Working Environment and Human Factors, Technical University of Dortmund (IfADo), Dortmund, Germany
| | - Jose V Castell
- Instituto de Investigación Sanitaria La Fe, Experimental Hepatology Unit, Valencia, Spain
| | - Leroy Elenschneider
- Fraunhofer Institute for Toxicology and Experimental Medicine Preclinical Pharmacology and In-Vitro Toxicology, Nikolai-Fuchs-Straße 1, 30625, Hannover, Germany
| | - Sylvia Escher
- Fraunhofer Institute for Toxicology and Experimental Medicine Preclinical Pharmacology and In-Vitro Toxicology, Nikolai-Fuchs-Straße 1, 30625, Hannover, Germany
| | - Patrick Guye
- InSphero AG, Wagistrasse 27, 8952, Schlieren, Switzerland
| | - Jan G Hengstler
- Leibniz Research Centre for Working Environment and Human Factors, Technical University of Dortmund (IfADo), Dortmund, Germany
| | - Ahmed Ghallab
- Leibniz Research Centre for Working Environment and Human Factors, Technical University of Dortmund (IfADo), Dortmund, Germany
- Department of Forensic Medicine and Toxicology, Faculty of Veterinary Medicine, South Valley University, Qena, 83523, Egypt
| | - Tanja Hansen
- Fraunhofer Institute for Toxicology and Experimental Medicine Preclinical Pharmacology and In-Vitro Toxicology, Nikolai-Fuchs-Straße 1, 30625, Hannover, Germany
| | - Marcel Leist
- In Vitro Toxicology and Biomedicine, Department Inaugurated, Doerenkamp-Zbinden Foundation, University of Konstanz, Konstanz, Germany
| | - Richard Maclennan
- Cyprotex Discovery, No 24 Mereside, Alderley Park, Cheshire, SK10 4TG, UK
| | | | - Laia Tolosa
- Instituto de Investigación Sanitaria La Fe, Unidad Hepatología Experimental, Valencia, Spain
| | - Tine Tricot
- Stem Cell Institute, Department of Development and Regeneration, KU Leuven, Herestraat 49, 3000, Leuven, Belgium
| | - Catherine Verfaillie
- Stem Cell Institute, Department of Development and Regeneration, KU Leuven, Herestraat 49, 3000, Leuven, Belgium
| | - Paul Walker
- Cyprotex Discovery, No 24 Mereside, Alderley Park, Cheshire, SK10 4TG, UK
| | - Bob van de Water
- Division of Drug Discovery and Safety, Leiden Academic Centre for Drug Research, Leiden University, PO Box 9503, 2300 RA, Leiden, The Netherlands
| | - Jos Kleinjans
- Department of Toxicogenomics, School of Oncology and Developmental Biology (GROW), Maastricht University, Maastricht, The Netherlands
| | - Florian Caiment
- Department of Toxicogenomics, School of Oncology and Developmental Biology (GROW), Maastricht University, Maastricht, The Netherlands.
| |
Collapse
|
25
|
Dusart P, Hallström BM, Renné T, Odeberg J, Uhlén M, Butler LM. A Systems-Based Map of Human Brain Cell-Type Enriched Genes and Malignancy-Associated Endothelial Changes. Cell Rep 2020; 29:1690-1706.e4. [PMID: 31693905 DOI: 10.1016/j.celrep.2019.09.088] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2019] [Revised: 04/09/2019] [Accepted: 09/27/2019] [Indexed: 02/02/2023] Open
Abstract
Changes in the endothelium of the cerebral vasculature can contribute to inflammatory, thrombotic, and malignant disorders. The importance of defining cell-type-specific genes and their modification in disease is increasingly recognized. Here, we develop a bioinformatics-based approach to identify normal brain cell-enriched genes, using bulk RNA sequencing (RNA-seq) data from 238 normal human cortex samples from 2 independent cohorts. We compare endothelial cell-enriched gene profiles with astrocyte, oligodendrocyte, neuron, and microglial cell profiles. Endothelial changes in malignant disease are explored using RNA-seq data from 516 lower-grade gliomas and 401 glioblastomas. Lower-grade gliomas appear to be an "endothelial intermediate" between normal brain and glioblastoma. We apply our method for the prediction of glioblastoma-specific endothelial biomarkers, providing potential diagnostic or therapeutic targets. In summary, we provide a roadmap of endothelial cell identity in normal and malignant brain, using a method developed to resolve bulk RNA-seq into constituent cell-type-enriched profiles.
Collapse
Affiliation(s)
- Philip Dusart
- Science for Life Laboratory, Department of Protein Science, Royal Institute of Technology (KTH), 171 21 Stockholm, Sweden; K.G. Jebsen Thrombosis Research and Expertise Centre, Department of Clinical Medicine, The Arctic University of Norway, 9019 Tromsø, Norway
| | - Björn Mikael Hallström
- Science for Life Laboratory, Department of Protein Science, Royal Institute of Technology (KTH), 171 21 Stockholm, Sweden
| | - Thomas Renné
- Institute for Clinical Chemistry and Laboratory Medicine, University Medical Centre Hamburg-Eppendorf, 20246 Hamburg, Germany
| | - Jacob Odeberg
- Science for Life Laboratory, Department of Protein Science, Royal Institute of Technology (KTH), 171 21 Stockholm, Sweden; K.G. Jebsen Thrombosis Research and Expertise Centre, Department of Clinical Medicine, The Arctic University of Norway, 9019 Tromsø, Norway; The University Hospital of North Norway (UNN), PB100, 9038 Tromsø, Norway; Department of Hematology, Karolinska University Hospital, 171 77 Stockholm, Sweden
| | - Mathias Uhlén
- Science for Life Laboratory, Department of Protein Science, Royal Institute of Technology (KTH), 171 21 Stockholm, Sweden
| | - Lynn Marie Butler
- Science for Life Laboratory, Department of Protein Science, Royal Institute of Technology (KTH), 171 21 Stockholm, Sweden; K.G. Jebsen Thrombosis Research and Expertise Centre, Department of Clinical Medicine, The Arctic University of Norway, 9019 Tromsø, Norway; Institute for Clinical Chemistry and Laboratory Medicine, University Medical Centre Hamburg-Eppendorf, 20246 Hamburg, Germany; Clinical Chemistry and Blood Coagulation Research, Department of Molecular Medicine and Surgery, Karolinska Institute, 171 76 Stockholm, Sweden.
| |
Collapse
|
26
|
Mukherjee S, Agarwal D, Zhang NR, Bhattacharya BB. Distribution-Free Multisample Tests Based on Optimal Matchings With Applications to Single Cell Genomics. J Am Stat Assoc 2020. [DOI: 10.1080/01621459.2020.1791131] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Somabha Mukherjee
- Department of Statistics, The Wharton School, University of Pennsylvania, Philadelphia, PA
| | - Divyansh Agarwal
- Department of Statistics, The Wharton School, University of Pennsylvania, Philadelphia, PA
| | - Nancy R. Zhang
- Department of Statistics, The Wharton School, University of Pennsylvania, Philadelphia, PA
| | | |
Collapse
|
27
|
Noé A, Cargill TN, Nielsen CM, Russell AJC, Barnes E. The Application of Single-Cell RNA Sequencing in Vaccinology. J Immunol Res 2020; 2020:8624963. [PMID: 32802896 PMCID: PMC7411487 DOI: 10.1155/2020/8624963] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2020] [Accepted: 07/09/2020] [Indexed: 02/06/2023] Open
Abstract
Single-cell RNA sequencing allows highly detailed profiling of cellular immune responses from limited-volume samples, advancing prospects of a new era of systems immunology. The power of single-cell RNA sequencing offers various opportunities to decipher the immune response to infectious diseases and vaccines. Here, we describe the potential uses of single-cell RNA sequencing methods in prophylactic vaccine development, concentrating on infectious diseases including COVID-19. Using examples from several diseases, we review how single-cell RNA sequencing has been used to evaluate the immunological response to different vaccine platforms and regimens. By highlighting published and unpublished single-cell RNA sequencing studies relevant to vaccinology, we discuss some general considerations how the field could be enriched with the widespread adoption of this technology.
Collapse
MESH Headings
- Animals
- Betacoronavirus/immunology
- COVID-19
- Cell Line
- Clinical Trials as Topic
- Coronavirus Infections/epidemiology
- Coronavirus Infections/immunology
- Coronavirus Infections/prevention & control
- Coronavirus Infections/virology
- Disease Models, Animal
- Drug Evaluation, Preclinical
- Host-Pathogen Interactions/genetics
- Host-Pathogen Interactions/immunology
- Humans
- Immunity, Cellular/genetics
- Immunity, Innate/genetics
- Immunogenicity, Vaccine
- Pandemics/prevention & control
- Pneumonia, Viral/epidemiology
- Pneumonia, Viral/immunology
- Pneumonia, Viral/prevention & control
- Pneumonia, Viral/virology
- RNA, Viral/isolation & purification
- RNA-Seq/methods
- SARS-CoV-2
- Single-Cell Analysis
- Vaccinology/methods
- Viral Vaccines/administration & dosage
- Viral Vaccines/immunology
Collapse
Affiliation(s)
- Andrés Noé
- The Jenner Institute, University of Oxford, Old Road Campus Research Building, Oxford OX3 7DQ, UK
| | - Tamsin N. Cargill
- Peter Medawar Building for Pathogen Research and Oxford NIHR Biomedical Research Centre, Nuffield Department of Medicine, University of Oxford, South Parks Road, Oxford OX1 3SY, UK
- Translational Gastroenterology Unit, John Radcliffe Hospital, Oxford OX3 9DU, UK
| | - Carolyn M. Nielsen
- The Jenner Institute, University of Oxford, Old Road Campus Research Building, Oxford OX3 7DQ, UK
| | | | - Eleanor Barnes
- Peter Medawar Building for Pathogen Research and Oxford NIHR Biomedical Research Centre, Nuffield Department of Medicine, University of Oxford, South Parks Road, Oxford OX1 3SY, UK
- Translational Gastroenterology Unit, John Radcliffe Hospital, Oxford OX3 9DU, UK
| |
Collapse
|
28
|
Ahmed R, Omidian Z, Giwa A, Cornwell B, Majety N, Bell DR, Lee S, Zhang H, Michels A, Desiderio S, Sadegh-Nasseri S, Rabb H, Gritsch S, Suva ML, Cahan P, Zhou R, Jie C, Donner T, Hamad ARA. A Public BCR Present in a Unique Dual-Receptor-Expressing Lymphocyte from Type 1 Diabetes Patients Encodes a Potent T Cell Autoantigen. Cell 2020; 177:1583-1599.e16. [PMID: 31150624 DOI: 10.1016/j.cell.2019.05.007] [Citation(s) in RCA: 79] [Impact Index Per Article: 19.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2018] [Revised: 12/10/2018] [Accepted: 05/02/2019] [Indexed: 12/17/2022]
Abstract
T and B cells are the two known lineages of adaptive immune cells. Here, we describe a previously unknown lymphocyte that is a dual expresser (DE) of TCR and BCR and key lineage markers of both B and T cells. In type 1 diabetes (T1D), DEs are predominated by one clonotype that encodes a potent CD4 T cell autoantigen in its antigen binding site. Molecular dynamics simulations revealed that this peptide has an optimal binding register for diabetogenic HLA-DQ8. In concordance, a synthetic version of the peptide forms stable DQ8 complexes and potently stimulates autoreactive CD4 T cells from T1D patients, but not healthy controls. Moreover, mAbs bearing this clonotype are autoreactive against CD4 T cells and inhibit insulin tetramer binding to CD4 T cells. Thus, compartmentalization of adaptive immune cells into T and B cells is not absolute, and violators of this paradigm are likely key drivers of autoimmune diseases.
Collapse
Affiliation(s)
- Rizwan Ahmed
- Department of Pathology, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Zahra Omidian
- Department of Pathology, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Adebola Giwa
- Department of Pediatrics, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Benjamin Cornwell
- Department of Pathology, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Neha Majety
- Department of Pathology, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - David R Bell
- Computational Biology Center, IBM Thomas J. Watson Research Center, Yorktown Heights, NY 10598, USA
| | - Sangyun Lee
- Computational Biology Center, IBM Thomas J. Watson Research Center, Yorktown Heights, NY 10598, USA
| | - Hao Zhang
- Department of Molecular Microbiology and Immunology, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD 21205, USA
| | - Aaron Michels
- Barbara Davis Center for Diabetes, University of Colorado, Aurora, CO 80045, USA
| | - Stephen Desiderio
- Department of Molecular Biology and Genetics and Institute for Cell Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | | | - Hamid Rabb
- Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Simon Gritsch
- Department of Pathology and Center for Cancer Research, Massachusetts General Hospital and Harvard Medical School, Boston, MA, 02114, USA
| | - Mario L Suva
- Department of Pathology and Center for Cancer Research, Massachusetts General Hospital and Harvard Medical School, Boston, MA, 02114, USA; Broad Institute of Harvard and MIT, Cambridge, MA 02142, USA
| | - Patrick Cahan
- Department of Molecular Biology and Genetics and Institute for Cell Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA; Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Ruhong Zhou
- Computational Biology Center, IBM Thomas J. Watson Research Center, Yorktown Heights, NY 10598, USA; Department of Chemistry, Columbia University, New York, NY 10027, USA.
| | - Chunfa Jie
- Department of Biochemistry and Nutrition, Des Moines University, Des Moines, IA 50312, USA
| | - Thomas Donner
- Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Abdel Rahim A Hamad
- Department of Pathology, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA; Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA.
| |
Collapse
|
29
|
Zhang MJ, Ntranos V, Tse D. Determining sequencing depth in a single-cell RNA-seq experiment. Nat Commun 2020; 11:774. [PMID: 32034137 PMCID: PMC7005864 DOI: 10.1038/s41467-020-14482-y] [Citation(s) in RCA: 52] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2018] [Accepted: 12/13/2019] [Indexed: 11/25/2022] Open
Abstract
An underlying question for virtually all single-cell RNA sequencing experiments is how to allocate the limited sequencing budget: deep sequencing of a few cells or shallow sequencing of many cells? Here we present a mathematical framework which reveals that, for estimating many important gene properties, the optimal allocation is to sequence at a depth of around one read per cell per gene. Interestingly, the corresponding optimal estimator is not the widely-used plug-in estimator, but one developed via empirical Bayes. For single-cell RNA-seq experiments the sequencing budget is limited, and how it should be optimally allocated to maximize information is not clear. Here the authors develop a mathematical framework to show that, for estimating many gene properties, the optimal allocation is to sequence at the depth of one read per cell per gene.
Collapse
Affiliation(s)
- Martin Jinye Zhang
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Vasilis Ntranos
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA.,Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, CA, USA
| | - David Tse
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA.
| |
Collapse
|
30
|
Lähnemann D, Köster J, Szczurek E, McCarthy DJ, Hicks SC, Robinson MD, Vallejos CA, Campbell KR, Beerenwinkel N, Mahfouz A, Pinello L, Skums P, Stamatakis A, Attolini CSO, Aparicio S, Baaijens J, Balvert M, Barbanson BD, Cappuccio A, Corleone G, Dutilh BE, Florescu M, Guryev V, Holmer R, Jahn K, Lobo TJ, Keizer EM, Khatri I, Kielbasa SM, Korbel JO, Kozlov AM, Kuo TH, Lelieveldt BP, Mandoiu II, Marioni JC, Marschall T, Mölder F, Niknejad A, Rączkowska A, Reinders M, Ridder JD, Saliba AE, Somarakis A, Stegle O, Theis FJ, Yang H, Zelikovsky A, McHardy AC, Raphael BJ, Shah SP, Schönhuth A. Eleven grand challenges in single-cell data science. Genome Biol 2020; 21:31. [PMID: 32033589 PMCID: PMC7007675 DOI: 10.1186/s13059-020-1926-6] [Citation(s) in RCA: 565] [Impact Index Per Article: 141.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2019] [Accepted: 01/02/2020] [Indexed: 02/08/2023] Open
Abstract
The recent boom in microfluidics and combinatorial indexing strategies, combined with low sequencing costs, has empowered single-cell sequencing technology. Thousands-or even millions-of cells analyzed in a single experiment amount to a data revolution in single-cell biology and pose unique data science problems. Here, we outline eleven challenges that will be central to bringing this emerging field of single-cell data science forward. For each challenge, we highlight motivating research questions, review prior work, and formulate open problems. This compendium is for established researchers, newcomers, and students alike, highlighting interesting and rewarding problems for the coming years.
Collapse
Affiliation(s)
- David Lähnemann
- Algorithms for Reproducible Bioinformatics, Genome Informatics, Institute of Human Genetics, University Hospital Essen, University of Duisburg-Essen, Essen, Germany
- Department of Paediatric Oncology, Haematology and Immunology, Medical Faculty, Heinrich Heine University, University Hospital, Düsseldorf, Germany
- Computational Biology of Infection Research Group, Helmholtz Centre for Infection Research, Braunschweig, Germany
| | - Johannes Köster
- Algorithms for Reproducible Bioinformatics, Genome Informatics, Institute of Human Genetics, University Hospital Essen, University of Duisburg-Essen, Essen, Germany
- Medical Oncology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, USA
| | - Ewa Szczurek
- Institute of Informatics, Faculty of Mathematics, Informatics and Mechanics, University of Warsaw, Warszawa, Poland
| | - Davis J. McCarthy
- Bioinformatics and Cellular Genomics, St Vincent’s Institute of Medical Research, Fitzroy, Australia
- Melbourne Integrative Genomics, School of BioSciences–School of Mathematics & Statistics, Faculty of Science, University of Melbourne, Melbourne, Australia
| | - Stephanie C. Hicks
- Department of Biostatistics, Johns Hopkins University, Baltimore, MD USA
| | - Mark D. Robinson
- Institute of Molecular Life Sciences and SIB Swiss Institute of Bioinformatics, University of Zürich, Zürich, Switzerland
| | - Catalina A. Vallejos
- MRC Human Genetics Unit, Institute of Genetics and Molecular Medicine, University of Edinburgh, Western General Hospital, Edinburgh, UK
- The Alan Turing Institute, British Library, London, UK
| | - Kieran R. Campbell
- Department of Statistics, University of British Columbia, Vancouver, Canada
- Department of Molecular Oncology, BC Cancer Agency, Vancouver, Canada
- Data Science Institute, University of British Columbia, Vancouver, Canada
| | - Niko Beerenwinkel
- Department of Biosystems Science and Engineering, ETH Zurich, Basel, Switzerland
- SIB Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Ahmed Mahfouz
- Leiden Computational Biology Center, Leiden University Medical Center, Leiden, The Netherlands
- Delft Bioinformatics Lab, Faculty of Electrical Engineering, Mathematics and Computer Science, Delft University of Technology, Delft, The Netherlands
| | - Luca Pinello
- Molecular Pathology Unit and Center for Cancer Research, Massachusetts General Hospital Research Institute, Charlestown, USA
- Department of Pathology, Harvard Medical School, Boston, USA
- Broad Institute of Harvard and MIT, Cambridge, MA USA
| | - Pavel Skums
- Department of Computer Science, Georgia State University, Atlanta, USA
| | - Alexandros Stamatakis
- Computational Molecular Evolution Group, Heidelberg Institute for Theoretical Studies, Heidelberg, Germany
- Institute for Theoretical Informatics, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | | | - Samuel Aparicio
- Department of Molecular Oncology, BC Cancer Agency, Vancouver, Canada
- Department of Pathology and Laboratory Medicine, University of British Columbia, Vancouver, Canada
| | - Jasmijn Baaijens
- Life Sciences and Health, Centrum Wiskunde & Informatica, Amsterdam, The Netherlands
| | - Marleen Balvert
- Life Sciences and Health, Centrum Wiskunde & Informatica, Amsterdam, The Netherlands
- Theoretical Biology and Bioinformatics, Science for Life, Utrecht University, Utrecht, The Netherlands
| | - Buys de Barbanson
- Center for Molecular Medicine, University Medical Center Utrecht, Utrecht, The Netherlands
- Oncode Institute, Utrecht, The Netherlands
- Quantitative biology, Hubrecht Institute, Utrecht, The Netherlands
| | - Antonio Cappuccio
- Institute for Advanced Study, University of Amsterdam, Amsterdam, The Netherlands
| | - Giacomo Corleone
- Department of Surgery and Cancer, The Imperial Centre for Translational and Experimental Medicine, Imperial College London, London, UK
| | - Bas E. Dutilh
- Theoretical Biology and Bioinformatics, Science for Life, Utrecht University, Utrecht, The Netherlands
- Centre for Molecular and Biomolecular Informatics, Radboud University Medical Center, Nijmegen, The Netherlands
| | - Maria Florescu
- Center for Molecular Medicine, University Medical Center Utrecht, Utrecht, The Netherlands
- Oncode Institute, Utrecht, The Netherlands
- Quantitative biology, Hubrecht Institute, Utrecht, The Netherlands
| | - Victor Guryev
- European Research Institute for the Biology of Ageing, University Medical Center Groningen, University of Groningen, Groningen, The Netherlands
| | - Rens Holmer
- Bioinformatics Group, Wageningen University, Wageningen, The Netherlands
| | - Katharina Jahn
- Department of Biosystems Science and Engineering, ETH Zurich, Basel, Switzerland
- SIB Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Thamar Jessurun Lobo
- European Research Institute for the Biology of Ageing, University Medical Center Groningen, University of Groningen, Groningen, The Netherlands
| | - Emma M. Keizer
- Biometris, Wageningen University & Research, Wageningen, The Netherlands
| | - Indu Khatri
- Department of Immunohematology and Blood Transfusion, Leiden University Medical Center, Leiden, The Netherlands
| | - Szymon M. Kielbasa
- Department of Biomedical Data Sciences, Leiden University Medical Center, Leiden, The Netherlands
| | - Jan O. Korbel
- Genome Biology Unit, European Molecular Biology Laboratory, Heidelberg, Germany
| | - Alexey M. Kozlov
- Computational Molecular Evolution Group, Heidelberg Institute for Theoretical Studies, Heidelberg, Germany
| | - Tzu-Hao Kuo
- Computational Biology of Infection Research Group, Helmholtz Centre for Infection Research, Braunschweig, Germany
| | - Boudewijn P.F. Lelieveldt
- PRB lab, Delft University of Technology, Delft, The Netherlands
- Division of Image Processing, Department of Radiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Ion I. Mandoiu
- Computer Science & Engineering Department, University of Connecticut, Storrs, USA
| | - John C. Marioni
- Cancer Research UK Cambridge Institute, Li Ka Shing Centre, University of Cambridge, Cambridge, UK
- Wellcome Trust Sanger Institute, Wellcome Genome Campus, Hinxton, UK
- European Molecular Biology Laboratory, European Bioinformatics Institute, Hinxton, UK
| | - Tobias Marschall
- Center for Bioinformatics, Saarland University, Saarbrücken, Germany
- Max Planck Institute for Informatics, Saarbrücken, Germany
| | - Felix Mölder
- Algorithms for Reproducible Bioinformatics, Genome Informatics, Institute of Human Genetics, University Hospital Essen, University of Duisburg-Essen, Essen, Germany
- Institute of Pathology, University Hospital Essen, University of Duisburg-Essen, Essen, Germany
| | - Amir Niknejad
- Computation molecular design, Zuse Institute Berlin, Berlin, Germany
- Mathematics Department, Mount Saint Vincent, New York, USA
| | - Alicja Rączkowska
- Institute of Informatics, Faculty of Mathematics, Informatics and Mechanics, University of Warsaw, Warszawa, Poland
| | - Marcel Reinders
- Leiden Computational Biology Center, Leiden University Medical Center, Leiden, The Netherlands
- Delft Bioinformatics Lab, Faculty of Electrical Engineering, Mathematics and Computer Science, Delft University of Technology, Delft, The Netherlands
| | - Jeroen de Ridder
- Center for Molecular Medicine, University Medical Center Utrecht, Utrecht, The Netherlands
- Oncode Institute, Utrecht, The Netherlands
| | - Antoine-Emmanuel Saliba
- Helmholtz Institute for RNA-based Infection Research, Helmholtz-Center for Infection Research, Würzburg, Germany
| | - Antonios Somarakis
- Division of Image Processing, Department of Radiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Oliver Stegle
- Genome Biology Unit, European Molecular Biology Laboratory, Heidelberg, Germany
- European Molecular Biology Laboratory, European Bioinformatics Institute, Hinxton, UK
- Division of Computational Genomics and Systems Genetics, German Cancer Research Center–DKFZ, Heidelberg, Germany
| | - Fabian J. Theis
- Institute of Computational Biology, Helmholtz Zentrum München–German Research Center for Environmental Health, Neuherberg, Germany
| | - Huan Yang
- Division of Drug Discovery and Safety, Leiden Academic Center for Drug Research–LACDR–Leiden University, Leiden, The Netherlands
| | - Alex Zelikovsky
- Department of Computer Science, Georgia State University, Atlanta, USA
- The Laboratory of Bioinformatics, I.M. Sechenov First Moscow State Medical University, Moscow, Russia
| | - Alice C. McHardy
- Computational Biology of Infection Research Group, Helmholtz Centre for Infection Research, Braunschweig, Germany
| | | | - Sohrab P. Shah
- Computational Oncology, Department of Epidemiology and Biostatistics, Memorial Sloan Kettering Cancer Center, New York, USA
| | - Alexander Schönhuth
- Life Sciences and Health, Centrum Wiskunde & Informatica, Amsterdam, The Netherlands
- Theoretical Biology and Bioinformatics, Science for Life, Utrecht University, Utrecht, The Netherlands
| |
Collapse
|
31
|
Molecular profiling of single neurons of known identity in two ganglia from the crab Cancer borealis. Proc Natl Acad Sci U S A 2019; 116:26980-26990. [PMID: 31806754 PMCID: PMC6936480 DOI: 10.1073/pnas.1911413116] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Single-cell transcriptional profiling has become a widespread tool in cell identification, particularly in the nervous system, based on the notion that genomic information determines cell identity. However, many cell-type classification studies are unconstrained by other cellular attributes (e.g., morphology, physiology). Here, we systematically test how accurately transcriptional profiling can assign cell identity to well-studied anatomically and functionally identified neurons in 2 small neuronal networks. While these neurons clearly possess distinct patterns of gene expression across cell types, their expression profiles are not sufficient to unambiguously confirm their identity. We suggest that true cell identity can only be determined by combining gene expression data with other cellular attributes such as innervation pattern, morphology, or physiology. Understanding circuit organization depends on identification of cell types. Recent advances in transcriptional profiling methods have enabled classification of cell types by their gene expression. While exceptionally powerful and high throughput, the ground-truth validation of these methods is difficult: If cell type is unknown, how does one assess whether a given analysis accurately captures neuronal identity? To shed light on the capabilities and limitations of solely using transcriptional profiling for cell-type classification, we performed 2 forms of transcriptional profiling—RNA-seq and quantitative RT-PCR, in single, unambiguously identified neurons from 2 small crustacean neuronal networks: The stomatogastric and cardiac ganglia. We then combined our knowledge of cell type with unbiased clustering analyses and supervised machine learning to determine how accurately functionally defined neuron types can be classified by expression profile alone. The results demonstrate that expression profile is able to capture neuronal identity most accurately when combined with multimodal information that allows for post hoc grouping, so analysis can proceed from a supervised perspective. Solely unsupervised clustering can lead to misidentification and an inability to distinguish between 2 or more cell types. Therefore, this study supports the general utility of cell identification by transcriptional profiling, but adds a caution: It is difficult or impossible to know under what conditions transcriptional profiling alone is capable of assigning cell identity. Only by combining multiple modalities of information such as physiology, morphology, or innervation target can neuronal identity be unambiguously determined.
Collapse
|
32
|
Jiang Y, Jiang Y, Wang S, Zhang Q, Ding X. Optimal sequencing depth design for whole genome re-sequencing in pigs. BMC Bioinformatics 2019; 20:556. [PMID: 31703550 PMCID: PMC6839175 DOI: 10.1186/s12859-019-3164-z] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2019] [Accepted: 10/16/2019] [Indexed: 12/30/2022] Open
Abstract
BACKGROUND As whole-genome sequencing is becoming a routine technique, it is important to identify a cost-effective depth of sequencing for such studies. However, the relationship between sequencing depth and biological results from the aspects of whole-genome coverage, variant discovery power and the quality of variants is unclear, especially in pigs. We sequenced the genomes of three Yorkshire boars at an approximately 20X depth on the Illumina HiSeq X Ten platform and downloaded whole-genome sequencing data for three Duroc and three Landrace pigs with an approximately 20X depth for each individual. Then, we downsampled the deep genome data by extracting twelve different proportions of 0.05, 0.1, 0.15, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8 and 0.9 paired reads from the original bam files to mimic the sequence data of the same individuals at sequencing depths of 1.09X, 2.18X, 3.26X, 4.35X, 6.53X, 8.70X, 10.88X, 13.05X, 15.22X, 17.40X, 19.57X and 21.75X to evaluate the influence of genome coverage, the variant discovery rate and genotyping accuracy as a function of sequencing depth. In addition, SNP chip data for Yorkshire pigs were used as a validation for the comparison of single-sample calling and multisample calling algorithms. RESULTS Our results indicated that 10X is an ideal practical depth for achieving plateau coverage and discovering accurate variants, which achieved greater than 99% genome coverage. The number of false-positive variants was increased dramatically at a depth of less than 4X, which covered 95% of the whole genome. In addition, the comparison of multi- and single-sample calling showed that multisample calling was more sensitive than single-sample calling, especially at lower depths. The number of variants discovered under multisample calling was 13-fold and 2-fold higher than that under single-sample calling at 1X and 22X, respectively. A large difference was observed when the depth was less than 4.38X. However, more false-positive variants were detected under multisample calling. CONCLUSIONS Our research will inform important study design decisions regarding whole-genome sequencing depth. Our results will be helpful for choosing the appropriate depth to achieve the same power for studies performed under limited budgets.
Collapse
Affiliation(s)
- Yifan Jiang
- National Engineering Laboratory for Animal Breeding, Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture, College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| | - Yao Jiang
- National Engineering Laboratory for Animal Breeding, Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture, College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| | - Sheng Wang
- National Engineering Laboratory for Animal Breeding, Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture, College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| | - Qin Zhang
- Shandong Provincial Key Laboratory of Animal Biotechnology and Disease Control and Prevention, College of Animal Science and Technology, Shandong Agricultural University, Taian, 271001 China
| | - Xiangdong Ding
- National Engineering Laboratory for Animal Breeding, Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture, College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| |
Collapse
|
33
|
Zeng W, Chen X, Duren Z, Wang Y, Jiang R, Wong WH. DC3 is a method for deconvolution and coupled clustering from bulk and single-cell genomics data. Nat Commun 2019; 10:4613. [PMID: 31601804 PMCID: PMC6787340 DOI: 10.1038/s41467-019-12547-1] [Citation(s) in RCA: 40] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2019] [Accepted: 09/10/2019] [Indexed: 12/11/2022] Open
Abstract
Characterizing and interpreting heterogeneous mixtures at the cellular level is a critical problem in genomics. Single-cell assays offer an opportunity to resolve cellular level heterogeneity, e.g., scRNA-seq enables single-cell expression profiling, and scATAC-seq identifies active regulatory elements. Furthermore, while scHi-C can measure the chromatin contacts (i.e., loops) between active regulatory elements to target genes in single cells, bulk HiChIP can measure such contacts in a higher resolution. In this work, we introduce DC3 (De-Convolution and Coupled-Clustering) as a method for the joint analysis of various bulk and single-cell data such as HiChIP, RNA-seq and ATAC-seq from the same heterogeneous cell population. DC3 can simultaneously identify distinct subpopulations, assign single cells to the subpopulations (i.e., clustering) and de-convolve the bulk data into subpopulation-specific data. The subpopulation-specific profiles of gene expression, chromatin accessibility and enhancer-promoter contact obtained by DC3 provide a comprehensive characterization of the gene regulatory system in each subpopulation.
Collapse
Affiliation(s)
- Wanwen Zeng
- Department of Statistics, Department of Biomedical Data Science, Bio-X Program, Stanford University, Stanford, CA, 94305, USA
- MOE Key Laboratory of Bioinformatics, Bioinformatics Division, Department of Automation, Beijing National Research Center for Information Science and Technology, Tsinghua University, Beijing, 100084, China
| | - Xi Chen
- Department of Statistics, Department of Biomedical Data Science, Bio-X Program, Stanford University, Stanford, CA, 94305, USA
| | - Zhana Duren
- Department of Statistics, Department of Biomedical Data Science, Bio-X Program, Stanford University, Stanford, CA, 94305, USA
| | - Yong Wang
- CEMS, NCMIS, MDIS, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing, 100080, China
- Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, Kunming, 650223, China
| | - Rui Jiang
- MOE Key Laboratory of Bioinformatics, Bioinformatics Division, Department of Automation, Beijing National Research Center for Information Science and Technology, Tsinghua University, Beijing, 100084, China.
| | - Wing Hung Wong
- Department of Statistics, Department of Biomedical Data Science, Bio-X Program, Stanford University, Stanford, CA, 94305, USA.
| |
Collapse
|
34
|
Stevenson K, Uversky VN. Single-cell RNA-Seq: a next generation sequencing tool for a high-resolution view of the individual cell. J Biomol Struct Dyn 2019; 38:3730-3735. [PMID: 31524088 DOI: 10.1080/07391102.2019.1659859] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Abstract
Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Kevin Stevenson
- Department of Molecular Medicine, Morsani College of Medicine, University of South Florida, Tampa, FL, USA
| | - Vladimir N Uversky
- Department of Molecular Medicine, Morsani College of Medicine, University of South Florida, Tampa, FL, USA.,USF Health Byrd Alzheimer's Research Institute, Morsani College of Medicine, University of South Florida, Tampa, FL, USA.,Laboratory of New Methods in Biology, Institute for Biological Instrumentation, Russian Academy of Sciences, Moscow Region, Russia
| |
Collapse
|
35
|
Saldaña-Meyer R, Rodriguez-Hernaez J, Escobar T, Nishana M, Jácome-López K, Nora EP, Bruneau BG, Tsirigos A, Furlan-Magaril M, Skok J, Reinberg D. RNA Interactions Are Essential for CTCF-Mediated Genome Organization. Mol Cell 2019; 76:412-422.e5. [PMID: 31522988 DOI: 10.1016/j.molcel.2019.08.015] [Citation(s) in RCA: 169] [Impact Index Per Article: 33.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2019] [Revised: 07/15/2019] [Accepted: 08/16/2019] [Indexed: 01/14/2023]
Abstract
The function of the CCCTC-binding factor (CTCF) in the organization of the genome has become an important area of investigation, but the mechanisms by which CTCF dynamically contributes to genome organization are not clear. We previously discovered that CTCF binds to large numbers of endogenous RNAs, promoting its self-association. In this regard, we now report two independent features that disrupt CTCF association with chromatin: inhibition of transcription and disruption of CTCF-RNA interactions through mutations of 2 of its 11 zinc fingers that are not required for CTCF binding to its cognate DNA site: zinc finger 1 (ZF1) or zinc finger 10 (ZF10). These mutations alter gene expression profiles as CTCF mutants lose their ability to form chromatin loops and thus the ability to insulate chromatin domains and to mediate CTCF long-range genomic interactions. Our results point to the importance of CTCF-mediated RNA interactions as a structural component of genome organization.
Collapse
Affiliation(s)
- Ricardo Saldaña-Meyer
- Department of Biochemistry and Molecular Pharmacology, NYULSoM, New York, NY, USA; Howard Hughes Medical Institute, NYU Langone Health, New York, NY, USA.
| | | | - Thelma Escobar
- Department of Biochemistry and Molecular Pharmacology, NYULSoM, New York, NY, USA; Howard Hughes Medical Institute, NYU Langone Health, New York, NY, USA
| | | | - Karina Jácome-López
- Departamento de Genética Molecular, Instituto de Fisiología Celular, Universidad Nacional Autónoma de México, Mexico City, Mexico
| | - Elphege P Nora
- Gladstone Institutes, San Francisco, CA 94158, USA; Roddenberry Center for Stem Cell Biology and Medicine at Gladstone, San Francisco, CA 94158, USA
| | - Benoit G Bruneau
- Gladstone Institutes, San Francisco, CA 94158, USA; Roddenberry Center for Stem Cell Biology and Medicine at Gladstone, San Francisco, CA 94158, USA; Department of Pediatrics, University of California, San Francisco, San Francisco, CA 94143, USA; Cardiovascular Research Institute, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Aristotelis Tsirigos
- Department of Pathology, NYULSoM, New York, NY, USA; Applied Bioinformatics Laboratories, NYU School of Medicine, New York, NY 10016, USA
| | - Mayra Furlan-Magaril
- Departamento de Genética Molecular, Instituto de Fisiología Celular, Universidad Nacional Autónoma de México, Mexico City, Mexico
| | - Jane Skok
- Department of Pathology, NYULSoM, New York, NY, USA
| | - Danny Reinberg
- Department of Biochemistry and Molecular Pharmacology, NYULSoM, New York, NY, USA; Howard Hughes Medical Institute, NYU Langone Health, New York, NY, USA.
| |
Collapse
|
36
|
Characterizing Human Cell Types and Tissue Origin Using the Benford Law. Cells 2019; 8:cells8091004. [PMID: 31470662 PMCID: PMC6770594 DOI: 10.3390/cells8091004] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2019] [Revised: 08/27/2019] [Accepted: 08/28/2019] [Indexed: 11/16/2022] Open
Abstract
Processing massive transcriptomic datasets in a meaningful manner requires novel, possibly interdisciplinary, approaches. One principle that can address this challenge is the Benford law (BL), which posits that the occurrence probability of a leading digit in a large numerical dataset decreases as its value increases. Here, we analyzed large single-cell and bulk RNA-seq datasets to test whether cell types and tissue origins can be differentiated based on the adherence of specific genes to the BL. Then, we used the Benford adherence scores of these genes as inputs to machine-learning algorithms and tested their separation accuracy. We found that genes selected based on their first-digit distributions can distinguish between cell types and tissue origins. Moreover, despite the simplicity of this novel feature-selection method, its separation accuracy is higher than that of the mean-expression level approach and is similar to that of the differential expression approach. Thus, the BL can be used to obtain biological insights from massive amounts of numerical genomics data—a capability that could be utilized in various biomedical applications, e.g., to resolve samples of unknown primary origin, identify possible sample contaminations, and provide insights into the molecular basis of cancer subtypes.
Collapse
|
37
|
Afik S, Raulet G, Yosef N. Reconstructing B-cell receptor sequences from short-read single-cell RNA sequencing with BRAPeS. Life Sci Alliance 2019; 2:2/4/e201900371. [PMID: 31451449 PMCID: PMC6709718 DOI: 10.26508/lsa.201900371] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2019] [Revised: 08/13/2019] [Accepted: 08/14/2019] [Indexed: 12/17/2022] Open
Abstract
BRAPeS is a software for B-cell receptor reconstruction in single cells from very short (25–30 bp) read lengths, which achieves similar success rates and accuracy as applying other methods on long reads. RNA sequencing of single B cells provides simultaneous measurements of the cell state and its antigen specificity as determined by the B-cell receptor (BCR). However, to uncover the latter, further reconstruction of the BCR sequence is needed. We present BRAPeS (“BCR Reconstruction Algorithm for Paired-end Single cells” ), an algorithm for reconstructing BCRs from short-read paired-end single-cell RNA sequencing. BRAPeS is accurate and achieves a high success rate even at very short (25 bp) read length, which can decrease the cost and increase the number of cells that can be analyzed compared with long reads. BRAPeS is publicly available at the following link: https://github.com/YosefLab/BRAPeS.
Collapse
Affiliation(s)
- Shaked Afik
- Center for Computational Biology, University of California, Berkeley, Berkeley, CA, USA
| | - Gabriel Raulet
- Department of Computer Science, University of California, Davis, Davis, CA, USA
| | - Nir Yosef
- Center for Computational Biology, University of California, Berkeley, Berkeley, CA, USA .,Department of Electrical Engineering and Computer Science, University of California, Berkeley, Berkeley, CA, USA.,Ragon Institute of Massachusetts General Hospital, Massachusetts Institute of Technology and Harvard, Cambridge, MA, USA.,Chan Zuckerberg Biohub, San Francisco, CA, USA
| |
Collapse
|
38
|
Abstract
Single-cell RNA-seq (scRNA-seq) has provided novel routes to investigate the heterogeneous populations of T cells and is rapidly becoming a common tool for molecular profiling and identification of novel subsets and functions. This chapter offers an experimental and computational workflow for scRNA-seq analysis of T cells. We focus on the analyses of scRNA-seq data derived from plate-based sorted T cells using flow cytometry and full-length transcriptome protocols such as Smart-Seq2. However, the proposed pipeline can be applied to other high-throughput approaches such as UMI-based methods. We describe a detailed bioinformatics pipeline that can be easily reproduced and discuss future directions and current limitations of these methods in the context of T cell biology.
Collapse
|
39
|
Singh M, Al-Eryani G, Carswell S, Ferguson JM, Blackburn J, Barton K, Roden D, Luciani F, Giang Phan T, Junankar S, Jackson K, Goodnow CC, Smith MA, Swarbrick A. High-throughput targeted long-read single cell sequencing reveals the clonal and transcriptional landscape of lymphocytes. Nat Commun 2019; 10:3120. [PMID: 31311926 PMCID: PMC6635368 DOI: 10.1038/s41467-019-11049-4] [Citation(s) in RCA: 166] [Impact Index Per Article: 33.2] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2018] [Accepted: 06/10/2019] [Indexed: 01/08/2023] Open
Abstract
High-throughput single-cell RNA sequencing is a powerful technique but only generates short reads from one end of a cDNA template, limiting the reconstruction of highly diverse sequences such as antigen receptors. To overcome this limitation, we combined targeted capture and long-read sequencing of T-cell-receptor (TCR) and B-cell-receptor (BCR) mRNA transcripts with short-read transcriptome profiling of barcoded single-cell libraries generated by droplet-based partitioning. We show that Repertoire and Gene Expression by Sequencing (RAGE-Seq) can generate accurate full-length antigen receptor sequences at nucleotide resolution, infer B-cell clonal evolution and identify alternatively spliced BCR transcripts. We apply RAGE-Seq to 7138 cells sampled from the primary tumor and draining lymph node of a breast cancer patient to track transcriptome profiles of expanded lymphocyte clones across tissues. Our results demonstrate that RAGE-Seq is a powerful method for tracking the clonal evolution from large numbers of lymphocytes applicable to the study of immunity, autoimmunity and cancer. Single cell RNA sequencing generates short reads from one end of a template, providing incomplete transcript coverage and limiting identification of diverse sequences such as antigen receptors. Here the authors combine long read nanopore sequencing with short read profiling of barcoded libraries to generate full-length antigen receptor sequences.
Collapse
Affiliation(s)
- Mandeep Singh
- Garvan Institute of Medical Research, Sydney, NSW, 2010, Australia.,St Vincent's Clinical School, Faculty of Medicine, UNSW, Sydney, NSW, 2010, Australia
| | - Ghamdan Al-Eryani
- Garvan Institute of Medical Research, Sydney, NSW, 2010, Australia.,St Vincent's Clinical School, Faculty of Medicine, UNSW, Sydney, NSW, 2010, Australia
| | - Shaun Carswell
- Garvan Institute of Medical Research, Sydney, NSW, 2010, Australia
| | - James M Ferguson
- Garvan Institute of Medical Research, Sydney, NSW, 2010, Australia
| | - James Blackburn
- Garvan Institute of Medical Research, Sydney, NSW, 2010, Australia.,St Vincent's Clinical School, Faculty of Medicine, UNSW, Sydney, NSW, 2010, Australia
| | - Kirston Barton
- Garvan Institute of Medical Research, Sydney, NSW, 2010, Australia.,St Vincent's Clinical School, Faculty of Medicine, UNSW, Sydney, NSW, 2010, Australia
| | - Daniel Roden
- Garvan Institute of Medical Research, Sydney, NSW, 2010, Australia.,St Vincent's Clinical School, Faculty of Medicine, UNSW, Sydney, NSW, 2010, Australia
| | - Fabio Luciani
- St Vincent's Clinical School, Faculty of Medicine, UNSW, Sydney, NSW, 2010, Australia.,Kirby Institute for Infection and Immunity, School of Medical Sciences, UNSW, Sydney, NSW, 2052, Australia
| | - Tri Giang Phan
- Garvan Institute of Medical Research, Sydney, NSW, 2010, Australia.,St Vincent's Clinical School, Faculty of Medicine, UNSW, Sydney, NSW, 2010, Australia
| | - Simon Junankar
- Garvan Institute of Medical Research, Sydney, NSW, 2010, Australia.,St Vincent's Clinical School, Faculty of Medicine, UNSW, Sydney, NSW, 2010, Australia
| | - Katherine Jackson
- Garvan Institute of Medical Research, Sydney, NSW, 2010, Australia.,St Vincent's Clinical School, Faculty of Medicine, UNSW, Sydney, NSW, 2010, Australia
| | - Christopher C Goodnow
- Garvan Institute of Medical Research, Sydney, NSW, 2010, Australia. .,St Vincent's Clinical School, Faculty of Medicine, UNSW, Sydney, NSW, 2010, Australia.
| | - Martin A Smith
- Garvan Institute of Medical Research, Sydney, NSW, 2010, Australia. .,St Vincent's Clinical School, Faculty of Medicine, UNSW, Sydney, NSW, 2010, Australia.
| | - Alexander Swarbrick
- Garvan Institute of Medical Research, Sydney, NSW, 2010, Australia. .,St Vincent's Clinical School, Faculty of Medicine, UNSW, Sydney, NSW, 2010, Australia.
| |
Collapse
|
40
|
Abstract
MOTIVATION Single-cell RNA sequencing (scRNA-seq) has revolutionized biological sciences by revealing genome-wide gene expression levels within individual cells. However, a critical challenge faced by researchers is how to optimize the choices of sequencing platforms, sequencing depths and cell numbers in designing scRNA-seq experiments, so as to balance the exploration of the depth and breadth of transcriptome information. RESULTS Here we present a flexible and robust simulator, scDesign, the first statistical framework for researchers to quantitatively assess practical scRNA-seq experimental design in the context of differential gene expression analysis. In addition to experimental design, scDesign also assists computational method development by generating high-quality synthetic scRNA-seq datasets under customized experimental settings. In an evaluation based on 17 cell types and 6 different protocols, scDesign outperformed four state-of-the-art scRNA-seq simulation methods and led to rational experimental design. In addition, scDesign demonstrates reproducibility across biological replicates and independent studies. We also discuss the performance of multiple differential expression and dimension reduction methods based on the protocol-dependent scRNA-seq data generated by scDesign. scDesign is expected to be an effective bioinformatic tool that assists rational scRNA-seq experimental design and comparison of scRNA-seq computational methods based on specific research goals. AVAILABILITY AND IMPLEMENTATION We have implemented our method in the R package scDesign, which is freely available at https://github.com/Vivianstats/scDesign. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Wei Vivian Li
- Department of Statistics, University of California, Los Angeles, CA, USA
| | - Jingyi Jessica Li
- Department of Statistics, University of California, Los Angeles, CA, USA
- Department of Human Genetics, University of California, Los Angeles, CA, USA
| |
Collapse
|
41
|
Lindeman I, Stubbington MJT. Antigen Receptor Sequence Reconstruction and Clonality Inference from scRNA-Seq Data. Methods Mol Biol 2019; 1935:223-249. [PMID: 30758830 DOI: 10.1007/978-1-4939-9057-3_15] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
In this chapter, we describe TraCeR and BraCeR, our computational tools for reconstruction of paired full-length antigen receptor sequences and clonality inference from single-cell RNA-seq (scRNA-seq) data. In brief, TraCeR reconstructs T-cell receptor (TCR) sequences from scRNA-seq data by extracting sequencing reads derived from TCRs by aligning the reads from each cell against synthetic TCR sequences. TCR-derived reads are then assembled into full-length recombined TCR sequences. BraCeR builds on the TraCeR pipeline and accounts for somatic hypermutations (SHM) and isotype switching. Here we discuss experimental design, use of the tools, and interpretation of the results.
Collapse
Affiliation(s)
- Ida Lindeman
- Wellcome Sanger Institute, Hinxton, Cambridge, UK
- KG Jebsen Coeliac Disease Research Centre and Department of Immunology, University of Oslo, Oslo, Norway
| | | |
Collapse
|
42
|
Abstract
During the last decade, ncRNAs have been investigated intensively and revealed their regulatory role in various biological processes. Worldwide research efforts have identified numerous ncRNAs and multiple RNA subtypes, which are attributed to diverse functionalities known to interact with different functional layers, from DNA and RNA to proteins. This makes the prediction of functions for newly identified ncRNAs challenging. Current bioinformatics and systems biology approaches show promising results to facilitate an identification of these diverse ncRNA functionalities. Here, we review (a) current experimental protocols, i.e., for Next Generation Sequencing, for a successful identification of ncRNAs; (b) sequencing data analysis workflows as well as available computational environments; and (c) state-of-the-art approaches to functionally characterize ncRNAs, e.g., by means of transcriptome-wide association studies, molecular network analyses, or artificial intelligence guided prediction. In addition, we present a strategy to cover the identification and functional characterization of unknown transcripts by using connective workflows.
Collapse
|
43
|
AutoImpute: Autoencoder based imputation of single-cell RNA-seq data. Sci Rep 2018; 8:16329. [PMID: 30397240 PMCID: PMC6218547 DOI: 10.1038/s41598-018-34688-x] [Citation(s) in RCA: 82] [Impact Index Per Article: 13.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2018] [Accepted: 10/15/2018] [Indexed: 12/22/2022] Open
Abstract
The emergence of single-cell RNA sequencing (scRNA-seq) technologies has enabled us to measure the expression levels of thousands of genes at single-cell resolution. However, insufficient quantities of starting RNA in the individual cells cause significant dropout events, introducing a large number of zero counts in the expression matrix. To circumvent this, we developed an autoencoder-based sparse gene expression matrix imputation method. AutoImpute, which learns the inherent distribution of the input scRNA-seq data and imputes the missing values accordingly with minimal modification to the biologically silent genes. When tested on real scRNA-seq datasets, AutoImpute performed competitively wrt., the existing single-cell imputation methods, on the grounds of expression recovery from subsampled data, cell-clustering accuracy, variance stabilization and cell-type separability.
Collapse
|
44
|
AlJanahi AA, Danielsen M, Dunbar CE. An Introduction to the Analysis of Single-Cell RNA-Sequencing Data. Mol Ther Methods Clin Dev 2018; 10:189-196. [PMID: 30094294 PMCID: PMC6072887 DOI: 10.1016/j.omtm.2018.07.003] [Citation(s) in RCA: 67] [Impact Index Per Article: 11.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
The recent development of single-cell RNA sequencing has deepened our understanding of the cell as a functional unit, providing new insights based on gene expression profiles of hundreds to hundreds of thousands of individual cells, and revealing new populations of cells with distinct gene expression profiles previously hidden within analyses of gene expression performed on bulk cell populations. However, appropriate analysis and utilization of the massive amounts of data generated from single-cell RNA sequencing experiments are challenging and require an understanding of the experimental and computational pathways taken between preparation of input cells and output of interpretable data. In this review, we will discuss the basic principles of these new technologies, focusing on concepts important in the analysis of single-cell RNA-sequencing data. Specifically, we summarize approaches to quality-control measures for determination of which single cells to include for further examination, methods of data normalization and scaling to overcome the relatively inefficient capture rate of mRNA from each cell, and clustering and visualization algorithms used for dimensional reduction of the data to a two-dimensional plot.
Collapse
Affiliation(s)
- Aisha A. AlJanahi
- Translational Stem Cell Biology Branch, NHLBI, NIH, Bethesda, MD, USA
- Department of Biochemistry and Molecular & Cellular Biology, Georgetown University Medical Center, Washington, DC, USA
| | - Mark Danielsen
- Department of Biochemistry and Molecular & Cellular Biology, Georgetown University Medical Center, Washington, DC, USA
| | - Cynthia E. Dunbar
- Translational Stem Cell Biology Branch, NHLBI, NIH, Bethesda, MD, USA
| |
Collapse
|
45
|
Venturi V, Thomas PG. The expanding role of systems immunology in decoding the T cell receptor repertoire. ACTA ACUST UNITED AC 2018; 12:37-45. [PMID: 31106281 DOI: 10.1016/j.coisb.2018.09.005] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
T cells play a crucial role in the immune system's defense against many infectious diseases, including persistent infections for which no effective vaccines currently exist. The T cell component of the adaptive immune system is highly complex involving a constantly evolving landscape of various inter-related T cell populations. These T cell populations are characterized by their phenotypic and functional properties as well as the collection, or repertoire, of T cell receptors (TCR) that mediate T cell recognition of antigenic peptides derived from pathogens. Understanding the various processes and factors that impact the development and evolution of the broader T cell repertoire available to recognize and respond to pathogens and the characteristics of antigen-experienced T cell repertoires associated with effective immune control of pathogens is critical to the rational design of T cell-based vaccines and therapies. In this article we discuss, using examples of recent research, the promise that systems immunology approaches, involving quantitative analysis and mathematical and computational modeling of immunological data, hold for decoding the complex TCR repertoire system in the current era of advancing technologies.
Collapse
Affiliation(s)
- Vanessa Venturi
- Infection Analytics Program, Kirby Institute for Infection and Immunity, UNSW Australia, Sydney, NSW, Australia
| | - Paul G Thomas
- Department of Immunology, St. Jude Children's Research Hospital, Memphis, TN, USA
| |
Collapse
|
46
|
Nguyen A, Khoo WH, Moran I, Croucher PI, Phan TG. Single Cell RNA Sequencing of Rare Immune Cell Populations. Front Immunol 2018; 9:1553. [PMID: 30022984 PMCID: PMC6039576 DOI: 10.3389/fimmu.2018.01553] [Citation(s) in RCA: 76] [Impact Index Per Article: 12.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2018] [Accepted: 06/22/2018] [Indexed: 11/29/2022] Open
Abstract
Single-cell RNA sequencing (scRNA-Seq) is transforming our ability to characterize cells, particularly rare cells that are often overlooked in bulk population analytical approaches. This has lead to the discovery of new cell types and cellular states that echo the underlying heterogeneity and plasticity in the immune system. Technologies for the capture, sequencing, and bioinformatic analysis of single cells are rapidly improving, and scRNA-Seq is now becoming much more accessible to non-specialized laboratories. Here, we describe our experiences in adopting scRNA-Seq to the study of rare immune cells in their microanatomical niches.
Collapse
Affiliation(s)
- Akira Nguyen
- Immunology Division, Garvan Institute of Medical Research, Darlinghurst, NSW, Australia
- St Vincent’s Clinical School, Faculty of Medicine, University of New South Wales, Darlinghurst, NSW, Australia
| | - Weng Hua Khoo
- Bone Biology Division, Garvan Institute of Medical Research, Darlinghurst, NSW, Australia
- School of Biotechnology and Biomolecular Sciences, University of New South Wales, Kensington, NSW, Australia
| | - Imogen Moran
- Immunology Division, Garvan Institute of Medical Research, Darlinghurst, NSW, Australia
- St Vincent’s Clinical School, Faculty of Medicine, University of New South Wales, Darlinghurst, NSW, Australia
| | - Peter I. Croucher
- St Vincent’s Clinical School, Faculty of Medicine, University of New South Wales, Darlinghurst, NSW, Australia
- Bone Biology Division, Garvan Institute of Medical Research, Darlinghurst, NSW, Australia
| | - Tri Giang Phan
- Immunology Division, Garvan Institute of Medical Research, Darlinghurst, NSW, Australia
- St Vincent’s Clinical School, Faculty of Medicine, University of New South Wales, Darlinghurst, NSW, Australia
| |
Collapse
|
47
|
Wang Z, Zhu L, Nguyen THO, Wan Y, Sant S, Quiñones-Parra SM, Crawford JC, Eltahla AA, Rizzetto S, Bull RA, Qiu C, Koutsakos M, Clemens EB, Loh L, Chen T, Liu L, Cao P, Ren Y, Kedzierski L, Kotsimbos T, McCaw JM, La Gruta NL, Turner SJ, Cheng AC, Luciani F, Zhang X, Doherty PC, Thomas PG, Xu J, Kedzierska K. Clonally diverse CD38 +HLA-DR +CD8 + T cells persist during fatal H7N9 disease. Nat Commun 2018; 9:824. [PMID: 29483513 PMCID: PMC5827521 DOI: 10.1038/s41467-018-03243-7] [Citation(s) in RCA: 84] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2017] [Accepted: 01/30/2018] [Indexed: 12/21/2022] Open
Abstract
Severe influenza A virus (IAV) infection is associated with immune dysfunction. Here, we show circulating CD8+ T-cell profiles from patients hospitalized with avian H7N9, seasonal IAV, and influenza vaccinees. Patient survival reflects an early, transient prevalence of highly activated CD38+HLA-DR+PD-1+ CD8+ T cells, whereas the prolonged persistence of this set is found in ultimately fatal cases. Single-cell T cell receptor (TCR)-αβ analyses of activated CD38+HLA-DR+CD8+ T cells show similar TCRαβ diversity but differential clonal expansion kinetics in surviving and fatal H7N9 patients. Delayed clonal expansion associated with an early dichotomy at a transcriptome level (as detected by single-cell RNAseq) is found in CD38+HLA-DR+CD8+ T cells from patients who succumbed to the disease, suggesting a divergent differentiation pathway of CD38+HLA-DR+CD8+ T cells from the outset during fatal disease. Our study proposes that effective expansion of cross-reactive influenza-specific TCRαβ clonotypes with appropriate transcriptome signatures is needed for early protection against severe influenza disease.
Collapse
MESH Headings
- ADP-ribosyl Cyclase 1/genetics
- ADP-ribosyl Cyclase 1/immunology
- CD8-Positive T-Lymphocytes/immunology
- CD8-Positive T-Lymphocytes/pathology
- CD8-Positive T-Lymphocytes/virology
- Clonal Selection, Antigen-Mediated/genetics
- Cohort Studies
- Critical Illness
- Gene Expression Regulation
- HLA-DR Antigens/genetics
- HLA-DR Antigens/immunology
- Hospitalization
- Humans
- Influenza A Virus, H7N9 Subtype/immunology
- Influenza A Virus, H7N9 Subtype/pathogenicity
- Influenza, Human/genetics
- Influenza, Human/immunology
- Influenza, Human/mortality
- Influenza, Human/virology
- Lymphocyte Activation
- Membrane Glycoproteins/genetics
- Membrane Glycoproteins/immunology
- Programmed Cell Death 1 Receptor/genetics
- Programmed Cell Death 1 Receptor/immunology
- Receptors, Antigen, T-Cell, alpha-beta/genetics
- Receptors, Antigen, T-Cell, alpha-beta/immunology
- Survival Analysis
- T-Lymphocyte Subsets/immunology
- T-Lymphocyte Subsets/pathology
- T-Lymphocyte Subsets/virology
- Transcriptome/immunology
Collapse
Affiliation(s)
- Zhongfang Wang
- Shanghai Public Health Clinical Center & Institutes of Biomedical Sciences, Key Laboratory of Medical Molecular Virology of Ministry of Education/Health, Shanghai Medical College, Fudan University, 201508, Shangai, China
- Department of Microbiology and Immunology, University of Melbourne, The Peter Doherty Institute for Infection and Immunity, Parkville, Melbourne, VIC, 3010, Australia
| | - Lingyan Zhu
- Shanghai Public Health Clinical Center & Institutes of Biomedical Sciences, Key Laboratory of Medical Molecular Virology of Ministry of Education/Health, Shanghai Medical College, Fudan University, 201508, Shangai, China
| | - Thi H O Nguyen
- Department of Microbiology and Immunology, University of Melbourne, The Peter Doherty Institute for Infection and Immunity, Parkville, Melbourne, VIC, 3010, Australia
| | - Yanmin Wan
- Shanghai Public Health Clinical Center & Institutes of Biomedical Sciences, Key Laboratory of Medical Molecular Virology of Ministry of Education/Health, Shanghai Medical College, Fudan University, 201508, Shangai, China
| | - Sneha Sant
- Department of Microbiology and Immunology, University of Melbourne, The Peter Doherty Institute for Infection and Immunity, Parkville, Melbourne, VIC, 3010, Australia
| | - Sergio M Quiñones-Parra
- Department of Microbiology and Immunology, University of Melbourne, The Peter Doherty Institute for Infection and Immunity, Parkville, Melbourne, VIC, 3010, Australia
| | - Jeremy Chase Crawford
- Department of Immunology, St. Jude Children's Research Hospital, Memphis, TN, 38105, USA
| | - Auda A Eltahla
- School of Medical Sciences and The Kirby Institute, UNSW Sydney, Sydney, NSW, 2052, Australia
| | - Simone Rizzetto
- School of Medical Sciences and The Kirby Institute, UNSW Sydney, Sydney, NSW, 2052, Australia
| | - Rowena A Bull
- School of Medical Sciences and The Kirby Institute, UNSW Sydney, Sydney, NSW, 2052, Australia
| | - Chenli Qiu
- Shanghai Public Health Clinical Center & Institutes of Biomedical Sciences, Key Laboratory of Medical Molecular Virology of Ministry of Education/Health, Shanghai Medical College, Fudan University, 201508, Shangai, China
| | - Marios Koutsakos
- Department of Microbiology and Immunology, University of Melbourne, The Peter Doherty Institute for Infection and Immunity, Parkville, Melbourne, VIC, 3010, Australia
| | - E Bridie Clemens
- Department of Microbiology and Immunology, University of Melbourne, The Peter Doherty Institute for Infection and Immunity, Parkville, Melbourne, VIC, 3010, Australia
| | - Liyen Loh
- Department of Microbiology and Immunology, University of Melbourne, The Peter Doherty Institute for Infection and Immunity, Parkville, Melbourne, VIC, 3010, Australia
| | - Tianyue Chen
- Shanghai Public Health Clinical Center & Institutes of Biomedical Sciences, Key Laboratory of Medical Molecular Virology of Ministry of Education/Health, Shanghai Medical College, Fudan University, 201508, Shangai, China
| | - Lu Liu
- Shanghai Public Health Clinical Center & Institutes of Biomedical Sciences, Key Laboratory of Medical Molecular Virology of Ministry of Education/Health, Shanghai Medical College, Fudan University, 201508, Shangai, China
| | - Pengxing Cao
- School of Mathematics and Statistics, University of Melbourne, Parkville, VIC, 3010, Australia
| | - Yanqin Ren
- Shanghai Public Health Clinical Center & Institutes of Biomedical Sciences, Key Laboratory of Medical Molecular Virology of Ministry of Education/Health, Shanghai Medical College, Fudan University, 201508, Shangai, China
| | - Lukasz Kedzierski
- Department of Microbiology and Immunology, University of Melbourne, The Peter Doherty Institute for Infection and Immunity, Parkville, Melbourne, VIC, 3010, Australia
| | - Tom Kotsimbos
- Department of Respiratory Medicine, Alfred Hospital Health and Department Medicine, Monash University, 99 Commercial Road, Melbourne, VIC, 3004, Australia
| | - James M McCaw
- School of Mathematics and Statistics, University of Melbourne, Parkville, VIC, 3010, Australia
| | - Nicole L La Gruta
- Department of Microbiology and Immunology, University of Melbourne, The Peter Doherty Institute for Infection and Immunity, Parkville, Melbourne, VIC, 3010, Australia
- Department of Biochemistry and Molecular Biology, Infection and Immunity Program, Biomedicine Discovery Institute, Monash University, Clayton, VIC, 3800, Australia
| | - Stephen J Turner
- Department of Microbiology and Immunology, University of Melbourne, The Peter Doherty Institute for Infection and Immunity, Parkville, Melbourne, VIC, 3010, Australia
- Department of Microbiology, Biomedical Discovery Institute, Monash University, Clayton, VIC, 3800, Australia
| | - Allen C Cheng
- Infection Prevention and Healthcare Epidemiology Unit, Alfred Health and School of Public Health and Preventive Medicine, 99 Commercial Road, Melbourne, VIC, 3004, Australia
| | - Fabio Luciani
- School of Medical Sciences and The Kirby Institute, UNSW Sydney, Sydney, NSW, 2052, Australia
| | - Xiaoyan Zhang
- Shanghai Public Health Clinical Center & Institutes of Biomedical Sciences, Key Laboratory of Medical Molecular Virology of Ministry of Education/Health, Shanghai Medical College, Fudan University, 201508, Shangai, China
| | - Peter C Doherty
- Department of Microbiology and Immunology, University of Melbourne, The Peter Doherty Institute for Infection and Immunity, Parkville, Melbourne, VIC, 3010, Australia
- Department of Immunology, St. Jude Children's Research Hospital, Memphis, TN, 38105, USA
| | - Paul G Thomas
- Department of Immunology, St. Jude Children's Research Hospital, Memphis, TN, 38105, USA
| | - Jianqing Xu
- Shanghai Public Health Clinical Center & Institutes of Biomedical Sciences, Key Laboratory of Medical Molecular Virology of Ministry of Education/Health, Shanghai Medical College, Fudan University, 201508, Shangai, China.
| | - Katherine Kedzierska
- Shanghai Public Health Clinical Center & Institutes of Biomedical Sciences, Key Laboratory of Medical Molecular Virology of Ministry of Education/Health, Shanghai Medical College, Fudan University, 201508, Shangai, China.
- Department of Microbiology and Immunology, University of Melbourne, The Peter Doherty Institute for Infection and Immunity, Parkville, Melbourne, VIC, 3010, Australia.
| |
Collapse
|