- Open Access
Advances in single-cell sequencing: insights from organ transplantation
Military Medical Research volume 8, Article number: 45 (2021)
Single-cell RNA sequencing (scRNA-seq) is a comprehensive technical tool to analyze intracellular and intercellular interaction data by whole transcriptional profile analysis. Here, we describe the application in biomedical research, focusing on the immune system during organ transplantation and rejection. Unlike conventional transcriptome analysis, this method provides a full map of multiple cell populations in one specific tissue and presents a dynamic and transient unbiased method to explore the progression of allograft dysfunction, starting from the stress response to final graft failure. This promising sequencing technology remarkably improves individualized organ rejection treatment by identifying decisive cellular subgroups and cell-specific interactions.
The demand for organs due to end-stage organ failure is permanently increasing. The best solution is autologous transplantation, enabled by patient-derived induced pluripotent stem cells. This approach is currently envisaged for single cells and basic cellular units such as islets. Tissue engineering is being optimized as one solution for larger organs. However, this option is not realistic presently. Xenotransplantation could soon become available clinically and an alternative backup for allotransplantation. However, both will depend on matching grafts and require immune suppression. The importance of understanding immune reactions, identifying cellular subtypes involved in graft acceptance and tolerance induction and identifying early indicators for rejection mechanisms require detailed immune cell profiling.
The immune system is a complex biological network comprising multiple layers of orchestrated genes, proteins and cells. In response to the challenge of pathogens or transplants, the immune system triggers various interactions between immune cells and other cells, provoking specific responses . Innate and adaptive immune cells interact to ensure tissue protection according to functional requirements. Disruption of immune system homeostasis causes transplant rejection [2, 3]. Although inhibition of acute rejection has improved significantly in the past two decades, long-term rejection and immunosuppression can lead to high morbidity and mortality, and chronic transplant rejection can cause irreversible damage to the graft [4, 5]. The most common clinical acute rejection is mainly mediated by cellular responses, while hyperacute rejection and chronic rejection are mainly mediated by humoral immunity. The transplant rejection mechanism is an immunological reaction that recognizes foreign molecules of the donor cells, triggering attack and destruction. Various immunological factors are involved in acute and chronic rejection, including human leukocyte antigen (HLA) mismatch, donor-specific antibodies and non-immune factors (e.g., donor age, infection, and immunosuppressive drug toxicity) [6, 7].
Assessing the immune reaction
Cellular and molecular assays to measure immune cell differentiation, cellular function and antigen specificity contribute to solving important problems in immune-mediated transplant rejection. Fluorescence-activated cell sorting (FACS) is used for post-transplantation immune surveillance to identify cells of the immune system by detecting surface marker expression or intracellular proteins, including cytokines. This tool has been supplemented by in situ histological tests [8, 9]. However, histological diagnosis may miss subtle alterations in individual patients, which are essential for rejective pathology . Therefore, the combination of transcript sets and histological diagnosis of tissue samples is used to predict antibody-mediated or cellular rejections [11, 12]. Additionally, bulk RNA sequencing of allograft biopsies is a method to determine gene-specific single-nucleotide variants of donors and recipients . Microarray technology has been used to study the pathogenesis of transplant rejection processes. For example, Roedder et al.  used microarray technology to determine a test set of 17 relevant genes to predict renal rejection. Each of these approaches provides valuable insights. However, they are subject to limitations due to the complex immune rejection response. Furthermore, extensive analysis cannot resolve phenotypic heterogeneity and distinguish the gene profile of donors and recipients in mixed cell populations . Even with the application of Cytometry by Time of Flight (CyTOF), it remains challenging to assess a truly unbiased dataset of a single cell, a process that requires single-cell resolution . Among single-cell profiling methods, scRNA-seq comprehensively measures the transcriptional expression of bulk cells [1, 17] and quantitatively analyzes all transcripts expressed in a single cell, providing an unbiased strategy to identify and characterize different cell populations [18,19,20].
This review discusses the development and application of scRNA-seq in organ transplantation. This cutting-edge technology will improve immunotherapies and help to predict recipient outcomes.
Sample harvest and tissue processing
The entire process starting with sample isolation to the final evaluation of scRNA-seq data is summarized in Fig. 1. The first step of scRNA-seq analysis is the dissociation of the graft tissue, which is obtained in most cases by biopsies. The currently used protocols were developed and improved over decades, and each has its strengths and weaknesses. For transplanted organs, laser capture microdissection, digestion or enzyme-related approaches, followed by density gradient centrifugation or fluorescence-activated cell sorting are used. Tissue dissociation is more difficult to implement for frozen samples, e.g., those of the liver . In particular, for liver samples, a different sensitivity to cell death caused by the dissociation step may result in bias because hepatocytes die very fast, while other cells become activated during tissue dissociation, indicating a transcriptional stress response. Macparland et al.  developed a milder approach to reduce the cell damage rates by maintaining the tissue at 4 °C during all steps, including collagenase perfusion. Wang et al.  used a hypothermic strategy for kidney preservation for up to 4 days for scRNA-seq analysis. Recently, Guillaumet-Adkins et al.  published an improved method for sample preservation, gradual freezing by 1 ℃ per minute to -80 ℃, that does not change the transcriptional profiles and makes cryopreserved cells and tissues applicable for scRNA-seq. However, frozen tissues cannot be processed like fresh tissues. Tissue dissociation leads to the loss of spatial and anatomical information for cells, and the precise location of each cell should not be ignored. To address this issue, RNA probes identifying cellular transcriptional organization  or spatial transcriptomic protocols  are both helpful alternative methods.
The next steps after tissue dissociation include RNA capture, reverse transcription, RNA sequencing, and library construction. Selection of a suitable sequencing method is challenging because several methods exist, such as CEL-seq2, Drop-seq, MARS-seq, MATQ-seq, Quartz-seq, SCRB-seq, Smart-seq, Smart-seq2, Drop-seq, FB5P-seq, SPLIT-seq, and DNBelab C4 [27,28,29,30]. We summarize the most commonly used scRNA-seq methods in Table 1 according to their capturing format, cDNA amplification, sequencing method, transcript coverage advantages and limitations. Microfluidic technologies for scRNA-seq involve droplet-based and plate-based technologies. Amplification is performed by PCR for Smart-seq [31, 32] and Quartz-seq [33, 34] and in vitro transcription, generating RNA in vitro, e.g., by InDrop  and CEL-seq [36, 37]. Drop-seq , InDrop , and CEL-seq [36, 37] incorporate unique molecular identifiers (UMIs) into cDNA. A UMI is a short sequence barcode to detect and quantify transcripts. These molecular barcodes uniquely tag each molecule in a sample library and reduce quantitative and error biases introduced by amplification. Smart-seq2, Quartz-seq, and MATQ-seq produce almost full-length sequencing data, while others (e.g., CEL-seq and Drop-seq) only capture the 3'-end sequence or 5'-sequence (e.g., FB5P-seq and STRT-seq) [10, 30, 39, 40]. Each platform provides multiple and specific but not completely comprehensive advantages in data capture. The reduction of mRNA amplification noise by CEL-seq2, InDrop, Drop-seq, MARS-seq, and SCRB-seq is a favorable feature that makes these platforms preferable. However, MARS-seq, SCRB-seq, and particularly Smart-seq2 platforms capture more genes using the same number of cells, making them preferable for relatively low quantities . Drop-seq analyzes thousands of individual cells simultaneously without losing the original transcript . Compared with other widely used single-cell RNA sequencing platforms (such as Smart-seq2), 10× Genomics Chromium is a more cost-effective and time-efficient system. This platform generates droplets and forms a single-cell suspension. Additionally, it can process many cells and detect even rare cell types or transcripts  by combining one of the following methods: InDrop for rare cell populations  or CEL-seq for complex tissues containing multiple cell populations [36, 37]. Smart-seq increases the thermal stability of DNA base pairs [31, 32]. MATQ-seq is implemented on low-abundance genes and noncoding and non-polyadenylated RNA . Quartz-seq is applied to detect the different cell cycle phases and transcriptome heterogeneity. SCRB-seq is used for heterogeneous populations [33, 34]. FB5P-seq  and T-cell receptor repertoire sequencing (TCR-seq)  are used to identify the repertoire and diversity of BCRs and TCRs. STRT-seq tracks the cell origin efficiency without quantitative bias against long transcripts. One of the common disadvantages is the limited throughput and read coverage, e.g., by Smart-seq 1 and 2 [31, 32]. Another is the requirement for skilled technicians, e.g., for Quartz-seq 1 and 2 [33, 34]. Because each method has weaknesses (Table 1), investigators must choose a platform according to their specific interests.
Data analysis after sequencing generally comprises quality control (QC), batch effect correction, normalization, data imputation, dimensionality adjustment, subsequent expression analysis and cell subpopulation identification . QC is required for high technical noise and low-quality data, which are often generated because of dead cells. Batch effects are caused by large-scale scRNA-seq datasets, samples prepared in different laboratories, even those based on the same protocol, and large data generated from separate technicians or different time points [49, 50]. They can cause systematic errors and differences among multiple transcriptional profiles. Thus, batch effect correction is vital to avoid this misinterpretation, and removing unwanted variation (RUV) is a good normalization strategy to adjust confounding technical effects . Data normalization is necessary to adjust technical biases. Two types of normalization are used: sample normalization, which adjusts within-sample differences, and gene profile normalization, which eliminates gene expression biases. Data imputation is an effective strategy to insert substituted values into dropouts, eliminating the influence of missing data. Because of technological developments in scRNA-seq and bulk data generation with high dimensionalities, computational bioinformatics analysis is essential to process raw data. Automatic annotation methods such as the "SingleR" package , online databases  and gene expression markers  can be used for cell marker identification. The commonly used "SingleR" labels new cells from the scRNA-seq dataset based on similarity to the reference dataset of samples with known labels.
Challenges during data analysis involve multiple biases in the entire procedure and high dimensional datasets . The possibility of low-quality data or dropouts in the scRNA-seq results, caused by low viability or dead cells, can lead to misinterpretation . An optimal method can avoid false results to enable correct transcriptional interpretation (Table 2). "Seurat" is an R package designed for QC, analysis and exploration of scRNA-seq data to reduce some biases [54, 57]. To distinguish technical biases from biological signals, Bayesian Inference for Single-cell ClUstering and ImpuTing (BISCUIT) provides a discernible advantage for graphical algorithms. BISCUIT imputes dropouts and improves both clustering and normalization . Using these tools, the data are clustered to reduce technical variations (amplification bias, sequencing depth, GC content, capture inefficiency, and RNA content variations ). Because scRNA-seq data contain many genes and cells in high dimensions, large-scale computational resources are required. T-distributed stochastic neighbor embedding (tSNE) is an approach for high-dimensional data because of the computation time and potentially manifold embeddings for the same datasets from run to run . PhenoGraph is also an algorithmic method used for partitions of high-dimensional data into distinct subgroups within complex tissues . Dimensional reduction and uniform manifold approximation and projection (UMAP) are alternative methods for data analysis developed to project the data into lower dimensions and visualize cell clusters with high reproducibility without losing the properties of the original data . Xu et al. described a new clustering algorithm of graph-based shared nearest neighbor (SNN)-Cliq implemented in Python and MATLAB software, considering the surrounding data points, including low-density region data and detecting more cell subtypes with high accuracy . For zero-inflated data comprising an excessive number of zeros, zero inflated factor analysis (ZIFA) is a novel approach to reduce dataset dimensions . "Harmony" is another R package with an efficient batch algorithm method to integrate multiple datasets and requires fewer computational resources. Korsunsky et al.  developed this method (https://github.com/immunogenomics/harmony) and applied it to large datasets integrated with spatial transcriptomics data.
However, these algorithms are most commonly used for static analysis. Another promising means, considering a continuous transition between different cellular states, are the machine learning approaches listed in Table 3, Monocle , Monocle2 , TSCAN , Wishbone , Slingshot , Diffusion pseudotime software  and Wanderlust  allow the construction of trajectories of cells in dynamic gene regulation and explain normal physiological and pathophysiologic alterations of cellular subgroups in specific locations of the human body. Saelens et al.  concluded that Slingshot, TSCAN and Monocle2 exhibited better trajectory identification. Therefore, the combination of personalized medicine and artificial intelligence will become applicable in the near future . This combination will provide a map of cells, considering temporal dynamics and spatial positioning by evaluating the pathological microenvironment, phase of the cell cycle and responses to clinical medication. For transplanted organs, selecting an optimized process for accurate subsequent analysis is highly recommended.
Applications of scRNA-seq for transplantation
Currently, the main difficulties for successful transplantation comprise, among others, best-matched donor selection and a reduction in lifelong immunosuppression . For transplant rejection, the latest advances in scRNA-seq provide an opportunity to fully reveal new cell types and states without result bias and RNA degradation . Snapshots of the single-cell transcriptome exhibit various stages of immune differentiation and activation, while these stages are rarely synchronized among cells . At single-cell resolution, it can describe immune cells, stromal cells and new cell subtypes that undergo immune rejection and further compare the unique characteristics of the signaling pathways between different cellular subgroups . Here, we describe the advantages of scRNA-seq in the fields of kidney, liver, lung and hematopoietic stem cell (HSC) transplantation and for immunological applications.
T cells play a crucial role in graft rejection. Most studies have focused on bulk methods based on biopsy samples without providing information about αβ chain pairing of T-cell receptors (TCRs) . This lack can underestimate actual library differences and fails to reflect that T cells with the same TCR can exert opposite biological functions. ScRNA-seq technology overcomes the abovementioned limitations and brings library analysis a higher diversity . ScRNA-seq detects T-cell subclones . For example, Morris et al.  monitored donor-reactive T cells in patients with tolerant and non-tolerant kidney transplantation. Donor reactivity has been detected in patients with tolerance. The decrease in T-cell clones in non-tolerant patients did not show a reduced number of donor-reactive T cells. Contrary to data from tolerant patients undergoing standard immunosuppression, an expansion of donor-reactive T-cell clones was observed in peripheral blood [81,82,83]. Antibody-mediated rejection (AMR) in the kidney is not easily identified by histological diagnoses. However, a transcriptional expression profile can strengthen the diagnosis. AMR injury is the most common driver of late allograft loss [12, 84]. Several groups have performed scRNA-seq of kidney allograft undergoing AMR, where monocytes, B cells, plasma cells and T cells invade into the tissue  and donor endothelial cells (ECs) are the primary targets of the recipients’ humoral immune response [86, 87]. Thus, scRNA-seq provides a comprehensive understanding of subtle mechanisms in conceptualizing heterogeneous kidney rejection. Macrophages and T cells, activated in the recipient, significantly differ from the original populations in either the donor or recipient. Some of these immune cells exist for only a few days after transplantation, but macrophages can persist for several years . Furthermore, immune cell populations residing in donor-derived tissues can be replaced by recipient cells, particularly during rejection . Liu et al.  presented a novel heterogeneous profile of immune cells based on allograft biopsies and matching healthy kidneys by scRNA-seq, integrating the key alterations of molecular functions, establishing therapeutic surveillance for kidney allograft rejection and improving allograft survival . They identified subclusters of cytotoxic T lymphocytes that exhibit a more proinflammatory role in renal allograft rejection, while activated B cells interacted with surrounding stromal cells, mostly emerging in allograft kidneys, leading to immune cell recruitment and an activated inflammatory response. Non-invasive urine or blood biomarkers are applicable for a limited group of pathologies , whereas invasive biopsies are used to profile non-circulating immune cells in transplant rejection, providing more diagnostic accuracy and prognostic biomarkers amenable to therapeutic tools [10, 85].
The interaction of immune cells and liver cells in a transplant setting is a key mechanism for liver tolerance induction . ScRNA-seq can improve the hepatic immune cellular map in interpreting the specific CD4+ T-cell subgroup from other T cells in liver transplant rejection and tolerance . Immune cells such as dendritic cells (DCs), T cells, and probably NK cells interact with liver sinusoidal ECs and hepatocytes, adding specific signaling molecules to generate a tolerogenic state. Apoptosis of infiltrating T cells may be critical for allograft tolerance [93, 94]. Despite possible tolerance induction after liver transplantation, human liver allografts are likely to be rejected without applying immunosuppressive drugs . Approximately 10%-30% of allograft recipients are diagnosed with acute cellular rejection , likely due to T-cell-related immune responses . Applying scRNA-seq may better elucidate the cellular immune responses in the liver allograft.
Using scRNA-seq to test single cells isolated from lung transplantation donors and lung fibrosis has revealed transcriptionally distinct populations of alveolar macrophages that express profibrotic genes in patients with pulmonary fibrosis [98, 99]. Some mesenchymal cell markers that play a role in Wnt/β-catenin signaling during lung regeneration and some previously described rare cell populations have been identified. These technologies can improve the diagnosis of patients with fibrosis after lung transplantation and can be used to identify patients most likely to benefit from targeted therapy and monitor their response during disease progression . Mould et al.  assessed the cell populations between healthy samples and found highly conserved cellular heterogeneity in bronchoalveolar lavage (BAL) cells. By dynamically comparing the lungs of donors and recipients, persistent donor resident memory T cells are correlated with better clinical outcomes , providing a potential therapeutic tool for extended allograft survival.
To study the rejection of transplanted HSCs, Dong et al.  used scRNA-seq to obtain a transcriptome-based classification of 28 hematopoietic cell types. According to this classification, most transplanted HSCs are dedicated to transcriptional immunophenotypical multipotent progenitors (tMPP1). Transcriptional analysis and functional evaluation showed that the proliferation of transplanted cells is accompanied by a gradually decreased percentage of HSCs . However, a balance between proliferation, differentiation and stem cell maintenance likely exists. Graft-versus-host disease (GVHD) is the main complication of allogeneic hematopoietic cell transplantation (HCT) . Acute GVHD is mediated by donor T cells . TCR-seq, as a type of scRNA-seq, can clarify how acute GVHD occurs. Although donor T-cell pools have highly similar TCRs, the TCR repertoire after HCT is very specific to recipients. TCR recombination is highly stochastic and may not depend on evaluating the most expanded TCR clones in any single recipient but on the complex polyclonal architecture. These results can be used to guide clinical decisions to prevent or treat acute GVHD . By analyzing skin and peripheral blood T cells using TCR-seq, host skin-resident T cells were found to have an unanticipatedly pathogenic impact on GVHD .
Immunological applications of scRNA-seq
ScRNA-seq identifies gene profiles of various cell populations. This technical tool avoids the weakness of bulk analysis, which is likely to miss cell-specific signatures . It also improves the understanding of immune cell ontogeny and interaction with stromal cells in a given organ. The advantages of scRNA are maximized by its combination with databases or improvements by other techniques [such as single-nucleus RNA sequencing (snRNA-seq)] . The Immune Cell Atlas (ICA) is an essential part of the international Human Cell Atlas initiative (https://www.humancellatlas.org/), which collects samples from patients who have undergone rejection response and evaluates different reaction stages using scRNA-seq technology. By visualizing the dynamic observation of cell processes, the subtle transcriptional differences among cell types can be qualitative. The rejection response, caused by specific gene regulation, provides robust target genes and molecular mechanisms to diagnose and treat transplant rejection and identify potential diagnostic markers [109, 110]. The limitations of scRNA-seq due to fresh tissue requirements, artifactual transcriptional biases and loss of fragile or low viability cells can be overcome by snRNA-seq, enabling the storage of frozen samples and analysis using a quality comparable to scRNA-seq [108, 111,112,113]. The latest developments in snRNA-seq, imaging technologies [such as single-molecule fluorescence in situ hybridization (FISH)], proteomics (MIBI-TOF), and genomics all together provide benefits to further investigate cellular functions .
Classical immune cell characterization in different organs has limitations due to technologies such as microscopy and high-affinity antibody labeling. Additionally, conventional transcriptome studies may miss some essential immune cell subtypes . ScRNA-seq is currently widely used in immunology to unravel the differentiation process of HSCs, resolve previously under recognized immune cellular heterogeneity, decipher the immune cellular repertoire and predict disease-related phenotypes [39, 115, 116]. ScRNA-seq can help identify HSC fate branch points during differentiation; for example, conventional dendritic cells (cDCs) rely on the abundance of Siglec‑H and Ly6C to determine the direction of cDC type 1 (cDC1) or cDC type 2 (cDC2) . The myeloid progenitor cells that produce mast cells and eosinophils or monocytes and macrophages depend on the presence of GATA1 . Jaitin et al.  performed scRNA-seq on DCs and found different states of bone marrow-derived and other immune tissue-derived cells. Using scRNA-seq, progenitor immune cells can be identified by analyzing transcriptional variations in immature bone marrow and mature resident immune cells of specific organs . Many new types of immune cells were identified using scRNA-seq, which can analyze cell differentiation and cell lineages, including innate lymphocytes and lung interstitial macrophages . Recent studies have shown the heterogeneity of hematopoietic progenitor cells with a mixed lymphoid phenotype using scRNA-seq. This technology can also be used to identify novel cell types in several diseases, for example, cancer diagnosis and efficacy evaluation. By analyzing circulating tumor cells (CTCs) from prostate cancer patients, Miyamoto et al.  found that CTCs are highly heterogeneous in gene expression. The special feature of CTCs is that they activate Wnt signaling, which increases resistance to drug therapy. By analyzing brain cells in Alzheimer’s disease, scRNA-seq also enables the identification of microglia and perivascular macrophages related to neurodegenerative diseases . ScRNA-seq also improves the diagnosis of disease heterogeneity, such as identifying specific B-cell receptor signaling pathways and gene expression patterns in non-Hodgkin’s lymphoma. Most adult B-cell lymphomas exhibit a B-cell phenotype at the germinal center (GC). By combining scRNA-seq, the mixed characteristics of B cells derived from GC and follicular lymphoma (FL) revealed unique transcription characteristics [122, 123]. In the future, organ transplantation single-cell sequencing will likely help further elucidate the pathogenesis of transplant rejection and initiate the development of clinical trials and the emergence of more effective drugs to reduce the immune response associated with transplant rejection, thereby improving the quality of life of patients and extending patient survival.
Taken together, scRNA-seq can accurately interpret gene expression data, recognize cell heterogeneity, including new cell types or subtypes, and take snapshots of gene expression during the transition of cells from one state to another. All data can be integrated to define the critical process of cell development and differentiation, reveal the key signaling pathways and understand the gene regulatory network that predicts immune function [124,125,126,127].
Availability of data and materials
Bayesian Inference for Single-cell ClUstering and ImpuTing
Conventional dendritic cell
Circulating tumor cell
Cytometry by Time of Flight
Fluorescence Activated Cell Sorting
Hematopoietic cell transplantation
Human leukocyte antigen
Hematopoietic stem cells
Immune Cell Atlas
In vitro transcription
Remove unwanted variation
Single-cell RNA sequencing
Shared nearest neighbor
Single-nucleus RNA sequencing
T-cell receptor repertoire sequencing
Transcriptional immunophenotypical multipotent progenitors
T-distributed stochastic neighbor embedding
Uniform manifold approximation and projection
Unique molecular identifiers
Zero inflated factor analysis
Hao S, Yan KK, Ding L, Qian C, Chi H, Yu J. Network approaches for dissecting the immune system. iScience. 2020;23(8):101354.
Stubbington MJT, Rozenblatt-Rosen O, Regev A, Teichmann SA. Single-cell transcriptomics to explore the immune system in health and disease. Science. 2017;358(6359):58–63.
Ingulli E. Mechanism of cellular rejection in transplantation. Pediatr Nephrol. 2010;25(1):61–74.
Menon MC, Keung KL, Murphy B, O’connell PJ. The use of genomics and pathway analysis in our understanding and prediction of clinical renal transplant injury. Transplantation. 2016;100(7):1405–14.
Claeys E, Vermeire K. Immunosuppressive drugs in organ transplantation to prevent allograft rejection: mode of action and side effects. J Immunol Sci. 2019;3(4):14–21.
Nankivell BJ, Borrows RJ, Fung CLS, O’connell PJ, Allen RDM, Chapman JR. The natural history of chronic allograft nephropathy. N Engl J Med. 2003;349(24):2326–33.
El-Zoghby ZM, Stegall MD, Lager DJ, Kremers WK, Amer H, Gloor JM, et al. Identifying specific causes of kidney allograft loss. Am J Transplant. 2009;9(3):527–35.
Maguire O, Tario JD Jr, Shanahan TC, Wallace PK, Minderman H. Flow cytometry and solid organ transplantation: a perfect match. Immunol Invest. 2014;43(8):756–74.
Colvin RB. Antibody-mediated renal allograft rejection: diagnosis and pathogenesis. J Am Soc Nephrol. 2007;18(4):1046–56.
Malone AF, Wu H, Humphreys BD. Bringing renal biopsy interpretation into the molecular age with single-cell RNA sequencing. Semin Nephrol. 2018;38(1):31–9.
Sellares J, Reeve J, Loupy A, Mengel M, Sis B, Skene A, et al. Molecular diagnosis of antibody-mediated rejection in human kidney transplants. Am J Transplant. 2013;13(4):971–83.
Halloran PF, Pereira AB, Chang J, Matas A, Picton M, De Freitas D, et al. Microarray diagnosis of antibody-mediated rejection in kidney transplant biopsies: an international prospective study (INTERCOM). Am J Transplant. 2013;13(11):2865–74.
Thareja G, Yang H, Hayat S, Mueller FB, Lee JR, Lubetzky M, et al. Single nucleotide variant counts computed from RNA sequencing and cellular traffic into human kidney allografts. Am J Transplant. 2018;18(10):2429–42.
Roedder S, Sigdel T, Salomonis N, Hsieh S, Dai H, Bestard O, et al. The kSORT assay to detect renal transplant patients at high risk for acute rejection: results of the multicenter AART study. PLoS Med. 2014;11(11):e1001759.
Higdon LE, Schaffert S, Khatri P, Maltzman JS. Single cell immune profiling in transplantation research. Am J Transplant. 2019;19(5):1278–87.
Roy AL. Transcriptional regulation in the immune system: one cell at a time. Front Immunol. 2019;10:1355.
Giladi A, Amit I. Single-cell genomics: a stepping stone for future immunology discoveries. Cell. 2018;172(1–2):14–21.
Neu KE, Tang Q, Wilson PC, Khan AA. Single-cell genomics: approaches and utility in immunology. Trends Immunol. 2017;38(2):140–9.
Chen HI, Jin Y, Huang Y, Chen Y. Detection of high variability in gene expression from single-cell RNA-seq profiling. BMC Genom. 2016;17(Suppl 7):508.
Clatworthy MR. How to find a resident kidney macrophage: the single-cell sequencing solution. J Am Soc Nephrol. 2019;30(5):715–6.
Malone AF, Humphreys BD. Single-cell transcriptomics and solid organ transplantation. Transplantation. 2019;103(9):1776–82.
Macparland SA, Liu JC, Ma XZ, Innes BT, Bartczak AM, Gage BK, et al. Single cell RNA sequencing of human liver reveals distinct intrahepatic macrophage populations. Nat Commun. 2018;9(1):4383.
Wang W, Penland L, Gokce O, Croote D, Quake SR. High fidelity hypothermic preservation of primary tissues in organ transplant preservative for single cell transcriptome analysis. BMC Genom. 2018;19(1):140.
Guillaumet-Adkins A, Rodriguez-Esteban G, Mereu E, Mendez-Lago M, Jaitin DA, Villanueva A, et al. Single-cell transcriptome conservation in cryopreserved cells and tissues. Genome Biol. 2017;18(1):45.
Codeluppi S, Borm LE, Zeisel A, La Manno G, Van Lunteren JA, Svensson CI, et al. Spatial organization of the somatosensory cortex revealed by osmFISH. Nat Methods. 2018;15(11):932–5.
Ståhl PL, Salmén F, Vickovic S, Lundmark A, Navarro JF, Magnusson J, et al. Visualization and analysis of gene expression in tissue sections by spatial transcriptomics. Science. 2016;353:78–82.
Ziegenhain C, Vieth B, Parekh S, Reinius B, Guillaumet-Adkins A, Smets M, et al. Comparative analysis of single-cell RNA sequencing methods. Mol Cell. 2017;65(4):631-43.e4.
Lafzi A, Moutinho C, Picelli S, Heyn H. Tutorial: guidelines for the experimental design of single-cell RNA sequencing studies. Nat Protoc. 2018;13(12):2742–57.
Munoz-Manchado AB, Bengtsson Gonzales C, Zeisel A, Munguba H, Bekkouche B, Skene NG, et al. Diversity of interneurons in the dorsal striatum revealed by single-cell RNA sequencing and PatchSeq. Cell Rep. 2018;24(8):2179-90.e7.
Attaf N, Cervera-Marzal I, Dong C, Gil L, Renand A, Spinelli L, et al. FB5P-seq: FACS-based 5-prime end single-cell RNA-seq for integrative analysis of transcriptome and antigen receptor repertoire in B and T cells. Front Immunol. 2020;11:216.
Picelli S, Faridani OR, Björklund ÅK, Winberg G, Sagasser S, Sandberg R. Full-length RNA-seq from single cells using Smart-seq2. Nat Protoc. 2014;9(1):171–81.
Picelli S, Björklund ÅK, Faridani OR, Sagasser S, Winberg G, Sandberg R. Smart-seq2 for sensitive full-length transcriptome profiling in single cells. Nat Methods. 2013;10(11):1096–8.
Sasagawa Y, Nikaido I, Hayashi T, Danno H, Uno KD, Imai T, et al. Quartz-Seq: a highly reproducible and sensitive single-cell RNA sequencing method, reveals non-genetic gene-expression heterogeneity. Genome Biol. 2013;14(4):R31.
Sasagawa Y, Danno H, Takada H, Ebisawa M, Tanaka K, Hayashi T, et al. Quartz-Seq2: a high-throughput single-cell RNA-sequencing method that effectively uses limited sequence reads. Genome Biol. 2018;19(1):29.
Klein AM, Mazutis L, Akartuna I, Tallapragada N, Veres A, Li V, et al. Droplet barcoding for single-cell transcriptomics applied to embryonic stem cells. Cell. 2015;161(5):1187–201.
Hashimshony T, Senderovich N, Avital G, Klochendler A, De Leeuw Y, Anavy L, et al. CEL-Seq2: sensitive highly-multiplexed single-cell RNA-Seq. Genome Biol. 2016;17:77.
Hashimshony T, Wagner F, Sher N, Yanai I. CEL-Seq: Single-cell RNA-Seq by multiplexed linear amplification. Cell Rep. 2012;2(3):666–73.
Macosko EZ, Basu A, Satija R, Nemesh J, Shekhar K, Goldman M, et al. Highly parallel genome-wide expression profiling of individual cells using nanoliter droplets. Cell. 2015;161(5):1202–14.
Chen G, Ning B, Shi T. Single-cell RNA-seq technologies and related computational data analysis. Front Genet. 2019;10:317.
Kolodziejczyk AA, Kim JK, Svensson V, Marioni JC, Teichmann SA. The technology and biology of single-cell RNA sequencing. Mol Cell. 2015;58(4):610–20.
Gupta RK, Kuznicki J. Biological and medical importance of cellular heterogeneity deciphered by single-cell RNA sequencing. Cells. 2020;9(8):E1751.
Zou A, Ramanathan S, Dale RC, Brilot F. Single-cell approaches to investigate B cells and antibodies in autoimmune neurological disorders. Cell Mol Immunol. 2021;18(2):294–306.
Sheng K, Cao W, Niu Y, Deng Q, Zong C. Effective detection of variation in single-cell transcriptomes using MATQ-seq. Nat Methods. 2017;14(3):267–70.
Woodsworth DJ, Castellarin M, Holt RA. Sequence analysis of T-cell repertoires in health and disease. Genome Med. 2013;5(10):98.
Salomon R, Kaczorowski D, Valdes-Mora F, Nordon RE, Neild A, Farbehi N, et al. Droplet-based single cell RNAseq tools: a practical guide. Lab Chip. 2019;19(10):1706–27.
Ni J, Hu C, Li H, Li X, Fu Q, Czajkowsky DM, et al. Significant improvement in data quality with simplified SCRB-seq. Acta Biochim Biophys Sin (Shanghai). 2020;52(4):457–9.
Islam S, Kjallquist U, Moliner A, Zajac P, Fan JB, Lonnerberg P, et al. Characterization of the single-cell transcriptional landscape by highly multiplex RNA-seq. Genome Res. 2011;21(7):1160–7.
Keren-Shaul H, Kenigsberg E, Jaitin DA, David E, Paul F, Tanay A, et al. MARS-seq2.0: an experimental and analytical pipeline for indexed sorting combined with single-cell RNA sequencing. Nat Protoc. 2019;14(6):1841–62.
Hicks SC, Townes FW, Teng M, Irizarry RA. Missing data and technical variability in single-cell RNA-sequencing experiments. Biostatistics. 2018;19(4):562–78.
Leek JT, Scharpf RB, Bravo HC, Simcha D, Langmead B, Johnson WE, et al. Tackling the widespread and critical impact of batch effects in high-throughput data. Nat Rev Genet. 2010;11(10):733–9.
Risso D, Ngai J, Speed TP, Dudoit S. Normalization of RNA-seq data using factor analysis of control genes or samples. Nat Biotechnol. 2014;32(9):896–902.
Aran D, Looney AP, Liu L, Wu E, Fong V, Hsu A, et al. Reference-based analysis of lung single-cell sequencing reveals a transitional profibrotic macrophage. Nat Immunol. 2019;20(2):163–72.
Zhang X, Lan Y, Xu J, Quan F, Zhao E, Deng C, et al. Cell Marker: a manually curated resource of cell markers in human and mouse. Nucleic Acids Res. 2019;47(D1):D721–8.
Liao J, Yu Z, Chen Y, Bao M, Zou C, Zhang H, et al. Single-cell RNA sequencing of human kidney. Sci Data. 2020;7(1):4.
Kharchenko PV, Silberstein L, Scadden DT. Bayesian approach to single-cell differential expression analysis. Nat Methods. 2014;11(7):740–2.
Ilicic T, Kim JK, Kolodziejczyk AA, Bagger FO, Mccarthy DJ, Marioni JC, et al. Classification of low quality cells from single-cell RNA-seq data. Genome Biol. 2016;17:29.
Liu Y, Hu J, Liu D, Zhou S, Liao J, Liao G, et al. Single-cell analysis reveals immune landscape in kidneys of patients with chronic transplant rejection. Theranostics. 2020;10(19):8851–62.
Prabhakaran S, Azizi E, Carr A, Pe’er D. Dirichlet process mixture model for correcting technical variation in single-cell gene expression data. JMLR Workshop Conf Proc. 2016;48:1070–9.
Bacher R, Kendziorski C. Design and computational analysis of single-cell RNA-sequencing experiments. Genome Biol. 2016;17:63.
van der Maaten L, Hinton G. Visualizing data using t-SNE. J Mach Learn Res. 2008;9:2579–605.
Levine JH, Simonds EF, Bendall SC, Davis KL, El Amir AD, Tadmor MD, et al. Data-driven phenotypic dissection of AML reveals progenitor-like cells that correlate with prognosis. Cell. 2015;162(1):184–97.
Becht E, Mcinnes L, Healy J, Dutertre CA, Kwok IWH, Ng LG, et al. Dimensionality reduction for visualizing single-cell data using UMAP. Nat Biotechnol. 2018;37:38–44.
Xu C, Su Z. Identification of cell types from single-cell transcriptomes using a novel clustering method. Bioinformatics. 2015;31(12):1974–80.
Pierson E, Yau C. ZIFA: Dimensionality reduction for zero-inflated single-cell gene expression analysis. Genome Biol. 2015;16:241.
Korsunsky I, Millard N, Fan J, Slowikowski K, Zhang F, Wei K, et al. Fast, sensitive and accurate integration of single-cell data with Harmony. Nat Methods. 2019;16(12):1289–96.
Trapnell C, Cacchiarelli D, Grimsby J, Pokharel P, Li S, Morse M, et al. The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells. Nat Biotechnol. 2014;32(4):381–6.
Qiu X, Hill A, Packer J, Lin D, Ma YA, Trapnell C. Single-cell mRNA quantification and differential analysis with Census. Nat Methods. 2017;14(3):309–15.
Ji Z, Ji H. TSCAN: pseudo-time reconstruction and evaluation in single-cell RNA-seq analysis. Nucleic Acids Res. 2016;44(13):e117.
Setty M, Tadmor MD, Reich-Zeliger S, Angel O, Salame TM, Kathail P, et al. Wishbone identifies bifurcating developmental trajectories from single-cell data. Nat Biotechnol. 2016;34(6):637–45.
Street K, Risso D, Fletcher RB, Das D, Ngai J, Yosef N, et al. Slingshot: cell lineage and pseudotime inference for single-cell transcriptomics. BMC Genom. 2018;19(1):477.
Haghverdi L, Büttner M, Wolf FA, Buettner F, Theis FJ. Diffusion pseudotime robustly reconstructs lineage branching. Nat Methods. 2016;13(10):845–8.
Bendall SC, Davis KL, Amir EAD, Tadmor MD, Simonds EF, Chen TJ, et al. Single-cell trajectory detection uncovers progression and regulatory coordination in human B cell development. Cell. 2014;157(3):714–25.
Saelens W, Cannoodt R, Todorov H, Saeys Y. A comparison of single-cell trajectory inference methods. Nat Biotechnol. 2019;37(5):547–54.
Dangi A, Yu S, Luo X. Emerging approaches and technologies in transplantation: the potential game changers. Cell Mol Immunol. 2019;16(4):334–42.
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.
See P, Lum J, Chen J, Ginhoux F. A single-cell sequencing guide for immunologists. Front Immunol. 2018;9:2425.
de Simone M, Rossetti G, Pagani M. Single cell T cell receptor sequencing: techniques and future challenges. Front Immunol. 2018;9:1638.
Zheng P, Tamaresis J, Thangavelu G, Xu L, You X, Blazar BR, et al. Recipient-specific T-cell repertoire reconstitution in the gut following murine hematopoietic cell transplant. Blood Adv. 2020;4(17):4232–43.
Stubbington MJT, Lönnberg T, Proserpio V, Clare S, Speak AO, Dougan G, et al. T cell fate and clonality inference from single-cell transcriptomes. Nat Methods. 2016;13(4):329–32.
Morris H, Dewolf S, Robins H, Sprangers B, Locascio SA, Shonts BA, et al. Tracking donor-reactive T cells: evidence for clonal deletion in tolerant kidney transplant patients. Sci Transl Med. 2015;7(272):272ra10.
DeWolf S, Sykes M. Alloimmune T cells in transplantation. J Clin Invest. 2017;127(7):2473–81.
Alachkar H, Mutonga M, Kato T, Kalluri S, Kakuta Y, Uemura M, et al. Quantitative characterization of T-cell repertoire and biomarkers in kidney transplant rejection. BMC Nephrol. 2016;17(1):181.
Yew PY, Alachkar H, Yamaguchi R, Kiyotani K, Fang H, Yap KL, et al. Quantitative characterization of T-cell repertoire in allogeneic hematopoietic stem cell transplant recipients. Bone Marrow Transplant. 2015;50(9):1227–34.
Halloran PF, Pereira AB, Chang J, Matas A, Picton M, De Freitas D, et al. Potential impact of microarray diagnosis of T cell-mediated rejection in kidney transplants: the INTERCOM study. Am J Transplant. 2013;13(9):2352–63.
Stewart BJ, Ferdinand JR, Clatworthy MR. Using single-cell technologies to map the human immune system: implications for nephrology. Nat Rev Nephrol. 2020;16(2):112–28.
Reinders ME, Rabelink TJ, Briscoe DM. Angiogenesis and endothelial cell repair in renal disease and allograft rejection. J Am Soc Nephrol. 2006;17(4):932–42.
Sis B, Jhangri GS, Bunnag S, Allanach K, Kaplan B, Halloran PF. Endothelial gene expression in kidney transplants with alloantibody indicates antibody-mediated damage despite lack of C4d staining. Am J Transplant. 2009;9(10):2312–23.
Malone AF, Wu H, Fronick C, Fulton R, Gaut JP, Humphreys BD. Harnessing expressed single nucleotide variation and single cell rna sequencing to define immune cell chimerism in the rejecting kidney transplant. J Am Soc Nephrol. 2020;31(9):1977–86.
Bellan C, Amato T, Carmellini M, Onorati M, D’amuri A, Leoncini L, et al. Analysis of the IgVH genes in T cell-mediated and antibody-mediated rejection of the kidney graft. J Clin Pathol. 2011;64(1):47–53.
Dangi A, Natesh NR, Husain I, Ji Z, Barisoni L, Kwun J, et al. Single cell transcriptomics of mouse kidney transplants reveals a myeloid cell pathway for transplant rejection. JCI Insight. 2020;5(20):e141321.
Varma E, Luo X, Muthukumar T. Dissecting the human kidney allograft transcriptome: single-cell RNA sequencing. Curr Opin Organ Transplant. 2021;26(1):43–51.
Lei H, Reinke P, Volk HD, Lv Y, Wu R. Mechanisms of immune tolerance in liver transplantation-crosstalk between alloreactive T cells and liver cells with therapeutic prospects. Front Immunol. 2019;10:2667.
Dai H, Zheng Y, Thomson AW, Rogers NM. Transplant tolerance induction: insights from the liver. Front Immunol. 2020;11:1044.
Feng S, Bucuvalas J. Tolerance after liver transplantation: Where are we? Liver Transpl. 2017;23(12):1601–14.
Thomson AW, Vionnet J, Sanchez-Fueyo A. Understanding, predicting and achieving liver transplant tolerance: from bench to bedside. Nat Rev Gastroenterol Hepatol. 2020;17(12):719–39.
Choudhary NS, Saigal S, Bansal RK, Saraf N, Gautam D, Soin AS. Acute and chronic rejection after liver transplantation: what a clinician needs to know. J Clin Exp Hepatol. 2017;7(4):358–66.
Feng S, Bucuvalas JC, Demetris AJ, Burrell BE, Spain KM, Kanaparthi S, et al. Evidence of chronic allograft injury in liver biopsies from long-term pediatric recipients of liver transplants. Gastroenterology. 2018;155(6):1838-51.e7.
Misharin AV, Morales-Nebreda L, Reyfman PA, Cuda CM, Walter JM, Mcquattie-Pimentel AC, et al. Monocyte-derived alveolar macrophages drive lung fibrosis and persist in the lung over the life span. J Exp Med. 2017;214(8):2387–404.
Mccubbrey AL, Barthel L, Mohning MP, Redente EF, Mould KJ, Thomas SM, et al. Deletion of c-FLIP from CD11bhi macrophages prevents development of bleomycin-induced lung fibrosis. Am J Respir Cell Mol Biol. 2018;58(1):66–78.
Reyfman PA, Walter JM, Joshi N, Anekalla KR, Mcquattie-Pimentel AC, Chiu S, et al. Single-cell transcriptomic analysis of human lung provides insights into the pathobiology of pulmonary fibrosis. Am J Respir Crit Care Med. 2018;199(12):1517–36.
Mould KJ, Moore CM, Mcmanus SA, Mccubbrey AL, Mcclendon JD, Griesmer CL, et al. Airspace macrophages and monocytes exist in transcriptionally distinct subsets in healthy adults. Am J Respir Crit Care Med. 2021;203(8):946–56.
Snyder ME, Finlayson MO, Connors TJ, Dogra P, Senda T, Bush E, et al. Generation and persistence of human tissue-resident memory T cells in lung transplantation. Sci Immunol. 2019;4(33):eaav5581.
Dong F, Hao S, Zhang S, Zhu C, Cheng H, Yang Z, et al. Differentiation of transplanted haematopoietic stem cells tracked by single-cell transcriptomic analysis. Nat Cell Biol. 2020;22(6):630–9.
Fujisaki J, Wu J, Carlson AL, Silberstein L, Putheti P, Larocca R, et al. In vivo imaging of Treg cells providing immune privilege to the haematopoietic stem-cell niche. Nature. 2011;474(7350):216–9.
Beilhack A, Schulz S, Baker J, Beilhack GF, Wieland CB, Herman EI, et al. In vivo analyses of early events in acute graft-versus-host disease reveal sequential infiltration of T-cell subsets. Blood. 2005;106(3):1113–22.
DiRienzo CG, Murphy GF, Jones SC, Korngold R, Friedman TM. T-cell receptor Valpha spectratype analysis of a CD4-mediated T-cell response against minor histocompatibility antigens involved in severe graft-versus-host disease. Biol Blood Marrow Transplant. 2006;12(8):818–27.
Divito SJ, Aasebo AT, Matos TR, Hsieh PC, Collin M, Elco CP, et al. Peripheral host T cells survive hematopoietic stem cell transplantation and promote graft-versus-host disease. J Clin Invest. 2020;130(9):4624–36.
Habib N, Avraham-Davidi I, Basu A, Burks T, Shekhar K, Hofree M, et al. Massively parallel single-nucleus RNA-seq with DroNc-seq. Nat Methods. 2017;14(10):955–8.
Babel N, Stervbo U, Reinke P, Volk HD. The identity card of T cells-clinical utility of T-cell receptor repertoire analysis in transplantation. Transplantation. 2019;103(8):1544–55.
Chen H, Ye F, Guo G. Revolutionizing immunology with single-cell RNA sequencing. Cell Mol Immunol. 2019;16(3):242–9.
Denisenko E, Guo BB, Jones M, Hou R, De Kock L, Lassmann T, et al. Systematic assessment of tissue dissociation and storage biases in single-cell and single-nucleus RNA-seq workflows. Genome Biol. 2020;21(1):130.
Wu H, Kirita Y, Donnelly EL, Humphreys BD. Advantages of single-nucleus over single-cell RNA sequencing of adult kidney: rare cell types and novel cell states revealed in fibrosis. J Am Soc Nephrol. 2019;30(1):23–32.
Slyper M, Porter CBM, Ashenberg O, Waldman J, Drokhlyansky E, Wakiro I, et al. A single-cell and single-nucleus RNA-Seq toolbox for fresh and frozen human tumors. Nat Med. 2020;26(5):792–802.
Roy AL, Conroy RS. Toward mapping the human body at a cellular resolution. Mol Biol Cell. 2018;29(15):1779–85.
Vegh P, Haniffa M. The impact of single-cell RNA sequencing on understanding the functional organization of the immune system. Brief Funct Genom. 2018;17(4):265–72.
Zeng T, Dai H. Single-cell rna sequencing-based computational analysis to describe disease heterogeneity. Front Genet. 2019;10:629.
Schlitzer A, Sivakamasundari V, Chen J, Sumatoh HR, Schreuder J, Lum J, et al. Identification of cDC1- and cDC2-committed DC progenitors reveals early lineage priming at the common DC progenitor stage in the bone marrow. Nat Immunol. 2015;16(7):718–28.
Pang WW, Price EA, Sahoo D, Beerman I, Maloney WJ, Rossi DJ, et al. Human bone marrow hematopoietic stem cells are increased in frequency and myeloid-biased with age. Proc Natl Acad Sci USA. 2011;108(50):20012–7.
Jaitin DA, Kenigsberg E, Keren-Shaul H, Elefant N, Paul F, Zaretsky I, et al. Massively parallel single-cell RNA-seq for marker-free decomposition of tissues into cell types. Science. 2014;343(6172):776–9.
Miyamoto DT, Zheng Y, Wittner BS, Lee RJ, Zhu H, Broderick KT, et al. RNA-Seq of single prostate CTCs implicates noncanonical Wnt signaling in antiandrogen resistance. Science. 2015;349(6254):1351–6.
Keren-Shaul H, Spinrad A, Weiner A, Matcovitch-Natan O, Dvir-Szternfeld R, Ulland TK, et al. A unique microglia type associated with restricting development of Alzheimer’s disease. Cell. 2017;169(7):1276-90.e17.
Myklebust JH, Brody J, Kohrt HE, Kolstad A, Czerwinski DK, Walchli S, et al. Distinct patterns of B-cell receptor signaling in non-Hodgkin lymphomas identified by single-cell profiling. Blood. 2017;129(6):759–70.
Young RM, Staudt LM. Targeting pathological B cell receptor signalling in lymphoid malignancies. Nat Rev Drug Discov. 2013;12(3):229–43.
Zheng GX, Terry JM, Belgrader P, Ryvkin P, Bent ZW, Wilson R, et al. Massively parallel digital transcriptional profiling of single cells. Nat Commun. 2017;8:14049.
Bjorklund AK, Forkel M, Picelli S, Konya V, Theorell J, Friberg D, et al. The heterogeneity of human CD127(+) innate lymphoid cells revealed by single-cell RNA sequencing. Nat Immunol. 2016;17(4):451–60.
Paul F, Arkin Y, Giladi A, Jaitin DA, Kenigsberg E, Keren-Shaul H, et al. Transcriptional heterogeneity and lineage commitment in myeloid progenitors. Cell. 2016;164(1–2):325.
Ye Y, Song H, Zhang J, Shi S. Understanding the biology and pathogenesis of the kidney by single-cell transcriptomic analysis. Kidney Dis (Basel). 2018;4(4):214–25.
The authors thank Seq Health Technology Co., Ltd. (Wuhan, China) for providing guidance on sequencing technologies.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Wang, Y., Wang, JY., Schnieke, A. et al. Advances in single-cell sequencing: insights from organ transplantation. Military Med Res 8, 45 (2021). https://doi.org/10.1186/s40779-021-00336-1
- Single-cell RNA sequencing
- Transplant rejection
- Immune cell interactions
- Transcriptional profiling
- 10× Genomics chromium