201
|
Thomas N, Best K, Cinelli M, Reich-Zeliger S, Gal H, Shifrut E, Madi A, Friedman N, Shawe-Taylor J, Chain B. Tracking global changes induced in the CD4 T-cell receptor repertoire by immunization with a complex antigen using short stretches of CDR3 protein sequence. ACTA ACUST UNITED AC 2014; 30:3181-8. [PMID: 25095879 PMCID: PMC4221123 DOI: 10.1093/bioinformatics/btu523] [Citation(s) in RCA: 68] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Motivation: The clonal theory of adaptive immunity proposes that immunological responses are encoded by increases in the frequency of lymphocytes carrying antigen-specific receptors. In this study, we measure the frequency of different T-cell receptors (TcR) in CD4 + T cell populations of mice immunized with a complex antigen, killed Mycobacterium tuberculosis, using high throughput parallel sequencing of the TcRβ chain. Our initial hypothesis that immunization would induce repertoire convergence proved to be incorrect, and therefore an alternative approach was developed that allows accurate stratification of TcR repertoires and provides novel insights into the nature of CD4 + T-cell receptor recognition. Results: To track the changes induced by immunization within this heterogeneous repertoire, the sequence data were classified by counting the frequency of different clusters of short (3 or 4) continuous stretches of amino acids within the antigen binding complementarity determining region 3 (CDR3) repertoire of different mice. Both unsupervised (hierarchical clustering) and supervised (support vector machine) analyses of these different distributions of sequence clusters differentiated between immunized and unimmunized mice with 100% efficiency. The CD4 + TcR repertoires of mice 5 and 14 days postimmunization were clearly different from that of unimmunized mice but were not distinguishable from each other. However, the repertoires of mice 60 days postimmunization were distinct both from naive mice and the day 5/14 animals. Our results reinforce the remarkable diversity of the TcR repertoire, resulting in many diverse private TcRs contributing to the T-cell response even in genetically identical mice responding to the same antigen. However, specific motifs defined by short stretches of amino acids within the CDR3 region may determine TcR specificity and define a new approach to TcR sequence classification. Availability and implementation: The analysis was implemented in R and Python, and source code can be found in Supplementary Data. Contact:b.chain@ucl.ac.uk Supplementary information:Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Niclas Thomas
- UCL CoMPLEX, UCL Division of Infection and Immunity, London WC1 6BT, UK, Weizmann Institute of Science, Rehovot 76000, Israel and UCL Computer Science, London WC1E 6BT, UK
| | - Katharine Best
- UCL CoMPLEX, UCL Division of Infection and Immunity, London WC1 6BT, UK, Weizmann Institute of Science, Rehovot 76000, Israel and UCL Computer Science, London WC1E 6BT, UK
| | - Mattia Cinelli
- UCL CoMPLEX, UCL Division of Infection and Immunity, London WC1 6BT, UK, Weizmann Institute of Science, Rehovot 76000, Israel and UCL Computer Science, London WC1E 6BT, UK
| | - Shlomit Reich-Zeliger
- UCL CoMPLEX, UCL Division of Infection and Immunity, London WC1 6BT, UK, Weizmann Institute of Science, Rehovot 76000, Israel and UCL Computer Science, London WC1E 6BT, UK
| | - Hilah Gal
- UCL CoMPLEX, UCL Division of Infection and Immunity, London WC1 6BT, UK, Weizmann Institute of Science, Rehovot 76000, Israel and UCL Computer Science, London WC1E 6BT, UK
| | - Eric Shifrut
- UCL CoMPLEX, UCL Division of Infection and Immunity, London WC1 6BT, UK, Weizmann Institute of Science, Rehovot 76000, Israel and UCL Computer Science, London WC1E 6BT, UK
| | - Asaf Madi
- UCL CoMPLEX, UCL Division of Infection and Immunity, London WC1 6BT, UK, Weizmann Institute of Science, Rehovot 76000, Israel and UCL Computer Science, London WC1E 6BT, UK
| | - Nir Friedman
- UCL CoMPLEX, UCL Division of Infection and Immunity, London WC1 6BT, UK, Weizmann Institute of Science, Rehovot 76000, Israel and UCL Computer Science, London WC1E 6BT, UK
| | - John Shawe-Taylor
- UCL CoMPLEX, UCL Division of Infection and Immunity, London WC1 6BT, UK, Weizmann Institute of Science, Rehovot 76000, Israel and UCL Computer Science, London WC1E 6BT, UK
| | - Benny Chain
- UCL CoMPLEX, UCL Division of Infection and Immunity, London WC1 6BT, UK, Weizmann Institute of Science, Rehovot 76000, Israel and UCL Computer Science, London WC1E 6BT, UK
| |
Collapse
|
202
|
Bashford-Rogers RJM, Palser AL, Idris SF, Carter L, Epstein M, Callard RE, Douek DC, Vassiliou GS, Follows GA, Hubank M, Kellam P. Capturing needles in haystacks: a comparison of B-cell receptor sequencing methods. BMC Immunol 2014; 15:29. [PMID: 25189176 PMCID: PMC4243823 DOI: 10.1186/s12865-014-0029-0] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2014] [Accepted: 07/15/2014] [Indexed: 01/20/2023] Open
Abstract
BACKGROUND Deep-sequencing methods are rapidly developing in the field of B-cell receptor (BCR) and T-cell receptor (TCR) diversity. These promise to revolutionise our understanding of adaptive immune dynamics, identify novel antibodies, and allow monitoring of minimal residual disease. However, different methods for BCR and TCR enrichment and amplification have been proposed. Here we perform the first systematic comparison between different methods of enrichment, amplification and sequencing for generating BCR and TCR repertoires using large sample numbers. RESULTS Resampling from the same RNA or cDNA pool results in highly correlated and reproducible repertoires, but resampling low frequency clones leads to stochastic variance. Repertoires generated by different sequencing methods (454 Roche and Illumina MiSeq) and amplification methods (multiplex PCR, 5' Rapid amplification of cDNA ends (5'RACE), and RNA-capture) are highly correlated, and resulting IgHV gene frequencies between the different methods were not significantly different. Read length has an impact on captured repertoire structure, and ultimately full-length BCR sequences are most informative for repertoire analysis as diversity outside of the CDR is very useful for phylogenetic analysis. Additionally, we show RNA-based BCR repertoires are more informative than using DNA. CONCLUSIONS Repertoires generated by different sequencing and amplification methods are consistent, but we show that read lengths, depths and error profiles should be considered in experimental design, and multiple sampling approaches could be employed to minimise stochastic sampling variation. This detailed investigation of immune repertoire sequencing methods is essential for informing basic and clinical research.
Collapse
Affiliation(s)
| | - Anne L Palser
- />Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA UK
| | - Saad F Idris
- />Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA UK
| | - Lisa Carter
- />Molecular Haematology and Cancer Biology Unit, UCL Institute of Child Health, London, WC1N 1EH UK
| | - Michael Epstein
- />Molecular Haematology and Cancer Biology Unit, UCL Institute of Child Health, London, WC1N 1EH UK
| | - Robin E Callard
- />Molecular Haematology and Cancer Biology Unit, UCL Institute of Child Health, London, WC1N 1EH UK
| | - Daniel C Douek
- />Human Immunology Section, Vaccine Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD 20892 USA
| | - George S Vassiliou
- />Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA UK
| | - George A Follows
- />Department of Hematology, Addenbrooke’s Hospital, Hills Road, Cambridge, CB2 0QQ UK
| | - Mike Hubank
- />Molecular Haematology and Cancer Biology Unit, UCL Institute of Child Health, London, WC1N 1EH UK
| | - Paul Kellam
- />Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA UK
- />Research Department of Infection, Division of Infection and Immunity, University College London, Gower Street, London, WC1E 6BT UK
| |
Collapse
|
203
|
Lavinder JJ, Hoi KH, Reddy ST, Wine Y, Georgiou G. Systematic characterization and comparative analysis of the rabbit immunoglobulin repertoire. PLoS One 2014; 9:e101322. [PMID: 24978027 PMCID: PMC4076286 DOI: 10.1371/journal.pone.0101322] [Citation(s) in RCA: 45] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2014] [Accepted: 06/04/2014] [Indexed: 11/18/2022] Open
Abstract
Rabbits have been used extensively as a model system for the elucidation of the mechanism of immunoglobulin diversification and for the production of antibodies. We employed Next Generation Sequencing to analyze Ig germline V and J gene usage, CDR3 length and amino acid composition, and gene conversion frequencies within the functional (transcribed) IgG repertoire of the New Zealand white rabbit (Oryctolagus cuniculus). Several previously unannotated rabbit heavy chain variable (VH) and light chain variable (VL) germline elements were deduced bioinformatically using multidimensional scaling and k-means clustering methods. We estimated the gene conversion frequency in the rabbit at 23% of IgG sequences with a mean gene conversion tract length of 59±36 bp. Sequencing and gene conversion analysis of the chicken, human, and mouse repertoires revealed that gene conversion occurs much more extensively in the chicken (frequency 70%, tract length 79±57 bp), was observed to a small, yet statistically significant extent in humans, but was virtually absent in mice.
Collapse
Affiliation(s)
- Jason J. Lavinder
- Department of Chemical Engineering, University of Texas at Austin, Austin, Texas, United States of America
- Institute for Cellular and Molecular Biology, University of Texas at Austin, Austin, Texas, United States of America
| | - Kam Hon Hoi
- Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas, United States of America
| | - Sai T. Reddy
- Department of Chemical Engineering, University of Texas at Austin, Austin, Texas, United States of America
- Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas, United States of America
| | - Yariv Wine
- Department of Chemical Engineering, University of Texas at Austin, Austin, Texas, United States of America
- Institute for Cellular and Molecular Biology, University of Texas at Austin, Austin, Texas, United States of America
| | - George Georgiou
- Department of Chemical Engineering, University of Texas at Austin, Austin, Texas, United States of America
- Institute for Cellular and Molecular Biology, University of Texas at Austin, Austin, Texas, United States of America
- Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas, United States of America
- Section of Molecular Genetics and Microbiology, University of Texas at Austin, Austin, Texas, United States of America
| |
Collapse
|
204
|
Abstract
The efficient recognition of pathogens by the adaptive immune system relies on the diversity of receptors displayed at the surface of immune cells. T-cell receptor diversity results from an initial random DNA editing process, called VDJ recombination, followed by functional selection of cells according to the interaction of their surface receptors with self and foreign antigenic peptides. Using high-throughput sequence data from the β-chain of human T-cell receptors, we infer factors that quantify the overall effect of selection on the elements of receptor sequence composition: the V and J gene choice and the length and amino acid composition of the variable region. We find a significant correlation between biases induced by VDJ recombination and our inferred selection factors together with a reduction of diversity during selection. Both effects suggest that natural selection acting on the recombination process has anticipated the selection pressures experienced during somatic evolution. The inferred selection factors differ little between donors or between naive and memory repertoires. The number of sequences shared between donors is well-predicted by our model, indicating a stochastic origin of such public sequences. Our approach is based on a probabilistic maximum likelihood method, which is necessary to disentangle the effects of selection from biases inherent in the recombination process.
Collapse
|
205
|
Antibody repertoire deep sequencing reveals antigen-independent selection in maturing B cells. Proc Natl Acad Sci U S A 2014; 111:E2622-9. [PMID: 24927543 DOI: 10.1073/pnas.1403278111] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Antibody repertoires are known to be shaped by selection for antigen binding. Unexpectedly, we now show that selection also acts on a non-antigen-binding antibody region: the heavy-chain variable (VH)-encoded "elbow" between variable and constant domains. By sequencing 2.8 million recombined heavy-chain genes from immature and mature B-cell subsets in mice, we demonstrate a striking gradient in VH gene use as pre-B cells mature into follicular and then into marginal zone B cells. Cells whose antibodies use VH genes that encode a more flexible elbow are more likely to mature. This effect is distinct from, and exceeds in magnitude, previously described maturation-associated changes in heavy-chain complementarity determining region 3, a key antigen-binding region, which arise from junctional diversity rather than differential VH gene use. Thus, deep sequencing reveals a previously unidentified mode of B-cell selection.
Collapse
|
206
|
Michaeli M, Tabibian-Keissar H, Schiby G, Shahaf G, Pickman Y, Hazanov L, Rosenblatt K, Dunn-Walters DK, Barshack I, Mehr R. Immunoglobulin gene repertoire diversification and selection in the stomach - from gastritis to gastric lymphomas. Front Immunol 2014; 5:264. [PMID: 24917868 PMCID: PMC4042156 DOI: 10.3389/fimmu.2014.00264] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2013] [Accepted: 05/20/2014] [Indexed: 01/06/2023] Open
Abstract
Chronic gastritis is characterized by gastric mucosal inflammation due to autoimmune responses or infection, frequently with Helicobacter pylori. Gastritis with H. pylori background can cause gastric mucosa-associated lymphoid tissue lymphoma (MALT-L), which sometimes further transforms into diffuse large B-cell lymphoma (DLBCL). However, gastric DLBCL can also be initiated de novo. The mechanisms underlying transformation into DLBCL are not completely understood. We analyzed immunoglobulin repertoires and clonal trees to investigate whether and how immunoglobulin gene repertoires, clonal diversification, and selection in gastritis, gastric MALT-L, and DLBCL differ from each other and from normal responses. The two gastritis types (positive or negative for H. pylori) had similarly diverse repertoires. MALT-L dominant clones (defined as the largest clones in each sample) presented higher diversification and longer mutational histories compared with all other conditions. DLBCL dominant clones displayed lower clonal diversification, suggesting the transforming events are triggered by similar responses in different patients. These results are surprising, as we expected to find similarities between the dominant clones of gastritis and MALT-L and between those of MALT-L and DLBCL.
Collapse
Affiliation(s)
- Miri Michaeli
- The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan, Israel
| | - Hilla Tabibian-Keissar
- The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan, Israel
- Department of Pathology, Sheba Medical Center, Ramat Gan, Israel
| | - Ginette Schiby
- Department of Pathology, Sheba Medical Center, Ramat Gan, Israel
| | - Gitit Shahaf
- The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan, Israel
| | - Yishai Pickman
- The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan, Israel
| | - Lena Hazanov
- The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan, Israel
| | | | - Deborah K. Dunn-Walters
- Division of Immunology, Infection, and Inflammatory Diseases, King’s College London School of Medicine, London, UK
| | - Iris Barshack
- Department of Pathology, Sheba Medical Center, Ramat Gan, Israel
- Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Ramit Mehr
- The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan, Israel
| |
Collapse
|
207
|
Murray D, Barnidge D. Characterization of immunoglobulin by mass spectrometry with applications for the clinical laboratory. Crit Rev Clin Lab Sci 2014; 50:91-102. [PMID: 24156651 DOI: 10.3109/10408363.2013.838206] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Studies monitoring immunoglobulin (Ig) antigen specificity have brought to light key Ig biomarkers for immunity, autoimmunity, cancer detection, and immune system function evaluation. A fundamentally new approach to the detection of Igs based on the primary structure of the Ig is beginning to emerge in the literature. This approach has only become feasible in light of advances in proteomics and rapid improvements in mass spectrometry (MS). Driven primarily by the development of Ig pharmaceuticals, Ig MS-based proteomic methods are revealing structural features which were previously unavailable with other characterization techniques. The task of adapting these techniques to clinical chemistry is in its infancy, but these methods have the potential to dramatically alter testing for Ig biomarkers. The purpose of this article is to review the advances that have been made in proteomic characterization of Igs by MS and the early attempts to apply these methods to clinical samples.
Collapse
Affiliation(s)
- David Murray
- Department of Laboratory Medicine and Pathology, Mayo Clinic , Rochester, MN , USA
| | | |
Collapse
|
208
|
Giraud M, Salson M, Duez M, Villenet C, Quief S, Caillault A, Grardel N, Roumier C, Preudhomme C, Figeac M. Fast multiclonal clusterization of V(D)J recombinations from high-throughput sequencing. BMC Genomics 2014; 15:409. [PMID: 24885090 PMCID: PMC4070559 DOI: 10.1186/1471-2164-15-409] [Citation(s) in RCA: 53] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2013] [Accepted: 05/08/2014] [Indexed: 11/28/2022] Open
Abstract
BACKGROUND V(D)J recombinations in lymphocytes are essential for immunological diversity. They are also useful markers of pathologies. In leukemia, they are used to quantify the minimal residual disease during patient follow-up. However, the full breadth of lymphocyte diversity is not fully understood. RESULTS We propose new algorithms that process high-throughput sequencing (HTS) data to extract unnamed V(D)J junctions and gather them into clones for quantification. This analysis is based on a seed heuristic and is fast and scalable because in the first phase, no alignment is performed with germline database sequences. The algorithms were applied to TR γ HTS data from a patient with acute lymphoblastic leukemia, and also on data simulating hypermutations. Our methods identified the main clone, as well as additional clones that were not identified with standard protocols. CONCLUSIONS The proposed algorithms provide new insight into the analysis of high-throughput sequencing data for leukemia, and also to the quantitative assessment of any immunological profile. The methods described here are implemented in a C++ open-source program called Vidjil.
Collapse
Affiliation(s)
- Mathieu Giraud
- />Laboratoire d’Informatique Fondamentale de Lille (LIFL, UMR CNRS 8022, Université Lille 1) and Inria Lille – Cité scientifique – Bâtiment M3, 59655 Villeneuve d’Ascq, France
| | - Mikaël Salson
- />Laboratoire d’Informatique Fondamentale de Lille (LIFL, UMR CNRS 8022, Université Lille 1) and Inria Lille – Cité scientifique – Bâtiment M3, 59655 Villeneuve d’Ascq, France
| | - Marc Duez
- />Laboratoire d’Informatique Fondamentale de Lille (LIFL, UMR CNRS 8022, Université Lille 1) and Inria Lille – Cité scientifique – Bâtiment M3, 59655 Villeneuve d’Ascq, France
- />SIRIC OncoLille, Lille, France
| | - Céline Villenet
- />Functional and Structural Genomic Platform, Université Lille 2, IFR 114 Lille, France
| | - Sabine Quief
- />Functional and Structural Genomic Platform, Université Lille 2, IFR 114 Lille, France
- />Lille Institute for Cancer Research (IRCL), Lille, France
| | - Aurélie Caillault
- />Department of Hematology, Biology and Pathology Center, Lille, France
| | - Nathalie Grardel
- />Department of Hematology, Biology and Pathology Center, Lille, France
| | - Christophe Roumier
- />Department of Hematology, Biology and Pathology Center, Lille, France
- />Inserm U-837, Cancer Research Institute, Lille, France
| | - Claude Preudhomme
- />Department of Hematology, Biology and Pathology Center, Lille, France
- />Inserm U-837, Cancer Research Institute, Lille, France
| | - Martin Figeac
- />Functional and Structural Genomic Platform, Université Lille 2, IFR 114 Lille, France
| |
Collapse
|
209
|
Galson JD, Pollard AJ, Trück J, Kelly DF. Studying the antibody repertoire after vaccination: practical applications. Trends Immunol 2014; 35:319-31. [PMID: 24856924 DOI: 10.1016/j.it.2014.04.005] [Citation(s) in RCA: 81] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2013] [Revised: 04/17/2014] [Accepted: 04/28/2014] [Indexed: 12/25/2022]
Abstract
Nearly all licensed vaccines have been developed to confer protection against infectious diseases by stimulating the production of antibodies by B cells, but the nature of a successful antibody response has been difficult to capture. Recent advances in next-generation sequencing (NGS) technology have allowed high-resolution characterization of the antibody repertoire, and of the changes that occur following vaccination. These approaches have yielded important insights into the B cell response, and have raised the possibility of using specific antibody sequences as measures of vaccine immunogenicity. Here, we review recent findings based on antibody repertoire sequencing, and discuss potential applications of these new technologies and of the analyses of the increasing volume of antibody sequence data in the context of vaccine development.
Collapse
Affiliation(s)
- Jacob D Galson
- Oxford Vaccine Group, Department of Paediatrics, University of Oxford and the NIHR Oxford Biomedical Research Centre, Oxford, UK.
| | - Andrew J Pollard
- Oxford Vaccine Group, Department of Paediatrics, University of Oxford and the NIHR Oxford Biomedical Research Centre, Oxford, UK.
| | - Johannes Trück
- Oxford Vaccine Group, Department of Paediatrics, University of Oxford and the NIHR Oxford Biomedical Research Centre, Oxford, UK
| | - Dominic F Kelly
- Oxford Vaccine Group, Department of Paediatrics, University of Oxford and the NIHR Oxford Biomedical Research Centre, Oxford, UK
| |
Collapse
|
210
|
Menzel U, Greiff V, Khan TA, Haessler U, Hellmann I, Friedensohn S, Cook SC, Pogson M, Reddy ST. Comprehensive evaluation and optimization of amplicon library preparation methods for high-throughput antibody sequencing. PLoS One 2014; 9:e96727. [PMID: 24809667 PMCID: PMC4014543 DOI: 10.1371/journal.pone.0096727] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2014] [Accepted: 04/10/2014] [Indexed: 11/18/2022] Open
Abstract
High-throughput sequencing (HTS) of antibody repertoire libraries has become a powerful tool in the field of systems immunology. However, numerous sources of bias in HTS workflows may affect the obtained antibody repertoire data. A crucial step in antibody library preparation is the addition of short platform-specific nucleotide adapter sequences. As of yet, the impact of the method of adapter addition on experimental library preparation and the resulting antibody repertoire HTS datasets has not been thoroughly investigated. Therefore, we compared three standard library preparation methods by performing Illumina HTS on antibody variable heavy genes from murine antibody-secreting cells. Clonal overlap and rank statistics demonstrated that the investigated methods produced equivalent HTS datasets. PCR-based methods were experimentally superior to ligation with respect to speed, efficiency, and practicality. Finally, using a two-step PCR based method we established a protocol for antibody repertoire library generation, beginning from inputs as low as 1 ng of total RNA. In summary, this study represents a major advance towards a standardized experimental framework for antibody HTS, thus opening up the potential for systems-based, cross-experiment meta-analyses of antibody repertoires.
Collapse
Affiliation(s)
- Ulrike Menzel
- Department of Biosystems Science and Engineering, ETH Zürich, Basel, Switzerland
| | - Victor Greiff
- Department of Biosystems Science and Engineering, ETH Zürich, Basel, Switzerland
| | - Tarik A Khan
- Department of Biosystems Science and Engineering, ETH Zürich, Basel, Switzerland
| | - Ulrike Haessler
- Department of Biosystems Science and Engineering, ETH Zürich, Basel, Switzerland
| | - Ina Hellmann
- Department of Biosystems Science and Engineering, ETH Zürich, Basel, Switzerland
| | - Simon Friedensohn
- Department of Biosystems Science and Engineering, ETH Zürich, Basel, Switzerland
| | - Skylar C Cook
- Department of Biosystems Science and Engineering, ETH Zürich, Basel, Switzerland
| | - Mark Pogson
- Department of Biosystems Science and Engineering, ETH Zürich, Basel, Switzerland
| | - Sai T Reddy
- Department of Biosystems Science and Engineering, ETH Zürich, Basel, Switzerland
| |
Collapse
|
211
|
de Costa D, Broodman I, Calame W, Stingl C, Dekker LJM, Vernhout RM, de Koning HJ, Hoogsteden HC, Smitt PAES, van Klaveren RJ, Luider TM, VanDuijn MM. Peptides from the variable region of specific antibodies are shared among lung cancer patients. PLoS One 2014; 9:e96029. [PMID: 24787687 PMCID: PMC4006902 DOI: 10.1371/journal.pone.0096029] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2013] [Accepted: 04/03/2014] [Indexed: 11/19/2022] Open
Abstract
Late diagnosis of lung cancer is still the main reason for high mortality rates in lung cancer. Lung cancer is a heterogeneous disease which induces an immune response to different tumor antigens. Several methods for searching autoantibodies have been described that are based on known purified antigen panels. The aim of our study is to find evidence that parts of the antigen-binding-domain of antibodies are shared among lung cancer patients. This was investigated by a novel approach based on sequencing antigen-binding-fragments (Fab) of immunoglobulins using proteomic techniques without the need of previously known antigen panels. From serum of 93 participants of the NELSON trial IgG was isolated and subsequently digested into Fab and Fc. Fab was purified from the digested mixture by SDS-PAGE. The Fab containing gel-bands were excised, tryptic digested and measured on a nano-LC-Orbitrap-Mass-spectrometry system. Multivariate analysis of the mass spectrometry data by linear canonical discriminant analysis combined with stepwise logistic regression resulted in a 12-antibody-peptide model which was able to distinguish lung cancer patients from controls in a high risk population with a sensitivity of 84% and specificity of 90%. With our Fab-purification combined Orbitrap-mass-spectrometry approach, we found peptides from the variable-parts of antibodies which are shared among lung cancer patients.
Collapse
Affiliation(s)
- Dominique de Costa
- Department of Pulmonology, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Ingrid Broodman
- Department of Clinical Chemistry, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Wim Calame
- StatistiCal BV, Wassenaar, The Netherlands
| | - Christoph Stingl
- Department of Neurology, Erasmus Medical Center, Rotterdam, The Netherlands
| | | | - René M. Vernhout
- Unit of Trials and Statistics, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Harry J. de Koning
- Unit of Public Health, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Henk C. Hoogsteden
- Department of Pulmonology, Erasmus Medical Center, Rotterdam, The Netherlands
| | | | - Rob J. van Klaveren
- Department of Pulmonology, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Theo M. Luider
- Department of Neurology, Erasmus Medical Center, Rotterdam, The Netherlands
- * E-mail:
| | | |
Collapse
|
212
|
Boutz DR, Horton AP, Wine Y, Lavinder JJ, Georgiou G, Marcotte EM. Proteomic identification of monoclonal antibodies from serum. Anal Chem 2014; 86:4758-66. [PMID: 24684310 PMCID: PMC4033631 DOI: 10.1021/ac4037679] [Citation(s) in RCA: 53] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Characterizing the in vivo dynamics of the polyclonal antibody repertoire in serum, such as that which might arise in response to stimulation with an antigen, is difficult due to the presence of many highly similar immunoglobulin proteins, each specified by distinct B lymphocytes. These challenges have precluded the use of conventional mass spectrometry for antibody identification based on peptide mass spectral matches to a genomic reference database. Recently, progress has been made using bottom-up analysis of serum antibodies by nanoflow liquid chromatography/high-resolution tandem mass spectrometry combined with a sample-specific antibody sequence database generated by high-throughput sequencing of individual B cell immunoglobulin variable domains (V genes). Here, we describe how intrinsic features of antibody primary structure, most notably the interspersed segments of variable and conserved amino acid sequences, generate recurring patterns in the corresponding peptide mass spectra of V gene peptides, greatly complicating the assignment of correct sequences to mass spectral data. We show that the standard method of decoy-based error modeling fails to account for the error introduced by these highly similar sequences, leading to a significant underestimation of the false discovery rate. Because of these effects, antibody-derived peptide mass spectra require increased stringency in their interpretation. The use of filters based on the mean precursor ion mass accuracy of peptide-spectrum matches is shown to be particularly effective in distinguishing between "true" and "false" identifications. These findings highlight important caveats associated with the use of standard database search and error-modeling methods with nonstandard data sets and custom sequence databases.
Collapse
Affiliation(s)
- Daniel R Boutz
- Center for Systems & Synthetic Biology, †Institute for Cellular and Molecular Biology, ⊥Department of Biomedical Engineering, §Department of Chemical Engineering, and ∥Department of Molecular Biosciences, University of Texas at Austin , Austin, Texas 78712, United States
| | | | | | | | | | | |
Collapse
|
213
|
Tkach KE, Barik D, Voisinne G, Malandro N, Hathorn MM, Cotari JW, Vogel R, Merghoub T, Wolchok J, Krichevsky O, Altan-Bonnet G. T cells translate individual, quantal activation into collective, analog cytokine responses via time-integrated feedbacks. eLife 2014; 3:e01944. [PMID: 24719192 PMCID: PMC3980879 DOI: 10.7554/elife.01944] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022] Open
Abstract
Variability within isogenic T cell populations yields heterogeneous ‘local’ signaling responses to shared antigenic stimuli, but responding clones may communicate ‘global’ antigen load through paracrine messengers, such as cytokines. Such coordination of individual cell responses within multicellular populations is critical for accurate collective reactions to shared environmental cues. However, cytokine production may saturate as a function of antigen input, or be dominated by the precursor frequency of antigen-specific T cells. Surprisingly, we found that T cells scale their collective output of IL-2 to total antigen input over a large dynamic range, independently of population size. Through experimental quantitation and computational modeling, we demonstrate that this scaling is enforced by an inhibitory cross-talk between antigen and IL-2 signaling, and a nonlinear acceleration of IL-2 secretion per cell. Our study reveals how time-integration of these regulatory loops within individual cell signaling generates scaled collective responses and can be leveraged for immune monitoring. DOI:http://dx.doi.org/10.7554/eLife.01944.001 The cells of the immune system face the challenge of removing viruses and other pathogens without endangering healthy tissues. Cells called T cells plays a variety of roles in the immune response: some T cells directly destroy infected cells, some recruit other cells called phagocytes to the site of infection, and some release small proteins called cytokines. These cytokines help cells to communicate with other cells and, therefore, to tailor the overall immune responses to deal with a particular pathogen. It is known that mammals are capable of adjusting the T cell response to match the overall severity of an infection. However, it is not clear how individual T cells coordinate their seemingly binary response—they are either activated when they recognize a pathogen, or they are not activated—into a response at the collective cell level that can be varied continuously over a wide range of values. Here, Tkach et al. show that T cell populations match their production of the cytokine interleukin 2 (IL-2) to the abundance of antigens—molecules released by the pathogen—over an unexpectedly large range of concentrations. Through a combination of experimental and computational analyses, Tkach et al. identified two novel IL-2 feedback loops that help to generate the correct quantity of cytokine, irrespective of the total number of T cells. Furthermore, this model can be used to estimate antigen quantities within diseased tissues. The work of Tkach et al. illustrates the potential of feedback integration in cell signalling and gene regulation as a mechanism to allow cellular populations to respond to environmental stimuli in a graded, collective fashion. DOI:http://dx.doi.org/10.7554/eLife.01944.002
Collapse
Affiliation(s)
- Karen E Tkach
- Program in Computational Biology and Immunology, Memorial Sloan Kettering Cancer Center, New York, United States
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
214
|
High-resolution antibody dynamics of vaccine-induced immune responses. Proc Natl Acad Sci U S A 2014; 111:4928-33. [PMID: 24639495 DOI: 10.1073/pnas.1323862111] [Citation(s) in RCA: 126] [Impact Index Per Article: 12.6] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Abstract
The adaptive immune system confers protection by generating a diverse repertoire of antibody receptors that are rapidly expanded and contracted in response to specific targets. Next-generation DNA sequencing now provides the opportunity to survey this complex and vast repertoire. In the present work, we describe a set of tools for the analysis of antibody repertoires and their application to elucidating the dynamics of the response to viral vaccination in human volunteers. By analyzing data from 38 separate blood samples across 2 y, we found that the use of the germ-line library of V and J segments is conserved between individuals over time. Surprisingly, there appeared to be no correlation between the use level of a particular VJ combination and degree of expansion. We found the antibody RNA repertoire in each volunteer to be highly dynamic, with each individual displaying qualitatively different response dynamics. By using combinatorial phage display, we screened selected VH genes paired with their corresponding VL library for affinity against the vaccine antigens. Altogether, this work presents an additional set of tools for profiling the human antibody repertoire and demonstrates characterization of the fast repertoire dynamics through time in multiple individuals responding to an immune challenge.
Collapse
|
215
|
Kodangattil S, Huard C, Ross C, Li J, Gao H, Mascioni A, Hodawadekar S, Naik S, Min-debartolo J, Visintin A, Almagro JC. The functional repertoire of rabbit antibodies and antibody discovery via next-generation sequencing. MAbs 2014; 6:628-36. [PMID: 24481222 DOI: 10.4161/mabs.28059] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
Abstract
To gain insight into the functional antibody repertoire of rabbits, the VH and VL repertoires of bone marrow (BM) and spleen (SP) of a naïve New Zealand White rabbit (NZW; Oryctolagus cuniculus) and that of lymphocytes collected from a NZW rabbit immunized (IM) with a 16-mer peptide were deep-sequenced. Two closely related genes, IGHV1S40 (VH1a3) and IGHV1S45 (VH4), were found to dominate (~90%) the VH repertoire of BM and SP, whereas, IGHV1S69 (VH1a1) contributed significantly (~40%) to IM. BM and SP antibodies recombined predominantly with IGHJ4. A significant proportion (~30%) of IM sequences recombined with IGHJ2. The VK repertoire was encoded by nine IGKV genes recombined with one IGKJ gene, IGKJ1. No significant bias in the VK repertoire of the BM, SP and IM samples was observed. The complementarity-determining region (CDR)-H3 and -L3 length distributions were similar in the three samples following a Gaussian curve with average length of 12.2 ± 2.4 and 11.1 ± 1.1 amino acids, respectively. The amino acid composition of the predominant CDR-H3 and -L3 loop lengths was similar to that of humans and mice, rich in Tyr, Gly, Ser and, in some specific positions, Asp. The average number of mutations along the IGHV/KV genes was similar in BM, SP and IM; close to 12 and 15 mutations for VH and VL, respectively. A monoclonal antibody specific for the peptide used as immunogen was obtained from the IM rabbit. The CDR-H3 sequence was found in 1,559 of 61,728 (2.5%) sequences, at position 10, in the rank order of the CDR-H3 frequencies. The CDR-L3 was found in 24 of 11,215 (0.2%) sequences, ranking 102. No match was found in the BM and SP samples, indicating positive selection for the hybridoma sequence. Altogether, these findings lay foundations for engineering of rabbit V regions to enhance their potential as therapeutics, i.e., design of strategies for selection of specific rabbit V regions from NGS data mining, humanization and design of libraries for affinity maturation campaigns.
Collapse
Affiliation(s)
| | | | | | - Jian Li
- CTI-Boston; Pfizer Inc.; Boston, MA USA
| | | | | | | | | | | | | | | |
Collapse
|
216
|
Weiss-Ottolenghi Y, Gershoni JM. Profiling the IgOme: meeting the challenge. FEBS Lett 2014; 588:318-25. [PMID: 24239539 PMCID: PMC7094557 DOI: 10.1016/j.febslet.2013.11.005] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2013] [Revised: 11/06/2013] [Accepted: 11/06/2013] [Indexed: 02/03/2023]
Abstract
The entire repertoire of antibodies in our serum, the IgOme, is a historical record of our past experiences and a reflection of our immune status at any given moment. Understanding the dynamics of the IgOme and how the diversity and specificities of serum antibodies change in response to disease and maintenance of homeostasis can directly impact the ability to design and develop novel vaccines, diagnostics and therapeutics. Here we review both direct and indirect methodologies that are being developed to map the complexity and specificities of the antibodies in polyclonal serum - the IgOme.
Collapse
Affiliation(s)
- Yael Weiss-Ottolenghi
- Department of Cell Research and Immunology, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv 69978, Israel
| | - Jonathan M Gershoni
- Department of Cell Research and Immunology, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv 69978, Israel.
| |
Collapse
|
217
|
The promise and challenge of high-throughput sequencing of the antibody repertoire. Nat Biotechnol 2014; 32:158-68. [PMID: 24441474 PMCID: PMC4113560 DOI: 10.1038/nbt.2782] [Citation(s) in RCA: 471] [Impact Index Per Article: 47.1] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2013] [Accepted: 12/04/2013] [Indexed: 12/16/2022]
Abstract
Georgiou and colleagues discuss rapidly evolving methods for high-throughput sequencing of the antibody repertoire, and how the resulting data may be applied to answer basic and translational research questions. Efforts to determine the antibody repertoire encoded by B cells in the blood or lymphoid organs using high-throughput DNA sequencing technologies have been advancing at an extremely rapid pace and are transforming our understanding of humoral immune responses. Information gained from high-throughput DNA sequencing of immunoglobulin genes (Ig-seq) can be applied to detect B-cell malignancies with high sensitivity, to discover antibodies specific for antigens of interest, to guide vaccine development and to understand autoimmunity. Rapid progress in the development of experimental protocols and informatics analysis tools is helping to reduce sequencing artifacts, to achieve more precise quantification of clonal diversity and to extract the most pertinent biological information. That said, broader application of Ig-seq, especially in clinical settings, will require the development of a standardized experimental design framework that will enable the sharing and meta-analysis of sequencing data generated by different laboratories.
Collapse
|
218
|
Jiang N, He J, Weinstein JA, Penland L, Sasaki S, He XS, Dekker CL, Zheng NY, Huang M, Sullivan M, Wilson PC, Greenberg HB, Davis MM, Fisher DS, Quake SR. Lineage structure of the human antibody repertoire in response to influenza vaccination. Sci Transl Med 2014; 5:171ra19. [PMID: 23390249 DOI: 10.1126/scitranslmed.3004794] [Citation(s) in RCA: 256] [Impact Index Per Article: 25.6] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
The human antibody repertoire is one of the most important defenses against infectious disease, and the development of vaccines has enabled the conferral of targeted protection to specific pathogens. However, there are many challenges to measuring and analyzing the immunoglobulin sequence repertoire, including that each B cell's genome encodes a distinct antibody sequence, that the antibody repertoire changes over time, and the high similarity between antibody sequences. We have addressed these challenges by using high-throughput long read sequencing to perform immunogenomic characterization of expressed human antibody repertoires in the context of influenza vaccination. Informatic analysis of 5 million antibody heavy chain sequences from healthy individuals allowed us to perform global characterizations of isotype distributions, determine the lineage structure of the repertoire, and measure age- and antigen-related mutational activity. Our analysis of the clonal structure and mutational distribution of individuals' repertoires shows that elderly subjects have a decreased number of lineages but an increased prevaccination mutation load in their repertoire and that some of these subjects have an oligoclonal character to their repertoire in which the diversity of the lineages is greatly reduced relative to younger subjects. We have thus shown that global analysis of the immune system's clonal structure provides direct insight into the effects of vaccination and provides a detailed molecular portrait of age-related effects.
Collapse
Affiliation(s)
- Ning Jiang
- Department of Bioengineering, Stanford University School of Medicine, Stanford, CA 94305, USA
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
219
|
Thomas N, Heather J, Pollara G, Simpson N, Matjeka T, Shawe-Taylor J, Noursadeghi M, Chain B. The immune system as a biomonitor: explorations in innate and adaptive immunity. Interface Focus 2014; 3:20120099. [PMID: 24427535 DOI: 10.1098/rsfs.2012.0099] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
The human immune system has a highly complex, multi-layered structure which has evolved to detect and respond to changes in the internal microenvironment of the body. Recognition occurs at the molecular or submolecular scale, via classical reversible receptor-ligand interactions, and can lead to a response with great sensitivity and speed. Remarkably, recognition is coupled to memory, such that responses are modulated by events which occurred years or even decades before. Although the immune system in general responds differently and more vigorously to stimuli entering the body from the outside (e.g. infections), this is an emergent property of the system: many of the recognition molecules themselves have no inherent bias towards external stimuli (non-self) but also bind targets found within the body (self). It is quite clear that the immune response registers pathophysiological changes in general. Cancer, wounding and chronic tissue injury are some obvious examples. Against this background, the immune system 'state' tracks the internal processes of the body, and is likely to encode information regarding both current and past disease processes. Moreover, the distributed nature of most immune responses (e.g. typically involving lymphoid tissue, non-lymphoid tissue, bone marrow, blood, extracellular interstitial spaces, etc.) means that many of the changes associated with immune responses are manifested systemically, and specifically can be detected in blood. This provides a very convenient route to sampling immune cells. We consider two different and complementary ways of querying the human immune 'state' using high-dimensional genomic screening methodologies, and discuss the potentials of these approaches and some of the technological and computational challenges to be overcome.
Collapse
Affiliation(s)
- Niclas Thomas
- CoMPLEX, University College London, Gower Street, London WC1E 6BT, UK
| | - James Heather
- Division of Infection and Immunity, University College London, Gower Street, London WC1E 6BT, UK
| | - Gabriel Pollara
- Division of Infection and Immunity, University College London, Gower Street, London WC1E 6BT, UK
| | - Nandi Simpson
- Division of Infection and Immunity, University College London, Gower Street, London WC1E 6BT, UK
| | - Theres Matjeka
- Division of Infection and Immunity, University College London, Gower Street, London WC1E 6BT, UK
| | - John Shawe-Taylor
- Department of Computer Science, University College London, Gower Street, London WC1E 6BT, UK
| | - Mahdad Noursadeghi
- Division of Infection and Immunity, University College London, Gower Street, London WC1E 6BT, UK
| | - Benjamin Chain
- Division of Infection and Immunity, University College London, Gower Street, London WC1E 6BT, UK
| |
Collapse
|
220
|
Parameswaran P, Liu Y, Roskin KM, Jackson KKL, Dixit VP, Lee JY, Artiles KL, Zompi S, Vargas MJ, Simen BB, Hanczaruk B, McGowan KR, Tariq MA, Pourmand N, Koller D, Balmaseda A, Boyd SD, Harris E, Fire AZ. Convergent antibody signatures in human dengue. Cell Host Microbe 2013; 13:691-700. [PMID: 23768493 DOI: 10.1016/j.chom.2013.05.008] [Citation(s) in RCA: 192] [Impact Index Per Article: 17.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2013] [Revised: 04/03/2013] [Accepted: 04/30/2013] [Indexed: 01/05/2023]
Abstract
Dengue is the most prevalent mosquito-borne viral disease in humans, and the lack of early prognostics, vaccines, and therapeutics contributes to immense disease burden. To identify patterns that could be used for sequence-based monitoring of the antibody response to dengue, we examined antibody heavy-chain gene rearrangements in longitudinal peripheral blood samples from 60 dengue patients. Comparing signatures between acute dengue, postrecovery, and healthy samples, we found increased expansion of B cell clones in acute dengue patients, with higher overall clonality in secondary infection. Additionally, we observed consistent antibody sequence features in acute dengue in the highly variable major antigen-binding determinant, complementarity-determining region 3 (CDR3), with specific CDR3 sequences highly enriched in acute samples compared to postrecovery, healthy, or non-dengue samples. Dengue thus provides a striking example of a human viral infection where convergent immune signatures can be identified in multiple individuals. Such signatures could facilitate surveillance of immunological memory in communities.
Collapse
Affiliation(s)
- Poornima Parameswaran
- Division of Infectious Diseases and Vaccinology, School of Public Health, University of California, Berkeley, Berkeley, CA 94720-3370, USA
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
221
|
Error analysis of deep sequencing of phage libraries: peptides censored in sequencing. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2013; 2013:491612. [PMID: 24416071 PMCID: PMC3876701 DOI: 10.1155/2013/491612] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/06/2013] [Accepted: 07/30/2013] [Indexed: 01/12/2023]
Abstract
Next-generation sequencing techniques empower selection of ligands from phage-display libraries because they can detect low abundant clones and quantify changes in the copy numbers of clones without excessive selection rounds. Identification of errors in deep sequencing data is the most critical step in this process because these techniques have error rates >1%. Mechanisms that yield errors in Illumina and other techniques have been proposed, but no reports to date describe error analysis in phage libraries. Our paper focuses on error analysis of 7-mer peptide libraries sequenced by Illumina method. Low theoretical complexity of this phage library, as compared to complexity of long genetic reads and genomes, allowed us to describe this library using convenient linear vector and operator framework. We describe a phage library as N × 1 frequency vector n = ||ni||, where ni is the copy number of the ith sequence and N is the theoretical diversity, that is, the total number of all possible sequences. Any manipulation to the library is an operator acting on n. Selection, amplification, or sequencing could be described as a product of a N × N matrix and a stochastic sampling operator (Sa). The latter is a random diagonal matrix that describes sampling of a library. In this paper, we focus on the properties of Sa and use them to define the sequencing operator (Seq). Sequencing without any bias and errors is Seq = Sa IN, where IN is a N × N unity matrix. Any bias in sequencing changes IN to a nonunity matrix. We identified a diagonal censorship matrix (CEN), which describes elimination or statistically significant downsampling, of specific reads during the sequencing process.
Collapse
|
222
|
Xiao M, Ponraj P, Chen W, Kessing B, Dimitrov DS. Deep sequencing and Circos analyses of antibody libraries reveal antigen-driven selection of Ig VH genes during HIV-1 infection. Exp Mol Pathol 2013; 95:357-63. [PMID: 24158018 PMCID: PMC3889869 DOI: 10.1016/j.yexmp.2013.10.004] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2013] [Accepted: 10/11/2013] [Indexed: 01/21/2023]
Abstract
The vast diversity of antibody repertoires is largely attributed to heavy chain (V(H)) recombination of variable (V), diversity (D) and joining (J) gene segments. We used 454 sequencing information of the variable domains of the antibody heavy chain repertoires from neonates, normal adults and an HIV-1-infected individual, to analyze, with Circos software, the VDJ pairing patterns at birth, adulthood and a time-dependent response to HIV-1 infection. Our comparative analyses of the Ig VDJ repertoires from these libraries indicated that, from birth to adulthood, VDJ recombination patterns remain the same with some slight changes, whereas some V(H) families are selected and preferentially expressed after long-term infection with HIV-1. We also demonstrated that the immune system responds to HIV-1 chronic infection by selectively expanding certain HV families in an attempt to combat infection. Our findings may have implications for understanding immune responses in pathology as well as for development of new therapeutics and vaccines.
Collapse
Affiliation(s)
- Madelyne Xiao
- Urbana High School, Ijamsville, Maryland 21754, USA
- BSP CCR Genetics Core, Science Applications International Corporation-Frederick, Inc., NCI-Frederick, Frederick, MD 21702, USA
| | - Prabakaran Ponraj
- Protein Interactions Group, Center for Cancer Research Cancer and Inflammation Program, Science Applications International Corporation-Frederick, Inc., NCI-Frederick, Frederick, MD 21702, USA
- Basic Research Program, Science Applications International Corporation-Frederick, Inc., NCI-Frederick, Frederick, MD 21702, USA
| | - Weizao Chen
- Protein Interactions Group, Center for Cancer Research Cancer and Inflammation Program, Science Applications International Corporation-Frederick, Inc., NCI-Frederick, Frederick, MD 21702, USA
| | - Bailey Kessing
- BSP CCR Genetics Core, Science Applications International Corporation-Frederick, Inc., NCI-Frederick, Frederick, MD 21702, USA
| | - Dimiter S. Dimitrov
- Protein Interactions Group, Center for Cancer Research Cancer and Inflammation Program, Science Applications International Corporation-Frederick, Inc., NCI-Frederick, Frederick, MD 21702, USA
| |
Collapse
|
223
|
Six A, Mariotti-Ferrandiz ME, Chaara W, Magadan S, Pham HP, Lefranc MP, Mora T, Thomas-Vaslin V, Walczak AM, Boudinot P. The past, present, and future of immune repertoire biology - the rise of next-generation repertoire analysis. Front Immunol 2013; 4:413. [PMID: 24348479 PMCID: PMC3841818 DOI: 10.3389/fimmu.2013.00413] [Citation(s) in RCA: 113] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2013] [Accepted: 11/12/2013] [Indexed: 01/09/2023] Open
Abstract
T and B cell repertoires are collections of lymphocytes, each characterized by its antigen-specific receptor. We review here classical technologies and analysis strategies developed to assess immunoglobulin (IG) and T cell receptor (TR) repertoire diversity, and describe recent advances in the field. First, we describe the broad range of available methodological tools developed in the past decades, each of which answering different questions and showing complementarity for progressive identification of the level of repertoire alterations: global overview of the diversity by flow cytometry, IG repertoire descriptions at the protein level for the identification of IG reactivities, IG/TR CDR3 spectratyping strategies, and related molecular quantification or dynamics of T/B cell differentiation. Additionally, we introduce the recent technological advances in molecular biology tools allowing deeper analysis of IG/TR diversity by next-generation sequencing (NGS), offering systematic and comprehensive sequencing of IG/TR transcripts in a short amount of time. NGS provides several angles of analysis such as clonotype frequency, CDR3 diversity, CDR3 sequence analysis, V allele identification with a quantitative dimension, therefore requiring high-throughput analysis tools development. In this line, we discuss the recent efforts made for nomenclature standardization and ontology development. We then present the variety of available statistical analysis and modeling approaches developed with regards to the various levels of diversity analysis, and reveal the increasing sophistication of those modeling approaches. To conclude, we provide some examples of recent mathematical modeling strategies and perspectives that illustrate the active rise of a "next-generation" of repertoire analysis.
Collapse
Affiliation(s)
- Adrien Six
- UPMC University Paris 06, UMR 7211, Immunology-Immunopathology-Immunotherapy (I3) , Paris , France ; CNRS, UMR 7211, Immunology-Immunopathology-Immunotherapy (I3) , Paris , France ; INSERM, UMR_S 959, Immunology-Immunopathology-Immunotherapy (I3) , Paris , France ; AP-HP, Hôpital Pitié-Salpêtrière, CIC-BTi Biotherapy , Paris , France ; AP-HP, Hôpital Pitié-Salpêtrière, Département Hospitalo-Universitaire (DHU), Inflammation-Immunopathology-Biotherapy (i2B) , Paris , France
| | - Maria Encarnita Mariotti-Ferrandiz
- UPMC University Paris 06, UMR 7211, Immunology-Immunopathology-Immunotherapy (I3) , Paris , France ; CNRS, UMR 7211, Immunology-Immunopathology-Immunotherapy (I3) , Paris , France ; INSERM, UMR_S 959, Immunology-Immunopathology-Immunotherapy (I3) , Paris , France ; AP-HP, Hôpital Pitié-Salpêtrière, Département Hospitalo-Universitaire (DHU), Inflammation-Immunopathology-Biotherapy (i2B) , Paris , France
| | - Wahiba Chaara
- UPMC University Paris 06, UMR 7211, Immunology-Immunopathology-Immunotherapy (I3) , Paris , France ; CNRS, UMR 7211, Immunology-Immunopathology-Immunotherapy (I3) , Paris , France ; INSERM, UMR_S 959, Immunology-Immunopathology-Immunotherapy (I3) , Paris , France ; AP-HP, Hôpital Pitié-Salpêtrière, CIC-BTi Biotherapy , Paris , France ; AP-HP, Hôpital Pitié-Salpêtrière, Département Hospitalo-Universitaire (DHU), Inflammation-Immunopathology-Biotherapy (i2B) , Paris , France
| | - Susana Magadan
- Institut National de la Recherche Agronomique, Unité de Virologie et Immunologie Moléculaires , Jouy-en-Josas , France
| | - Hang-Phuong Pham
- UPMC University Paris 06, UMR 7211, Immunology-Immunopathology-Immunotherapy (I3) , Paris , France ; CNRS, UMR 7211, Immunology-Immunopathology-Immunotherapy (I3) , Paris , France
| | - Marie-Paule Lefranc
- IMGT®, The International ImMunoGeneTics Information System®, Institut de Génétique Humaine, UPR CNRS 1142, Université Montpellier 2 , Montpellier , France
| | - Thierry Mora
- Laboratoire de Physique Statistique, UMR8550, CNRS and Ecole Normale Supérieure , Paris , France
| | - Véronique Thomas-Vaslin
- UPMC University Paris 06, UMR 7211, Immunology-Immunopathology-Immunotherapy (I3) , Paris , France ; CNRS, UMR 7211, Immunology-Immunopathology-Immunotherapy (I3) , Paris , France ; INSERM, UMR_S 959, Immunology-Immunopathology-Immunotherapy (I3) , Paris , France ; AP-HP, Hôpital Pitié-Salpêtrière, Département Hospitalo-Universitaire (DHU), Inflammation-Immunopathology-Biotherapy (i2B) , Paris , France
| | - Aleksandra M Walczak
- Laboratoire de Physique Théorique, UMR8549, CNRS and Ecole Normale Supérieure , Paris , France
| | - Pierre Boudinot
- Institut National de la Recherche Agronomique, Unité de Virologie et Immunologie Moléculaires , Jouy-en-Josas , France
| |
Collapse
|
224
|
Yin L, Hou W, Liu L, Cai Y, Wallet MA, Gardner BP, Chang K, Lowe AC, Rodriguez CA, Sriaroon P, Farmerie WG, Sleasman JW, Goodenow MM. IgM Repertoire Biodiversity is Reduced in HIV-1 Infection and Systemic Lupus Erythematosus. Front Immunol 2013; 4:373. [PMID: 24298273 PMCID: PMC3828670 DOI: 10.3389/fimmu.2013.00373] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2013] [Accepted: 10/30/2013] [Indexed: 12/25/2022] Open
Abstract
Background: HIV-1 infection or systemic lupus erythematosus (SLE) disrupt B cell homeostasis, reduce memory B cells, and impair function of IgG and IgM antibodies. Objective: To determine how disturbances in B cell populations producing polyclonal antibodies relate to the IgM repertoire, the IgM transcriptome in health and disease was explored at the complementarity determining region 3 (CDRH3) sequence level. Methods: 454-deep pyrosequencing in combination with a novel analysis pipeline was applied to define populations of IGHM CDRH3 sequences based on absence or presence of somatic hypermutations (SHM) in peripheral blood B cells. Results: HIV or SLE subjects have reduced biodiversity within their IGHM transcriptome compared to healthy subjects, mainly due to a significant decrease in the number of unique combinations of alleles, although recombination machinery was intact. While major differences between sequences without or with SHM occurred among all groups, IGHD and IGHJ allele use, CDRH3 length distribution, or generation of SHM were similar among study cohorts. Antiretroviral therapy failed to normalize IGHM biodiversity in HIV-infected individuals. All subjects had a low frequency of allelic combinations within the IGHM repertoire similar to known broadly neutralizing HIV-1 antibodies. Conclusion: Polyclonal expansion would decrease overall IgM biodiversity independent of other mechanisms for development of the B cell repertoire. Applying deep sequencing as a strategy to follow development of the IgM repertoire in health and disease provides a novel molecular assessment of multiple points along the B cell differentiation pathway that is highly sensitive for detecting perturbations within the repertoire at the population level.
Collapse
Affiliation(s)
- Li Yin
- Department of Pathology, Immunology and Laboratory Medicine, College of Medicine, University of Florida , Gainesville, FL , USA ; Florida Center for AIDS Research, University of Florida , Gainesville, FL , USA
| | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
225
|
Abstract
The evolutionary emergence of vertebrates was accompanied by major morphological and functional innovations, including the development of an adaptive immune system. Vertebrate adaptive immunity is based on the clonal expression of somatically diversifying antigen receptors on lymphocytes. This is a common feature of both the jawless and jawed vertebrates , although these two groups of extant vertebrates employ structurally different types of antigen receptors and principal mechanisms for their somatic diversification . These observations suggest that the common vertebrate ancestor must have already possessed a complex immune system, including B- and T-like lymphocyte lineages and primary lymphoid organs, such as the thymus, but possibly lacked the facilities for somatic diversification of antigen receptors. Interestingly, memory formation, previously considered to be a defining feature of adaptive immunity, also occurs in the context of innate immune responses and can even be observed in unicellular organisms, attesting to the convergent evolutionary history of distinct aspects of adaptive immunity.
Collapse
Affiliation(s)
- Thomas Boehm
- Max Planck Institute of Immunobiology and Epigenetics, 79108 Freiburg, Germany; ,
| | | |
Collapse
|
226
|
IgG variable region and VH CDR3 diversity in unimmunized mice analyzed by massively parallel sequencing. Mol Immunol 2013; 57:274-83. [PMID: 24211535 DOI: 10.1016/j.molimm.2013.09.008] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2013] [Revised: 09/11/2013] [Accepted: 09/28/2013] [Indexed: 11/24/2022]
Abstract
Most antigen-specific mouse antibodies have been derived by hybridoma technology, predominantly through use of the Balb/c strain. Much of the Balb/c germline repertoire of variable genes (V regions) is known. However, there is little information about the background expressed repertoire of IgG antibodies in mice, which reflects the baseline against which antigen-specific antibodies are generated through immunization. To assess this baseline repertoire, RNA was isolated from splenic B-cells enriched for expression of IgG from three mice. The RNA was individually amplified with three distinct PCR primer sets for comprehensive recovery of the heavy and light chain variable regions. Each PCR product was independently subjected to deep sequencing using 454 pyro-sequencing technology and analysed for redundancy, open reading frame, germline representation, and CDR3 sequence of the heavy chain variable region (VH CDR3) within and across the primer sets and mice. A highly skewed abundance of heavy and light chain variable gene usage was observed for all three primers in all three mice. While showing considerable overlap, there were differences among these profiles indicative of primer bias and animal-to-animal variation. VH CDR3 sequences were likewise highly skewed indicating that the heavy chain genes profiles substantially reflected individual antibodies. This observation was confirmed through analysis of randomly selected complete heavy chain variable sequences. However, there was very little redundancy in VH CDR3 sequences across the different mice. We conclude that the background IgG repertoire in young, unimmunized mice is highly skewed within individual mice and is diverse among them, a pattern similar to that observed in highly immunized mice.
Collapse
|
227
|
Busse CE, Czogiel I, Braun P, Arndt PF, Wardemann H. Single-cell based high-throughput sequencing of full-length immunoglobulin heavy and light chain genes. Eur J Immunol 2013; 44:597-603. [PMID: 24114719 DOI: 10.1002/eji.201343917] [Citation(s) in RCA: 92] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2013] [Revised: 08/27/2013] [Accepted: 09/19/2013] [Indexed: 11/09/2022]
Abstract
Single-cell PCR and sequencing of full-length Ig heavy (Igh) and Igk and Igl light chain genes is a powerful tool to measure the diversity of antibody repertoires and allows the functional assessment of B-cell responses through direct Ig gene cloning and the generation of recombinant mAbs. However, the current methodology is not high-throughput compatible. Here we developed a two-dimensional bar-coded primer matrix to combine Igh and Igk/Igl chain gene single-cell PCR with next-generation sequencing for the parallel analysis of the antibody repertoire of over 46 000 individual B cells. Our approach provides full-length Igh and corresponding Igk/Igl chain gene-sequence information and permits the accurate correction of sequencing errors by consensus building. The use of indexed cell sorting for the isolation of single B cells enables the integration of flow cytometry and Ig gene sequence information. The strategy is fully compatible with established protocols for direct antibody gene cloning and expression and therefore advances over previously described high-throughput approaches to assess antibody repertoires at the single-cell level.
Collapse
Affiliation(s)
- Christian E Busse
- Research Group Molecular Immunology, Max Planck Institute for Infection Biology, Berlin, Germany
| | | | | | | | | |
Collapse
|
228
|
Genomic and systems approaches to translational biomarker discovery in immunological diseases. Drug Discov Today 2013; 19:133-9. [PMID: 24126144 DOI: 10.1016/j.drudis.2013.10.002] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2013] [Revised: 09/13/2013] [Accepted: 10/04/2013] [Indexed: 02/07/2023]
Abstract
The high failure rate of new therapeutic mechanisms tested in clinical development has spurred an upsurge in research dedicated to discovering biomarker readouts that can improve decision-making. Increasingly, systems biology and genomic technologies, such as transcriptional profiling, are being leveraged to aid in the discovery of biomarker readouts. For inflammatory and immunological diseases, such as rheumatoid arthritis (RA) and asthma, progress has been made in developing biomarkers to monitor disease activity, prediction of response to therapy, and pharmacodynamic (PD) measurements. In this review, we discuss recent successes and challenges in these endeavors, highlighting the importance of human clinical studies of standard-of-care treatments in control subjects and patients with disease as the most direct path toward identifying useful translational biomarkers for clinical development.
Collapse
|
229
|
Robins H. Immunosequencing: applications of immune repertoire deep sequencing. Curr Opin Immunol 2013; 25:646-52. [DOI: 10.1016/j.coi.2013.09.017] [Citation(s) in RCA: 129] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2013] [Revised: 10/01/2013] [Accepted: 09/30/2013] [Indexed: 01/25/2023]
|
230
|
Colombo M, Cutrona G, Reverberi D, Bruno S, Ghiotto F, Tenca C, Stamatopoulos K, Hadzidimitriou A, Ceccarelli J, Salvi S, Boccardo S, Calevo MG, De Santanna A, Truini M, Fais F, Ferrarini M. Expression of immunoglobulin receptors with distinctive features indicating antigen selection by marginal zone B cells from human spleen. Mol Med 2013; 19:294-302. [PMID: 23877718 PMCID: PMC4344459 DOI: 10.2119/molmed.2013.00069] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2013] [Accepted: 07/11/2013] [Indexed: 01/21/2023] Open
Abstract
Marginal zone (MZ) B cells, identified as surface (s)IgM(high)sIgD(low)CD23(low/-)CD21(+)CD38(-) B cells, were purified from human spleens, and the features of their V(D)J gene rearrangements were investigated and compared with those of germinal center (GC), follicular mantle (FM) and switched memory (SM) B cells. Most MZ B cells were CD27(+) and exhibited somatic hypermutations (SHM), although to a lower extent than SM B cells. Moreover, among MZ B-cell rearrangements, recurrent sequences were observed, some of which displayed intraclonal diversification. The same diversifying sequences were detected in very low numbers in GC and FM B cells and only when a highly sensitive, gene-specific polymerase chain reaction was used. This result indicates that MZ B cells could expand and diversify in situ and also suggested the presence of a number of activation-induced cytidine deaminase (AID)-expressing B cells in the MZ. The notion of antigen-driven expansion/selection in situ is further supported by the VH CDR3 features of MZ B cells with highly conserved amino acids at specific positions and by the finding of shared ("stereotyped") sequences in two different spleens. Collectively, the data are consistent with the notion that MZ B cells are a special subset selected by in situ antigenic stimuli.
Collapse
Affiliation(s)
- Monica Colombo
- Direzione Scientifica, Istituto di Ricovero e Cura a Carattere Scientifico (IRCCS) Azienda Ospedaliera Universitaria (AOU) San Martino-IST, Genoa, Italy
| | - Giovanna Cutrona
- Division of Molecular Diagnostics, IRCCS AOU San Martino-IST, Genoa, Italy
| | | | - Silvia Bruno
- Department Experimental Medicine, University of Genoa, Genoa, Italy
| | - Fabio Ghiotto
- Department Experimental Medicine, University of Genoa, Genoa, Italy
| | - Claudya Tenca
- Department Experimental Medicine, University of Genoa, Genoa, Italy
| | - Kostas Stamatopoulos
- Institute of Agrobiotechnology, Center for Research and Technology Hellas, Thessaloniki, Greece
- Hematology Department and Hematopoietic Cell Transplantation (HCT) Unit, G. Papanicolau Hospital, Thessaloniki, Greece
| | | | - Jenny Ceccarelli
- Direzione Scientifica, Istituto di Ricovero e Cura a Carattere Scientifico (IRCCS) Azienda Ospedaliera Universitaria (AOU) San Martino-IST, Genoa, Italy
| | - Sandra Salvi
- Division of Pathology, IRCCS AOU San Martino-IST, Genoa, Italy
| | - Simona Boccardo
- Division of Pathology, IRCCS AOU San Martino-IST, Genoa, Italy
| | - Maria Grazia Calevo
- Division of Epidemiology and Statistic, IRCCS Istituto G. Gaslini, Genoa, Italy
| | | | - Mauro Truini
- Division of Pathology, IRCCS AOU San Martino-IST, Genoa, Italy
| | - Franco Fais
- Department Experimental Medicine, University of Genoa, Genoa, Italy
| | - Manlio Ferrarini
- Direzione Scientifica, Istituto di Ricovero e Cura a Carattere Scientifico (IRCCS) Azienda Ospedaliera Universitaria (AOU) San Martino-IST, Genoa, Italy
| |
Collapse
|
231
|
Mathonet P, Ullman CG. The application of next generation sequencing to the understanding of antibody repertoires. Front Immunol 2013; 4:265. [PMID: 24062737 PMCID: PMC3769625 DOI: 10.3389/fimmu.2013.00265] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2013] [Accepted: 08/19/2013] [Indexed: 01/01/2023] Open
Abstract
In the decade since the human genome sequence was declared complete, the development of next generation sequencing (NGS) or "deep" sequencing to deliver cost-effective genomic sequencing has influenced advances beyond its primary application and changed the research landscape in many other areas. This review will survey recent applications of NGS which have broadened the understanding of natural antibody repertoires (the "antibodyome") and how these evolve in response to viral infection. We will also report examples where deep sequencing of binding populations, derived from both natural and synthetic repertoires, have been used to benefit antibody engineering. This knowledge will ultimately lead to the design of more effective biological drugs and vaccines.
Collapse
|
232
|
Lakhani KR, Boudreau KJ, Loh PR, Backstrom L, Baldwin C, Lonstein E, Lydon M, MacCormack A, Arnaout RA, Guinan EC. Prize-based contests can provide solutions to computational biology problems. Nat Biotechnol 2013; 31:108-11. [PMID: 23392504 DOI: 10.1038/nbt.2495] [Citation(s) in RCA: 90] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
233
|
Michaeli M, Barak M, Hazanov L, Noga H, Mehr R. Automated analysis of immunoglobulin genes from high-throughput sequencing: life without a template. J Clin Bioinforma 2013; 3:15. [PMID: 23977981 PMCID: PMC3846367 DOI: 10.1186/2043-9113-3-15] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2013] [Accepted: 08/22/2013] [Indexed: 12/20/2022] Open
Abstract
Background Immunoglobulin (that is, antibody) and T cell receptor genes are created through somatic gene rearrangement from gene segment libraries. Immunoglobulin genes are further diversified by somatic hypermutation and selection during the immune response. Studying the repertoires of these genes yields valuable insights into immune system function in infections, aging, autoimmune diseases and cancers. The introduction of high throughput sequencing has generated unprecedented amounts of repertoire and mutation data from immunoglobulin genes. However, common analysis programs are not appropriate for pre-processing and analyzing these data due to the lack of a template or reference for the whole gene. Results We present here the automated analysis pipeline we created for this purpose, which integrates various software packages of our own development and others’, and demonstrate its performance. Conclusions Our analysis pipeline presented here is highly modular, and makes it possible to analyze the data resulting from high-throughput sequencing of immunoglobulin genes, in spite of the lack of a template gene. An executable version of the Automation program (and its source code) is freely available for downloading from our website: http://immsilico2.lnx.biu.ac.il/Software.html.
Collapse
Affiliation(s)
- Miri Michaeli
- The Mina & Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat-Gan 52900, Israel.
| | | | | | | | | |
Collapse
|
234
|
Singh V, Stoop MP, Stingl C, Luitwieler RL, Dekker LJ, van Duijn MM, Kreft KL, Luider TM, Hintzen RQ. Cerebrospinal-fluid-derived immunoglobulin G of different multiple sclerosis patients shares mutated sequences in complementarity determining regions. Mol Cell Proteomics 2013; 12:3924-34. [PMID: 23970564 DOI: 10.1074/mcp.m113.030346] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
B lymphocytes play a pivotal role in multiple sclerosis pathology, possibly via both antibody-dependent and -independent pathways. Intrathecal immunoglobulin G in multiple sclerosis is produced by clonally expanded B-cell populations. Recent studies indicate that the complementarity determining regions of immunoglobulins specific for certain antigens are frequently shared between different individuals. In this study, our main objective was to identify specific proteomic profiles of mutated complementarity determining regions of immunoglobulin G present in multiple sclerosis patients but absent in healthy controls. To achieve this objective, we purified immunoglobulin G from the cerebrospinal fluid of 29 multiple sclerosis patients and 30 healthy controls and separated the corresponding heavy and light chains via SDS-PAGE. Subsequently, bands were excised, trypsinized, and measured with high-resolution mass spectrometry. We sequenced 841 heavy and 771 light chain variable region peptides. We observed 24 heavy and 26 light chain complementarity determining regions that were solely present in a number of multiple sclerosis patients. Using stringent criteria for the identification of common peptides, we found five complementarity determining regions shared in three or more patients and not in controls. Interestingly, one complementarity determining region with a single mutation was found in six patients. Additionally, one other patient carrying a similar complementarity determining region with another mutation was observed. In addition, we found a skew in the κ-to-λ ratio and in the usage of certain variable heavy regions that was previously observed at the transcriptome level. At the protein level, cerebrospinal fluid immunoglobulin G shares common characteristics in the antigen binding region among different multiple sclerosis patients. The indication of a shared fingerprint may indicate common antigens for B-cell activation.
Collapse
Affiliation(s)
- Vaibhav Singh
- Department of Neurology, Erasmus University Medical Center, Rotterdam CA 3000, The Netherlands
| | | | | | | | | | | | | | | | | |
Collapse
|
235
|
Genetic measurement of memory B-cell recall using antibody repertoire sequencing. Proc Natl Acad Sci U S A 2013; 110:13463-8. [PMID: 23898164 DOI: 10.1073/pnas.1312146110] [Citation(s) in RCA: 176] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Annual influenza vaccinations aim to protect against seasonal infections, and vaccine strain compositions are updated every year. This protection is based on antibodies that are produced by either newly activated or memory B cells recalled from previous encounters with influenza vaccination or infection. The extent to which the B-cell repertoire responds to vaccination and recalls antibodies has so far not been analyzed at a genetic level--which is to say, at the level of antibody sequences. Here, we developed a consensus read sequencing approach that incorporates unique barcode labels on each starting RNA molecule. These labels allow one to combine multiple sequencing reads covering the same RNA molecule to reduce the error rate to a desired level, and they also enable accurate quantification of RNA and isotype levels. We validated this approach and analyzed the differential response of the antibody repertoire to live-attenuated or trivalent-inactivated influenza vaccination. Additionally, we analyzed the antibody repertoire in response to repeated yearly vaccinations with trivalent-inactivated influenza vaccination. We found antibody sequences that were present in both years, providing a direct genetic measurement of B-cell recall.
Collapse
|
236
|
Network properties derived from deep sequencing of human B-cell receptor repertoires delineate B-cell populations. Genome Res 2013; 23:1874-84. [PMID: 23742949 PMCID: PMC3814887 DOI: 10.1101/gr.154815.113] [Citation(s) in RCA: 100] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
The adaptive immune response selectively expands B- and T-cell clones following antigen recognition by B- and T-cell receptors (BCR and TCR), respectively. Next-generation sequencing is a powerful tool for dissecting the BCR and TCR populations at high resolution, but robust computational analyses are required to interpret such sequencing. Here, we develop a novel computational approach for BCR repertoire analysis using established next-generation sequencing methods coupled with network construction and population analysis. BCR sequences organize into networks based on sequence diversity, with differences in network connectivity clearly distinguishing between diverse repertoires of healthy individuals and clonally expanded repertoires from individuals with chronic lymphocytic leukemia (CLL) and other clonal blood disorders. Network population measures defined by the Gini Index and cluster sizes quantify the BCR clonality status and are robust to sampling and sequencing depths. BCR network analysis therefore allows the direct and quantifiable comparison of BCR repertoires between samples and intra-individual population changes between temporal or spatially separated samples and over the course of therapy.
Collapse
|
237
|
Forsyth CM, Juan V, Akamatsu Y, DuBridge RB, Doan M, Ivanov AV, Ma Z, Polakoff D, Razo J, Wilson K, Powers DB. Deep mutational scanning of an antibody against epidermal growth factor receptor using mammalian cell display and massively parallel pyrosequencing. MAbs 2013; 5:523-32. [PMID: 23765106 PMCID: PMC3906306 DOI: 10.4161/mabs.24979] [Citation(s) in RCA: 54] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
We developed a method for deep mutational scanning of antibody complementarity-determining regions (CDRs) that can determine in parallel the effect of every possible single amino acid CDR substitution on antigen binding. The method uses libraries of full length IgGs containing more than 1000 CDR point mutations displayed on mammalian cells, sorted by flow cytometry into subpopulations based on antigen affinity and analyzed by massively parallel pyrosequencing. Higher, lower and neutral affinity mutations are identified by their enrichment or depletion in the FACS subpopulations. We applied this method to a humanized version of the anti-epidermal growth factor receptor antibody cetuximab, generated a near comprehensive data set for 1060 point mutations that recapitulates previously determined structural and mutational data for these CDRs and identified 67 point mutations that increase affinity. The large-scale, comprehensive sequence-function data sets generated by this method should have broad utility for engineering properties such as antibody affinity and specificity and may advance theoretical understanding of antibody-antigen recognition.
Collapse
|
238
|
Deep mRNA sequencing analysis to capture the transcriptome landscape of zebrafish embryos and larvae. PLoS One 2013; 8:e64058. [PMID: 23700457 PMCID: PMC3659048 DOI: 10.1371/journal.pone.0064058] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2012] [Accepted: 04/09/2013] [Indexed: 11/19/2022] Open
Abstract
Transcriptome analysis is a powerful tool to obtain large amount genome-scale gene expression profiles. Despite its extensive usage to diverse biological problems in the last decade, transcriptomic researches approaching the zebrafish embryonic development have been very limited. Several recent studies have made great progress in this direction, yet the large gap still exists, especially regarding to the transcriptome dynamics of embryonic stages from early gastrulation onwards. Here, we present a comprehensive analysis about the transcriptomes of 9 different stages covering 7 major periods (cleavage, blastula, gastrula, segmentation, pharyngula, hatching and early larval stage) in zebrafish development, by recruiting the RNA-sequencing technology. We detected the expression for at least 24,065 genes in at least one of the 9 stages. We identified 16,130 genes that were significantly differentially expressed between stages and were subsequently classified into six clusters. Each revealed gene cluster had distinct expression patterns and characteristic functional pathways, providing a framework for the understanding of the developmental transcriptome dynamics. Over 4000 genes were identified as preferentially expressed in one of the stages, which could be of high relevance to stage-specific developmental and molecular events. Among the 68 transcription factor families active during development, most had enhanced average expression levels and thus might be crucial for embryogenesis, whereas the inactivation of the other families was likely required by the activation of the zygotic genome. We discussed our RNA-seq data together with previous findings about the Wnt signaling pathway and some other genes with known functions, to show how our data could be used to advance our understanding about these developmental functional elements. Our study provides ample information for further study about the molecular and cellular mechanisms underlying vertebrate development.
Collapse
|
239
|
Mazutis L, Gilbert J, Ung WL, Weitz DA, Griffiths AD, Heyman JA. Single-cell analysis and sorting using droplet-based microfluidics. Nat Protoc 2013; 8:870-91. [PMID: 23558786 PMCID: PMC4128248 DOI: 10.1038/nprot.2013.046] [Citation(s) in RCA: 825] [Impact Index Per Article: 75.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
We present a droplet-based microfluidics protocol for high-throughput analysis and sorting of single cells. Compartmentalization of single cells in droplets enables the analysis of proteins released from or secreted by cells, thereby overcoming one of the major limitations of traditional flow cytometry and fluorescence-activated cell sorting. As an example of this approach, we detail a binding assay for detecting antibodies secreted from single mouse hybridoma cells. Secreted antibodies are detected after only 15 min by co-compartmentalizing single mouse hybridoma cells, a fluorescent probe and single beads coated with anti-mouse IgG antibodies in 50-pl droplets. The beads capture the secreted antibodies and, when the captured antibodies bind to the probe, the fluorescence becomes localized on the beads, generating a clearly distinguishable fluorescence signal that enables droplet sorting at ∼200 Hz as well as cell enrichment. The microfluidic system described is easily adapted for screening other intracellular, cell-surface or secreted proteins and for quantifying catalytic or regulatory activities. In order to screen ∼1 million cells, the microfluidic operations require 2-6 h; the entire process, including preparation of microfluidic devices and mammalian cells, requires 5-7 d.
Collapse
Affiliation(s)
- Linas Mazutis
- School of Engineering and Applied Sciences (SEAS), Harvard University, Cambridge, Massachusetts, USA
| | | | | | | | | | | |
Collapse
|
240
|
Doorenspleet ME, Klarenbeek PL, de Hair MJH, van Schaik BDC, Esveldt REE, van Kampen AHC, Gerlag DM, Musters A, Baas F, Tak PP, de Vries N. Rheumatoid arthritis synovial tissue harbours dominant B-cell and plasma-cell clones associated with autoreactivity. Ann Rheum Dis 2013; 73:756-62. [PMID: 23606709 DOI: 10.1136/annrheumdis-2012-202861] [Citation(s) in RCA: 87] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
OBJECTIVE To identify potential autoreactive B-cell and plasma-cell clones by quantitatively analysing the complete human B-cell receptor (BCR) repertoire in synovium and peripheral blood in early and established rheumatoid arthritis (RA). METHODS The BCR repertoire was screened in synovium and blood of six patients with early RA (ERA) (<6 months) and six with established RA (ESRA) (>20 months). In two patients, the repertoires in different joints were compared. Repertoires were analysed by next-generation sequencing from mRNA, generating >10 000 BCR heavy-chain sequence reads per sample. For each clone, the degree of expansion was calculated as the percentage of the total number of reads encoding the specific clonal sequence. Clones with a frequency ≥ 0.5% were considered dominant. RESULTS Multiple dominant clones were found in inflamed synovium but hardly any in blood. Within an individual patient, the same dominant clones were detected in different joints. The majority of the synovial clones were class-switched; however, the fraction of clones that expressed IgM was higher in ESRA than ERA patients. Dominant synovial clones showed autoreactive features: in ERA in particular the clones were enriched for immunoglobulin heavy chain gene segment V4-34 (IGHV4-34) and showed longer CDR3 lengths. Dominant synovial clones that did not encode IGHV4-34 also had longer CDR3s than peripheral blood. CONCLUSIONS In RA, the synovium forms a niche where expanded--potentially autoreactive--B cells and plasma cells reside. The inflamed target tissue, especially in the earliest phase of disease, seems to be the most promising compartment for studying autoreactive cells.
Collapse
Affiliation(s)
- M E Doorenspleet
- Department of Clinical Immunology & Rheumatology, Academic Medical Center of the University of Amsterdam, , Amsterdam, The Netherlands
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
241
|
Baum PD, Venturi V, Price DA. Wrestling with the repertoire: the promise and perils of next generation sequencing for antigen receptors. Eur J Immunol 2013; 42:2834-9. [PMID: 23108932 DOI: 10.1002/eji.201242999] [Citation(s) in RCA: 41] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Next generation sequencing technologies are revolutionizing the study of immune repertoires. These methods provide a previously unimaginable amount of sequence data, unfortunately accompanied by numerous challenges associated with error correction and interpretation that remain to be solved. For antigen receptors, these challenges will require dedicated solutions beyond those developed for genome sequencing, which may differ depending on the sequencing technology used and the purpose of the experiment. Many investigators are developing such methods, based on different sequencing platforms, but critical details of protocol and performance are proprietary. The field will move forward when these methods are shared and standardized, and when the accuracy, sensitivity and reproducibility of various sequencing, and analytic methods are evaluated using standardized samples in comparative experiments.
Collapse
Affiliation(s)
- Paul D Baum
- Department of Medicine, University of California, San Francisco, CA 94114, USA.
| | | | | |
Collapse
|
242
|
Mining the antibodyome for HIV-1-neutralizing antibodies with next-generation sequencing and phylogenetic pairing of heavy/light chains. Proc Natl Acad Sci U S A 2013; 110:6470-5. [PMID: 23536288 DOI: 10.1073/pnas.1219320110] [Citation(s) in RCA: 127] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Next-generation sequencing of antibody transcripts from HIV-1-infected individuals with broadly neutralizing antibodies could provide an efficient means for identifying somatic variants and characterizing their lineages. Here, we used 454 pyrosequencing and identity/divergence grid sampling to analyze heavy- and light-chain sequences from donor N152, the source of the broadly neutralizing antibody 10E8. We identified variants with up to 28% difference in amino acid sequence. Heavy- and light-chain phylogenetic trees of identified 10E8 variants displayed similar architectures, and 10E8 variants reconstituted from matched and unmatched phylogenetic branches displayed significantly lower autoreactivity when matched. To test the generality of phylogenetic pairing, we analyzed donor International AIDS Vaccine Initiative 84, the source of antibodies PGT141-145. Heavy- and light-chain phylogenetic trees of PGT141-145 somatic variants also displayed remarkably similar architectures; in this case, branch pairings could be anchored by known PGT141-145 antibodies. Altogether, our findings suggest that phylogenetic matching of heavy and light chains can provide a means to approximate natural pairings.
Collapse
|
243
|
Ravn U, Didelot G, Venet S, Ng KT, Gueneau F, Rousseau F, Calloud S, Kosco-Vilbois M, Fischer N. Deep sequencing of phage display libraries to support antibody discovery. Methods 2013; 60:99-110. [PMID: 23500657 DOI: 10.1016/j.ymeth.2013.03.001] [Citation(s) in RCA: 87] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2012] [Revised: 02/18/2013] [Accepted: 03/01/2013] [Indexed: 11/18/2022] Open
Abstract
The use of next generation sequencing (NGS) for the analysis of antibody sequences both in phage display libraries and during in vitro selection processes has become increasingly popular in the last few years. Here, our methods developed for DNA preparation, sequencing and data analysis are presented. A key parameter has also been to develop new software designed for high throughput antibody sequence analysis that is used in combination with publicly available tools. As an example of our methods, we provide data from the extensive analysis of five scFv libraries generated using different heavy chain CDR3 diversification strategies. The results not only confirm that the library designs were correct but also reveal differences in quality not easily identified by standard DNA sequencing approaches. The very large number of reads permits extensive sequence coverage after the selection process. Furthermore, as samples can be multiplexed, costs decrease and more information is gained per NGS run. Using examples of results obtained post phage display selections against two antigens, frequency and clustering analysis identified novel antibody fragments that were then shown to be specific for the target antigen. In summary, the methods described here demonstrate how NGS analysis enhances quality control of complex antibody libraries as well as facilitates the antibody discovery process.
Collapse
Affiliation(s)
- Ulla Ravn
- NovImmune S.A., 14 Chemin des Aulx, CH-1228 Plan-les-Ouates, Switzerland
| | | | | | | | | | | | | | | | | |
Collapse
|
244
|
Birnbaum ME, Dong S, Garcia KC. Diversity-oriented approaches for interrogating T-cell receptor repertoire, ligand recognition, and function. Immunol Rev 2013; 250:82-101. [PMID: 23046124 DOI: 10.1111/imr.12006] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Molecular diversity lies at the heart of adaptive immunity. T-cell receptors and peptide-major histocompatibility complex molecules utilize and rely upon an enormous degree of diversity at the levels of genetics, chemistry, and structure to engage one another and carry out their functions. This high level of diversity complicates the systematic study of important aspects of T-cell biology, but recent technical advances have allowed for the ability to study diversity in a comprehensive manner. In this review, we assess insights gained into T-cell receptor function and biology from our increasingly precise ability to assess the T-cell repertoire as a whole or to perturb individual receptors with engineered reagents. We conclude with a perspective on a new class of high-affinity, non-stimulatory peptide ligands we have recently discovered using diversity-oriented techniques that challenges notions for how we think about T-cell receptor signaling.
Collapse
Affiliation(s)
- Michael E Birnbaum
- Department of Molecular and Cellular Physiology, Program in Immunology, Stanford University School of Medicine, CA, USA
| | | | | |
Collapse
|
245
|
Fillatreau S, Six A, Magadan S, Castro R, Sunyer JO, Boudinot P. The astonishing diversity of Ig classes and B cell repertoires in teleost fish. Front Immunol 2013; 4:28. [PMID: 23408183 PMCID: PMC3570791 DOI: 10.3389/fimmu.2013.00028] [Citation(s) in RCA: 109] [Impact Index Per Article: 9.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2012] [Accepted: 01/24/2013] [Indexed: 12/17/2022] Open
Abstract
With lymphoid tissue anatomy different than mammals, and diverse adaptations to all aquatic environments, fish constitute a fascinating group of vertebrate to study the biology of B cell repertoires in a comparative perspective. Fish B lymphocytes express immunoglobulin (Ig) on their surface and secrete antigen-specific antibodies in response to immune challenges. Three antibody classes have been identified in fish, namely IgM, IgD, and IgT, while IgG, IgA, and IgE are absent. IgM and IgD have been found in all fish species analyzed, and thus seem to be primordial antibody classes. IgM and IgD are normally co-expressed from the same mRNA through alternative splicing, as in mammals. Tetrameric IgM is the main antibody class found in serum. Some species of fish also have IgT, which seems to exist only in fish and is specialized in mucosal immunity. IgM/IgD and IgT are expressed by two different sub-populations of B cells. The tools available to investigate B cell responses at the cellular level in fish are limited, but the progress of fish genomics has started to unravel a rich diversity of IgH and immunoglobulin light chain locus organization, which might be related to the succession of genome remodelings that occurred during fish evolution. Moreover, the development of deep sequencing techniques has allowed the investigation of the global features of the expressed fish B cell repertoires in zebrafish and rainbow trout, in steady state or after infection. This review provides a description of the organization of fish Ig loci, with a particular emphasis on their heterogeneity between species, and presents recent data on the structure of the expressed Ig repertoire in healthy and infected fish.
Collapse
Affiliation(s)
- Simon Fillatreau
- Deutsches Rheuma-Forschungszentrum, Leibniz Institute Berlin, Germany
| | | | | | | | | | | |
Collapse
|
246
|
DeKosky BJ, Ippolito GC, Deschner RP, Lavinder JJ, Wine Y, Rawlings BM, Varadarajan N, Giesecke C, Dörner T, Andrews SF, Wilson PC, Hunicke-Smith SP, Willson CG, Ellington AD, Georgiou G. High-throughput sequencing of the paired human immunoglobulin heavy and light chain repertoire. Nat Biotechnol 2013; 31:166-9. [PMID: 23334449 DOI: 10.1038/nbt.2492] [Citation(s) in RCA: 327] [Impact Index Per Article: 29.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2012] [Accepted: 01/02/2013] [Indexed: 02/07/2023]
Abstract
Each B-cell receptor consists of a pair of heavy and light chains. High-throughput sequencing can identify large numbers of heavy- and light-chain variable regions (V(H) and V(L)) in a given B-cell repertoire, but information about endogenous pairing of heavy and light chains is lost after bulk lysis of B-cell populations. Here we describe a way to retain this pairing information. In our approach, single B cells (>5 × 10(4) capacity per experiment) are deposited in a high-density microwell plate (125 pl/well) and lysed in situ. mRNA is then captured on magnetic beads, reverse transcribed and amplified by emulsion V(H):V(L) linkage PCR. The linked transcripts are analyzed by Illumina high-throughput sequencing. We validated the fidelity of V(H):V(L) pairs identified by this approach and used the method to sequence the repertoire of three human cell subsets-peripheral blood IgG(+) B cells, peripheral plasmablasts isolated after tetanus toxoid immunization and memory B cells isolated after seasonal influenza vaccination.
Collapse
Affiliation(s)
- Brandon J DeKosky
- Department of Chemical Engineering, University of Texas at Austin, Austin, Texas, USA
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
247
|
Meier J, Roberts C, Avent K, Hazlett A, Berrie J, Payne K, Hamm D, Desmarais C, Sanders C, Hogan KT, Archer KJ, Manjili MH, Toor AA. Fractal organization of the human T cell repertoire in health and after stem cell transplantation. Biol Blood Marrow Transplant 2013; 19:366-77. [PMID: 23313705 DOI: 10.1016/j.bbmt.2012.12.004] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2012] [Accepted: 12/12/2012] [Indexed: 11/18/2022]
Abstract
T cell repertoire diversity is generated in part by recombination of variable (V), diversity (D), and joining (J) segments in the T cell receptor β (TCR) locus. T cell clonal frequency distribution determined by high-throughput sequencing of TCR β in 10 stem cell transplantation (SCT) donors revealed a fractal, self-similar frequency distribution of unique TCR bearing clones with respect to V, D, and J segment usage in the T cell repertoire of these individuals. Further, ranking of T cell clones by frequency of gene segment usage in the observed sequences revealed an ordered distribution of dominant clones conforming to a power law, with a fractal dimension of 1.6 and 1.8 in TCR β DJ and VDJ containing clones in healthy stem cell donors. This self-similar distribution was perturbed in the recipients after SCT, with patients demonstrating a lower level of complexity in their TCR repertoire at day 100 followed by a modest improvement by 1 year post-SCT. A large shift was observed in the frequency distribution of the dominant T cell clones compared to the donor, with fewer than one third of the VDJ-containing clones shared in the top 4 ranks. In conclusion, the normal T cell repertoire is highly ordered with a TCR gene segment usage that results in a fractal self-similar motif of pattern repetition across levels of organization. Fractal analysis of high-throughput TCR β sequencing data provides a comprehensive measure of immune reconstitution after SCT.
Collapse
MESH Headings
- Antilymphocyte Serum/pharmacology
- Antilymphocyte Serum/therapeutic use
- Clone Cells
- Fractals
- Hematologic Neoplasms/immunology
- Hematologic Neoplasms/pathology
- Hematologic Neoplasms/therapy
- High-Throughput Nucleotide Sequencing
- Humans
- Myeloablative Agonists/pharmacology
- Myeloablative Agonists/therapeutic use
- Receptors, Antigen, T-Cell, alpha-beta/chemistry
- Receptors, Antigen, T-Cell, alpha-beta/genetics
- Receptors, Antigen, T-Cell, alpha-beta/immunology
- Stem Cell Transplantation
- T-Lymphocytes/classification
- T-Lymphocytes/immunology
- T-Lymphocytes/pathology
- Transplantation Chimera/immunology
- Transplantation Conditioning
- Transplantation, Homologous
Collapse
Affiliation(s)
- Jeremy Meier
- Bone Marrow Transplant Program, Department of Internal Medicine, Virginia Commonwealth University, Richmond, VA 23298, USA
| | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
248
|
Castro R, Jouneau L, Pham HP, Bouchez O, Giudicelli V, Lefranc MP, Quillet E, Benmansour A, Cazals F, Six A, Fillatreau S, Sunyer O, Boudinot P. Teleost fish mount complex clonal IgM and IgT responses in spleen upon systemic viral infection. PLoS Pathog 2013; 9:e1003098. [PMID: 23326228 PMCID: PMC3542120 DOI: 10.1371/journal.ppat.1003098] [Citation(s) in RCA: 116] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2012] [Accepted: 11/09/2012] [Indexed: 01/12/2023] Open
Abstract
Upon infection, B-lymphocytes expressing antibodies specific for the intruding pathogen develop clonal responses triggered by pathogen recognition via the B-cell receptor. The constant region of antibodies produced by such responding clones dictates their functional properties. In teleost fish, the clonal structure of B-cell responses and the respective contribution of the three isotypes IgM, IgD and IgT remain unknown. The expression of IgM and IgT are mutually exclusive, leading to the existence of two B-cell subsets expressing either both IgM and IgD or only IgT. Here, we undertook a comprehensive analysis of the variable heavy chain (VH) domain repertoires of the IgM, IgD and IgT in spleen of homozygous isogenic rainbow trout (Onchorhynchus mykiss) before, and after challenge with a rhabdovirus, the Viral Hemorrhagic Septicemia Virus (VHSV), using CDR3-length spectratyping and pyrosequencing of immunoglobulin (Ig) transcripts. In healthy fish, we observed distinct repertoires for IgM, IgD and IgT, respectively, with a few amplified μ and τ junctions, suggesting the presence of IgM- and IgT-secreting cells in the spleen. In infected animals, we detected complex and highly diverse IgM responses involving all VH subgroups, and dominated by a few large public and private clones. A lower number of robust clonal responses involving only a few VH were detected for the mucosal IgT, indicating that both IgM(+) and IgT(+) spleen B cells responded to systemic infection but at different degrees. In contrast, the IgD response to the infection was faint. Although fish IgD and IgT present different structural features and evolutionary origin compared to mammalian IgD and IgA, respectively, their implication in the B-cell response evokes these mouse and human counterparts. Thus, it appears that the general properties of antibody responses were already in place in common ancestors of fish and mammals, and were globally conserved during evolution with possible functional convergences.
Collapse
Affiliation(s)
- Rosario Castro
- Virologie et Immunologie Moléculaires, INRA, Jouy-en-Josas, France
| | - Luc Jouneau
- Virologie et Immunologie Moléculaires, INRA, Jouy-en-Josas, France
| | - Hang-Phuong Pham
- UPMC Univ Paris 06, UMR 7211, “Integrative Immunology” Team, Paris, France; CNRS, UMR 7211, “Immunology, Immunopathology, Immunotherapy,” Paris, France
| | - Olivier Bouchez
- UMR INRA 0444 Laboratoire de Génétique Cellulaire, GeT-PlaGe Core Facility, Castanet Tolosan, France
| | - Véronique Giudicelli
- IMGT, the International ImMunoGeneTics Information System, Laboratoire d'ImmunoGénétique Moléculaire LIGM, IGH, UPR CNRS 1142 and Université Montpellier 2, Montpellier, France
| | - Marie-Paule Lefranc
- IMGT, the International ImMunoGeneTics Information System, Laboratoire d'ImmunoGénétique Moléculaire LIGM, IGH, UPR CNRS 1142 and Université Montpellier 2, Montpellier, France
| | - Edwige Quillet
- Génétique Animale et Biologie Intégrative, INRA, Jouy-en-Josas, France
| | | | - Frédéric Cazals
- INRIA Sophia-Antipolis - Méditerranée, Algorithms-Biology-Structure, Sophia-Antipolis, France
| | - Adrien Six
- UPMC Univ Paris 06, UMR 7211, “Integrative Immunology” Team, Paris, France; CNRS, UMR 7211, “Immunology, Immunopathology, Immunotherapy,” Paris, France
| | - Simon Fillatreau
- Deutsches RheumaForschungszentrum, a Leibniz Institute, Berlin, Germany
| | - Oriol Sunyer
- Department of Pathobiology, School of Veterinary Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
| | - Pierre Boudinot
- Virologie et Immunologie Moléculaires, INRA, Jouy-en-Josas, France
| |
Collapse
|
249
|
Yin L, Liu L, Sun Y, Hou W, Lowe AC, Gardner BP, Salemi M, Williams WB, Farmerie WG, Sleasman JW, Goodenow MM. High-resolution deep sequencing reveals biodiversity, population structure, and persistence of HIV-1 quasispecies within host ecosystems. Retrovirology 2012; 9:108. [PMID: 23244298 PMCID: PMC3531307 DOI: 10.1186/1742-4690-9-108] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2012] [Accepted: 11/20/2012] [Indexed: 02/05/2023] Open
Abstract
Background Deep sequencing provides the basis for analysis of biodiversity of taxonomically similar organisms in an environment. While extensively applied to microbiome studies, population genetics studies of viruses are limited. To define the scope of HIV-1 population biodiversity within infected individuals, a suite of phylogenetic and population genetic algorithms was applied to HIV-1 envelope hypervariable domain 3 (Env V3) within peripheral blood mononuclear cells from a group of perinatally HIV-1 subtype B infected, therapy-naïve children. Results Biodiversity of HIV-1 Env V3 quasispecies ranged from about 70 to 270 unique sequence clusters across individuals. Viral population structure was organized into a limited number of clusters that included the dominant variants combined with multiple clusters of low frequency variants. Next generation viral quasispecies evolved from low frequency variants at earlier time points through multiple non-synonymous changes in lineages within the evolutionary landscape. Minor V3 variants detected as long as four years after infection co-localized in phylogenetic reconstructions with early transmitting viruses or with subsequent plasma virus circulating two years later. Conclusions Deep sequencing defines HIV-1 population complexity and structure, reveals the ebb and flow of dominant and rare viral variants in the host ecosystem, and identifies an evolutionary record of low-frequency cell-associated viral V3 variants that persist for years. Bioinformatics pipeline developed for HIV-1 can be applied for biodiversity studies of virome populations in human, animal, or plant ecosystems.
Collapse
Affiliation(s)
- Li Yin
- Department of Pathology, Immunology and Laboratory Medicine, College of Medicine, University of Florida, 2033 Mowry Road, PO Box 103633, Gainesville, FL 32610-3633, USA.
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
250
|
Rubelt F, Sievert V, Knaust F, Diener C, Lim TS, Skriner K, Klipp E, Reinhardt R, Lehrach H, Konthur Z. Onset of immune senescence defined by unbiased pyrosequencing of human immunoglobulin mRNA repertoires. PLoS One 2012; 7:e49774. [PMID: 23226220 PMCID: PMC3511497 DOI: 10.1371/journal.pone.0049774] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2012] [Accepted: 10/12/2012] [Indexed: 12/15/2022] Open
Abstract
The immune system protects us from foreign substances or pathogens by generating specific antibodies. The variety of immunoglobulin (Ig) paratopes for antigen recognition is a result of the V(D)J rearrangement mechanism, while a fast and efficient immune response is mediated by specific immunoglobulin isotypes obtained through class switch recombination (CSR). To get a better understanding on how antibody-based immune protection works and how it changes with age, the interdependency between these two parameters need to be addressed. Here, we have performed an in depth analysis of antibody repertoires of 14 healthy donors representing different gender and age groups. For this task, we developed a unique pyrosequencing approach, which is able to monitor the expression levels of all immunoglobulin V(D)J recombinations of all isotypes including subtypes in an unbiased and quantitative manner. Our results show that donors have individual immunoglobulin repertoires and cannot be clustered according to V(D)J recombination patterns, neither by age nor gender. However, after incorporating isotype-specific analysis and considering CSR information into hierarchical clustering the situation changes. For the first time the donors cluster according to age and separate into young adults and elderly donors (>50). As a direct consequence, this clustering defines the onset of immune senescence at the age of fifty and beyond. The observed age-dependent reduction of CSR ability proposes a feasible explanation why reduced efficacy of vaccination is seen in the elderly and implies that novel vaccine strategies for the elderly should include the "Golden Agers".
Collapse
Affiliation(s)
- Florian Rubelt
- Max Planck Institute for Molecular Genetics, Berlin, Germany
- Faculty of Biology, Chemistry, and Pharmacy, Freie Universität Berlin, Berlin, Germany
| | - Volker Sievert
- Max Planck Institute for Molecular Genetics, Berlin, Germany
| | - Florian Knaust
- Max Planck Institute for Molecular Genetics, Berlin, Germany
| | - Christian Diener
- Max Planck Institute for Molecular Genetics, Berlin, Germany
- Theoretische Biophysik, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Theam Soon Lim
- Max Planck Institute for Molecular Genetics, Berlin, Germany
- Institute for Research in Molecular Medicine, Universiti Sains Malaysia, Penang, Malaysia
| | - Karl Skriner
- Department of Rheumatology and Clinical Immunology, Charité – Universitätsmedizin Berlin, Berlin, Germany
| | - Edda Klipp
- Theoretische Biophysik, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Richard Reinhardt
- Max Planck Institute for Molecular Genetics, Berlin, Germany
- Max Planck Genome Centre Cologne, Max Planck Institute for Plant Breeding Research, Cologne, Germany
| | - Hans Lehrach
- Max Planck Institute for Molecular Genetics, Berlin, Germany
| | - Zoltán Konthur
- Max Planck Institute for Molecular Genetics, Berlin, Germany
- * E-mail:
| |
Collapse
|