{"id":13788,"date":"2025-12-16T06:00:00","date_gmt":"2025-12-16T11:00:00","guid":{"rendered":"https:\/\/cov19longhaulfoundation.org\/?p=13788"},"modified":"2025-12-16T09:38:50","modified_gmt":"2025-12-16T14:38:50","slug":"heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin-of-severe-covid-19-cases","status":"publish","type":"post","link":"https:\/\/cov19longhaulfoundation.org\/?p=13788","title":{"rendered":"Heightened Innate Immune and Inflammatory Gene Expression in the Skin of Severe COVID-19 Cases"},"content":{"rendered":"\n<p class=\"has-small-font-size\"><strong>Authors<\/strong>\u00a0Fernandes IG,\u00a0Ferreira FM,\u00a0De Sousa ESA\u00a0<a href=\"https:\/\/orcid.org\/0000-0001-8861-635X\" target=\"_blank\" rel=\"noreferrer noopener\"><\/a>,\u00a0<a href=\"https:\/\/www.dovepress.com\/author_profile.php?id=2271883\" target=\"_blank\" rel=\"noreferrer noopener\">P<\/a>et.al., <a href=\"https:\/\/orcid.org\/0000-0002-7911-1824\" target=\"_blank\" rel=\"noreferrer noopener\"><\/a><strong>DOI<\/strong>\u00a0<a href=\"https:\/\/doi.org\/10.2147\/JIR.S516737\">https:\/\/doi.org\/10.2147\/JIR.S516737<\/a><\/p>\n\n\n\n<p>Correspondence: Maria Notomi Sato, Laboratory of Medical Investigation in Dermatology and Immunodeficiencies (LIM-56), Department of Dermatology, Faculdade de Medicina, Universidade de Sao Paulo, Sao Paulo, SP, 05403-000, Brazil, Tel +55 11 30617499, Email\u00a0marisato@usp.br<br><br><strong>Purpose:<\/strong>\u00a0Cutaneous manifestations of SARS-CoV-2 infection exhibit significant variability, yet the role of innate immune responses in the skin of COVID-19 patients remains poorly understood. In this study, we investigated the transcriptomic profile of skin samples from patients who succumbed to COVID-19.<\/p>\n\n\n\n<p><strong>Patients and Methods:<\/strong>\u00a0Skin autopsies from COVID-19 patients with post-mortem time of less than 20\u00a0hours were obtained from University of S\u00e3o Paulo Medical School Hospital and healthy skin samples, were submitted to RNA sequencing analysis. Validation of differentially\u00a0expressed genes (DEGs)\u00a0was performed by real-Time PCR.<\/p>\n\n\n\n<p><strong>Results:<\/strong>\u00a0Our analysis revealed markedly elevated expression of type I interferon (IFN)-inducible antiviral factors, antioxidant enzymes, and components of several cytokine-signaling pathways in COVID-19 skin samples compared to healthy controls. SARS-CoV-2 infection robustly induced numerous interferon-stimulated genes (ISGs), IFITs, IRF, S100 family with a notable enrichment of those associated with antiviral and inflammatory responses. Moreover, the presence of counter-regulatory factors such as SOCS3 and NFKBIA indicate the involvement of anti-inflammatory mechanisms in the skin. Furthermore, deconvolution data indicated increased presence of macrophages other nucleated cells in vessels in the skin of COVID-19 patients, highlighting the involvement of innate immune mechanisms.<br><strong>Conclusion:<\/strong>\u00a0The results revealed cutaneous alterations in the expression of genes associated with innate immunity and inflammation factors. This suggests that, unlike tissues with viral tropism, the skin is enriched with antiviral factors to defend against SARS-CoV-2. This information could be useful for developing specific antiviral therapies.<\/p>\n\n\n\n<h2 class=\"wp-block-heading\">Introduction<\/h2>\n\n\n\n<p>The cutaneous manifestations associated with COVID-19 vary widely and can be broadly classified into two main categories: inflammatory reactions (eg, morbilliform&nbsp;(measles-like)&nbsp;rash, vascular&nbsp;rashes&nbsp;and&nbsp;vesicular rash) and vascular-origin lesions (eg, chilblain-like rashes, petechiae\/purpura, and livedo racemosa-like patterns).<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0001\">1<\/a>,<a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0002\">2<\/a><\/sup>&nbsp;Vesicular eruptions appear early in the course of the disease (15% before other symptoms). The pseudo\u2010chilblain pattern frequently appears late in the evolution of the COVID\u201019 disease (59% after other symptoms), while the rest tend to appear with other symptoms of COVID\u201019.<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0003\">3<\/a><\/sup>&nbsp;Despite numerous reports of diverse cutaneous manifestations linked to SARS-CoV-2 infection, dermatologists treating COVID-19 patients at the University of S\u00e3o Paulo Medical School Hospital\u2014a major reference center and one of the largest university hospitals in Latin America confirmed that was not possible to attribute with certainty that any cutaneous lesions had a direct correlation to the SARS-CoV-2 infection in 86 patients confirmed with COVID-19.<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0004\">4<\/a><\/sup>&nbsp;Instead, the dermatological presentations observed mirrored those commonly encountered prior to the pandemic, with fungal infections being the most frequent, followed by drug eruptions and viral infections. This underscores the need for further molecular investigations to unravel the immunological changes occurring in the skin during COVID-19. In same line, it has been verified a lack of association between chilblains and SARS-CoV-2 infection.<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0005\">5<\/a><\/sup><\/p>\n\n\n\n<p>SARS-CoV-2 enters host cells via the angiotensin-converting enzyme 2 (ACE2) receptor and TMPRSS2 protease, which are expressed in multiple tissues, including the lungs, kidneys, colon, heart, testicles, and skin.<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0006\">6<\/a><\/sup>&nbsp;Despite this, studies have found no virally induced cytopathic alterations or intranuclear inclusions in the skin,<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0007\">7<\/a><\/sup>&nbsp;nor detectable SARS-CoV-2 viral load in skin biopsy specimens.<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0008\">8<\/a>,<a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0009\">9<\/a><\/sup>&nbsp;Nevertheless, upon viral entry, after the latent period, depending on the variant, such as Delta, the latency period was 3.2 days.<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0010\">10<\/a><\/sup>&nbsp;After this, a cascade of activation of innate immunity is rapidly activated, including the recruitment of dendritic cells, macrophages, and monocytes. These immune cells adhere to infected cells and release type I interferons (IFNs), which act as danger signals, along with pro-inflammatory and antiviral cytokines.<\/p>\n\n\n\n<p>Understanding the molecular immunological events in the skin during SARS-CoV-2 infection requires a focus on innate immunity, which involves pattern recognition receptors (PRRs), Toll-like receptor (TLR) signaling, and cytokine-mediated antiviral responses. To determine the extent to which the skin is activated in response to SARS-CoV-2. If activation does occur, it is crucial to identify the controlling factors.<\/p>\n\n\n\n<p>In this study, we conducted a transcriptomic analysis of skin samples from patients with acute COVID-19 as depicted in&nbsp;<a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#f0001\" id=\"ref-f0001\">Figure 1A<\/a>. Our findings reveal an integrated gene expression profile encompassing antiviral, antioxidant, and innate immune responses in the skin, shedding light on the complex host-pathogen interactions in this compartment.<\/p>\n\n\n\n<p><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#ref-f0001\"><\/a><\/p>\n\n\n\n<figure class=\"wp-block-table\"><table class=\"has-fixed-layout\"><tbody><tr><td><a href=\"https:\/\/www.dovepress.com\/article\/fulltext_file\/516737\/aW1n\/JIR_A_516737_O_F0001g.jpg\"><\/a><\/td><td><strong>Figure 1<\/strong>&nbsp;(<strong>A<\/strong>) Schematic workflow of the study; (<strong>B<\/strong>) Unsupervised z-score heatmap of normalized gene expression counts, based on 150 down and 150 upregulated genes. The control group (CT) is represented in green, and the COVID-19 group (CV) in pink. Red shades indicate z-scores above the row average (higher expression), while blue shades indicate z-scores below the row average (lower expression).<\/td><\/tr><\/tbody><\/table><\/figure>\n\n\n\n<h2 class=\"wp-block-heading\">Materials and Methods<\/h2>\n\n\n\n<h3 class=\"wp-block-heading\">Patients\u2019 Enrollment<\/h3>\n\n\n\n<p>Skin samples used in this study were obtained from patients who succumbed to severe SARS-CoV-2 infection while hospitalized in the Intensive Care Unit of Hospital das Cl\u00ednicas, S\u00e3o Paulo. Post-mortem biopsies were collected in the thigh region using a 5 mm punch, by the Death Verification Service. All patients were diagnosed with COVID-19 through nasopharyngeal detection of SARS-CoV-2 RNA using reverse transcriptase polymerase chain reaction (RT-PCR).<\/p>\n\n\n\n<p>For RNA sequencing, we included skin samples from eight COVID-19 patients (4 men and 4 women), aged 32\u201370 years, who died between April and June 2020 (<a href=\"https:\/\/www.dovepress.com\/article\/supplementary_file\/516737\/ID_516737_supplementary_file.zip\">Supplementary Table 1<\/a>). Additional skin samples from COVID-19 patients (n=10, 3 men and 7 women, aged 26\u201376 years) were used for real-time PCR analysis. At the time, vaccines against SARS-CoV-2 were unavailable. Post-mortem skin biopsies were collected within 20&nbsp;hours of death from the same anatomical region for all patients. Post-mortem skin biopsies were taken up to 20&nbsp;hours after death due to complications arising from a diagnosis of SARS-CoV-2. The exclusion criteria were cancer, HIV-1 infection, and others skin diseases.<\/p>\n\n\n\n<p>As controls, skin biopsies were obtained from healthy donors (n=5) prior to the COVID-19 pandemic. These control samples were collected from Caucasian individuals of both sexes, aged 35\u201370 years. The autopsy was carried out with the authorization of the relatives, who signed an informed consent form approved by the National Research Ethics Committee (CONEP). The study was approved by the ethics committee of the HC-FMUSP (CAAE 30364720.0.0000.0068) and conducted following the Declaration of Helsinki.<\/p>\n\n\n\n<h3 class=\"wp-block-heading\">Transcriptomic Analysis<\/h3>\n\n\n\n<p>Total RNA was extracted from autopsy-derived skin specimens using the RNeasy Plus Mini Kit (Qiagen, Hilden, Germany), following the manufacturer\u2019s instructions. Tissue homogenization was performed with a TissueRuptor (Qiagen). RNA integrity was assessed using the TapeStation system (Agilent, Santa Clara, CA, USA).<\/p>\n\n\n\n<p>Transcriptomic profiling was carried out using the Illumina TruSeq Stranded mRNA Library Prep Kit (Illumina, San Diego, CA, USA) at the FMUSP Large-Scale Sequencing Laboratory (SELA). Sequencing adapters were removed, and quality control of raw reads was performed using FastQC.<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0011\">11<\/a><\/sup>&nbsp;Reads were aligned to the reference human genome (GRCh38.p13, version 107) using the SUBREAD package<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0012\">12<\/a><\/sup>&nbsp;and quantified with FEATURECOUNTS.<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0013\">13<\/a><\/sup><\/p>\n\n\n\n<p>Differentially expressed genes (DEGs) were identified using the DESeq2 package from R\/Bioconductor,<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0014\">14<\/a><\/sup>&nbsp;with statistical significance determined at an adjusted p-value \u2264 0.05 (Benjamini-Hochberg correction). Functional enrichment analyses were conducted with the FGSEA package<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0015\">15<\/a><\/sup>&nbsp;using KEGG<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0016\">16<\/a><\/sup>&nbsp;and GO Biological Process<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0017\">17<\/a><\/sup>&nbsp;databases. Protein-protein interaction networks were constructed using Cytoscape<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0018\">18<\/a><\/sup>&nbsp;with data from the InnateDB database.<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0019\">19<\/a><\/sup><\/p>\n\n\n\n<h3 class=\"wp-block-heading\">Gene Expression by Real-Time PCR<\/h3>\n\n\n\n<p>Total RNA was quantified using a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific). Complementary DNA (cDNA) was synthesized using the iScript\u2122 Reverse Transcription Kit (Bio-Rad, Hercules, CA, USA). Real-time PCR was performed with SYBR Green (Applied Biosystems, Waltham, MA, USA) and primers specific to the target genes (primer sequences listed in&nbsp;<a href=\"https:\/\/www.dovepress.com\/article\/supplementary_file\/516737\/ID_516737_supplementary_file.zip\">Supplementary Table 2<\/a>). Amplification was carried out on a 7500 Real-Time PCR System (Applied Biosystems), and data analysis was conducted using 7500 software v2.0.6 (Applied Biosystems) based on the \u0394\u0394Ct method.<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0020\">20<\/a><\/sup><\/p>\n\n\n\n<h3 class=\"wp-block-heading\">Statistical Analysis<\/h3>\n\n\n\n<p>Statistical comparisons between groups were performed using the Mann\u2013Whitney test, with significance set at p \u2264 0.05. Detection of outlier samples, normalization, statistical tests and identification of DEGs were carried out using the package R DESEQ2 from the Bioconductor repository. The p-values were corrected at using the Benjamini-Hochberg method. Genes were considered differentially expressed for p-adj values \u2264 0.05 and absolute fold-change values, FC \u2265 1.5.<\/p>\n\n\n\n<h2 class=\"wp-block-heading\">Results<\/h2>\n\n\n\n<h3 class=\"wp-block-heading\">Characteristics of the COVID-19 Patients and Healthy Controls<\/h3>\n\n\n\n<p>A group of 8 patients who died due to COVID-19 complications (CV) from April to June 2020, composed of 4 men and 4 women, aged between 32 to 76 years (mean = 55,25), was evaluated (<a href=\"https:\/\/www.dovepress.com\/article\/supplementary_file\/516737\/ID_516737_supplementary_file.zip\">Supplementary Table 1<\/a>). All cases configure severity and death, and comorbidities were identified in almost all individuals. As a control group, we used skin samples from 5&nbsp;healthy donors (HD), composed of 2 men and 3 women, aged between 29 and 70 years (mean = 62). HD skin samples were collected in 2019, before the COVID-19 pandemic. No comorbidities or infectious diseases were present in this group.<\/p>\n\n\n\n<h3 class=\"wp-block-heading\">Transcriptomic Profile of Skin Samples<\/h3>\n\n\n\n<p>Skin samples were processed for RNA extraction, followed by library construction and sequencing. Computational analyses were performed to identify differentially expressed genes (DEGs), enriched molecular pathways, and gene interaction networks. Gene expression of the most important DEGs were validated using real-time PCR (<a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#f0001\" id=\"ref-f0001\">Figure 1A<\/a>).<\/p>\n\n\n\n<p>Unsupervised z-score heatmap of normalized counts revealed a distinct transcriptomic profile between COVID-19 (CV) and healthy control (CT) skin samples (<a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#f0001\" id=\"ref-f0001\">Figure 1B<\/a>). This showed 150 upregulated and 150 downregulated genes. The more important genes related to innate immunity were labelled. A total list of genes can be found in the&nbsp;<a href=\"https:\/\/www.dovepress.com\/article\/supplementary_file\/516737\/ID_516737_supplementary_file.zip\">Supplementary Table 3<\/a>. A volcano plot demonstrated with total genes, that 1188 upregulated and 803 downregulated genes in the CV group compared to the CT group (<a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#f0002\" id=\"ref-f0002\">Figure 2<\/a>). Principal component analysis (PCA) showed that the great majority of the samples grouped closely within their respective cohorts, except for one perhaps mild CV sample (<a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#f0002\" id=\"ref-f0002\">Figure 2<\/a>).<\/p>\n\n\n\n<p><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#ref-f0002\"><\/a><\/p>\n\n\n\n<figure class=\"wp-block-table\"><table class=\"has-fixed-layout\"><tbody><tr><td><a href=\"https:\/\/www.dovepress.com\/article\/fulltext_file\/516737\/aW1n\/JIR_A_516737_O_F0002g.jpg\"><\/a><\/td><td><strong>Figure 2<\/strong>&nbsp;Principal component analysis (PCA) showing the explained variance projected onto the two largest orthogonal principal components; Volcano plot depicting the -log10(adjusted p-value) as a function of log2 fold change (log2FC). Upregulated genes are shown as red circles, downregulated genes as blue circles, and non-significant genes as grey circles. Vertical dashed lines denote fold change thresholds (absolute log<sub>2<\/sub>FC &gt; 1.5), while the horizontal dashed line represents the adjusted p-value significance cutoff (p \u2264 0.05).<\/td><\/tr><\/tbody><\/table><\/figure>\n\n\n\n<p>Enrichment molecular pathway analysis revealed upregulation of 25 pathways associated with inflammatory responses in the CV group compared to control group. These pathways included \u201cregulation of inflammatory response\u201d, \u201cresponse to molecules of bacterial origin\u201d, \u201cleukocyte chemotaxis\u201d, \u201csignaling by interleukins\u201d, \u201cresponse to oxidative stress\u201d, \u201cmatrix organization\u201d, and the \u201cJAK-STAT signaling pathway\u201d, all of which are linked to innate immunity and inflammation (<a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#f0003\" id=\"ref-f0003\">Figure 3A<\/a>). Moreover, pathways analyzed included response to oxidative stress as \u201creactive oxygen species metabolic response\u201d and \u201cresponse to oxidative stress\u201d. This finding was in line with downregulation of pathways related to controlling of inflammatory process, such \u201cpositive regulation of signaling\u201d, \u201cregulation of cell differentiation\u201d, \u201cnegative regulation of response to stimulus\u201d, \u201chomeostatic process\u201d (<a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#f0003\" id=\"ref-f0003\">Figure 3B<\/a>).<\/p>\n\n\n\n<p><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#ref-f0003\"><\/a><\/p>\n\n\n\n<figure class=\"wp-block-table\"><table class=\"has-fixed-layout\"><tbody><tr><td><a href=\"https:\/\/www.dovepress.com\/article\/fulltext_file\/516737\/aW1n\/JIR_A_516737_O_F0003g.jpg\"><\/a><\/td><td><strong>Figure 3<\/strong>&nbsp;(<strong>A<\/strong>) Enrichment analysis of the MSIG molecular pathways database. The normalized enrichment score (NES) is displayed at the center of each circle, with shades of red indicating upregulated pathways. The circle size represents the number of genes in each pathway. (<strong>B<\/strong>) Over-representation analysis obtained with the downregulated genes against the MSIG molecular pathways database. The percentage of overlapping genes is shown at the center of each circle, with colors reflecting the degree of overlap. The circle size corresponds to the number of overlapping genes.<\/td><\/tr><\/tbody><\/table><\/figure>\n\n\n\n<p>Over-representation analysis revealed downregulated molecular pathways in CV samples included those related to the skin matrix, such as \u201cmatrisome\u201d, \u201ccollagen-containing extracellular matrix\u201d, \u201cepidermis development\u201d, \u201cregulation of cell differentiation\u201d, and \u201cconnective tissue development\u201d (<a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#f0003\" id=\"ref-f0003\">Figure 3B<\/a>).<\/p>\n\n\n\n<p>The gene interaction network revealed several upregulated DEGs (pink nodes) associated with antiviral responses, including interferon-stimulated genes (ISGs) such as IFIT1, IFIT2, IFIT3, and interferon-regulatory factors IRF3 and IRF7. Additional antiviral genes included members of the S100 protein family, such as S100P, S100A7, S100A8, and S100A9. Cytokine and cytokine-related factors were also prominent, including STAT2, STAT3, SOCS3, NF-\u03baB2, RELB, NFKBIA, MyD88, IL-6, CCL20, and CXCL2 (<a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#f0004\" id=\"ref-f0004\">Figure 4<\/a>). Some factors were downregulated (blue nodes), such as TNF, a multifunctional proinflammatory cytokine, and MAVS (mitochondrial antiviral signaling protein), which regulates IFN-\u03b2 expression and contributes to antiviral innate immunity (<a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#f0004\" id=\"ref-f0004\">Figure 4<\/a>).&nbsp;<a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#f0005\" id=\"ref-f0005\">Figure 5<\/a>&nbsp;shows the selected DEGs for validation in skin samples, which showed upregulation of S100P, GPX3, SOD2, SOCS3, SERPINE1, and CXCL2, alongside downregulation of TNF-\u03b1 in the skin of COVID-19 patients compared to the CT group.<\/p>\n\n\n\n<p><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#ref-f0004\"><\/a><\/p>\n\n\n\n<figure class=\"wp-block-table\"><table class=\"has-fixed-layout\"><tbody><tr><td><a href=\"https:\/\/www.dovepress.com\/article\/fulltext_file\/516737\/aW1n\/JIR_A_516737_O_F0004g.jpg\"><\/a><\/td><td><strong>Figure 4<\/strong>&nbsp;Network of Genes Involved in COVID-19 Skin Infection Pathways. The network visualization illustrates genes associated with COVID-19-related pathways in the skin. Each circle represents a gene, with the color intensity indicating the log<sub>2<\/sub>FC as shades of red represent upregulated genes (FC \u2265 1.5), and shades of blue represent downregulated genes (FC \u2264 \u22121.5). Only genes with an adjusted p-value &lt; 0.05 were included. Genes directly involved in COVID-19-specific pathways are highlighted with magenta borders. The size of each circle corresponds to the number of connections (degree) the gene has within the network, indicating its relative importance.<\/td><\/tr><\/tbody><\/table><\/figure>\n\n\n\n<p><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#ref-f0005\"><\/a><\/p>\n\n\n\n<figure class=\"wp-block-table\"><table class=\"has-fixed-layout\"><tbody><tr><td><a href=\"https:\/\/www.dovepress.com\/article\/fulltext_file\/516737\/aW1n\/JIR_A_516737_O_F0005g.jpg\"><\/a><\/td><td><strong>Figure 5<\/strong>&nbsp;Selecting innate immune factors to validated genes in the skin of COVID-19 patients via qRT-PCR. Gene expression levels of S100P, S100A8, S100A9, MMP14, CXCL2, CCL20, TNF-\u03b1, SERPINE1, GPX3, SOD2, SOCS3, and ADAMTS1 were assessed in skin autopsies. Data are presented as medians. Statistical significance was determined using the Mann\u2013Whitney&nbsp;<em>U<\/em>-test, with significance levels indicated as p &lt; 0.05 (*), p &lt; 0.01 (**), and p &lt; 0.001 (***). Green dots represent healthy controls (CT), while pink dots represent COVID-19 patients (CV).<\/td><\/tr><\/tbody><\/table><\/figure>\n\n\n\n<p>In addition, deconvolution analysis of immune cell types revealed an increased presence of macrophages and reduced dendritic cells in the skin tissue of COVID-19 patients, likely reflecting their migratory status (<a href=\"https:\/\/www.dovepress.com\/article\/supplementary_file\/516737\/ID_516737_supplementary_file.zip\">Supplementary Figure 1<\/a>). Eosinophils were only present inside the dermal perivascular vessel, not in the tissue infiltration. There is mild lymphohistiocytic infiltration in skin tissue (see&nbsp;<a href=\"https:\/\/www.dovepress.com\/article\/supplementary_file\/516737\/ID_516737_supplementary_file.zip\">Supplementary Figure 2<\/a>).<\/p>\n\n\n\n<h2 class=\"wp-block-heading\">Discussion<\/h2>\n\n\n\n<p>The gene expression profile of the skin in COVID-19 patients revealed significant upregulation of innate immune factors including antiviral components, antioxidant defenses, and inflammatory markers. Although the skin is not a primary target organ for SARS-CoV-2, it engages protective mechanisms that may prevent severe cutaneous lesions. Network analysis highlighted the involvement of innate antiviral immune factors such as Mx1, APOBEC3G, IRF3, IRF7, IRF9, DDX58, IFIT3, and IFIT1. These genes are associated with antiviral responses triggered by interferon signaling, indicating a robust cutaneous antiviral response. Additionally, interferon-induced transmembrane (IFITM) proteins, known to be upregulated in multiple tissues such as the lung, gut, heart, and brain,<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0021\">21<\/a><\/sup>&nbsp;were also expressed in the skin, further underscoring its protective role. Despite the absence of detectable viral load in most cases (present in only two samples), all skin samples exhibited upregulation of antiviral genes.<\/p>\n\n\n\n<p>In terms of innate immunity, the differentially expressed genes found in COVID-19 skin samples included members of the S100 protein family, such as S100A8, S100A9, and S100A7. These proteins, particularly S100A8\/A9, are known alarmins that mediate host pro-inflammatory responses and have been identified as biomarkers of COVID-19 severity.<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0022\">22<\/a><\/sup>&nbsp;Consistent with findings in nasal swabs, the upregulation of S100 genes was associated with disease severity.<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0023\">23<\/a><\/sup>&nbsp;In the skin, S100A8\/A9 acts as a ligand for TLR4 and RAGE, promoting NF-\u03baB signaling and inducing pro-inflammatory cytokines such as IL-6, which was notably upregulated in COVID-19 skin samples. This contrasts with the downregulation of TNF-\u03b1, a cytokine often implicated in the cytokine storm of severe COVID-19 cases.<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0024\">24<\/a><\/sup>&nbsp;The reduced expression of TNF-\u03b1 in the skin may represent a counter-regulatory mechanism to limit excessive inflammation, possibly mediated by the upregulation of SOCS3, a known negative feedback inhibitor of cytokine signaling, and also includes NFKBIA which encodes a protein known as I\u03baB\u03b1 (NF-\u03baB alpha inhibitor). These findings align with prior evidence of cross-regulation between TNF-\u03b1 and type I IFNs triggering,<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0025\">25<\/a><\/sup>&nbsp;where TNF-\u03b1 has been shown to suppress type I IFN production and pDC activity, while IFN-\u03b2 can inhibit TNF-\u03b1 secretion. Given that TNF-\u03b1 may suppress an important antiviral factor such as type I IFN, the cutaneous reduction in TNF-\u03b1 is timely.<\/p>\n\n\n\n<p>Over-representation analysis obtained with the downregulated genes related with impaired development of the extracellular epidermis and connective tissue and influence the dynamic repair of skin tissue in skin of severe cases of COVID-19. However, the upregulation of CXCL2 plays a role in the stages of haemostasis, inflammation, proliferation and remodeling in the skin. This finding suggests that, rather than acting as an eosinophil\/neutrophil chemoattractant, CXCL2 acts more in skin remodeling. This is particularly relevant given that skin infiltration in samples from patients with SARS-CoV-2 was primarily composed of lymphocytes and histiocytes (macrophages).<\/p>\n\n\n\n<p>The number of downregulated genes related to innate immunity in the skin was much smaller than the number of upregulated genes. Such as, TNF-\u03b1, mitochondrial antiviral signaling protein (MAVS) and WNT2. MAVS plays a key role in the antiviral response. MAVS is activated in response to viral infection by RIG-I-like receptors that recognize viral RNA, mediating the activation of NF-\u03baB and interferon-regulatory factors, as well as the induction of interferons, in response to viral infection. It is required for interferon induction in dendritic cells (DCs).<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0026\">26<\/a><\/sup>&nbsp;This observation aligns with previous reports that SARS-CoV-2 nucleocapsid protein inhibits type I IFN responses by interacting with MAVS, thereby suppressing host antiviral defenses. However, it remains unclear whether the downregulation of MAVS in the skin reflects a specific viral evasion strategy. The Wnt\/\u03b2-catenin signaling pathway inhibits the replication of SARS-CoV-2 and other pathogenic RNA viruses in vitro.<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0027\">27<\/a><\/sup>&nbsp;Furthermore, the Wnt\/\u03b2-catenin signaling pathway strongly suppresses peroxisome biogenesis. These are intracellular organelles containing enzymes for beta-oxidation, and they are depleted from cells and tissues during SARS-CoV-2 infection. This depletion may be expected to impair the IFN response.<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0028\">28<\/a><\/sup><\/p>\n\n\n\n<p>Despite the downregulation of MAVS, our findings showed that some first-line defense antioxidants, such as superoxide dismutase (SOD) and glutathione peroxidase (GPX), are upregulated in the skin. These genes play critical roles in mitigating oxidative stress, with SOD2 converting reactive oxygen species (ROS) to hydrogen peroxide and GPX3 metabolizing hydrogen peroxide and lipid hydroperoxides toxins.<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0029\">29<\/a><\/sup>&nbsp;This robust antioxidant response highlights the importance of oxidative stress mitigation in SARS-CoV-2 pathogenesis and suggests potential therapeutic targets.<\/p>\n\n\n\n<p>Another notable upregulated DEG was SERPINE1, which encodes plasminogen activator inhibitor 1 (PAI-1). Upregulation of PAI-1 may counteract the thrombosis frequently observed in COVID-19 by modulating fibrinolysis. SERPINE 1 is required to stimulate keratinocyte migration during the repair of cutaneous injuries.<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0030\">30<\/a><\/sup>&nbsp;A polymorphism in the human SERPINE1 gene increases influenza A virus susceptibility in vitro.<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0031\">31<\/a><\/sup>&nbsp;Furthermore, the ability of SERPINE1 to inhibit TMPRSS2&nbsp;has been demonstrated, resulting in a reduction in SARS-CoV-2 entry and infection.<sup><a href=\"https:\/\/www.dovepress.com\/heightened-innate-immune-and-inflammatory-gene-expression-in-the-skin--peer-reviewed-fulltext-article-JIR#cit0032\">32<\/a><\/sup>&nbsp;We observed that protease inhibitors such as SERPINE 1, was observed in heatmap as DEGs, in the expressions in skin from COVID-19, have shown to be valuable component of the antiviral innate immunity.<\/p>\n\n\n\n<p>Analysis of immune cell composition revealed an increased presence of macrophages in the skin of COVID-19 patients, along with a reduction in DCs. The decreased number of DCs in the skin may reflect their migration to the draining lymph nodes, where antigen-specific adaptive immunity occurs. Conversely, eosinophils were only present in the perivascular dermal vessels and not in the tissue infiltrate, which was found to be mild lymphohistiocytic infiltration of the skin.<\/p>\n\n\n\n<h2 class=\"wp-block-heading\">Conclusion<\/h2>\n\n\n\n<p>SARS-CoV-2 triggers a significant innate immune response in the skin, characterised by the upregulation of antiviral, inflammatory and antioxidant pathways. Several genes related to innate immunity are upregulated, including members of the S100 family, IFITs, ISGs, pro-inflammatory cytokines, and chemokines, indicating a robust cutaneous antiviral response but also with inflammatory potential. In this line, a counter-regulatory response is induced, by genes such as SOCS3, which controls cytokine levels, and NFKBIA, which controls NF-\u03baB activity. These genes may, partly, control the inflammatory status in the skin. The severe COVID-19 patient showed intense innate immunity enriched with antiviral factors. This could explain why they did not exhibit significant cutaneous lesions. Further investigation into the interaction of immune pathways in the skin could lead to the development of new therapeutic strategies for controlling the inflammation associated with the disease.<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Authors\u00a0Fernandes IG,\u00a0Ferreira FM,\u00a0De Sousa ESA\u00a0,\u00a0Pet.al., DOI\u00a0https:\/\/doi.org\/10.2147\/JIR.S516737 Correspondence: Maria Notomi Sato, Laboratory of Medical Investigation in Dermatology and Immunodeficiencies (LIM-56), Department of Dermatology, Faculdade de Medicina, Universidade de Sao Paulo, Sao [&hellip;]<\/p>\n","protected":false},"author":2,"featured_media":14164,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[797,253,533],"tags":[],"class_list":["post-13788","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-epidermis","category-inflammation-2","category-skin-dermatology"],"_links":{"self":[{"href":"https:\/\/cov19longhaulfoundation.org\/index.php?rest_route=\/wp\/v2\/posts\/13788","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/cov19longhaulfoundation.org\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/cov19longhaulfoundation.org\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/cov19longhaulfoundation.org\/index.php?rest_route=\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/cov19longhaulfoundation.org\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=13788"}],"version-history":[{"count":3,"href":"https:\/\/cov19longhaulfoundation.org\/index.php?rest_route=\/wp\/v2\/posts\/13788\/revisions"}],"predecessor-version":[{"id":13804,"href":"https:\/\/cov19longhaulfoundation.org\/index.php?rest_route=\/wp\/v2\/posts\/13788\/revisions\/13804"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/cov19longhaulfoundation.org\/index.php?rest_route=\/wp\/v2\/media\/14164"}],"wp:attachment":[{"href":"https:\/\/cov19longhaulfoundation.org\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=13788"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/cov19longhaulfoundation.org\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=13788"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/cov19longhaulfoundation.org\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=13788"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}