*Article* **Correlates with Vaccine Protective Capacity and COVID-19 Disease Symptoms Identified by Serum Proteomics in Vaccinated Individuals**

**Margarita Villar 1,2 , José Miguel Urra 3,4, Sara Artigas-Jerónimo <sup>1</sup> , Lorena Mazuecos <sup>1</sup> , Marinela Contreras <sup>1</sup> , Rita Vaz-Rodrigues <sup>1</sup> , Francisco J. Rodríguez-del-Río 5 , Christian Gortázar <sup>1</sup> and José de la Fuente 1,6,\***


**Abstract:** In the last two years, the coronavirus disease 19 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been a scientific and social challenge worldwide. Vaccines have been the most effective intervention for reducing virus transmission and disease severity. However, genetic virus variants are still circulating among vaccinated individuals with different disease symptomatology. Understanding the protective- or disease-associated mechanisms in vaccinated individuals is relevant to advances in vaccine development and implementation. To address this objective, serum-protein profiles were characterized by quantitative proteomics and data-analysis algorithms in four cohorts of uninfected and SARS-CoV-2-infected vaccinated individuals with asymptomatic, non-severe, and severe disease symptomatology. The results show that immunoglobulins were the most overrepresented proteins in infected cohorts when compared to PCR-negative individuals. The immunoglobulin profile varied between different infected cohorts and correlated with protective- or disease-associated capacity. Overrepresented immunoglobulins in PCR-positive individuals correlated with protective response against SARS-CoV-2, other viruses, and thrombosis in asymptomatic cases. In non-severe cases, correlates of protection against SARS-CoV-2 and HBV together with risk of myasthenia gravis and allergy and autoantibodies were observed. Patients with severe symptoms presented risk for allergy, chronic idiopathic thrombocytopenic purpura, and autoantibodies. The analysis of underrepresented immunoglobulins in PCR-positive compared to PCR-negative individuals identified vaccine-induced protective epitopes in various coronavirus proteins, including the spike receptor-binding domain RBD. Non-immunoglobulin proteins were associated with COVID-19 symptoms and biological processes. These results evidence host-associated differences in response to vaccination and the possibility of improving vaccine efficacy against SARS-CoV-2.

**Keywords:** COVID-19; proteomic; vaccine; immunology; biomarker

### **1. Introduction**

Millions of deaths have been reported worldwide associated with coronavirus disease 19 (COVID-19), a pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (https://www.google.com/search?client=safari&rls=en&q=COVID-19 +worldwide+cases&ie=UTF-8&oe=UTF-8#colocmid=/m/02j71&coasync=0).

**Citation:** Villar, M.; Urra, J.M.; Artigas-Jerónimo, S.; Mazuecos, L.; Contreras, M.; Vaz-Rodrigues, R.; Rodríguez-del-Río, F.J.; Gortázar, C.; de la Fuente, J. Correlates with Vaccine Protective Capacity and COVID-19 Disease Symptoms Identified by Serum Proteomics in Vaccinated Individuals. *Molecules* **2022**, *27*, 5933. https://doi.org/ 10.3390/molecules27185933

Academic Editors: Giovanni N. Roviello and Caterina Vicidomini

Received: 1 August 2022 Accepted: 7 September 2022 Published: 13 September 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

Furthermore, the global number of deaths caused by COVID-19 may be up to four times this figure [1]. Vaccines using different platforms have been developed as the most safe and effective intervention for reducing SARS-CoV-2-virus transmission and disease severity [2,3]. However, genetic variants of the coronavirus are still circulating among vaccinated individuals with different disease symptomatology [4,5]. Understanding the protective- or disease host-associated mechanisms in vaccinated individuals is relevant to advances in vaccine development and implementation [6].

To address this challenge, understanding the protective- or disease-associated mechanisms in vaccinated individuals is relevant to advances in vaccine development and implementation. Previous proteomics studies, e.g., [7,8], have addressed response to infection and vaccination, but our study addressed for the first time the immune response to COVID-19 vaccination in uninfected and SARS-CoV-2-infected individuals with asymptomatic, non-severe, and severe disease symptomatology. In this study, serum-protein profiles were characterized by previously validated quantitative proteomics [7] in four cohorts of uninfected and SARS-CoV-2-infected vaccinated individuals with asymptomatic, non-severe, and severe disease symptomatology. The results evidence host-associated differences in response to vaccination and the possibility of advances in vaccine development and implementation against SARS-CoV-2.

#### **2. Materials and Methods**

#### *2.1. Study Design with Serum Samples from Different Cohorts*

The study design is described in Figure 1. A retrospective case-control study was conducted in patients suffering from COVID-19 and healthy controls sampled at the University General Hospital of Ciudad Real (HGUCR), Spain [8]. Individuals were confirmed as SARS-CoV-2-infected by reverse transcriptase-polymerase chain reaction (RT-PCR) and sampled between November and December 2021 (Table 1). In this study with individuals vaccinated against COVID-19, vaccine administration, clinical symptoms, and laboratory determinations associated with COVID-19 were obtained from patients' medical records to create cohorts of PCR– and PCR+ asymptomatic, non-severe, and severe individuals (Table 1). Patient symptoms can be found in Table 1. Blood samples were drawn in a vacutainer tube without anticoagulant. The tube remained at rest for 15–30 min at room temperature (RT) for clotting. Subsequently, the tube was centrifuged at 1500× *g* for 10 min at RT to remove the clot and obtain serum. Serum samples were heat-inactivated for 30 min at 56 ◦C and conserved at −20 ◦C until used for analysis. The use of samples and individuals' data was approved by the Ethical and Scientific Committees (University Hospital of Ciudad Real C-352 and SESCAM C-73).

**Table 1.** Data on serum samples included in the analysis.



```
Table 1. Cont.
```
(\*) Days between first vaccine dose and serum sampling are shown in parentheses.

1

**Figure 1.** Experimental design and rationale. The experimental design used in our study was based on sera collected from vaccinated individuals (mostly with Pfizer and Moderna-BioNTech) and SARS-CoV-2 PCR-negative (PCR–) or infected PCR-positive (PCR+) and with asymptomatic, non-severe, and severe COVID-19 symptomatology (Table 1). Two PCR+ cases with severe symptoms and in the hospital intensive-care unit (ICU) were included only as reference for selected proteins. Sera were collected between 156 and 298 days after first vaccine-dose administration and subjected to SWATH-MS quantitative proteomics to characterize serum-protein profiles in different cohorts. The proteomics results were then translated into the identification of correlates with protective- or disease-associated capacity and vaccine-induced protective epitopes.

#### *2.2. Serum Proteomics*

The methodology and algorithms for serum proteomics were as previously described [8]. Serum samples from PCR– controls and PCR+ COVID-19 asymptomatic, non-severe, and severe individuals (n = 5 each) were individually analyzed. Two PCR+ cases with severe symptoms and in the hospital ICU were included only as reference for the selected proteins. Protein concentration in samples was determined using the BCA Protein Assay with BSA (Sigma-Aldrich) as standard. Protein-serum samples (100 µg per sample) were trypsin digested using the FASP Protein Digestion Kit (Expedeon Ltd., UK) and sequencing-grade trypsin (Promega, Madison, WI, USA) following the manufacturer's recommendations. The resulting tryptic peptides were desalted onto OMIX Pipette tips C18 (Agilent Technologies, Santa Clara, CA, USA), dried down, and stored at −20 ◦C until massspectrometry analysis. The desalted protein digests were resuspended in 2% acetonitrile, 5% acetic acid in water, and analyzed by reverse-phase liquid chromatography coupled

with mass spectrometry (RP-LC-MS/MS) using an EkspertTM nanoLC 415 system coupled online with a 6600 TripleTOF® mass spectrometer (AB SCIEX; Framingham, MA, USA) through Information-Dependent Acquisition (IDA) followed by Sequential Windowed data-independent Acquisition of the Total High-resolution Mass Spectra (SWATH-MS). The peptides were concentrated in a 0.1 × 20 mm C18 RP precolumn (Thermo Fisher Scientific, Waltham, MA, USA) with a flow rate of 2 µL/min for 10 min in solvent A. Then, peptides were separated in a 0.075 × 250 mm C18 RP column (New Objective, Woburn, MA, USA) with a flow rate of 300 nl/min. Elution was carried out in a 120 min gradient from 5% B to 30% B followed by a 15 min gradient from 30% B to 60% B (Solvent A: 0.1% formic acid in water, solvent B: 0.1% formic acid in acetonitrile) and directly injected into the mass spectrometer for analysis. For IDA experiments, the mass spectrometer was set to scanning full spectra from 350 *m*/*z* to 1400 *m*/*z* (250 ms accumulation time), followed by up to 50 MS/MS scans (100–1500 *m*/*z*). Candidate ions with a charge state between +2 and +5 and counts per second above a minimum threshold of 100 were isolated for fragmentation. One MS/MS spectrum was collected for 100 ms before adding those precursor ions to the exclusion list for 15 s (mass spectrometer operated by Analyst® TF 1.6, ABSciex®). Dynamic background subtraction was turned off. Data were acquired in high-sensitivity mode with rolling collision energy on and a collision energy spread of 5. An equal amount of the five samples for each experimental group joined together as a representative mixed sample of each of the 4 experimental groups, which were used for the generation of the reference spectral-ion library as part of SWATH-MS analysis. A total amount of 4 µg protein digests for each mixed sample was injected. For SWATH quantitative analysis, 50 independent samples (2 technical replicates from each of the 5 biological replicates for each of the 4 experimental groups and 5 technical replicates from each of the 2 biological replicates in the case of ICU samples) (6 µg each) were subjected to the cyclic data-independent acquisition (DIA) of mass spectra using the SWATH variable-window calculator (V 1.0, AB SCIEX) and the SWATH acquisition-method editor (AB SCIEX), similar to established methods [8]. A set of 50 overlapping windows was constructed (containing 1 *m*/*z* for the window overlap), covering the precursor mass range of 400–1250 *m*/*z*. For these experiments, a 50 ms survey scan (350–1400 *m*/*z*) was acquired at the beginning of each cycle, and SWATH-MS/MS spectra were collected from 100–1500 *m*/*z* for 70 ms in high-sensitivity mode, resulting in a cycle time of 3.6 s. Collision energy for each window was determined according to the calculation for a charge +2 ion centered upon the window with a collision-energy spread of 15. To create a spectral library of all the detectable peptides in the samples, the IDA MS raw files were combined and subjected to database searches in unison using Protein-Pilot software v. 5.0.1 (AB SCIEX) with the Paragon algorithm. Spectra identification was performed by searching against the Uniprot human-proteome database (79,038 entries in January 2022) with the following parameters: iodoacetamide cysteine alkylation, trypsin digestion, identification focus on biological modification, and thorough ID as search effort. The detected protein threshold was set at 0.05. To assess the quality of identifications, an independent False Discovery Rate (FDR) analysis with the target-decoy approach provided by Protein PilotTM was performed. Positive identifications were considered when identified proteins reached a 1% global FDR. The mass-spectrometry proteomics data were deposited in the Proteome Xchange Consortium via the PRIDE partner repository with the dataset identifier PXD031969 and 10.6019/PXD031969.

#### *2.3. Quality Control of Proteomics Data*

The quality of the proteomics data was controlled at multiple levels. First, a rat-ileum digest was used for the evaluation of instrument performance. Buffer A samples were run as blanks every two injections to prevent carryover. Two technical replicates were injected for each sample. For validation of serum-proteomics data, protein representation for previously identified selected biomarkers for COVID-19 and proteomics studies were used to show correlation with disease severity. An enrichment analysis was conducted using the Coronascape COVID database (https://metascape.org/COVID; [8]) to identify proteins found in our study as differentially represented in response to COVID-19 and reported in previous COVID-19 omics datasets.

#### *2.4. Data Analysis*

For SWATH processing, up to 10 peptides with 7 transitions per protein were automatically selected by the SWATH Acquisition MicroApp 2.0 in the PeakView 2.2 software with the following parameters: 15 ppm ion library tolerance, 5 min XIC extraction window, 0.01 Da XIC width, and considering only peptides with at least 99% confidence and excluding those that were shared or contained modifications. However, to ensure reliable quantitation, only proteins with 3 or more peptides available for quantitation were selected for XIC peak-area extraction and exported for analysis in the MarkerView 1.3 software (AB SCIEX). Global normalization according to the total area sums of all detected proteins in the samples was conducted (Supplementary Materials Data File S1). The Student's *t*-test (*p* < 0.05) was used to perform two-sample comparisons between the averaged area sums of all the transitions derived for each protein across the 10 replicate runs for each group under comparison to identify proteins that were significantly represented between groups (Supplementary Materials Data File S1). Protein representation was compared between groups by One-way ANOVA test followed by post-hoc Bonferroni and Holm multiple comparisons tests (*p* < 0.05; https://astatsa.com/OneWay\_Anova\_with\_TukeyHSD/), and relative intensity was compared between PCR- and PCR+ cohorts by Welch's unpaired *t*-test (*p* < 0.05; https://www.graphpad.com/quickcalcs/ttest1/?Format=C) [8]. Data were separately analyzed for overrepresented and underrepresented proteins using the Metascape gene annotation and analysis resource (https://metascape.org/gp/ index.html#/main/step1). The analytical algorithm developed using Protein BLAST sequence alignment against non-redundant protein database (nr) using compositional matrix adjustment (https://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastp&PAGE\_TYPE= BlastSearch&LINK\_LOC=blasthome) and Paratome (http://www.ofranlab.org) was used for the identification of antigen-binding regions and vaccine-induced antibody-protective epitopes and correlates of identified proteins with protective- or disease-associated capacity (Supplementary Materials Data Files S2 and S4).

#### *2.5. Antibody Neutralization Test*

Antibody titers specific for the neutralization of SARS-CoV-2 virus were determined with a cPass SARS-CoV-2 neutralization antibody-detection kit (Genscript, Piscataway, NJ, USA) following the manufacturer0 s instructions. Briefly, 100 µL of the positive and negative controls and serum samples at 1:100 dilution and previously incubated with HRP conjugated RBD during 30 min at 37 ◦C were added to the 96-microwell plate coated with RBD-SARS-CoV-2 protein and incubated for 15 min at 37 ◦C. After washing four times with 260 µL/well of wash buffer, 100 µL/well of chromogen-substrate solution were added and incubated for 15 min at RT. Finally, the colorimetric reaction was stopped with 50 µL/well of stop solution and the absorbance was measured in a spectrophotometer (Thermo Fisher Scientific) at O.D. 450 nm. Results were evaluated by calculating the ratio between the O.D. of the sample and the O.D. of the calibrator using the following formula: Inhibition = (1 − O.D. value of sample/O.D. value of negative control) × 100.

#### *2.6. Antibody Levels against HBV and Zika Virus*

Individual sera from all cohorts included in the study were characterized for antibody levels against HBV and Zika virus using pathogen-specific ELISA tests. The ELISA kits are designed for the detection of antibodies to hepatitis B-virus surface antigen: Hepatitis B surface antigen Ab ELISA kit (AB-KA0287; Biogen Científica, Madrid, Spain) or anti-Zika Virus Non-structural Protein, Anti-Zika virus IgG ELISA kit (ab221844; Abcam, Cambridge, UK) in human serum. Antibody levels were compared between PCR– and PCR+ cohorts by Chi-squared test (*p* < 0.05). A Spearman's Rho (*rs*) correlation analysis was conducted

between virus cross-reactive Ig levels and proteomic-protein relative intensity (*p* < 0.05; https://www.socscistatistics.com/tests/spearman/default2.aspx).

#### *2.7. Human Autoantibody General Survey Microarray*

Reactive autoantibodies were characterized with individual sera from all cohorts using the Human Autoantigens General Survey Antigen Microarray (GeneCopoeia's OmicsArray; Rockville, MD, USA), a protein microarray enabling powerful detection of autoantibodies associated with many diseases, including rheumatoid arthritis, muscular dystrophy, systemic lupus erythematosus, and type-1 diabetes. The array carries 120 superior-quality purified proteins spotted onto nitrocellulose filters, which are adhered to glass slides. Antigens known to be associated with specific autoimmune diseases are chosen based on a thorough review of peer-reviewed publications. Data interpretation: NSI, net signal intensity, averaged fluorescent-signal intensity for each antigen subtracted by local background and negative control signal; NSI-nor, NSI normalized to internal Ig controls; SBR, averaged signal-to-background ratio of each antigen; SBR-nor, SBR normalized to internal Ig controls; SNR, signal-to-noise ratio that represents the significance of the signal above the background (SNR ≥ 3 means the signal is significantly higher than the background). Additional information is in Supplementary Materials Data File S3.

#### **3. Results and Discussion**

#### *3.1. Characterization of Immunoglobulin Protein Profiles and Correlation with Protective- or Disease-Associated Capacity*

The experimental design used in our study was based on sera collected from SARS-CoV-2 PCR-negative (PCR–) or infected PCR-positive (PCR+) vaccinated individuals (mostly with Pfizer and Moderna-BioNTech) and with asymptomatic, non-severe, and severe COVID-19 symptomatology (Figure 1, Table 1). Two PCR+ cases with severe symptoms and in the hospital intensive-care unit (ICU) were included only as reference for selected proteins. Sera were collected in November–December 2021 before the appearance of the SARS-CoV-2 Omicron variant in Spain, and thus most PCR+ cases were probably of the Delta virus variant (Table 1). As expected, 91% of the individuals' serum samples (20 out of 22 except for samples C5 and S1; Table 1) showed neutralization antibodies against SARS-CoV-2. Nevertheless, the characterization of antibody- and non-antibodymediated immune response to vaccination is important to understanding the response to COVID-19 vaccines. Proteomics analysis identified significantly dysregulated proteins in vaccinated and infected cohorts (asymptomatic, n = 134; non-severe, n = 117; severe, n = 230) when compared to vaccinated PCR– individuals (Figure 2A, Supplementary Materials Data File S1). More than 55% of the dysregulated proteins identified in infected cohorts when compared to PCR– individuals were immunoglobulins (Igs), mostly overrepresented in asymptomatic and severely infected cohorts (Figure 2A,B). This result is associated with response to vaccination, as supported by the finding using a similar serum-proteomics approach in healthy, unvaccinated individuals and with different COVID-19 symptomatology, in whom only 32% (60/189) of the dysregulated proteins were Igs [8].

An analytical workflow was developed to characterize selected Ig proteins identified as significantly dysregulated in vaccinated and infected cohorts when compared to vaccinated PCR– individuals (Supplementary Materials Data File S2). The use of the Paratome web server (http://www.ofranlab.org) allowed for the identification of antigen-binding regions in identified Ig light- or heavy-chain variable regions, including but not limited to complementarity-determining regions (CDRs) [9]. As has been discussed, heavy-chain complementarity-determining region 3 (HCDR3) is necessary, but insufficient for specific antibody binding [10].

**Figure 2.** Serum-protein profiles and correlation of Ig overrepresented in vaccinated infected cohorts with protective- or disease-associated capacity. (**A**) Number of identified overrepresented or underrepresented Ig and non-Ig proteins in infected asymptomatic, non-severe, and severe cohorts when compared to PCR– individuals. (**B**) Change in Ig-protein levels in PCR+ cohorts. PCR+/− Log fold-change relative intensity was compared between groups by One-way ANOVA test followed by post-hoc Bonferroni and Holm multiple comparisons tests (*p* < 0.05). (**C**) Heatmap of PCR+/− Log fold-change relative intensity (Z-scored original value) for Ig proteins overrepresented in infected cohorts. Correlates with protective- or disease-associated capacity are shown for Igs highly overrepresented in PCR+ individuals. **Figure 2.** Serum-protein profiles and correlation of Ig overrepresented in vaccinated infected cohorts with protective- or disease-associated capacity. (**A**) Number of identified overrepresented or underrepresented Ig and non-Ig proteins in infected asymptomatic, non-severe, and severe cohorts when compared to PCR– individuals. (**B**) Change in Ig-protein levels in PCR+ cohorts. PCR+/− Log fold-change relative intensity was compared between groups by One-way ANOVA test followed by post-hoc Bonferroni and Holm multiple comparisons tests (*p* < 0.05). (**C**) Heatmap of PCR+/− Log fold-change relative intensity (Z-scored original value) for Ig proteins overrepresented in infected cohorts. Correlates with protective- or disease-associated capacity are shown for Igs highly overrepresented in PCR+ individuals.

An analytical workflow was developed to characterize selected Ig proteins identified as significantly dysregulated in vaccinated and infected cohorts when compared to vaccinated PCR– individuals (Supplementary Materials Data File S2). The use of the Paratome web server (http://www.ofranlab.org) allowed for the identification of antigen-binding regions in identified Ig light- or heavy-chain variable regions, including but not limited to complementarity-determining regions (CDRs) [9]. As has been discussed, heavychain complementarity-determining region 3 (HCDR3) is necessary, but insufficient for specific antibody binding [10]. Focusing on Igs highly overrepresented in PCR+ individuals, the results showed differences between infected cohorts (Figure 2C, Supplementary Materials Data File S2). In vaccinated and infected asymptomatic cases, predictive models associated the Igs with protection against SARS-CoV-2, Zika virus, rotavirus, Hepatitis B virus (HBV), and thrombosis. However, in cases with COVID-19 symptoms, Igs associated with protection against SARS-CoV-2 and HBV were only identified in non-severe cases, whereas Igs Focusing on Igs highly overrepresented in PCR+ individuals, the results showed differences between infected cohorts (Figure 2C, Supplementary Materials Data File S2). In vaccinated and infected asymptomatic cases, predictive models associated the Igs with protection against SARS-CoV-2, Zika virus, rotavirus, Hepatitis B virus (HBV), and thrombosis. However, in cases with COVID-19 symptoms, Igs associated with protection against SARS-CoV-2 and HBV were only identified in non-severe cases, whereas Igs associated with autoantibodies and risk of allergic reactions and diseases such as myasthenia gravis and chronic idiopathic thrombocytopenic purpura (ITP) were identified only in non-severe and severe patients (Figure 2C, Supplementary Materials Data File S2). Additionally, autoantibodies were identified in all PCR+ cohorts. The analysis of mass-spectra relative intensity for selected overrepresented Ig proteins corroborated in individual samples the cohort-dependent results (Figure 3A). Then, an independent analysis of predicted biomarkers was used to validate these results for Zika virus, HBV, and a human-autoantibody general survey (Figure 3B–D).

man-autoantibody general survey (Figures 3B–D).

**Figure 3.** Validation of serum Ig proteins overrepresented in vaccinated infected individuals. (**A**) Changes in Ig-protein mass-spectra relative intensity for selected overrepresented Ig in individual samples in all PCR– and PCR+ cohorts. Relative intensity was compared between PCR– and PCR+ cohorts by Welch's unpaired t-test (\* *p* < 0.05, \*\* *p* < 0.005; n = 5 biological replicates). Samples from individuals in the ICU were only included as reference for selected proteins. (**B**) Independent analysis of predicted biomarker A0A5C2FZ03 for anti-Zika virus Ig levels. Antibody levels were compared between PCR– and PCR+ cohorts by Chi-squared test (*p* = 0.004 for asymptomatic cohort). A Spearman's Rho correlation analysis was conducted between virus cross-reactive Ig levels and proteomics-protein relative intensity (*p* < 0.05). (**C**) Independent analysis of predicted biomarkers A0A5C2FZ03 and A0A5C2GPZ0 for anti-HBV Ig levels in asymptomatic and non-severe cohorts, respectively. A Spearman's Rho correlation analysis was conducted between HBV cross-reactive Ig levels and proteomics-protein relative intensity (*p* < 0.05). (**D**) A human-autoantibody general survey identified proteins reactive to IgM and IgG autoantibodies with significant NSI-nor values in PCR– and PCR+ cohorts. Proteins were distributed based on the highest NSI-nor average IgM+IgG value. **Figure 3.** Validation of serum Ig proteins overrepresented in vaccinated infected individuals. (**A**) Changes in Ig-protein mass-spectra relative intensity for selected overrepresented Ig in individual samples in all PCR– and PCR+ cohorts. Relative intensity was compared between PCR– and PCR+ cohorts by Welch's unpaired t-test (\* *p* < 0.05, \*\* *p* < 0.005; n = 5 biological replicates). Samples from individuals in the ICU were only included as reference for selected proteins. (**B**) Independent analysis of predicted biomarker A0A5C2FZ03 for anti-Zika virus Ig levels. Antibody levels were compared between PCR– and PCR+ cohorts by Chi-squared test (*p* = 0.004 for asymptomatic cohort). A Spearman's Rho correlation analysis was conducted between virus cross-reactive Ig levels and proteomics-protein relative intensity (*p* < 0.05). (**C**) Independent analysis of predicted biomarkers A0A5C2FZ03 and A0A5C2GPZ0 for anti-HBV Ig levels in asymptomatic and non-severe cohorts, respectively. A Spearman's Rho correlation analysis was conducted between HBV cross-reactive Ig levels and proteomics-protein relative intensity (*p* < 0.05). (**D**) A human-autoantibody general survey identified proteins reactive to IgM and IgG autoantibodies with significant NSI-nor values in PCR– and PCR+ cohorts. Proteins were distributed based on the highest NSI-nor average IgM+IgG value.

associated with autoantibodies and risk of allergic reactions and diseases such as myasthenia gravis and chronic idiopathic thrombocytopenic purpura (ITP) were identified only in non-severe and severe patients (Figure 2C, Supplementary Materials Data File S2). Additionally, autoantibodies were identified in all PCR+ cohorts. The analysis of massspectra relative intensity for selected overrepresented Ig proteins corroborated in individual samples the cohort-dependent results (Figure 3A). Then, an independent analysis of predicted biomarkers was used to validate these results for Zika virus, HBV, and a hu-

In agreement with serum-protein profiles and correlation of Igs overrepresented in vaccinated infected cohorts with protective-associated capacity (Figures 2C and 3A), the results showed that the anti-Zika-virus IgG levels were significantly higher in the asymptomatic cohort only (Figure 3B). Furthermore, a significant positive correlation between Zika-virus cross-reactive Ig levels and relative intensity of associated A0A5C2FZ03 protein was obtained (Figure 3B, Supplementary Materials Data File S2). For HBV, the sensitivity of the surface-antigen ELISA (positive values > 0.034) did not allow for identification of positive individuals in PCR+ cohorts and allowed identification of only one individual in the PCR–
