Activation of NLRP3 Inflammasome by Virus-Like Particles of Human Polyomaviruses in Macrophages

Viral antigens can activate phagocytes, inducing inflammation, but the mechanisms are barely explored. The aim of this study is to investigate how viral oligomeric proteins of different structures induce inflammatory response in macrophages. Human THP-1 cell line was used to prepare macrophages that were treated with filamentous nucleocapsid-like particles (NLPs) of paramyxoviruses and spherical virus-like particles (VLPs) of human polyomaviruses. The effects of viral proteins on cell viability, pro-inflammatory cytokines’ production, and NLRP3 inflammasome activation were investigated. Filamentous NLPs did not induce inflammation while spherical VLPs mediated inflammatory response followed by NLRP3 inflammasome activation. Inhibitors of cathepsins and K+ efflux decreased IL-1β release and cell death, indicating a complex inflammasome activation process. A similar activation pattern was observed in primary human macrophages. Single-cell RNAseq analysis of THP-1 cells revealed several cell activation states different in inflammation-related genes. This study provides new insights into the interaction of viral proteins with immune cells and suggests that structural properties of oligomeric proteins may define cell activation pathways.


INTRODUCTION
The components of innate immunity such as macrophages play a key role in the onset and progression of inflammatory and age-related diseases (1,2). Macrophages are considered as a potential target for treatment of many diseases; therefore, molecular mechanisms related to their activation are currently a top priority in research (3). Besides, macrophages are known to recognize the structural properties of the activation agents and influence the inflammatory response via inflammasome activation (4,5). Inflammasomes are intracellular protein complexes representing important components of the innate immune system (6). The best-described representative is NLRP3 inflammasome. It contains three major components-nucleotide-binding and oligomerization domain-like receptor, apoptosisassociated speck-like protein containing CARD (ASC), and pro-caspase-1. NLRP3 inflammasome assembly results in IL-1b release and inflammatory cell death-pyroptosis (7). Endogenous and external factors can trigger its assembly. Activation of NLRP3 inflammasome is associated with various diseases, including gout and Alzheimer's disease (8). Our previous study showed that NLRP3 inflammasome is activated by amyloid beta (Ab), both oligomers and protofibrils (9). Inflammasome activation was also induced by a-synuclein (10) and tau oligomers (11). However, the most explored inflammasome triggers are polymeric nanoparticles (12), cholesterol crystals (13), and airborne pollutants such as silica (14). The structural properties of synthetic nanoparticles determine the outcome of cell activation (15). Previously, we have demonstrated that spherical oligomeric proteins of viral origin induce inflammatory responses in macrophages, but the mechanism and molecular components of this process are not yet confirmed (16). In addition, another study demonstrated that the nucleocapsid (N) protein of Zika virus activated the inflammasome (17). Repetitive and latent viral infections are potential agents of inflammation (18). This frame of reference draws attention to the latent viral infections and their role in inflammation.
It is well known that some viruses lie dormant within the host after an acute infection. According to the World Health Organization, 50%-80% of the world population is infected by polyomaviruses (PyVs) during childhood (19). The effects of these viruses and their antigens on innate immunity throughout life are not fully understood. For example, herpes virus in its latent form can cause multisymptom illnesses with a broad range of simultaneous symptoms, such as cognitive disorders, depression, and fatigue (18). However, the mechanism of multisymptomatic illness is unclear. Investigation of inflammatory responses induced by proteins of viruses causing acute and latent infection may indicate whether these viral proteins alone are capable of activating the immune system.
Inflammasome activation by latent viruses was demonstrated in different viral infections, and it could be related to their deleterious effects. So far, the identified viral triggers of the inflammasome are viral nucleic acids and viroporins, which make pores in cell membranes and disturb ion homeostasis and induce mitochondrial stress or ROS production (20). For example, IFI16 inflammasome is activated by Epstein-Barr virus DNA during its latency (21). IFI16 also senses Kaposi's sarcoma-associated herpesvirus genome (22). This could be a potential mechanism of inflammation associated with Kaposi's sarcoma. In a mouse model, it was demonstrated that mouse cytomegalovirus-induced activation of NLRP3/NLRP1 may contribute to pulmonary fibrosis caused by latent cytomegalovirus infection (23). However, the ability of viral surface antigens and other structural proteins to activate the inflammasome is still insufficiently studied. Another question is whether viral proteins or other factors during their latency are produced, which could mediate inflammation. For example, herpes simplex virus-1 (HSV-1) can induce Ab production in infected brain cells in its latent phase (24). Furthermore, these amyloids promote neurotoxicity and neuroinflammation. HSV-1 is known to interact with AD-related genes and proteins (25). There is an increasing number of evidence that Ab has a protective function by inhibiting viral replication and virus attachment to neurons (26). In response to HSV-1, neurons start to produce Ab to fight infection; however, for genetically susceptible persons, Ab accumulation can lead to AD (27,28). The impact of virus reactivation to neuroinflammation and neurodegeneration in mouse models reveals the role of the immune system in disease progression (29,30). Studies show that HSV-1 reactivation promotes cognitive impairment (31,32). However, processes during viral latency period are less investigated. It was demonstrated that inflammatory cytokines are produced even during viral latency (33,34).
It is important to investigate how the components of innate and adaptive immunity are activated, including the inflammasome, interferons, and other inflammatory molecules, in viral infections or during vaccination. VLPs are already established as carriers and adjuvants in vaccines. Components of vaccines can activate the inflammasome as it was reported for the ISCOMATRIX adjuvant (35). Studying how different viral proteins especially VLPs stimulate the immune cells would allow the selection of better tools for vaccination.
Inflammatory reactions induced by synthetic polymeric particles in macrophages vary depending on their size and shape (4,5). Similarly, viral proteins may also determine cellular response depending on their structural properties. Our study was aimed to investigate how recombinant viral proteins induce inflammatory response in macrophages focusing on inflammasome activation. The VLPs derived from major capsid protein VP1 of human PyVs-Karolinska Institute (KI) PyV and Merkel cell (MC) PyV, self-assembling to spherical particles about 20-60 nm in diameter (36)-were selected as typical representatives of spherical oligomeric proteins. The NLPs of measles and mumps viruses were chosen as a model of oligomeric proteins forming filamentous rod-shaped structures (37,38). We showed that inflammatory responses and activation of NLRP3 inflammasome in macrophages depend on the structural properties of viral proteins. This study provides new insights into the ability of multimeric viral antigens to induce inflammatory response in innate immune cells.

Human Macrophage Cell Culture
Human macrophage cell cultures were prepared by differentiation of THP-1 cells (39). The cells were seeded in the 24-well plate at a density of 0.125 × 10 6 /well using RPMI 1640 medium supplemented with 10% FBS and 1% P/S and differentiated to macrophages using 100 ng/ml of PMA. After 48 h of differentiation, the medium was replaced with the fresh medium without PMA and cells were left to rest for another 24 h. After the rest period, the cells are differentiated into macrophages and were used in experiments with viral proteins. These macrophage-like cells were used in the study and called THP-1 macrophages.
THP-1 macrophages were washed once with serum-free RPMI and treated with viral proteins for 24 h. Viral proteins were prepared in PBS, so control when PBS was added instead of viral proteins was used. As a positive control, the inflammasome inducer nigericin was used at 10 mM concentration. MCC950, which selectively inhibits the NLRP3 inflammasome, was used at 1 mM concentration and added 30 min before the treatment. Inhibitors of cathepsins were used with the following concentrations: CA-074 Me at 2 mM and 10 mM, and K777 at 15 mM, added 30 min before the treatment. Another inhibitor, glibenclamide, which blocks K + efflux, was used at 50 μM concentration, 30 min before the treatment. After incubation, cell culture supernatants were collected and stored at −20°C for further cytokine analysis. Supernatants for LDH assay were used instantly.
Primary human macrophages were purchased from Lonza (#4W-700). Human macrophages were derived from CD14 + human monocytes of one donor. Cryopreserved cells were thawed and cultured for 2 days in RPMI 1640 medium supplemented with 10% FBS and 1% P/S before treatment. The cells were treated as THP-1 macrophages in serum-free RPMI. Cells for FLICA assay were plated in an Ibidi 96-well m-plate and for ELISA in a TPP 24-well plate.

VLP Production, Purification, and Analysis
Purification of VLPs of recombinant PyV VP1 proteins and electron microscopy were carried out as described previously (36). Briefly, Saccharomyces cerevisiae yeast biomass after the induction of recombinant proteins synthesis was mechanically homogenized in DB450 buffer (450 mM NaCl, 1 mM CaCl 2 , 0.25 M L-Arginine, and 0.001% Triton X-100 in 10 mM Tris/HCl buffer, pH 7.2) with 2 mM PMSF and EDTA-free Complete Protease Inhibitors Cocktail tablets (Roche Diagnostics, Mannheim, Germany), and its supernatant was transferred onto 30%-60% sucrose gradient. After overnight centrifugation (at 4°C) at 100,000×g (Beckman Coulter Optima L-90 ultracentrifuge), collected 2-ml fractions were analyzed by SDS-PAGE. The mixture of fractions containing PyV VP1 proteins diluted in DB150 buffer (150 mM NaCl, 1 mM CaCl 2 , 0.25 M L-Arginine, and 0.001% Triton X-100 in 10 mM Tris/HCl buffer, pH 7.2) and VLPs were concentrated by ultracentrifugation at 100,000×g for 4 h (at 4°C). Thereafter, pellets containing VP1 were subjected to ultracentrifugation overnight on CsCl gradient (1.23-1.46 g/ml density) at 4°C. One-milliliter fractions of formed gradient were collected and analyzed by SDS-PAGE. Positive fractions were pooled, diluted in DB150 buffer, and concentrated as described above. The isolated VLPs were dissolved in PBS, dialyzed, and stored in PBS with 50% glycerol. The VLP formation was verified by examination of the purified proteins using Morgagni-268 electron microscope (FEI, Inc., Hillsboro, OR, USA). The protein samples were placed on 400-mesh carbon-coated palladium grids (Agar Scientific, Stansted, UK) and stained with 2% aqueous uranyl acetate.

Cell Cytotoxicity Assays
Cell cytotoxicity was measured using lactate dehydrogenase (LDH) release assay (LDH cytotoxicity detection Kit). A quantity of 50 ml of cell supernatants was used to perform the cytotoxicity assay according to the manufacturer's protocol. Briefly, 50 ml of supernatant was mixed with freshly prepared LDH reagent and incubated at 37°C. After 30 min, absorbance was measured at 490 nm using a Multiskan GO microplate spectrophotometer (Thermo Fisher Scientific Oy, Finland).
To determine cell viability, PI/Hoechst nuclear staining was also used. Nuclei were stained with 1.25 mg/ml PI and 1 mg/ml Hoechst33342 in cell culture medium for 30 min. The cells were washed with PBS. The fluorescence signal was measured by taking photos automatically with an EVOS FL Auto fluorescence microscope. Images were taken using a 20× objective. Viability was quantified according to a ratio of PI (dead cells) and Hoechst (all cells), expressed in percentages.

Quantitation of Cytokines in Cell Culture Supernatants
ELISA kits for the measurement of human cytokine-IL-6, IL-10, IL-1b, and TNF-a-levels in cell culture supernatants were used (#88-7066, #88-7106, #88-7261, and #88-7346, Thermo Fisher Scientific, USA). ELISA kits are based on the sandwich immunoassay technique. Supernatants were used diluted up to 1:600. All procedures were performed according to the manufacturer's protocols. In the last step, 3,3',5,5'tetramethylbenzidine (TMB) substrate solution was added to each well. The plates were monitored for 15 min for color development, the reaction in wells was stopped with 3.6% H 2 SO 4 solution, and the wells were read at 450 nm with reference wavelength at 620 nm using a Multiskan GO microplate spectrophotometer. A standard curve was generated from cytokine standard, and the cytokine concentration in the samples was calculated.

Detection of ASC Speck Formation in THP-1 Macrophages
THP-1 monocytes expressing ASC fused with green fluorescent protein (GFP) was used for ASC speck formation. The cells were cultivated in RPMI 1640 medium supplemented with 10% FBS, 1% P/S, 100 mg/ml normocin, and 100 mg/ml selective antibiotic zeocin. THP-1-ASC-GFP cells were differentiated to macrophages as origin THP-1 using RPMI 1640 medium supplemented with 10% FBS, 1% P/S, and PMA. THP-1-ASC-GFP macrophages were treated with viral proteins in serum-free RPMI. Thirty minutes before the treatment termination, Hoechst33342 was added to stain cell nuclei. After 24 h of treatment, cell culture medium was replaced with FluoroBrite DMEM. ASC specks were analyzed with an EVOS FL Auto fluorescence microscope (Life Technologies, USA) by taking photos with 20× objective. Cells were counted according to the number of nuclei. ASC speck number per cell was counted using the image processing program ImageJ.

Measurement of Active Caspase-1
Active caspase-1 was detected using Fluorochrome Labeled Inhibitors of Caspases (FLICA) assay according to the manufacturer's protocol. Briefly, FLICA reagent (FAM-YVAD-FMK-caspase-1 inhibitor probe) was added after the 15-h treatment with viral proteins and incubated for 1 h. Cells were washed three times and stained with Hoechst33342 at 1 mg/ml and PI at 1.25 mg/ml. After washing, cells were analyzed directly by an EVOS FL Auto fluorescence microscope using a 20× or 40× objective.

SDS-PAGE and Western Blot Analysis
To determine cleaved caspase-1 p20 fragment, Western blot assay (WB) was performed. After cell treatment with viral proteins, supernatant was collected and centrifuged at 600×g for 10 min to remove cellular debris. Then, the protein content of the supernatants was concentrated 10× using centrifugal filters with a 10-kDa cutoff (#UFC501096, Amicon, Merck). The concentrated samples were boiled in a reducing sample buffer and separated in 4%-12% polyacrylamide gel (#NW04122BOX, Thermo Fisher Scientific) electrophoresis (PAGE) in MES SDS running buffer (#B0002 Invitrogen, Thermo Fisher Scientific). The proteins from the SDS-PAGE gel were blotted onto a 0.2-μm nitrocellulose (NC) membrane (#LC2000, Invitrogen, Thermo Fisher Scientific) by wet transfer. The membrane was blocked with 5% BSA in PBS for 1 h at RT and rinsed with TBST. The membrane was then incubated with primary antibodies in TBST with 1% BSA overnight at 4°C. The primary antibodies against human caspase-1 (clone Bally-1, #AG-20B-0048-C100, RRID: AB_2490257, Adipogene) were used at 1:1,000 dilution. Thereafter, the membrane was incubated with secondary antibodies Goat Anti-Mouse IgG (H+L)-HRP Conjugate (Bio-Rad) diluted 1:5,000 in TBST with 1% BSA for 1 h at RT. The horseradish peroxidase (HRP) enzymatic reaction was developed using chemiluminescent substrate (#34094, Thermo Fisher Scientific).

Immunocytochemistry for Studying the Uptake of VLPs by Macrophages
Cells were stained in IbidiTreat 96-well μ-plates. After the treatment, cells were washed with PBS and fixed in 4% PFA dissolved in PBS for 15 min and permeabilized with 0.1% Triton X-100 prepared in PBS for 10 min. Blocking solution-PBS containing 2% BSA-was applied for 30 min followed by two washing steps. The primary antibodies rabbit polyclonal anti-CD68 (1:100; #25747-1-AP, RRID : AB_2721140, Proteintech) and mouse anti-PyV VP1 VLPs (monoclonal antibodies of hybridoma supernatant at dilution 1:2) were added to the blocking solution and incubated overnight. The following secondary antibodies were respectively used: goat anti-rabbit (1:1,000) and goat anti-mouse (1:1,000). The secondary antibodies were applied for 2 h followed by two washing steps. Hoechst33342 was used for nuclear staining at 1 mg/ml for 30 min in PBS. The images were taken using a 40× objective. CD68 was used as a macrophage and lysosomal marker. The experiment was imaged using an EVOS FL Auto fluorescence microscope (Thermo Fisher Scientific, USA). Acquired images were processed using ImageJ (Wayne Rusband; National Institute of Health, Bethesda, MD, USA).

Microscopy
All images were taken with fluorescence microscope EVOS FL Auto (#AMAFD1000, Thermo Fisher Scientific). EVOS imaging systems use LED light cubes, which combine bright LED illumination with excitation and emission filters into single components. PI signal was detected in RFP light cube; Hoechst33342 signal-in DAPI light cube (AMEP#4650); GFP and AlexaFluor 488 signal-in GFP light cube (AMEP#4651); PI signal-in RFP light cube (AMEP#4652); AlexaFluor 594 signalin TxRed light cube (AMEP#4655). The following objective lenses were used 20× (#AMEP4682) and 40× (#AMEP4683). Acquisition software EVOS ® FL Auto v1.6 was used. Images were prepared with ImageJ program. Firstly, appropriate pseudo-color was added at images of 8-bit format. Then, the images were converted to RGB color format. Finally, brightness and contrast were adjusted equally to all images per channel. Images (fluorescent signal and object counting) were analyzed with ImageJ. The final figures were arranged using Adobe Photoshop without any brightness/contrast and color manipulations.

Cell Preparation for Single-Cell RNA Sequencing
THP-1 macrophages for RNA sequencing analysis were treated for 15 h with PyV-derived VLPs (20 μg/ml). After the treatment, cells were washed with PBS and detached with TrypLE reagent after 15 min of incubation at a cell culture incubator (37°C, 5% CO 2 ). After centrifugation, cells were washed twice with 1× RNase-free PBS (cat#AM9624, Invitrogen). After final centrifugation, cells were resuspended in RNase-free PBS and kept on ice till cell preparation for RNA sequencing analysis

Single-Cell RNA Sequencing
Single-cell RNA sequencing (scRNAseq) was performed using a modified version of the inDrops method (42,43). While the original version is based on the CEL-Seq protocol (44) that involves linear cDNA amplification by in vitro transcription, the modified version relies on template switching and cDNA amplification by PCR, as in the Smart-seq protocol (45). First, single-cell transcriptomes were barcoded in 1-nl droplets by coencapsulating: (i) barcoding hydrogel beads; (ii) the reversetranscription/lysis (RT/lysis) mix; and (iii) the cell suspension. Cell encapsulation was performed using the microfluidic device described in (43) on an Onyx platform (Droplet Genomics). Maxima H-minus (Thermo Fisher, cat. no. EP0751) was used for reverse transcription with template switching. RT was performed for 1 h at 42°C followed by heat inactivation for 5 min at 85°C. The emulsion was then broken and the pooled material was taken through library preparation for Illumina sequencing, which involves the following steps: (i) cDNA amplification (Terra PCR direct Polymerase, Takara, cat. no. 639270),

ScRNAseq Raw Data Processing
The solo-in-drops pipeline (https://github.com/jsimonas/soloin-drops), which is a wrapper around STARsolo (https:// github.com/alexdobin/STAR/blob/master/docs/STARsolo.md), was used to obtain cells × genes expression matrices. STARsolo expects as input two fastq files per library, one containing the barcode information and the other the transcript information. Meanwhile, the custom scRNAseq protocol used in the current study outputs 3 fastq files per library: barcode half 1, barcode half 2, and transcript. Solo-in-drops prepares the data for compatibility with STARsolo. STAR (version 2.7.6a) was run with the following parameters: -soloType CB_UMI_Simple, -soloUMIfiltering MultiGeneUMI, and -soloCBmatchWLtype 1MM_multi_pseudocounts. Homo sapiens (human) genome assembly GRCh38 (hg38) was used as the reference.

ScRNAseq Count Data Analysis
Jupyter notebooks with commented code for scRNAseq data analyses, including data filtering, normalization, visualization, clustering, cell population annotation, differential gene expression (DGE) analysis, cell cycle scoring, and plots used in The Scanpy toolbox was used for data analysis. Transcriptomes with fewer than 900 total counts were excluded. No filtering on the fraction of mitochondrial counts was performed. UMAP (46) was used for data visualization in 2D. Leiden clustering was used to divide the graph into clusters, which were annotated based on their gene expression profiles. Gene Ontology (GO) gene set enrichment analysis ( Figures 8D and 10E) was performed using the tool available online (http://geneontology.org) using Panther GO-Slim Biological Processes as the annotation data set. For DGE analyses ( Figure 9C, 10E), the Mann-Whitney U test was used to test for significance, and the Benjamini-Hochberg procedure was used to correct for multiple hypothesis testing.

Statistical Analysis
All statistical analyses were performed with GraphPad Prism 9.2.0 (GraphPad Software, Inc., La Jolla, CA). The data in the figures are represented as individual data points from at least 6 independent experiments using box plots (showing minimum, first quartile, median, third quartile, and maximum) or bar graphs. Independent experiments referred to as N means the number of independent cell culture preparations and n means the number of technical repeats. Normality test was carried out to test if the values come from a Gaussian distribution. Statistical comparisons of vehicle controls versus treatment were performed with one-way ANOVA in conjunction with a Tukey's multiple comparison test or Student's t-test. A Kruskal-Wallis test with Dunn's post hoc test was used for non-parametric data. Differences with p value less than 0.05 were considered to be statistically significant: *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001. To emphasize non-significant results, ns was used.

RESULTS
We investigated the ability of viral multimeric proteins to induce inflammatory response in macrophages focusing on inflammasome activation. For this study, we have selected structurally diverse viral proteins: filamentous NLPs of measles and mumps viruses (Paramyxoviridae family) that usually cause acute infections, and spherical VLPs derived from VP1 of PyVs (Polyomaviridae family) that generally induce latent viral infections. We also investigated whether viral oligomeric proteins of diverse structure can determine different patterns of cell activation.

NLPs Did Not Activate THP-1 Macrophages Despite Their Interaction With the Cells
First at all, we investigated the uptake of viral proteins by THP-1 macrophages. The cells were treated with recombinant NLPs of measles and mumps viruses, as well as recombinant PyV VLP, for 24 h to observe the uptake of multimeric viral proteins. The NLPs and VLPs were immunostained with the respective monoclonal antibodies. We also immunostained the cells for the macrophage and lysosomal marker CD68. The uptake of both NLPs and VLPs was detected microscopically demonstrating their interaction with THP-1 macrophages ( Figure 1A).
To investigate macrophage activation by multimeric viral proteins, we started with recombinant NLPs of measles and mumps viruses forming filamentous structures. These NLPs are long rod-shaped structures, about 20 nm in diameter, mimicking the nucleocapsids of native viruses (40,41). The NLPs were not cytotoxic to THP-1 macrophages according to propidium iodide (PI) and Hoechst nuclear staining assay (Figures 1B, C) and lactate dehydrogenase (LDH) release ( Figure 1D). Treatment of THP-1 macrophages with NLPs did not cause any inflammatory response according to TNF-a release data ( Figure 1E). The NLPs also did not activate the inflammasome as no change in IL-1b secretion ( Figure 1F) was detected. We concluded that NLPs have no effect on inflammatory response.

PyV-Derived VLPs Induced Inflammatory Response Followed by NLRP3 Inflammasome Activation in THP-1 Macrophages
Next, we investigated macrophage activation by PyV-derived VLPs of spherical structures that are similar to native viruses in their shape and size. In order to evaluate the effects of diverse VLPs, we have selected PyV-derived VLPs of different sizes, ranging 20-60 nm in diameter (36). KIPyV proteins VP1 form heterogeneous VLPs, ranging from 20 to 60 nm, while MCPyV proteins VP1 form more homogeneous VLPs, 45-50 nm in diameter. Thus, we investigated VLPs of different sizes. Treatment of human THP-1 macrophages with PyV-derived VLPs induced cell activation and inflammatory response according to TNF-a and IL-6 release (Figures 2A, B). Therefore, we assumed that VLPs activated the NF-kb signaling pathway, which could be a priming step for further inflammasome activation. In addition, we assayed the cells for anti-inflammatory cytokine IL-10 secretion and did not detect its release after VLP treatment (detected IL-10 values were below the assay detection limit). PyV-derived VLPs induced not only secretion of inflammatory cytokines but also provoked cell death as LDH release was increased ( Figure 2D) and dead cells were detected using PI staining ( Figure 2E). We also demonstrated that VLPs of different sizes induced cell activation signal of different strength. Large-sized (45-50 nm in diameter) homogeneous MCPyV VLPs induced higher inflammatory response as compared to heterogeneous KIPyV VLPs ( Figure 2).
Then, we assessed whether PyV-derived VLPs induce activation of NLRP3 inflammasome in THP-1 macrophages. To test this, we used a small-molecule MCC950, which specifically inhibits NLRP3 inflammasome. PyV-derived VLPs induced cell death (Figures 2D, E) and IL-1b release, showing inflammasome activation ( Figure 2C). Pre-treatment with MCC950 inhibitor before adding the VLPs significantly decreased this activation signal ( Figures 2C-E). The inhibitor MCC950 did not influence the secretion level of other inflammatory cytokine TNF-a showing the specificity of inhibitor to inflammasome activation ( Figure 2A). These data demonstrate that VLPs are potent inflammatory agents triggering inflammatory cytokine secretion and pyroptotic cell death to a similar extent as classical NLRP3 inflammasome activator nigericin (Supplementary Figure 1).
In the next step, we investigated the formation of ASC specks that would indicate the assembly of the inflammasome. To detect ASC specks, we used THP-1-ASC-GFP macrophages. We identified ASC specks after treating cells with PyV-derived VLPs ( Figures 3A, B). However, MCC950 did not fully inhibit ASC speck formation ( Figure 3B), which might be explained by ASC speck release from the inflammasome-activated cells. Secreted ASC specks can promote further maturation of inflammatory cytokines (47). In addition, ASC specks could induce a subsequent inflammatory response by activating the surrounding cells.
Next, we performed time-lapse experiments to observe the dynamics of cell activation by PyV-derived VLPs. We took the microscopy images of THP-1-ASC-GFP within 3-24 h fluorescence of treatment with VLPs. The cell activation started to be seen at 12 h, and after 24 h, ASC specks were evidently observed (Supplementary Figure 2A). We also took the images after 14-24 h treatment with VLPs ( Figure 3D). There were no differences in green intensity and ASC speck formation within this time interval. We measured LDH release in wild-type THP-1 cells treated with VLPs for 15-23 h. Again, we did not find considerable differences within the 15-23 h treatment period ( Figure 3C). We concluded that THP-1 macrophages are activated at about 15 h after VLP addition and their activation pattern remains constant from this time point.
There was a different pattern of cell activation by VLPs and nigericin used as a positive control. In the first hours of nigericin treatment, we have observed caspase-1 expression and ASC speck formation (Supplementary Figures 2B, 3C). However, the number of ASC speck-positive cells decreased in time and no ASC specks were detected 24 h after nigericin treatment. In contrast, VLPs induced the formation of ASC specks after 15 h ( Figure 3D), and the signal did not change within the 15-24 h period of VLP treatment.
To prove VLP-induced activation of the inflammasome, we assayed for active caspase-1, which is a major component of the inflammasome cascade, converting pro-IL-1b to its mature form.
Having the time-lapse experiment data, we collected cell culture supernatant for caspase-1 assay after 15 h-from the steady cell activation time point in order not to lose the active enzyme. The cleaved caspase-1 p20 fragment was detected by WB after THP-1 treatment with VLPs ( Figures 4A, B and Supplementary Figure 3). NLRP3 inhibitor MCC950 reduced generation of the activated caspase-1. Then, we analyzed caspase-1 activation at a single-cell level using FLICA reagent. It contains FAM-YVAD-FMK caspase-1 inhibitor probe, which covalently binds to only activated caspase-1. We revealed caspase-1 activation in some THP-1 macrophages after treatment with VLPs ( Figures 4C, D). Inflammasome inhibitor MCC950 reduced the level of activated caspase-1. We also stained VLP-treated cells with PI to identify dead cells. Activated caspase-1 colocalized with dead cells confirming the pyroptotic cell death and inflammasome activation.

The Mechanism of VLP-Mediated Inflammasome Activation Is Related to Lysosomal Damage
We further investigated the mechanism of inflammasome activation by PyV-derived VLPs. It is likely that accumulation of phagocytosed VLPs in lysosomes can damage them and induce the release of cathepsins, in particular cathepsin B (CtsB), which is one of the activators of NLRP3 inflammasome. We observed a significant decrease of IL-1b secretion after treatment of THP-1 cells with CtsB inhibitor Ca-074 Me ( Figure 5A). However, CtsB did not reduce cell death (Figures 5B, C). The results were the same even at higher (10 μM) Ca-074 Me concentration ( Figures 5D, E). This suggests that VLPs may induce the cytotoxicity by a different mechanism next to lysosomal damage-associated cell death. Then, we treated THP-1 cells with the pan-cathepsins inhibitor K777. We observed a significant decrease in cell death ( Figure 5F) according to LDH assay. Staining of cell nuclei with PI confirmed that the K777 inhibitor prevented cell death ( Figure 5G). Next, we measured IL-1b release and found that K777 significantly suppressed IL-1b secretion ( Figure 5H), but not to control baseline. This indicates that K777 is a partial inhibitor of inflammasome activation induced by VLPs.

PyV-Derived VLPs Induced Inflammasome Activation via K + Efflux
Lysosomal damage can induce further processes, which consolidate inflammasome activation by VLPs as it was shown  with inhaled particles, like silica and polystyrene nanoparticles (5). The VLPs might also trigger different cell activation pathways leading to inflammasome activation. To prove this assumption, we investigated an alternative inflammasome activation pathway using the K + efflux inhibitor, glibenclamide. It significantly blocked cell death ( Figure 6A). We also observed a significant decrease in IL-1b release after glybenclamide pretreatment ( Figure 6B). Taken together, these results reveal the complexity of inflammatory responses leading to inflammasome activation induced by PyV-derived VLPs in THP-1 macrophages. It is likely that several different mechanisms are involved in this process.

Inflammasome Activation by PyV-Derived VLPs Was Confirmed in Primary Human Macrophages
As cell lines may misrepresent real cell activation pattern, we used primary human macrophages derived from monocytes of peripheral blood mononuclear cells to confirm cell activation profile observed in THP-1 cells treated with PyV-derived VLPs. Activation of primary human macrophages with VLPs revealed a similar cell activation pattern to that observed in THP-1 macrophages-a significant increase in activated caspase-1 ( Figures 7A, B), cell death ( Figure 7C), and TNF-a and IL-1b release ( Figures 7D, E). Moreover, NLRP3 inflammasome inhibitor MCC950 significantly reduced primary human macrophage activation and rescued from cell death ( Figures 7A-C). As in experiments with THP-1 macrophages, not all primary human macrophages were activated by the VLPs, and the percentage of dead cells did not exceed 20% (Figures 2E and 7C). According to FLICA assay indicating caspase-1 activation ( Figures 7A, B), only a part of primary human macrophages was activated by VLPs. According to PI staining, primary macrophages having activated caspase-1 also were dead ( Figure 7A) proving the pyroptotic cell death. These experiments indicate a similar pattern of PyV VLPinduced inflammation followed by inflammasome activation both in THP-1 cell line and in primary human macrophages. We found that in some THP-1 cells, the inflammasome is activated, and in some, it is not based on data obtained from caspase-1 and ASC speck formation assays. To identify differences between these cells, we performed single-cell RNA sequencing (ScRNAseq) analysis. This analysis was also used to reveal whether KIPyV and MCPyV VLPs induce different cell activation states since we observed different cell activation levels according to inflammatory cytokines' production.

ScRNAseq Reveals an Overall Conserved Gene Expression Response to PyV-Derived VLPs
In vitro differentiated THP-1 macrophages were subjected to stimulation by VLPs of KIPyV or MCPyV (and unstimulated control-PBS) for 15 h followed by single-cell RNAseq ( Figure 8A). Upon removing transcriptomes with <900 total counts, 32,415 cells were retained, with >10,000 cells per condition. The mitochondrial gene count filter, conventionally used as a signature of dead cells (48), was omitted to enable a comparison of the viability observed from scRNAseq data. Uniform Manifold Approximation and Projection (UMAP) visualization of scRNAseq data revealed an overall similar population structure between KIPyV and MCPyV VLP-treated cells, and an apparent difference compared to the unstimulated control ( Figure 8B). A bulk-like DGE analysis identified 43 upregulated and 12 downregulated genes in KIPyV VLP case vs. the control ( Figure 8C and Supplementary Table 1 genes compared to the control were the same regardless of the VLPs used for stimulation, and the 5 most enriched terms in a GO gene set enrichment analysis of these common genes were associated with immune cell activation processes (migration, chemotaxis, response to TNF, and IL-1) ( Figure 8D and Supplementary Figure 4, Supplementary Table 2).

ScRNAseq Reveals Multiple Subpopulations of THP-1 Macrophages and Changes in Their Abundance Upon VLP Stimulation
The UMAP visualization of scRNAseq data revealed a continuous structure with no clear boundaries between gene expression states. Upon interactive exploration (49) of gene expression patterns and Leiden clustering results at different resolution (i.e., different number of clusters), we chose to report on 9 populations ( Figure 9A) of variable abundance across the 3 conditions ( Figures 9B, C). For each cell population, we identified the 20 most enriched genes (FDR < 0.05, Mann-Whitney U test). Hierarchical clustering of the expression of these genes across all 9 populations revealed both distinct and overlapping gene expression signatures and the transcriptional relationship between populations (Figures 9D, E). The major split of the population dendrogram ( Figure 9D) separates activated states (AI, AII, HA, and IR) enriched after VLP stimulation ( Figure 9F) from SS, RS, MA, and PA. The population of dead cells cluster with the latter but form a distant branch. Consistently with cell viability assays ( Figures 2D, E), the fraction of dead cells, characterized by a high expression of mitochondrial genes ( Figure 9E and Supplementary Figure 5A; Supplementary Table 3), increases with VLP stimulation from 2% to 5%-6% ( Figure 9C). None of the states was completely quiescent, i.e., with all cells in the G0 stage, as suggested by a classification into cell cycle stages based on gene expression, although different populations showed a variable fraction of cells in G1/M/G2/S stages ( Figure 9F).
The cell cluster most uniformly represented across the three conditions (PBS, KIPyV, and MCPyV VLPs) was called steady state (SS) (Figures 9A, F). Among its enriched genes are CD52, LMO4, CHI3L1, and COX5B ( Figure 9E and Supplementary Figure 5B; Table S3). The function of most of the proteins encoded by these genes is unclear. CD52 is a glycoprotein of an unknown function. A recent study showed the ability of soluble CD52 to suppress inflammatory cytokine production by inhibiting TLR-induced NF-kB activation in macrophages and dendritic cells (50). COX5B encodes the 5B subunit of Cytochrome C Oxidase, an essential mitochondrial respiratory chain enzyme (51). An increased expression of COX5B may be related to an increased cellular respiration. LMO4-encoded protein may play a role as an oncogene targeting the TGF-b signaling pathway (52). Overall, the phenotype of SS population is unclear as the functions of most enriched genes in macrophages are unknown. In addition, GO enrichment analysis did not show any relation to known GO terms ( Figure 9E). Therefore, this population could be described as a new possible phenotype of the THP-1 cell line.
One cell cluster, the relative abundance of which changed upon cell activation, was called resting state (RS) ( Figure 9A). This cell population did not express inflammatory molecules, and the enriched genes were related to usual cellular processes. Even so, this cluster distinguished oneself by high enrichment of ARHGAP18, TMEM158, TGM2, OXR, and FN1, which have vital or still unknown functions ( Figure 9E Table 3). ARHGAP18 encodes Rho GTPase Activating Protein 18, which is essential for actin remodeling; thus, it is important for cell migration and controls cell shape (53). TMEM158 encodes Transmembrane Protein 158, whose biological function is unclear, although its involvement in the activation of Ras pathway was shown. It was reported that TMEM158 enhanced proliferation and migration of cancer cells (54,55). For example, another transmembrane protein, TMEM119, was identified as a marker of microglia, representing microglia with homeostatic properties (56,57). Noticeably, TMEM158 expression decreased after activation with VLPs addressing its presence in ramified macrophages (Supplementary Figure 5I). TGM2 encodes Transglutaminase 2 implicated in cell death pathways (58). OXR1 encodes Oxidation Resistance Protein 1, which is involved in protection from oxidative stress and is important for lysosomal function (59,60). FN1 encodes fibronectin 1 known for its function in cell adhesion, motility, and maintenance of cell shape (61,62). Therefore, expression of these genes indicates restful cells.
Another cell cluster that changed upon cell activation was called metabolically active (MA) ( Figure 9A). It had a higher expression of mitochondrial and ribosomal genes, for example, MT-CO3, MT-CYB, RPL37A, and RPS23 ( Figure 9E and Supplementary Figure 5D). The MA cluster did not show expression of genes related to inflammatory cell activation, but the observed gene expression profile indicates increased cellular respiration that is related to high cell metabolic activity (63). Enrichment in ribosomal genes characterized protein production indicating metabolic cell activity.
In the control condition (PBS) we found a cluster in which the cells seemed to be activated a priori. It was named prone to activation (PA) ( Figure 9A). Higher expression of SPP1, CTSC, TREM2, S100A4, and IL1B was detected in some cells of this cluster ( Figure 9E; Supplementary Figure 5E   Supplementary Table 3). It is possible that these cells showed a delayed response to phorbol 12-myristate 13-acetate (PMA) used for THP-1 differentiation. This cell cluster disappeared after treatment with VLPs. CTSC encodes lysosomal protease Cathepsin C. As other cathepsins, CTSC is important for cargo degradation (64). CTSC was also shown to be necessary for activating serine proteases since its absence altered extracellular IL-1b activation (65). In addition, upregulation of CTSC mediates macrophage polarization to inflammatory phenotype (66,67). S100A4 encodes calcium-binding protein S100A4 involved in inflammatory reactions (68,69). SPP1 encodes Secreted Phosphoprotein 1 that is chemotactic, induces IFN-g and IL-12 production, and promotes cell survival (70). TREM2 encodes Triggering Receptor Expressed on Myeloid Cells 2 that triggers secretion of inflammatory molecules (71,72), although the antiinflammatory effect of TREM2 was demonstrated in macrophages lacking TREM2 as toll-like receptor stimulation induced higher pro-inflammatory cytokine secretion (73). Soluble TREM2 was identified as an activator of inflammatory response (74). It enhances phagocytosis as its loss impairs cellular uptake of various substrates, such as cellular debris (75), bacteria (76), or amyloid-beta aggregates (77). Overall, PA cells are differently activated compared to other THP-1 cell populations (AI, AII, HA, and IF), and they distinguish themselves as having inflammation-related and phagocytic cell properties. The cluster of cells formed upon VLP stimulation was defined as a high activation (HA) state ( Figure 9A). This cluster was more abundant in MCPyV VLP-treated cells than in KIPyV VLP-treated cells ( Figure 9F). This cell population was characterized with the highest expression of IL1B and chemokine genes, such as CXCL1, CXCL3, CXCL8, CCL3, CCL20, CCL4L2, and CSTB ( Figure 9E Table 3). It is known that IL1B is expressed at extremely high levels in myeloid-derived cells in response to microbial invasion and tissue injury (78). The product of IL1B is pro-inflammatory cytokine IL-1b, a key mediator of inflammation and one of the main indicators of NLRP3 inflammasome assembly (79). IL-1b induces the production of chemokines and proteases to attract other immune cells to the infection site. Furthermore, secretion of chemokines is reduced in NLRP3-deficient mice demonstrating the importance of inflammasome activation in chemotaxis (80). IL-1b is a product of inflammasome activation; thus, factors stimulating the inflammasome also recruit immune cells. CXCL1, CXCL3, and CXCL8 encode members of the CXC subfamily of chemokines, the chemoattractants for neutrophils (81). CXCL8, also known as IL-8, has no homologs in rats or mice and is a significant component of inflammation-mediated processes as it attracts neutrophils, basophils, and T cells to the site of infection and promotes endothelial cell migration and proliferation (82). HA state is also rich in genes of other chemoattractants CCL3 and CCL20. CCL3, also known as macrophage inflammatory protein 1 alpha, is induced by NF-kB signaling pathway and triggers inflammatory reactions (83). CCL20 acts as a ligand for C-C chemokine receptor CCR6 that induces a strong chemotactic response and plays an important role at skin and mucosal surfaces under homeostatic and inflammatory conditions (84). HA cluster was also highly enriched in CSTB, which encodes cystatin B, a cysteine protease inhibitor known to interact with cathepsin B (85) and considered to protect from cathepsin leakage from damaged lysosomes (86). VLPs. Overall, highly enriched genes of HA state showed a strong inflammatory response in VLP-treated macrophages.
Other cell clusters formed upon cell activation were named the Activation I (AI) and Activation II (AII) states ( Figure 9A). AI cluster was more characteristic for KIPyV than MCPyV VLP-treated cells and AII state, and vice versa ( Figure 9F). Both clusters had enriched inflammation-related genes similarly to HA state but at lower levels ( Figure 9E). AI state was characterized by high IL1B, CXCL8, CCL3L1, and CCL3 expression ( Figure 9E and Supplementary Figure 5F; Supplementary Table 3). AII had relatively low expression of IL1B but the enrichment in CXCL8, CCL3L1, and CCL3 was similar to AI. In general, the same genes were enriched in HA, AI, and AII states, although the expression levels were different. Furthermore, HA and AII were different from AI state as HA and AII had enrichment in STC1, MMP8, ATP2B1, FTH1, and SNX9 genes. Interestingly, enrichment in some of these genes, like ATP2B1 and SNX9, was similar to RS cluster. STC1 is known to encode Stanniocalcin 1, a secreted glycoprotein involved in inflammation and carcinogenesis (87) that possibly serves as a negative mediator of inflammation (88). Interestingly, HA state, which was highly enriched in pro-inflammatory cytokines, also expressed high levels of STC1. MMP8 is a Matrix Metalloproteinase-8 that cleaves collagen, some cell adhesion proteins, growth factors, and chemokines (89), and promotes polarization of macrophages into alternatively activated (M2) macrophages (90). Therefore, AII and HA states describe a group of activated macrophages expressing inflammationrelated genes and a couple of inflammation suppressors. SNX9 encodes Sorting Nexin 9 involved in intracellular trafficking (91) and regulates clathrin-mediated endocytosis (92). It also plays a role in inflammatory reactions (93) and regulation of the micropinocytosis-endocytosis pathway (94). In innate immune cells, this endocytosis pathway may function for the delivery of antigens to their respective intracellular pattern recognition receptors (95). Since the endocytosis pathway of polyomavirus VLPs is unclear, SNX9 may contribute to intracellular recognition of VLPs. ATP2B1 encodes Plasma Membrane Calcium ATPase 1 (PMCA1) important in maintaining cytosolic Ca 2+ for physiological cell functions (96). Ca 2+ mobilization is critical for NLRP3 inflammasome activation (97). Thus, calcium pumps could be implicated in inflammasome signaling detected after VLP treatment. FTH1 encodes the heavy subunit of ferritin, which is essential to store iron inside cells (98). Macrophages are central players in iron metabolism as they recycle senescent erythrocytes and modulate iron availability as part of host protective mechanisms (99). FTH1 is crucial in protection against iron-induced oxidative stress (100). In summary, HA and AII clusters next to inflammation-related genes also expressed higher levels of genes encoding ion channels and proteins involved in endocytosis contrary to AI state.  Figure 6A). A population abundance foldchange analysis emphasized the overall similar phenotype of the two VLP-treatment conditions ( Figure 10F), while differences in abundance relatively to control exceeded 30× ( Figures 10D, E). However, a notable exception was the interferon-response population (IR), highly enriched after MCPyV VLP stimulation compared to KIPyV VLP and absent in control cells ( Figures 10D, E and Supplementary Figure 6B). This suggests differences in pathogen-associated molecular patterns derived from human PyVs. IR population was characterized by the expression of CXCL10, ISG15, IFITIM3, IFIT1-3, IFI6, CXCL11, and other genes ( Figure 9E and Supplementary Figure 5H; Supplementary Table 3), known to be involved in inhibition of viral invasion. We assume that macrophage engagement with viral antigens itself mediates antiviral activity in macrophages without external stimulation, such as lipopolysaccharide.

While Most Population and Gene
In detail, most of IR cluster genes are induced by interferon (IFN) and related to cell response during viral infection. IFITMs encode Interferon-Induced Transmembrane Proteins that limit viral infection via entry pathway (101). In interferon-stimulated cells, IFITM3 is localized on endocytic vesicles that fuse with incoming virus particles and enhance their trafficking to lysosomes (102). A recent study revealed that IFITM3 is able to inhibit SARS-CoV-2 infection in vitro (103). IFITs encode Interferon-Induced Proteins with Tetratricopeptide repeats that block viral infection by interacting with factors responsible for virus replication (104). IFI6 encodes Interferon-a Inducible Protein 6. The exact function of this protein is still unclear. The family of IFI6 genes encodes mitochondrial proteins implicated in inhibition of apoptosis (105). IFI6 was also shown to inhibit DNA replication in hepatitis B virus (106). ISG15 (Interferon Stimulated Gene 15) encodes ubiquitin-like protein that is induced by type I interferons and involved in various processes related to host antiviral response (107). Its upregulation was detected in macrophages after infection by vaccinia and SARS-CoV-2 viruses (108,109). This protein is likely to be an important factor implicated in cytokine storms triggered by SARS-CoV-2. ISG15 deletion in bone marrowderived macrophages induced mitochondrial dysfunction and altered nitric oxide production exposing its importance in regulating mitochondrial function (110). Other genes rich in IR state encode CXCL10 and CXCL11, but the latter expression was detected only in several cells. CXCL10 is known as early interferon response gene encoding CXCL10 with inflammationrelated pleiotropic effects (111). CXCL11 participates in inflammatory reactions and acts as a chemoattractant of activated T cells (112). CXCL10 and CXCL11 are thought to be key mediators of the cytokine storm in immune response to SARS-CoV-2 infection (113). Overall, high expression of IRspecific genes indicates cellular response to interferon induced in THP-1 cells by VLP treatment.

DISCUSSION
In our study, we focused on the inflammatory response of human macrophages treated with viral oligomeric proteins. Cellular response to nanoparticles has been explored previously, covering different mechanisms. It was demonstrated that polymeric nanoparticles induced NLRP3 inflammasome activation dependent on cathepsin B in macrophages (114). Another study showed that oligomeric proteins, Ab fibrils, cause lysosomal damage and induce inflammasome activation (115). In addition, tau oligomers and a-synuclein fibrils were shown to trigger an inflammatory response in microglia (11,116). The latter reports demonstrated inflammasome activation in macrophages by endogenous pathogenic protein oligomers. However, there are limited data on molecular mechanisms and how oligomeric proteins activate NLRP3 inflammasome. To address this question, we investigated the ability of exogenous oligomeric proteins of viral origin to induce inflammatory responses in human macrophages.
In our study, we have used different viral oligomeric proteinsfilamentous NLPs of measles and mumps viruses and spherical VLPs of KIPyV and MCPyV. We did not observe any inflammatory response to NLPs, although NLRP3 inflammasome activation was demonstrated to be induced by N protein of SARS-CoV-2 (117). Interestingly, N protein of SARS-CoV-2 inhibits the cleavage of gasdermin D, which forms pores in the membrane, reducing IL-1b secretion and pyroptosis (118). In our study, only PyV-derived VLPs induced inflammatory response of human macrophages followed by NLRP3 inflammasome activation. Specific NLRP3 inflammasome inhibitor MCC950 blocked the detected activation signal proving that PyV-derived VLPs are a trigger of NLRP3 inflammasome. In addition, VLPs induced secretion of other inflammatory cytokines TNF-a and IL-6, suggesting the engagement of NF-kB signaling pathway in the activated macrophages.
NLRP3 inflammasome can be activated by different mechanisms, such as changes in intracellular ion concentration, mitochondrial, or lysosomal damage followed by cathepsins release (119). First, we studied inflammasome activation mechanism related to lysosomal damage since PyV-derived VLPs are phagocytosed particles. The VLP-induced inflammasome activation signal was significantly reduced by cathepsin B and pan-cathepsin inhibitors. However, the inhibitory effect was not complete compared to control. Cathepsin B inhibitor reduced only IL-1b release and had no effect on cell death. In the case of MCPyVderived VLPs, IL-1b secretion and cell death did not drop to control baseline even using pan-cathepsin inhibitor. In the case of KIPyVderived VLPs, pan-cathepsin inhibitor completely reduced cell death, contrary to IL-1b release. This suggests that VLP size and possibly other structural features of VLPs may define cell activation profile. Viral capsids of KIPyV and MCPyV are structurally different; thus, they may interact with different cellular receptors. In addition, the hemagglutination ability was demonstrated only for MCPyV, indicating different glycoproteins on KIPyV and MCPyV capsid surface (120,121).
There are controversial data on particle-induced inflammasome activation. Some of them show that nanoparticles induce NLRP3 inflammasome activation via phagosomal destabilization (122). Other studies reveal that different nanoparticles induce different mechanisms depending on the composition and structure of the particles (123). For example, cholesterol crystals activate NLRP3 inflammasome; however, inhibitors of cathepsins reduce only IL-1b release and do not change the level of cell death. The role of cathepsins in inflammasome activation is not fully understood. It was recently shown that cathepsins can induce activation of poreforming protein gasdermin D (124). We assumed that PyV-derived VLPs induced inflammasome activation via lysosomal damage and investigated the underlying mechanisms, in particular the mechanism of K + ion efflux. We demonstrated that a specific inhibitor of K + ion efflux significantly reduced cell death and IL-1b release, suggesting the complexity of macrophage activation by PyV-derived VLPs. In addition, studies on the possible outcome of lysosomal leakage indicate that effector molecules such as cathepsins released after the permeabilization of lysosomal membrane may activate the inflammasome by several mechanisms, including K + ion efflux and intracellular Ca 2 + -initiated signaling events (125). Cathepsins can also induce synthesis of IL-1b precursor, leading to a higher release of inflammatory cytokines and involvement of cathepsins in NF-kB signaling (126). This demonstrates that inflammasome activation by phagocytosed particles is beyond the classical mechanism.
Moreover, we have observed a different pattern of cell activation by VLPs and nigericin-a small molecule that promotes K+ efflux and activates NLRP3 inflammasome. In contrast to nigericin that induced caspase-1 expression and ASC speck formation within the first hours of treatment, VLPs induced the formation of ASC specks after 15 h of treatment. These findings might be explained by the fact that nigericin makes pores in cell membranes finally leading to cell lysis, while VLPs induce inflammasome activation via a phagocytosisrelated pathway.
Finally, we proved the inflammatory response observed in human macrophage cell line THP-1 to PyV-derived VLPs using primary human macrophages. The same cellular events were identified in both the THP-1 cell line and the primary human macrophages. Activated caspase-1 was detected in dead primary macrophages, suggesting the pyroptotic cell death. As in THP-1 experiments, only a part of the primary macrophages was activated by the VLPs. Thus, large-sized multimeric protein particles, VLPs, are capable of inducing a similar cell activation pattern in both macrophage cell line and primary human macrophages.
From in vitro studies, it is well known that strong inducers of inflammation, such as lipopolysaccharides, can induce high cell activation states. However, other cell activators usually mediate modest cellular response. For example, in previous research, we showed lower inflammasome activation by Ab oligomers compared to the classical NLRP3 inflammasome inducer nigericin (9). However, when Ab-activated cells were observed microscopically, it was clear that not all of them exhibited similar level of inflammasome activation. It raises a question on the activation state of separate cells. The presence of heterogeneous cell populations was demonstrated several times by other groups. For example, microglia are at different cell activation states depending on the distance from amyloid-beta (Ab) plaque in the brain of Alzheimer's disease patients (127). In addition, microglia are heterogeneous in healthy brain, and cells enriched in inflammation-related genes are present throughout the lifespan and increase in the aged brain (128). ScRNAseq analysis of peripheral blood mononuclear cell culture revealed that monocytes respond differently to lipopolysaccharide (129). One cell population expressed pro-inflammatory alarmins and chemokines to attract immune cells while another was enriched in later inflammation-related and anti-inflammatory genes to either enhance or terminate inflammatory reaction. In our study, IL-1b and TNF-a release data showed overall cell activation by PyV-derived VLPs; however, caspase-1, ASC speck, and cell viability assays allowed the identification of differently activated cells. Therefore, we performed a single-cell analysis to identify gene expression differences in cells, which primarily were expected to be homogeneous-THP-1 cell culture model. ScRNAseq analysis showed a highly divergent cell activation pattern after treatment with PyV-derived VLPs. Four activated THP-1 cell populations varying in expression of inflammationrelated genes, among them IL1B and chemokine genes, were identified. IL-1 family cytokines are known to induce the production of chemokines (130). The stress response and prolonged inflammation also lead to recruitment of immune cells (131). Therefore, increased chemokine expression in cells, enriched in IL-1b and other inflammation mediators, indicates a strong inflammatory response. After VLP treatment, one cell cluster was highly enriched in inflammation-related genes (HA state), including IL-1b. This population was also enriched in CSTB, which encodes cystatin B, an inhibitor of cathepsin B. It possibly copes with cathepsins released after lysosomal leakage (86). Therefore, CSTB expression verifies lysosomal damage induced by VLPs. This cell cluster may represent a population of inflammasome-activated cells, which we revealed according to ASC speck formation, cell viability, and caspase-1 assays.
Interestingly, two populations of activated cells, including the HA cell cluster, also expressed negative mediators of inflammation next to inflammation-related genes. This may indicate a prolonged cell activation state when self-protecting genes were switched on to avoid hyperinflammation. It is known that inflammatory response follows the resolution phase when anti-inflammatory molecules are secreted (132).
Comparing KIPyV and MCPyV VLP-induced cellular response, we found higher expression of interferon responserelated genes after MCPyV VLP treatment. This extends the known differences of KIPyV and MCPyV VLPs. Among interferon-response genes were CXCL10 and IGS15, involved in host antiviral response and found to induce cytokine storm in SARS-CoV-2-infected patients (108,113). Interferoninduced proteins participate in inflammatory reactions and attract different cells of innate and adaptive immunity. Interferon-stimulated genes can be triggered not only by viral antigens directly but also by cellular stress response; thus, their encoded proteins may have a broader biological function (133).
Summarizing, single-cell analysis revealed the presence of heterogeneous cell populations in the in vitro cell culture model, THP-1 macrophages, which might be expected to be a homogeneous cell culture. Although the molecular mechanisms behind different cell activation states are unknown, our study shows that some cells respond to the activating agent by inflammatory reactions while other cells remain unaffected. Assuming the inflammasome activation by protein oligomers as a result of phagocytosis-related process, we suppose that some cells are capable of degrading their cargo and some not, which may lead to lysosomal disruption.
Next to the inflammasome, several other mechanisms are implicated in the response of immune cells to intact viruses or viral components. One example is toll-like receptors (TLRs), such as TLR2 and TLR4, which recognize a variety of microbial molecules, among them viral proteins (134). After TLR activation, NF-kB signaling cascade is initiated to induce expression of pro-inflammatory cytokines. It was demonstrated that polyomavirus VP1 pentamers bind to TLR2 on antigenpresenting cells and induce IL-12 secretion (135). In addition, our results show that PyV-derived VLPs induced NF-kB activation. Interestingly, TLR2 activation is also induced by asynuclein oligomers (136). TLR activation by human papillomavirus-derived VLPs leading to production of proinflammatory cytokines was also demonstrated (137,138). Therefore, VLPs may trigger TLR2-mediated NF-kB activation, followed by pro-inflammatory cytokine production. Similar processes were induced by surface glycoproteins of other viruses, for example, gH/gL and gB of herpes simplex virus (139,140), and gB and gH of human cytomegalovirus (141). Non-structural viral proteins were also demonstrated to activate TLRs. Non-structural protein 1 (NS1) of dengue virus activated TLR4 (142) while both NS3 and capsid proteins of hepatitis C virus activated TLR2, inducing pro-inflammatory cytokine release (143). NLRP3 inflammasome activation in macrophages by hepatitis C capsid protein was also demonstrated (144), suggesting that the same viral protein can mediate several cell signaling cascades.
Inactivated viruses and recombinant viral proteins are broadly used for vaccination (145). However, the mechanisms of the immune response and especially activation of the innate immunity components induced by viral proteins are barely investigated. Therefore, the results of our study on PyVderived VLPs as potent inducers of the inflammatory response in macrophages as compared to filamentous NLPs that are incapable to trigger the inflammation would broaden the understanding of the interaction of viral proteins with innate immune cells. Our study suggests that not all viral proteins can trigger the inflammatory response depending on the structural properties.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are publicly available. This data can be found here: https://www.ncbi.nlm.nih. gov/geo/query/acc.cgi?acc=GSE190666.

AUTHOR CONTRIBUTIONS
AL designed experiments, analyzed the data, and wrote the manuscript. ID and KM contributed to the experiments and manuscript writing. MN was responsible for VLP preparation. AŠperformed single-cell RNA sequencing. RŽwas responsible for single-cell RNA sequencing analysis. IK-K and AG helped to design the experiments and were responsible for manuscript revision. AŽwas responsible for the final manuscript revision. All authors contributed to the article and approved the submitted version.