Joshua A. Weinstein

High-Throughput Sequencing of the Zebrafish Antibody Repertoire

Antibody Repertoire Revealed Antibodies produced by B cells protect us against infection by a wide array of pathogens. Such wide-ranging responses are possible because the specific domain of the antibody that binds to the invader is highly variable owing to the somatic recombination of an inherited set of gene segments that encode the antibody molecule. Although antibody diversity within an individual organism is well-established, the specific antibody repertoire of an individual organism has not been characterized. Using high-throughput sequencing technology, Weinstein et al. (p. 807) characterized the diversity of the antigen-binding domain of the antibody heavy chain in zebrafish. The antibody repertoire of individual fish covered at least 50% of the possible gene combinations. Although the specific gene combinations differed between fish, a similar frequency distribution of the repertoire was observed. Unexpected examples of evolutionary conversion were also seen, with the same antibody observed in different animals. Sequencing of immunoglobulin messenger RNA characterizes the diversity of the antibody repertoire in individual zebrafish. Despite tremendous progress in understanding the nature of the immune system, the full diversity of an organism’s antibody repertoire is unknown. We used high-throughput sequencing of the variable domain of the antibody heavy chain from 14 zebrafish to analyze VDJ usage and antibody sequence. Zebrafish were found to use between 50 and 86% of all possible VDJ combinations and shared a similar frequency distribution, with some correlation of VDJ patterns between individuals. Zebrafish antibodies retained a few thousand unique heavy chains that also exhibited a shared frequency distribution. We found evidence of convergence, in which different individuals made the same antibody. This approach provides insight into the breadth of the expressed antibody repertoire and immunological diversity at the level of an individual organism.

Lineage Structure of the Human Antibody Repertoire in Response to Influenza Vaccination

High-throughput sequencing suggests age-related differences in the changes in the human antibody repertoire after influenza vaccination. Antibodies Act Their Age One of the main advantages of the immune system in fighting infection is its ability to diversify—antibody and T cell receptor genes physically rearrange, creating a repertoire of potential responses that can be called upon and expanded if needed. However, the very diversity of this repertoire is what makes immune responses hard to study. Although we can know how any individual B cell or antibody responds to stimulation, getting the big picture is much more difficult. Adding in another variable, such as time, further complicates things. Now, Jiang et al. use high-throughput long read sequencing to characterize the human antibody response after influenza vaccination. People of age are thought to have altered immune systems compared to younger individuals. However, why exactly how the antibody repertoire changes with age remains unclear. By analyzing more than 5 million antibody heavy chain sequences, the authors were able to compare isotype diversity, lineage structure, and mutational activity in differently aged populations. They found that the elderly have fewer lineages, with reduced diversity, compared with younger subjects; however, the antibodies present before vaccination had higher levels of mutation. These same techniques can be used to study individual-specific immune responses and may aid in optimizing vaccination in the future. 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.

Genetic measurement of memory B-cell recall using antibody repertoire sequencing

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.

Determinism and stochasticity during maturation of the zebrafish antibody repertoire

It is thought that the adaptive immune system of immature organisms follows a more deterministic program of antibody creation than is found in adults. We used high-throughput sequencing to characterize the diversifying antibody repertoire in zebrafish over five developmental time points. We found that the immune system begins in a highly stereotyped state with preferential use of a small number of V (variable) D (diverse) J (joining) gene segment combinations, but that this stereotypy decreases dramatically as the zebrafish mature, with many of the top VDJ combinations observed in 2-wk-old zebrafish virtually disappearing by 1 mo. However, we discovered that, in the primary repertoire, there are strong correlations in VDJ use that increase with zebrafish maturity, suggesting that VDJ recombination involves a level of deterministic programming that is unexpected. This stereotypy is masked by the complex diversification processes of antibody maturation; the variation and lack of correlation in full repertoires between individuals appears to be derived from randomness in clonal expansion during the affinity maturation process. These data provide a window into the mechanisms of VDJ recombination and diversity creation and allow us to better understand how the adaptive immune system achieves diversity.

A Multiscale Computational Approach to Dissect Early Events in the Erb Family Receptor Mediated Activation, Differential Signaling, and Relevance to Oncogenic Transformations

We describe a hierarchical multiscale computational approach based on molecular dynamics simulations, free energy-based molecular docking simulations, deterministic network-based kinetic modeling, and hybrid discrete/continuum stochastic dynamics protocols to study the dimer-mediated receptor activation characteristics of the Erb family receptors, specifically the epidermal growth factor receptor (EGFR). Through these modeling approaches, we are able to extend the prior modeling of EGF-mediated signal transduction by considering specific EGFR tyrosine kinase (EGFRTK) docking interactions mediated by differential binding and phosphorylation of different C-terminal peptide tyrosines on the RTK tail. By modeling signal flows through branching pathways of the EGFRTK resolved on a molecular basis, we are able to transcribe the effects of molecular alterations in the receptor (e.g., mutant forms of the receptor) to differing kinetic behavior and downstream signaling response. Our molecular dynamics simulations show that the drug sensitizing mutation (L834R) of EGFR stabilizes the active conformation to make the system constitutively active. Docking simulations show preferential characteristics (for wildtype vs. mutant receptors) in inhibitor binding as well as preferential enhancement of phosphorylation of particular substrate tyrosines over others. We find that in comparison to the wildtype system, the L834R mutant RTK preferentially binds the inhibitor erlotinib, as well as preferentially phosphorylates the substrate tyrosine Y1068 but not Y1173. We predict that these molecular level changes result in preferential activation of the Akt signaling pathway in comparison to the Erk signaling pathway for cells with normal EGFR expression. For cells with EGFR over expression, the mutant over activates both Erk and Akt pathways, in comparison to wildtype. These results are consistent with qualitative experimental measurements reported in the literature. We discuss these consequences in light of how the network topology and signaling characteristics of altered (mutant) cell lines are shaped differently in relationship to native cell lines.

Correlation of Gene Expression and Genome Mutation in Single B-Cells

High-throughput measurement of gene-expression and immune receptor repertoires have recently become powerful tools in the study of adaptive immune response. However, despite their now-widespread use, both tend to discard cell identity by treating cell populations in bulk, and therefore lose the correlation between genetic variability and gene-expression at the single cell level. In order to recover this information, we developed a method to simultaneously measure gene expression profiles and genome mutations in single cells. We applied this method by quantifying the relationships between gene expression and antibody mutation in ensembles of individual B-cells from immunized mice. The results reveal correlations reflecting the manner in which information propagates between a B-cell’s antigen receptors, its gene expression, and its mutagenic machinery, and demonstrate the power of this approach to illuminate both heterogeneity and physiology in cell populations.