Jason J. Lavinder

High-throughput sequencing of the paired human immunoglobulin heavy and light chain repertoire

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 (VH and VL) 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 × 104 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 VH:VL linkage PCR. The linked transcripts are analyzed by Illumina high-throughput sequencing. We validated the fidelity of VH:VL 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.

High-Throughput Thermal Scanning: A General, Rapid Dye-Binding Thermal Shift Screen for Protein Engineering

The low stability of natural proteins often limits their use in therapeutic, industrial, and research applications. The scale and throughput of methods such as circular dichroism, fluorescence spectroscopy, and calorimetry severely limit the number of variants that can be examined. Here we demonstrate a high-throughput thermal scanning (HTTS) method for determining the approximate stabilities of protein variants at high throughput and low cost. The method is based on binding to a hydrophobic dye akin to ANS, which fluoresces upon binding to molten globules and thermal denaturation intermediates. No inherent properties of the protein, such as enzymatic activity or presence of an intrinsic fluorophore, are required. Very small sample sizes are analyzed using a real-time PCR machine, enabling the use of high-throughput purification. We show that the apparent T(M) values obtained from HTTS are approximately linearly related to those from CD thermal denaturation for a series of four-helix bundle hydrophobic core variants. We demonstrate similar results for a small set of TIM barrel variants. This inexpensive, general, and scaleable approach enables the search for conservative, stable mutants of biotechnologically important proteins and provides a method for statistical correlation of sequence-stability relationships.

Identification and characterization of the constituent human serum antibodies elicited by vaccination.

Significance Most vaccines confer immunity by eliciting long-term production of antibodies that bind to and neutralize the vaccine antigen. Remarkably, very little is known regarding the identities, sequence diversity, relative concentrations, or binding functionalities of the mAbs that comprise the serum repertoire elicited by vaccination. Here, we have delineated the constituent antibodies of the human serum IgG repertoire after vaccination and examined their relationship to the antibody V gene repertoire encoded by circulating B cells. The results detail the molecular composition and characteristics of the vaccine-specific serum antibody repertoire and demonstrate differences between the end-point response (the serum antibodies) and the peripheral B cells responding to the vaccine. Most vaccines confer protection via the elicitation of serum antibodies, yet more than 100 y after the discovery of antibodies, the molecular composition of the human serum antibody repertoire to an antigen remains unknown. Using high-resolution liquid chromatography tandem MS proteomic analyses of serum antibodies coupled with next-generation sequencing of the V gene repertoire in peripheral B cells, we have delineated the human serum IgG and B-cell receptor repertoires following tetanus toxoid (TT) booster vaccination. We show that the TT+ serum IgG repertoire comprises ∼100 antibody clonotypes, with three clonotypes accounting for >40% of the response. All 13 recombinant IgGs examined bound to vaccine antigen with Kd ∼ 10−8–10−10 M. Five of 13 IgGs recognized the same linear epitope on TT, occluding the binding site used by the toxin for cell entry, suggesting a possible explanation for the mechanism of protection conferred by the vaccine. Importantly, only a small fraction (<5%) of peripheral blood plasmablast clonotypes (CD3−CD14−CD19+CD27++CD38++CD20−TT+) at the peak of the response (day 7), and an even smaller fraction of memory B cells, were found to encode antibodies that could be detected in the serological memory response 9 mo postvaccination. This suggests that only a small fraction of responding peripheral B cells give rise to the bone marrow long-lived plasma cells responsible for the production of biologically relevant amounts of vaccine-specific antibodies (near or above the Kd). Collectively, our results reveal the nature and dynamics of the serological response to vaccination with direct implications for vaccine design and evaluation.

Molecular deconvolution of the monoclonal antibodies that comprise the polyclonal serum response.

We have developed and validated a methodology for determining the antibody composition of the polyclonal serum response after immunization. Pepsin-digested serum IgGs were subjected to standard antigen-affinity chromatography, and resulting elution, wash, and flow-through fractions were analyzed by bottom-up, liquid chromatography–high-resolution tandem mass spectrometry. Identification of individual monoclonal antibodies required the generation of a database of IgG variable gene (V-gene) sequences constructed by NextGen sequencing of mature B cells. Antibody V-gene sequences are characterized by short complementarity determining regions (CDRs) of high diversity adjacent to framework regions shared across thousands of IgGs, greatly complicating the identification of antigen-specific IgGs from proteomically observed peptides. By mapping peptides marking unique VH CDRH3 sequences, we identified a set of V-genes heavily enriched in the affinity chromatography elution, constituting the serum polyclonal response. After booster immunization in a rabbit, we find that the antigen-specific serum immune response is oligoclonal, comprising antibodies encoding 34 different CDRH3s that group into 30 distinct antibody VH clonotypes. Of these 34 CDRH3s, 12 account for ∼60% of the antigen-specific CDRH3 peptide mass spectral counts. For comparison, antibodies with 18 different CDRH3s (12 clonotypes) were represented in the antigen-specific IgG fraction from an unimmunized rabbit that fortuitously displayed a moderate titer for BSA. Proteomically identified antibodies were synthesized and shown to display subnanomolar affinities. The ability to deconvolute the polyclonal serum response is likely to be of key importance for analyzing antibody responses after vaccination and for more completely understanding adaptive immune responses in health and disease.

Molecular-level analysis of the serum antibody repertoire in young adults before and after seasonal influenza vaccination

Molecular understanding of serological immunity to influenza has been confounded by the complexity of the polyclonal antibody response in humans. Here we used high-resolution proteomics analysis of immunoglobulin (referred to as Ig-seq) coupled with high-throughput sequencing of transcripts encoding B cell receptors (BCR-seq) to quantitatively determine the antibody repertoire at the individual clonotype level in the sera of young adults before and after vaccination with trivalent seasonal influenza vaccine. The serum repertoire comprised between 40 and 147 clonotypes that were specific to each of the three monovalent components of the trivalent influenza vaccine, with boosted pre-existing clonotypes accounting for ∼60% of the response. An unexpectedly high fraction of serum antibodies recognized both the H1 and H3 monovalent vaccines. Recombinant versions of these H1 + H3 cross-reactive antibodies showed broad binding to hemagglutinins (HAs) from previously circulating virus strains; several of these antibodies, which were prevalent in the serum of multiple donors, recognized the same conserved epitope in the HA head domain. Although the HA-head-specific H1 + H3 antibodies did not show neutralization activity in vitro, they protected mice against infection with the H1N1 and H3N2 virus strains when administered before or after challenge. Collectively, our data reveal unanticipated insights regarding the serological response to influenza vaccination and raise questions about the added benefits of using a quadrivalent vaccine instead of a trivalent vaccine.

Direct single-molecule observation of a protein living in two opposed native structures

Biological activity in proteins requires them to share the energy landscape for folding and global conformational motions, 2 key determinants of function. Although most structural studies to date have focused on fluctuations around a single structural basin, we directly observe the coexistence of 2 symmetrically opposed conformations for a mutant of the Rop-homodimer (Repressor of Primer) in single-molecule fluorescence resonance energy transfer (smFRET) measurements. We find that mild denaturing conditions can affect the sensitive balance between the conformations, generating an equilibrium ensemble consisting of 2 equally occupied structural basins. Despite the need for large-scale conformational rearrangement, both native structures are dynamically and reversibly adopted for the same paired molecules without separation of the constituent monomers. Such an ability of some proteins or protein complexes to switch between conformations by thermal fluctuations and/or minor environmental changes could be central to their ability to control biological function.

Protein stability by number: high-throughput and statistical approaches to one of protein science's most difficult problems.

Most proteins are only barely stable, which impedes research, complicates therapeutic applications, and makes proteins susceptible to pathologically destabilizing mutations. Our ability to predict the thermodynamic consequences of even single point mutations is still surprisingly limited, and established methods of measuring stability are slow. Recent advances are bringing protein stability studies into the high-throughput realm. Some methods are based on inferential read-outs such as activity, proteolytic resistance or split-protein fragment reassembly. Other methods use miniaturization of direct measurements, such as intrinsic fluorescence, H/D exchange, cysteine reactivity, aggregation and hydrophobic dye binding (DSF). Protein engineering based on statistical analysis (consensus and correlated occurrences of amino acids) is promising, but much work remains to understand and implement these methods.

Next-generation sequencing and protein mass spectrometry for the comprehensive analysis of human cellular and serum antibody repertoires

Recent developments of high-throughput technologies are enabling the molecular-level analysis and bioinformatic mining of antibody-mediated (humoral) immunity in humans at an unprecedented level. These approaches explore either the sequence space of B-cell receptor repertoires using next-generation deep sequencing (BCR-seq), or the amino acid identities of antibody in blood using protein mass spectrometry (Ig-seq), or both. Generalizable principles about the molecular composition of the protective humoral immune response are being defined, and as such, the field could supersede traditional methods for the development of diagnostics, vaccines, and antibody therapeutics. Three key challenges remain and have driven recent advances: (1) incorporation of innovative techniques for paired BCR-seq to ascertain the complete antibody variable-domain VH:VL clonotype, (2) integration of proteomic Ig-seq with BCR-seq to reveal how the serum antibody repertoire compares with the antibody repertoire encoded by circulating B cells, and (3) a demand to link antibody sequence data to functional meaning (binding and protection).

Systematic Characterization and Comparative Analysis of the Rabbit Immunoglobulin Repertoire

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.

Proteomic identification of monoclonal antibodies from serum.

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.