High-throughput targeted long-read single cell sequencing reveals the clonal and transcriptional landscape of lymphocytes

Nature Communications

Volumn 10, Issue 1, 2019, Pages

  Singh, Mandeep ^a,b   Al Eryani, Ghamdan ^a,b   Carswell, Shaun ^a   Ferguson, James M ^a   Blackburn, James ^a,b   Barton, Kirston ^a,b   Roden, Daniel ^a,b   Luciani, Fabio ^b,c   Giang Phan, Tri ^a,b   Junankar, Simon ^a,b   Jackson, Katherine ^a,b   Goodnow, Christopher C ^a,b   Smith, Martin A ^a,b   Swarbrick, Alexander ^a,b  

^a GARVAN INSTITUTE OF MEDICAL RESEARCH   (Australia)

^b ST VINCENT S HOSPITAL   (Australia)

^c UNIVERSITY OF NEW SOUTH WALES   (Australia)

Author keywords

[No Author keywords available]

Indexed keywords

B LYMPHOCYTE RECEPTOR; COMPLEMENTARY DNA; LYMPHOCYTE ANTIGEN RECEPTOR; T LYMPHOCYTE RECEPTOR; TRANSCRIPTOME;

ANTIGEN; CANCER; CELL; GENE EXPRESSION; GENETIC ANALYSIS; IMMUNE SYSTEM; RNA;

ARTICLE; AUTOIMMUNITY; CELL SELECTION; CELL SUSPENSION; CLONAL EVOLUTION; CLONING; DNA LIBRARY; EXON; FLOW CYTOMETRY; GENE EXPRESSION; GENE SEQUENCE; GENETIC DISTANCE; GERMLINE MUTATION; HIGH THROUGHPUT SEQUENCING; HUMAN; HUMAN CELL; HUMAN TISSUE; LYMPH NODE; LYMPHOCYTE; POINT MUTATION; POLYMERASE CHAIN REACTION; PRIMARY TUMOR; QUALITY CONTROL; REGULATORY SEQUENCE; RNA SEQUENCE; SEQUENCE ALIGNMENT; SEQUENCE ANALYSIS; SINGLE CELL ANALYSIS; SOMATIC MUTATION; TRANSCRIPTION REGULATION; TRIPLE NEGATIVE BREAST CANCER; GENE EXPRESSION PROFILING; GENETICS; IMMUNOLOGY; METABOLISM; PRIMARY CELL CULTURE; PROCEDURES;

CLONAL EVOLUTION; GENE EXPRESSION PROFILING; HIGH-THROUGHPUT NUCLEOTIDE SEQUENCING; HUMANS; LYMPHOCYTES; PRIMARY CELL CULTURE; RECEPTORS, ANTIGEN, B-CELL; RECEPTORS, ANTIGEN, T-CELL; SEQUENCE ANALYSIS, RNA; SINGLE-CELL ANALYSIS;

EID: 85069532270     PISSN: None     EISSN: 20411723     Source Type: Journal    
DOI: 10.1038/s41467-019-11049-4     Document Type: Article

References (67)

1. Market, E. & Papavasiliou, F. N. V(D)J recombination and the evolution of the adaptive immune system. PLoS Biol. 1, E16, 10.1371/journal.pbio.0000016 (2003).
2. Wang, E. T. et al. Alternative isoform regulation in human tissue transcriptomes. Nature 456, 470–476, 10.1038/nature07509 (2008).
3. Bassing, C. H., Swat, W. & Alt, F. W. The mechanism and regulation of chromosomal V(D)J recombination. Cell 109(Suppl), S45–55 (2002).
4. Chaudhuri, J. & Alt, F. W. Class-switch recombination: interplay of transcription, DNA deamination and DNA repair. Nat. Rev. Immunol. 4, 541–552, 10.1038/nri1395 (2004).
5. Alt, F. W. et al. Synthesis of secreted and membrane-bound immunoglobulin mu heavy chains is directed by mRNAs that differ at their 3’ ends. Cell 20, 293–301 (1980).
6. Willis, N. A., Rass, E. & Scully, R. Deciphering the code of the cancer. Genome: mechanisms of chromosome rearrangement. Trends Cancer 1, 217–230, 10.1016/j.trecan.2015.10.007 (2015).
7. Calis, J. J. & Rosenberg, B. R. Characterizing immune repertoires by high throughput sequencing: strategies and applications. Trends Immunol. 35, 581–590, 10.1016/j.it.2014.09.004 (2014).
8. Laydon, D. J., Bangham, C. R. & Asquith, B. Estimating T-cell repertoire diversity: limitations of classical estimators and a new approach. Philos. Trans. R. Soc. Lond. Ser. B 10.1098/rstb.2014.0291 (2015).
9. Papalexi, E., Satija, R. & Single-cell, R. N. A. sequencing to explore immune cell heterogeneity. Nat. Rev. Immunol. 18, 35–45, 10.1038/nri.2017.76 (2018).
10. Stubbington, M. J. T., Rozenblatt-Rosen, O., Regev, A. & Teichmann, S. A. Single-cell transcriptomics to explore the immune system in health and disease. Science 358, 58–63, 10.1126/science.aan6828 (2017).
11. Picelli, S. et al. Smart-seq2 for sensitive full-length transcriptome profiling in single cells. Nat. Methods 10, 1096–1098, 10.1038/nmeth.2639 (2013).
12. Afik, S. et al. Targeted reconstruction of T cell receptor sequence from single cell RNA-seq links CDR3 length to T cell differentiation state. Nucleic Acids Res. 45, e148, 10.1093/nar/gkx615 (2017).
13. Eltahla, A. A. et al. Linking the T cell receptor to the single cell transcriptome in antigen-specific human T cells. Immunol. Cell Biol. 94, 604–611, 10.1038/icb.2016.16 (2016).
14. Rizzetto, S. et al. B-cell receptor reconstruction from single-cell RNA-seq with VDJPuzzle. Bioinformatics 10.1093/bioinformatics/bty203 (2018).
15. Stubbington, M. J. T. et al. T cell fate and clonality inference from single-cell transcriptomes. Nat. Methods 13, 329–332, 10.1038/nmeth.3800 (2016).
16. Upadhyay, A. A. et al. BALDR: a computational pipeline for paired heavy and light chain immunoglobulin reconstruction in single-cell RNA-seq data. Genome Med. 10, 20, 10.1186/s13073-018-0528-3 (2018).
17. Rizzetto, S. et al. Impact of sequencing depth and read length on single cell RNA sequencing data of T. Cells Sci. Rep. 7, 12781, 10.1038/s41598-017-12989-x (2017).
18. Ziegenhain, C. et al. Comparative analysis of single-cell RNA sequencing methods. Mol. Cell 65, 631–643.e634, 10.1016/j.molcel.2017.01.023 (2017).
19. Cao, J. et al. Comprehensive single-cell transcriptional profiling of a multicellular organism. Science 357, 661–667, 10.1126/science.aam8940 (2017).
20. Fan, H. C., Fu, G. K. & Fodor, S. P. Expression profiling. Combinational labeling of single cells for gene expression cytometry. Science 347, 1258367, 10.1126/science.1258367 (2015).
21. Gierahn, T. M. et al. Seq-Well: portable, low-cost RNA sequencing of single cells at high throughput. Nat. Methods 14, 395–398, 10.1038/nmeth.4179 (2017).
22. Han, X. et al. Mapping the mouse cell atlas by microwell-seq. Cell 173, 1307, 10.1016/j.cell.2018.05.012 (2018).
23. Klein, A. M. et al. Droplet barcoding for single-cell transcriptomics applied to embryonic stem cells. Cell 161, 1187–1201, 10.1016/j.cell.2015.04.044 (2015).
24. Macosko, E. Z. et al. Highly parallel genome-wide expression profiling of individual cells using nanoliter droplets. Cell 161, 1202–1214, 10.1016/j.cell.2015.05.002 (2015).
25. Rosenberg, A. B. et al. Single-cell profiling of the developing mouse brain and spinal cord with split-pool barcoding. Science 360, 176–182, 10.1126/science.aam8999 (2018).
26. Zheng, G. X. et al. Massively parallel digital transcriptional profiling of single cells. Nat. Commun. 8, 14049, 10.1038/ncomms14049 (2017).
27. Azizi, E. et al. Single-cell map of diverse immune phenotypes in the breast tumor microenvironment. Cell, 10.1016/j.cell.2018.05.060 (2018).
28. Byrne, A. et al. Nanopore long-read RNAseq reveals widespread transcriptional variation among the surface receptors of individual B cells. Nat. Commun. 8, 16027, 10.1038/ncomms16027 (2017).
29. Carlson, C. S. et al. Using synthetic templates to design an unbiased multiplex PCR assay. Nat. Commun. 4, 2680, 10.1038/ncomms3680 (2013).
30. Shugay, M. et al. Towards error-free profiling of immune repertoires. Nat. Methods 11, 653–655, 10.1038/nmeth.2960 (2014).
31. Chapman, C. J., Zhou, J. X., Gregory, C., Rickinson, A. B. & Stevenson, F. K. VH and VL gene analysis in sporadic Burkitt’s lymphoma shows somatic hypermutation, intraclonal heterogeneity, and a role for antigen selection. Blood 88, 3562–3568 (1996).
32. Yanagi, Y., Chan, A., Chin, B., Minden, M. & Mak, T. W. Analysis of cDNA clones specific for human T cells and the alpha and beta chains of the T-cell receptor heterodimer from a human T-cell line. Proc. Natl Acad. Sci. USA 82, 3430–3434 (1985).
33. Koren, S. et al. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 27, 722–736, 10.1101/gr.215087.116 (2017).
34. Vaser, R., Sovic, I., Nagarajan, N. & Sikic, M. Fast and accurate de novo genome assembly from long uncorrected reads. Genome Res. 27, 737–746, 10.1101/gr.214270.116 (2017).
35. Loman, N. J., Quick, J. & Simpson, J. T. A complete bacterial genome assembled de novo using only nanopore sequencing data. Nat. Methods 12, 733–735, 10.1038/nmeth.3444 (2015).
36. Sale, J. E. & Neuberger, M. S. TdT-accessible breaks are scattered over the immunoglobulin V domain in a constitutively hypermutating B cell line. Immunity 9, 859–869 (1998).
37. Maki, R. et al. The role of DNA rearrangement and alternative RNA processing in the expression of immunoglobulin delta genes. Cell 24, 353–365 (1981).
38. Yuan, D. & Tucker, P. W. Regulation of IgM and IgD synthesis in B lymphocytes. I. Changes in biosynthesis of mRNA for mu- and delta-chains. J. Immunol. 132, 1561–1565 (1984).
39. Roskin, K. M. et al. IgH sequences in common variable immune deficiency reveal altered B cell development and selection. Sci. Transl. Med. 7, 302ra135, 10.1126/scitranslmed.aab1216 (2015).
40. Lamson, G. & Koshland, M. E. Changes in J chain and mu chain RNA expression as a function of B cell differentiation. J. Exp. Med. 160, 877–892 (1984).
41. Hsu, S. M., Raine, L. & Nayak, R. N. Medullary carcinoma of breast: an immunohistochemical study of its lymphoid stroma. Cancer 48, 1368–1376 (1981).
42. Zhao, Y., Niu, C. & Cui, J. Gamma-delta (gammadelta) T cells: friend or foe in cancer development? J. Transl. Med. 16, 3, 10.1186/s12967-017-1378-2 (2018).
43. Joachims, M. L., Chain, J. L., Hooker, S. W., Knott-Craig, C. J. & Thompson, L. F. Human alpha beta and gamma delta thymocyte development: TCR gene rearrangements, intracellular TCR beta expression, and gamma delta developmental potential–differences between men and mice. J. Immunol. 176, 1543–1552 (2006).
44. Reantragoon, R. et al. Antigen-loaded MR1 tetramers define T cell receptor heterogeneity in mucosal-associated invariant T cells. J. Exp. Med. 210, 2305–2320, 10.1084/jem.20130958 (2013).
45. Van Rhijn, I. et al. A conserved human T cell population targets mycobacterial antigens presented by CD1b. Nat. Immunol. 14, 706–713, 10.1038/ni.2630 (2013).
46. Cheuk, S. et al. CD49a expression defines tissue-resident CD8(+) T cells poised for cytotoxic function in human skin. Immunity 46, 287–300, 10.1016/j.immuni.2017.01.009 (2017).
47. Savas, P. et al. Clinical relevance of host immunity in breast cancer: from TILs to the clinic. Nat. Rev. Clin. Oncol. 13, 228–241, 10.1038/nrclinonc.2015.215 (2016).
48. Cha, E. et al. Improved survival with T cell clonotype stability after anti-CTLA-4 treatment in cancer patients. Sci. Transl. Med. 6, 238ra270, 10.1126/scitranslmed.3008211 (2014).
49. Tumeh, P. C. et al. PD-1 blockade induces responses by inhibiting adaptive immune resistance. Nature 515, 568–571, 10.1038/nature13954 (2014).
50. Ji, R. R. et al. An immune-active tumor microenvironment favors clinical response to ipilimumab. Cancer Immunol. Immunother.: 61, 1019–1031, 10.1007/s00262-011-1172-6 (2012).
51. Gee, M. H. et al. Antigen identification for orphan T cell receptors expressed on tumor-infiltrating lymphocytes. Cell 172, 549–563.e516, 10.1016/j.cell.2017.11.043 (2018).
52. Saikia, M. et al. Simultaneous multiplexed amplicon sequencing and transcriptome profiling in single cells. Nat. Methods 16, 59–62, 10.1038/s41592-018-0259-9 (2019).
53. Peterson, V. M. et al. Multiplexed quantification of proteins and transcripts in single cells. Nat. Biotechnol. 35, 936–939, 10.1038/nbt.3973 (2017).
54. Shahi, P., Kim, S. C., Haliburton, J. R., Gartner, Z. J. & Abate, A. R. Abseq: Ultrahigh-throughput single cell protein profiling with droplet microfluidic barcoding. Sci. Rep. 7, 44447, 10.1038/srep44447 (2017).
55. Stoeckius, M. et al. Simultaneous epitope and transcriptome measurement in single cells. Nat. Methods 14, 865–868, 10.1038/nmeth.4380 (2017).
56. Kebschull, J. M. & Zador, A. M. Sources of PCR-induced distortions in high-throughput sequencing data sets. Nucleic Acids Res. 43, e143, 10.1093/nar/gkv717 (2015).
57. Wick, R. R., Judd, L. M. & Holt, K. E. Deepbinner: demultiplexing barcoded Oxford Nanopore reads with deep convolutional neural networks. PLoS Comput. Biol. 14, e1006583, 10.1371/journal.pcbi.1006583 (2018).
58. Gupta, I. et al. Single-cell isoform RNA sequencing characterizes isoforms in thousands of cerebellar cells. Nat. Biotechnol., 10.1038/nbt.4259 (2018).
59. Picelli, S. et al. Full-length RNA-seq from single cells using Smart-seq2. Nat. Protoc. 9, 171–181, 10.1038/nprot.2014.006 (2014).
60. Lefranc, M. P. et al. IMGT(R), the international ImMunoGeneTics information system(R) 25 years on. Nucleic Acids Res. 43, D413–422, 10.1093/nar/gku1056 (2015).
61. Satija, R., Farrell, J. A., Gennert, D., Schier, A. F. & Regev, A. Spatial reconstruction of single-cell gene expression data. Nat. Biotechnol. 33, 495–502, 10.1038/nbt.3192 (2015).
62. Tirosh, I. et al. Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq. Science 352, 189–196, 10.1126/science.aad0501 (2016).
63. Li, H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics, 10.1093/bioinformatics/bty191 (2018).
64. Ye, J., Ma, N., Madden, T. L. & Ostell, J. M. IgBLAST: an immunoglobulin variable domain sequence analysis tool. Nucleic acids Res. 41, W34–40, 10.1093/nar/gkt382 (2013).
65. Camacho, C. et al. BLAST+: architecture and applications. BMC Bioinforma. 10, 421, 10.1186/1471-2105-10-421 (2009).
66. Shannon, P. et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504, 10.1101/gr.1239303 (2003).
67. Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079, 10.1093/bioinformatics/btp352 (2009).