Viral Capsid, Antibody, and Receptor Interactions: Experimental Analysis of the Antibody Escape Evolution of Canine Parvovirus.

Canine parvovirus (CPV) is a small nonenveloped single-stranded DNA virus that causes serious diseases in dogs worldwide. The original strain of the virus (CPV-2) emerged in dogs during the late 1970s due to a host range switch of a virus similar to the feline panleukopenia virus that infected another host. The virus that emerged in dogs had altered capsid receptor and antibody binding sites, with some changes affecting both functions. Further receptor and antibody binding changes arose when the virus became better adapted to dogs or to other hosts. Here, we used in vitro selection and deep sequencing to reveal how two antibodies with known interactions select for escape mutations in CPV. The antibodies bound two distinct epitopes, and one largely overlapped the host receptor binding site. We also generated mutated antibody variants with altered binding structures. Viruses were passaged with wild-type (WT) or mutated antibodies, and their genomes were deep sequenced during the selective process. A small number of mutations were detected only within the capsid protein gene during the first few passages of selection, and most sites remained polymorphic or were slow to go to fixation. Mutations arose both within and outside the antibody binding footprints on the capsids, and all avoided the transferrin receptor type 1 binding footprint. Many selected mutations matched those that have arisen in the natural evolution of the virus. The patterns observed reveal the mechanisms by which these variants have been selected in nature and provide a better understanding of the interactions between antibody and receptor selections. IMPORTANCE Antibodies protect animals against infection by many different viruses and other pathogens, and we are gaining new information about the epitopes that induce antibody responses against viruses and the structures of the bound antibodies. However, less is known about the processes of antibody selection and antigenic escape and the constraints that apply in this system. Here, we used an in vitro model system and deep genome sequencing to reveal the mutations that arose in the virus genome during selection by each of two monoclonal antibodies or their mutated variants. High-resolution structures of each of the Fab:capsid complexes revealed their binding interactions. The wild-type antibodies or their mutated variants allowed us to examine how changes in antibody structure influence the mutational selection patterns seen in the virus. The results shed light on the processes of antibody binding, neutralization escape, and receptor binding, and they likely have parallels for many other viruses.

Computational Analysis of Mutations in the Receptor-Binding Domain of SARS-CoV-2 Spike and Their Effects on Antibody Binding.

Currently, SARS-CoV-2 causing coronavirus disease 2019 (COVID-19) is responsible for one of the most deleterious pandemics of our time. The interaction between the ACE2 receptors at the surface of human cells and the viral Spike (S) protein triggers the infection, making the receptor-binding domain (RBD) of the SARS-CoV-2 S-protein a focal target for the neutralizing antibodies (Abs). Despite the recent progress in the development and deployment of vaccines, the emergence of novel variants of SARS-CoV-2 insensitive to Abs produced in response to the vaccine administration and/or monoclonal ones represent a potential danger. Here, we analyzed the diversity of neutralizing Ab epitopes and assessed the possible effects of single and multiple mutations in the RBD of SARS-CoV-2 S-protein on its binding affinity to various antibodies and the human ACE2 receptor using bioinformatics approaches. The RBD-Ab complexes with experimentally resolved structures were grouped into four clusters with distinct features at sequence and structure level. The performed computational analysis indicates that while single amino acid replacements in RBD may only cause partial impairment of the Abs binding, moreover, limited to specific epitopes, the variants of SARS-CoV-2 with multiple mutations, including some which were already detected in the population, may potentially result in a much broader antigenic escape. Further analysis of the existing RBD variants pointed to the trade-off between ACE2 binding and antigenic escape as a key limiting factor for the emergence of novel SAR-CoV-2 strains, as the naturally occurring mutations in RBD tend to reduce its binding affinity to Abs but not to ACE2. The results provide guidelines for further experimental studies aiming to identify high-risk RBD mutations that allow for an antigenic escape.

Antibody VH domain sequence analysis by a bioinformatics approach based on electronic amino acid properties may help to predict paratop location.

Gene as the basic functional unit of DNA encodes information about the product such as protein. The majority of proteins realize function through protein-protein interactions involving short protein motifs. However, some proteins such as antibodies are established by the rearrangement of several (V-D-J) gene segments with the potential addition of nontemplated nucleotides that may change information encoded by the respective gene segment used. Antibody VH domain sequence analysis by ISM bioinformatics approach that is based on amino acids physicochemical features, enable to distinguish the contribution of the information encoded by VH gene or generated during VDJ gene recombination for antibody-antigen interaction. The data presented in this report revealed the significance of CDRH3 for the interaction of antibody specific for immunogenic molecules while CDRH3 contribution is minor for antibody interaction with nonimmunogenic molecules such as haptens and native mammalian dsDNA. Thus, paratopes might be located in the CDRH3 or VH regions.

Sequence analysis of the DA domain of glycoprotein E2 of pestiviruses isolated from beef cattle in Southern Brazil.

The envelope glycoprotein E2 of pestiviruses is a major target for neutralizing antibodies. In this study, we analyzed the E2 DA domain of 43 pestiviruses from Southern Brazil. The isolates were identified as Bovine viral diarrhea virus (BVDV) subtypes 1a and 1b or BVDV-2b. Compared to reference strains, the BVDV-1 and -2 isolates had four and two mutations in the DA domain, respectively. All BVDV-2 isolates had a deletion of residues 724 and 725. All mutated amino acids in the BVDV isolates had the same aa substitution, and all were in previously identified antibody binding sites. It is possible that an immunity-mediated selection is acting on the pestiviruses circulating in Southern Brazil.

DutaFabs are engineered therapeutic Fab fragments that can bind two targets simultaneously.

We report the development of a platform of dual targeting Fab (DutaFab) molecules, which comprise two spatially separated and independent binding sites within the human antibody CDR loops: the so-called H-side paratope encompassing HCDR1, HCDR3 and LCDR2, and the L-side paratope encompassing LCDR1, LCDR3 and HCDR2. Both paratopes can be independently selected and combined into the desired bispecific DutaFabs in a modular manner. X-ray crystal structures illustrate that DutaFabs are able to bind two target molecules simultaneously at the same Fv region comprising a VH-VL heterodimer. In the present study, this platform is applied to generate DutaFabs specific for VEGFA and PDGF-BB, which show high affinities, physico-chemical stability and solubility, as well as superior efficacy over anti-VEGF monotherapy in vivo. These molecules exemplify the usefulness of DutaFabs as a distinct class of antibody therapeutics, which is currently being evaluated in patients.

Grabbing the Bull by Both Horns: Bovine Ultralong CDR-H3 Paratopes Enable Engineering of 'Almost Natural' Common Light Chain Bispecific Antibodies Suitable For Effector Cell Redirection.

A subset of antibodies found in cattle comprises ultralong CDR-H3 regions of up to 70 amino acids. Interestingly, this type of immunoglobulin usually pairs with the single germline VL gene, V30 that is typically very conserved in sequence. In this work, we have engineered ultralong CDR-H3 common light chain bispecific antibodies targeting Epidermal Growth Factor Receptor (EGFR) on tumor cells as well as Natural Cytotoxicity Receptor NKp30 on Natural Killer (NK) cells. Antigen-specific common light chain antibodies were isolated by yeast surface display by means of pairing CDR-H3 diversities following immunization with a single V30 light chain. After selection, EGFR-targeting paratopes as well as NKp30-specific binders were combined into common light chain bispecific antibodies by exploiting the strand-exchange engineered domain (SEED) technology for heavy chain heterodimerization. Biochemical characterization of resulting bispecifics revealed highly specific binding to the respective antigens as well as simultaneous binding to both targets. Most importantly, engineered cattle-derived bispecific common light chain molecules elicited potent NK cell redirection and consequently tumor cell lysis of EGFR-overexpressing cells as well as robust release of proinflammatory cytokine interferon-γ. Taken together, this data is giving clear evidence that bovine bispecific ultralong CDR-H3 common light chain antibodies are versatile for biotechnological applications.

Bi-paratopic and multivalent VH domains block ACE2 binding and neutralize SARS-CoV-2.

Neutralizing agents against SARS-CoV-2 are urgently needed for the treatment and prophylaxis of COVID-19. Here, we present a strategy to rapidly identify and assemble synthetic human variable heavy (VH) domains toward neutralizing epitopes. We constructed a VH-phage library and targeted the angiotensin-converting enzyme 2 (ACE2) binding interface of the SARS-CoV-2 Spike receptor-binding domain (Spike-RBD). Using a masked selection approach, we identified VH binders to two non-overlapping epitopes and further assembled these into multivalent and bi-paratopic formats. These VH constructs showed increased affinity to Spike (up to 600-fold) and neutralization potency (up to 1,400-fold) on pseudotyped SARS-CoV-2 virus when compared to standalone VH domains. The most potent binder, a trivalent VH, neutralized authentic SARS-CoV-2 with a half-maximal inhibitory concentration (IC50) of 4.0 nM (180 ng ml-1). A cryo-EM structure of the trivalent VH bound to Spike shows each VH domain engaging an RBD at the ACE2 binding site, confirming our original design strategy.

Conformational changes in a Vernier zone region: Implications for antibody dual specificity.

Understanding the determinants of antibody specificity is one of the challenging tasks in antibody development. Monospecific antibodies are still dominant in approved antibody therapeutics but there is a significant body of work to show that multispecific antibodies can increase the overall therapeutic effect. Dual-specific or "Two-in-One" antibodies can bind to two different antigens separately with the same antigen-binding site as opposed to bispecifics, which simultaneously bind to two different antigens through separate antigen-binding units. These nonstandard dual-specific antibodies were recently shown to be promising for new antibody-based therapeutics. Here, we physicochemically and structurally analyzed six different antibodies of which two are monospecific and four are dual-specific antibodies derived from monospecific templates to gain insight about dual-specificity determinants. These dual-specific antibodies can target both human epidermal growth factor receptor 2 and vascular endothelial growth factor at different binding affinities. We showed that a particular region of clustered Vernier zone residues might play key roles in gaining dual specificity. While there are minimal intramolecular interactions between a certain Vernier zone region, namely LV4 and LCDR1 of monospecific template, there is a significant structural change and consequently close contact formation between LV4-LCDR1 loops of derived dual-specific antibodies. Although Vernier zone residues were previously shown to be important for humanization applications, they are mostly underestimated in the literature. Here, we also aim to resurrect Vernier zone residues for antibody engineering efforts.

HIV-1 gp120-CD4-Induced Antibody Complex Elicits CD4 Binding Site-Specific Antibody Response in Mice.

Elicitation of broadly neutralizing Ab (bNAb) responses toward the conserved HIV-1 envelope (Env) CD4 binding site (CD4bs) by vaccination is an important goal for vaccine development and yet to be achieved. The outcome of previous immunogenicity studies suggests that the limited accessibility of the CD4bs and the presence of predominant nonneutralizing determinants (nND) on Env may impede the elicitation of bNAbs and their precursors by vaccination. In this study, we designed a panel of novel immunogens that 1) preferentially expose the CD4bs by selective elimination of glycosylation sites flanking the CD4bs, and 2) minimize the nND immune response by engineering fusion proteins consisting of gp120 Core and one or two CD4-induced (CD4i) mAbs for masking nND epitopes, referred to as gp120-CD4i fusion proteins. As expected, the fusion proteins possess improved antigenicity with retained affinity for VRC01-class, CD4bs-directed bNAbs and dampened affinity for nonneutralizing Abs. We immunized C57BL/6 mice with these fusion proteins and found that overall the fusion proteins elicit more focused CD4bs Ab response than prototypical gp120 Core by serological analysis. Consistently, we found that mice immunized with selected gp120-CD4i fusion proteins have higher frequencies of germinal center-activated B cells and CD4bs-directed memory B cells than those inoculated with parental immunogens. We isolated three mAbs from mice immunized with selected gp120-CD4i fusion proteins and found that their footprints on Env are similar to VRC01-class bNAbs. Thus, using gp120-CD4i fusion proteins with selective glycan deletion as immunogens could focus Ab response toward CD4bs epitope.

Primary structure of human neutrophil antigens 1a and 1b.

BACKGROUND: Neutrophil specific Fcγ receptor IIIb (CD16b) is a low-affinity IgG receptor. Its polymorphic variants are associated with human neutrophil antigens (HNA). HNA-1a and HNA-1b differ in four amino acids. Immunization can lead to the production of alloantibodies. The exact contribution of four amino acid exchanges for the formation of HNA-1a, -1b epitopes is currently unknown. STUDY DESIGN AND METHODS: Permutation of each polymorphic amino acid from wild-type CD16b cDNA constructs was performed and expressed on HEK293 cells. All 16 receptor variants were produced and tested against 19 well-characterized HNA antisera in an antigen capture assay. RESULTS: Analyzing the reaction pattern revealed that anti-HNA-1a antibodies can bind whenever asparagine (N) is present in position 65, regardless of the three other positions (CD16b *N**). Anti-HNA-1b antibodies can bind when serine (S) is present in position 36 (CD16b S***), when N is present in position 82 (CD16b **N*), or both (CD16b S*N*). CD16b variants with N65 and S36 and/or N82 (such as CD16b SNN*) bind both, anti-HNA-1a and anti-HNA-1b alloantibodies. If these specific amino acids are missing (as in CD16b RSD*), no antibodies will bind. CONCLUSION: Whereas the primary structure of HNA-1a and HNA-1b usually differs in four amino acids, epitope composition is not "antithetical". N65 alone determines the presence of HNA-1a, and S36 and/or N82 determine the presence of HNA-1b. Amino acid 106 does not participate in epitope formation. Our findings are of specific relevance when a HNA-1 phenotype is predicted from a genotype.

Kappa-on-Heavy (KoH) bodies are a distinct class of fully-human antibody-like therapeutic agents with antigen-binding properties.

We describe a Kappa-on-Heavy (KoH) mouse that produces a class of highly diverse, fully human, antibody-like agents. This mouse was made by replacing the germline variable sequences of both the Ig heavy-chain (IgH) and Ig kappa (IgK) loci with the human IgK germline variable sequences, producing antibody-like molecules with an antigen binding site made up of 2 kappa variable domains. These molecules, named KoH bodies, structurally mimic naturally existing Bence-Jones light-chain dimers in their variable domains and remain wild-type in their antibody constant domains. Unlike artificially diversified, nonimmunoglobulin alternative scaffolds (e.g., DARPins), KoH bodies consist of a configuration of normal Ig scaffolds that undergo natural diversification in B cells. Monoclonal KoH bodies have properties similar to those of conventional antibodies but exhibit an enhanced ability to bind small molecules such as the endogenous cardiotonic steroid marinobufagenin (MBG) and nicotine. A comparison of crystal structures of MBG bound to a KoH Fab versus a conventional Fab showed that the KoH body has a much deeper binding pocket, allowing MBG to be held 4 Å further down into the combining site between the 2 variable domains.

Effects of a remote mutation from the contact paratope on the structure of CDR-H3 in the anti-HIV neutralizing antibody PG16.

PG16 is a broadly neutralizing antibody to the human immunodeficiency virus (HIV). A crystal structure of PG16 revealed that the unusually long 28-residue complementarity determining region (CDR) H3 forms a unique subdomain, referred to as a "hammerhead", that directly contacts the antigen. The hammerhead apparently governs the function of PG16 while a previous experimental assay showed that the mutation of TyrH100Q to Ala, which does not directly contact the antigen, decreased the neutralization ability of PG16. However, the molecular mechanism by which a remote mutation from the hammerhead or contact paratope affects the neutralization potency has remained unclear. Here, we performed molecular dynamics simulations of the wild-type and variants (TyrH100Q to Ala, and TyrH100Q to Phe) of PG16, to clarify the effects of these mutations on the dynamics of CDR-H3. Our simulations revealed that the structural rigidity of the CDR-H3 in PG16 is attributable to the hydrogen bond interaction between TyrH100Q and ProH99, as well as the steric support by TyrH100Q. The loss of both interactions increases the intrinsic fluctuations of the CDR-H3 in PG16, leading to a conformational transition of CDR-H3 toward an inactive state.

A comprehensive search of functional sequence space using large mammalian display libraries created by gene editing.

The construction of large libraries in mammalian cells allows the direct screening of millions of molecular variants for binding properties in a cell type relevant for screening or production. We have created mammalian cell libraries of up to 10 million clones displaying a repertoire of IgG-formatted antibodies on the cell surface. TALE nucleases or CRISPR/Cas9 were used to direct the integration of the antibody genes into a single genomic locus, thereby rapidly achieving stable expression and transcriptional normalization. The utility of the system is illustrated by the affinity maturation of a PD-1-blocking antibody through the systematic mutation and functional survey of 4-mer variants within a 16 amino acid paratope region. Mutating VH CDR3 only, we identified a dominant "solution" involving substitution of a central tyrosine to histidine. This appears to be a local affinity maximum, and this variant was surpassed by a lysine substitution when light chain variants were introduced. We achieve this comprehensive and quantitative interrogation of sequence space by combining high-throughput oligonucleotide synthesis with mammalian display and flow cytometry operating at the multi-million scale.

CD20-Mimotope Peptides: A Model to Define the Molecular Basis of Epitope Spreading.

Antigen-mimicking peptide (mimotope)-based vaccines are one of the most promising forms of active-immunotherapy. The main drawback of this approach is that it induces antibodies that react poorly with the nominal antigen. The aim of this study was to investigate the molecular basis underlying the weak antibody response induced against the naïve protein after peptide vaccination. For this purpose, we analyzed the fine specificity of monoclonal antibodies (mAb) elicited with a 13-mer linear peptide, complementary to theantigen-combining site of the anti-CD20 mAb, Rituximab, in BALB/c mice. Anti-peptide mAb competed with Rituximab for peptide binding. Even so, they recognized a different antigenic motif from the one recognized by Rituximab. This explains their lack of reactivity with membrane (naïve) CD20. These data indicate that even on a short peptide the immunogenic and antigenic motifs may be different. These findings highlight an additional mechanism for epitope spreading and should be taken into account when designing peptides for vaccine purposes.

Defining the structural basis for human alloantibody binding to human leukocyte antigen allele HLA-A*11:01.

Our understanding of the conformational and electrostatic determinants that underlie targeting of human leukocyte antigens (HLA) by anti-HLA alloantibodies is principally based upon in silico modelling. Here we provide a biochemical/biophysical and functional characterization of a human monoclonal alloantibody specific for a common HLA type, HLA-A*11:01. We present a 2.4 Å resolution map of the binding interface of this antibody on HLA-A*11:01 and compare the structural determinants with those utilized by T-cell receptor (TCR), killer-cell immunoglobulin-like receptor (KIR) and CD8 on the same molecule. These data provide a mechanistic insight into the paratope-epitope relationship between an alloantibody and its target HLA molecule in a biological context where other immune receptors are concomitantly engaged. This has important implications for our interpretation of serologic binding patterns of anti-HLA antibodies in sensitized individuals and thus, for the biology of human alloresponses.

SAMURAI (Solid-phase Assisted Mutagenesis by Uracil Restriction for Accurate Integration) for antibody affinity maturation and paratope mapping.

Mutagenesis libraries are essential for combinatorial protein engineering. Despite improvements in gene synthesis and directed mutagenesis, current methodologies still have limitations regarding the synthesis of complete antibody single-chain variable fragment (scFv) genes and simultaneous diversification of all six CDRs. Here, we describe the generation of mutagenesis libraries for antibody affinity maturation using a cell-free solid-phase technique for annealing of single-strand mutagenic oligonucleotides. The procedure consists of PCR-based incorporation of uracil into a wild-type template, bead-based capture, elution of single-strand DNA, and in vitro uracil excision enzyme based degradation of the template DNA. Our approach enabled rapid (8 hours) mutagenesis and automated cloning of 50 position-specific alanine mutants for mapping of a scFv antibody paratope. We further exemplify our method by generating affinity maturation libraries with diversity introduced in critical, nonessential, or all CDR positions randomly. Assessment with Illumina deep sequencing showed less than 1% wild-type in two libraries and the ability to diversify all CDR positions simultaneously. Selections of the libraries with bacterial display and deep sequencing evaluation of the selection output showed that diversity introduced in non-essential positions allowed for a more effective enrichment of improved binders compared to the other two diversification strategies.

Improvement of antibody functionality by structure-guided paratope engraftment.

Broadly neutralizing antibodies (bNAbs) represent a promising alternative to antiretroviral drugs for HIV-1 prevention and treatment. Selected antibodies to the CD4-binding site bolster envelope trimer binding via quaternary contacts. Here, we rationally engraft a new paratope, i.e., the extended heavy-chain framework region 3 (FR3) loop of VRC03, which mediates quaternary interaction, onto several potent bNAbs, enabling them to reach an adjacent gp120 protomer. The interactive quaternary surface is delineated by solving the crystal structure of two FR3 loop-chimeric antibodies. Chimerization enhances the neutralizing activity of several potent bNAbs against a majority of global HIV-1 strains. Compared to unmodified antibodies, chimeric antibodies display lower autoreactivity and prolonged in vivo half-life in huFcRn mice and rhesus macaques. Thus, paratope engraftment may be used to expand the epitope repertory of natural antibodies, improving their functionality for disease prevention and treatment.

Designing antibody against highly conserved region of dengue envelope protein by in silico screening of scFv mutant library.

Dengue being one of the deadliest diseases of tropical regions, enforces to put continuous efforts for the development of vaccine and effective therapeutics. Most of the antibodies generated during dengue infection are non-neutralizing and cause antibody dependent enhancement. Hence, making a potent neutralizing antibody against all four dengue serotypes could be very effective for the treatment. However, designing a single antibody for all serotypes is difficult due to variation in protein sequences. Therefore, the objective is to identify conserved region of dengue envelope protein and then develop an antibody against that conserved region. Before advancing to the development of such an antibody, it is desirable to validate the interactions between antibody and dengue envelope protein. In silico analysis of such interactions provides a good platform to find out a suitable region to design and construct an antibody against it by analyzing antigen-antibody interaction before synthesizing the antibody. In this study, two highly conserved regions of dengue envelope protein were identified and an scFv was constructed against it. Both scFv and FuBc proteins were expressed in bacterial expression system and binding efficiency was analyzed by SPR analysis with KD value 2.3 µM. In order to improve binding efficiency, an in silico scFv mutant library was created which was virtually screened for higher binding efficiency. Six mutants with high binding efficiency were selected for further analysis. The binding ability of these mutants were predicted using simulation analysis which shows these mutations were stabilizing scFv-FuBc complex.

Anti-PD1 'SHR-1210' aberrantly targets pro-angiogenic receptors and this polyspecificity can be ablated by paratope refinement.

Monoclonal anti-programmed cell death 1 (PD1) antibodies are successful cancer therapeutics, but it is not well understood why individual antibodies should have idiosyncratic side-effects. As the humanized antibody SHR-1210 causes capillary hemangioma in patients, a unique toxicity amongst anti-PD1 antibodies, we performed human receptor proteome screening to identify nonspecific interactions that might drive angiogenesis. This screen identified that SHR-1210 mediated aberrant, but highly selective, low affinity binding to human receptors such as vascular endothelial growth factor receptor 2 (VEGFR2), frizzled class receptor 5 and UL16 binding protein 2 (ULBP2). SHR-1210 was found to be a potent agonist of human VEGFR2, which may thereby drive hemangioma development via vascular endothelial cell activation. The v-domains of SHR-1210's progenitor murine monoclonal antibody 'Mab005' also exhibited off-target binding and agonism of VEGFR2, proving that the polyspecificity was mediated by the original mouse complementarity-determining regions (CDRs), and had survived the humanization process. Molecular remodelling of SHR-1210 by combinatorial CDR mutagenesis led to deimmunization, normalization of binding affinity to human and cynomolgus PD1, and increased potency in PD1/PD-L1 blockade. Importantly, CDR optimization also ablated all off-target binding, rendering the resulting antibodies fully PD1-specific. As the majority of changes to the paratope were found in the light chain CDRs, the germlining of this domain drove the ablation of off-target binding. The combination of receptor proteome screening and optimization of the antibody binding interface therefore succeeded in generating novel, higher-potency, specificity-enhanced therapeutic IgGs from a single, clinically sub-optimal progenitor. This study showed that highly-specific off-target binding events might be an under-appreciated phenomenon in therapeutic antibody development, but that these unwanted properties can be fully ameliorated by paratope refinement.

Tracing Antibody Repertoire Evolution by Systems Phylogeny.

Antibody evolution studies have been traditionally limited to either tracing a single clonal lineage (B cells derived from a single V-(D)-J recombination) over time or examining bulk functionality changes (e.g., tracing serum polyclonal antibody proteins). Studying a single B cell disregards the majority of the humoral immune response, whereas bulk functional studies lack the necessary resolution to analyze the co-existing clonal diversity. Recent advances in high-throughput sequencing (HTS) technologies and bioinformatics have made it possible to examine multiple co-evolving antibody monoclonal lineages within the context of a single repertoire. A plethora of accompanying methods and tools have been introduced in hopes of better understanding how pathogen presence dictates the global evolution of the antibody repertoire. Here, we provide a comprehensive summary of the tremendous progress of this newly emerging field of systems phylogeny of antibody responses. We present an overview encompassing the historical developments of repertoire phylogenetics, state-of-the-art tools, and an outlook on the future directions of this fast-advancing and promising field.