Dandelion utilizes single cell adaptive immune receptor repertoire to explore lymphocyte developmental origins

Assessment of single-cell gene expression (scRNA-seq) and adaptive immune receptor sequencing (scVDJ-seq) has been invaluable in studying lymphocyte biology. Here, we introduce Dandelion, a computational pipeline for scVDJ-seq analysis. It enables the application of standard V(D)J analysis workflows to single-cell datasets, delivering improved V(D)J contig annotation and the identification of non-productive and partially spliced contigs. We devised a novel strategy to create an adaptive immune receptor feature space that can be used for both differential V(D)J usage analysis and pseudotime trajectory inference. The application of Dandelion improved the alignment of human thymic development trajectories of double positive T cells to mature single-positive CD4/CD8 T cells, with important new predictions of factors regulating lineage commitment. Dandelion analysis of other cell compartments provided novel insights into the origins of human B1 cells and ILC/NK cell development, illustrating the power of our approach. Dandelion is an open access resource (https://www.github.com/zktuong/dandelion) that will enable future discoveries.

[1]  S. Teichmann,et al.  Integrating T-cell receptor and transcriptome for large-scale single-cell immune profiling analysis , 2021, bioRxiv.

[2]  D. Nemazee,et al.  Asymmetrical forward and reverse developmental trajectories determine molecular programs of B cell antigen receptor editing , 2022, Science Immunology.

[3]  J. Park,et al.  Cross-tissue immune cell analysis reveals tissue-specific features in humans , 2022, Science.

[4]  D. Jaffe,et al.  Functional antibodies exhibit light chain coherence , 2022, bioRxiv.

[5]  Daniel B. Burkhardt,et al.  Multimodal single cell data integration challenge: results and lessons learned , 2022, bioRxiv.

[6]  S. Teichmann,et al.  Mapping the developing human immune system across organs , 2022, bioRxiv.

[7]  S. Teichmann,et al.  Differential abundance testing on single-cell data using k-nearest neighbor graphs , 2021, Nature Biotechnology.

[8]  Michael J. T. Stubbington,et al.  Integrating T cell receptor sequences and transcriptional profiles by clonotype neighbor graph analysis (CoNGA) , 2021, Nature Biotechnology.

[9]  X. Liu,et al.  TRUST4: immune repertoire reconstruction from bulk and single-cell RNA-seq data , 2021, Nature Methods.

[10]  K. McNagny,et al.  ILC-You in the Thymus: A Fresh Look at Innate Lymphoid Cell Development , 2021, Frontiers in Immunology.

[11]  Frances E. Muldoon,et al.  Single-cell multi-omics analysis of the immune response in COVID-19 , 2021, Nature Medicine.

[12]  Mohammad M. Karimi,et al.  The order and logic of CD4 versus CD8 lineage choice and differentiation in mouse thymus , 2021, Nature Communications.

[13]  T. Mora,et al.  Population variability in the generation and selection of T-cell repertoires , 2020, PLoS Comput. Biol..

[14]  Cédric R. Weber,et al.  Platypus: an open-access software for integrating lymphocyte single-cell immune repertoires with transcriptomes , 2020, bioRxiv.

[15]  M. Messing,et al.  Abortive γδTCR rearrangements suggest ILC2s are derived from T-cell precursors. , 2020, Blood advances.

[16]  H. Hosokawa,et al.  How transcription factors drive choice of the T cell fate , 2020, Nature reviews. Immunology.

[17]  Do Young Hyeon,et al.  Single-cell multiomics: technologies and data analysis methods , 2020, Experimental & Molecular Medicine.

[18]  Zlatko Trajanoski,et al.  Scirpy: A Scanpy extension for analyzing single-cell T-cell receptor sequencing data , 2020, bioRxiv.

[19]  S. Teichmann,et al.  Immunology in the Era of Single-Cell Technologies. , 2020, Annual review of immunology.

[20]  S. Teichmann,et al.  A cell atlas of human thymic development defines T cell repertoire formation , 2020, Science.

[21]  S. Chatterjee.  A New Coefficient of Correlation , 2019, Journal of the American Statistical Association.

[22]  Nicholas L. Bormann,et al.  scRepertoire: An R-based toolkit for single-cell immune receptor analysis. , 2020, F1000Research.

[23]  David McDonald,et al.  Decoding human fetal liver haematopoiesis , 2019, Nature.

[24]  Fabian J Theis,et al.  Predicting antigen specificity of single T cells based on TCR CDR3 regions , 2019, bioRxiv.

[25]  B. Povinelli,et al.  Discovery of a CD10 negative B-progenitor in human fetal life identifies unique ontogeny-related developmental programs. , 2019, Blood.

[26]  A. Krug,et al.  What Makes a pDC: Recent Advances in Understanding Plasmacytoid DC Development and Heterogeneity , 2019, Front. Immunol..

[27]  Yvan Saeys,et al.  A comparison of single-cell trajectory inference methods , 2019, Nature Biotechnology.

[28]  J. Wren,et al.  Suppression of ILC2 differentiation from committed T cell precursors by E protein transcription factors , 2019, The Journal of experimental medicine.

[29]  Andrew J. Hill,et al.  The single cell transcriptional landscape of mammalian organogenesis , 2019, Nature.

[30]  Kong-Peng Lam,et al.  BCR-dependent lineage plasticity in mature B cells , 2019, Science.

[31]  D. Pe’er,et al.  Characterization of cell fate probabilities in single-cell data with Palantir , 2019, Nature Biotechnology.

[32]  M. Farrar,et al.  B-1a cells acquire their unique characteristics by bypassing the pre-BCR selection stage , 2017, Nature Communications.

[33]  Michael I. Jordan,et al.  Deep Generative Modeling for Single-cell Transcriptomics , 2018, Nature Methods.

[34]  Syed Ahmad Chan Bukhari,et al.  AIRR Community Standardized Representations for Annotated Immune Repertoires , 2018, Front. Immunol..

[35]  Thierry Mora,et al.  Detecting T-cell receptors involved in immune responses from single repertoire snapshots , 2018, bioRxiv.

[36]  Samuel L. Wolock,et al.  Scrublet: computational identification of cell doublets in single-cell transcriptomic data , 2018, bioRxiv.

[37]  Paul Hoffman,et al.  Integrating single-cell transcriptomic data across different conditions, technologies, and species , 2018, Nature Biotechnology.

[38]  M. Busslinger,et al.  Control of B-1a cell development by instructive BCR signaling. , 2018, Current opinion in immunology.

[39]  T. Hughes,et al.  The Human Transcription Factors , 2018, Cell.

[40]  R. Satija,et al.  Single-cell RNA sequencing to explore immune cell heterogeneity , 2017, Nature Reviews Immunology.

[41]  Fabian J Theis,et al.  SCANPY: large-scale single-cell gene expression data analysis , 2018, Genome Biology.

[42]  Syed Ahmad Chan Bukhari,et al.  Adaptive Immune Receptor Repertoire Community recommendations for sharing immune-repertoire sequencing data , 2017, Nature Immunology.

[43]  Syed Ahmad Chan Bukhari,et al.  Reproducibility and Reuse of Adaptive Immune Receptor Repertoire Data , 2017, Front. Immunol..

[44]  M. Krangel,et al.  Tcrd Rearrangement Redirects a Processive Tcra Recombination Program to Expand the Tcra Repertoire. , 2017, Cell reports.

[45]  Fabian J Theis,et al.  Diffusion pseudotime robustly reconstructs lineage branching , 2016, Nature Methods.

[46]  Steven H. Kleinstein,et al.  Change-O: a toolkit for analyzing large-scale B cell immunoglobulin repertoire sequencing data , 2015, Bioinform..

[47]  T. Mora,et al.  Inferring processes underlying B-cell repertoire diversity , 2015, bioRxiv.

[48]  G. Yaari,et al.  Automated analysis of high-throughput B-cell sequencing data reveals a high frequency of novel immunoglobulin V gene segment alleles , 2015, Proceedings of the National Academy of Sciences.

[49]  Tiago P. Peixoto,et al.  The graph-tool python library , 2014 .

[50]  Malay Mandal,et al.  Orchestrating B cell lymphopoiesis through interplay of IL-7 receptor and pre-B cell receptor signalling , 2013, Nature Reviews Immunology.

[51]  D. Kerby.  The Simple Difference Formula: An Approach to Teaching Nonparametric Correlation1: , 2014 .

[52]  Ning Ma,et al.  IgBLAST: an immunoglobulin variable domain sequence analysis tool , 2013, Nucleic Acids Res..

[53]  Li Wu,et al.  Convergent differentiation: myeloid and lymphoid pathways to murine plasmacytoid dendritic cells. , 2013, Blood.

[54]  E. Montecino-Rodriguez,et al.  B-1 B cell development in the fetus and adult. , 2012, Immunity.

[55]  M. Yoder,et al.  Embryonic day 9 yolk sac and intra-embryonic hemogenic endothelium independently generate a B-1 and marginal zone progenitor lacking B-2 potential , 2011, Proceedings of the National Academy of Sciences.

[56]  Eugene V Koonin,et al.  Complex selection on 5' splice sites in intron-rich organisms. , 2009, Genome research.

[57]  Qingzhao Zhang,et al.  A differentiation pathway for B1 cells in adult bone marrow , 2009, Proceedings of the National Academy of Sciences.

[58]  A. Singer,et al.  Lineage fate and intense debate: myths, models and mechanisms of CD4- versus CD8-lineage choice , 2008, Nature Reviews Immunology.

[59]  J. Kaye,et al.  Development of all CD4 T lineages requires nuclear factor TOX. , 2008, The Journal of experimental medicine.

[60]  E. Montecino-Rodriguez,et al.  Identification of a B-1 B cell–specified progenitor , 2006, Nature Immunology.

[61]  Yifan Hu,et al.  Efficient, High-Quality Force-Directed Graph Drawing , 2006 .

[62]  M. Busslinger,et al.  Derivation of 2 categories of plasmacytoid dendritic cells in murine bone marrow. , 2005, Blood.

[63]  Xiaolong Liu,et al.  The zinc finger protein cKrox directs CD4 lineage differentiation during intrathymic T cell positive selection , 2005, Nature Immunology.

[64]  S. Habu,et al.  Dual functions of Runx proteins for reactivating CD8 and silencing CD4 at the commitment process into CD8 thymocytes. , 2005, Immunity.

[65]  Yi Zhang,et al.  The zinc finger transcription factor Th-POK regulates CD4 versus CD8 T-cell lineage commitment , 2005, Nature.

[66]  P. Leder,et al.  Plasmacytoid dendritic cells activate lymphoid-specific genetic programs irrespective of their cellular origin. , 2004, Immunity.

[67]  S. Pai,et al.  Critical roles for transcription factor GATA-3 in thymocyte development. , 2003, Immunity.

[68]  J. Alberola-Ila,et al.  GATA-3 expression is controlled by TCR signals and regulates CD4/CD8 differentiation. , 2003, Immunity.

[69]  Li Wu,et al.  The Lymphoid Past of Mouse Plasmacytoid Cells and Thymic Dendritic Cells1 , 2003, The Journal of Immunology.

[70]  D. Littman,et al.  Differential Requirements for Runx Proteins in CD4 Repression and Epigenetic Silencing during T Lymphocyte Development , 2002, Cell.

[71]  H. Le Hir,et al.  The exon–exon junction complex provides a binding platform for factors involved in mRNA export and nonsense‐mediated mRNA decay , 2001, The EMBO journal.

[72]  Michael J. Bevan,et al.  Itm2a Is Induced during Thymocyte Selection and T Cell Activation and Causes Downregulation of Cd8 When Overexpressed in Cd4+Cd8+ Double Positive Thymocytes , 1999, The Journal of experimental medicine.

[73]  F. Alt,et al.  Assembly of Productive T Cell Receptor (cid:100) Variable Region Genes Exhibits Allelic Inclusion , 1998 .

[74]  Marie-Paule Lefranc,et al.  IMGT, the international ImMunoGeneTics database. , 1997, Nucleic acids research.

[75]  C. Janeway.  Immunobiology: The Immune System in Health and Disease , 1996 .

[76]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[77]  K. Rajewsky,et al.  A critical role of λ5 protein in B cell development , 1992, Cell.

[78]  J. Kearney,et al.  An embryonic source of Ly1 but not conventional B cells. , 1991, International immunology.

[79]  L. Herzenberg,et al.  Toward a layered immune system , 1989, Cell.

[80]  M. Woodruff,et al.  Immune Response , 1969, Nature.