ORIGINAL RESEARCH article
Enhanced Expression of Autoantigens During SARS-CoV-2 Viral Infection
- 1Sharjah Institute of Medical Research, University of Sharjah, Sharjah, United Arab Emirates
- 2Department of Clinical Sciences, College of Medicine, University of Sharjah, Sharjah, United Arab Emirates
- 3Department of Rheumatology, Ministry of Health and Prevention, Dubai, United Arab Emirates
- 4Division of Surgery and Interventional Science, University College London, London, United Kingdom
- 5Meakins-Christie Laboratories, Research Institute of the McGill University Healthy Center, McGill University, Montreal, QC, Canada
- 6Prince Abdullah Ben Khaled Celiac Disease Research Chair, Department of Pediatrics, Faculty of Medicine, King Saud University, Riyadh, Saudi Arabia
Immune homeostasis is disturbed during severe viral infections, which can lead to loss of tolerance to self-peptides and result in short- or long-term autoimmunity. Using publicly available transcriptomic datasets, we conducted an in-silico analyses to evaluate the expression levels of 52 autoantigens, known to be associated with 24 autoimmune diseases, during SAR-CoV-2 infection. Seven autoantigens (MPO, PRTN3, PADI4, IFIH1, TRIM21, PTPRN2, and TSHR) were upregulated in whole blood samples. MPO and TSHR were overexpressed in both lung autopsies and whole blood tissue and were associated with more severe COVID-19. Neutrophil activation derived autoantigens (MPO, PRTN3, and PADI4) were prominently increased in blood of both SARS-CoV-1 and SARS-CoV-2 viral infections, while TSHR and PTPRN2 autoantigens were specifically increased in SARS-CoV-2. Using single-cell dataset from peripheral blood mononuclear cells (PBMCs), we observed an upregulation of MPO, PRTN3, and PADI4 autoantigens within the low-density neutrophil subset. To validate our in-silico analysis, we measured plasma protein levels of two autoantigens, MPO and PRTN3, in severe and asymptomatic COVID-19. The protein levels of these two autoantigens were significantly upregulated in more severe COVID-19 infections. In conclusion, the immunopathology and severity of COVID-19 could result in transient autoimmune activation. Longitudinal follow-up studies of confirmed cases of COVID-19 could determine the enduring effects of viral infection including development of autoimmune disease.
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus causing coronavirus disease 2019 (COVID-19), appeared first in Wuhan, China, in December 2019 and has since rapidly spread globally (1, 2). COVID-19 severity level ranged from asymptomatic infection to life-threatening disease (3–6). Severe COVID-19 disease has been associated with innate immune dysregulation, early immunosuppression, lymphopenia, and cytokine storm (7–10).
Multiple factors are involved in the development of autoimmunity, including genetics, age, and environment (11). Between the environmental triggers, viral infections, particularly those resulting in low interferon production, as it is the case with SARS-CoV-2 infection, have long been associated with induction of autoimmunity (11). Similar to many severe viral infections, SARS-CoV-2 could trigger the autoimmune reaction through multiple mechanism including molecular mimicry, epitope spreading, bystander activation, and persistence of latent virus (12–15). Aforementioned mechanisms could be understood through examination of homology between various antigens of SARS-CoV-2 and self-antigens (16). Development of cross-reactive epitopes are then dependent on viral strain as well as host genetic susceptibility, including human leucocyte antigen (HLA) polymorphism.
Alternatively, viral induced inflammation and dysregulation of innate and adaptive immune system, could both directly and indirectly lead to loss of tolerance to self-peptides and result in short- or long-term autoimmunity (17, 18). Emerging reports are associating COVID-19 infection with worsening, relapse or de novo induction of several autoimmune diseases including systemic lupus erythematosus (19–21), Guillain-Barré syndrome (22, 23), vasculitis (24, 25), and multiple sclerosis (26). Further, high prevalence of antinuclear antibodies (ANA) and lupus anticoagulant were reported in severe hospitalized SARS-CoV-2 infection, which could hint at induction of autoimmunity (27).
Although the concept of autoimmunity had been explored in previous viral infection (28), its relevance to COVID-19 respiratory infection deserves more attention; especially that immune derangement during SARS-CoV-2 infection could potentially trigger relapse and induction of many new cases. Therefore, the aim of current study was to utilize publicly available transcriptomic COVID-19 data to evaluate autoimmune activation during COVID-19 infection through measuring the gene expression levels of 52 autoantigens, known to be associated with 24 different autoimmune diseases.
Material and Methods
For the purpose of this study, we used a list of 52 autoantigens established by Burbelo et al. (29) for diagnosis of 24 different autoimmune disease including Hashimoto’s thyroiditis, ANCA-associated vasculitis, rheumatoid arthritis (RA), and Systemic lupus erythematosus (SLE) (Table 1). The expression of these autoantigens in the lungs and whole blood of COVID-19 patients was then determined in-silico using publicly available datasets. We validated autoantigen with mRNA expression equal or more than 1.5 LogFC change in leucocyte isolated from COIVD-19. These datasets were publicly available at National Center for Biotechnology Information Gene Expression Omnibus (NCIB GEO, http://www.ncbi.nlm.nih.gov/geo) and the European Bioinformatics Institute (EMBL-EBI, https://www.ebi.ac.uk). Moreover, single cell transcriptomic datasets of sorted neutrophils were used. In addition, the expression of autoantigens following COVID-19 infection was compared to that following infection with three respiratory viruses: SARS-CoV-1, IAV, and RSV.
RNA-sequencing platforms were used for COVID-19 studies, while microarray platforms were used for older datasets of SARS-CoV-1, IAV, and RSV (Table 2). For the COVID-19 lung autopsies dataset (PRJNA646224) (32), the authors have extracted RNA from Formalin fixed paraffin embedded (FFPE) tissues from 9 COVID-19 fatal cases, and 10 SARS-CoV-2-uninfected individuals who undertook biopsy as part of routine clinical care for lung cancer. For this lung autopsy datasets, we used processed sequencing data provided by Wu Meng et al. (32). The authors used DESeq2 to identify differentially expressed genes between the cases and controls. Benjamini-Horchberg correction was used for multiple testing (37).
For COVID-19 whole blood transcriptomic dataset, we used processed sequencing data deposited under project number EGAS00001004503 (33). In this study, Aschenbrenner et al. extracted the RNA from whole blood of 39 COVID-19 patients and 10 healthy controls and analyzed it using NovaSeq 6000 (33). The authors used DESeq2 to identify differentially expressed genes between the cases and controls. Independent hypothesis weighting was used for multiple testing correction (33). Transcriptomic datasets of leucocytes isolated from 10 controls, 51 moderate COVID-19, 49 severe COVID-19 (12 non-critical and 37 critical) were used to validate the results of autoantigens with mRNA expression levels of equal or greater than 1.5 log fold change (GSE157103) (34). For leukocytes study, we processed the raw data using the Bioconductor package limma-voom (38), and presented the results as log2 counts per million (LogCPM). Independent student t-test (39) was used to compare between the independent groups. In addition, for SLE whole blood transcriptomic study, we used processed data provided by Panousis et al. (35). In this study, the authors used DESeq2 to identify differentially expressed genes between 79 active SLE and 58 controls. Benjamini-Horchberg correction was used for multiple testing (35).
Transcriptomic datasets of whole blood isolated from RSV and IAV infected patients (GSE17156) (31) and PBMCs isolated from SARS-CoV-1 infected patients (GSE1739) (30) were analyzed. In both studies, blood was obtained during peak of patient’s symptoms, and processed by authors for RNA extraction and hybridization following Affymetrix protocol. After quality check, we normalized, and log transformed the raw Affymetrix data. Microarray data (CEL files) were pre-processed in our study with Robust Multi-Array Average (RMA) technique using R software (40). The probe set with the largest interquartile range (IQR) of expression values were selected to represent the gene. For RNA-seq study, we processed the data using the Bioconductor package limma-voom (38), and presented the results as logCPM. Log-transformed normalized intensities were also used in Linear Models for MicroArray data (LIMMA) analyses to identify differentially expressed genes between diseased and control groups. We used the default Benjamini-Horchberg correction for multiple testing. Raw data from different studies were never mixed or combined. For each study, the LogFC was obtained separately by analyzing data of diseased and controls. Statistical analyses were performed using R software (v 3.0.2) and Prism (v8; GraphPad Software). For all analyses, p-values <0.05 were considered significant.
For the single cell dataset, transcriptomic datasets of sorted neutrophils were used. Wilk, AJ, et al. (36) performed single sequencing on PBMC from seven COVID-19 patients, and six healthy controls. The details of sample isolation, sequencing, and data processing are available at NCBI GEO, and the study protocols (36). Briefly, PBMCs (GSE150728) were isolated from blood via standard Ficoll-Paque density gradient centrifugation. Authors performed the single-cell RNA-seq library preparation using a Nextera XT DNA library preparation kit (Illumina FC-131-1096) with 1 ng of pooled library and dual-index primers. Sequencing was performed on a NovaSeq S2 instrument (Illumina; Chan Zuckerberg Biohub) (36). Differentially expressed genes (DEGs) were calculated by comparing gene expression of individual COVID-19 samples with gene expression of all healthy controls using Seurat’s implementation of the Wilcoxon rank-sum test. Only DEGs with a two-sided p value <0.05 adjusted for multiple comparisons by Bonferroni’s correction were selected. The investigators clustered neutrophils in to two clusters, low-density neutrophils and canonical neutrophils. The novel cell population of low-density neutrophils was significantly increased only in patients with ARDS.
Gene ontology enrichments analyses was performed using DisGeNET (41) and Gene Ontology biology process databases. Metascape.org (42) was used to identify the enrichment in DisGeNET (41). Terms with a p-value <0.01, a minimum count of 3, and an enrichment factor >1.5 were collected and grouped into clusters based on their membership similarities. The top few enriched clusters (one term per cluster) were presented. Gene Ontology biology process database was accessed through Enrichr open source, available as a gene set enrichment analysis web server (43, 44). GO biological processes were ranked according to combined score. This score was computed in Enrichr by taking the log of the p-value from the Fisher exact test and multiplying that by the z score of the deviation from the expected rank (43, 44).
The plasma level of MPO and PRTN3 (PR3) for seven non-infected controls, eight severe, and eight asymptomatic COVID-19 patients was determined using commercially available human ELISA kit (MPO, Cat # ab119605 and PR3, Cat # ab226902, Abcam, Cambridge, MA, USA). The plasma used in our study was obtained from COVID-19 patients recruited from Rashid Hospital. Plasma was isolated from blood via standard Ficoll-Paque density gradient centrifugation (Sigma, Histopaque-10771). Assays were preformed strictly following the manufacturer’s instructions. Each sample was assayed in duplicate, and values were expressed as the mean of 2 measures per sample. One-way analysis of variance (ANOVA) and post hoc Tukey multiple comparison analyses were applied.
MPO and TSHR autoantigens Are Associated With More Severe COVID-19 in Lung Autopsy and Whole Blood
Using publicly available transcriptomic datasets, we have determined the expression levels of 52 autoantigens, known to be associated with 24 different autoimmune diseases. The list of these genes and their associated autoimmune disease is presented in Table 1. The datasets used in this study are presented in Table 2. Expression levels of autoantigens were determined in lung autopsies and whole blood (Figure 1). For lung, RNA-sequencing data was obtained from 9 deceased COVID-19 patients and 10 negative controls (PRJNA646224) (Figure 1A). For whole blood, RNA-sequencing data was extracted from 10 controls, 19 mild COVID-19, 20 severe COVID-19 (Figures 1B, C). MPO and TSHR were the only two autoantigens overexpressed in lung autopsies and were shared between the lung and blood of severe COVID-19 patients. In whole blood study, seven autoantigens were upregulated in blood of severe COVID-19 patients (MPO, PRTN3, PADI4, IFIH1, TRIM21, PTPRN2, and TSHR) (Figure 1B), while four autoantigens (MPO, PRTN3, IFIH1, and PADI4) were increased in blood of mild patients (Figure 1C). The upregulation of MPO and PRTN3 was significantly higher in severe versus mild COVID-19 (MPO, 3.4 ± 0.3 vs 1.9±0.3 log fold change [LogFC]; p-value = 0.004, and PRTN3 (2.9 ± 0.5 vs 1.6 ± 0.5 LogFC; p-value = 0.01 as presented in Figure 1C). We also looked at the levels of the selected seven autoantigens during active systemic lupus erythematosus (SLE) autoimmune disease using whole blood transcriptomic data deposited by Panousis et al. (35). Of the seven autoantigens, six markers were also upregulated during active SLE, however, MPO (3.4 ± 0.3 vs 0.9 ± 0.2 LogFC; p-value <0.0001), PRTN3 (2.9 ± 0.5 vs 1.2 ± 0.4 LogFC; p-value = 0.007), and PADI4 (1.8 ± 0.3 vs 0.4 ± 0.2 LogFC; p-value = 0.04) levels were significantly higher in severe COVID-19 blood (Figure 1D).
Figure 1 Gene expression of autoantigens in lung autopsies and whole blood of COVID-19 patients. (A) Enhanced expression of MPO and TSHR autoantigens in lung autopsies (n = 9 COVID-19 vs n = 10 controls, Dataset: PRJNA646224). (B) For whole blood COVID-19 dataset (EGAS00001004503), seven autoantigens were upregulated in severe COVID-19 vs controls comparison (C), four autoantigens were upregulated in mild COVID-19 vs controls. The upregulation of MPO and PRTN3 was higher in severe versus mild COVID-19. The sample size presented in B–C EGAS00001004503 dataset was as following; n = 10 controls, n = 19 mild COVID-19, and n = 20 severe COVID-19. Results are presented as log fold change of gene expression with adjusted p-values of less than 0.05. The log fold changes were compared between severe and mild COVID-19 groups. (D) Upregulation of autoantigens in whole blood of both severe COVID-19 and active systemic lupus erythematosus patients. The following datasets were used; EGAS00001004503 (n = 20 severe COVID-19 vs n = 10 controls), and Panousis et al. (35) (n = 79 active SLE vs n = 58 controls). Unpaired student t-test was used to compare between each two independent groups. ∗p < 0.05, ∗∗p < 0.01, ∗∗∗∗p < 0.0001.
We next screened the expression of these seven genes in a second COVID-19 dataset of leukocyte transcriptome (GSE157103) (Figure 2). Out of the seven genes, four autoantigens (MPO, PRTN3, PADI4, and TSHR) were upregulated in leukocytes of COVID-19 patients, while three autoantigens (MPO, PRTN3, and PADI4) were significantly higher in severe versus moderate COVID-19. We have presented the results of leukocyte transcriptomic in Figure 2.
Figure 2 Upregulation of MPO, PRTN3, PADI4, and TSHR mRNA levels in leukocytes of severe COVID-19 patients. The COVID-19 leukocyte data set (GSE157103) analyzed included n = 10 controls, n = 51 moderate COVID-19, and n = 49 severe COVID-19 (12 non-critical and 37 critical). An increase in gene expression of MPO, PRTN3, PADI4, and TSHR in severe COVID-19 compared to moderate and healthy controls. IFIH1 was increased in moderate but not severe COVID-19. Unpaired student t-test was used to compare between the independent groups (mRNA normalized expression between different groups). ∗p < 0.05, ∗∗p < 0.01, ∗∗∗p < 0.001, ∗∗∗∗p < 0.0001. Results are presented as mean (± SEM) of mRNA expression.
We next used gene ontology databases to determine general diseases associated with these seven genes (Figure 3A). DisGeNET database was pooled using publicly available metascape.org tool showing that beside autoimmune disease, these genes were associated with other inflammatory and fibrotic conditions affecting different organs including lung (Bronchiectasis; MPO, PRTN3, and TRIM21) and kidney (Glomerulonephritis; MPO, PRTN3, TRIM21, and PADI4). In addition, the GO biological process database revealed the association of these seven genes with interferon alpha signaling, neutrophil activation, and regulation of cytokine production (Figure 3B).
Figure 3 Gene ontology enrichments identified in the DisGeNET and Gene Ontology biology process databases. (A) Summary of enrichment of DisGeNet category of Metascape. Terms with a p-value <0.01, a minimum count of 3, and an enrichment factor >1.5 (the enrichment factor is the ratio between the observed counts and the counts expected by chance) are collected and grouped into clusters based on their membership similarities. The top few enriched clusters (one term per cluster) were presented. (B) Top 10 categories of the GO biological processes. GO analysis associated with the top seven differentially expressed autoantigens was performed using Enrichr. Bar lines represent cumulative score for the 10 top-ranked categories of GO biological processes. Combined score is computed by taking the log of the p-value from the Fisher exact test and multiplying that by the z score of the deviation from the expected rank.
Upregulation of Autoantigens in Low-Density Neutrophils During COVID-19 Infection
After establishing an overall upregulation of autoantigens in lung autopsies and whole blood COVID-19 patients, we next determined whether the observed increase in autoantigens is reflected on the main inflammatory cells regulating COVID-19 severity. We extracted the neutrophil and lymphocyte counts from PRJNA646224 study and compared the neutrophil to lymphocyte ratio between 16 mild and 15 severe COVID-19 patients. As presented in Figure 4A, severe COVID-19 had significant higher neutrophil to lymphocyte ratio. Next, A single cell dataset of immune cells isolated from peripheral blood mononuclear cells (PBMCs) (GSE150728) of COVID-19 severe patients were used (36). Between the different immune cells three autoantigens, MPO, PRTN3, and PADI4, were significantly enriched within the low-density neutrophil subset (Figure 4B). Moreover, the counts of low-density neutrophil and canonical neutrophil were increased during severe COVID-19 infection (Figure 4C). To validate our in-silico analysis, we next measured plasma protein levels of two autoantigens of MPO and PRTN3 in severe and asymptomatic COVID-19. The protein levels were estimated using human ELISA assays. The results revealed an increase of MPO (mean 36,787 ± 1,961 vs mean 29,007 ± 1,860 pg/ml; p-value = 0.038) and PRTN3 (mean 151.5 ± 38 vs mean 14.77 ± 1.4 ng/ml; p-value = 0.001) in severe compared to asymptomatic COVID-19 infection (Figure 5). Level of these proteins were not different in asymptomatic and non-infected controls (Figure 5).
Figure 4 Neutrophil to lymphocyte ratio and gene expression isolated from severe COVID-19 patients. (A) Neutrophil to lymphocyte ratio from whole blood of 16 mild and 15 severe COIVD-19 patients were compared (PRJNA646224). This ratio was significantly higher in severe COVID-19 (B) Single-cell RNA sequencing was performed on PBMCs from six severe COVID-19 patients and seven healthy controls (GSE150728). Differential expression of MPO, PRTN3, and PADI4 within the low-density neutrophils cluster. (C) Increased counts of low-density neutrophil and canonical neutrophil during severe COVID-19 infection. The count value represents the absolute cell number within the PBMCs. Unpaired student t-test was used to compare between the independent groups. ∗p < 0.05, ∗∗p < 0.01, ∗∗∗p < 0.001.
Figure 5 Elevated plasma MPO and PRTN3 levels in severe COIVD-19 patients. An increase in the plasma levels of MPO and PRTN3 in severe compared to asymptomatic COVID-19 patients. The protein levels were estimated using human ELISA assays (n = non-infected controls, n = 8 asymptomatic COVID-19, and severe n = 8 COVID-19). ∗p < 0.05, ∗∗p < 0.01, ∗∗∗p < 0.001.
Prominent Autoantigens Upregulation in Coronavirus Infections Relative to Other Viral Infections
We next compared the profile of autoantigen upregulation observed during SARS-CoV-2 to that detected during other respiratory viral infection. To do that, we used transcriptomic microarrays and RNA-sequencing data from blood of SARS-CoV-1, influenza A virus (IAV), and respiratory syncytial virus (RSV) infected patients at the peak of disease. For each condition, differential expression and LogFC were obtained by comparing the normalized gene expression of the infected group versus healthy donors (Figure 6A). None of autoantigens were upregulated more than one LogFC in IAV and RSV, while five autoantigens in SARS-CoV-1 and seven autoantigens in SARS-CoV-2 were upregulated more than one LogFC (Figure 6A). MPO, PRTN3 (PR3), and PADI4 were the top shared autoantigens appearing in both coronavirus respiratory infections, with an increase in expression of more than 1.5 LogFC. We then intersected the differentially expressed genes in all four respiratory infections to obtain the shared signatures (Figure 6B). Interestingly, TSHR and PTPRN2 autoantigens were specifically increased in SARS-CoV-2 (Figure 6B). TSHR was also overexpressed in COVID-19 lung autopsies. Two genes (IFIH1 and TRIM21) were upregulated in three viral infections; however, their expression was higher in SARS-CoV-2 compared to IAV, and RSV (Figure 6B).
Figure 6 Expression of autoantigens during COVID-19 and other viral infections. The following datasets were used to compare expression of autoantigens between SARS-CoV and other viral infections; GSE17156 (n = 17 IAV vs n = 17 controls), GSE17156 (n = 20 RSV vs n = 20 controls), GSE1739 (n = 10 SARS-CoV-1 vs n = 4 controls), and EGAS00001004503 (n = 39 COVID-19 vs n = 10 controls). (A) Log fold change in the expression of autoantigens following infection with SARS-CoV-1, SARS-CoV-2, IAV, and RSV relative to healthy controls. (B) Intersection of upregulated autoantigen signatures in four different respiratory viral infections, IAV, RSV, SARS-CoV-1, and SARS-CoV-2. Three top upregulated autoantigens (PADI4, MPO, and PRTN3) are shared between SARS-COV-1 and SARS-COV-2. For all analyses, p < 0.05 was considered significant. IAV, influenza A virus; RSV, Respiratory syncytial virus.
Herein we report the common autoantigens during COVID-19 viral infection. Fifty-two autoantigens, which are linked with 24 different autoimmune diseases, were tested (29). Seven autoantigens were found to be elevated in COVID-19. The level of three autoantigens, MPO, PRTN3, and PADI4, were higher in the blood of severe compared to mild COVID-19.
Autoimmune disease could be triggered by genetic and environmental factors; viral infections had been known as a major environmental cause of transient autoimmunity that could potentially lead to relapse or induction of de novo autoimmune disorders (11). These autoimmune disorders emerge weeks post viral infection; hence sensitive serological tests are needed to determine the cause-effect relationship between SARS-CoV-2 infection and autoimmune disease diagnosis (45–47).
Interestingly, MPO and TSHR were increased in both lung autopsies and whole blood of severe COVID-19 patients. Comparison of COVID-19 blood transcriptomic with IAV and RSV revealed that MPO, PRTN3, and PADI4 were selectively upregulated in coronavirus infections, SARS-CoV-1 and SARS-CoV-2, while TSHR and PTPRN2 autoantigens were distinctive to SARS-CoV-2 infection. These two genes (TSHR and PTPRN2) did not increase in the mild COVID-19 and hence they were associated with more severe infection.
Following the results obtained through gene ontology enrichment analyses, single cell transcriptomics of PBMCs revealed the significant increase in MPO, PRTN3, and PADI4 mRNA levels within low-density neutrophils. In addition, analyses of cell counts provided by Wilk et al. study showed significant increase in both low-density and canonical neutrophils during severe COVID-19 infection compared to controls (Figure 4C) (36).
The top seven autoantigens upregulated in the blood of severe COVID-19 patients, were associated with a wide range of vascular and inflammatory autoimmune disorders (Figure 3A). Vasculitis, a shared condition between fatal COVID-19 and vascular autoimmune diseases such as Anti-neutrophil cytoplasmic autoantibody (ANCA) vasculitis, is featured by elevation of MPO and PRTN3 levels (25, 48). Wilk et al. single cell data identified a distinct group of low-density neutrophil; these immune cells were only detected in severe COVID-19 complicated with acute respiratory distress syndrome (ARDS) (36). These cells had a significantly high level of MPO, PRTN3, and PADI4 indicating that they could be a major source of the observed increase in blood level of these autoantigens during severe COVID-19. Following SARS-CoV-2 infection, neutrophil-derived extracellular traps (NETs) formation by neutrophils, NETosis, may therefore lead to burst of autoantigens including PADI4, MPO, and PRTN3 in the context of immunostimulatory molecules.
Confirming previous findings (49, 50), increased MPO level reflected COVID-19 severity. Neutrophil activation leading to net formation or NETosis is observed in both viral infections and autoimmune diseases (47, 51, 52). NETosis could then be considered as a common pathological modulator of viral infection and autoimmune disease (47, 52, 53). NETosis autoantigens, MPO, PRTN3, and PADI4, were markedly increased in SARS-CoV-1 and SARS-CoV-2 but they did not appear in IAV and RSV. Supporting our findings, a recent investigation by FP Veras et al. showed viable SARS-CoV-2 ability to directly induce the release of NETs by healthy neutrophils (54). While, in another study, influenza A virus infection did not affect MPO release (55). MPO is a peroxidase enzyme responsible for intracellular catalytic reactions between hydrogen peroxidase and chlorides to form hypochlorous acid (56, 57). NETosis leads to extracellular burst of chromatic, histones, and neutrophil granules containing MPO, PRTN3, and PADI4. An exaggerated increase in the level of these antigens during COVID-19 infection could lead to breaking the autoimmune tolerance and the recognition of autoantigens by immune sentinel cells.
NETosis induces inflammation and DNA deployment that could trigger autoimmunity. While NET produced by neutrophils, their clearance is achieved by macrophages efferocytosis. During ARDS, neutrophil count and lifespan is significantly increased, and the ability of macrophages to engulf NETs and apoptotic cells is significantly decreased, prolonging the lung injury induced by neutrophil blast. Pharmacologic treatment could be used to enhance NETs clearance. Macrophage NETs efferocytosis could be restored by AMP-activated protein kinase (AMPK) activator such as Metformin or application of neutralizing antibody against HMGB1 (58). Therefore, such medications could be considered to reduce blood levels of these autoantigens, and hence lower chance of triggering autoimmunity.
TSHR is expressed by thyroid epithelial cells, and various extra-thyroidal tissue including the adipose, peripheral blood cells, and fibrocytes (59). Derived from monocytes, human fibrocytes express both thyroglobulin and thyrotropin receptor (60). They are increased during lung injury and have both the inflammatory characters of macrophages and the tissue remodeling features of fibroblasts (61). Chronic inflammatory conditions such as autoimmunity, cardiovascular disease, and asthma promote differentiation of immune cells to circulating fibrocytes and their accumulation at the site of injury (61, 62).
TSHR is targeted by autoantibodies during Graves’ disease (63). SARS-CoV-2 infection has been connected with the initiation and relapse phases of Grave disease (45, 46). Of note, this disease has emerged in some patients during COVID-19 recovery period. Patients diagnosed were negative for naso-pharyngeal swab PCR test but were positive for both IgM and IgG SARS-CoV-2 antibodies. This suggests that the observed increase of autoantigens during SARS-CoV-2 infection may trigger autoimmunity. This could lead to initiation or relapse of autoimmune disorders as a long-term COVID-19 outcome. In fact, evidence of transient autoimmunity has been reported among long COVID-19 outcomes by several studies (64–66). Follow-up data from survivors of viral infections have shown appearance of autoimmune disorders within weeks to months after recovery (67). In addition, TSHR elevation was only detected in severe COVID-19, which could suggest higher chance for appearance of autoimmune disorders post severe COVID-19 infection.
The majority of COVID-19 patients are expected to show one or more residual symptoms months after recovering from the infection (66, 68–70). In a post-acute COVID-19 follow-up study of 179 confirmed cases in Italy, fatigue and dyspnea persisted in around half of recovered individuals, while joint pain and chest pain lingered in 20–30% of recovered patients (68). Similarly, a 6-month follow-up study of 1,733 COVID-19 hospitalized patients from China reported lasting of fatigue and muscle weakness in more than half of patients, while patients who were more severely ill during their hospital stay had more persisting long-term symptoms (71).
Emerging case reports have shown that SARS-CoV-2 induces long-term immune-inflammatory abnormalities (66, 72). Schenker et al. reported a 65-year-old female patient with de novo reactive arthritis and cutaneous vasculitis 10 days after recovery of all COVID-19-related symptoms (72). In a more severe case, a 31-year female was presented with fatal multisystem inflammatory syndrome (MIS) 2 weeks post recovery from COVID-19 (65). This increase in incidence of autoimmunity was also reported among children, where SARS-CoV-2 epidemic was associated with a 30-fold increase in Kawasaki-like disease (71).
Long-COVID-19 persisting symptoms involve immune-mediated inflammatory disease and neurological abnormalities that could suggest possibility of triggering pre-existing or de novo autoimmune reactions weeks or month after COVID-19 recovery (37, 73). Previous studies have shown that autoantigen gene upregulation is often followed by an increase in the respective autoantibody level (29, 48, 74). Although the increase in autoantigen expression observed in this study could only trigger short lasting autoimmunity, follow-up longitudinal studies are needed to establish the long-term enduring effects of SARS-CoV-2 infection in developing autoimmune diseases.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material.
The studies involving human participants were reviewed and approved by Dubai Health Authority. The patients/participants provided their written informed consent to participate in this study.
RaH, QH, NSS-A, FSS-A, conceived and designed the experiments. NSS-A, FSS-A, RaH, RiH, SH and SA analyzed the data. All authors contributed to the article and approved the submitted version.
This research has been financially supported by Tissue Injury and Repair (TIR) group operational grant (Grant code: 150317); COVID-19 research grant (CoV19-0307) and (CoV19-0308), Seed grant (Grant code: 2001090275); and by collaborative research grant (Grant code: 2001090278) to RH, University of Sharjah, UAE; and by a Sandooq Al Watan Applied Research & Development grant to RH (SWARD-S20-007); and by Al Jalila Foundation Seed Grant (AJF202019); and by Prince Abdullah Ben Khalid Celiac Disease Research Chair, under the Vice Deanship of Research Chairs, King Saud University, Riyadh, Kingdom of Saudi Arabia.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
ANA, Antinuclear antibodies; ANCA, Anti-neutrophil cytoplasmic autoantibody; ARDS, Acute respiratory distress syndrome; COVID-19, Coronavirus disease of 2019; DEGs, Differentially expressed genes; EMBL-EBI, European Bioinformatics Institute; FFPE, Formalin fixed paraffin embedded; HLA, Human leukocyte antigen; IAV, Influenza A virus; IQR, Interquartile range; LIMMA, Linear Models for MicroArray data; Log CPM, Log2 counts per million; MPO, Myeloperoxidase; NETs, Neutrophil-derived extracellular traps; NCIB GEO, National Center for Biotechnology Information Gene Expression Omnibus; PBMC, Peripheral blood mononuclear cells; RA, Rheumatoid arthritis; RMA, Robust Multi-Array Average; RSV, Respiratory syncytial virus; SARS-CoV-2, Severe acute respiratory syndrome coronavirus 2; SLE, Systemic lupus erythematosus.
3. Chen H, Guo J, Wang C, Luo F, Yu X, Zhang W, et al. Clinical Characteristics and Intrauterine Vertical Transmission Potential of COVID-19 Infection in Nine Pregnant Women: A Retrospective Review of Medical Records. Lancet (2020) 395(10226):809–15. doi: 10.1016/S0140-6736(20)30360-3
4. Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, et al. Clinical Features of Patients Infected With 2019 Novel Coronavirus in Wuhan, China. Lancet (2020) 395(10223):497–506. doi: 10.1016/S0140-6736(20)30183-5
5. Rothe C, Schunk M, Sothmann P, Bretzel G, Froeschl G, Wallrauch C, et al. Transmission of 2019-Ncov Infection From an Asymptomatic Contact in Germany. N Engl J Med (2020) 382(10):970–1. doi: 10.1056/NEJMc2001468
8. Blanco-Melo D, Nilsson-Payant BE, Liu W-C, Uhl S, Hoagland D, Møller R, et al. Imbalanced Host Response to SARS-CoV-2 Drives Development of COVID-19. Cell (2020) 181(5):1036–45.e1039. doi: 10.1016/j.cell.2020.04.026
9. Lee JS, Park S, Jeong HW, Ahn JY, Choi SJ, Lee H, et al. Immunophenotyping of COVID-19 and Influenza Highlights the Role of Type I Interferons in Development of Severe COVID-19. Sci Immunol (2020) 5(49):eabd1554. doi: 10.1126/sciimmunol.abd1554
12. Fujinami RS, von Herrath MG, Christen U, Whitton JL. Molecular Mimicry, Bystander Activation, or Viral Persistence: Infections and Autoimmune Disease. Clin Microbiol Rev (2006) 19(1):80–94. doi: 10.1128/CMR.19.1.80-94.2006
13. Reddy V, Desai A, Krishna SS, Vasanthapuram R. Molecular Mimicry Between Chikungunya Virus and Host Components: A Possible Mechanism for the Arthritic Manifestations. PLoS Negl Trop Dis (2017) 11(1):e0005238. doi: 10.1371/journal.pntd.0005238
16. Megremis S, Walker TDJ, He X, Ollier WER, Chinoy H, Hampson L, et al. Antibodies Against Immunogenic Epitopes With High Sequence Identity to SARS-CoV-2 in Patients With Autoimmune Dermatomyositis. Ann Rheumatic Dis (2020) 79(10):1383. doi: 10.1136/annrheumdis-2020-217522
17. Canas CA. The Triggering of post-COVID-19 Autoimmunity Phenomena Could be Associated With Both Transient Immunosuppression and an Inappropriate Form of Immune Reconstitution in Susceptible Individuals. Med Hypotheses (2020) 145:110345. doi: 10.1016/j.mehy.2020.110345
18. Bastard P, Rosen LB, Zhang Q, Michailidis E, Hoffmann HH, Zhang Y, et al. Autoantibodies Against Type I Ifns in Patients With Life-Threatening COVID-19. Science (2020) 370(6515). doi: 10.1126/science.abd4585
19. Mantovani Cardoso E, Hundal J, Feterman D, Magaldi J. Concomitant New Diagnosis of Systemic Lupus Erythematosus and COVID-19 With Possible Antiphospholipid Syndrome. Just a Coincidence? A Case Report and Review of Intertwining Pathophysiology. Clin Rheumatol (2020) 39(9):2811–5. doi: 10.1007/s10067-020-05310-1
20. Slimani Y, Abbassi R, El Fatoiki F-Z, Barrou L, Chiheb S. Systemic Lupus Erythematosus and Varicella-Like Rash Following COVID-19 in a Previously Healthy Patient. J Med Virol (2020) 93(2):1184–7. doi: 10.1002/jmv.26513
23. Alberti P, Beretta S, Piatti M, Karantzoulis A, Piatti ML, Santoro P, et al. Guillain-Barré Syndrome Related to COVID-19 Infection. Neurol - Neuroimmunol Neuroinflamm (2020) 7(4):e741. doi: 10.1212/NXI.0000000000000741
25. Vacchi C, Meschiari M, Milic J, Marietta M, Tonelli R, Alfano G, et al. Covid-19-associated Vasculitis and Thrombotic Complications: From Pathological Findings to Multidisciplinary Discussion. Rheumatology (2020) 59(12):e147–50. doi: 10.1093/rheumatology/keaa581
26. Zanin L, Saraceno G, Panciani PP, Renisi G, Signorini L, Migliorati K, et al. Sars-CoV-2 Can Induce Brain and Spine Demyelinating Lesions. Acta Neurochir (2020) 162(7):1491–4. doi: 10.1007/s00701-020-04374-x
27. Gazzaruso C, Carlo Stella N, Mariani G, Nai C, Coppola A, Naldani D, et al. High Prevalence of Antinuclear Antibodies and Lupus Anticoagulant in Patients Hospitalized for SARS-CoV2 Pneumonia. Clin Rheumatol (2020) 39(7):2095–7. doi: 10.1007/s10067-020-05180-7
29. Burbelo PD, Iadarola MJ, Alevizos I, Sapio MR. Transcriptomic Segregation of Human Autoantigens Useful for the Diagnosis of Autoimmune Diseases. Mol Diagn Ther (2016) 20(5):415–27. doi: 10.1007/s40291-016-0211-6
30. Reghunathan R, Jayapal M, Hsu LY, Chng HH, Tai D, Leung BP, et al. Expression Profile of Immune Response Genes in Patients With Severe Acute Respiratory Syndrome. BMC Immunol (2005) 6:2. doi: 10.1186/1471-2172-6-2
31. Zaas AK, Chen M, Varkey J, Veldman T, Hero AO 3rd, Lucas J, et al. Gene Expression Signatures Diagnose Influenza and Other Symptomatic Respiratory Viral Infections in Humans. Cell Host Microbe (2009) 6(3):207–17. doi: 10.1016/j.chom.2009.07.006
32. Wu M, Chen Y, Xia H, Wang C, Tan CY, Cai X, et al. Transcriptional and Proteomic Insights Into the Host Response in Fatal COVID-19 Cases. Proc Natl Acad Sci (2020) 117(45):28336. doi: 10.1073/pnas.2018030117
33. Aschenbrenner AC, Mouktaroudi M, Krämer B, Oestreich M, Antonakos N, Nuesch-Germano M, et al. Disease Severity-Specific Neutrophil Signatures in Blood Transcriptomes Stratify COVID-19 Patients. Genome Med (2021) 13(1):7. doi: 10.1186/s13073-020-00823-5
34. Overmyer KA, Shishkova E, Miller IJ, Balnis J, Bernstein MN, Peters-Clarke TM, et al. Large-Scale Multi-Omic Analysis of COVID-19 Severity. Cell Syst (2020) 12(1):23–40. doi: 10.1016/j.cels.2020.10.003
35. Panousis NI, Bertsias GK, Ongen H, Gergianaki I, Tektonidou MG, Trachana M, et al. Combined Genetic and Transcriptome Analysis of Patients With SLE: Distinct, Targetable Signatures for Susceptibility and Severity. Ann Rheumatic Dis (2019) 78(8):1079. doi: 10.1136/annrheumdis-2018-214379
36. Wilk AJ, Rustagi A, Zhao NQ, Roque J, Martínez-Colón GJ, McKechnie JL, et al. A Single-Cell Atlas of the Peripheral Immune Response in Patients With Severe COVID-19. Nat Med (2020) 26(7):1070–6. doi: 10.1038/s41591-020-0944-y
38. Ritchie ME, Phipson B, Wu DI, Hu Y, Law CW, Shi W, et al. Limma Powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies. Nucleic Acids Res (2015) 43(7):e47–7. doi: 10.1093/nar/gkv007
41. Piñero J, Bravo À, Queralt-Rosinach N, Gutiérrez-Sacristán A, Deu-Pons J, Centeno E, et al. DisGeNET: A Comprehensive Platform Integrating Information on Human Disease-Associated Genes and Variants. Nucleic Acids Res (2017) 45(D1):D833–9. doi: 10.1093/nar/gkw943
42. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape Provides a Biologist-Oriented Resource for the Analysis of Systems-Level Datasets. Nat Commun (2019) 10(1):1523. doi: 10.1038/s41467-019-09234-6
43. Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, et al. Enrichr: Interactive and Collaborative HTML5 Gene List Enrichment Analysis Tool. BMC Bioinf (2013) 14:128. doi: 10.1186/1471-2105-14-128
44. Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, et al. Enrichr: A Comprehensive Gene Set Enrichment Analysis Web Server 2016 Update. Nucleic Acids Res (2016) 44(W1):W90–7. doi: 10.1093/nar/gkw377
45. Mateu-Salat M, Urgell E, Chico A. Sars-COV-2 as a Trigger for Autoimmune Disease: Report of Two Cases of Graves’ Disease After COVID-19. J Endocrinol Invest (2020) 43(10):1527–8. doi: 10.1007/s40618-020-01366-7
48. Jones BE, Herrera CA, Agosto-Burgos C, Starmer J, Bass WA, Poulton CJ, et al. ANCA Autoantigen Gene Expression Highlights Neutrophil Heterogeneity Where Expression in Normal-Density Neutrophils Correlates With ANCA-induced Activation. Kidney Int (2020) 98(3):744–57. doi: 10.1016/j.kint.2020.04.037
50. Chen T, Wu D, Chen H, Yan W, Yang D, Chen G, et al. Clinical Characteristics of 113 Deceased Patients With Coronavirus Disease 2019: Retrospective Study. BMJ (Clin Res ed.) (2020) 368:m10911. doi: 10.1136/bmj.m1091
52. Turunen S, Huhtakangas J, Nousiainen T, Valkealahti M, Melkko J, Risteli J, et al. Rheumatoid Arthritis Antigens Homocitrulline and Citrulline are Generated by Local Myeloperoxidase and Peptidyl Arginine Deiminases 2, 3 and 4 in Rheumatoid Nodule and Synovial Tissue. Arthritis Res Ther (2016) 18(1):239. doi: 10.1186/s13075-016-1140-9
53. Khandpur R, Carmona-Rivera C, Vivekanandan-Giri A, Gizinski A, Yalavarthi S, Knight JS, et al. Nets are a Source of Citrullinated Autoantigens and Stimulate Inflammatory Responses in Rheumatoid Arthritis. Sci Transl Med (2013) 5(178):178ra40. doi: 10.1126/scitranslmed.3005580
54. Veras FP, Pontelli MC, Silva CM, Toller-Kawahisa JE, de Lima M, Nascimento DC, et al. Sars-CoV-2-triggered Neutrophil Extracellular Traps Mediate COVID-19 Pathology. J Exp Med (2020) 217(12). doi: 10.1084/jem.20201129
55. Pang G, Clancy R, Cong M, Ortega M, Zhigang R, Reeves G. Influenza Virus Inhibits Lysozyme Secretion by Sputum Neutrophils in Subjects With Chronic Bronchial Sepsis. Am J Respir Crit Care Med (2000) 161(3 Pt 1):718–22. doi: 10.1164/ajrccm.161.3.9812047
57. Wallace FJ, Clancy RL, Cripps AW. An Animal Model Demonstration of Enhanced Clearance of Nontypable Haemophilus Influenzae From the Respiratory Tract After Antigen Stimulation of Gut-Associated Lymphoid Tissue. Am Rev Respir Dis (1989) 140(2):311–6. doi: 10.1164/ajrccm/140.2.311
58. Grégoire M, Uhel F, Lesouhaitier M, Gacouin A, Guirriec M, Mourcin F, et al. Impaired Efferocytosis and Neutrophil Extracellular Trap Clearance by Macrophages in ARDS. Eur Respir J (2018) 52(2):1702590. doi: 10.1183/13993003.02590-2017
60. Fernando R, Atkins S, Raychaudhuri N, Lu Y, Li B, Douglas RS, et al. Human Fibrocytes Coexpress Thyroglobulin and Thyrotropin Receptor. Proc Natl Acad Sci (2012) 109(19):7427. doi: 10.1073/pnas.1202064109
64. Bastida C, Also MA, Pericas JM, Letang E, Tuset M, Miró JM. [Rhabdomyolysis and Severe Hepatotoxicity Due to a Drug-Drug Interaction Between Ritonavir and Simvastatin. Could We Use the Most Cost-Effective Statin in All Human Immunodeficiency Virus-Infected Patients?]. Enferm Infecc Microbiol Clin (2014) 32(9):579–82. doi: 10.1016/j.eimc.2014.03.014
66. Sollini M, Ciccarelli M, Cecconi M, Aghemo A, Morelli P, Gelardi F, et al. Vasculitis Changes in COVID-19 Survivors With Persistent Symptoms: An [18F]FDG-PET/CT Study. Eur J Nucl Med Mol Imaging (2021) 48(5):1460–6. doi: 10.1007/s00259-020-05084-3
69. Zhao Y-M, Shang Y-m, Song W-b, Li Q-q, Xie H, Xu Q-f, et al. Follow-Up Study of the Pulmonary Function and Related Physiological Characteristics of COVID-19 Survivors Three Months After Recovery. EClinicalMedicine (2020) 25. doi: 10.1016/j.eclinm.2020.100463
70. van den Borst B, Peters JB, Brink M, Schoon Y, Bleeker-Rovers CP, Schers H, et al. Comprehensive Health Assessment Three Months After Recovery From Acute COVID-19. Clin Infect Dis (2020). doi: 10.1093/cid/ciaa1750
71. Verdoni L, Mazza A, Gervasoni A, Martelli L, Ruggeri M, Ciuffreda M, et al. An Outbreak of Severe Kawasaki-Like Disease at the Italian Epicentre of the SARS-CoV-2 Epidemic: An Observational Cohort Study. Lancet (2020) 395(10239):1771–8. doi: 10.1016/S0140-6736(20)31103-X
Keywords: autoantigen, bioinformatics, SARS-CoV-2, COVID-19, autoimmune disease, neutrophil, respiratory viral infection, lung autopsies
Citation: Saheb Sharif-Askari N, Saheb Sharif-Askari F, Ahmed SBM, Hannawi S, Hamoudi R, Hamid Q and Halwani R (2021) Enhanced Expression of Autoantigens During SARS-CoV-2 Viral Infection. Front. Immunol. 12:686462. doi: 10.3389/fimmu.2021.686462
Received: 26 March 2021; Accepted: 27 May 2021;
Published: 30 June 2021.
Edited by:Zhiwei Wu, Nanjing University, China
Reviewed by:Laura Cook, Peter Doherty Institute for Infection and Immunity, Australia
Marcelo A. Soares, National Cancer Institute (INCA), Brazil
Copyright © 2021 Saheb Sharif-Askari, Saheb Sharif-Askari, Ahmed, Hannawi, Hamoudi, Hamid and Halwani. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rabih Halwani, email@example.com