Investigation of Immunoreactivity Profiles and Epitope Landscape in Divergent COVID-19 Trajectories and SARS-CoV-2 Variants

Abstract

This study aimed to elucidate the complexity of the humoral immune response in COVID-19 patients with varying disease trajectories using a SARS-CoV-2 whole proteome peptide microarray chip. The microarray, containing 5347 peptides spanning the entire SARS-CoV-2 proteome and key variants of concern, was used to analyze IgG responses in 10 severe-to-recovered, 9 nonsevere-to-severe cases, and 10 control case (5 pre-pandemic and 5 SARS-CoV-2-negative) plasma samples. We identified 1151 IgG-reactive peptides corresponding to 647 epitopes, with 207 peptides being cross-reactive across 124 epitopes. Nonstructural protein 3 (nsp3) exhibited the highest number of total and unique epitopes, followed by the spike protein. nsp12 had the most number of cross-reactive epitopes. Peptides from the spike protein and nsps 2, 3, 5, and 13 were notably associated with recovery. Additionally, specific mutations in SARS-CoV-2 variants were found to alter peptide immunoreactivity, with some mutations (e.g., G142D, L452R, and N501Y) enhancing and others (e.g., R190S and E484 K) reducing immune recognition. These findings have critical implications for the development of diagnostics, vaccines, and therapeutics. Understanding the distribution of epitopes and the impact of viral mutations on antigenicity provides insights into immune evasion mechanisms, informing strategies for controlling COVID-19 and future coronavirus outbreaks.

Keywords: COVID-19; IgG immune response; Kruskal–Wallis test; Mann–Whitney U test; SARS-CoV-2; cross-reactivity; epitope; mutation; peptide microarray; z-score.