The mutation patterns in B-cell immunoglobulin receptors reflect the influence of selection acting at multiple time-scales

During the several-week course of an immune response, B cells undergo a process of clonal expansion, somatic hypermutation of the immunoglobulin (Ig) genes and affinity-dependent selection. Over a lifetime, each B cell may participate in multiple rounds of affinity maturation as part of different immune responses. These two time-scales for selection are apparent in the structure of B-cell lineage trees, which often contain a ‘trunk’ consisting of mutations that are shared across all members of a clone, and several branches that form a ‘canopy’ consisting of mutations that are shared by a subset of clone members. The influence of affinity maturation on the B-cell population can be inferred by analysing the pattern of somatic mutations in the Ig. While global analysis of mutation patterns has shown evidence of strong selection pressures shaping the B-cell population, the effect of different time-scales of selection and diversification has not yet been studied. Analysis of B cells from blood samples of three healthy individuals identifies a range of clone sizes with lineage trees that can contain long trunks and canopies indicating the significant diversity introduced by the affinity maturation process. We here show that observed mutation patterns in the framework regions (FWRs) are determined by an almost purely purifying selection on both short and long time-scales. By contrast, complementarity determining regions (CDRs) are affected by a combination of purifying and antigen-driven positive selection on the short term, which leads to a net positive selection in the long term. In both the FWRs and CDRs, long-term selection is strongly dependent on the heavy chain variable gene family.

[1]  L. Wysocki,et al.  Di- and trinucleotide target preferences of somatic mutagenesis in normal and autoreactive B cells. , 1996, Journal of immunology.

[2]  Steven H. Kleinstein,et al.  Detecting selection in immunoglobulin sequences , 2011, Nucleic Acids Res..

[3]  D. Allen,et al.  Antibody engineering for the analysis of affinity maturation of an anti‐hapten response. , 1988, The EMBO journal.

[4]  L. Wysocki,et al.  Sequence-specific targeting of two bases on both DNA strands by the somatic hypermutation mechanism. , 2003, Molecular immunology.

[5]  Thomas L. Rothstein,et al.  The role of clonal selection and somatic mutation in autoimmunity , 1987, Nature.

[6]  Steven H. Kleinstein,et al.  Estimating Hypermutation Rates from Clonal Tree Data 1 , 2003, The Journal of Immunology.

[7]  D. Pisetsky,et al.  Structure and function of anti-DNA autoantibodies derived from a single autoimmune mouse. , 1987, Proceedings of the National Academy of Sciences of the United States of America.

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

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

[10]  Uri Hershberg,et al.  Taking Advantage: High-Affinity B Cells in the Germinal Center Have Lower Death Rates, but Similar Rates of Division, Compared to Low-Affinity Cells1 , 2009, The Journal of Immunology.

[11]  Steven H. Kleinstein,et al.  Quantifying selection in high-throughput Immunoglobulin sequencing data sets , 2012, Nucleic acids research.

[12]  Garnett Kelsoe,et al.  Predicted and inferred waiting times for key mutations in the germinal centre reaction: Evidence for stochasticity in selection , 1998, Immunology and cell biology.

[13]  U. Storb,et al.  The molecular basis of somatic hypermutation of immunoglobulin genes. , 1996, Current opinion in immunology.

[14]  Mohamed Uduman,et al.  Integrating B Cell Lineage Information into Statistical Tests for Detecting Selection in Ig Sequences , 2014, The Journal of Immunology.

[15]  Steven H. Kleinstein,et al.  Models of Somatic Hypermutation Targeting and Substitution Based on Synonymous Mutations from High-Throughput Immunoglobulin Sequencing Data , 2013, Front. Immunol..

[16]  P. Hodgkin,et al.  The clonal selection theory: 50 years since the revolution , 2007, Nature Immunology.

[17]  A. Bergman,et al.  V-region mutation in vitro, in vivo, and in silico reveal the importance of the enzymatic properties of AID and the sequence environment , 2009, Proceedings of the National Academy of Sciences.

[18]  L. Staudt,et al.  Generation of antibody diversity in the immune response of BALB/c mice to influenza virus hemagglutinin. , 1984, Proceedings of the National Academy of Sciences of the United States of America.

[19]  C. Kocks,et al.  Stepwise intraclonal maturation of antibody affinity through somatic hypermutation. , 1988, Proceedings of the National Academy of Sciences of the United States of America.

[20]  I. Maclennan,et al.  Mechanism of antigen-driven selection in germinal centres , 1989, Nature.

[21]  S. Sinha,et al.  Problems in using statistical analysis of replacement and silent mutations in antibody genes for determining antigen‐driven affinity selection , 2005, Immunology.

[22]  Steven H. Kleinstein,et al.  B cells populating the multiple sclerosis brain mature in the draining cervical lymph nodes , 2014, Science Translational Medicine.

[23]  C. Berek,et al.  Maturation of the immune response in germinal centers , 1991, Cell.

[24]  M. Kimura,et al.  On the probability of fixation of mutant genes in a population. , 1962, Genetics.

[25]  A Tramontano,et al.  Conformations of the third hypervariable region in the VH domain of immunoglobulins. , 1998, Journal of molecular biology.

[26]  Y. Louzoun,et al.  Rep‐Seq: uncovering the immunological repertoire through next‐generation sequencing , 2012, Immunology.

[27]  V. Giudicelli,et al.  IMGT unique numbering for immunoglobulin and T cell receptor variable domains and Ig superfamily V-like domains. , 2003, Developmental and comparative immunology.

[28]  Yoram Louzoun,et al.  Multi Step Selection in Ig H Chains is Initially Focused on CDR3 and Then on Other CDR Regions , 2013, Front. Immunol..

[29]  Uri Hershberg,et al.  Improved methods for detecting selection by mutation analysis of Ig V region sequences. , 2008, International immunology.

[30]  David A. Hafler,et al.  pRESTO: a toolkit for processing high-throughput sequencing raw reads of lymphocyte receptor repertoires , 2014, Bioinform..