Phage display demonstrates durable differences in serological profile by route of inoculation in primary infections of non-human primates with Dengue Virus 1

Natural dengue virus (DENV) infections occur by mosquito bite but how the inoculation route affects the humoral immune response is unknown. We serologically profiled 20 non-human primates (NHP) from a prior study of DENV1 infection where animals were inoculated by mosquito (N = 10) or subcutaneous injection (N = 10). Using a comprehensive, densely tiled and highly redundant pan-flavivirus programmable phage library containing 91,562 overlapping 62 amino acid peptides, we produced a high-resolution map of linear peptide sequences enriched during DENV seroconversion. Profiles in mosquito-inoculated and subcutaneously-inoculated animals were similar up to 90 days after primary infection, but diverged at 1 year with differences in sero-reactivity in the Envelope (E; residues 215–406; p < 0.08), and Nonstructural-3 (NS3; residues 549–615; p < 0.05) proteins in mosquito-inoculated versus subcutaneously-inoculated animals. Within the E protein, residues 339–384 in domain III accounted for > 99% of the observed sero-reactivity difference. Antibody breadth did not vary by mode of inoculation. The differential reactivity to E domain III seen by phage display validated orthogonally by ELISA, but did not correlate with late neutralization titers. Serological profiling of humoral immune responses to DENV infection in NHP by programmable phage display demonstrated durable differences in sero-reactivity by route of inoculation.


Results
Inoculation of NHP with DENV1 by mosquito or by subcutaneous injection induces similar humoral responses. The programmable phage library used here consists of 91,562 overlapping 62 amino acid peptides derived from 74 flaviviruses (see "Methods"). Next generation sequencing (NGS) of the packaged phage display library confirmed > 96% representation of designed peptides, with 96% of them present at ≤ 1 read per 100,000 reads (Supplemental Figure 1). The experimental procedure for immunoprecipitations is depicted in Fig. 1A. For the 10 non-human primates (NHP) inoculated by mosquito, the number of enriched peptides relative to the pre-infection baseline sample at 7, 35, 90 and 365 days post-infection (p.i.) ranged from 22 to 344, 111 to 416, 22 to 540, and 28 to 501. For the 10 NHP inoculated by subcutaneous injection of DENV1 the number of enriched peptides at 7, 35, 90 and 365 days p.i. ranged from 26 to 93, 88 to 930, 64 to 789 and 31 to 786. Among enriched peptides, there was no bias towards peptides that were more abundant in the phage library (Supplemental Figure 2).
The proportion of DENV1-specific peptides, defined as peptides that were ≥ 90% similar to DENV1 but not to any other DENV serotype, peaked at 35 days p.i. and decreased at 90 and 365 days p.i. (Fig. 1B, left). The difference in the median proportions of DENV1-specific peptides for each of these timepoints as compared to 7 days p.i. was statistically significant (Mann-Whitney U test) and this was true regardless of the route of inoculation. At 7, 35, 90 and 365 days p.i., the median proportions of DENV1-specific peptides in the mosquitoinoculated group were 6 Cross-reactive peptides, defined here as peptides that were ≥ 90% similar to DENV1 as well as at least one other DENV sero-type or as peptides that were < 90% similar to DENV1, showed a different pattern, decreasing from day 7 p.i. as the proportion of DENV1-specific antibodies increased. As with DENV1-specific peptides, this pattern did not differ by route of inoculation. In the mosquito-inoculated group, the median proportion of cross-reactive peptides at 7 Primary infection by mosquito induces differences in serological profile compared to subcutaneous inoculation at one year post-infection. Minimal sero-reactivity was observed at 7 days postinfection, but by 35 days post-infection, the median antibody coverage for each group of 10 animals enriched for 17.3% (mosquito) and 14.3% (subcutaneous) of the DENV1 proteome respectively (Fig. 2). This coverage focused on the envelope (E), nonstructural-1 (NS1), nonstructural-3 (NS3) and nonstructural-5 (NS5) proteins. In the NHP inoculated by mosquito, the E, NS1, NS3 and NS5 proteins represented 44.2%, 20 At 365 days post-infection, the median antibody coverage targeted 16.8% and 7.4% of the DENV1 proteome in NHP inoculated by mosquito and subcutaneous routes. Of the enriched regions, the E, NS1, NS3 and NS5 proteins represented 46.5%, 13.4%, 14.6%, and 25.5% in mosquito-inoculated NHP versus 62.1%, 8.3%, 0% and 29.6% in subcutaneously inoculated NHP. The area under the curve (AUC) was calculated as a measure of antibody breadth and magnitude and showed distinct patterns at this late timepoint by route of inoculation (Fig. 2, bottom). We observed differences in the AUC spanning residues 549-615 of the NS3 protein with a larger AUC in the mosquito-inoculated group (log 10 [AUC msq − AUC sc ] = 4.03, p < 0.05). We obtained similar results in the E protein, with a region spanning residues 1-91 (domain I) having a larger AUC in the subcutaneously-inoculated group (log 10 [AUC sc − AUC msq ] = 5.62, p < 0.075) and another region spanning residues 215-406 (domain III) where the AUC was larger in the mosquito-infected group (log 10 [AUC msq − AUC sc ] = 5.02, p < 0.075). For the latter region, 97.6% of this AUC difference (log 10 [AUC msq − AUC sc ] = 4.90, p < 0.1) was accounted for by residues 339-384 in domain III (Supplemental Figure 5). While this region was present in both mosquito-and subcutaneously infected groups at 35 and 90 days post-infection, by 365 days post-infection, as indicated by AUC difference, it had waned substantially in the subcutaneously infected group. A multiple sequence alignment of reference sequences of the E protein of each of the four DENV serotypes showed that 34.7% of the positions within this region are 100% conserved across the 4 serotypes (Supplemental Figure 3). Route of inoculation does not affect breadth of antibodies targeting ED III. Multiple sequence alignment of the peptide sequences enriched in mosquito-inoculated and subcutaneously inoculated NHP were done against the DENV1 reference sequence (Fig. 3B). At 7 days p.i. at 38/46 of the positions in the targeted region (E 339-384 ), more than 90% of the aligned peptides covering each position were identical to the reference sequence regardless of the route of inoculation. At both 35 and 90 days p.i., 42/46 positions were > 90% conserved in the mosquito group and 41/46 in the subcutaneous group. By 365 days p.i., the proportions of > 90% conserved residues were 43/46 and 42/46 in the mosquito and subcutaneously-inoculated groups. In both groups, only 2 positions, Alanine 368 and Valine 379 demonstrated < 70% sequence conservation relative to the reference sequence.
MHC Class I epitopes map to the region of differential sero-reactivity in EDIII. The immune epitope database (IEDB) was queried for validated MHC Class I, MHC Class II and B-cell linear epitopes contained in the differentially sero-reactive region in EDIII. A total of 11 MHC Class I epitopes were found to map to this region, but no MHC Class II or B-cell linear epitopes (Fig. 3C).
EDIII targeting did not correlate with late neutralization titer. The ELISA data from 365 days p.i.
did not correlate with late (13 month) neutralization titers for either the mosquito-inoculated or the subcutaneously inoculated group, with Pearson correlation coefficients of − 0.34 and 0.34 for each group, respectively. We next determined whether the phage display results for the targeted region of EDIII correlated with the ELISA data, calculating the Pearson correlation between the sum of the median enrichments spanning residues 339-384 and the ELISA OD 450 ratios. The correlation coefficients were 0.84 and 0.52 for the mosquito and subcutaneously inoculated NHP. Because the ELISA data for the mosquito group correlated well with phage display results, we next examined correlations between the summed enrichments of all 45 amino acid windows across the DENV1 proteome for the 365 days p.i. profiles and the 13 month neutralization titers (Fig. 4). The maximum correlation we observed in the mosquito-infected group was 0.91 for a 45 amino-acid window centered at amino acid residue 2618 in the NS5 protein, although several other amino acid windows also achieved a correlation of ≥ 0.9, predominantly in the NS3 and NS5 proteins. At this same timepoint, the maximum correlation observed for any 45 amino acid window and the 13 month neutralization titers in the subcutaneously infected group was 0.62.

Discussion
The primary goal of this work was the determination of a high-resolution map of sero-reactivity after DENV infection in non-human primates using a comprehensive set of tiled linear 62 amino acid peptides displayed on phage. Secondarily, we sought to determine whether the route of inoculation resulted in differences in this map. While serial heterotypic DENV infection by mosquito bite induces long-term immunity, it is unknown whether this route of inoculation affects the humoral immune response. In this study, serological profiling of primary DENV1 infections by programmable phage display demonstrated statistically significant differences in immune profiles in NHP inoculated by mosquito bite versus subcutaneous injection. Sero-reactivity to the DENV E, NS1, NS3 and NS5 proteins was observed at 35 and 90 days post-primary infection and did not vary by route. In contrast, we observed minimal sero-reactivity at 7 days post-primary infection, suggesting that the signals measured by PhIP-Seq in this study reflect IgG and not IgM responses. At one year post-infection, however, reactivity to the E and NS3 proteins demonstrated statistically significant increases in mosquito-inoculated animals versus subcutaneously-inoculated animals. The sustained differential reactivity to the E protein centered on residues 339-384 in domain III (ED III), a result that orthogonally validated by ELISA. The targeted region showed a moderate degree of sequence conservation across DENV serotypes. There was no difference in antibody breadth overall, as measured by sequence diversity, by route of inoculation. Although none of the reported validated MHC Class II or B-cell linear epitopes mapped to the targeted EDIII region, multiple validated MHC Class I epitopes did. ELISA results for the targeted ED III region did not correlate with 13-month neutralization titers, but further analysis of phage data suggested several regions in the NS3 and NS5 proteins that might do so in animals inoculated by mosquito. There are multiple potential causes for the late differences in sero-reactivity observed after mosquito and subcutaneous infection with DENV that we observed. Mosquito saliva, introduced during biting, increases DENV infectivity and enhances viral replication and could have resulted in a higher more immunogenic inoculum [15][16][17] . Arguing against this possibility, however, is that there was in fact delayed onset of viremia in   www.nature.com/scientificreports/ Given the large number of uncharacterized proteins in mosquito saliva, it is plausible that one or more of these could function as an adjuvant [22][23][24] . The virus stock used to infect NHP in the parent study was produced in Vero cells and is another potential explanation of our findings. For mosquito infections, the same virus stock was injected into the thoracic cavity of mosquitoes which were then used to infect NHP 14 days later. During this time, multiple viral replication cycles would have occurred in the mosquitoes. Post-translational modifications (PTM) of the DENV E protein, which contains two well-documented glycosylation sites, play an important role in viral attachment and entry 25,26 . PTM machinery in mammalian and insect cells is different, with mammalian cells capable of producing more complex structures 27 . These PTM differences have not been reported to affect viral infectivity, but could affect the innate immune response and thus shape the adaptive immune response to DENV infection 28 . While mosquito-derived virus would only be present for early rounds of infection, we nevertheless cannot rule out the possibility that the differences in sero-reactivity that we observed were attributable to the fact that the virus that was introduced in mosquito-inoculated NHP was different than that of subcutaneously inoculated NHP.
The DENV E protein is the primary target of the neutralizing antibodies (nAb) which are crucial for immunity to DENV, with ED III mediating viral attachment through its lateral ridge 29 . Anti-EDIII antibodies, however, constitute only a small proportion (< 10%) of serum DENV nAb in humans and are not required for DENV neutralization 6,30,31 . It is thus unsurprising that we observed no correlation between ED III ELISA results and 13-month neutralization titers, although our phage display results do suggest sero-reactivity to parts of NS3 and NS5 could be useful as candidate biomarkers. While the functional significance of the differential sero-reactivity to EDIII we observe remains to be determined, it is notable that the region we identified does contain multiple reported Class I epitopes, suggesting the potential immunological relevance of the findings we report here. It is surprising that we identified no MHC Class II epitopes, but this may be due to biased database representation, with roughly three times more Class I epitopes (1195) as Class II epitopes (428).
Using linear peptides displayed on bacteriophage to understand humoral immune responses has some limitations. DENV nAbs are thought to be against conformational epitopes defined by non-contiguous residues in the E protein and cannot be represented in 62 amino acid peptides 31,32 . Phage displayed peptides also lack eukaryotic PTM that could be important for antibody recognition. Nevertheless, using this approach to characterize the full spectrum of DENV sero-reactivity in the context of natural infection can provide valuable insights about DENV immunity. Finally, while NHP are the closest model to approximate DENV replication in humans, the findings we report here represent a single group of animals and need to be verified to determine if they generalize.
There are three DENV vaccine candidates currently that have either completed or are in the midst of phase III clinical trials and serologic diagnostics for DENV remain challenging. The results reported here show that in an NHP infection model, infection by mosquito can result in very different humoral immune profiles than infection by subcutaneous injection and that these differences can be long-lasting. While the precise causes for these differences and their functional consequences remain to be determined, the findings reported here suggest that natural immunity to DENV infection is influenced by vector delivery and should be investigated further.

Dengue virus cell culture and non-human primate infections. Culture of the DENV1 and DENV2
strains and NHP infections have previously been described 14 . Phage display library design and construction. All protein sequences present in the National Center for Biotechnology Information (NCBI) GenBank database as of November 4, 2017 were downloaded for 74 flaviviruses known to infect humans. For each virus a multiple sequence alignment was performed and the aligned sets of sequences divided into sets of overlapping 62 amino acid peptides. Consecutive peptides overlapped by 54 amino acids. The full set of peptides was collapsed on 98% sequence similarity, resulting in a final set of 91,562 flavivirus-specific peptides. Nucleic acid sequences coding for each peptide were designed using randomized sets of codons such that overlapping regions of consecutive peptides did not have the same nucleic acid sequence. Common 5′ (GCA GGA GTA GCT GGT GTT GTG, coding for AGVAGVV) and 3′ (TGA TAA GCA TAT GCC ATG GCCTC) linker sequences were appended to each oligonucleotide coding sequence, all sequenced outputted to a FASTA format file and sent to Agilent, Inc. for synthesis. All bioinformatic analysis was done using R and R studio.
Lyophilized oligonucleotides were resuspended in 10 mM Tris-HCl-1 mM EDTA, pH 8.0 at 2 nM. The oligonucleotide stock was diluted to 0.2 nM and used as the template for PCR. The forward primer was TAG TTA AGC GGA ATT CAG CAG GAG TAG CTG GTG TTG TG and the reverse primer was ATC CTG AGC TAA GCT TGA GGC CAT GGC ATA TGC TTA TCA. Cycling conditions used were: 98 °C × 30 s followed by 20 cycles of 98 °C × 5 s, 70 °C × 15 s and 72 °C × 10 s and terminating with 72 °C × 2 min. A total of 4 independent PCR reactions (50ul reactions) were pooled together and cleaned using Ampure XP beads at 1×. Bead clean-up followed the manufacturer's protocol. Bead-cleaned PCR products were quantified using Qubit reagent. One microgram (1ug) of PCR was digested with 1ul of EcoR I-HF (New England Biolabs) and 1ul Hind III-HF (New England Biolabs) in 1× Cutsmart buffer and a final reaction volume of 50ul for 1H at 37 °C. Restriction digests were heat inactivated for 10 min at 65 °C, cleaned with 1X Ampure XP beads, resuspended in a final volume of 40ul nuclease-free water and quantified by Qubit. Cut and uncut template was visualized using an Agilent Bioanalyzer high sensitivity DNA chip to confirm 100% digestion. Digested product (0.06 pmol) was ligated into pre-EcoR I/Hind III digested T7 phage vector arms (EMD Millipore, Inc) with 1ul T4 DNA ligase (NEB), 0.5 ul T4 DNA ligase buffer, 1ul of pre-digested T7 vector arms (0.02 pmol) and nuclease-free water to a final reaction volume of 5ul. Four (4) ligations were set up in parallel, incubated at 16 °C × 12H in a thermal cycler and heat inactivated for 10 min at 65 °C. The entire ligation reaction (5ul) was mixed with 25ul of T7 phage packaging extract (EMD Millipore, Inc) Non-human primate infections and samples. Infection of rhesus macaques and collection of samples has previously been described. Serum samples were shipped from WRAIR to San Francisco on dry ice, thawed and randomly aliquoted into plates so that samples from the same individual were not contiguous in the plate. All serum samples were diluted 1:1 with 2× sample storage buffer (40% glycerol, 40 mM HEPES pH 7.3, 0.04% NaN 3 , PBS) and stored in a 4 °C cold-room after thawing. Samples from mosquito-infected and subcutaneously infected subjects were distributed between two plates.
Phage immunoprecipitation sequencing. All immunoprecipitation experiments were conducted in 2 ml deep-well polypropylene plates (Genesee Scientific, Inc). Plates were blocked with 3% BSA in TBST overnight with overhead rotation on a Rotator Genie prior to use for immunoprecipitations. After removal of blocking buffer, 1 ml of phage library (2 × 10 10 pfu) was decanted with a multichannel pipet by hand. Serum samples were added directly to each well using an Integra ViaFlow 96 multichannel pipettor, 1ul per well. After addition of serum, plates were sealed with a rubber sealing mat (Genesee Scientific, Inc) and rotated overhead overnight on a Rotator Genie. Plates were removed from the rotator the following morning and spun at 1500RPM in an Eppendorf 5810R centrifuge for 5 min. The sealing mat was carefully removed, and 40ul of a mix of Protein A (20ul) and Protein G (20ul) Dynabeads (Thermo-Fisher, Inc) resuspended in Tris/NP-40 was added to each well by multichannel pipet by hand. Protein A and Protein G beads were washed 3 × with overhead rotation × 5 min in TNP-40 prior to addition to phage-library/serum incubations. Immunoprecipitations were carried out by overhead rotation for 1H at 4 °C, after which plates were spun at 400 RCF × 1 min at 4 °C in the 5810R centrifuge. The sealing mat was removed, the plate placed on a magnet compatible with 96-well plates and the beads allowed to collect on the magnet and the liquid aspirated. 900ul of RIPA buffer was added to each well, the sealing mat replaced and plates rotated overhead at 4 °C for 5 min. After rotation, plates were spun at 400RCF in the 5810R centrifuge, the sealing mat removed, beads collected on the magnet, wash buffer aspirated and a fresh 900ul of RIPA added. This process was repeated for a total of 5 washes. After the fifth wash, 50ul of phage elution buffer was added to the beads, gently agitated by hand and incubated for 3 min on ice. After incubation, beads were collected and inoculated into 1 ml of OD 600 0.5 E coli strain BLT-5403 in a fresh deep-well plate, sealed with a gas-permeable plate sealer and incubated with agitation at 37 °C at 750 rpm using an Infors shaking incubator until all wells showed 100% lysis (1.5-2H). After lysis, 100ul of 5 M NaCl was added to each well, mixed and spun at 3220 RPM for 10 min. 500ul of spun lysate was added to each well of a blocked deep 96-well plate and an additional 600ul of LB-Carbenicillin added to each well, followed by 1ul of each serum sample to each well. Plates were sealed with a fresh sealing mat and incubated overnight with rotation as already described. The next day, immunoprecipitations and washes were carried out as already described, with the exception that the lysate was not used for an additional round of selection but for sequencing library preparation.
Lysates were heated to 70 °C for 15 min to liberate phage genomes. For each immunoprecipitation, 2ul of lysate was used as the template for a PCR reaction to add Illumina P5 and P7 sequences. Forward and reverse primers for the PCR reaction were AAT GAT ACG GCG ACC ACC GAG ATC TACACNNNNNNNNATG GGC CAC GGT GGT CTT CGCCC and CAA GCA GAA GAC GGC ATA CGA GAT NNNNNNNNGGG TTA ACT AGT TAC TCG AGT GCG GCCG, with 96-different barcodes. PCR reactions were done in 25ul volumes with 2 ul lysate, 1 ul of each primer at 10uM, 0.5 ul of 10 mM dNTP, 5 ul of 5× Phusion buffer, 0.25 ul pf Phusion Taq DNA polymerase and nuclease free water to the final volume. PCR conditions were identical to those already listed for library cloning. After PCR, 5 ul of each reaction was pooled, cleaned up as already described with Ampure XP beads and quantified using the KAPA Biosystems Illumina quantification kit. All samples were subsequently submitted to the Chan-Zuckerberg Biohub sequencing facility and sequenced on a NextSeq 550 high output (400 million reads) cartridge using custom read 1 and indexing primers. A single-ended 150 base pair sequencing strategy (sequencing 150/186 base pairs of inserts) was used, with a goal read-depth of 1.5-2 million reads/sample. Bioinformatic analysis of PhIP-Seq results. Demultiplexed sequencing data was aligned to the library oligonucleotide sequences using the Bowtie2 aligner. All alignments that were not full-length, perfect alignments were discarded. Perfect alignments were translated and only those that were perfect matches to the input library retained. This set of sequences was tabulated for each of the immunoprecipitation reactions, resulting in a read-count per peptide. Peptide read-counts were converted to reads per 100,000 reads (rp100k) for each sample by dividing by the total number of reads in a sample and multiplying by 100,000.
All samples were organized by animal in temporal order, starting with baseline-pre-infection samples. To determine which hits were significant, we calculated enrichment scores for all identified peptides relative to a reference sample for each animal. For samples that were 7, 35, 90 and 365 days post-primary infection, the For any peptide not present in the baseline sample, the median rp100k from the reference sample was used to calculate enrichments. The fit of several probability distributions with enrichment data was examined, including the normal, negative binomial and Poisson. The normal distribution was found to have the best fit, particularly after log transformation of enrichments. In brief, a set of enrichments for a sample was log-transformed and a z-score calculated for each peptide. Peptides with a z-score > 3 were considered as significantly enriched compared to the reference sample.
All peptides in the library were aligned against reference genomes for all four DENV serotypes using blastp. These data were then used to classify peptides as either DENV1 type-specific or cross-reactive. DENV1 typespecific peptides were those that were ≥ 90% homologous across their full length to DENV1 alone, while crossreactive peptides were those that were at least 50% homologous to at least 2 DENV strains. Using these data, the proportion of significantly enriched peptides for each timepoint for each animal by each route of infection was determined.
For generating peptide coverage maps, a different strategy was used. Some regions of the flavivirus genome are highly conserved across multiple viruses (e.g. the fusion loop sequence of the envelope protein) and linear peptide sequences recognized by antibodies, such as those represented in the phage display library used here are short. To account for both of these factors, all peptides in the library were decomposed into a series of overlapping 10mers, with each 10mer aligned to the DENV1 reference genome using blastp and only those that were ≥ 90% homologous to the DENV1 reference genome retained for downstream analysis. For each sample, each significant peptide was mapped to its set of DENV1 homologous 10mers. Using these alignment data and enrichment scores, a cumulative enrichment score was calculated at each position of the DENV1 proteome. As before, these data were computed for each animal, at each timepoint by each route of infection. The median cumulative enrichment art each position was computed and used to generate coverage plots for each timepoint by each route of infection.

Statistical analysis of phage immunoprecipitation sequencing results.
Area under the curve (AUC) was calculated using the median cumulative enrichment data generated as described above. These were calculated for all 10 amino-acid windows across the DENV1 proteome but also for specific regions of interest in the E, NS1, NS3 and NS5 proteins. Permutation testing was used to test for the statistical significance of the differences of AUCs in the mosquito-inoculated and subcutaneously inoculated groups. Briefly, all 20 animals were randomly assigned to the mosquito-inoculated or the subcutaneously inoculated group AUCs calculated and the difference between mosquito and subcutaneously groups calculated. This process was repeated for a total of 1000 iterations, generating the null distribution for AUC msq − AUC sc across the DENV1 proteome. The measurements at each distribution were used to model the cumulative distribution function of AUC msq − AUC sc and to estimate P(X > x) at that position.

ED III ELISAs.
Envelope protein residues 301-394 were PCR amplified from a gene block of the Entire DENV1 E protein reference sequence ordered from IDT, Inc. An internal Hind III site was removed by PCR and a second PCR performed to add a 5′ EcoR I site, a 3′ FLAG tag (DYKDDDDK) followed by a 3′ Hind III site. The PCR product was digested in a 50ul reaction as previously described and full digestion confirmed by Agilent Bioanalyzer. Digested product was cloned into an expression vector (OriGene pCMV6), sequenced confirmed by Sanger sequencing (Quintara Bio, Inc.) and CaPO 4 transfected into 50% confluent HEK293 cells. Transfected cells were lysed with RIPA buffer, to which a protease inhibitor cocktail pill was added (Roche, Inc). Lysates were spun at 13.2 K RPM in an Eppendorf microcentrifuge at 4 °C for 15 min, supernatants aspirated, aliquoted in 200-400ul aliquots and stored at − 80C until use.
Clear polystyrene plates for ELISAs (R & D Biosystems, Inc) were adsorbed with anti-FLAG antibody (Cell Signalling, Inc.) diluted 1:2500, 50ul per well overnight at 4 °C in a sealed plate. Plates were washed 3 times with PBS + 0.05% Tween-20 (PBST), 200ul of 2% BSA + PBST added to each well and plates incubated at room temperature for 2H. Blocking buffer was removed, sample diluted 1:100 in blocking buffer added to each well, the plate sealed and incubated overnight at 4 °C. After sample incubation, all wells were washed 4 times with 300ul PBST per well by multichannel pipet, 50ul of anti-human IgG-HRP conjugated detection antibody (Sigma Chemical, #) added, the plates sealed and incubated at room temperature in the dark for 1H. Plates were washed 5× with 300ul PBST per well as already described and 50ul of TMP Substrate/Peroxidase solution added. Reactions were stopped using 2 N H 2 SO 4 and OD450 measured on a plate reader. For each sample series, the ratio of the OD 450 to pre-infection OD 450 was calculated. Median ratios were compared using a one-sided Mann-Whitney U test to determine the statistical significance of the observed differences.

Analysis of antibody binding specificities.
For each animal at each timepoint by each route of infection, the significantly enriched peptides that spanned residues 347-377 of the E protein were identified. Multiple sequence alignments were performed using the package DECIPHER and the proportion of aligned sequences at that position whose amino acid residue matched that of the reference sequence was calculated.
Correlation analysis with neutralization data. We used the R function cor to calculate Pearson correlation coefficients between log transformed ELISA data and 13 month geometric mean neutralization titers (GMT) against DENV1 that had been determined for each of the 20 animals. For the analysis of correlation between phage data and the GMT data, we calculated the sum of all enrichments for all possible 45 amino-acid windows for each animal for both the mosquito and subcutaneously inoculated groups. We then correlated these Scientific Reports | (2021) 11:10823 | https://doi.org/10.1038/s41598-021-90318-z www.nature.com/scientificreports/ vectors of summed enrichments with GMTs and plotted the resulting correlation coefficients for each position. All of these analyses were conducted using R.
Ethics statement. The work described here was approved by the University of California, San Francisco Institutional Review Board and classified as 'Not Human Subjects Research' . All methods were carried out in accordance with relevant guidelines and regulations.