Abstract
Population genetic variability in immune system genes can often underlie variability in immune responses to pathogens. Cytotoxic T-lymphocytes are emerging as critical determinants of both severe acute respiratory syndrome coronavirus 2 infection severity and long-term immunity, after either recovery or vaccination. A hallmark of coronavirus disease 2019 is its highly variable severity and breadth of immune responses between individuals. To address the underlying mechanisms behind this phenomenon, we analyzed the proteolytic processing of S1 spike glycoprotein precursor antigenic peptides across ten common allotypes of endoplasmic reticulum aminopeptidase 1 (ERAP1), a polymorphic intracellular enzyme that can regulate cytotoxic T-lymphocyte responses by generating or destroying antigenic peptides. We utilized a systematic proteomic approach that allows the concurrent analysis of hundreds of trimming reactions in parallel, thus better emulating antigen processing in the cell. While all ERAP1 allotypes were capable of producing optimal ligands for major histocompatibility complex class I molecules, including known severe acute respiratory syndrome coronavirus 2 epitopes, they presented significant differences in peptide sequences produced, suggesting allotype-dependent sequence biases. Allotype 10, previously suggested to be enzymatically deficient, was rather found to be functionally distinct from other allotypes. Our findings suggest that common ERAP1 allotypes can be a major source of heterogeneity in antigen processing and through this mechanism contribute to variable immune responses in coronavirus disease 2019.