Abstract
Introduction:
Neurological complications are increasingly recognized as a significant consequence of COVID-19; however, time-resolved, brain-specific characterization of transcriptional alterations underlying SARS-CoV-2–associated neuroinflammation and neuronal injury remain limited. We hypothesized that brain transcriptional responses evolve dynamically during acute SARS-CoV-2 infection, resulting in temporal transcriptional programs.
Methods:
KRT18-hACE2 transgenic mice were intranasally inoculated with SARS-CoV-2. Brain was harvested at 4 and 6 days post-infection (dpi) for analyses.
Results:
Immunohistochemical analyses confirmed a broad spectrum of viral neurotropism and gliotropism, accompanied by an increased apoptotic burden, particularly in cortical neurons (ClCas3/SATB2+). Robust activation of myeloid cells (Iba1+/CD68+) provided evidence of neuroinflammation. Cytokine/chemokine profiling demonstrated pronounced upregulation of inflammatory mediators (CXCL10, IL-12p40, CCL12), alongside reduced CX3CL1, suggesting impaired neuron–microglia communication. Whole-transcriptome and gene ontology analyses uncovered stage-dependent molecular programs, with early alterations at 4 dpi enriched in protein ubiquitination, vesicle trafficking, and synaptic processes, followed by intensified innate immune activation and engagement of chromosomal maintenance pathways at 6 dpi. In parallel, pronounced suppression of mitochondrial function at 6 dpi, pointing to energy exhaustion and transcriptional-translational discordance, as supported by digital PCR and a substantial reduction in COXIV protein levels.
Discussion:
These findings provide a time-resolved molecular landscape of SARS-CoV-2–induced neuroinflammation and metabolic stress, highlighting CNS vulnerability during severe infection and suggesting pathways potentially relevant to COVID-19-associated sequelae.
1 Introduction
COVID-19, caused by SARS-CoV-2, has become a global health crisis, placing unprecedented strain on healthcare systems and economies worldwide (1, 2). Although primarily a respiratory disease, accumulating evidence indicates that SARS-CoV-2 can affect central nervous system (CNS) (3, 4), where it is capable of triggering neuroinflammation (5, 6) and contributing to persistent neurological and psychiatric sequelae (7). Clinical manifestations span a broad spectrum, ranging from anosmia and ageusia to severe neurological conditions such as encephalitis, seizures, and cerebrovascular events (8, 9). These understandings mark SARS-CoV-2 as a significant, though relatively indirect, neurotropic virus, underscored by its capacity to perturb CNS homeostasis. Importantly, converging clinical, pathological, and experimental data suggest that direct viral neurotropism is not the principal driver of most COVID-19-associated neurological sequelae; rather, indirect mechanisms, including cytokine dysregulation, endothelial and blood-brain barrier disruption, microthrombosis, and hypoxia, appear to play dominant pathogenic roles (10–13), with direct CNS infection contributing in a subset of cases, potentially acting as a triggering factor (14). Despite this growing body of evidence, the molecular and transcriptional programs underlying virus-driven neuroinflammation remain incompletely characterized, particularly in a time-resolved manner.
hACE2 transgenic mouse models, which permit SARS-CoV-2 entry via the ACE2 receptor, are widely used to investigate CNS pathogenesis and have provided evidence for viral neuroinvasion, neuroinflammation, and neuronal injury (15–19). SARS-CoV-2 infection progresses through distinct phases, each characterized by dominant viral and host processes, including induction of inflammatory gene programs (19), elevated cytokine and chemokine expression in the brain (17), and progressive manipulation of host cellular machinery (20, 21). Notably, epidemiological and biomarker studies have indicated that SARS−CoV−2 infection is associated with an increased risk of subsequent neurodegenerative diagnoses and with sustained markers of neuronal and astroglial injury in long−COVID cohorts (22, 23), suggesting a potential link between neuro-COVID and neurodegenerative-like brain changes that extends beyond viral burden alone, although causality and the precise mechanistic cascade remain to be fully elucidated. Temporal interrogation of inflammatory and transcriptional dynamics is essential for understanding the evolution of SARS-CoV-2 associated neuropathology. In this context, previous comparative analyses on beta over original variant has reported time-dependent alterations in systemic and CNS immune responses (24). However, a time-resolved characterization of brain-specific transcriptional programs at key disease stages remains limited, and such insights are critical for understanding SARS-CoV-2 neuropathogenesis and for identifying therapeutic targets aimed at mitigating both acute and long-term neurological complications.
Here, we aimed to characterize the neuropathological and transcriptional responses in the brains of KRT18-hACE2 transgenic mice following intranasal SARS-CoV-2 infection, focusing on two critical stages, 4- and 6-days post-infection (dpi), which correspond to peak viral replication and heightened immune activation preceding lethal diseases in this model. By integrating cytokine/chemokine profiling with whole-transcriptome profiling, gene ontology analysis and absolute quantification of selected targets by digital PCR, we delineate time-dependent neuroinflammatory and metabolic remodeling. Our analyses reveal stage-specific transcriptional programs and RNA–protein discordance in key inflammatory and mitochondrial pathways, highlighting complex regulatory mechanisms underlying SARS-CoV-2–induced neuroinflammation. Collectively, these findings provide a descriptive, time-resolved molecular landscape that informs mechanisms of acute neuro-COVID and lays the groundwork for the development of stage-specific therapeutic strategies aimed at mitigating both acute and long-term neurological complications.
2 Materials and methods
2.1 Animals and viral infection
All work involving live SARS-CoV-2 was conducted in a Biosafety Level 3 (BSL-3) containment facility, in accordance with approved safety protocols and with the use of appropriate personal protective equipment. Animal procedures were approved by the Institutional Animal Care and Use Committee (IACUC) of the Korea Research Institute of Chemical Technology (permit no. 2021-8A-04-01, protocol no. 8A-M6). Male K18-hACE2 mice (no. 034860) [B6.Cg-Tg(K18-ACE2)2Prlmn/J] were purchased from the Jackson Laboratory. Mice were housed in groups under standard laboratory conditions with free access to food and water and were randomly divided into groups before infection.
SARS-CoV-2 wild-type strain (BetaCoV/Korea/KCDC03/2020) (GISAID accession number EPI_ISL_407193) was provided by the National Culture Collection for Pathogens (25) (NCCP; numbers 43326 and 43382) of the Korea Disease Control and Prevention Agency. For SARS-CoV-2 infection, mice were anesthetized with isoflurane and injected intranasally with 50 μL of saline containing 2 × 103 PFU of SARS-CoV-2 wild-type strain (25 μL per nostril). Control mice were mock-infected with vehicles. This inoculum represents a sublethal respiratory dose but is the minimum required to induce reproducible and consistent neuroinvasion and CNS pathology in K18-hACE2 mice. Weight loss, individual clinical scores and mortality were monitored daily after infection and reported in our previous study (24).
2.2 Tissue collection
Mice were sacrificed at 4- and 6- day post-infection (dpi). Uninfected and SARS-CoV-2-infected mice were anesthetized with isoflurane. Whole brain was perfused with 10 mL of phosphate-buffered saline (PBS) to remove blood. For immunofluorescence, brains were fixed in 10% neutral buffered formalin at room temperature for 2 to 5 days. For microarray and cytokine/chemokine arrays, brains were homogenized in cold PBS using a FastPrep-24 homogenizer (MP Biomedicals) for three cycles of 15 s on and 20 s off. Tissue debris was removed by brief centrifugation.
2.3 Immunohistochemistry
Fixed brain samples were embedded in paraffin blocks, and 5-μm-thick sagittal sections of hemisphere were generated using a rotary microtome (Roundfin RD-315, Shenyang, China), mounted on glass slides, and allowed to dry. Prior to staining, sections were deparaffinized and rehydrated by sequential immersion in three changes of xylene, followed by xylene:100% ethanol (1:1, v/v), 100%, 95%, and 80% ethanol (10 min each), and finally distilled water. Antigen retrieval was performed by boiling the sections in citrate buffer (10 mM sodium citrate, 0.05% Tween-20, pH 6.0) for 20 min. After cooling to room temperature, samples were blocked with 5% bovine serum albumin (BSA) for 30 min and incubated with primary antibodies overnight at 4°C. The following day, sections were washed three times with PBS (10 min each), incubated with fluorescently labeled secondary antibodies diluted in 1% BSA for 1 h at RT, and counterstained with 4′,6-diamidino-2-phenylindole (DAPI, 100ng/mL). Finally, slides were washed twice (10 min each), mounted, and imaged using K1 Fluo confocal microscopy (Nanoscope Systems, Inc., Korea). Images were acquired with a 40× objective lens and a resolution of 2048 × 2048 pixels. For each marker, acquisition parameters were optimized using control sections and kept constant for 6dpi group.
Three mice per group (n = 3) were analyzed, with three brain sections per mouse examined, focusing on cortical regions. Quantification was performed manually in a blinded manner across a defined number of high-power fields per section. Marker-positive and double-positive cells were identified based on fluorescence intensity above background and cellular morphology, and results were expressed as the percentage of positive cells relative to DAPI-positive nuclei. Representative images were selected to reflect the quantified results.
Primary antibodies used in this study included anti-NP (Cell Signaling Technology, 33717), anti-SATB2 (Abcam, ab34735), anti-MAP2 (Abcam, ab183830), anti-GFAP (Abcam, ab7260), anti-O4 (R&D system, MAB1326), anti-Iba1 (Cell Signaling Technology, 58970), anti-CD68 (Cell Signaling Technology, 97778), anti-Cleaved Caspase 3 (Cell Signaling Technology, 9661), all used at a dilution ratio of 1:500.
Secondary antibodies used in this study included goat anti-rabbit Alexa Fluor 488-conjugated IgG (Abcam, ab1500770), goat anti-rabbit Alexa Fluor 555-conjugated IgG (Abcam, ab1500780), goat anti-mouse Alexa Fluor 488-conjugated IgG (Abcam, ab150113), goat anti-mouse Alexa Fluor 555-conjugated IgG (Abcam, ab150114), all used at a dilution of 1:300.
2.4 TUNEL assay
DNA fragmentation was detected on rehydrated paraffin-embedded sections using a one-step TUNEL in situ apoptosis detection kit (E-CK-A325; Elabscience), according to the manufacturer’s instructions. The TUNEL assay was performed on the same set of tissue sections used for immunohistochemical analyses.
2.5 Proteome profile analysis
Five mice per group (n=5) were analyzed. Supernatant of brain homogenates were treated with 1% Triton X-100 to inactivate SARS-CoV-2 (15). Initial screening of 111 inflammatory proteins was performed using Proteome Profiler Mouse XL Cytokine Array (ARY028, R&D System), according to the manufacturer’s instructions. Chemiluminescent signals were detected and quantified using Quick Spots Tool (Western Vision software, Version 22.0.1b) and presented as relative fold-change in signal intensity compared with uninfected control.
2.6 Total RNA extraction & cDNA synthesis
The control group consisted of five mice (n=5), the S4d and S6d groups each comprised six mice (n=6). Total RNA was isolated from homogenized tissues using Maxwell RSC SimplyRNA Tissue Kit (Promega) according to the manufacturer’s instructions. RNA concentration and purity were assessed using a NanoDrop spectrophotometer (Nano400A; Allsheng). Samples that passed quality control (A260/A280 > 1.6) were used for cDNA synthesis using iScript Reverse Transcription Supermix for RT–qPCR (Bio-Rad) and a Mastercycler (Eppendorf). The resulting cDNA samples were used for whole-transcriptome analysis and digital PCR.
2.7 Whole-transcriptome analysis
Whole-transcriptome analysis was performed as in prior study using the Affymetrix GeneChip Mouse Gene 2.0 ST array (no. 902463) and GeneChip WT Pico kit (no. 902623) (Thermo Fisher Scientific) according to the manufacturer’s instructions. The raw data were summarized and normalized using the Robust Multi-Average method in Affymetrix Power Tools, followed by differential expression analysis. Statistical significance was determined using an independent t test and fold change. The null hypothesis was that no difference exists between the groups. The false-discovery rate was controlled by adjusting the P value using the Benjamini-Hochberg algorithm.
2.8 Differential expression analysis
We determined the statistical significance of the DEGs by performing pairwise contrasts between the 4dpi or 6dpi groups compared to the uninfected group. Genes with p-adj <0.05 (t-test) and ∣log2FC| >1 were defined as differentially expressed genes. The resulting upregulated and downregulated gene lists for each experimental comparison were then used for Gene Ontology (GO) enrichment analysis via the ShinyGo web tool (http://geneontology.org). All data analyses and visualization were performed in R 3.3.2 (www.r-project.org).
2.9 Digital PCR
Absolute quantification of copy number of each target gene in each sample was performed using Absolute Q™ DNA Digital PCR Master Mix (ABI, Thermo Fisher Scientific) and QuantStudio Absolute Q Digital PCR System (ABI) following manufacture’s procedure. The primers used for dPCR were provided in Supplementary Table 1.
2.10 Statistical analysis
All the statistical analyses were performed via GraphPad Prism 9 software. The data were tested for normal distribution and within-group variance prior to statistical analysis. The data obtained from at least three independent experiments were averaged and are presented as the means ± SDs as stated in the figure legends. In all cases, a p value less than 0.05 was considered significant: *P < 0.05, **P < 0.01, ***P < 0.001.
3 Results
3.1 Neuroinvasion and cellular tropism of SARS-CoV-2 in the brain of KRT18-hACE2 transgenic mice
Previous studies, including our own, have shown that K18-hACE2 mice develop severe clinical disease and mortality around 5–6 days after SARS-CoV-2 infection, with peak viral burden and widespread brain distribution occurring at 6–7 days post-infection (dpi) (26, 27). Therefore, 6dpi was chosen to capture pronounce inflammation and severe pathology status. Accordingly, to specifically interrogate the neuroinvasive phase, we first confirmed SARS-CoV-2 neuroinvasion in the brains of intranasally infected K18-hACE2 mice at 6 dpi (Figure 1A). IHC analysis revealed co-expression of the viral nucleocapsid protein (NP) with MAP2 (Figure 1B) and SATB2 (Supplementary Figure S1A), indicating infection of mature neurons and cortical neurons, respectively. In addition, NP+ GFAP+ astrocytes (Figure 1C) and NP+ O4+ oligodendrocytes (Figure 1D) were detected, demonstrating viral infection of multiple glial populations.
Figure 1

Neuroinvasion and cellular tropism of SARS-CoV-2 in a mouse model. (A) Schematic illustration of experimental design. KRT18-hACE2 transgenic mice were intranasally infected with SARS-CoV-2 and brain tissue was harvested six days post-infection (6 dpi). (B-D) Immunofluorescence staining of brain tissue from uninfected and infected KRT18-hACE2 transgenic mice. Each panel shows co-localization of DAPI (blue, nuclear stain) with SARS-CoV-2 nucleocapsid protein (NP, green) and specific (B) mature neuron MAP2 marker (red), (C) astrocyte GFAP marker (red), (D) oligodendrocyte O4 marker (red).
Together, these results confirm direct neuroinvasion by SARS-CoV-2 at a defined acute disease stage and reveal broad cellular tropism across neuronal and glial lineages, providing a cellular framework for subsequent analyses of transcriptional and functional disruption in the infected brain.
3.2 SARS-CoV-2–induced neuropathologic alterations in the brain of KRT18-hACE2 transgenic mice
The identified neuro- and glio-tropic effects of SARS-CoV-2 prompted us to investigate neurodegenerative changes in the brain of KRT18-hACE2 mice following infection. To this end, we first assessed neuroinflammation as an early histopathological indicator. IHC analysis of 6 dpi brain revealed a marked increase in Iba1-positive cells (Figure 2A) and elevated CD68 expression (Figure 2B) compared with control, indicating activation of brain-associated myeloid cells, comprising both CNS-resident and infiltrating immune populations, and enhanced phagocytic activity. Consistently, infected brains exhibited an increased proportion of apoptotic cells, as demonstrated by TUNEL staining (Figures 2C, D), accompanied by a substantial increase in cortical neuronal apoptosis, evidenced by ClCas3/SATB2 double-positive cells (Supplementary Figures S1B, C).
Figure 2

Neuropathological changes and inflammatory mediators in K18-hACE mouse brain following SARS-CoV-2 infection. (A, B) Immunofluorescence staining of brain tissue from uninfected and SARS-CoV-2 infected mice showing colocalization of DAPI (blue, nuclear stain) with (A) Iba-1 (green) and (B) CD68 (red). Lower panels show single staining for Iba-1 or CD68. Scale bar: 50 µm. (C) TUNEL assay of brain tissue from uninfected and SARS-CoV-2 infected mice. (D) Quantification of panel C showing the percentage of TUNEL+ cells in each group. Student t-test was performed to determine the significance. **p < 0.01, ***p < 0.001).
Collectively, these results demonstrate that SARS-CoV-2 elicits robust immune responses in the brain and is associated with neurodegeneration, especially upper-layer cortical neurons, underscoring the contribution of inflammatory and apoptotic pathways to the neuropathological consequences of infection.
3.3 Enhanced immune response in SARS-CoV-2 infected brain of KRT18-hACE2 transgenic mice
To further characterize the inflammatory state of the brain at 6dpi, we performed proteomic profiling using a cytokine/chemokine targeting 111 analytes (Figure 3A). The results demonstrated a pronounce pro-inflammatory signatures, marked by the significant upregulation of the chemokines CCL5, CCL12, CXCL9, and CXCL10. In parallel, enhanced expressions of immune-associated proteins, including coagulation factor III, IL-12p40, and the acute-phase proteins Fetuin-A, were observed (Figures 3B, C). In contrast, expression of CX3CL1, a chemokine critical for maintaining the neuronal – microglial axis (28, 29), was reduced. This reciprocal pattern of increased pro-inflammatory mediators and decreased homeostatic signaling molecules suggests a shift toward an inflammatory state of the infected brain at 6 dpi. Furthermore, the elevated RBP4 expression, a key mediator of retinol signaling in the brain (30), plausibly pointing to metabolic disorders.
Figure 3

Cytokine and chemokine expression changes in SARS-CoV-2 infected brains. (A) Cytokine/chemokine antibody array showing protein expression in the brains of uninfected and SARS-CoV-2 infected mice. Each numbered spot corresponds to a specific cytokine or chemokine. (B) Quantification of the relative density of the spots in (A) Bars represent the relative density of each spot, indicating the expression level of the corresponding cytokine or chemokine. One-way ANOVA was performed to calculate significance. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001). (C) Volcano plot of differentially expressed cytokines/chemokines from (B) The axes represent the difference of mean density (fold change, x-axis) and statistical significance (-log10(q-value), y-axis). Labels for specific cytokines/chemokines of interest are shown.
These data, therefore, demonstrate that SARS-CoV-2 infection is associated with an enhanced immune and inflammatory response in the brain at the protein level and suggest the association of metabolic signaling to neuropathological alterations.
3.4 Stage-dependent antiviral and immune transcriptional responses to SARS-CoV-2 infection in the brain of KRT18-hACE2 transgenic mice
Our prior studies (24) and published literatures (26, 27) demonstrated that 6 dpi represents fully developed acute disease phase with maximal inflammation and imminent mortality, whereas 4 dpi corresponds to peak viral replication and early inflammatory responses in K18-hACE2 mice, we asked whether transcriptomic alterations also occur in a time-resolved manner. To examine the temporal dynamics underlying the development of immune and inflammatory responses during infection, RNA microarray analyses were performed at 4 dpi, representing early immune activation, and at 6 dpi, capturing severe pathology and mortality (Figure 4A). Consistent with the proteomic profiling, most immune-related transcripts were more prominently upregulated at 6dpi, indicating a stage-dependent escalation of immune activation (Figure 4B). Specifically, transcripts associated with interleukins, chemokines, TNF signaling, inflammasomes, and interferon pathways were markedly elevated post-infection, particularly at 6 dpi, highlight the progressive immune activation in the brain. Principal component analysis (PCA) further supported this observation, showing distinct stage-specific clustering that separated uninfected from infected groups at both time points (Figure 4C). Volcano plots for diffentially expressed genes (DEGs) analysis relative to the uninfected group further demonstrated a substantially broader range of up- and downregulated transcriptional signatures in 6dpi compared to 4 dpi (Figures 4D, E). Notably, transcripts that were significantly altered at 6 dpi showed no significant changes at 4 dpi, indicating extensive dynamic remodeling of the host responses over the course of infection. Identification of the top 100 DEGs across three groups (Figure 4F) emphasized clear expression differences between 4dpi and 6dpi, further supporting a time-dependent pathological transition in the KRT18-hACE2 mouse brain following SARS-CoV-2 invasion.
Figure 4

Stage-dependent antiviral responses to SARS-CoV-2 infection in the brain of KRT18-hACE2 transgenic mice. (A) Schematic overview of experimental design. (B) Heat map demonstrating RNA level of each gene in each experimental condition (Uninfected, 4dpi, 6dpi). (C) PCA plot showing distinct clustering patterns among RNA microarray profiles for each experimental group. (D, E) Volcano plot displaying the log2FC values of the differentially expressed genes in the (D) 4dpi and (E) 6dpi group versus the control group. (F) Heatmap displaying differential RNA expression profiles across Uninfected, 4 dpi and 6 dpi groups, highlighting distinct transcriptional changes over the course of SARS-CoV-2 infection.
3.5 Early-stage transcriptional remodeling in the brain at 4 dpi
To further dissect the stage-dependent immune response to SARS-CoV-2 infection, we performed gene ontology (GO) enrichment analysis of DEGs at each stage. At 4dpi, upregulated biological processes (BP) terms were predominantly enriched for pathways related to protein biogenesis and turnover, as well as synaptic organization and plasticity (Figure 5A). MA plot analysis supported these enrichments, highlighting the top 10 marker genes within each selected pathway such as proteasome/ubiquitin, endosomal transport (Supplementary Figure S2A). Enrichment of macroautophagy-related terms was also observed. Upregulation was also notable in macroautophagy-related pathways, suggesting potential engagement of autophagy in response to viral infection and organelle stress. Additionally, GO terms related to synaptic vesicle cycling and TOR pathway were enriched among upregulated genes, indicating transcriptional changes affecting pathways involved in neuronal function at 4dpi. These transcription signatures are primarily associated with subcellular domains mediating ubiquitin-proteasome degradation, transcriptional regulation, particularly within neuronal compartments (Figure 5B), with prominent ubiquitination-related enzymatic activities (Figure 5C). Together, these findings indicate early transcriptional changes in pathways associated with protein homeostasis and vesicular processes in infected brain. Importantly, the high distribution within vesicular system together with promoted GTPase associated activities suggests the involvement of intracellular trafficking processes required for SARS-CoV-2 egress (31, 32), where GTPases play a critical role (33).
Figure 5

Comparative gene ontology analysis of brains at 4 dpi relative to uninfected controls. (A-C) Top 15 significantly upregulated terms in each category: (A) biological process (BP), (B) cellular component (CC), (C) molecular function (MF); (D-F) Top 15 significantly downregulated terms in each category: (D) BP, (E) CC, (F) MF. The results illustrate enriched pathways and molecular features associated with SARS-CoV-2-induced transcriptional alterations.
In contrast, downregulated GO terms at 4dpi were enriched for pathways governing cytoarchitecture, including neuronal patterning, regionalization and specialization (Figure 5D). The most strongly suppressed genes were associated with extracellular matrix (ECM) organization, cell aggregation, and morphogenetic programmed cell death (Supplementary Figure S2B). Downregulation of genes controlling local protein synthesis at presynapse sites was also detected, suggesting alterations in transcriptional program linked to synaptic plasticity. These downregulated transcripts were predominantly associated with ECM-related compartments (Figure 5E). Notably, GO terms related to retinoid pathways, a critical antiviral mechanism known to interfer with spike-mediated cell entry (34), were significantly suppressed at 4dpi (Figure 5F).
Collectively, these GO enrichment analyses reveal stage-specific transcriptional remodeling in the brain at 4dpi characterized by coordinated changes in gene sets associated with protein homeostasis, vesicular organization, and neuronal structural programs during the neurotropic phase of SARS-CoV-2 at 4 dpi.
3.6 Late-stage transcriptional remodeling in the brain at 6 dpi
At 6 dpi, transcriptomic signatures were predominantly enriched for innate immune processes in brain (Figure 6A; Supplementary Figure S2C). GO analysis revealed significant enrichment of genes sets associated with endolysosomal and phagocytic compartments, and activation of multiple RNA regulators localized to stress granules, P-bodies, and RNP granules (Figure 6B). These enrichments indicate broad transcriptional representation of pathways typically linked to antiviral and inflammatory responses.
Figure 6

Comparative gene ontology analysis of brains at 6 dpi relative to uninfected controls. (A-C) Top 15 significantly upregulated terms in each category: (A) BP, (B) CC, (C) MF. (D-F) Top 15 significantly downregulated terms in each category: (D) BP, (E) CC, (F) MF. The results reflect enriched pathways and molecular features associated with SARS-CoV-2-induced transcriptional alterations.
In parallel, molecular function GO terms related to ubiquitin-like (UBL) ligases and transferases were significantly enriched (Figure 6C), reflecting transcriptional changes in pathways involved in protein modification and turnover. Additional enrichment was observed for pathways associated with nucleocytoplasmic and nuclear transport, and membrane-associated compartments, suggesting altered expression of genes involved in nuclear substructure such as PML bodies and nuclear specks, which are broadly associated with transcriptional regulation and nuclear organization.
Conversely, GO analysis revealed significant downregulation in multiple processes spanning different steps of cellular respiration at 6dpi (Figure 6D). Pathways involved in purine (ribo)nucleotide and ribose phosphate biosynthesis, tightly linked to ATP synthesis, were among the most suppressed. Consistent with these findings, downregulated genes were enriched in mitochondrial and mitochondria-associated cellular components (Figure 6E), including transcripts corresponding to multiple stages of mitochondrial respiration (Supplementary Figure S1D). Molecular function analysis further highlighted reduced representation of enzymatic activities essential for core metabolic processes, along with altered expression of genes associated with ECM organization (Figure 6F).
Together, these GO enrichment analyses indicate that late-stage SARS-CoV-2 infection is associated with coordinated transcriptional changes characterized by enrichment of immune-related pathways alongside suppression of metabolic and mitochondrial gene programs in the brain.
3.7 Time-resolved remodeling of immune- and mitochondria- associated genes
To further validate the stage-dependent alterations in immune response–associated genes identified in the K18-hACE2 mouse brain following SARS-CoV-2 infection, we measured absolute gene copy number in each sample through digital PCR (dPCR) analysis. We first quantified viral burden at 4 and 6 dpi by targeting the viral NP gene and found that viral RNA was not increased at 6 dpi compared with 4 dpi (Figure 7A), which is consistent with western blot result for SARS-CoV-2 NP protein (Supplementary Figure S3). These results indicate that transcriptional changes at 6dpi are not simply driven by higher viral load.
Figure 7

dPCR quantification of absolute copy numbers of selected genes in brain tissue at 4 and 6 dpi. The left panel displays representative fluorescence amplitude plots showing the distribution of individual droplets. The right panels present quantification of transcript copy concentrations (copies/µL) for each target gene across samples: (A) SARS-CoV-2 and Interleukins, (B) Chemokines, (C) Interferons and TNF, (D) The top five upregulated mitochondria-associated genes identified by RNA microarray; and (E) The top five downregulated mitochondria-associated genes identified by RNA microarray. Statistical significance was determined using Student’s t-test (* p<0.05, **p < 0.01, ***p < 0.001).
dPCR analysis of key interleukin genes, including IL-1β, IL-10, and IL-6, revealed a marked increase in transcript copy numbers at 6 dpi compared with 4 dpi, whereas IL-2 transcripts exhibited no significant change (Figure 7B). In contrast, chemokines involved in monocyte chemoattraction (CCL2) and neutrophil recruitment (CXCL1, CXCL2) were significantly elevated at 4 dpi but showed a substantial decline by 6 dpi. Meanwhile, CCL3 and CXCL10, both broadly recognized as interferon-responsive chemokines, exhibited a sustained and progressive increase over the course of infection. Consistently, interferon transcripts also increased in a time-dependent manner, reaching significantly higher levels at 6 dpi compared with 4 dpi (Figure 7C).
Together, these transcriptional patterns suggest a shift in immune state between 4 and 6 dpi, suggesting attenuation of early innate inflammatory recruitment signals, while interferon-associated and adaptive immune responses are enhanced.
Based on our GO enrichment analyses, impaired energy metabolism and mitochondrial pathways emerged as prominent features of the 6-dpi brain transcriptome. To validate these pathway-level findings, we filtered mitochondria-associated genes and selected the top five up- and down-regulated genes at 6 dpi compared with uninfected controls, ranked by adjusted p-value and log2 fold change, for dPCR analysis. Consistent with RNA microarray results, the top five upregulated mitochondrial signatures displayed time-dependent increases in transcript abundance (Figure 7D). Gene ontology annotation of these genes converged on processes related to mitophagy, mitochondrial fission, and mitochondrial organization, supporting progressive mitochondrial alterations associated with time-resolved infection.
Conversely, dPCR analysis of the five most downregulated mitochondrial genes revealed increased absolute transcript abundance at 6 dpi compared with uninfected controls. This apparent discrepancy, may reflect a compensatory transcriptional response to mitochondrial dysfunction, as exemplified by Ndufb7 and Ndufs6 that encode Complex I subunits, and is consistent with prior reports indicating disrupted mitochondrial integrity following SARS-CoV-2 infection (35). Despite Complex I-associated transcript increases, COXIV protein, a marker for Complex IV in the mitochondrial electron transport chain, was significantly reduced at both 4- and 6-dpi compared to the control (Supplementary Figure S3), suggesting complex post-infection mitochondrial remodeling, potentially including impaired assembly or stability of oxidative phosphorylation complexes at the protein level. Collectively, these data validated several key pathways significantly altered in 4- and 6-dpi, and reveal a transcriptional–translational discrepancy that may reflect compensatory host responses aimed at preserving cellular homeostasis while mounting antiviral defenses.
4 Discussion
In this study, we characterize the neuroinvasive potential of SARS-CoV-2 and the associated temporal transcriptional remodeling in the brains of intranasally infected a KRT18-hACE2 transgenic mice. Our findings demonstrate direct viral effects, evidenced by viral infiltration of the brain with both neurotropism and gliotropism, culminating in increased apoptosis, particularly among cortical projection neurons. In parallel, viral invasion elicits secondary immune responses, including robust inflammation marked by upregulation of chemokines involved in immune cell recruitment (e.g., CCL5, CXCL10), and activation of myeloid cells, accompanied by a shift in host defense program toward interferon-associated and adaptive immune response.
Mechanistic insight was strengthened through an integrated analytical strategy combining cytokine/chemokine profiling, RNA microarray with GO enrichment analysis, and dPCR validation of selected inflammatory and mitochondria-associated genes. This layered approach enables unbiased identification of inflammation-driven transcriptional programs and mitochondrial stress pathways, followed by absolute quantification of key immune mediators. Such analyses are particular informative in the K18-hACE2 model, given its pronounced neuroinflammatory phenotype, and allowed delineation of the temporal landscape of innate immune activation and host defense, as well as identification of potentially vulnerable molecular nodes.
Cytokine/chemokine array analysis revealed marked upregulation of immune modulators including CXCL9, CXCL10, CCL12, and IL-12p40, consistent with enhanced immune cell recruitment to infected sites. These findings are concordant with RNA microarray data at 6 dpi, which demonstrated strong activation of innate antiviral defense pathways, and align with previous reports identifying these molecules as indicator of neuroinflammation and biomarkers associated with COVID-19 disease progression (36–38). In contrast, CX3CL1 (fractalkine), a neuron-specific chemokine that mediate neuron-microglia communication and has been reported to be elevated in COVID-19 patients (39), was significantly reduced, suggesting potential disruption of protective neuron–microglia crosstalk and a maladaptive response at 6 dpi.
Noteworthily, the combined analyses revealed transcription – translation divergence across several pathways at 6 dpi, suggesting altered post-transcriptional or post-translational regulation during infection. This observation is consistent with our GO analysis showing sustained enrichment of ubiquitin-associated functions. Notably, elevated RBP4 proteins levels observed at 6dpi, whereas the retinoid signaling was among the most downregulated molecular function at 4 dpi, suggesting temporal and multi−layered modulation of retinoid-associated processes. Given prior reports implicating RBP4 in antiviral responses (34) and inflammation (40), these changes may reflect adaptive remodeling of host pathways under infectious stress. Likewise, IL-1β, TNF- α, and IL-6, cytokines frequently linked to SARS-CoV-2-associated immunopathology (6, 41, 42), exhibited time-dependent upregulations at transcript level but were not detected in cytokine profiling of 6 dpi brain tissue. As these cytokines are subjected to stringent post-transcriptional control (43, 44) and rapid turnover (45, 46), this RNA-protein discordance may indicate regulatory restraint of inflammatory signaling at later stages of infection rather than an absence of inflammatory activation. Given that chronic elevation of TNF-α and IL-6 beyond the acute inflammatory phase is known to disrupt brain homeostasis (10), impair neurogenesis (45, 47), and compromise blood–brain barrier integrity (10, 48), the observed discordance between RNA and protein levels may reflect a host regulatory mechanism that constrains excessive inflammatory signaling while maintaining antiviral defense. Importantly, such RNA–protein mismatches under severe infectious conditions highlight a regulatory layer that may shape downstream neurological outcomes. Further studies delineating this dysregulated yet restrained inflammatory signaling may contribute not only to a better understanding of COVID-19–associated neuropathology but also provide broader insight into neurodegenerative and neuroinflammatory disorders in which inflammation plays a central pathogenic role.
Temporal GO analysis further highlighted stage-specific biological alterations in host responses. At 4 dpi, transcriptional changes were primarily associated with cellular architecture, specialization, synaptogenesis, and intracellular trafficking, with enrichment of ubiquitin- and GTPase-associated functions. These enrichments, while not experimentally validated at the protein level, are consistent with prior reports indicating that viral infection can perturb vesicular trafficking and synaptic transport through modulation of protein ubiquitination (49). Given the established role of small GTPases in vesicle docking and neurotransmitter release, such alterations may contribute to disrupted neuronal communication and synaptic homeostasis (50). In parallel, enrichments of GTPase-associated pathways at this stage may also reflect virus-driven intracellular trafficking during the post-replication egress phase (51), in line with previous observations of SARS-CoV-2-associated synaptic dysfunction (52).
By 6 dpi, innate immune pathways, DNA repair and transcriptional regulatory processes dominated the transcriptional landscape, accompanied by increased numbers of TUNEL+ cells, collectively indicating heightened cellular stress and engagement of mechanisms involved in maintaining chromosome integrity and transcriptional fidelity. In contrast to early-stage alterations affecting cellular structure and synaptic function, downregulated pathways at 6 dpi were most prominently associated with energy metabolism, supported by a marked reduction in COXIV protein expression, suggesting progressive metabolic compromise and cellular exhaustion at later stages of infection. Additionally, the persistent enrichment of ubiquitin-associated pathways at this stage may reflect either continued viral exploitation of ubiquitin and UBL systems to facilitate replication, trafficking and immune evasion (53–55), or host-driven targeting of viral proteins for degradation as part of antiviral defense mechanisms (56). Notably, we also observed enrichment of ECM-associated pathways, raising the possibility that SARS-CoV-2 infection may influence ECM remodeling in brain-resident cells. Given the established relationship between ECM dysregulation and disease severity in COVID-19 lung pathology (57), analogous alterations in the brain could plausibly contribute to impaired tissue integrity and neurological manifestations, although further investigation is required.
Collectively, our findings delineate a transcriptional landscape charactered by discordant RNA-protein regulation of immune pathway-associated factors, together with evidence of metabolic and energy exhaustion at later stage of infectious, within which ubiquitin-associated pathways appears prominently engaged. The KRT18-hACE2 mouse model enables efficient and consistent viral entry to CNS and induces a robust, high penetrance neuroinflammatory and innate immune response (58, 59). Therefore, our results aim to capture the maximal transcriptional and inflammatory responses of brain tissue, rather than to model mild or resolving infection. Insights from this fulminant phenotype are particularly relevant for understanding the upper bounds of virus-driven neuroinflammation, which may be applicable to rare but clinically devastating manifestations such as severe acute neuro-COVID and COVID-associated encephalopathy, conditions that have been linked to pronounced systemic and CNS inflammatory signatures (60, 61). Importantly, since post-acute neurological sequelae can also arise following mild infection, defining pathways activated under extreme inflammatory pressure provides critical insight into mechanisms that, if incompletely resolved, may predispose to persistent neurological dysfunction. These findings establish a mechanistic framework that can be tested and refined in milder models and human systems to advance understanding of long COVID–associated neurological sequelae.
The current study has several limitations. Our layered analyses delineate pathway signatures within a global landscape of transcriptional alterations at 4 and 6 dpi; however, functional validation of the identified metabolic impairments remains limited. Accordingly, these findings should be regarded as primarily descriptive and hypothesis-generating. Further mechanistic studies incorporating cell-type–specific markers, lineage-resolved approaches, and region-specific analyses will be required to determine whether these upregulations primarily reflect viral exploitation or host-driven defense mechanism, whether they exhibit stage-dependent differences, and to delineate the respective contributions of resident and infiltrating immune cells during SARS-CoV-2–associated neuroinflammation. In addition, evaluation of sex-dependent differences in neuroinflammatory and metabolic responses may provide important biological context, particularly in light of well-documented gender discrepancy in neuroimmune regulation.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Ethics statement
The animal study was approved by the Institutional Animal Care and Use Committee (IACUC) of the Korea Research Institute of Chemical Technology (permit no. 2021-8A-04-01). The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
D-GA: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Validation, Visualization, Writing – review & editing. NM: Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing. D-JJ: Formal analysis, Software, Visualization, Writing – review & editing. B-SM: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This work was supported by the Standardization research for the practical application of alternative animal testing funded by the Ministry of Health & Welfare and the Ministry of Food and Drug Safety (RS-2024–00332014 to B-SM), a grant of the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health & Welfare, Republic of Korea (RS-2025–02303821 to B-SM), the Basic Science Research Program through the National Research Foundation (NRF) funded by the Ministry of Education (RS-2024–00455950 to B-SM) and the 2023 Yeungnam University Research Grant (223A580045).
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2026.1716597/full#supplementary-material
References
1
Zhou P Yang X-L Wang X-G Hu B Zhang L Zhang W et al . A pneumonia outbreak associated with a new coronavirus of probable bat origin. nature. (2020) 579:270–3. doi: 10.1038/s41586-020-2012-7
2
Zhu N Zhang D Wang W Li X Yang B Song J et al . A novel coronavirus from patients with pneumonia in China, 2019. New Engl J Med. (2020) 382:727–33. doi: 10.1056/NEJMoa2001017
3
Bauer L Laksono BM de Vrij FMS Kushner SA Harschnitz O van Riel D . The neuroinvasiveness, neurotropism, and neurovirulence of SARS-CoV-2. Trends Neurosci. (2022) 45:358–68. doi: 10.1016/j.tins.2022.02.006
4
Song E Zhang C Israelow B Lu-Culligan A Prado AV Skriabine S et al . Neuroinvasion of SARS-CoV-2 in human and mouse brain. J Exp Med. (2021) 218(3):e20202135. doi: 10.3892/etm.2021.10943
5
Yachou Y El Idrissi A Belapasov V Ait Benali S . Neuroinvasion, neurotropic, and neuroinflammatory events of SARS-CoV-2: understanding the neurological manifestations in COVID-19 patients. Neurol Sci. (2020) 41:2657–69. doi: 10.1007/s10072-020-04575-3
6
Oh S-J Kumari P Auroni TT Stone S Pathak H Elsharkawy A et al . Upregulation of Neuroinflammation-associated genes in the brain of SARS-CoV-2-infected mice. Pathogens. (2024) 13:528. doi: 10.3390/pathogens13070528
7
Orsini A Corsi M Santangelo A Riva A Peroni D Foiadelli T et al . Challenges and management of neurological and psychiatric manifestations in SARS-CoV-2 (COVID-19) patients. Neurol Sci. (2020) 41:2353–66. doi: 10.1007/s10072-020-04544-w
8
Helms J Kremer S Merdji H Clere-Jehl R Schenck M Kummerlen C et al . Neurologic features in severe SARS-coV-2 infection. N Engl J Med. (2020) 382:2268–70. doi: 10.1056/NEJMc2008597
9
Spudich S Nath A . Nervous system consequences of COVID-19. science. (2022) 375:267–9. doi: 10.1126/science.abm2052
10
Yang J Ran M Li H Lin Y Ma K Yang Y et al . New insight into neurological degeneration: Inflammatory cytokines and blood–brain barrier. Front Mol Neurosci. (2022) 15:1013933. doi: 10.3389/fnmol.2022.1013933
11
Li C Verduzco-Gutierrez M . Neurologic and neuromuscular sequelae of COVID-19. Phys Med Rehabil Clinics North America. (2023) 34(3). doi: 10.1016/j.pmr.2023.04.002
12
Pang Z Tang A He Y Fan J Yang Q Tong Y et al . Neurological complications caused by SARS-CoV-2. Clin Microbiol Rev. (2024) 37:e00131–24. doi: 10.1128/cmr.00131-24
13
Freitas NL Deus JVC Sampaio K Gomes YCP Torres RC Brandão CO et al . Central nervous system and systemic inflammatory networks associated with acute neurological outcomes in COVID-19. Sci Rep. (2025) 15:24154. doi: 10.1038/s41598-025-08632-9
14
Wesselingh R . Prevalence, pathogenesis and spectrum of neurological symptoms in COVID-19 and post-COVID-19 syndrome: a narrative review. Med J Aust. (2023) 219:230–6. doi: 10.5694/mja2.52063
15
Winkler ES Bailey AL Kafai NM Nair S McCune BT Yu J et al . SARS-CoV-2 infection of human ACE2-transgenic mice causes severe lung inflammation and impaired function. Nat Immunol. (2020) 21:1327–35. doi: 10.1038/s41590-020-0778-2
16
Jiang R-D Liu M-Q Chen Y Shan C Zhou Y-W Shen X-R et al . Pathogenesis of SARS-CoV-2 in transgenic mice expressing human angiotensin-converting enzyme 2. Cell. (2020) 182:50–8. doi: 10.1016/j.cell.2020.05.027
17
Kumari P Rothan HA Natekar JP Stone S Pathak H Strate PG et al . Neuroinvasion and encephalitis following intranasal inoculation of SARS-CoV-2 in K18-hACE2 mice. Viruses. (2021) 13:132. doi: 10.3390/v13010132
18
Golden JW Cline CR Zeng X Garrison AR Carey BD Mucker EM et al . Human angiotensin-converting enzyme 2 transgenic mice infected with SARS-CoV-2 develop severe and fatal respiratory disease. JCI Insight. (2020) 5:e142032. doi: 10.1172/jci.insight.142032
19
Rothan HA Kumari P Stone S Natekar JP Arora K Auroni TT et al . SARS-CoV-2 infects primary neurons from human ACE2 expressing mice and upregulates genes involved in the inflammatory and necroptotic pathways. Pathogens. (2022) 11:257. doi: 10.3390/pathogens11020257
20
Sui L Li L Zhao Y Zhao Y Hao P Guo X et al . Host cell cycle checkpoint as antiviral target for SARS-CoV-2 revealed by integrative transcriptome and proteome analyses. Signal Transduct Target Ther. (2023) 8:21. doi: 10.1038/s41392-022-01296-1
21
Singh KK Chaubey G Chen JY Suravajhala P . Decoding SARS-CoV-2 hijacking of host mitochondria in COVID-19 pathogenesis. Am J Physiology-Cell Physiol. (2020) 319:C258–67. doi: 10.1152/ajpcell.00224.2020
22
Rahmati M Yon DK Lee SW Soysal P Koyanagi A Jacob L et al . New-onset neurodegenerative diseases as long-term sequelae of SARS-CoV-2 infection: a systematic review and meta-analysis. J Med Virol. (2023) 95:e28909. doi: 10.1002/jmv.28909
23
Seo D Choi Y Jeong E Bang S Lee J-S Jang I-H et al . Distinct brain alterations and neurodegenerative processes in cognitive impairment associated with post-acute sequelae of COVID-19. Nat Commun. (2025) 16:10552. doi: 10.1038/s41467-025-65597-z
24
Lee S Yoon GY Lee SJ Kwon Y-C Moon HW Kim Y-J et al . Immunological and pathological peculiarity of severe acute respiratory syndrome coronavirus 2 beta variant. Microbiol Spectr. (2022) 10:e02371–22. doi: 10.1128/spectrum.02371-22
25
Kim J-M Chung Y-S Jo HJ Lee N-J Kim MS Woo SH et al . Identification of coronavirus isolated from a patient in korea with COVID-19. Osong Public Health Res Perspect. (2020) 11:3–7. doi: 10.24171/j.phrp.2020.11.1.02
26
Oladunni FS Park J-G Pino PA Gonzalez O Akhter A Allué-Guardia A et al . Lethality of SARS-CoV-2 infection in K18 human angiotensin-converting enzyme 2 transgenic mice. Nat Commun. (2020) 11:6122. doi: 10.1038/s41467-020-19891-7
27
Moreau GB Burgess SL Sturek JM Donlan AN Petri WA Jr. Mann BJ . Evaluation of K18-hACE2 mice as a model of SARS-CoV-2 infection. Am J Trop Med Hygiene. (2020) 103:1215. doi: 10.4269/ajtmh.20-0762
28
Verge GM Milligan ED Maier SF Watkins LR Naeve GS Foster AC . Fractalkine (CX3CL1) and fractalkine receptor (CX3CR1) distribution in spinal cord and dorsal root ganglia under basal and neuropathic pain conditions. Eur J Neurosci. (2004) 20:1150–60. doi: 10.1111/j.1460-9568.2004.03593.x
29
Bivona G Iemmolo M Ghersi G . CX3CL1 pathway as a molecular target for treatment strategies in Alzheimer’s disease. Int J Mol Sci. (2023) 24:8230. doi: 10.3390/ijms24098230
30
Steinhoff JS Lass A Schupp M . Biological functions of RBP4 and its relevance for human diseases. Front Physiol. (2021) 12:659977. doi: 10.3389/fphys.2021.659977
31
Scherer KM Mascheroni L Carnell GW Wunderlich LC Makarchuk S Brockhoff M et al . SARS-CoV-2 nucleocapsid protein adheres to replication organelles before viral assembly at the Golgi/ERGIC and lysosome-mediated egress. Sci Adv. (2022) 8:eabl4895. doi: 10.1126/sciadv.abl4895
32
Zhang H Zhang H . Entry, egress and vertical transmission of SARS-CoV-2. J Mol Cell Biol. (2021) 13:168–74. doi: 10.1093/jmcb/mjab013
33
Sergio MC Ricciardi S Guarino AM Giaquinto L De Matteis MA . Membrane remodeling and trafficking piloted by SARS-CoV-2. Trends Cell Biol. (2024) 34:785–800. doi: 10.1016/j.tcb.2023.12.006
34
Tong L Wang L Liao S Xiao X Qu J Wu C et al . A retinol derivative inhibits SARS-CoV-2 infection by interrupting spike-mediated cellular entry. MBio. (2022) 13:e01485–22. doi: 10.1128/mbio.01485-22
35
Shin HJ Lee W Ku KB Yoon GY Moon H-W Kim C et al . SARS-CoV-2 aberrantly elevates mitochondrial bioenergetics to induce robust virus propagation. Signal Transduct Target Ther. (2024) 9:125. doi: 10.1038/s41392-024-01836-x
36
Lorè NI De Lorenzo R Rancoita PM Cugnata F Agresti A Benedetti F et al . CXCL10 levels at hospital admission predict COVID-19 outcome: hierarchical assessment of 53 putative inflammatory biomarkers in an observational study. Mol Med. (2021) 27:129. doi: 10.1186/s10020-021-00390-4
37
Callahan V Hawks S Crawford MA Lehman CW Morrison HA Ivester HM et al . The pro-inflammatory chemokines CXCL9, CXCL10 and CXCL11 are upregulated following SARS-CoV-2 infection in an AKT-dependent manner. Viruses. (2021) 13:1062. doi: 10.3390/v13061062
38
Michael BD Dunai C Needham EJ Tharmaratnam K Williams R Huang Y et al . Para-infectious brain injury in COVID-19 persists at follow-up despite attenuated cytokine and autoantibody responses. Nat Commun. (2023) 14:8487. doi: 10.1038/s41467-023-42320-4
39
Hsu R-J Yu W-C Peng G-R Ye C-H Hu S Chong PCT et al . The role of cytokines and chemokines in severe acute respiratory syndrome coronavirus 2 infections. Front Immunol. (2022) 13:832394. doi: 10.3389/fimmu.2022.832394
40
Rubin LP Ross AC Stephensen CB Bohn T Tanumihardjo SA . Metabolic effects of inflammation on vitamin A and carotenoids in humans and animal models. Adv Nutr. (2017) 8:197–212. doi: 10.3945/an.116.014167
41
Vanderheiden A Hill JD Jiang X Deppen B Bamunuarachchi G Soudani N et al . Vaccination reduces central nervous system IL-1β and memory deficits after COVID-19 in mice. Nat Immunol. (2024) 25:1158–71. doi: 10.1038/s41590-024-01868-z
42
Fekete R Simats A Bíró E Pósfai B Cserép C Schwarcz AD et al . Microglia dysfunction, neurovascular inflammation and focal neuropathologies are linked to IL-1-and IL-6-related systemic inflammation in COVID-19. Nat Neurosci. (2025) 28(3):558–76. doi: 10.1038/s41593-025-01871-z
43
Mino T Takeuchi O . Post-transcriptional regulation of cytokine mRNA controls the initiation and resolution of inflammation. Biotechnol Genet Eng Rev. (2013) 29:49–60. doi: 10.1080/02648725.2013.801236
44
Uchida Y Chiba T Kurimoto R Asahara H . Post-transcriptional regulation of inflammation by RNA-binding proteins via cis-elements of mRNAs. J Biochem. (2019) 166:375–82. doi: 10.1093/jb/mvz067
45
Kummer KK Zeidler M Kalpachidou T Kress M . Role of IL-6 in the regulation of neuronal development, survival and function. Cytokine. (2021) 144:155582. doi: 10.1016/j.cyto.2021.155582
46
Frankola KA Greig NH Luo W Tweedie D . Targeting TNF-alpha to elucidate and ameliorate neuroinflammation in neurodegenerative diseases. CNS Neurol Disorders-Drug Targets (Formerly Curr Drug Targets-CNS Neurol Disorders). (2011) 10:391–403. doi: 10.2174/187152711794653751
47
Belarbi K Jopson T Tweedie D Arellano C Luo W Greig NH et al . TNF-α protein synthesis inhibitor restores neuronal function and reverses cognitive deficits induced by chronic neuroinflammation. J Neuroinflamm. (2012) 9:23. doi: 10.1186/1742-2094-9-23
48
Zhou M-Y Xie F Zhao Y Wang X Sun Z-W Qian L-J . Stress-induced regional characteristics of the blood-brain barrier. Brain Res Bull. (2025) 231:111548. doi: 10.1016/j.brainresbull.2025.111548
49
Wang Q Peng W Yang Y Wu Y Han R Ding T et al . Proteome and ubiquitinome analyses of the brain cortex in K18-hACE2 mice infected with SARS-CoV-2. Iscience. (2024) 27:110602. doi: 10.1016/j.isci.2024.110602
50
Partiot E Hirschler A Colomb S Lutz W Claeys T Delalande F et al . Brain exposure to SARS-CoV-2 virions perturbs synaptic homeostasis. Nat Microbiol. (2024) 9:1189–206. doi: 10.1038/s41564-024-01657-2
51
Hou W Wang S Wu H Xue L Wang B Wang S et al . Small GTPase—a key role in host cell for coronavirus infection and a potential target for coronavirus vaccine adjuvant discovery. Viruses. (2022) 14:2044. doi: 10.3390/v14092044
52
Ngo BT Marik P Kory P Shapiro L Thomadsen R Iglesias J et al . The time to offer treatments for COVID-19. Expert Opin Investigat Drugs. (2021) 30:505–18. doi: 10.1080/13543784.2021.1901883
53
Zhao M Zhang M Yang Z Zhou Z Huang J Zhao B . Role of E3 ubiquitin ligases and deubiquitinating enzymes in SARS-CoV-2 infection. Front Cell Infect Microbiol. (2023) 13:1217383. doi: 10.3389/fcimb.2023.1217383
54
Zheng Y Yang H Zhang X Gao C . Regulation of SARS-CoV-2 infection and antiviral innate immunity by ubiquitination and ubiquitin-like conjugation. Biochim Biophys Acta (BBA)-Gene Regul Mech. (2023) 1866:194984. doi: 10.1016/j.bbagrm.2023.194984
55
Viswanathan K Früh K DeFilippis V . Viral hijacking of the host ubiquitin system to evade interferon responses. Curr Opin Microbiol. (2010) 13:517–23. doi: 10.1016/j.mib.2010.05.012
56
Yu M Li J Gao W Li Z Zhang W . Multiple E3 ligases act as antiviral factors against SARS-CoV-2 via inducing the ubiquitination and degradation of ORF9b. J Virol. (2024) 98:e01624–23. doi: 10.1128/jvi.01624-23
57
Huang J-J Wang C-W Liu Y Zhang Y-Y Yang N-B Yu Y-C et al . Role of the extracellular matrix in COVID-19. World J Clin cases. (2023) 11:73. doi: 10.12998/wjcc.v11.i1.73
58
Carossino M Kenney D O’Connell AK Montanaro P Tseng AE Gertje HP et al . Fatal neurodissemination and SARS-CoV-2 tropism in K18-hACE2 mice is only partially dependent on hACE2 expression. Viruses. (2022) 14:535. doi: 10.3390/v14030535
59
Dong W Mead H Tian L Park J-G Garcia JI Jaramillo S et al . The K18-human ACE2 transgenic mouse model recapitulates non-severe and severe COVID-19 in response to an infectious dose of the SARS-CoV-2 virus. J Virol. (2022) 96:e00964–21. doi: 10.1128/JVI.00964-21
60
Etter MM Martins TA Kulsvehagen L Pössnecker E Duchemin W Hogan S et al . Severe Neuro-COVID is associated with peripheral immune signatures, autoimmunity and neurodegeneration: a prospective cross-sectional study. Nat Commun. (2022) 13:6777. doi: 10.1038/s41467-022-34068-0
61
Duve K Petakh P Kamyshnyi O . COVID-19-associated encephalopathy: connection between neuroinflammation and microbiota-gut-brain axis. Front Microbiol. (2024) 15:1406874. doi: 10.3389/fmicb.2024.1406874
Summary
Keywords
chemokines, cytokines, immune signaling, inflammation, neuroinvasion, SARS-CoV-2
Citation
Ahn D-G, Mai NTQ, Jeong D-J and Moon B-S (2026) Neuroinflammatory and transcriptional dynamics during SARS-CoV-2 infection in KRT18-hACE2 mouse brain. Front. Immunol. 17:1716597. doi: 10.3389/fimmu.2026.1716597
Received
30 September 2025
Revised
28 December 2025
Accepted
09 January 2026
Published
02 February 2026
Volume
17 - 2026
Edited by
Guirong Wang, Upstate Medical University, United States
Reviewed by
Fabrícia Lima Fontes-Dantas, Rio de Janeiro State University, Brazil
Akinkunmi Lawal, Upstate Medical University, United States
Updates
Copyright
© 2026 Ahn, Mai, Jeong and Moon.
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: Byoung-San Moon, bsmoon@yu.ac.kr
†These authors have contributed equally to this work
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.