SONAR: A High-Throughput Pipeline for Inferring Antibody Ontogenies from Longitudinal Sequencing of B Cell Transcripts

The rapid advance of massively parallel or next-generation sequencing technologies has made possible the characterization of B cell receptor repertoires in ever greater detail, and these developments have triggered a proliferation of software tools for processing and annotating these data. Of especial interest, however, is the capability to track the development of specific antibody lineages across time, which remains beyond the scope of most current programs. We have previously reported on the use of techniques such as inter- and intradonor analysis and CDR3 tracing to identify transcripts related to an antibody of interest. Here, we present Software for the Ontogenic aNalysis of Antibody Repertoires (SONAR), capable of automating both general repertoire analysis and specialized techniques for investigating specific lineages. SONAR annotates next-generation sequencing data, identifies transcripts in a lineage of interest, and tracks lineage development across multiple time points. SONAR also generates figures, such as identity–divergence plots and longitudinal phylogenetic “birthday” trees, and provides interfaces to other programs such as DNAML and BEAST. SONAR can be downloaded as a ready-to-run Docker image or manually installed on a local machine. In the latter case, it can also be configured to take advantage of a high-performance computing cluster for the most computationally intensive steps, if available. In summary, this software provides a useful new tool for the processing of large next-generation sequencing datasets and the ontogenic analysis of neutralizing antibody lineages. SONAR can be found at https://github.com/scharch/SONAR, and the Docker image can be obtained from https://hub.docker.com/r/scharch/sonar/.

[1]  G. Yaari,et al.  Practical guidelines for B-cell receptor repertoire sequencing analysis , 2015, Genome Medicine.

[2]  S. Plotkin Correlates of Protection Induced by Vaccination , 2010, Clinical and Vaccine Immunology.

[3]  IV FrederickA.Matsen,et al.  Consistency of VDJ Rearrangement and Substitution Parameters Enables Accurate B Cell Receptor Sequence Annotation , 2015, PLoS Comput. Biol..

[4]  Thomas B Kepler,et al.  Reconstructing a B-cell clonal lineage. I. Statistical inference of unobserved ancestors , 2013, F1000Research.

[5]  Xun Xu,et al.  IMonitor: A Robust Pipeline for TCR and BCR Repertoire Analysis , 2015, Genetics.

[6]  C. Cheah,et al.  Breakthrough therapies in B-cell non-Hodgkin lymphoma. , 2016, Annals of oncology : official journal of the European Society for Medical Oncology.

[7]  Chaim A. Schramm,et al.  Co-evolution of a broadly neutralizing HIV-1 antibody and founder virus , 2013, Nature.

[8]  Hanlee P. Ji,et al.  Next-generation DNA sequencing , 2008, Nature Biotechnology.

[9]  Stephen R. Quake,et al.  Genetic measurement of memory B-cell recall using antibody repertoire sequencing , 2013, Proceedings of the National Academy of Sciences.

[10]  D. Higgins,et al.  Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega , 2011, Molecular systems biology.

[11]  Kwan-Yuet Ho,et al.  HTJoinSolver: Human immunoglobulin VDJ partitioning using approximate dynamic programming constrained by conserved motifs , 2015, BMC Bioinformatics.

[12]  Mark M. Davis,et al.  Lineage Structure of the Human Antibody Repertoire in Response to Influenza Vaccination , 2013, Science Translational Medicine.

[13]  Mario Roederer,et al.  Focused Evolution of HIV-1 Neutralizing Antibodies Revealed by Structures and Deep Sequencing , 2011, Science.

[14]  Eline T. Luning Prak,et al.  The analysis of clonal expansions in normal and autoimmune B cell repertoires , 2015, Philosophical Transactions of the Royal Society B: Biological Sciences.

[15]  Dong Xie,et al.  BEAST 2: A Software Platform for Bayesian Evolutionary Analysis , 2014, PLoS Comput. Biol..

[16]  Matthew R. Pocock,et al.  The Bioperl toolkit: Perl modules for the life sciences. , 2002, Genome research.

[17]  J. Felsenstein,et al.  A Hidden Markov Model approach to variation among sites in rate of evolution. , 1996, Molecular biology and evolution.

[18]  Dennis R. Burton,et al.  Clonify: unseeded antibody lineage assignment from next-generation sequencing data , 2016, Scientific Reports.

[19]  Young Do Kwon,et al.  Maturation and Diversity of the VRC01-Antibody Lineage over 15 Years of Chronic HIV-1 Infection , 2015, Cell.

[20]  L. Staudt,et al.  B-cell receptor signaling in diffuse large B-cell lymphoma. , 2015, Seminars in hematology.

[21]  Simon J. Henderson,et al.  Monoclonal antibody therapeutics: history and future. , 2012, Current opinion in pharmacology.

[22]  Robert C. Edgar,et al.  MUSCLE: multiple sequence alignment with high accuracy and high throughput. , 2004, Nucleic acids research.

[23]  Robert C. Edgar,et al.  BIOINFORMATICS APPLICATIONS NOTE , 2001 .

[24]  Mark M. Davis,et al.  Human responses to influenza vaccination show seroconversion signatures and convergent antibody rearrangements. , 2014, Cell host & microbe.

[25]  M. Metzker Sequencing technologies — the next generation , 2010, Nature Reviews Genetics.

[26]  Thomas B Kepler,et al.  Reconstructing a B-cell clonal lineage , 2016 .

[27]  Chaim A. Schramm,et al.  Maturation Pathway from Germline to Broad HIV-1 Neutralizer of a CD4-Mimic Antibody , 2016, Cell.

[28]  K. Katoh,et al.  MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability , 2013, Molecular biology and evolution.

[29]  Jonathan R. McDaniel,et al.  Structures of HIV-1-Env V1V2 with broadly neutralizing antibodies reveal commonalities that enable vaccine design , 2015, Nature Structural &Molecular Biology.

[30]  J. Mullikin,et al.  Somatic Populations of PGT135–137 HIV-1-Neutralizing Antibodies Identified by 454 Pyrosequencing and Bioinformatics , 2012, Front. Microbio..

[31]  P. Lipsky,et al.  Characterization of the Human Ig Heavy Chain Antigen Binding Complementarity Determining Region 3 Using a Newly Developed Software Algorithm, JOINSOLVER , 2004, The Journal of Immunology.

[32]  S. Quake,et al.  The promise and challenge of high-throughput sequencing of the antibody repertoire , 2014, Nature Biotechnology.

[33]  Joaquín Dopazo,et al.  ETE: a python Environment for Tree Exploration , 2010, BMC Bioinformatics.

[34]  Elizabeth Ernestina Godoy-Lozano,et al.  Reconstructing and mining the B cell repertoire with ImmunediveRsity , 2015, mAbs.

[35]  William D. Lees,et al.  Utilities for High-Throughput Analysis of B-Cell Clonal Lineages , 2015, Journal of immunology research.

[36]  M. Egholm,et al.  Measurement and Clinical Monitoring of Human Lymphocyte Clonality by Massively Parallel V-D-J Pyrosequencing , 2009, Science Translational Medicine.

[37]  Lynn Morris,et al.  Effects of Darwinian Selection and Mutability on Rate of Broadly Neutralizing Antibody Evolution during HIV-1 Infection , 2016, PLoS Comput. Biol..

[38]  M. Shlomchik,et al.  From T to B and back again: positive feedback in systemic autoimmune disease , 2001, Nature Reviews Immunology.

[39]  Baoshan Zhang,et al.  Mining the antibodyome for HIV-1–neutralizing antibodies with next-generation sequencing and phylogenetic pairing of heavy/light chains , 2013, Proceedings of the National Academy of Sciences.

[40]  Young Do Kwon,et al.  Enhanced Potency of a Broadly Neutralizing HIV-1 Antibody In Vitro Improves Protection against Lentiviral Infection In Vivo , 2014, Journal of Virology.

[41]  D. Koller,et al.  High-resolution antibody dynamics of vaccine-induced immune responses , 2014, Proceedings of the National Academy of Sciences.

[42]  Chaim A. Schramm,et al.  Developmental pathway for potent V1V2-directed HIV-neutralizing antibodies , 2014, Nature.

[43]  Bartek Wilczynski,et al.  Biopython: freely available Python tools for computational molecular biology and bioinformatics , 2009, Bioinform..

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

[45]  Ron Unger,et al.  IgTree: creating Immunoglobulin variable region gene lineage trees. , 2008, Journal of immunological methods.

[46]  J. Mascola,et al.  Safety, pharmacokinetics and neutralization of the broadly neutralizing HIV‐1 human monoclonal antibody VRC01 in healthy adults , 2015, Clinical and experimental immunology.

[47]  Joseph G. Jardine,et al.  HIV-1 broadly neutralizing antibody precursor B cells revealed by germline-targeting immunogen , 2016, Science.

[48]  Lynn Morris,et al.  New Member of the V1V2-Directed CAP256-VRC26 Lineage That Shows Increased Breadth and Exceptional Potency , 2015, Journal of Virology.

[49]  Ning Ma,et al.  BLAST+: architecture and applications , 2009, BMC Bioinformatics.

[50]  Tongqing Zhou,et al.  De novo identification of VRC01 class HIV-1–neutralizing antibodies by next-generation sequencing of B-cell transcripts , 2013, Proceedings of the National Academy of Sciences.

[51]  B. Haynes,et al.  Optimization of the Solubility of HIV-1-Neutralizing Antibody 10E8 through Somatic Variation and Structure-Based Design , 2016, Journal of Virology.

[52]  Enrico Macii,et al.  VDJSeq-Solver: In Silico V(D)J Recombination Detection Tool , 2015, PloS one.

[53]  Patrice Duroux,et al.  IMGT®, the international ImMunoGeneTics information system® 25 years on , 2014, Nucleic Acids Res..

[54]  Thomas B. Kepler,et al.  Reconstructing a B-Cell Clonal Lineage. II. Mutation, Selection, and Affinity Maturation , 2014, Front. Immunol..

[55]  Thierry Mora,et al.  The Past, Present, and Future of Immune Repertoire Biology – The Rise of Next-Generation Repertoire Analysis , 2013, Front. Immunol..

[56]  Lynn Morris,et al.  Viral variants that initiate and drive maturation of V1V2-directed HIV-1 broadly neutralizing antibodies , 2015, Nature Medicine.

[57]  K. P. Murphy,et al.  Janeway's immunobiology , 2007 .