Different adjuvanted pediatric HIV envelope vaccines induced distinct plasma antibody responses despite similar B cell receptor repertoires in infant rhesus macaques

Abstract

Different HIV vaccine regimens elicit distinct plasma antibody responses in both human and nonhuman primate models. Previous studies in human and non-human primate infants showed that adjuvants influenced the quality of plasma antibody responses induced by pediatric HIV envelope vaccine regimens. We recently reported that use of the 3M052-SE adjuvant and longer intervals between vaccinations are associated with higher magnitude of antibody responses in infant rhesus macaques. However, the impact of different adjuvants in HIV vaccine regimens on the developing infant B cell receptor (BCR) repertoire has not been studied. This study evaluated whether pediatric HIV envelope vaccine regimens with different adjuvants induced distinct antigen-specific memory B cell repertoires and whether specific immunoglobulin (Ig) immunogenetic characteristics are associated with higher magnitude of plasma antibody responses in vaccinated infant rhesus macaques. We utilized archived preclinical pediatric HIV vaccine studies PBMCs and tissue samples from 19 infant rhesus macaques immunized either with (i) HIV Env protein with a squalene adjuvant, (ii) MVA-HIV and Env protein co-administered using a 3-week interval, (iii) MVA-HIV prime/ protein boost with an extended 6-week interval between immunizations, or (iv) with HIV Env administered with 3M-052-SE adjuvant. Frequencies of vaccine-elicited HIV Env-specific memory B cells from PBMCs and tissues were similar across vaccination groups (frequency range of 0.06–1.72%). There was no association between vaccine-elicited antigen-specific memory B cell frequencies and plasma antibody titer or avidity. Moreover, the epitope specificity and Ig immunogenetic features of vaccine-elicited monoclonal antibodies did not differ between the different vaccine regimens. These data suggest that pediatric HIV envelope vaccine candidates with different adjuvants that previously induced higher magnitude and quality of plasma antibody responses in infant rhesus macaques were not driven by distinct antigen-specific memory BCR repertoires.

About the Authors:

Stella J. Berendam

Roles Data curation, Formal analysis, Visualization, Writing – original draft

¶‡ SJB and PKM shared first authors. MAM and GGF shared last authors.

Affiliation: Duke Human Vaccine Institute, Duke University Medical Center, Durham, North Carolina, United States of America

ORCID logo https://orcid.org/0000-0001-9622-2238

Papa K. Morgan-Asiedu

Roles Formal analysis, Investigation

¶‡ SJB and PKM shared first authors. MAM and GGF shared last authors.

Affiliation: Duke Human Vaccine Institute, Duke University Medical Center, Durham, North Carolina, United States of America

ORCID logo https://orcid.org/0000-0001-5622-524X

Riley J. Mangan

Roles Investigation

Affiliation: Duke Human Vaccine Institute, Duke University Medical Center, Durham, North Carolina, United States of America

Shuk Hang Li

Roles Data curation, Formal analysis

Affiliation: Duke Human Vaccine Institute, Duke University Medical Center, Durham, North Carolina, United States of America

Holly Heimsath

Roles Investigation

Affiliation: Duke Human Vaccine Institute, Duke University Medical Center, Durham, North Carolina, United States of America

Kan Luo

Roles Methodology

Affiliation: Duke Human Vaccine Institute, Duke University Medical Center, Durham, North Carolina, United States of America

Alan D. Curtis II

Roles Investigation, Methodology

Affiliation: Department of Microbiology and Immunology, Children’s Research Institute and Center for AIDS Research, School of Medicine, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America

Joshua A. Eudailey

Roles Investigation

Affiliations Duke Human Vaccine Institute, Duke University Medical Center, Durham, North Carolina, United States of America, Department of Pediatrics, Weill Cornell College of Medicine, New York City, New York, United States of America

Christopher B. Fox

Roles Resources

Affiliations Infectious Disease Research Institute (IDRI), Seattle, Washington State, United States of America, Department of Global Health, University of Washington, Seattle, Washington State, United States of America

Mark A. Tomai

Roles Resources

Affiliation: 3M Center, 3 M Drug Delivery Systems, St. Paul, Minnesota, United States of America

Bonnie Phillips

Roles Investigation

Affiliation: Department of Microbiology and Immunology, Children’s Research Institute and Center for AIDS Research, School of Medicine, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America

Hannah L. Itell

Roles Investigation

Affiliation: Duke Human Vaccine Institute, Duke University Medical Center, Durham, North Carolina, United States of America

Erika Kunz

Roles Investigation

Affiliation: Duke Human Vaccine Institute, Duke University Medical Center, Durham, North Carolina, United States of America

Michael Hudgens

Roles Formal analysis, Methodology

Affiliation: Department of Biostatistics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America

Kenneth Cronin

Roles Methodology

Affiliation: Duke Human Vaccine Institute, Duke University Medical Center, Durham, North Carolina, United States of America

Kevin Wiehe

Roles Formal analysis, Methodology, Writing – review & editing

Affiliation: Duke Human Vaccine Institute, Duke University Medical Center, Durham, North Carolina, United States of America

S. Munir Alam

Roles Methodology

Affiliation: Duke Human Vaccine Institute, Duke University Medical Center, Durham, North Carolina, United States of America

Koen K. A. Van Rompay

Roles Project administration, Supervision, Writing – review & editing

Affiliation: California National Primate Research Center, University of California at Davis, Davis, California, United States of America

ORCID logo https://orcid.org/0000-0002-7375-1337

Kristina De Paris

Roles Conceptualization, Funding acquisition, Supervision, Writing – review & editing

Affiliation: Department of Microbiology and Immunology, Children’s Research Institute and Center for AIDS Research, School of Medicine, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America

ORCID logo https://orcid.org/0000-0001-5199-3004

Sallie R. Permar

Roles Conceptualization, Funding acquisition, Supervision, Writing – review & editing

Affiliations Duke Human Vaccine Institute, Duke University Medical Center, Durham, North Carolina, United States of America, Department of Pediatrics, Weill Cornell College of Medicine, New York City, New York, United States of America

M. Anthony Moody

Roles Conceptualization, Supervision, Writing – review & editing

¶‡ SJB and PKM shared first authors. MAM and GGF shared last authors.

Affiliation: Duke Human Vaccine Institute, Duke University Medical Center, Durham, North Carolina, United States of America

Genevieve G. Fouda

Roles Conceptualization, Supervision, Writing – review & editing

* E-mail: [email protected]

¶‡ SJB and PKM shared first authors. MAM and GGF shared last authors.

Affiliation: Duke Human Vaccine Institute, Duke University Medical Center, Durham, North Carolina, United States of America

ORCID logo https://orcid.org/0000-0002-0960-1019

Introduction

In 2019, 85% of the estimated 1.3 million pregnant women living with HIV-1 globally received antiretroviral drugs to prevent transmission to their children [1]. While the implementation of antiretroviral prophylaxis has significantly decreased the global frequency of mother-to-child transmission (MTCT) of HIV-1, issues of maternal adherence to antiretroviral therapy (ART) [2, 3], development of ART-resistant viruses [4], and insufficient coverage of ART in some of the hardest-hit areas globally have limited the effectiveness of ART [5]. Furthermore, women with acute HIV-1 infection in late pregnancy or during the breastfeeding period are less likely to be diagnosed and receive treatment to prevent MTCT [6]. Thus, despite advancements in therapy, breast milk transmission still accounts for approximately 50% of pediatric HIV infections [7, 8]. Additional prevention strategies, such as a pediatric HIV-1 vaccine, are therefore critically needed to eradicate breast milk transmission of HIV-1.

Early efforts in HIV-1 vaccine development focused on the humoral arm of the immune system as other vaccines that successfully prevent viral diseases relied on antibodies for protection [9]. However, early-phase HIV-1 vaccine studies using recombinant HIV-1 envelope (Env) proteins showed no efficacy [10–12] with the exception of the RV144 trial that demonstrated moderate vaccine efficacy of 61% during the first year and an overall efficacy of 31% at 3.5 years after vaccination [13, 14]. Interestingly, non-neutralizing IgG that targeted the HIV-1 Env variable loops 1 and 2 (V1V2) were identified as correlates of protection from the RV144 study; meanwhile HIV-1 Env-specific IgA plasma antibodies were associated with lack of protection [15]. This finding reinvigorated the optimism that an effective HIV-1 vaccine is attainable. A phase 2b/3 study (HVTN 702 or Uhambo) in South Africa, which used a pox vector prime-protein boost vaccine regimen similar to RV144, albeit with distinct vaccine strains and a different adjuvant, did not reproduce the results from the RV144 study and showed no efficacy [16]. Nevertheless, these studies demonstrated that HIV-1 vaccine immunogens could induce robust levels of V1V2 IgG antibodies as well as polyfunctional CD4+ T cell responses as surrogate of possible protection [15, 17].

To date, only a few previous pediatric HIV-1 vaccine trials have been conducted and these trials demonstrated that immunization with recombinant subunit HIV gp120 vaccines (PACTG 230) or with canarypox vectors expressing HIV antigens (PACTG 326, HPTN 027) are safe and immunogenic [18–21]. Importantly, HIV-1 infant vaccination with MF-59 adjuvanted HIV gp120 was able to generate robust and durable Env-specific IgG responses including anti-V1V2 IgG responses with low levels of Env-specific IgA responses [22]. Similarly, neonatal rhesus macaques are capable of developing robust virus-specific humoral and cellular immune responses after SIV vaccination [23, 24]. These clinical and preclinical studies demonstrate the feasibility of initiating HIV-1 immunization within the first few days of life.

Our previous studies in non-human primate infants have indicated that several factors can modulate the quality of the vaccine-elicited antibody response in infants [25, 26]. Notably, we observed that extending the interval between immunizations and the use of toll-like receptor (TLR) agonist adjuvants can enhance the magnitude and breadth of the cellular and humoral responses [26, 27]. However, the influence of distinct HIV vaccine regimens and adjuvants on the developing infant B cell repertoire has not been studied and the relationship between the magnitude and quality of HIV-1 vaccine-elicited responses and immunogenetic characteristics of the infant B cell repertoire remains unclear. Taking advantage of archived plasma samples obtained from completed preclinical studies [25, 26], we assessed the antigen-specific B cell repertoire in HIV vaccinated infant rhesus macaques. Our results indicated that while the magnitude and quality of plasma antibody responses induced by pediatric HIV vaccine regimens with different adjuvants, the B cell repertoire profiles were not distinct between different vaccine regimens in infant rhesus macaques.

Methods

Animals

A total of 19 Simian Immunodeficiency Virus (SIV)-negative and type D retrovirus-negative newborn Indian-origin rhesus macaques (Macaca mulatta) were hand reared in the nursery of the California National Primate Research Center (CNPRC, Davis, CA) as previously described [25]. Animals were reared in accordance with the American Association for Accreditation of Laboratory Animal Care Standards, the guidelines of the Guide for the Care and Use of Laboratory Animals of the Institute for Laboratory Research, National Research Council, and the International Guiding Principles for Biomedical Research Involving Animals. All protocols were assessed and approved by the University of California at Davis Institutional Animal Care and Use Committee prior to beginning the study. Animals were randomly assigned to groups and anesthetized for vaccinations and sample collection as previously reported [25, 26].

Immunization regimens and study designs

Immunization regimens and archived samples used in this study for all 4 animal groups are summarized in Fig 1 and detailed study designs were previously published [25, 26]. Briefly, infants in group 1 (HIV Env+adjuvant, n = 5) were vaccinated at 0, 3 and 6 weeks of age intramuscularly (IM) with 5 x 108 infectious units [IU] of MVA/SIV gag/pol, with 15 μg of C.1086 gp120 administered IM in Span85-Tween 80-squalene (STS) adjuvant, and with 200 μg of C.1086 gp120 administered intranasally (IN) in Toll-like receptor 7 and 8 (TLR7/8) agonist, R848 adjuvant [25]. Infants in group 2 (Co-administration, n = 5) were vaccinated at 0, 3 and 6 weeks of age with similar regimen as group 1 with addition of the MVA-HIV Env (5 x 108 IU, IM). Infants from these two groups were followed for 19 and 15 weeks, respectively, after which they were euthanized to analyze vaccine-induced tissue responses. Infants in group 3 (Extended Interval, n = 4) received the same vaccine regimen as the co-administration group but were immunized at 0, 6 and 12 weeks and then boosted at 32 before euthanasia at week 35 [25]. Infants in group 4 (3M-052-SE, n = 4) were immunized at 0, 2 and 6 weeks of age with a combination of HIV.C.1086C gp120 and TV1 gp120 (15μg each, IM) in TLR7/8-based adjuvant 3M-052 formulated in stable emulsion (3M-052-SE) and was challenged with SHIV at week 15 [26].

[Figure omitted. See PDF.]

Fig 1. Animal study design and immunization schedule.

Infant rhesus macaques from four immunization schedules were included. HIV Env+adjuvant group was immunized with MVA-SIV gag/pol and HIV Envelope (Env) 1086c adjuvanted with 15 μg STS, intramuscular + 200 μg R848, intranasal. Co-administration group was immunized with MVA-SIV gag/pol, HIV Env 1086c adjuvanted with 15 μg STS, intramuscular + 200 μg R848, intranasal), and with MVA-HIV Env. Extended interval group was immunized with similar immunogens as Co-administration group with longer immunization intervals. 3M-052-SE group was immunized with HIV Env 1086c/TV1 bivalent adjuvanted with 15 μg/15 μg+3M-052-SE, intramuscular. Plus symbol (+) denoted necropsy, IM denoted intramuscular, IN denoted intranasal, and IU denoted infectious unit.

https://doi.org/10.1371/journal.pone.0256885.g001

Collection and processing of blood and tissue specimens

Whole blood, plasma, and peripheral blood mononuclear cells (PBMCs) were collected before each immunization and thereafter biweekly throughout the study as previously described [25, 26]. Archived peripheral blood and mononuclear cells (PBMCs) reported in this study were collected at week 10–12 for the HIV Env+adjuvant, the Co-administrationthe 3M-052-SE groups. There were no archived PBMC samples available at week 10 or 12 for the Extended Interval group and no PBMCs were sorted for this group. Archived tissue samples including spleen, lymph nodes (LNs; axillary, mesenteric, submandibular, cervical, submental, and retropharyngeal), and intestinal tissues (colon and ileum) were used for preparation of mononuclear cell (MNC) suspensions at necropsies for HIV Env+adjuvant, coadministration, and extended interval groups. Meanwhile archived LN biopsies from week 10 were used for preparation of MNC suspensions in the 3M-052-SE group. The 3M-052-SE group was part of a challenge study and therefore no necropsy tissues were available for antigen-specific B cell sorting. Summary of archived samples available for antigen-specific B cell sorting in this study are listed in Table 1.

[Figure omitted. See PDF.]

Table 1. Frequency of Env-specific memory B cells in sorted tissues and PBMCs of envelope-vaccinated infant monkeys.

https://doi.org/10.1371/journal.pone.0256885.t001

Single-cell flow cytometry sorting of antigen-specific memory B cells

Due to limited availability of archived samples, antigen-specific memory B cell sorting was performed in 15 of 19 infants (Table 1). Summary of percent antigen-specific memory B cells from archived samples of individual infants and sample types are listed in Table 1 and representative flow cytometry analyses for different sample types from each vaccine group are listed in S1A–S1J Fig. Briefly, PBMCs and tissue MNC suspensions were treated with 5 μM Chk2 inhibitor II or 2-[4-(4-chlorophenoxy)phenyl]-1H-benzimidazole-5-carboxamide (Sigma) prepared in final volume of 1% bovine serum albumin (Sigma-Aldrich) in 1X phosphate buffered saline (PBS, Sigma). PBMCs and MNC suspensions were blocked with 6.25 μg/ml anti-human CD4 antibody (BD Biosciences) at 4°C for 15 min followed by staining with a panel of fluorochrome-conjugated antibodies to identify antigen-specific memory B cells as described by gating strategy (S1A–S1J Fig). Briefly, lymphocytes were gated on singlets and live cells based on Aqua vital dye (Invitrogen), exclusion of T cells (CD3-PerCP-Cy5.5, clone SP34-2, BD Biosciences) and monocytes/macrophages (CD14-BV570, clone M5E2, BioLegend; CD16-Phycoerythrin-Cy7, clone 3G8, BD Biosciences), followed by selection of memory B cells by positive expression of CD20 and CD27 (CD20-FITC, clone 2H7; CD27-APC-Cy7, clone O323, both BioLegend) but not immunoglobulin D (IgD-PE, Southern Biotech) with double specificity to both HIV C.1086 Env hooks (BV421-gp120 C.1086 and AF647-gp120 C.1086, both generated in-house).

Polymerase chain reaction (PCR) amplification of immunoglobulin (Ig) VH and VL genes

The sorted single-cell antigen specific memory B cells VH and VL genes were amplified by nested PCR as previously described [28–30] followed by Sanger sequencing of the purified nested PCR products. Sequences were analyzed using a custom bioinformatics pipeline and were annotated with immunogenetic information using the Cloanalyst software package (https://www.bu.edu/computationalimmunology/research/software/) [31]. Identification of Ig subtypes and functional Ig heavy and light chains were first performed using human Ig sequence database as previously described [28] and recombinant monoclonal antibodies (mAbs) were generated based on this analysis. Subsequently the immunogenetic characteristics of the recombinant mAbs were reanalyzed using a rhesus Ig sequence database once it became available.

Statistical analysis

Due to the small number of previously archived samples used in this study, the statistical analysis is mostly descriptive. Antibody binding data were analyzed by a non-parametric Mann-Whitney U test with p-value of <0.05 considered as significant. Correlation analysis was performed using Spearman’s rank correlation test, with p-value of <0.05 considered as significant. All statistical analyses were performed using GraphPad Prism version 9.0.0 for Windows, GraphPad Software, San Diego, CA, USA.

Results

TLR agonist adjuvants showed enhanced pediatric HIV Env antibody responses despite similar antigen-specific memory B cell frequencies in PBMC and tissues across vaccination groups

We have evaluated different vaccine regimens and immunization schedules (Fig 1, see details in methods) to optimize HIV Env-specific antibody responses in infant rhesus macaques [25, 26]. Previously, we observed that while the peak immunogenicity across different vaccine groups was comparable the magnitude and quality of vaccine-induced HIV Env-specific responses was enhanced by increasing the timing of vaccination interval from 3 to 6 weeks [25] and by using the TLR7/8 agonist adjuvant 3M-052-SE when compared to alum or to the TLR4 ligand glucopyranosyl lipid formulated in SE (GLA-SE) [26]. Since these analyses were conducted in two separate studies [25, 26], here we compare the responses at peak immunogenicity among all vaccine groups. Overall, infants vaccinated with 3M-052-SE-adjuvanted vaccine developed the highest magnitude of HIV Env C.1086 gp120-specific plasma IgG antibody concentrations across all four groups at peak immunogenicity (Fig 2A). Additionally, at peak immunogenicity, infant plasma antibody from all four groups developed varying levels of cross clade Env gp120 responses as well as broad heterologous epitope specificities and breadth (S2A and S2B Fig). Infant plasma IgG antibody from group 4 (3M-052-SE) also demonstrated higher avidity strengths against the tested antigens (HIV Env 1086d7gp120 K160N and gp70 ConC V3) compared to other immunization groups (S2C Fig).

[Figure omitted. See PDF.]

Fig 2. Vaccine-induced Env binding IgG responses against HIV 1086c gp120 K160N does not correlate with frequency of antigen-specific memory B cells in infants.

(A) Binding antibody level was significantly higher in the 3M-052-SE group as measured by ELISA at peak immunogenicity. (B) Binding antibody levels at peak immunogenicity was not correlated with antigen-specific memory B cell frequencies in all vaccine groups. Statistical analyses were performed with GraphPad Prism, * denoted significant p-values of <0.05 by a non-parametric Mann-Whitney test. Correlation analysis was performed using Spearman’s rank correlation test with GraphPad Prism, with p-value of <0.05 considered as significant.

https://doi.org/10.1371/journal.pone.0256885.g002

Single-cell flow cytometry sorting of antigen-specific memory B cells indicated low level of frequencies across all vaccination groups

To determine whether the different vaccine regimens induced distinct frequency of antigen-specific memory B cells, we characterized antigen-specific memory B cells from PBMC and tissues at selected time points (Fig 1). We were able to obtain Env-specific memory B cells (CD3-CD16-CD14-CD20+CD27+IgD-, double positive for HIV Env C.1086 gp120) from a total of 15 infants, including 3 of 5 infants in group 1 (HIV Env+adjuvant), 5 of 5 infants in group 2 (Co-administration), 3 of 5 infants in group 3 (Extended Interval), and from 4 of 4 infants in group 4 (3M-052-SE) (Table 1). The frequency of Env-specific memory B cells were low across all vaccine groups (0.06–1.72%) with 0.07–1.72% range in group 1, 0.06–0.97% in group 2, 0.08–0.49% in group 3, and 0.06–0.11% in group 4 (Table 1). The frequency of antigen-specific memory B cells were not correlated with the Env-specific IgG responses against the HIV Env 1086c gp120 K160N, although the small sample size could be a major limitation of the analysis (Fig 2B).

B cell receptor (BCR) repertoires analysis using the human Ig-gene sequence database led to identification of functional heavy- and light-chain pairs, which were used to produce a total of 39 monoclonal antibodies (mAbs). The human Ig-gene sequence database was utilized as the rhesus Ig-gene database was unavailable at the time when the initial analyses were performed [28]. The 39 mAbs came from 13 specimens including 3 specimens from group 1 (2 infants), 3 specimens from group 2 (2 infants), 4 specimens from group 3 (2 infants), and 3 specimens from group 4 (3 infants) for final B cell repertoire analyses (S1 Table). No tissue specimen was available for group 4 as these infants were part of a challenge study. Overall, the frequencies of Env-specific memory B cells in PBMCs and tissues did not differ between the vaccine groups. The percent frequency of antigen-specific memory B cells ranged from 0.06 to 1.72%, with the highest frequency observed in the PBMCs of 2 infants in group 1 (HIV Env+adjuvant) (Table 1). Notably, despite higher magnitude and quality of Env-specific plasma antibody responses in group 3 (Extended Interval) compared to group 2 (Co-administration) [25], the frequencies of Env-specific memory B cells did not differ between these groups. Similarly, despite higher binding magnitude and avidity strength of Env-specific plasma antibody responses in group 4 (3M-052-SE) when compared to others, the frequencies of Env-specific memory B cells at least in the PBMCs of vaccinated infants in this group did not differ to the other groups. Altogether, these data suggest that the size of HIV Env vaccine-elicited memory B cell pool is not directly related to magnitude or quality of plasma Env-specific antibody responses in infant rhesus macaques.

Epitope specificities and immunogenetics of vaccine-elicited Env-specific memory BCR repertoires did not differ across vaccination groups

We next characterized the immunogenetics and epitope specificities of the HIV Env vaccine-elicited mAbs. In our initial analysis, the identification of functional heavy- and light-chains and their immunogenetic characteristics was conducted using a previously validated bioinformatic method for rhesus BCR repertoire characterization using the human Ig sequence database [28, 30] (S1 Table). Notably, the majority of identified functional heavy- and light-chain Ig pairs were of IgG isotype (31 functional pairs) with IgA and IgM isotypes represented in 4 and 3 functional pairs, respectively (S3 Fig and S1 Table). Most mAbs were specific to the HIV Env V3 loop (14 functional Ig pairs) followed by V1V2 loop-specific mAbs (5 functional pairs). Only one functional pair targeting the CD4 binding site was identified, and we were unable to determine the epitope specificities of 18 functional heavy and light chains pairs. These data are consistent with our prior observations which demonstrated that most of the polyclonal plasma antibody response in these vaccinated infant macaques was directed against the HIV Env V3 region [25, 26].

Recent advances in genome sequencing and detailed characterization of rhesus Ig loci provided better understanding of allelic diversity in rhesus Ig genes [32–34]. Moreover, the availability of rhesus Ig gene libraries [32] provided the opportunity to reanalyze the 39 functional Ig pairs that we previously identified using the human Ig sequence database. This secondary analysis confirmed that 26 of 39 Ig pairs were indeed functional (Table 2) with 5 Ig pairs in group 1 (HIV Env+adjuvant), 7 Ig pairs in group 2 (Co-administration), 8 Ig pairs in group 3 (Extended Interval), and 6 Ig pairs in group 4 (3M-052-SE) (Fig 3). Similar to earlier analysis using the human Ig sequence database, the majority of identified functional heavy- and light-chain Ig pairs are of IgG isotype (21 functional pairs). Based on the rhesus Ig sequence database, the variable heavy chain (VH) gene usage was largely restricted to the VH3 and VH4 gene families, across all the vaccination groups, and there was no apparent difference in VH usage in the 3M-052-SE group as compared to the other vaccination groups (Table 2). The majority of the Ig pairs confirmed as functional in the secondary analysis using rhesus Ig-gene database were against undefined epitopes (12 Ig pairs), whereas 6 pairs were specific to the V3 loop and 5 pairs targeted the V1V2 loop. this suggests that an important proportion of the vaccine-elicited mAbs were against non-linear conformational epitopes on the HIV Env.

[Figure omitted. See PDF.]

Fig 3. Analyses of epitope specificity and immunogenetic characteristics of the Env-specific functional heavy- and light-chains of 26 vaccine-elicited mAbs in infants using rhesus Ig-gene database.

Reanalysis using newly developed software based on rhesus macaque Ig sequences confirmed that 26 of 39 heavy- and light-chain pairs previously identified using human Ig sequences were functional. Overall, the epitope specificities, VH gene family usage, and isotype distribution were similar across different vaccina groups. Epitope specificity, VH gene family usage, and isotype distribution of identified functional heavy and light chains are displayed in concentric circles. The number of mAbs per group is displayed in the center.

https://doi.org/10.1371/journal.pone.0256885.g003

[Figure omitted. See PDF.]

Table 2. Immunogenetic characteristics of isolated envelope (Env)-reactive mAbs of Env-vaccinated infant monkeys based on rhesus macaque immunoglobulin database analysis.

https://doi.org/10.1371/journal.pone.0256885.t002

Frequency of somatic hypermutation and heavy chain complementarity-determining region 3 (HCDR3) length in HIV Env-specific Ig pairs did not differ across vaccination groups

We also sought to define whether increased somatic hypermutation and affinity maturation in Env-specific Ig pairs could potentially contributed to distinct plasma antibody responses induced by different pediatric HIV vaccine regimens [25, 26]. Based on the preliminary Ig sequence analysis using the human Ig gene database, the frequency of somatic hypermutation (SHM) was comparable across the vaccination groups (S4 Fig) with range of SMH of 5.31–10.03% in group 1 (HIV Env+adjuvant), 4.52–10.15% in group 2 (Co-administration), 6.05–11.08% in group 3 (Extended Interval), and 4.15–7.03% in group 4 (3M-052-SE) (S1 Table). Additionally, the HCDR3 of Env-reactive functional pairs were comparable across all vaccination groups, with HCDR3 length range of 10–23 aa. The highest median HCDR3 length was observed in the HIV Env+adjuvant group (median length of 23 aa) and the lowest was observed in the 3M-052-SE group (median length of 13 aa). The HCDR3 median length in group 2 (Co-administration) and 3 (Extended Interval) was 14 aa and 17 aa, respectively.

Comparable SHM frequency across the groups was confirmed with the reanalysis of the human Ig sequence database using the rhesus Ig sequence database (Fig 4 and Table 2), albeit mutation rates were slightly lower compared to the initial analysis using the human Ig sequence database. The range of SHM was 3.82–7.9% in group 1 (HIV Env+adjuvant), 4.17–7.99% in the group 2 (Co-administration), 1.39–13.75% in group 3 (Extended Interval), and 0–5.15% in group 4 (3M-052-SE) (Table 2). The HCDR3 length of Env-reactive functional pairs was also comparable across all vaccination groups. Similar to human Ig sequence, the highest median HCDR3 region was observed in group 1 (HIV Env+adjuvant) and the lowest median HCDR3 region was observed in group 4 (3M-052-SE). Interestingly, the median HCDR3 length in group 2 (Co-administration) was lower based on human Ig sequence database (median HCDR3 length of 14 aa) than the rhesus Ig sequence (median HCDR3 length of 17 aa). Meanwhile, the median HCDR3 length for group 3 (Extended Interval) was lower with the rhesus IgG sequence (median HCDR3 length of 15 aa) than in the human Ig sequence (median HCDR3 length of 17 aa). These findings highlight the limitations of using human database to analyze the BCR repertoire in rhesus.

[Figure omitted. See PDF.]

Fig 4. Frequency of somatic hypermutation and heavy chain complementarity-determining region 3 (HCDR3) length of vaccine-elicited Env-reactive functional heavy- and light-chains identified using rhesus Ig sequence database.

(A) Analysis of percent somatic hypermutation frequency and HCDR3 lengths for Env-reactive heavy and light chains pairs (39 mAb pairs) from infant antigen-specific B cells based on human immunoglobulin (Ig) sequence database. (B) Analysis of percent somatic hypermutation frequency and HCDR3 lengths for Env-reactive heavy- and light-chains functional pairs (26 mAb pairs) from vaccinated infant antigen-specific memory B cells using rhesus Ig-gene database. Horizontal lines indicated median values of individual groups. Corresponding functional heavy and light chains isolated from individual infants are denoted by symbols.

https://doi.org/10.1371/journal.pone.0256885.g004

Altogether, our data suggest that the magnitude and quality of vaccine-elicited plasma Env-specific antibody responses administered with different adjuvants are not related to the size of the antigen-specific memory B cell pool or to the immunogenetics characteristics of the vaccine-elicited Ig pairs.

Discussion

Based on previous observations that different pediatric HIV vaccine regimens induced distinct plasma antibody responses [25, 26], the goal of this study was to investigate whether these responses were driven by distinct memory BCR repertoire characteristics including distinct Ig gene usage, rates of SHM and HCDR3 lengths. We utilized samples from 2 completed immunization studies [25, 26], in which newborn infant rhesus macaques were immunized with four distinct vaccine regimens (Fig 1). Our results indicate that although the different vaccine regimens induced distinct plasma antibody responses, there were no significant differences in the B cell repertoires.

The development of a safe and effective pediatric HIV-1 vaccine to eliminate postnatal infant HIV-1 infections will probably require the use of novel adjuvants such as Toll-like receptor (TLR) agonists. Indeed, recent studies have demonstrated that TLR7/8 adjuvantation can enhance vaccine responses in early life [35–37]. Notably, our group previously reported that an HIV vaccine adjuvanted with 3M-052-SE induces superior plasma antibody levels than other adjuvanted HIV regimens in infant rhesus macaques [26]. We showed that antibody levels against autologous and heterologous envelope proteins, and observed that overall, the magnitude of the vaccine-elicited antibody response was higher in the 3M-052-SE group (Fig 2A and S2A–S2C Fig). However, the impact of these different pediatric HIV Env vaccine regimens on the developing infant antigen-specific B cell BCR repertoires is still unclear.

Previous study investigating the BCR repertoire in HIV immunized adult monkeys have reported a preferential usage of the VH4 and VH3 families [28]. Our initial analysis of the BCR repertoire in the vaccinated infant rhesus macaques was conducted using the same human Ig sequence database and bioinformatics methods, in which we similarly observed that 69% of the infant Env-specific functional heavy and light chains in this study used VH4 genes and 26% use VH3 genes (S1 Table). However, we found that the Env-specific functional heavy and light chains in infant rhesus macaques were slightly shorter and had lower mutation rates than in adults. For example, in adult rhesus macaques immunized with a pox prime/protein boost regimen the median HCDR3 length after the fifth immunization was 16 aa and the SHM rate was 9.3%. Meanwhile, in our study the overall median HCDR3 length across all vaccination groups was 14 aa and the SMH rate was 6.3%. The lower levels of SHM rate in infant rhesus macaques as compared to adult rhesus macaques is in accordance with the observation that SHM rate increases with age. Moreover, broadly neutralizing antibodies (bNAbs) isolated from HIV-infected pediatric patients (1-year post-infection) also have lower levels of SHM rate than HIV-infected adult bNAbs directed against the same epitopes [38].

Due to their close genetic similarity to humans, rhesus macaques are valuable animal model for studies of infectious diseases including the understanding of vaccine-elicited immune responses. However, despite their wide usage as a human surrogate model system, many aspects of the rhesus macaque immune system are still under characterized and poorly annotated. Recent advances in high-throughput NGS and specialized computational methods provide tools to the scientific community to compare rhesus repertoires of heavy and light-chains to humans in order to better understand how they may perform as a model system for B- and T-cell mediated immunity in humans [39–41]. Characterization of human and rhesus BCR repertoires showed that the frequency of V- and J-gene segment usage and HCDR3 lengths between human and rhesus were in concordance with one another [39, 42]. However, comparative analyses of different Ig subtypes (IgM, IgG, IgK and IgL) sequences revealed significant differences in the overall BCR repertoires. Rhesus macaques have higher diversity of BCR repertoires with different family gene usage and slight difference in the frequencies of HCDR3 lengths within the IgM+ BCR repertoires, likely due to gene family usage in the class-switched (IgG) compartment [39]. However, importantly, the low abundance of long CDRH3s in rhesus IgM+ B cells did not impede their expansion into the IgG+ B cells, and in rhesus IgG+ B cell frequency was comparable to human IgG+ B cells. Thus, given the complexity of gene recombination, high diversity in rhesus BCR repertoires, and the close genetic relationship between rhesus and humans, it is likely that rhesus B-cell compartment recapitulates its human counterpart and is poised to respond to antigen in similar manners.

Historically, a limitation of BCR repertoire analysis in the nonhuman primate (NHP) model has been the lack of rhesus macaque Ig sequence database. Recent development of bioinformatics tools that enable comparison of Ig sequences from immunized animals to a NHP reference database allows for a more accurate characterization of the BCR repertoire in response to vaccination or infection. We found some differences in BCR repertoire characteristics including VH usage, SMH rate, and Ig subclass/isotypes when the same Ig pairs were analyzed using human Ig and rhesus Ig sequence databases (Figs 3 and 4, Table 2, S3 and S4 Figs, S1 Table). Notably only 26 of 39 pairs identified as functional based on the human Ig sequence database were confirmed to be functional using the rhesus Ig sequence database. This could be due to high sequence homology using the rhesus Ig sequence database compared to human Ig sequence database. These findings highlight the importance of species-specific database for comprehensive understanding of the BCR repertoire and antibody maturation in response to vaccinations and/or infections using the rhesus macaque model.

Interestingly, the observed SHM rate in infants across all vaccinations groups in this study (4.7%) is only slightly higher than the observed SHM rate in human HIV vaccine trials such as GSK PRO HIV-002 (3.8%) and the RV144 trial (2.4%). HIV-1-infected infants have been shown to develop neutralization breadth earlier than HIV-infected adults. Plasma antibody responses in HIV-infected infants neutralized a panel of diverse HIV-1 viruses, including more difficult to neutralize cross clade variants [43]. Additionally, these responses were observed as early as 1 to 2 years post-infection. Furthermore, bNAbs isolated from pediatric HIV cases appeared to have lower levels of SHM when compared to bNAbs isolated from HIV-infected adults [38]. Altogether, these data suggest that induction of HIV-1-specific plasma antibody neutralization can be achieved in children without prolonged extensive SHM and affinity maturation.

Conclusion

Our results suggest that the high plasma antibody magnitude and functionality achieved with 3M-052-SE adjuvantation is not accompanied by distinct B cell repertoire characteristics. Nevertheless, our study has several limitations including, (1) small sample size, (2) distinct vaccine regimens and timing used, and (3) the low number of functional heavy and light chain pairs identified from the different vaccination groups. Nevertheless, this study provides novel insights into (1) the B cell repertoires in vaccinated infant rhesus macaques, which has not been analyzed, and (2) significant increase of antibody magnitude observed in our previous studies in the 3M-052-SE group was not related to distinct immunogenetics characteristics of vaccine-induced antibodies. Further investigation of the mechanism by which the 3M-052-SE adjuvant leads to enhanced immune responses in the setting of the developing early life immune system are warranted. Additionally, it will be important to evaluate whether the enhanced immunogenicity of HIV vaccines with 3M-052-SE adjuvantation in infants is associated with protection from oral virus exposure.

Supporting information

S1 Fig.

A-J. Representative flow cytometry analyses for single cell HIV Env-specific memory B cell sorting from PBMC and tissues of infants by vaccination groups. Percentage of memory B cells (CD20+CD27+IgD-), which were reactive to HIV 1086c gp120 K160N tagged with two colors-BV421 and AF647 (double positive cells) indicated in the upper right quadrant. Gates were drawn based on isotype controls and fluorescent minus one (FMO) controls.

https://doi.org/10.1371/journal.pone.0256885.s001

(PDF)

S2 Fig. Characterization of vaccination-induced plasma envelope IgG responses in infant rhesus macaques.

(A) Cross clade HIV gp120-specific IgG responses at peak immunogenicity by luminex assay (BAMA). (B) Variable and conserved epitopes-specific IgG responses at peak immunogenicity by BAMA. (C) The 3M-052 group antibody responses had higher avidity strength against most of the tested antigens. Avidity 1/k off (the inverse of the dissociation rate) was plotted as a measure of the strength of binding and avidity scores which take into consideration the magnitude are also shown. Statistical analyses were performed with GraphPad Prism, * denoted significant p-values of <0.05 by a non-parametric Mann-Whitney test.

https://doi.org/10.1371/journal.pone.0256885.s002

(PDF)

S3 Fig. Analyses of epitope specificity and immunogenetic characteristics of the Env-specific functional heavy- and light-chains of 39 vaccine-elicited mAbs in infants using human Ig-gene database.

Initial analysis with human immunoglobulin (Ig) database indicated a total of 39 heavy- and light-chain pairs isolated from antigen-specific memory B cells across different vaccine groups. Epitope specificity, VH gene family usage, and isotype distribution of identified functional heavy- and light-chain pairs were similar across vaccine groups. Epitope specificity, VH gene family usage, and isotype distribution of identified functional heavy and light chains are displayed in concentric circles. The number of mAbs per group is displayed in the center.

https://doi.org/10.1371/journal.pone.0256885.s003

(PDF)

S4 Fig. Frequency of somatic hypermutation and heavy chain complementarity-determining region 3 (HCDR3) length of vaccine-elicited Env-reactive functional heavy- and light-chains identified using rhesus Ig sequence database.

Analysis of percent somatic hypermutation frequency and HCDR3 lengths for Env-reactive heavy and light chains pairs (39 mAb pairs) from infant antigen-specific B cells based on human immunoglobulin (Ig) sequence database. Horizontal lines indicated median values of individual groups. Corresponding functional heavy and light chains isolated from individual infants are denoted by symbols.

https://doi.org/10.1371/journal.pone.0256885.s004

(PDF)

S1 Table. Immunogenetic characteristics of isolated envelope (Env)-reactive mAbs of Env-vaccinated infant monkeys based on human immunoglobulin database analysis.

A total of 39 pairs of potentially Env-reactive mAbs were isolated from the four vaccination groups across several anatomic compartments. Frequency of gene usage, percent somatic hypermutation, and complementarity-region 3 (CDR3) length are displayed for the heavy and light chains for each mAb along with the isotype and epitope specificity.

https://doi.org/10.1371/journal.pone.0256885.s005

(PDF)

S1 Data.

https://doi.org/10.1371/journal.pone.0256885.s006

(PDF)

Acknowledgments

We would like to thank Carolyn Weinbaum for her contribution in animal study and sample coordination. Flow cytometry and single cell sorting was performed at the Duke Human Vaccine Institute (DHVI) Flow Cytometry Facility, Duke University, Durham, NC.

Citation: Berendam SJ, Morgan-Asiedu PK, Mangan RJ, Li SH, Heimsath H, Luo K, et al. (2021) Different adjuvanted pediatric HIV envelope vaccines induced distinct plasma antibody responses despite similar B cell receptor repertoires in infant rhesus macaques. PLoS ONE 16(12): e0256885. https://doi.org/10.1371/journal.pone.0256885

References

1. Organization WH. HIV/AIDS 2019 [https://www.who.int/news-room/fact-sheets/detail/hiv-aids.

2. Haas AD M M, Egger M, Tenthani L, Tweya H, Jahn A, Gadabu OJ, et al. Adherence to Antiretroviral Therapy During and After Pregnancy: Cohort Study on Women Receiving Care in Malawi’s Option B+ Program. Clin Infect Dis. 2016;63(9):1227–35. pmid:27461920

3. Tenthani L H A, Tweya H, Jahn A, van Oosterhout JJ, Chimbwandira F, Chirwa Z, et al. Retention in care under universal antiretroviral therapy for HIV-infected pregnant and breastfeeding women (’Option B+’) in Malawi. AIDS. 2014;28(4):589–98. pmid:24468999

4. Stadeli KM R D. Rates of emergence of HIV drug resistance in resource-limited settings: a systematic review. Antivir Ther. 2015;18(1):115–23.

5. UNAIS. Prevention Gap Report. Joint United Nations Programme on HIV/AIDS, Geneva, Switzerland. 2016.

6. Drake AL W A, Richardson B, John-Stewart G. Incident HIV during pregnancy and postpartum and risk of mother-to-child HIV transmission: a systematic review and meta-analysis. PloS Med. 2014;11. pmid:24586123

7. WHO. Global Update on the Health Sector Response to HIV. World Health Organization, Geneva, Switzerland. 2014.

8. Siberry GK. Preventing and managing HIV infection in infants, children, and adolescents in the United States. Pediatr Rev. 2014;35(7):268–86. pmid:24986927

9. Gao Y, McKay PF, Mann JFS. Advances in HIV-1 Vaccine Development. Viruses. 2018;10(4). pmid:29614779

10. Dolin R, Graham BS, Greenberg SB, Tacket CO, Belshe RB, Midthun K, et al. The safety and immunogenicity of a human immunodeficiency virus type 1 (HIV-1) recombinant gp160 candidate vaccine in humans. NIAID AIDS Vaccine Clinical Trials Network. Ann Intern Med. 1991;114(2):119–27. pmid:1984386

11. Wright PF, Lambert JS, Gorse GJ, Hsieh RH, McElrath MJ, Weinhold K, et al. Immunization with envelope MN rgp120 vaccine in human immunodeficiency virus-infected pregnant women. J Infect Dis. 1999;180(4):1080–8. pmid:10479134

12. Belshe RB, Stevens C, Gorse GJ, Buchbinder S, Weinhold K, Sheppard H, et al. Safety and immunogenicity of a canarypox-vectored human immunodeficiency virus Type 1 vaccine with or without gp120: a phase 2 study in higher- and lower-risk volunteers. J Infect Dis. 2001;183(9):1343–52. pmid:11294665

13. Rerks-Ngarm S, Pitisuttithum P, Nitayaphan S, Kaewkungwal J, Chiu J, Paris R, et al. Vaccination with ALVAC and AIDSVAX to prevent HIV-1 infection in Thailand. N Engl J Med. 2009;361(23):2209–20. pmid:19843557

14. Robb ML, Rerks-Ngarm S, Nitayaphan S, Pitisuttithum P, Kaewkungwal J, Kunasol P, et al. Risk behaviour and time as covariates for efficacy of the HIV vaccine regimen ALVAC-HIV (vCP1521) and AIDSVAX B/E: a post-hoc analysis of the Thai phase 3 efficacy trial RV 144. Lancet Infect Dis. 2012;12(7):531–7. pmid:22652344

15. Haynes BF, Gilbert PB, McElrath MJ, Zolla-Pazner S, Tomaras GD, Alam SM, et al. Immune-correlates analysis of an HIV-1 vaccine efficacy trial. N Engl J Med. 2012;366(14):1275–86. pmid:22475592

16. Release. NIoHN. Experimental HIV Vaccine Regimen Ineffective in Preventing HIV. 3 February 2020.

17. Lin L, Finak G, Ushey K, Seshadri C, Hawn TR, Frahm N, et al. COMPASS identifies T-cell subsets correlated with clinical outcomes. Nat Biotechnol. 2015;33(6):610–6. pmid:26006008

18. Johnson DC, McFarland EJ, Muresan P, Fenton T, McNamara J, Read JS, et al. Safety and immunogenicity of an HIV-1 recombinant canarypox vaccine in newborns and infants of HIV-1-infected women. J Infect Dis. 2005;192(12):2129–33. pmid:16288378

19. McFarland EJ, Johnson DC, Muresan P, Fenton T, Tomaras GD, McNamara J, et al. HIV-1 vaccine induced immune responses in newborns of HIV-1 infected mothers. AIDS. 2006;20(11):1481–9. pmid:16847402

20. Kaleebu P, Njai HF, Wang L, Jones N, Ssewanyana I, Richardson P, et al. Immunogenicity of ALVAC-HIV vCP1521 in infants of HIV-1-infected women in Uganda (HPTN 027): the first pediatric HIV vaccine trial in Africa. J Acquir Immune Defic Syndr. 2014;65(3):268–77. pmid:24091694

21. Kintu K, Andrew P, Musoke P, Richardson P, Asiimwe-Kateera B, Nakyanzi T, et al. Feasibility and safety of ALVAC-HIV vCP1521 vaccine in HIV-exposed infants in Uganda: results from the first HIV vaccine trial in infants in Africa. J Acquir Immune Defic Syndr. 2013;63(1):1–8. pmid:23221981

22. Fouda GG, Cunningham CK, McFarland EJ, Borkowsky W, Muresan P, Pollara J, et al. Infant HIV type 1 gp120 vaccination elicits robust and durable anti-V1V2 immunoglobulin G responses and only rare envelope-specific immunoglobulin A responses. J Infect Dis. 2015;211(4):508–17. pmid:25170104

23. Van Rompay KK, Greenier JL, Cole KS, Earl P, Moss B, Steckbeck JD, et al. Immunization of newborn rhesus macaques with simian immunodeficiency virus (SIV) vaccines prolongs survival after oral challenge with virulent SIVmac251. J Virol. 2003;77(1):179–90. pmid:12477823

24. Van Rompay KK, Abel K, Lawson JR, Singh RP, Schmidt KA, Evans T, et al. Attenuated poxvirus-based simian immunodeficiency virus (SIV) vaccines given in infancy partially protect infant and juvenile macaques against repeated oral challenge with virulent SIV. J Acquir Immune Defic Syndr. 2005;38(2):124–34. pmid:15671796

25. Phillips B, Fouda GG, Eudailey J, Pollara J, Curtis AD 2nd, Kunz E, et al. Impact of Poxvirus Vector Priming, Protein Coadministration, and Vaccine Intervals on HIV gp120 Vaccine-Elicited Antibody Magnitude and Function in Infant Macaques. Clin Vaccine Immunol. 2017;24(10).

26. Phillips B, Van Rompay KKA, Rodriguez-Nieves J, Lorin C, Koutsoukos M, Tomai M, et al. Adjuvant-Dependent Enhancement of HIV Env-Specific Antibody Responses in Infant Rhesus Macaques. J Virol. 2018;92(20). pmid:30089691

27. Dowling DJ, van Haren SD, Scheid A, Bergelson I, Kim D, Mancuso CJ, et al. TLR7/8 adjuvant overcomes newborn hyporesponsiveness to pneumococcal conjugate vaccine at birth. JCI Insight. 2017;2(6):e91020. pmid:28352660

28. Luo K, Liao HX, Zhang R, Easterhoff D, Wiehe K, Gurley TC, et al. Tissue memory B cell repertoire analysis after ALVAC/AIDSVAX B/E gp120 immunization of rhesus macaques. JCI Insight. 2016;1(20):e88522. pmid:27942585

29. Liao HX, Levesque MC, Nagel A, Dixon A, Zhang R, Walter E, et al. High-throughput isolation of immunoglobulin genes from single human B cells and expression as monoclonal antibodies. J Virol Methods. 2009;158(1–2):171–9. pmid:19428587

30. Wiehe K, Easterhoff D, Luo K, Nicely NI, Bradley T, Jaeger FH, et al. Antibody light-chain-restricted recognition of the site of immune pressure in the RV144 HIV-1 vaccine trial is phylogenetically conserved. Immunity. 2014;41(6):909–18. pmid:25526306

31. Kepler TB. Reconstructing a B-cell clonal lineage. I. Statistical inference of unobserved ancestors. F1000Res. 2013;2:103. pmid:24555054

32. Ramesh A, Darko S, Hua A, Overman G, Ransier A, Francica JR, et al. Structure and Diversity of the Rhesus Macaque Immunoglobulin Loci through Multiple De Novo Genome Assemblies. Front Immunol. 2017;8:1407. pmid:29163486

33. Cirelli KM, Carnathan DG, Nogal B, Martin JT, Rodriguez OL, Upadhyay AA, et al. Slow Delivery Immunization Enhances HIV Neutralizing Antibody and Germinal Center Responses via Modulation of Immunodominance. Cell. 2019;177(5):1153–71 e28. pmid:31080066

34. Vazquez Bernat N, Corcoran M, Nowak I, Kaduk M, Castro Dopico X, Narang S, et al. Rhesus and cynomolgus macaque immunoglobulin heavy-chain genotyping yields comprehensive databases of germline VDJ alleles. Immunity. 2021;54(2):355–66 e4. pmid:33484642

35. Holbrook BC, Aycock ST, Machiele E, Clemens E, Gries D, Jorgensen MJ, et al. An R848 adjuvanted influenza vaccine promotes early activation of B cells in the draining lymph nodes of non-human primate neonates. Immunology. 2018;153(3):357–67. pmid:28940186

36. Holbrook BC, D’Agostino RB Jr., Tyler Aycock S, Jorgensen MJ, Hadimani MB, Bruce King S, et al. Adjuvanting an inactivated influenza vaccine with conjugated R848 improves the level of antibody present at 6months in a nonhuman primate neonate model. Vaccine. 2017;35(45):6137–42. pmid:28967521

37. Holbrook BC, Kim JR, Blevins LK, Jorgensen MJ, Kock ND, D’Agostino RB Jr., et al. A Novel R848-Conjugated Inactivated Influenza Virus Vaccine Is Efficacious and Safe in a Neonate Nonhuman Primate Model. J Immunol. 2016;197(2):555–64. pmid:27279374

38. Simonich CA, Williams KL, Verkerke HP, Williams JA, Nduati R, Lee KK, et al. HIV-1 Neutralizing Antibodies with Limited Hypermutation from an Infant. Cell. 2016;166(1):77–87. pmid:27345369

39. Vigdorovich V, Oliver BG, Carbonetti S, Dambrauskas N, Lange MD, Yacoob C, et al. Repertoire comparison of the B-cell receptor-encoding loci in humans and rhesus macaques by next-generation sequencing. Clin Transl Immunology. 2016;5(7):e93. pmid:27525066

40. Brochu HN, Tseng E, Smith E, Thomas MJ, Jones AM, Diveley KR, et al. Systematic Profiling of Full-Length Ig and TCR Repertoire Diversity in Rhesus Macaque through Long Read Transcriptome Sequencing. J Immunol. 2020;204(12):3434–44. pmid:32376650

41. Zhang W, Li X, Wang L, Deng J, Lin L, Tian L, et al. Identification of Variable and Joining Germline Genes and Alleles for Rhesus Macaque from B Cell Receptor Repertoires. J Immunol. 2019;202(5):1612–22. pmid:30700589

42. Dai K, He L, Khan SN, O’Dell S, McKee K, Tran K, et al. Rhesus Macaque B-Cell Responses to an HIV-1 Trimer Vaccine Revealed by Unbiased Longitudinal Repertoire Analysis. mBio. 2015;6(6):e01375–15. pmid:26530382

43. Goo L, Chohan V, Nduati R, Overbaugh J. Early development of broadly neutralizing antibodies in HIV-1-infected infants. Nat Med. 2014;20(6):655–8. pmid:24859529