Mounting data suggest that environmental pollution due to airborne fine particles (AFPs) increases the occurrence and severity of respiratory virus infection in humans. However, it is unclear whether and how interactions with AFPs alter viral infection and distribution. We report synergetic effects between various AFPs and the H1N1 virus, regulated by physicochemical properties of the AFPs. Unlike infection caused by virus alone, AFPs facilitated the internalization of virus through a receptor-independent pathway. Moreover, AFPs promoted the budding and dispersal of progeny virions, likely mediated by lipid rafts in the host plasma membrane. Infected animal models demonstrated that AFPs favored penetration of the H1N1 virus into the distal lung, and its translocation into extrapulmonary organs including the liver, spleen, and kidney, thus causing severe local and systemic disorders. Our findings revealed a key role of AFPs in driving viral infection throughout the respiratory tract an beyond. These insights entail stronger air quality management and air pollution reduction policies.
SIGN UP FOR THE SCIENCE ADVANCES eTOC
Get the latest table of contents from Science Advances delivered right to you!
INTRODUCTION
Exposure to airborne fine particles (AFPs) has been linked to increased morbidity and mortality of respiratory viruses (1), including the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (2, 3). AFPs have also been shown to promote the spread of Streptococcus pneumoniae from the nasopharynx to the lung in mice and enhance the replication of human polyomavirus (4, 5). Furthermore, a few studies to date have identified diverse viral nucleic acids in ambient AFPs, suggesting that AFPs are likely important vehicles for pathogen transmission (6, 7). These observations are of substantial importance, because smaller-sized or ultrafine AFPs are recognized to cross the air-blood barrier and accumulate in distant organs, such as the knee joint and even brai in humans and animals (8–11). Infection of extrapulmonary organs by highly pathogenic influenza A viruses (IAVs), e.g., the H5N1 and H7N9 strains (12, 13), have also been detected, despite the viral exposure and transport pathway warranting detailed investigation.
In the presence of AFPs, viruses could incur more severe infection by disrupting epithelial barriers or altering immune homeostasis (14–16). More recently, AFPs were reported to elevate the expression of the SARS-CoV-2 receptor, angiotensin-converting enzyme 2, and enhance virus entry via the receptor-dependent route (17, 18). While it is known that physicochemical properties such as size of inhaled AFPs can decide where they accumulate along the airway (19), how AFPs interact with viruses and whether the interactions influence viral infectivity, replication, release, and dispersal are still unclear (20, 21). Understanding the impact of AFPs on virus infection and biodistribution could have considerable implications for air ollution and disease management.
In this regard, we aimed to shed light on the synergistic outcomes between H1N1 viruses and different types of AFPs in vitro and in vivo. To gain a general understanding of the AFP properties responsible for respiratory virus infection, we deliberately selected four representative AFPs with diverse sources and compositions in the current study. First, particulate matter with a diameter smaller than 2.5 μm (PM2.5) was sampled from the atmospheric environment. Second, dust particles with silicon dioxide (SiO2) as the main constituent were selected, standing for AFPs originating from sandstorms, mineral dust, and so on (22, 23). Third, biochar particles were prepared through biomass burning, representing AFPs released from biomass fuel consumption, wildfires, and so forth (24, 25). Last, carbon black particles were produced from the combustion of natural gas, representing AFPs emitted from fossil fuel combustion in traffic and industrial activities (26, 27. Our selected AFPs exhibited a broad range of physicochemical properties, including wide variation in surface area and pore volume, differential morphology, and varied surface hydrophobicity. Here, our data unveiled that AFP-virus interactions largely relied on the AFP physicochemical properties, which essentially altered the infection and distribution of loaded viruses. To a differential extent, all AFPs were found to reinforce virus uptake, replication, and release, and to drive viral transport to extrapulmonary organs including the liver, spleen, and kidney. As a consequence, AFPs were uncovered to worsen the severity of viral infection in mice. Together, our findings highlight the remarkable contribution of AFPs to virus transmission and infection, and provide an avenue to further scrutinize the behavior and fate of viruses upon interaction with and loading onto AFPs.
RESULTS
Active association between AFPs and H1N1 viruses
First, we studied the interaction of human influenza viru A/Puerto Rico/8/34 (H1N1 PR8) with naturally sampled PM2.5, silica-based particles (i.e., dust), and carbon-based particles (i.e., biochar and carbon black) (as illustrated in Fig. 1A). Transmission electron microscopy (TEM) and scanning EM (SEM) images showed the size ranges of PM2.5 (300 to 2000 nm), dust (500 to 6000 nm), biochar (200 to 3200 nm), and carbon black (20 to 90 nm) (fig. S1, A and B). As characterized by the Brunauer-Emmett-Teller analysis (fig. S1C), the surface area and micropore volume of AFPs ranked in the order: carbon black > PM2.5> biochar > dust. X-ray photoelectron spectroscopy analysis of surface chemistry demonstrated that PM2.5, dust, and biochar bore relatively abundant oxygen-containing functional groups (C─O and O═C─O), while carbon black contained the least of all studied AFPs (fig. S1, D and E). In addition, the hydrophilicity and hydrophobicity of these particles were determined through the dye partitioning method with a hydrophilic dye (Nile Blue, NB and a hydrophobic dye (Rose Bengal, RB), respectively. Our data suggested that carbon black was the most hydrophobic of the four AFPs, whereas PM2.5, dust, and biochar appeared to be more hydrophilic (fig. S2, A and B).
Fig. 1. AFPs interact with and load H1N1 viruses.
(A) Schematic diagram showing the preparation of virus-laden airborne fine particles (AFPs). Here, AFPs [50 μg/ml of each; particulate matter with a diameter smaller than 2.5 μm (PM2.5), dust, biochar, and carbon black] were incubated with 1.0 × 104 plaque-forming unit (PFU) H1N1 PR8 viruses for 1 hour and washed three times with phosphate-buffered saline (PBS) before examination of morphology, viral load, and infectivity. (B) Transmission electron microscopy (TEM; top panel) and scanning EM (SEM; bottom panel) images show that H1N1 PR8 viruses (denoted by red squares) are adsorbed into the surface of PM2.5, dust, and biochar but are enveloped by carbon black particles. (C) Quantitative real-time polymerase chain reacton (qRT-PCR) analysis of viral loads on AFPs at 20, 50, and 100 μg/ml (n = 3). Viral load for each AFP was obtained by converting cycle threshold (Ct) values using the standard curves (fig. S3). (D and E) Correlation analysis between (D) surface area and (E) pore volume of AFPs and viral loads on AFPs (50 μg/ml). Pearson correlation (R2) and statistical significance (P value) are shown. *P < 0.05.
As seen in the negative-stained TEM images (Fig. 1B), the H1N1 PR8 virus was spherical and contained condensed ribonucleoprotein inside a well-resolved lipid envelope decorated with microspike projections, as reported (28). Even after extensive washing, intact viruses could still be visualized on all AFPs, indicating robust and sustained interaction between H1N1 PR8 viruses and AFPs. On closer inspection, we observed that virus particles were mainly adsorbed into the micropores of PM2.5, dust, and biochar but were wrapped by or embedded between carbon black particles (Fig. 1B), hinting at a dfferent association mode between carbon black and viruses.
To determine the viral loads on different AFPs, we measured the viral matrix protein 1 (M1) gene using quantitative real-time polymerase chain reaction (qRT-PCR), as established (29–31). The viral load for each AFP sample was obtained by converting the cycle threshold (Ct) values (namely, the cycle number at which the PCR product exceeds the set threshold) using the standard curve generated from known virus titers (fig. S3). All AFPs showed concentration-dependent viral loads (P < 0.05; Fig. 1C). Of note, PM2.5 adsorbed much more viruses than the other three AFPs with two- to sevenfold (P < 0.05; Fig. 1C), and biochar ranked second in viral loading, showing approximately two- to fourfold greater viral loads than carbon black and dust, respectively (P < 0.05; Fig. 1C).
Further correlation analyses revealed that viral loading closely correlated to the physicochemical properties of AFPs (Fig. 1, D and E). PM2.5 with higher surfacearea and larger pore volumes displayed higher viral loading than dust and biochar (P < 0.05). Otherwise, consistent with the different patterns of AFP-virus interaction illustrated in Fig. 1B, carbon black was not clustered with the other three AFPs in the correlation analyses (Fig. 1, D and E), implying a discrepant association mode for carbon black that only to a less extent was decided by surface area and pore volume, but likely to a greater extent by other factors including size, morphology, and hydrophobicity. This discrepancy warrants more detailed further investigations. Collectively, these results suggested an active association between AFPs and H1N1 viruses that was regulated by the AFP physicochemical properties, and of them, the area and volume of the contact interface dictated virus loading onto AFPs (including PM2.5, dust, and biochar), except for carbon black with a different particle-virus interaction pattern.
AFP-borne H1N1 viruses remain infectious
To determine whetherviruses loaded onto AFPs remain alive and infectious, we examined the viral hemagglutinin (HA) activity upon interaction with AFPs using the classical hemagglutination assay (32). HA is a glycoprotein on the virus surface that initiates infection, and hemagglutination is observed when the viruses cross-link with red blood cells (RBCs). Similar to H1N1 PR8 viruses (virus control), all AFP-virus complexes agglutinated RBCs (Fig. 2A). Approximately 70 to 90% of the viral HA activity stayed after 4 hours of incubation (Fig. 2B), demonstrating that AFPs had only a mild effect on viral HA activity.
Fig. 2. AFP-borne viruses remain infectious in vitro.
(A) Viral hemagglutinin (HA) activity assay involves coincubation of viruses and standardized concentration of red blood cells (RBCs; top illustration). Experimental data (bottom panel) showed that RBCs were agglutinated after incubation with solo viruses (virus control) or AFP-virus complexes for 4 hours. (B) Relative HA activity of AFP-borneviruses after 4 hours of incubation (n = 5). AFP concentration is 50 μg/ml, and virus titer is 1.0 × 106 PFU. (C) Numbers of plaques formed in Madin-Darby canine kidney (MDCK) cells treated with virus-laden AFPs (20, 50, and 100 μg/ml; n = 3). (D) Representative images of viral plaques. MDCK cells were infected with solo viruses (at 200 PFU) or virus-laden AFPs (50 μg/ml) and stained with crystal violet following a 48-hour infection. Yellow arrowheads denote viral plaques. (E) Schematic showing internalized (37°C for 1 hour), replicated, and released viruses (8, 24, and 48 hours post infection, hpi) in AFP-borne virus-infected target cells. (F) Quantification of internalized, replicated, and released viruses in A549 cells infected with solo viruses (at 200 PFU) through qRT-PCR assay (n = 4). (G to I) Quantification of (G) internalized, (H) replicated, and (I) released viruses in AFP-borne virus-infected A549 cells (n = 4). (J) Flow cytometry data (left) and corresponding quantitative aalysis (right) of A549 cells infected with enhanced green fluorescent protein (EGFP)–tagged viruses bound onto AFPs at 24 hpi (n = 4). SSC-H, side scatter height. (K) Correlation analysis between replicated viruses and EGFP-positive cell percentage in viable cells at 8, 24, and 48 hpi. Pearson correlation (R2) and statistical significance (P value) are shown. *P < 0.05; **P < 0.001.
We further looked at the infectivity of AFP-borne viruses using the standard plaque assay (6). Madin-Darby canine kidney (MDCK) monolayer cells infected with solo viruses (namely, the free virus control) or AFP-borne viruses formed plaques, indicating that AFP-borne viruses remained infectious (Fig. 2, C and D). The number of virus plaques elevated with increasing concentration of virus-laden AFPs (P < 0.05; Fig. 2C). Consistent with the viral loads shown in Fig. 1C, PM2.5-borne viruses gave rise to two- to fivefold more plaques than the other three AFPs at the same concentration (P < 0.05; Fig. 2, C and D) Together, these data qualitatively and quantitatively indicated that viruses bound to AFPs remained infectious.
To verify the infectivity of AFP-borne viruses, we used lung epithelial A549 cells, a sialic acid (SA) receptor-positive cell line commonly used for IAV infection analyses (33), to study the viral lifecycle, including internalization, replication, and release under different conditions (Fig. 2E). Consistent with the IAV infection pattern reported previously (34, 35), the solo viruses showed effective entry into target cells after a 1-hour infection (Fig. 2F). For the single cycle (8 hours post infection, hpi) and multicycle (24 and 48 hpi) virus replication kinetics (31), the levels of intracellular viruses (indicative of newly synthesized virions) increased progressively over time (P < 0.001; Fig. 2F). Moreover, the titers of extracellular viruses (suggestive of released mature virus progeny) displayed a similar incremental trend corresponding to the intracellular virus titrs (P < 0.05). Notably, similar results were demonstrated in primary human small airway epithelial cells (hSAECs; fig. S4A).
When cells were incubated with AFP-borne viruses at 37°C for 1 hour, we found that viruses adsorbed onto AFPs (especially PM2.5) had already entered cells (P < 0.001; Fig. 2G). The contents of intracellular viruses and extracellular viruses quantified after a single cycle and multiple cycles of infection manifested progressive increases (P < 0.05; Fig. 2, H and I). The increase in viral titers revealed that AFP-borne viruses were replicating within A549 cells and releasing mature viral progeny over time. Overall, PM2.5 induced significantly more virus entry, replication, and subsequent release than the other AFPs (P < 0.05; Fig. 2, G to I), in agreement with the greater viral loading as described in Fig. 1C. Consistent findings were also confirmed in hSAECs (fig. S4, B to D).
Flow cytometry analysis of A549 cells infected with enhanced green fluorescent protein (GFP)–tagged viruses or AFP-borne viruses substantiated intracellular accumulation of newly synthesized virions (EGFP-positive cells) in a time-dependent manner (P < 0.001; Fig. 2J and fig. S5). Consistent with qRT-PCR data, flow cytometry data also unveiled that PM2.5 induced a greater proportion of EGFP-positive cells than the other three AFPs with two- to threefold (P < 0.001; Fig. 2J and fig. S5B). The correlation analysis indicated that the percentage of EGFP-positive viable cells was significantly positively correlated with the contents of newly replicated viruses (P < 0.001; Fig. 2K). Together, our results uncovered that AFPs delivered viruses to host cells for effective internalization, replication, and release. PM2.5, which had higher viral loads, consistently elicited greater effects than dust, biochar, and carbon black.
AFPs are implicated in the initial transport of viruses into target cells
To decipher how AFPs participate in the lifecycle of H1N1 viruses, we immunostained he viruses with an anti-HA antibody (red) and observed their interactions with AFPs and target cells using laser scanning confocal microscopy (LSCM). High-magnification LSCM images showed that AFP-borne viruses (red) readily entered cells at sites where AFPs (black) were localized close to the plasma membrane at 1 hpi (yellow arrowheads in Fig. 3A), hinting that AFPs likely acted as a Trojan horse for viral entry into host cells at the early stages of infection. Likewise, similar results were determined in hSAECs (fig. S6). To interrogate whether AFP-borne virus infection is receptor-mediated, we treated A549 cells with neuraminidase (NA) enzyme to remove the surface SA groups that are indispensable for the attachment and entry of IAVs. In agreement with a previous study (36), the removal of SA completely blocked solo virus infection (P < 0.001; Fig. 3B). AFP-borne viruses continued to infect NA-treated cells by 20 to 59% relative to that without NA treatment (P < 0.05), suggesting tha AFPs aided the cell entry of loaded viruses independent of the viral receptors.
Fig. 3. AFPs transport viruses into target cells.
(A) Representative merged images of bright-field and fluorescence (top panel) and confocal immunostaining images (bottom panel) of A549 cells upon incubation with solo viruses (at 200 PFU) or virus-laden AFPs (50 μg/ml) for 1 hour. AFPs close to plasma membrane (yellow arrowheads), viruses (white arrowheads), and cell membranes (white dashed lines) are shown. Viral HA protein is stained with an anti-HA antibody (red), and cell nuclei are stained with 4′,6-diamidino-2-phenylindole (DAPI; blue). (B) Quantification of the replicated viruses in A549 cells with and without neuraminidase (NA) pretreatment following 24 hours of infection with virus-laden AFPs (50 μg/ml; n = 3). (C) Quantification of replicated viruses in NA-pretreated A549 cells followed by a 24-hour infection with virus-laden AFPs (50 μg/ml) in the presence of different inhibitors (n = 3). (D) Shematic showing that virus-laden AFPs are mainly internalized through clathrin-dependent endocytosis, followed by lipid raft–mediated endocytosis. *P < 0.05; **P < 0.001, relative to untreated control, or as indicated.
To further investigate the entry mechanisms of AFP-borne viruses, we used various inhibition modes to hinder the possible cellular uptake pathways, as established previously (37, 38). The chlorpromazine treatment decreased by 76% for PM2.5-borne, 51% for dust-borne, 57% for biochar-borne, and 74% for carbon black–borne virus infection, compared to cells without inhibitor treatment (P < 0.05; Fig. 3C), suggesting that AFP-borne viruses were internalized primarily by clathrin-mediated endocytosis. Methyl-β-cyclodextrin also significantly repressed viral infection with a weaker effect than chlorpromazine (P < 0.05; Fig. 3C), while cytochalasin D elicited negligible reductions, indicating that virus-laden AFPs were also taken up through lipid raft–mediated endocytosis but no to any degree via macropinocytosis. Carbon black may rely a little more on the micropinocytosis uptake pathway than the other three AFPs (P < 0.05; Fig. 3, C and D).
AFPs favor the budding and release of progeny virus
After multiple cycles of viral replication, LSCM imaging analysis revealed that the progeny virions were randomly distributed in the cytoplasm of cells infected with solo viruses, as reported (31, 34). Differently in cells infected with AFP-borne viruses, the progeny virions tended to accumulate around AFP localized sites (white arrowheads in Fig. 4A). Lipid rafts have been proved to serve as the preferred sites for viral assembly and budding at the plasma membrane, including for IAVs (39), Ebola virus (40), herpesviruses (41), human respiratory syncytial virus (42), and human immunodeficiency virus (43). Thus, we performed immunofluorescent staining using an anti-human flotillin 1 antibody (green) to visualize lipid rafts in target cells. The LSCM images showed that lipd rafts were enriched near the AFP-attached plasma membrane and were colocalized with viral HA proteins (Fig. 4A). AFPs, especially the hydrophobic carbon black (fig. S2), bearing a high affinity toward lipids, attached to the plasma membrane and formed lipid raft–enriched domains, whereby the progeny virions were postulated to take advantage of the lipid rafts for apical plasma membrane targeting and budding (Fig. 4B). The mechanism underlying AFP-aided progeny virion budding was also verified in hSAECs (fig. S7).
Fig. 4. AFPs favor progeny virion budding.
(A) Representative merged images of bright-field and fluorescence (top panel) and confocal immunostaining images (bottom panel) of A549 cells upon incubation with solo viruses (at 200 PFU) or virus-laden AFPs (50 μg/ml) for 24 hours. Intracellular virions are stained with an anti-HA antibody (red). Lipid rafts are stained with an anti-flotillin 1 antibody (green). Cell nuclei are stained with DAPI (blue). Yellow arrowheads denote te location of AFPs, and white arrowheads indicate intracellular virions. Cells are indicated through dashed lines. (B) Schematic illustrating AFPs in directing the budding and release of progeny virions. (C) Relative ratio of extracellular virus counts to intracellular virus counts at 24 hpi (n = 4). Virus concentration in virus control is the same as that loaded onto each AFP. Data are presented as percent changes relative to the according virus control. (D) Representative merged images of bright-field and fluorescence (top panel) and confocal fluorescence images (bottom panel) of A549 cells treated with AFPs (20 μg/ml) for 24 hours. Red fluorescence (Evans blue) is indicative of the injuries of cell plasma membrane. Cell nuclei are stained with DAPI in blue. *P < 0.05; **P < 0.001, relative to virus control.
Furthermore, when the ratio of extracellular to intracellular virus counts was evaluated, we found that PM2.5, dust, biochar, and carbon black enhanced the release of viral progey by 14, 23, 87, and 174%, respectively, compared to their according virus control groups (P < 0.05; Fig. 4C). The larger difference for carbon black compared to the other three AFPs (P < 0.05) was consistent with the greater mass of carbon black adsorbed onto the cell plasma membrane, as described in Figs. 3A and 4A. Of note, particles with greater hydrophobicity exhibit higher adsorption capacity onto the plasma membrane (44, 45), and carbon black was the most hydrophobic of the four AFPs (fig. S2), which fundamentally accounted for the highest particle-membrane affinity. Moreover, the injury of host cell membrane has been reported to enhance the release of hepatitis A virus (46), poliovirus (47), and dengue virus (48). On the basis of the principle that dyes tend to accumulate in cells with reduced membrane integrity (49), we looked into the cell plasma membrane injuries upon AFP exposure. In analogy to the above results, all AFPs, particularly carbon black, yielded marked red fluorscence with a three- to sevenfold increase relative to untreated control, indicating active association of these particles with and resultant damage to the plasma membrane (P < 0.001; Fig. 4D and fig. S8). Together, our results indicated that AFPs reinforced the intracellular delivery of H1N1 viruses, and the particle-plasma membrane association even favored the budding and release of progeny virions.
AFPs alter H1N1 distribution along the respiratory tract
Afterwards, we further examined the infection of AFP-borne viruses in vivo by administering BALB/c mice with virus-laden AFPs (200 μg/kg body weight) through the intranasal route. This exposure dose of AFPs is equal to about 12 μg/m3 in ambient air, which is lower than the exposure levels of inhaled particles reported in most studies on the health implications of air pollution (50, 51). Viral titers in the nasal cavity, trachea, and lung corresponding to the upper, middle, and lower respiratory tracts were assessed at 1, 3, and 8 das post infection (dpi) (Fig. 5A). As shown in Fig. 5B, viruses were evenly distributed in the nasal cavity and lung of mice infected with solo viruses at 1 dpi, as reported previously (12, 29). This finding is supported by the spatial expression pattern of specific virus receptors, such as α-2,6 SA-linked receptors, which are more abundant in the lung and nasal cavity than those in the trachea (52, 53). Unlike solo viruses, 82% of PM2.5-borne, 98% of dust-borne, and 76% of biochar-borne viruses that were administered ended up in the lung (Fig. 5B). In contrast, 83% of carbon black–borne viruses still remained in the nasal cavity, with only a small amount reaching the lung at 1 dpi. A similar biodistribution pattern along the airway was identified for Au-doped AFPs that were tracked using inductively coupled plasma mass spectrometry (Fig. 5C). These results collectively uncovered that different AFP types yielded different viral biodistribution.
Fig. 5. AFPs alter the distribution of viuses along the respiratory tract.
(A) Timeline and experimental design for evaluating viral titers in the nasal cavity, trachea, and lung in BALB/c mice that were intranasally infected with virus-laden AFPs (200 μg/kg body weight). i.n., intranasal. (B) Percentage of viruses detected in the nasal cavity (red), trachea (green), and lung (blue) through SYBR green–based qRT-PCR analysis at 1, 3, and 8 days post infection (dpi; n = 5). Data are presented as the proportion of virus titers at each site relative to the whole titers in the respiratory tract. (C) Inductively coupled plasma mass spectrometry analysis of Au content in the nasal cavity (red), trachea (green), and lung (blue) of mice intranasally administrated with 1.0 mg/kg body weight Au-doped AFPs for 2 hours (n = 4). Data are presented as the proportion of Au content at each site relative to the whole mass in the respiratory tract. (D) Representative nasal cavity sections from mice treated with 1.0 mg/kg body weight virus-ladenAFPs for 2 hours. Inserts show magnified views of smaller red squares. Blank control is PBS only. (E) Representative hematoxylin and eosin (H&E)–stained lung sections from mice treated with PBS (blank control), solo viruses, or virus-laden AFPs (200 μg/kg body weight) at 1 dpi. AFP accumulation (black arrowheads) is observed. (F) Representative confocal images of lung sections obtained from mice as described in (E). Viruses are stained with an anti-HA antibody (green), and cell nuclei are stained with DAPI (blue).
In the following, histological examination revealed that particles were visualized to accumulate in the turbinate of nasal cavity from mice upon exposure to virus-laden carbon black (Fig. 5D), but were found in the lung alveoli after exposure to virus-laden PM2.5, dust, and biochar at 1 dpi (black arrowheads in Fig. 5E). Analogously, no particles were readily visualized in the lung after 1 day of exposure to virus-laden carbon black (Fig. 5E), in support of the above results Fig. 5, B to D). To verify these findings, the virus localization in the lung was also qualitatively and quantitatively corroborated through in situ immunofluorescent staining using an anti-HA antibody (Fig. 5F), which correlated with the virus distribution pattern (Fig. 5B). Among the AFP physicochemical properties, the surface hydrophobicity essentially determined the interaction with and adsorption of AFPs to the plasma membrane (44, 45). Particles with greater hydrophobicity are more easily entrapped and immobilized in the respiratory mucus (54). To this end, carbon black particles were the most hydrophobic of the four AFPs (fig. S2), which contributed to understanding a differential distribution pattern of carbon black along the respiratory tract in comparison to the other three AFPs.
Similar to the reported features of IAV replication (55, 56), viral contents in the lung for solo viruses and all AFP-virus complexes increased rapidly from 1 to 3 dpi before declining steadily betwen 3 and 8 dpi (Fig. 5F and fig. S9, A and B). Nevertheless, at the terminal point of infection (8 dpi), abundant viruses still remained in the lung (fig. S9B), consistent with the robust replication ability of IAVs in the lower respiratory tract (53). Of note, AFP-borne virus infection still induced greater proportion of virus distribution in the lung than that in the virus control group at 3 and 8 dpi (Fig. 5B). Consequently, sustained inflammatory cell infiltration and histological damage were visualized (fig. S9, A and B).
AFPs function to spread the H1N1 virus to extrapulmonary organs
Although IAV infections largely attack and produce adverse outcomes along the respiratory tract (53, 55, 56), previous studies have shown that respiratory viruses can also reach distant organs and even cause lethal effects (12, 13). Furthermore, we recently unveiled that 125I-radiolabeled AFP particles could enter the circulation from the lung and thereafter accumulate in the liver, spleen, kidney, kne joint, and even the brain (8, 9). Inspired by these findings, we investigated whether AFP-virus complexes also spread to extrapulmonary sites. To realize this, H1N1 PR8 viruses labeled with an anti-HA antibody (green) and red fluorescent SiO2 particles with a similar composition to dust were used.
Whole-mount sections of the lung revealed that solo viruses tended to infect the tracheal epithelial cells (white arrowheads in Fig. 6A) at 1 dpi, as reported previously (56). On their own, SiO2 particles translocated from distal bronchioles into the alveolar parenchyma (Fig. 6A, middle panel), in agreement with our previous reports (8, 9). This intrusion into the interstitial compartment suggested that SiO2 particles could traverse the blood-air barrier (8, 9). Notably, with virus-laden SiO2 particles, we found virus signals (green) overlapping with particle signals (red) in the interstitial compartments (yellow arrowheads in Fig. 6A) and notably also in distant organs including the liver,spleen, and kidney (Fig. 6B). Constant findings were also verified in polystyrene particle–borne virus infection (fig. S10A). In addition, virus localization was validated in these organs from mice upon infections of all AFP-borne viruses using immunofluorescent staining of viral HA protein (fig. S10B). Collectively, these results substantiated that virus-laden AFPs could traverse the blood-air barrier and be transported to extrapulmonary organs through the blood circulation.
Fig. 6. AFPs drive the spread of viruses into extrapulmonary organs.
(A) Representative confocal images of lung sections obtained from mice treated with solo viruses (top panel), SiO2 particles (middle panel), or virus-laden SiO2 particles (200 μg/kg body weight; bottom panel) for 1 day. SiO2 control group received the same mass of particles as the virus-laden SiO2 group. Virus concentration in virus control is equal to that loaded onto SiO2 particles. (B) Representative immunofluorescent staining of the liver, sleen, and kidney sections obtained from mice treated with virus-laden SiO2 particles for 3 days. The right panel indicates magnified views of red squares. Viruses are stained with an anti-HA antibody (green), and cell nuclei are stained with DAPI (blue). (C) Relative changes of virus titers in the blood and distant organs (liver, spleen, and kidney) from mice infected with virus-laden AFPs (200 μg/kg body weight) or solo viruses at 1, 3, and 8 dpi (n = 5). Virus titers are quantified through TaqMan qRT-PCR and presented as percent changes relative to the corresponding virus control at 1 dpi. Ct value >38 was considered nondetectable. (D) Proportion of virus-positive organs in animals infected with virus-laden AFPs (200 μg/kg body weight) in contrast to the same titers of solo viruses at 1 and 3 dpi, as determined by TaqMan qRT-PCR (n = 10). *P < 0.05; **P < 0.001, relative to the corresponding virus control at 1 dpi, or as indicated.
To elaborate on how AFPs drive the virus spread follwing infection, we quantified the viral titers in the blood and distant organs using a more sensitive method via the determination of viral polymerase acidic (PA) protein gene by TaqMan qRT-PCR (fig. S10C), as previously reported (29). A quick and notable increase of the viral titers was observed in the blood and extrapulmonary organs (liver, spleen, and kidney) following the infection with virus-laden PM2.5, dust, and biochar at 1 dpi (P < 0.05; Fig. 6C). Except for carbon black, it remained in the upper airway and was slower in driving viral extrapulmonary transport (Fig. 6C). Because of the viral transport and replication, virus titers at each site increased rapidly at 3 dpi for AFP-borne virus infection (P < 0.05; Fig. 6C). AFP-borne virus infection showed higher proportions (up to threefold) of virus-positive organs than the virus control at 1 and 3 dpi (Fig. 6D). The viral titers markedly dropped to nearly nondetectable levels in the blood at 8 dpi for all infections (P < 0.001) ut still remained at high levels in the extrapulmonary organs (Fig. 6C). Of note, no virus was detected in the heart and brain in this experiment. These findings suggested that AFPs aided the transport of viruses to extrapulmonary organs.
AFP-borne viruses cause more severe phenotypes than AFPs alone or solo viruses in mice
With these encouraging findings, we looked further into the adverse outcomes both systemically and locally in mice challenged by diverse infections. Neutrophils are rapidly mobilized from the bone marrow upon viral infection, leading to an increase in the peripheral circulation (57). In the presence of inflammatory signals, circulating lymphocytes are recruited to the infected respiratory tract, resulting in lower peripheral counts (58). Thereby, as a classical parameter to assess the inflammatory state in virus infection, the ratio of neutrophil/lymphocyte has been reported to indicate the occurrence, development, and outcomes of infectious diseases (59, 60). In ageement with the accumulation of inflammatory cells in the lung (Fig. 5E and fig. S9), complete blood count analysis validated that PM2.5-, dust-, and biochar-borne viruses induced 44 to 105% increase in the neutrophil counts and 5 to 7% decrease in the lymphocyte counts in the peripheral blood, compared to untreated blank control at 1 dpi (P < 0.05; Fig. 7, A to C). Otherwise, carbon black–borne viruses did not induce substantial changes in these parameters early on in the infection (Fig. 7D).
Fig. 7. Virus-laden AFPs cause greater systemic immune responses in vivo than virus infection alone.
Relative changes of neutrophil counts (left panel), lymphocyte counts (middle panel), and the ratio of neutrophil/lymphocyte (right panel) in the peripheral blood of mice treated with (A) PM2.5-, (B) dust-, (C) biochar-, and (D) carbon black–bearing viruses at 200 μg/kg body weight or the same titers of solo viruses at 1, 3, and 8 dpi (n = 5). Change was calculated relative to blank control (PBS) (E) Relative changes in body weight of mice treated with PBS (blank control, green line), AFPs alone (orange line), virus control (purple line), and virus-laden AFPs (red line) over the 8-day infection period. Change in body weight was compared with initial weight (n = 5 mice per group). Virus titers in virus control are the same as those loaded onto according AFPs. *P < 0.05; **P < 0.001, relative to blank control, or as indicated.
During the later stage of infection, i.e., at 3 and 8 dpi, AFP-borne viruses and solo viruses, especially the former, induced systemic inflammation. We observed a significant increase in neutrophil counts and a significant drop in lymphocyte counts in the AFP-borne virus treatments (P < 0.001), while changes in the peripheral blood were marginal in mice treated with AFPs only (Fig. 7, A to D). Consistently, systemic inflammation was more pronounced in PM2.5-borne virus infection than in the other three AFP-borne virus infections (P < 0.05; Fig. 7, A to D).Moreover, AFP-borne virus infection resulted in a 14 to 65% increase in the ratio of neutrophil/lymphocyte compared to virus control, indicating that the AFP-virus complex elicited more severe inflammation (P < 0.05; Fig. 7, A to D). Notably, because of the enhanced transport of AFP-virus complexes to distant organs, local inflammation was observed in the hematoxylin and eosin (H&E)–stained images of the liver, spleen, and kidney (figs. S11 to S13). Together, these data revealed that AFPs markedly altered biodistribution and provoked inflammatory responses to the viruses in the lung and other sites.
In addition to the observed inflammation, all virus infections caused a marked reduction in body weight, which is one of the common symptoms after IAV infection and is believed to be caused by multiple pathologies, including fever and decreased food/water intake (61, 62). PM2.5-, dust-, and biochar-borne viruses afflicted the infected animals with quick body weight loss at the early stage o infection, with PM2.5 quickly showing the greatest weight loss (P < 0.001; Fig. 7E). However, carbon black–borne viruses displayed delayed body weight loss during 1 to 3 dpi (Fig. 7E). This difference should be ascribed to the retention of carbon black particles in the upper respiratory tract and slower transport into the lung. Overall, our data manifested that AFP-borne virus infections (to a divergent degree) resulted in more marked body weight reduction than their according virus control (Fig. 7E). The more severe phenotypes in animals induced by AFP-borne viruses in comparison to virus alone should reside in (at least partially) the AFP-aided transportation of viruses deep into the lung and on to distant organs, as well as the synergistic toxicities by AFPs and viruses. Together, our results unearthed that, besides acting as a carrier, AFPs could markedly drive the biodistribution and inflammatory responses of the viruses along the respiratory tract and afterwards to the distant stes, inducing detrimental effects that were different from (even more severe than) those produced by AFPs alone or solo viruses.
DISCUSSION
The implication and contribution of airborne particles in pathogen spread and infection have been studied for years. Under realistic exposure scenarios, respiratory viruses are reported to be mostly transported by aerosols and AFPs (6, 7). However, exactly how AFPs interact with viruses and involve in the viral infection process is largely unknown. The previous research mainly focused on the indirect effects of AFPs on viral infection, such as the damage of defense networks against pathogens or the promotion of viral receptor expression (14–18). Despite these progresses, it is crucial to understand the impacts of AFPs per se on virus loading, as well as the influence of AFPs on the viral infection process, including viral entry, replication, and release. Here, four typical AFPs with diverse sources and compositions were selected, comprising naturaly sampled PM2.5, silica-based dust, biomass-derived biochar, and fossil fuel–derived carbon black. Under influence of various physicochemical properties, the AFP-virus interactions markedly affected viral lifecycle and tissue distribution and ultimately elicited divergent pathogenesis and outcomes of respiratory infection. For example, AFPs (especially the hydrophobic carbon black) promoted the formation of lipid raft–enriched fractions on the plasma membrane, where the progeny virions targeted and took advantage of the lipid rafts for budding.
The translocation and biodistribution of AFPs have been verified to dictate their toxic effects (10, 11). Recent studies from our and other groups have continued to elaborate on the extrapulmonary translocation of exogenous fine particles (8–11). As evidenced by previous reports, particles with fine size could escape from the body’s protection system and reach deep into the alveolus regions of the lung, and then damage and penetrate the air-bloo barrier (8, 9). The particles that translocated into the blood circulation could distribute to other secondary organs, including the liver, spleen, and kidney (8–11). Inhaled AFPs have been found in diverse human tissues, such as the blood, heart, and even brain, highlighting the pathway of particle translocation from the bloodstream to distant organs (63, 64). Likewise, not restricted to the respiratory tract, IAV-associated extrapulmonary infection has been identified, particularly for infection with highly pathogenic virus strains (12, 13). In this context, AFPs provide intriguing possibilities to alter the viral biodistribution and the occurrence and progression of IAV infection. AFPs could act as a Trojan horse to load viruses and transport them into distant organs. The relatively easier localization of AFP-virus complexes in organs with large-volume blood suggested the tendentious accumulation of particles, as supported by previous findings (8, 9). Because of the detrimental heath effects, AFP-induced physical damage and immune disorders could reinforce the extrapulmonary translocation of viruses and also hinder viral clearance (65–67). Besides, RBCs could be a vehicle to transport both viruses and particles as well (68, 69), which was indicated as another potential mechanism in driving the persistent delivery of viruses to distant organs, especially the spleen as the graveyard of RBCs.
In summary, our study uncovered active interactions between AFPs and the H1N1 PR8 virus, which were largely regulated by the physicochemical properties of the AFPs. On the one hand, the surface area and pore volume of the AFPs were the main determinants in governing their virus loading capacity, except for carbon black with a different particle-virus interaction pattern. On the other hand, the surface hydrophobicity of AFPs determined the particle-plasma membrane interactions, resulting in enhanced virus budding and release in target cells and altered viral distribution along he respiratory tract in animals. Notably, we found that AFPs provided a receptor-independent path for viral entry into and budding from host cells. Moreover, these particles drove the loaded viruses deep into the lower respiratory tract and on to distant organs, such as the liver, spleen, and kidney. As a result, AFP-borne viruses (in particular, PM2.5-borne viruses) yielded more severe infection than the virus alone, characterized by local tissue injuries, systemic inflammation, and body weight loss. The overall findings are summarized in Fig. 8. Our results offer detailed insights into the impact of air pollution on respiratory infection, providing further urgency for prioritizing air quality management and air pollution reduction policies.
Fig. 8. Schematic illustration of the interaction between AFPs and H1N1 virus and the spread of the virus driven by AFPs.
AFPs such as PM2.5, dust, biochar, and carbon black adsorb different amounts of viruses to form infectious AFP-virus complexs (top panel). AFP-virus complexes infect the host cells through a receptor-independent pathway, and AFPs direct the internalization and budding of the viruses (middle panel). Different AFP-borne viruses are distributed at different locations along the respiratory tract and transported into distant extrapulmonary organs (bottom panel).
MATERIALS AND METHODS
Ethics statement
All animal experimental designs and protocols were approved by the Animal Ethics Committee at the Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences. The ethical approval number is AEWC-RCEES-2020001. All virus infection protocols were conducted in strict accordance with the standards of the Institute of Microbiology, Chinese Academy of Sciences.
Viral amplification and quantification
The mouse-adapted human H1N1 PR8 virus and the recombinant H1N1 PR8 virus with NA gene segment harboring EGFP gene (EGFP-PR8 virus) were respectively propagated and purified, as described previously (70). Viraltiters were assayed through the plaque assay.
H1N1 PR8 virus and EGFP-PR8 virus were quantified by qRT-PCR analysis. In brief, viral genomic RNAs were extracted via the Viral RNA Rapid Extraction Kit (Zomanbio, China) according to the manufacturer’s instructions. Extracted RNAs were reversely transcribed into cDNAs using Script RT All-in-One Mix (StarLighter, China). Gene expression levels were determined using an Mx3005P qRT-PCR instrument (Bio-Rad, USA). The primer sequences are listed in table S1.
The standard curves for viral quantitation were generated by a 10-fold serial dilution of virus stocks against the Ct values, with a range from 1.0 × 101 to 1.0 × 107 plaque-forming unit (PFU) viruses for SYBR green–based qRT-PCR analysis, and 1.0 × 101.7 to 1.0 × 106 PFU viruses for TaqMan qRT-PCR analysis. The viral load of each sample was converted using Ct values and the standard curves, as described previously (29–31).
Preparation and characterization of AFPs and virus-laden AFPs
PM2. samples were collected from Beijing, China. PM2.5 specimens loaded onto the polytetrafluoroethylene filter membranes were eluted in ultrapure water by ultrasonic extraction and then freeze-dried. Dust samples were purchased from the Powder Technology Inc. (USA). Biochar samples were prepared via pyrolyzing biomass from corn straw under oxygen-limited conditions at 500°C, as described previously (25). The biochar powder was crushed and passed through a 3-μm sieve. Carbon black samples were purchased from the Degussa Inc. (USA). Red fluorescent SiO2 particles (with a particle size of 2 μm) and red fluorescent polystyrene particles (with a particle size 2 μm) were purchased from the Baseline Inc. (China). The synthesis and characterization of Au-doped AFPs were performed, as described previously (8, 9).
To prepare virus-laden AFPs, AFPs (20, 50, and 100 μg/ml) were coincubated with 1.0 × 104 PFU viruses for 1 hour at room temperature. The mixtures were washed three times with cold phosphte-buffered saline (PBS) by centrifugation at 1000g for 6 min to remove any free viruses. The loaded viruses on AFPs were quantified through the SYBR green–based qRT-PCR analysis.
The morphological characteristics of AFPs and virus-laden AFPs were determined by TEM (H-7500, Hitachi, Japan) and SEM (SU-8020, Hitachi, Japan). The size distribution of AFPs was assessed by measuring more than 300 particles in TEM images. The surface area and micropore volume of AFPs were measured via a surface area and porosity analyzer (ASAP2460, Micrometrics, USA). The surface functional groups of AFPs were quantified using an x-ray photoelectron spectrometer (ESCALAB 250Xi, Thermo Fisher Scientific, USA).
The hydrophilicity and hydrophobicity of AFPs were determined using the dye partitioning method. In brief, a series of AFP suspensions at different concentrations were coincubated with NB (20 μg/ml; Sigma-Aldrich, USA) or RB (20 μg/ml; Sigma-Aldrich, USA) solution for 3 hours at room temperature. Aftercentrifugation at 16,000g for 1 hour, the NB and RB molecules in the supernatants were recorded at 620 and 549 nm with an ultraviolet-visible spectrometer (Epoch, BioTek, USA). The partitioning quotient (PQ) was calculated by the ratio of the dye bound onto the particle surface (Dbound) to free dye molecules in the liquid phase (Dfree); i.e., PQ = Dbound/Dfree. The linear relationship between PQ and AFP surface area was analyzed. The total surface area of the particles dispersed in the suspension was calculated from the hydrodynamic size of these particles determined by the Zetasizer Nano ZS (Malvern Panalytical, UK) with the assumption of a spherical shape, as previously described (71).
HA test
To test the viral HA activity, AFPs (50 μg/ml) were first coincubated with 1.0 × 106 PFU viruses for the designed times at room temperature. Then, 50 μl of mixtures were incubated with equal volumes of 2% chicken erythrocytes (Beijing Merial Vital Laboratory Animal Technology Co. Ltd., China) fr 30 min at room temperature in the V bottom of 96-well plates. High HA activity occurs when viral HA protein cross-links with RBCs and forms a diffuse reddish solution. RBCs precipitate at the bottom of the well when few or no viruses bind to them. Agglutinated erythrocytes were quantified by examining the absorbance of hemoglobin at 575 nm using a plate reader (Thermo Fisher Scientific, USA). The same amounts of viruses were used as the virus control.
Cell culture
A549 cells and MDCK cells were purchased from the Shanghai Cell Bank of Type Culture Collection (China) and were cultured in Dulbecco’s modified Eagle’s medium (DMEM) (Corning, USA) supplemented with 10% fetal bovine serum (NEWZERUM, New Zealand) and 1% penicillin/streptomycin (Gibco, USA) in a cell incubator with 5% CO2 and at 37°C. hSAECs were purchased from the iCell Bioscience Inc. (China) and cultured in an airway epithelial cell basal medium according to the manufacturer’s protocol.
Cell membrane integrity assay
To asay the membrane permeability, cells were treated with AFPs (20 μg/ml) for 24 hours and stained with a cell membrane integrity assay kit, following the manufacturer’s protocol (BestBio Biotechnology Co. Ltd., China).
Plaque assay
To test the titers of virus stock, MDCK cells in confluent monolayers were infected with 10-fold serially diluted viruses. To evaluate the infectivity of viruses, cells were incubated with virus-laden AFPs (20, 50, and 100 μg/ml) or solo viruses. After the infection at 37°C for 1 hour, cells were overlaid with DMEM supplemented with 1% low melting agarose (GenStar, China), 1% bovine serum albumin (Solarbio, China) and [N-tosyl-l-phenylalanyl chloromethyl ketone]–treated trypsin (1 μg/ml; Sigma-Aldrich, USA). At 48 hpi, cells were stained with 0.25% crystal violet solution (Solarbio, China), following fixing with 4% paraformaldehyde (Biorigin, China). Viral titers or infectivity were determined by the number of plaques.
Viral infection of pulmonary epithelial cels
For viral internalization assay, target cells were incubated with virus-laden AFPs (50 μg/ml) or solo viruses at 37°C for 1 hour. To investigate viral replication and release, cells were further cultured for 8, 24, and 48 hours. To verify the receptor-independent infection pattern of AFP-borne viruses, cells were treated with NA enzyme (1 U/ml; Solarbio, China) at 37°C for 2 hours, followed by infection for 24 hours.
To investigate the uptake mechanism of virus-laden AFPs, cells were pretreated with NA enzyme (1 U/ml) at 37°C for 2 hours and immediately incubated with diverse inhibitors in complete cell culture medium at 37°C for 1 hour, followed by infection with virus-laden AFPs (50 μg/ml) for 24 hours. Chlorpromazine (5 μg/ml; Solarbio, China) was used to inhibit clathrin-mediated endocytosis; methyl-β-cyclodextrin (20 μg/ml; Macklin Biochemical Co. Ltd., China) was used to inhibit lipid raft–mediated endocytosis; amiloride (50 μg/ml; Solarbio, China) was used to block micropinoctosis; cytochalasin D (2.5 μg/ml; Aladdin Biochemical Technology Co. Ltd., China) was used to inhibit actin polymerization and macropinocytosis, as established previously (37, 38).
Intracellular viruses were quantified via SYBR green–based qRT-PCR analysis following the extraction of total cellular RNAs, and progeny viruses in the corresponding culture medium were quantified according to the procedure of viral quantification.
Flow cytometry analysis
Target cells were infected with EGFP-PR8 virus-laden AFPs or solo viruses for 8, 24, and 48 hours. Afterwards, live cells were distinguished using the Fixable Viability Dye 780 (Thermo Fisher Scientific, USA) staining assay. The percentage of EGFP-positive cells within gated Fixable Viability Dye–negative cells was determined with Attune NxT Flow Cytometer (Thermo Fisher Scientific, USA).
Immunofluorescence assay
Target cells were infected with AFP-borne viruses or solo viruses and then were extensively washed with PBS and fixed with methanl for 15 min at −20°C. After washing, these cell specimens were blocked with blocking buffer (PBS containing 5% goat serum and 0.3% Triton X-100) for 1 hour at room temperature, and further incubated with mouse monoclonal anti-HA antibody (1:200; Abcam, UK) and/or rabbit monoclonal anti-flotillin 1 antibody (1:500; Abcam, UK) at 4°C overnight. Cells were sequentially stained with CoraLite 594–conjugated goat anti-mouse immunoglobulin G (IgG; 1:500; ProteinTech, USA) and/or fluorescein isothiocyanate–conjugated goat anti-rabbit IgG (1:500; ProteinTech, USA) for 1 hour, followed by staining with 4′,6-diamidino-2-phenylindole (DAPI; 1 μg/ml; Thermo Fisher Scientific, USA) for 15 min at room temperature. After thoroughly washing, fluorescent signals were visualized through LSCM (Leica, Germany).
Animal experimentation
Seven-week-old female BALB/c mice were purchased from the Vital River Laboratory Animal Technology Co. Ltd. (Beijing, China) and were cared and bred in the biosafety level 2 acility at the Institute of Microbiology, Chinese Academy of Sciences. Mice were intranasally administrated with virus-laden AFPs or solo viruses after being anesthetized by pentobarbital sodium. For in situ fluorescence imaging, mice were intranasally administered with virus-laden SiO2 particles or virus-laden polystyrene particles, and the frozen sections of organs were prepared and examined at 1, 3, and 8 dpi. Solo viruses and particles alone served as controls.
Peripheral blood was analyzed with a hematology analyzer (Nihon Kohden, Japan). Viral titers in the blood, nasal cavity, trachea, lung, liver, spleen kidney, brain, and heart were quantified via qRT-PCR analysis following the extraction of total tissue RNAs. H&E staining of various tissues was performed following the standard protocols (72, 73). For the immunofluorescence assay, frozen sections of the diverse tissues were prepared, and the viruses were visualized with an anti-HA antibody (1:200; Abcam, UK). All sections werescanned by a digital slide scanner (KFBIO, China).
Statistical analysis
Experimental results were presented as the means ± standard error of the mean. Statistical analysis was assessed with the independent t test or one-way analysis of variance (ANOVA) test using GraphPad Prism 8 software and Statistical Package for the Social Sciences software. In the current study, P values (P < 0.05 and P < 0.001) define the level of statistical significance.
Acknowledgments
We are grateful to J. Jin (Minzu University of China) for providing the PM2.5 samples and to Q. Ma (Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences) for providing the dust samples. We acknowledge L. Chen (Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences) for sharing the recombinant H1N1 viruses. We thank A. L. Chun of Science StoryLab for reading and editing the manuscript.
Funding: This work was supported by grants from the National Natural Science Foundation of China (gran numbers: 22150006, 21920102007, 22076210, 92043301, and 22021003), the Youth Innovation Promotion Association of Chinese Academy of Sciences (2022042), the National Key Research and Development Program of China (grant numbers: 2018YFA0901101 and 2021YFE0101500), and the International Collaboration Key Grant from the Chinese Academy of Sciences (grant number: 121311KYSB20190010).
Author contributions: Conceptualization: Z.D., J.M., M.F., and S.L. Methodology: Z.D., J.Q., Q.R., Q.S., X.D., G. Li., Y.Y.Z., Y.Q., and Y.L. Investigation: Z.D., J.M., X.D., Y.Q., G. Liu., M.F., and S.L. Supervision: J.M., M.F., I.L., and S.L. Writing: Z.D., J.M., X.D., Y.Y.Z., Y.Q., I.L., M.F., and S.L.
Competing interests: The authors declare that they have no competing interests.
Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials.