Ikechukwu B. Jacob1,2Akinkunmi O. Lawal1,2Salma S. Mahmoud1Emerson M. Kopsack1Erin S. Reynolds2Qinghe Meng1Hongkuan Fan3Paul T. Massa2,4Saravanan Thangamani2Hongpeng Jia5Guirong Wang1,2*- 1Department of Surgery, The State University of New York (SUNY) Upstate Medical University, Syracuse, NY, United States
- 2Department of Microbiology and Immunology, The State University of New York (SUNY) Upstate Medical University, Syracuse, NY, United States
- 3Department of Pathology and Laboratory Medicine, Medical University of South Carolina, Charleston, SC, United States
- 4Department of Neurology, The State University of New York (SUNY) Upstate Medical University, Syracuse, NY, United States
- 5Department of Surgery, Johns-Hopkins University, Baltimore, MD, United States
Introduction: COVID-19 remains a significant threat to public health globally. Infection in some susceptible individuals causes life-threatening acute lung injury (ALI/ARDS) and/or death. Human surfactant protein A (SP-A) is a C-type lectin expressed in the lung and other mucosal tissues, and it plays a critical role in host defense against various pathogens. The human SP-A genes (SFTPA1 and SFTPA2) are highly polymorphic and comprise several common genetic variants, i.e., SP-A1 (variants 6A2, 6A4) and SP-A2 (variants 1A0, 1A3). Here, we elucidated the differential antiviral and immunoregulatory roles of SP-A variants in response to SARS-CoV-2 infection in vivo.
Methods: Six genetically-modified mouse lines, expressing both hACE2 (SARS-CoV-2 receptor) and individual SP-A variants: (hACE2/6A2 (6A2), hACE2/6A4 (6A4), hACE2/1A0 (1A0), and hACE2/1A3 (1A3), one SP-A knockout (hACE2/SP-A KO (KO) and one hACE2/mouse SP-A (K18) mice, were challenged intranasally with 103 PFU SARS-CoV-2 or MEM medium (Sham).
Results: Infected KO and 1A0 mice had more weight loss and mortality compared to other mouse lines. Relative to other infected mouse lines, a more severe ALI was observed in KO, 1A0, and 6A2 mice. Reduced viral titers were generally observed in the lungs of infected SP-A mice relative to KO mice. Transcriptomic analysis revealed an upregulation in genes that play central roles in immune responses such as MyD88, Stat3, IL-18, and Jak2 in the lungs of KO and 1A0 mice. However, Mapk1 was significantly downregulated in 6A2 versus 1A0 mice. Analysis of biological pathways identified those involved in lung host defense and innate immunity, including pathogen-induced cytokine, NOD1/2, and Trem1 signaling pathways. Consistent with the transcriptomic data, levels of cytokines and chemokines such as G-CSF, IL-6, and IL-1β were comparatively higher in the lungs and sera of KO and 1A0 mice with the highest mortality rate. Furthermore, we observed the complexity of COVID-19, such as the difference between lung and systemic immune response to viral infection and of viral load and mortality among SP-A variants in this model.
Conclusion: These findings demonstrate that human SP-A variants differentially modulate SARS-CoV-2-induced lung injury and disease severity by differentially inhibiting viral infectivity and regulating immune-related gene expressions.
Introduction
COVID-19 is an infectious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (1, 2). Severe disease is characterized by a persistence of viral RNA and a robust influx of proinflammatory cytokines into the lungs, leading to acute lung injury (ALI) and acute respiratory distress syndrome (ARDS) (1, 2). Barely five years since first reported, there have been more than 7 million deaths due to COVID-19 globally (3). Although most people recover from acute infection without requiring hospitalization, some develop severe morbidity and mortality even in the absence of known risk factors (4). While Genome-Wide Association Studies (GWAS) have widely demonstrated differential disease outcomes in the general population, with multiple studies implicating genetic variations in host innate immunity, the reasons for the differential disease phenotypes are still unclear (4–8).
Human surfactant protein A (SP-A) is a C-type lectin (collectin) expressed by alveolar epithelial type 2 cells (AT2) of the lungs and other mucosal tissues. It plays critical roles in innate host defense, immunomodulation, and surfactant physiology in the lung and other mucosal surfaces (9). SP-A can bind glycoconjugates on microbes and interact with host cell receptors to inhibit infectivity and maintain lung homeostasis (9). Our data showed that SP-A can bind to SARS-CoV-2 spike protein (S), inhibit viral infectivity in lung epithelial cells, and protect against severe lung injury in humanized transgenic (hTG) mice (10). Other investigators have previously demonstrated the importance of collectins in attenuating viral infectivity and disease severity in the context of SARS-CoV-2, Influenza A virus (IAV), Human Immunodeficiency Virus (HIV), and Respiratory Syncytial Virus (RSV) infections (10–19).
The human SP-A protein is encoded by two highly functional and polymorphic genes (SFTPA1: SP-A1 and SFTPA2: SP-A2). Several genetic variants (alleles) of these genes have been identified in the general population, such as SP-A1 (variants: 6A2, 6A4) and SP-A2 (variants: 1A0, 1A3) (20). Although there are only a few amino acid differences between SP-A1 and SP-A2 variants, SP-A2 variants have a higher capacity to bind carbohydrates, enhance phagocytosis of invading pathogens by macrophages and induce higher proinflammatory cytokine production in immune cells under various stressed conditions, indicative of the potential of SP-A variants to differentially regulate immune function in the context of infection (21–24). Furthermore, SNPs resulting in an amino acid change from arginine to tryptophan in the 6A4 variant have been demonstrated to alter its structure and function and enhance susceptibility to chronic pulmonary diseases, including tuberculosis (20, 25–27). Previous population-based studies found that differences in collectin genes such as those of the Mannose-binding lectin (MBL), surfactant protein D (SP-D), and SP-A play important roles in disease severity following viral and bacterial infections (28–32); while dysfunctions in surfactant-related genes were shown to associate with severe COVID-19 (33). Importantly, recent GWAS studies have strongly implicated SP-A and SP-D genes in the differential COVID-19 susceptibility and severity observed in the human population (34, 35). Although the latter studies point to the potential host defense role of SP-A genes during infection, mechanistic pre-clinical studies of variant-specific functions of human SP-A genes in response to SARS-CoV-2 infection are lacking. Therefore, an adequate animal model to study the roles of SP-A genetic variants on host defense and lung function following SARS-CoV-2 challenge is vital.
In this study, we have generated and characterized double-humanized transgenic (double-hTG) mice expressing both hACE2 (SARS-CoV-2 cognate receptor) and single genetic variants of human SP-A (6A2, 6A4, 1A0 or 1A3). Using these mice, we explored whether there are differences in COVID-19 severity (mortality and ALI) and the underlying mechanisms of disease modulation by SP-A genetic variants post-SARS-CoV-2 (Delta variant) challenge. Our findings revealed distinct antiviral and immunomodulatory roles of SP-A genetic variants in COVID-19 symptomology among hTG mouse lines. The data indicate that SP-A genetic variants differentially modulate severe lung injury and survival rate in double-hTG mice by differentially regulating immune genes and inflammatory cytokine production locally and systemically. The data support ongoing efforts at producing personalized surfactant-based therapies for COVID-19 patients (36, 37).
Materials and methods
Mouse models
Mice breeding and maintenance were carried out at SUNY Upstate Medical University’s animal core facility. The animals were housed in a temperature-controlled room at 22°C under specific pathogen-free conditions. Mice between 8-16 weeks of age of both sexes were used in this study. The animal experiments were approved by SUNY Upstate Medical University Institutional Animal Care and Use Committee and conducted according to the National Institutes of Health and ARRIVE guidelines on using laboratory animals. The initial K18 (hACE2/mSP-A) mice were obtained from Jackson Laboratories (Bar Harbor, ME); these mice are susceptible to SARS-CoV-2 infection since they express the hACE2 transgene. To examine the role of human SP-A variants in the context of SARS-CoV-2 infection, we generated transgenic mice that carry both hACE2 and SP-A transgenes (double humanized transgenic or double-hTG) in which the hACE2 and individual single-gene variants of SP-A1 (6A2 and 6A4) and SP-A2 (1A0 and 1A3) are expressed and the mouse SP-A gene has been deleted by breeding with previously characterized hTG SP-A mice (38). To do this, we initially crossed K18 mice that are hemizygous for the hACE2 transgene with our SP-A knockout (KO) mice to remove the mSP-A gene. Genotyping was performed after each filial generation. Subsequently, we crossed the hACE2/SP-A KO mice with mice carrying individual variants of either SP-A1 (6A2 or 6A4) or SP-A2 (1A0 or 1A3) genes to generate a double-hTG mouse model with hACE2 and single-gene variants of human SP-A. We examined the presence of hACE2 and SP-A genes in the mice by PCR genotyping and analyzed SP-A expression in the lung by immunoblotting. The double-hTG mice (hACE2/6A2 (6A2), hACE2/6A4 (6A4), hACE2/1A0 (1A0) and hACE2/1A3 (1A3) alongside hACE2/SP-A KO (KO) and hACE2/mSP-A (K18) mice were subsequently used for SARS-CoV-2 (delta) challenge studies.
Mouse infection and sample collection
SARS-CoV-2- and mock-challenged mice (6 groups: 6A2, 6A4, 1A0, 1A3, K18, and KO) were anesthetized with isoflurane and infected intranasally (i.n.) with 30 µl (15 µl/nose) of virus solution containing 1× 103 PFU of SARS-CoV-2 in 1X MEM media. Control (Sham) mice were inoculated with 30 µl of 1X MEM. After viral infection, mice were observed daily for morbidity (body weight) and mortality (survival). Mice showing >25% decrease in their initial body weight were defined as reaching the experimental endpoint and euthanized. Mice were sacrificed by anesthesia and exsanguination on days 2, 4, and 6 post-infection (pi) to obtain lung samples for viral load analysis by RT-qPCR, immunohistochemistry (IHC), and plaque assay. Cytokine and gene-expression changes in the lung were analyzed, and systemic cytokines in sera were also determined.
Tissue processing
Mice were anesthetized with isoflurane and exsanguinated. Then, we perfused each mouse with 3 mL NBF. Then the lungs were inflation-fixed by tracheal instillation of 1 mL NBF. Mouse whole bodies were dissected along the abdomen (n = 7-11 mice/group) and fixed in NBF before histological and immunohistochemical analyses. Lungs from a subset of mice were harvested (n = 6-7 mice/group) and frozen immediately at -80°C for further processing. The lungs were weighed and homogenized in cold PBS and then centrifuged for 5 mins at 21130x g and supernatants were collected for viral load and serological analyses using RT-qPCR, plaque assay, and ELISA.
Histopathological analysis
Fixed lung tissues were embedded as previously described (39). 5 µm sections of individual tissue samples were cut and stained with Hematoxylin and Eosin (H&E) and then microscopically examined. Lung injury was scored using a 0-2 scale by counting the number of neutrophils in the alveolar space, neutrophils in the interstitial space, presence of hyaline membranes, proteinaceous debris filling the airspaces, and alveolar septal thickening as described by (40). Histopathological analyses were blinded and carried out by two independent pathologists.
Immunohistochemistry
Samples were analyzed by IHC as described previously by (39), with some modifications. Sections were deparaffinized in an oven for 1 h at 60°C and then rehydrated sequentially using 100%, 95%, 80%, 70% and 50% ethanol. The sections were rinsed twice in double distilled water and boiled in 10 mM citrate (pH 6.0) for 9 mins for antigen retrieval followed by cooling at room temperature for 20 mins. The sections were blocked with 10% BSA in PBS for 1 h at room temperature to prevent nonspecific background staining. This was followed by overnight incubation of the slides at 4°C with SARS-CoV-2 nucleocapsid protein (NP) rabbit antibody (cell Signaling, #33336S; 1:100). The next day, after several washes, the sections were incubated with biotinylated goat anti-rabbit antibody for 1 h (1:200), followed by staining using ABC/HRP complex (Vectastain ABC kit, peroxidase rat IgG PK-4004, and peroxidase rabbit IgG PK-4001; Vector Laboratories). Stained sections were visualized using 3’3-diaminobenzidine (DAB) (SK-4100; Vector Laboratories) and counterstained with hematoxylin (H-3404-100) (QS counterstain). Images were acquired using the Nikon Eclipse TE2000-U microscope (Nikon Corporation, Tokyo, Japan) and the number of SARS-CoV-2 NP-stained lesions was counted using Image J software.
Western blot analysis
Lungs from each group were homogenized in RIPA buffer (ThermoFisher Scientific, Rockford, IL) containing a cocktail of protease and phosphatase inhibitors (Roche, Indianapolis, IN). Total protein in lung homogenates was evaluated by the BCA assay kit (ThermoFisher Scientific), following the manufacturer’s instructions. Thirty micrograms of protein were subjected to SDS-PAGE on a 10% gel in reducing conditions and transferred to PVDF membranes (Bio-Rad). To block nonspecific interactions, The blots were incubated in a buffer with 5% non-fat milk for 30 mins before further incubation with SP-A polyclonal antibody (1:1000). Membranes were developed and observed by incubating with goat anti-rabbit IgG HRP-conjugated secondary antibody (Bio-Rad) and by using the ECL Western Blotting substrate (ThermoFisher Scientific).
Measurement of viral loads
Total RNA was isolated from homogenized lung tissues using the Quick-RNA extraction miniprep kit (# R1055 Zymo Research, CA, USA) following the manufacturer’s instructions and the concentration of RNA was measured using the nanodrop machine (Thermo Scientific). Real-time qPCR (RT-qPCR) was performed using the AB StepOnePlus Detection System (Applied Biosystems, Foster City, CA) using the one-step kit RT-PCR Master Mix Reagents (#64471423, Biorad), prepared according to the manufacturer’s instructions. Purified RNA (100 ng) was added into a 10 µl total volume of real-time PCR mix buffer containing forward/reverse primer pairs (forward, AGCCTCTTCTCGTTCCTCATCAC; reverse, CCGCCATTGCCAGCCATTC; each 500 nM) targeting SARS-CoV-2 N1 gene and a probe (250 nM, FAM) and other reagents provided by the manufacturer. The one-step RT-qPCR was performed through a cycle of reverse transcription at 55°C for 10 mins followed by 40 cycles of amplification at 95°C for 3 mins, 95°C for 15 s, and 55°C for 1 min. SARS-CoV-2 copy numbers in lung tissues were quantified using SARS-CoV-2 RNA standards as described by Winkler et al. (41).
Viral titer was also quantified using supernatants from lung homogenates by plaque assay. In brief, Confluent monolayers of Vero E6 in 24-well plates (1 × 105 cells/well seeded) were infected with 10-fold serial dilutions of supernatants from each mouse group. The cells were cultured for 1 h with intermittent rocking for viral adsorption. After which, the unbound virus was removed by washing with PBS, overlayed with 2% methylcellulose and cultured for another three days at 37°C in a humidified incubator. After plaque development, cells were fixed with 10% formalin for 60 mins and stained with 0.05% (w/v) crystal violet in 20% methanol and plaques were quantified as PFU/ml.
Transcriptomics analyses
The expressions of 84 innate and adaptive immune genes in the lungs of SARS-CoV-2 infected mice and Sham were determined using RT2 Profiler™ PCR array kit (Cat#: PAMM-052ZC-24, Qiagen) (42). In brief, total RNA was isolated from homogenized lung tissues using the Quick-RNA extraction miniprep kit (# R1055 Zymo Research, CA, USA) following the manufacturer’s instructions and RNA concentration and purity was determined by spectrophotometry using the nanodrop machine (Thermo Scientific). cDNA was synthesized from 500 ng of total RNA using the iScript Reverse Transcription Supermix for RT-qPCR (Cat#: 1708841, Biorad). Real-time PCR was performed using the RT2 Profiler™ PCR array system following the manufacturer’s recommendation in the AB StepOnePlus Detection System (Applied Biosystems, Foster City, CA). The levels of expression of the mRNA of each gene in the different groups were normalized to the housekeeping gene (HSP90). Data was exported into an Excel spreadsheet and analyzed with Qiagen’s PCR analysis web-based software (GeneGlobe Data Analysis Center). Relative gene expression changes were calculated by the 2−(averageΔΔCt) method. Fold change (FC) values were compared relative to the sham group or compared between two infected groups. Gene expression FCs are considered significant when FC > 2.0, p-value > 0.05). We further used specific primers to determine the expression of some significantly upregulated genes in infected hTG mice by RT-qPCR and then normalized to Gapdh expression (Supplementary Table 1 shows the primer list).
Signaling pathway association analysis
Differentially expressed genes (DEGs) were functionally analyzed using the Qiagen’s Ingenuity Pathway Analysis (IPA, QIAGEN, Redwood City, CA) software and the analyses included canonical pathway analysis, gene networks and upstream regulators. Canonical pathways were scored by analyzing the ratio of the number of genes that map to a specific pathway in the IPA database. The observed gene networks show functional relationships among DEGs based on database associations. IPA was used to gain additional insights into the immunoregulatory roles of each SP-A humanized mouse in response to SARS-CoV-2 infection.
ELISA of lung tissue cytokines
Levels of TNF-α (Lot# 227846-003), IL-1β (#372944-002), IL-6 (#330384-008), IL-10 (#355429-004), MCP-1 (#175884-002) and IFN-γ (#300277-008) in mouse lungs were measured by enzyme-linked immunosorbent assay (ELISA) (43). All measurements were in duplicate, according to the manufacturer’s instructions and the fold changes represented in the heat map among the infected mouse groups are relative to those observed in the mock-infected group.
Serum cytokine and chemokine analysis: multiplex cytokine assay
Inflammatory cytokines and chemokines in sera were analyzed using a Bioplex Mouse cytokine 23-plex panel magnetic bead Luminex assay (Biorad, 23-plex, Group 1, #64436152) following the manufacturer’s instructions. Twenty-five microliters of serum samples were diluted in 75 µl Bio-plex sample diluent and incubated with magnetic capture beads. Then, the plate was washed, and detection antibodies were added followed by streptavidin incubation. The analytes were measured using a Luminex 200 Analyzer (Luminex) and quantitated using a standard curve. Samples from each mouse were analyzed in duplicate (50 µl per well). All immunoassays were performed in the BSL-3 facility and fold changes represented in the heat map among the infected mouse groups are relative to those observed in the mock-infected group.
Statistical analysis
All experimental data are presented as mean ± standard error and statistically analyzed using GraphPad Prism 8.0 (GraphPad Software, San Diego, CA, USA). Comparisons between two independent groups were performed using Student’s t-test or multiple groups using one-way ANOVA. Survival analysis was performed with Kaplan-Meyer survival curves and evaluated statistically by the log-rank test. Results were considered statistically significant when P<0.05.
Results
Generation and characterization of double-hTG mice
To study the differential role of SP-A genetic variants in response to SARS-CoV-2 infection in the lungs, we generated double-hTG mice expressing both hACE2 (SARS-CoV-2 cognate receptor) and individual SP-A genetic variants derived from either SP-A1 (6A2 and 6A4) or SP-A2 (1A0 and 1A3) genes. This was done by crossing K18 mice expressing hACE2 and mSP-A with SP-A KO mice to generate mice that are positive for hACE2 and eliminate the mSP-A gene (Figure 1A). These KO mice, positive for hACE2 (hACE2/KO) were then bred with our previously characterized SP-A hTG mice carrying single-gene variants 6A2, 6A4, 1A0, or 1A3 to generate mice that carry both hACE2 and individual SP-A variant. The representative genotyping analysis of hACE2, SP-A1, SP-A2, and mSP-A genes from each mouse line is shown in Figure 1B. As expected, all the mice expressed hACE2, but KO mice had no human or mouse SP-A by PCR analysis. SP-A protein expression in the lungs of KO, K18, 6A2, 6A4, 1A0 and 1A3 mice was analyzed by Western blotting using an anti-SP-A antibody. Except for KO mice, SP-A was expressed in all SP-A hTG mice as expected (Figure 1C).
Figure 1. Characterization of double-humanized transgenic mice. Double-humanized transgenic (double-hTG) mice expressing both hACE2 and human SP-A1 (6A2 or 6A4), or SP-A2 (1A0 or 1A3) genetic variants were generated and characterized using the breeding strategy shown in Panel (A) We also generated mice with hACE2 but SP-A deficient (hACE2/KO) and those with hACE2 and mSP-A (K18). (A) Breeding scheme of double-hTG mice generated by crossing K18 mice with SP-A KO mice to generate hACE2/KO mice. hACE2/KO mice were subsequently crossed with previously characterized hTG SP-A mice. (B) Genotyping analysis of double-hTG, KO and K18 mice by PCR. (C) Western blotting analysis of SP-A expression in the lung of double-hTG mice showing hSP-A expressions in the lungs of double-hTG and K18 but not in mice deficient in SP-A (KO).
Significant body-weight loss and higher mortality rate in SP-A KO and 1A0 mice
Several population-based studies have shown that human SP-A variants may be either detrimental or protective depending on the disease (20, 44). Thus, to gain insights into the roles of human SP-A variants with regard to SARS-CoV-2-induced morbidity and mortality, we infected four different SP-A hTG mouse groups: 6A2, 6A4, 1A0 and 1A3, one K18 and one KO mice (Figure 2A). We observed that weight loss was particularly severe in KO mice and those carrying the 1A0 variant, particularly on days 5 and 6 post-viral challenge. On day 6 pi, KO, 1A0 and 6A4 mice had lost more than 13%, 12%, and 10% of their initial weights, respectively (Figure 2B). Although body-weight changes were observed early on day 3 pi. in KO and not in 1A0 mice, there was a drastic decrease in body weight in 1A0 mice by days 5 and 6 pi., suggestive of a differential response of the different SP-A variants to viral-challenge. Importantly, compared to KO mice on day 6 post challenge, mice carrying the 6A2 and 1A3 variants did not significantly lose body weight (P<0.05). The results until day 14 pi. showed that most of the mouse groups except KO, gradually regained their body weights until the end of the study. However, none of the KO mice survived beyond day 12 pi. Other clinically relevant characteristics observed in infected mice were ruffled fur, lethargy/unresponsiveness to stimuli, shivering, hunched posture, mucus around eyes (observed in the most severe mice), and labored breathing.
Figure 2. Body-weight loss and mortality in infected humanized mice. A total of six transgenic mouse lines, four of them carrying hACE2 alongside distinct human SP-A variants, i.e. hACE2/6A2 (6A2), hACE2/6A4 (6A4), hACE2/1A0 (1A0), and hACE2/1A3 (1A3), and the other two lines, i.e. hACE2/SP-A KO (KO) and hACE2/mSP-A (K18), were infected intranasally with 103 PFU/mouse of SARS-CoV-2 (Delta variant) or Sham (30 µl of MEM medium) on day 0; and animals were monitored daily for survival and body weight changes until 14 dpi or sacrificed at the indicated timepoints as shown in (A). (B, C) showed the body-weight changes and animal survival of each group, respectively. (n= 18-20 mice/group). Symbols in (B) represents mean ± SEM. *P< 0.05, **P< 0.01 (relative to sham), #<0.05, ##<0.01 (relative to infected KO), n=5 independent experiments. ns, not significant. Color corresponds to each mouse line.
On day 5 pi, some of the 1A0 and KO mice had reached mortality as defined by a loss of >25% of starting weight and were euthanized as a humane experimental endpoint (Figure 2C). By 6 dpi, about 45% of KO and 27% of 1A0 mice had died or reached the study endpoint for euthanasia. Strikingly, 100% survival was observed in mice carrying the 6A2 variant, closely followed by 88% and 82% survival observed in 6A4 and 1A3 mice, respectively, on day 6 pi. By day 14 pi. the results showed that mice carrying the 6A2 and 1A3 variants also had the highest survival rate (53.3% and 36.6% respectively) compared to 6A4 and 1A0 mice (26.6% and 12% respectively). KO mice, however, did not survive beyond day 12 pi. These findings suggest that mice carrying SP-A2 (1A0) genetic variant and those deficient in SP-A had more severe morbidity and higher mortality rates compared to those carrying SP-A1 (6A2 and 6A4) and SP-A2 (1A3).
Human SP-A variants differentially attenuate SARS-CoV-2-induced ALI in mice
As an important host defense protein and innate immune modulator, variations in the human SP-A gene have been associated with differences in susceptibilities and severities to acute and persistent pulmonary diseases (45–47). Therefore, we examined ALI features in mice carrying individual SP-A variants following SARS-CoV-2 infection. Lungs were collected on day 6 pi.; and the severity of lung injury in our mouse groups was evaluated in H&E-stained sections as recommended by the American Thoracic Society (40). In this study, we observed the presence of inflammatory cells in the alveolar space and interstitial (black arrows, Figure 3A, Supplementary Figure 1), proteinaceous debris within the airways, and a thickening of alveolar walls resulting from a loss of alveolar integrity and fluid accumulation (red arrows, Figure 3A) in infected lungs. As expected, none of these characteristics was found in the lung of sham mice (Figure 3A).
Figure 3. Differential lung injury after SARS-CoV-2 infection of humanized mouse lines. (A) Representative histological images of lungs from each mouse group on day 6 p.i. SARS-CoV-2 induced a more severe lung damage in KO mice characterized by an infiltration of inflammatory cells in alveolar spaces and interstitial (black arrows), accumulation of proteinaceous debris, and edema leakage of protein-rich fluid in the alveolar space (red arrows) compared to K18 and hTG SP-A1 6A2 (representative images are shown in Panel A). (B) Semi-quantitative histological assessment of lung injury from each mouse group showed significantly higher scores of lung injury in all infected groups compared to sham group, and differential lung injury scores were observed in mice carrying unique SP-A genetic variant. Symbols in (B) represents mean ± SEM of individual animals. N = 7-12 mice/group. *P< 0.05, **P< 0.01 (relative to sham), #<0.05, ##<0.01 (between two infected mouse groups from 5 independent experiments). Scale bars= 50 µm.
Interestingly, the markers of severe ARDS were more overt and severe upon scoring in KO mice relative to other groups, suggestive of the role of SP-A in host defense and maintaining alveolar integrity post-SARS-CoV-2 infection (KO vs hTG, P<0.01). However, a comparison of lung injury scores among the variants showed that 6A4 and 1A3 mice had the most attenuated ALI compared to those carrying 6A2 and 1A0 variants (the order of ALI score: 1A0>6A2>1A3>6A4) (Figure 3B). Further analysis of SP-A1 genetic variants showed that lung injury score trended higher in 6A2 relative to 6A4 mice while a comparison of SP-A2 gene-specific variants revealed a significantly higher lung injury score in 1A0 compared to 1A3 mice. These observations suggest that human SP-A variants 6A4 and 1A3 may be more protective against SARS-CoV-2-induced ALI compared to 6A2 and 1A0. Together, the results highlight two distinct outcomes of SARS-CoV-2-induced pathology in our experimental mouse groups based on body-weight changes and survival, and lung injury scores.
Human SP-A variants differentially attenuate SARS-CoV-2 load in mice
Given that viral burden has been associated with COVID-19 severity, we next quantified and compared viral load 6 dpi in infected mouse lungs from each group by RT-qPCR, IHC, and plaque assay (48, 49). IHC staining of viral N proteins revealed higher levels of viral antigens in the alveolar septum and mononuclear cells in the lungs of SP-A deficient mice relative to the other hTG mice on day 6 pi (P<0.01). However, we observed comparatively reduced N proteins in the lungs of mice carrying the 6A2 and 1A0 variants on Day 6 p.i. (Figures 4A, B). Relative to KO mice with the highest viral genome copy numbers, the other groups (except K18) had significantly reduced copy numbers (p<0.05). Among the variants, the 1A3 variant had the highest average viral RNA copies (Log10 8.6) (the order of viral RNA copies: 1A3>6A4>6A2>1A0) (Figure 4C). Furthermore, on day 6 pi, we were only able to determine infectious viral particles more consistently in KO, K18, and 1A3 mice (~103 PFU/mL) while other SP-A hTG lines had mostly cleared the virus and infectious virus was therefore not observed at the limit of detection (LOD) (Figure 4D). To address the innate antiviral role of human SP-A, infectious virus titer was quantified in lungs collected on day 2 pi since there was viral clearance in most SP-A hTG lungs by day 6 pi. We observed a significant decrease in viral titer in K18, 6A2, and 1A0 mice relative to KO mice, with 1A0 mice having the lowest viral burden in the lungs; suggestive of the more profound antiviral role of 1A0 and 6A2 variants in response to SARS-CoV-2 infection in vivo (Figure 4D). Further analysis of lungs collected on day 4 pi showed a decline in viral titer compared to those on day 2 pi. The results suggest that SP-A may attenuate SARS-CoV-2 severity through its inhibitory effect on virus infectivity, while also highlighting the potential role of other distinct functions of SP-A in modulating disease severity, besides its antiviral function.
Figure 4. Differential viral load in the lung of SARS-CoV-2 infected humanized mouse lines. Viral load in the lungs of SARS-CoV-2-infected mice were determined using IHC, RT-qPCR and plaque assay. (A) IHC analysis was performed for SARS-CoV-2 nucleocapsid protein (NP) in the lungs of infected and sham groups on day 6 pi. Viral antigen was observed in all infected groups but not in sham. The arrows indicate the presence of viral antigens in the alveolar septum and cells. (B) Quantification of NP-positive lesions in each lung section reflected a significantly higher number of viral antigens in infected KO relative to hTG SP-A lines. Significant differences in viral antigens were observed among hTG SP-A mouse lines. (n = 7 mice/group) (C) Lungs were collected from infected and sham mice on day 6 pi and viral load was determined by RT-qPCR quantification of the N gene. The horizontal dotted line shows the limit of detection (LOD) observed in sham controls. The data indicate higher viral copy number in KO mice compared to SP-A hTG mice. Significant differences in genome copies were observed among hTG SP-A mouse lines. (n = 3-7 mice/group) (D) Plaque assay quantification of infectious virus in lungs collected on days 2, 4, and 6 pi. We observed peak viral load on day 2 pi. Significant differences in infectious virus titer were observed among hTG SP-A mouse lines on day 2, but not on day 4 and 6. The dots below the LOD in Figure 4D (6 dpi) represent the number of mice without any detectable virus *P< 0.05 **P< 0.01 (relative to sham), #<0.05, ##<0.01 (between two infected groups) from 3-4 independent experiments. Scale bars= 50µm. Symbols in (B–D) represent mean ± SEM of individual animals.
Human SP-A variants differentially modulate innate and adaptive immune gene expressions in the lungs of SARS-CoV-2 infected humanized mice
Given that human SP-A can interact with both viral and host proteins to modulate immune responses, we assessed the role of SP-A in modulating immune gene expression post-viral challenge in this mouse model (50). We used a microarray-based system that allowed us to study the transcriptional changes of a focused panel of 84 innate and adaptive immune genes in the lungs from infected and sham groups. In the hierarchical clustering of genes illustrating gene-expression changes in infected vs mock-challenged mice, we observed that KO and 1A0 mice (these animals showed the most severe disease) had the greatest numbers of differentially expressed genes (DEGs) such as Il18, NLRP3, Jak2, Stat3, Irak1 and Il1r1 (Figure 5A); reflecting a more enhanced activation of inflammatory genes in these groups. Of the 84 genes studied, 17 and 15 genes were significantly upregulated in infected KO and 1A0 relative to mock-challenged mice respectively (DEGs were ranked by >2.0-fold change, p<0.05) (Figure 5A). Supplementary Table 2 shows all the DEGs in infected and sham groups.
Figure 5. Differential transcriptional signatures in innate and adaptive immune genes in the lung of SARS-CoV-2 infected humanized mouse lines. A total of 84 innate and adaptive immune genes were examined in the lung tissues of six infected mouse lines, including sham mice. (A) Hierarchical clustering of genes illustrating changes in gene expression compared between infected mice and sham mice (n=3 mice/group). Of the 84 genes studied, 17 and 15 genes were significantly upregulated in infected KO and 1A0 mice relative to mock-challenged mice, respectively. (B) Volcano plots of differentially expressed genes (DEGs>2.0-fold change, p<0.05) between infected K18 mice relative to mock-challenged mice. (C) Volcano plots of some significant DEGs between two infected groups (Top panels: KO vs 6A2; KO vs 6A4; KO vs 1A0; KO vs 1A3; 6A2 vs 6A4. Bottom Panels: 6A2 vs 1A0; 6A2 vs 1A3; 6A4 vs 1A0; 6A4 vs 1A3 & 1A0 vs 1A3). The data revealed an upregulation in immune genes such as MyD88, Stat3, & Icam1 in KO vs 6A2; IL17a, Csf2 and Casp1 in KO vs 1A0 and Mapk1 in 6A2 vs 1A0 mice. (D) RT-qPCR validation of DEGs in SARS-CoV-2 infected hTG Mice. Results showed significant differences in several gene expressions between infected mice and sham, as well as among hTG SP-A mouse lines. *P< 0.05, **P< 0.01 (relative to sham), #<0.05, ##<0.01 (between two infected mouse groups).
The analysis further revealed a significant upregulation of genes such as Nod2, Ifnb1, and Irf7 involved in interferon induction and antiviral response in infected K18 mice relative to sham (Figure 5B, Supplementary Table 2). Further analysis revealed that the antiviral molecule, Mx1 was significantly upregulated in K18 mice and appears to be a key innate host defense mediator in this line relative to all other groups as previously observed in critical COVID-19 patients (51). Strikingly, we further observed the upregulation of Myeloid Differentiation Primary Response 88 (Myd88), a unique signaling adaptor protein shown to activate inflammatory cytokines through the NF-kB and MAPK1 pathways following SARS-CoV-2 infection in KO mice compared to other lines (Figure 5C - Top panel; Supplementary Table 2). Other significant genes upregulated in KO mice relative to SP-A hTG and K18 mice are Mx1, Icam1, Stat3, Nod2, Casp1, and Mapk1. The more DEGs involved in cytokine induction and macrophage activation such as Mapk1, Irak1, and Il18 observed in KO and 1A0 mice on day 6 pi reflects a more sustained activation of inflammatory genes and, thus, worse pathology in these lines.
The Volcano plots of the two SP-A1 variants revealed the upregulation of CD86, a costimulatory molecule for T-cell mediated immunity activation, and a downregulation of Lyz2 (a marker for macrophage activation) in infected 6A2 mice (moderate mortality) vs 6A4 mice (severe mortality), suggestive of the important role of T cell-mediated immune responses in virus clearance. In contrast with the low number of DEGs observed in SP-A1, SP-A2 genotype 1A0 (critical disease) had several significantly upregulated proinflammatory immune genes such as Il18, Mapk1, Irak1, Nfkb1, and Ticam1 compared to 1A3 (moderate disease). Interestingly, Mapk1 was the only significant DEG downregulated in infected 6A2 mice relative to 1A0 mice. In 1A3 mice with moderate disease, 2 of 4 upregulated genes relative to sham were involved in antiviral state (Mx1 and STAT1), neutrophil/monocyte recruitment (CCL2), and anti-inflammatory response (IL10). Thus, mice with critically severe disease (KO and 1A0) had more significant upregulation of immune gene expression and the higher activation of these genes may have contributed to the severe morbidity and mortality observed in these groups. Thus, differences in the capacity of human SP-A genetic variants in inflammatory gene modulation in the context of SARS-CoV-2 infection could result in differential disease severity post-SARS-CoV-2 infection (Figures 5B, C, Supplementary Table 2). Furthermore, we validated by RT-qPCR a select group of essential DEGs in disease development and observed remarkable differences in their expression levels (Figure 5D).
Human SP-A variants differentially modulate inflammatory signaling pathways in SARS-CoV-2 infected humanized mice
We analyzed the DEGs using the Ingenuity Pathway Analysis (IPA) tool to identify the signaling pathways implicated in each hTG mouse line. The DEGs were defined as those genes with fold change ≥ 2 and p-value < 0.01. Data analysis of the respective genes that mapped to a given canonical pathway obtained using a ratio of the number of significant DEGs to those that map to a given pathway revealed the upregulation of major inflammatory response pathways such as, “Pathogen-induced cytokine storm”, “neuroinflammation”, “role of pattern recognition receptors in recognition of bacteria and viruses”, “Trem1” and “macrophage classical activation” signaling in infected hTG lungs relative to sham lungs (Figure 6A). These findings support the robust induction of cytokines and disease development in infected mice compared to sham mice. A comparison of the pathways specifically upregulated in mice with the most severe mortality (KO) compared to those with moderate mortality (6A2) showed that most genes mapped to Trem1 (p-value = 1.26E-7) and IL33 (p-value = 1.76E-6) signaling pathways. IL-33 activation has been reported in advanced COVID-19, contributing to severe lung inflammation by activating several proinflammatory cytokines such as IL-6 and TNF-α through the MyD88-MAPK-NF-kB signaling axis (52–54). As shown in Figures 6A, B, other pathways identified in KO mice relative to SP-A hTG lines with moderate disease are those related to the activation of inflammatory responses as described above. Several unique pathways activated in 1A0 vs 1A3 are neuroinflammation (p-value = 9.87E-21), “Toll-like receptor” (p-value = 1.22E-19), “pathogen-induced cytokine storm” (p-value = 1.32E-17) signaling pathways, indicating that the severe mortality and ALI observed in KO and 1A0 mice (groups with the most severe disease) may be due to an over-activation of immune genes resulting in a dysregulated cytokine storm in these mouse groups compared to those carrying the 6A2 and 1A3 (moderate disease) variants.
Figure 6. Top canonical pathways identified from functional analysis of DEGs in the lung of infected mouse lines. Top canonical pathways were identified using ingenuity pathway analysis (IPA). Canonical pathways were identified using only the significant up- and down-regulated genes from our microarray studies and ranked according to their –log (p-value) (Orange bars). Color intensity highlights the magnitude of fold change expressions. (A) Top one to five canonical pathways identified in SARS-CoV-2 infected double-hTG mice relative to sham mice. (B) Top five canonical pathways identified between two infected mouse groups: KO vs 6A2, KO vs 1A0, 6A2 vs 1A0, and 1A0 vs 1A3.
Differential cytokine signatures in the lungs of SARS-CoV-2 infected humanized mice
Since the most severe COVID-19 groups (KO and 1A0) showed significant upregulation of multiple immune-related genes in the lungs, we evaluated the inflammatory signatures in mice carrying individual SP-A single-gene variants by quantifying the level of some cytokines in the lungs of infected and mock-challenged mice. We observed that IL-1β, IL-6, and IL-10 were significantly higher in KO, 1A0, and 6A4 compared to 6A2 and 1A3 mice (Figures 7A, B); but there was no significant difference in the levels of these cytokines in infected 6A2 and 1A3 mice relative to sham. The latter supports the role of elevated cytokines as major contributors to severe disease observed in the most severe groups. TNF-α was increased in all infected lines compared to sham, however, its level trended lower in 6A2 relative to other double-hTG lines. A Comparison of cytokine profiles between the two SP-A1 genotypes showed that IL-10 and IFN-γ were significantly higher in 6A4 mice (severe disease) compared to 6A2 mice (moderate disease). Data analysis of SP-A2 gene variants revealed higher levels of IL-1β, IL-6, and IL-10 in 1A0 mice (severe disease) vs 1A3 mice (moderate disease). These observations also correlated with the significant upregulation of immune genes in the lungs of KO mice and 1A0-carrying mice and consistent with cytokine profiling in COVID-19 patients which showed higher IL-6, IL-10, and TNF-α as correlates of severe COVID-19 (1, 55, 56).
Figure 7. Differential cytokine expressions in the lungs of SARS-CoV-2 infected humanized mouse lines. (A) Bar graphs of absolute concentrations of cytokines and chemokines using ELISA in the lungs of infected and sham mice on day 6 pi. (B) Heatmap representation of cytokine and chemokines levels in the lungs of SARS-CoV-2 infected mice normalized to levels in sham mice. Symbols in (A) represents mean ± SEM of individual animals. *P< 0.05, **P< 0.01 (relative to sham), #<0.05, ##<0.01 (between infected mouse groups), n= 6-10 mice/group.
Differential cytokine signatures in sera of SARS-CoV-2 infected SP-A humanized mice
In severe COVID-19 patients, pulmonary and extrapulmonary manifestations such as acute kidney, intestinal, and brain injuries are observed (55). These extrapulmonary manifestations are often associated with increased levels of systemic cytokines and fatal SARS-CoV-2 infection (55, 57). Therefore, we further measured the levels of 23 cytokines in the sera of SARS-CoV-2-infected and sham mice using a multiplex assay. Data analysis revealed a similar concentration of cytokines such as IL-1α, IL-2, GM-CSF, and IL-13 between infected and sham mice irrespective of the individual SP-A genotype (Figures 8A, B). Interestingly, as observed in human COVID-19 patients with fatal disease (KO and 1A0) and those with severe disease (6A4) had overall, significantly increased systemic levels of G-CSF, IL-6, IL-10, TNF-α, IL-12 (p40), IL-17A, CCL5, and CCL11 compared to those with moderate disease (6A2 and 1A3) (55, 58, 59). However, we found that serum levels of the T-cell associated cytokine (IFN-γ) were markedly elevated in mice carrying 6A2 and 1A3 variants (moderate disease) relative to KO and 1A0 (severe disease) mice. This finding is in line with previous reports on the importance of IFN-γ signaling and T-cell mediated immune responses in alleviating severe COVID-19 (60, 61).
Figure 8. Differential cytokine signatures in sera of SARS-CoV-2 infected humanized mouse lines. (A) Bar graphs of absolute levels of some cytokines obtained from multiplex assay in the sera of infected mice relative to sham on day 6 pi. (B) Heatmap representation of cytokine and chemokines levels quantified using a multiplex platform in sera of SARS-CoV-2 infected mice on day 6 pi normalized to levels in sham mice. Symbols in (A) represents mean ± SEM of individual animals. *P< 0.05, **P< 0.01 (relative to sham), #<0.05, ##<0.01 (between infected mouse groups), n= 3-5 mice/group.
Discussion
Besides identifying those susceptible to a given disease, the goal of precision medicine is to tailor treatment and prevention strategies using the unique gene set of the individual (62). Although most SARS-CoV-2 cases are mild, often requiring no hospitalization, some patients develop fatal respiratory distress in the absence of known risk factors like smoking and pre-existing medical conditions (63–65). In-depth GWAS analyses have linked these COVID-19 inter-individual variabilities to genetic variations in host proteins such as those of the blood group system, SARS-CoV-2 entry receptors, and innate and adaptive immune responses as drivers of severe disease observed among apparently healthy individuals (34, 35, 66–70).
As important as GWAS studies are at finding associations, pre-clinical studies are needed to determine the causal roles of these genes in human disease. Herein, we demonstrated that mice carrying single-gene variants of human SP-A1 and SP-A2 differentially modulate COVID-19 severity in response to SARS-CoV-2 infection. The results showed that (1) individual variants of human SP-A genes alleviate SARS-CoV-2 morbidity and mortality via their differential antiviral effects, (2) immune genes are differentially activated in mice carrying unique human SP-A variants, and (3) cytokine levels vary systemically and in the lungs in the mouse lines following SARS-CoV-2 infection and their levels directly correlate with severe disease. These observations mechanistically elucidate the roles of human SP-A genetic variants in the observed differences in COVID-19 symptomatology in the general population for the first time.
To decipher the innate immune roles of human SP-A genetic variants, we determined lung injury and animal mortality following SARS-CoV-2 infection and evaluated the established correlates of severe disease - viral load and inflammatory cytokine induction in mice (49, 57). In this report, we observed two distinct patterns of COVID-19 severity in our hTG mouse lines, defined by survival rate and ALI severity, respectively. First, we observed that SP-A is important in preventing severe COVID-19 since SP-A deficient mice had the highest mortality rate and ALI scores compared to the other infected lines with either mouse SP-A (K18) or the respective human SP-A genes. This supports previous studies that observed higher viral titers in animals and patients with severe morbidity and mortality and reinforces the significant role of host innate pattern recognition molecules such as collectins in alleviating severe COVID-19 by inhibiting viral infectivity as was previously shown in the context of bacterial-induced pneumonia, RSV, flu, and recently, SARS-CoV-2 (10, 13, 48, 57, 71–76). Importantly, since the 6A4 and 1A3 mice had a more attenuated ALI, we speculate that there could be a small subset of individuals particularly protected from severe lung damage after SARS-CoV-2 infection (77, 78). In addition, individuals predominantly carrying the 6A2 and 1A3 alleles may have better survival rates following infection than those carrying 1A0 and 6A4 alleles as previously observed with better lung function and mechanics in 1A3 vs 1A0 mice post-bacterial infection (47). The profound ALI and death in the SP-A2 (1A0) line was, however, intriguing and deviates from previous reports of high viral titers as drivers of severe disease given the significantly reduced viral titers observed in this line (13, 76). The finding however aligns with a previous report of an association between the 1A0 variant and an increased risk of acute respiratory failure after IAV infection in humans and an animal model of pneumonia (28, 47). Therefore, besides its antiviral function, we speculated that SP-A may modulate disease severity in response to SARS-CoV-2 infection via its other effector functions.
It has been established that COVID-19 is a cytokine-release syndrome caused by a dysregulated and over-activated immune response to infection. We speculate that while human surfactant protein (SP) genes play vital host defense functions against SARS-CoV-2 infectivity, their capacity to modulate inflammation in response to an infection may be even more crucial for disease modulation as was previously observed (23, 29, 43, 76, 79). Collectins such as SP-A, SP-D, and MBL largely serve as first line of defense against invading microbes by neutralizing, opsonizing, and facilitating microbial phagocytosis to enhance antigen presentation and stimulate the adaptive immune response (80). Importantly, from both immunological and biophysical standpoints, SP-A genetic variants differ in several instances: (1) the capacity for macrophages to phagocytose bacteria was higher in the presence of the 1A0 variant compared to 6A2; (2) the 1A0 variant was shown to activate inflammatory cytokine production more than 6A2 variant; (3) Structurally, SP-A2 gene-specific variants (1A0 and 1A3) are more stable and may better withstand oxidation by reactive oxygen radicals during inflammation compared to SP-A1 gene-specific variants (6A2 and 6A4); and (4) The proteomics and miRNAs of alveolar macrophages are different between SP-A1 and SP-A2 genes (22, 23, 47, 79). Thus, the severe ALI and mortality observed in 1A0 mice may be due to some of the latter functional capacities of the SP-A2 (1A0) gene.
Furthermore, in this study we observed the differential role of human SP-A variants in lung and systemic immune response to SARS-CoV-2 infection, and different effect of human SP-A variants on the direct antiviral function and the immune modulation after viral infection in the humanized murine COVID-19 model, which can explain the observation, i.e. the lower viral load in the lung and higher levels of proinflammatory cytokines with severe morbidity and mortality in the SP-A2 (1A0) mice. These demonstrated that to the complex genetic variation of human SP-A genes also contributes the complexity of COVID-19 pathology and individual susceptibility to SARS-COV-2 infection.
The transcriptomics data supported our hypothesis and revealed enhanced activation of immune genes in SP-A deficient and 1A0 mice compared to the other mouse lines. It appears that in the absence of SP-A, viral replication is unimpeded, resulting in an overstimulation and recruitment of immune cells, resulting in severe disease. However, the high mortality observed in the 1A0 line may be due to more viral aggregation, internalization by alveolar macrophages, enhanced innate immune signaling, and greater presentation of antigenic peptides in the lung environment to induce robust immune cell infiltration into the lungs. Insufficient, or delayed type 1 responses resulting in prolonged hyperinflammation, cellular influx, cell death, and lung fibrosis are some hallmarks of severe COVID-19. A delicate balance must thus be struck by innate and adaptive immune responses against the virus and modulation of these inflammatory responses to prevent sustained and hyperinflammatory tissue damage in the context of COVID-19. It is essential to state that our transcriptomic data were performed using samples collected 6 dpi, which is indicative of a sustained innate pro-inflammatory immune response in KO and 1A0 mice and consistent with persistent inflammation observed in severe patients (81, 82).
The importance of T-cell-mediated immunity in preventing COVID-19 death was particularly highlighted with the upregulation of Cd86, a T-cell co-stimulatory molecule in 6A2 mice with moderate mortality compared to 6A4 mice with severe mortality. SARS-CoV-2, upon TLR-2 ligation, activates the Myd88 adaptor protein, subsequently signaling either through Mapk1 or NF-kB pathways (83). Interestingly, SP-A can activate NF-kB and induce cytokine production by directly interacting with TLR-2 and TLR-4 or by differentially regulating miRNA, which targets both proinflammatory and anti-inflammatory cytokine production (84–87). Moreover, SP-A was shown to reduce Mapk1 phosphorylation and activation (88). The downregulation of Mapk1 in 6A2 mice relative to 1A0 is an interesting observation as it may impact cytokine production, which may have contributed to the higher survival rate recorded in mice carrying the 6A2 gene. Therefore, sequence changes in 6A2 could have caused functional differences in any of the latter mechanisms, impacting Mapk1 activation and, potentially, cytokine levels in 6A2 mice.
A recurring canonical signaling pathway identified from our dataset is Trem1, an activating receptor mainly expressed in macrophages and neutrophils (89). Activation of Trem1 by soluble molecules, microbial ligands, inflammasome components, and cytokines through TLR ligation causes downstream activation of inflammatory cytokines and amplifies innate immune responses (89, 90). Studies have shown that NLRP3 inhibition decreased Trem1 activation, resulting in diminished proinflammatory cytokine levels and disease severity (90, 91). We further observed enrichment of Trem1 in infected KO compared to 6A2 mice, suggestive of the potential role of the 6A2 variant in attenuating severe disease by downregulating downstream inflammatory cytokine pathways. Further studies should address the underlying mechanism involved in the differential modulation of inflammatory responses via Trem1 signaling by SP-A variants. Do SP-A variants bind to Trem1 to interfere with subsequent activation by proinflammatory cytokines and viral structural proteins in myeloid cells (89, 91). It remains to be elucidated whether SP-A variants attenuate inflammatory cytokine production by directly binding to Trem1 or by indirectly interacting with TLRs to reduce cytokine production and mitigate the positive feedback loop of cytokines on Trem1 activation.
To explain the differences in disease severity, we speculated that 1A0 mice may be better at inducing inflammatory responses than other SP-A genetic variants, given the high number of DEGs observed from our transcriptomic data. Our immunoassay of inflammatory cytokines in the lungs strongly supports this speculation where we observed higher levels of these cytokines in KO and 1A0 mice with the highest mortality compared to 6A4, 6A2, and 1A3 with high to moderate mortality rates (Mortality Rate: KO>1A0>6A4>1A3>K18>6A2). Certain cytokines directly correlate with severe disease or death, such as IFN-γ, TNF-α, IL-1β, IL-6, and IL-10 (55, 56). The higher levels of IL-6 and IL-10 in KO and 1A0 mice in our study strongly correlate with previous clinical studies (56, 57, 82). We also speculate that variants 6A2 and 1A3 may be better than 1A0 at attenuating the levels of proinflammatory cytokines in macrophages stimulated with microbial or immune molecules as previously observed in functional studies using native SP-A and SP-A single-gene variants in vitro (23, 79, 92).
Cytokine expression and release in the context of SARS-CoV-2 infection is not only associated with ALI/ARDS but also with the development of multiple organ dysfunction (MOD) in COVID-19 patients (55). MOD has been linked with broad ACE2 expression, active viral replication in extrapulmonary organs, and high levels of systemic inflammatory cytokines (93). We therefore compared the levels of inflammatory cytokines in sera of individual hTG mice and observed that 6A2 and 1A3 mice with significantly higher survival rates had remarkably reduced levels of most cytokines except IFN-γ, a major inducer of cell-mediated immune responses, which are vital for clearing virally infected cells. Indeed, the higher systemic levels of IFN-γ observed in 6A2 mice support the double-edge sword theory of IFN-γ in COVID-19 pathology in which robust levels in blood may limit virus spread to other tissues, but increased levels in lungs may drive lung pathology as observed in the most severe mouse groups (61, 94). This finding warrants further studies on the role of human SP-A genetic variants in modulating extrapulmonary disease manifestations of COVID-19.
SP-A exhibits functional and quantitative differences that are influenced by the unique SP-A haplotype, sex, and age of an individual, which ultimately influences health outcomes following infection due to transcriptomic and proteomic differences (46, 47, 84). One study showed enhanced miRNA activation in male mice, significantly impacting signaling molecules and proinflammatory cytokine induction (84). We therefore infected mice of similar ages in this study and observed that ALI severity and mortality trended higher in male K18, 6A4, and 1A0 compared to female mice. However, due to sample size limitation, a definitive contribution of sex in SP-A function cannot be ascertained and future studies should delineate the intricate interplay of age, sex, and the unique SP-A variant in the context of SARS-CoV-2 infection. Viral strain differences in pathogenicity have been observed in animals and humans, and we speculate that the observed disease in the mouse groups may vary if SARS-CoV-2 variants other than Delta are used (95).
Changes in the levels of pulmonary surfactant proteins (SP-A and SP-D) and MBL have been shown to correlate with severe disease (96). Our previous data revealed a significant reduction in SP-A levels among severe COVID-19 patients while a dysregulation in SP genes was observed in SARS-CoV-2 infected lung tissues (10, 33). The decrease in collectin levels post-primary infection should be further addressed in a co-infection model and emphasizes the need to develop lung surfactant treatments for people particularly predisposed to severe COVID-19.
For the first time, we demonstrated that human SP-A genetic variants differentially attenuate ALI and affect survival rates post SARS-CoV-2 infection, by reducing viral burden, and differentially modulating viral-induced inflammatory responses. Here, we have demonstrated that compared to SP-A2 (1A0), mice carrying the SP-A1 gene- specific variant (6A2) have the highest survival rate, followed by SP-A2 (1A3) and SP-A1 (6A4) and survival rate inversely correlated with lung cytokine levels. SP-A genetic variants (particularly 6A2 and 1A3) seem more capable of attenuating the robust induction of cytokines observed in severe COVID-19 and SP-A gene variants are protective (albeit differentially) in response to SARS-CoV-2 infection by distinctively activating immune genes and signaling pathways (Figure 9). We, therefore, speculate that morbidity and mortality observed in some infected individuals may be explained by the unique and predominant SP-A genotypes in such individuals. Understanding the contributions of human SP-A genetic variants in response to SARS-CoV-2 infection is vital to developing novel gene/variant-specific surfactant-based therapies for COVID-19 patients.
Figure 9. A model of the differential role of human SP-A genetic variants against SARS-CoV-2-induced ALI and mortality. SP-A variants differentially inhibit viral load in the lung after SARS-CoV-2 inoculation. ALI and disease severity are characterized by enhanced expressions of immune genes as well as higher levels of cytokines and chemokines in the lungs. Infected SP-A KO and 1A0 mice exhibited higher levels of innate and adaptive immune gene expression and elevated levels of pro-inflammatory cytokines and chemokines compared to infected K18, 6A2, 6A4 and 1A3 mice.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The animal study was approved by SUNY Upstate Medical University Institutional Animal Care and Use Committee. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
IJ: Conceptualization, Data curation, Formal analysis, Writing – review & editing, Investigation, Methodology, Software, Visualization, Writing – original draft. AL: Investigation, Methodology, Visualization, Writing – review & editing. SM: Methodology, Writing – review & editing. EK: Methodology, Writing – review & editing. ER: Methodology, Writing – review & editing. QM: Methodology, Writing – review & editing, Resources. HF: Methodology, Writing – review & editing. PM: Writing – review & editing, Methodology, Supervision. ST: Methodology, Writing – review & editing, Resources, Supervision. HJ: Methodology, Writing – review & editing, Funding acquisition. GW: Funding acquisition, Methodology, Project administration, Writing – review & editing, Conceptualization, Formal analysis, Resources, Supervision, Validation, Visualization.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by NIH R01HL136706, R21AI171574, and the National Science Foundation research award (1722630) (to GW) and by the 1R01AI148446-01A1 (to HJ).
Acknowledgments
We thank the staff of the Upstate Vector Biocontainment Laboratory (BSL-3) for their technical support and SUNY Upstate Pathology Research Core (SUNY SPORE) for their histological services. We would like to thank Dr. Joanna Floros for her generous support with the biomaterials used for this project and her scientific perspective on this study. Vero E6 cells were obtained through BEI Resources, NIAID, and NIH.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2025.1462278/full#supplementary-material
References
1. Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. (2020) 395:497–506. doi: 10.1016/S0140-6736(20)30183-5
2. Shuai H, Chan JF, Hu B, Chai Y, Yuen TT, Yin F, et al. Attenuated replication and pathogenicity of SARS-CoV-2 B.1.1.529 Omicron. Nature. (2022) 603:693–9. doi: 10.1038/s41586-022-04442-5
3. WHO COVID-19 dashboard (2024). Available online at: https://covid19.who.int/ (Accessed January 23, 2024).
4. Anastassopoulou C, Gkizarioti Z, Patrinos GP, Tsakris A. Human genetic factors associated with susceptibility to SARS-CoV-2 infection and COVID-19 disease severity. Hum Genomics. (2020) 14:40. doi: 10.1186/s40246-020-00290-4
5. Verma S, Verma S, Khan FH, Siddiqi Z, Raza ST, Abbas M, et al. Genetic polymorphisms of IL6 gene -174G > C and -597G > A are associated with the risk of COVID-19 severity. Int J Immunogenet. (2023) 50:5–11. doi: 10.1111/iji.12605
6. Udomsinprasert W, Nontawong N, Saengsiwaritt W, Panthan B, Jiaranai P, Thongchompoo N, et al. Host genetic polymorphisms involved in long-term symptoms of COVID-19. Emerg Microbes Infect. (2023) 12:2239952. doi: 10.1080/22221751.2023.2239952
7. Elnagdy MH, Magdy A, Eldars W, lgamal M, El-Nagdy AH, Salem O, et al. Genetic association of ACE2 and TMPRSS2 polymorphisms with COVID-19 severity; a single centre study from Egypt. Virol J. (2024) 21:27. doi: 10.1186/s12985-024-02298-x
8. Uslu K, Ozcelik F, Zararsiz G, ldem V, Cephe A, Sahin IO, et al. Deciphering the host genetic factors conferring susceptibility to severe COVID-19 using exome sequencing. Genes Immun. (2024) 25:14–42. doi: 10.1038/s41435-023-00232-9
9. Watson A, Madsen J, Clark HW. SP-A and SP-D: dual functioning immune molecules with antiviral and immunomodulatory properties. Front Immunol. (2020) 11:622598. doi: 10.3389/fimmu.2020.622598
10. Jacob IB, Gemmiti A, Xiong W, Reynolds E, Nicholas B, Thangamani S, et al. Human surfactant protein A inhibits SARS-CoV-2 infectivity and alleviates lung injury in a mouse infection model. Front Immunol. (2024) 15:1370511. doi: 10.3389/fimmu.2024.1370511
11. Funk CJ, Wang J, Ito Y, Travanty EA, Voelker DR, Holmes KV, et al. Infection of human alveolar macrophages by human coronavirus strain 229E. J Gen Virol. (2012) 93:494–503. doi: 10.1099/vir.0.038414-0
12. Gaiha GD, Dong T, Palaniyar N, Mitchell DA, Reid KB, Clark HW. Surfactant protein A binds to HIV and inhibits direct infection of CD4+ cells, but enhances dendritic cell-mediated viral transfer. J Immunol. (2008) 181:601–9. doi: 10.4049/jimmunol.181.1.601
13. LeVine AM, Gwozdz J, Stark J, Bruno M, Whitsett J, Korfhagen T. Surfactant protein-A enhances respiratory syncytial virus clearance. vivo J Clin Invest. (1999) 103:1015–21. doi: 10.1172/JCI5849
14. Benne CA, Kraaijeveld CA, van Strijp JA, Brouwer E, Harmsen M, Verhoef J, et al. Interactions of surfactant protein A with influenza A viruses: binding and neutralization. J Infect Dis. (1995) 171:335–41. doi: 10.1093/infdis/171.2.335
15. Jacob IB, Gemmiti A, Xiong W, Reynolds E, Nicholas B, Thangamani S, et al. Human surfactant protein A alleviates SARS-coV-2 infectivity in human lung epithelial cells. bioRxiv. (2023). doi: 10.1101/2023.04.03.535215
16. Kase T, Suzuki Y, Kawai T, Sakamoto T, Ohtani K, Eda S, et al. Human mannan-binding lectin inhibits the infection of influenza A virus without complement. Immunology. (1999) 97:385–92. doi: 10.1046/j.1365-2567.1999.00781.x
17. Ip WK, Chan KH, Law HK, Tso GH, Kong EK, Wong WH, et al. Mannose-binding lectin in severe acute respiratory syndrome coronavirus infection. J Infect Dis. (2005) 191:1697–704. doi: 10.1086/429631
18. Hsieh MH, Beirag N, Murugaiah V, Chou YC, Kuo WS, Kao HF, et al. Human surfactant protein D binds spike protein and acts as an entry inhibitor of SARS-coV-2 pseudotyped viral particles. Front Immunol. (2021) 12:641360. doi: 10.3389/fimmu.2021.641360
19. Madan T, Biswas B, Varghese PM, Subedi R, Pandit H, Idicula-Thomas S, et al. A recombinant fragment of human surfactant protein D binds spike protein and inhibits infectivity and replication of SARS-coV-2 in clinical samples. Am J Respir Cell Mol Biol. (2021) 65:41–53. doi: 10.1165/rcmb.2021-0005OC
20. Silveyra P, Floros J. Genetic complexity of the human surfactant-associated proteins SP-A1 and SP-A2. Gene. (2013) 531:126–32. doi: 10.1016/j.gene.2012.09.111
21. Oberley RE, Snyder JM. Recombinant human SP-A1 and SP-A2 proteins have different carbohydrate-binding characteristics. Am J Physiol Lung Cell Mol Physiol. (2003) 284:L871–881. doi: 10.1152/ajplung.00241.2002
22. Mikerov AN, Wang G, Umstead TM, Zacharatos M, Thomas NJ, Phelps DS, et al. Surfactant protein A2 (SP-A2) variants expressed in CHO cells stimulate phagocytosis of Pseudomonas aeruginosa more than do SP-A1 variants. Infect Immun. (2007) 75:1403–12. doi: 10.1128/IAI.01341-06
23. Wang G, Phelps DS, Umstead TM, Floros J. Human SP-A protein variants derived from one or both genes stimulate TNF-alpha production in the THP-1 cell line. Am J Physiol Lung Cell Mol Physiol. (2000) 278:L946–954. doi: 10.1152/ajplung.2000.278.5.L946
24. Floros J, Phelps DS. Is the role of lung innate immune molecules, SP-A1 and SP-A2, and of the alveolar macrophage being overlooked in COVID-19 diverse outcomes? Pneumon. (2020) 33(2):1–5.
25. Garcia-Verdugo I, Wang G, Floros J, Casals C. Structural analysis and lipid-binding properties of recombinant human surfactant protein a derived from one or both genes. Biochemistry. (2002) 41:14041–53. doi: 10.1021/bi026540l
26. Floros J, Lin HM, Garcia A, Salazar MA, Guo X, DiAngelo S, et al. Surfactant protein genetic marker alleles identify a subgroup of tuberculosis in a Mexican population. J Infect Dis. (2000) 182:1473–8. doi: 10.1086/jid.2000.182.issue-5
27. Seifart C, Lin HM, Seifart U, Plagens A, DiAngelo S, von Wichert P, et al. Rare SP-A alleles and the SP-A1-6A(4) allele associate with risk for lung carcinoma. Clin Genet. (2005) 68:128–36. doi: 10.1111/j.1399-0004.2005.00470.x
28. Herrera-Ramos E, López-Rodríguez M, Ruíz-Hernández JJ, Horcajada JP, Borderias L, Lerma E, et al. Surfactant protein A genetic variants associate with severe respiratory insufficiency in pandemic influenza A virus infection. Crit Care. (2014) 18:R127. doi: 10.1186/cc13934
29. Thomas NJ, DiAngelo S, Hess JC, Fan R, Ball MW, Geskey JM, et al. Transmission of surfactant protein variants and haplotypes in children hospitalized with respiratory syncytial virus. Pediatr Res. (2009) 66:70–3. doi: 10.1203/PDR.0b013e3181a1d768
30. Zhang H, Zhou G, Zhi L, Yang H, Zhai Y, Dong X, et al. Association between mannose-binding lectin gene polymorphisms and susceptibility to severe acute respiratory syndrome coronavirus infection. J Infect Dis. (2005) 192:1355–61. doi: 10.1086/491479
31. Hultstrom M, Frithiof R, Grip J, Lindelof L, Rooijackers O, Pigazzini S, et al. Genetic determinants of mannose-binding lectin activity predispose to thromboembolic complications in critical COVID-19. Nat Immunol. (2022) 23:861–4. doi: 10.1038/s41590-022-01227-w
32. Garcia-Laorden MI, Rodriguez de Castro F, Sole-Violan J, Rajas O, Blanquer J, Borderias L, et al. Influence of genetic variability at the surfactant proteins A and D in community-acquired pneumonia: a prospective, observational, genetic study. Crit Care. (2011) 15:R57. doi: 10.1186/cc10030
33. Islam A, Khan MA. Lung transcriptome of a COVID-19 patient and systems biology predictions suggest impaired surfactant production which may be druggable by surfactant therapy. Sci Rep. (2020) 10:19395. doi: 10.1038/s41598-020-76404-8
34. Initiative C-HG. A first update on mapping the human genetic architecture of COVID-19. Nature. (2022) 608:E1–E10. doi: 10.1038/s41586-022-04826-7
35. Initiative C-HG. A second update on mapping the human genetic architecture of COVID-19. Nature. (2023) 621:E7–E26. doi: 10.1038/s41586-023-06355-3
36. Piva S, DiBlasi RM, Slee AE, Jobe AH, Roccaro AM, Filippini M, et al. Surfactant therapy for COVID-19 related ARDS: a retrospective case-control pilot study. Respir Res. (2021) 22:20. doi: 10.1186/s12931-020-01603-w
37. Li D, Wang X, Liao Y, Wang S, Shan J, Ji J. Insights gained into the treatment of COVID19 by pulmonary surfactant and its components. Front Immunol. (2022) 13:842453. doi: 10.3389/fimmu.2022.842453
38. Wang G, Guo X, Diangelo S, Thomas NJ, Floros J. Humanized SFTPA1 and SFTPA2 transgenic mice reveal functional divergence of SP-A1 and SP-A2: formation of tubular myelin in vivo requires both gene products. J Biol Chem. (2010) 285:11998–2010. doi: 10.1074/jbc.M109.046243
39. Oladunni FS, Park JG, Pino PA, Gonzalez O, Akhter A, Allue-Guardia A, et al. Lethality of SARS-CoV-2 infection in K18 human angiotensin-converting enzyme 2 transgenic mice. Nat Commun. (2020) 11:6122. doi: 10.1038/s41467-020-19891-7
40. Matute-Bello G, Downey G, Moore BB, Groshong SD, Matthay MA, Slutsky AS, et al. An official American Thoracic Society workshop report: features and measurements of experimental acute lung injury in animals. Am J Respir Cell Mol Biol. (2011) 44:725–38. doi: 10.1165/rcmb.2009-0210ST
41. Winkler ES, Bailey AL, Kafai NM, Nair S, McCune BT, Yu J, et al. SARS-CoV-2 infection of human ACE2-transgenic mice causes severe lung inflammation and impaired function. Nat Immunol. (2020) 21:1327–35. doi: 10.1038/s41590-020-0778-2
42. Liao WY, Li HJ, Chang MY, Tang AC, Hoffman AS, Hsieh PC. Comprehensive characterizations of nanoparticle biodistribution following systemic injection in mice. Nanoscale. (2013) 5:11079–86. doi: 10.1039/c3nr03954d
43. Yang F, Yang Y, Ruan F, Chen X, Guo J, Abdel-Razek O, et al. Regulatory roles of human surfactant protein B variants on genetic susceptibility to pseudomonas aeruginosa pneumonia-induced sepsis. Shock. (2020) 54:507–19. doi: 10.1097/SHK.0000000000001494
44. Silveyra P, Floros J. Genetic variant associations of human SP-A and SP-D with acute and chronic lung injury. Front Biosci (Landmark Ed). (2012) 17:407–29. doi: 10.2741/3935
45. Gandhi CK, Chen C, Wu R, Yang L, Thorenoor N, Thomas NJ, et al. Association of SNP-SNP interactions of surfactant protein genes with pediatric acute respiratory failure. J Clin Med. (2020) 9(4):1183. doi: 10.3390/jcm9041183
46. Thorenoor N, Phleps SD, Kala P, Ravi R, Phelps AF, Umstead TM. Impact of surfactant protein-A variants on survival in aged mice in response to klebsiella pneumoniae infection and ozone: serendipity in action. Microorganisms. (2020) 8(9):1276. doi: 10.3390/microorganisms8091276
47. Thorenoor N, Zhang X, Umstead TM, Scott Halstead E, Phelps DS, Floros J. Differential effects of innate immune variants of surfactant protein-A1 (SFTPA1) and SP-A2 (SFTPA2) in airway function after Klebsiella pneumoniae infection and sex differences. Respir Res. (2018) 19:23. doi: 10.1186/s12931-018-0723-1
48. Knudtzen FC, Jensen TG, Lindvig SO, Rasmussen LD, Madsen LW, Hoegh SV, et al. SARS-CoV-2 viral load as a predictor for disease severity in outpatients and hospitalised patients with COVID-19: A prospective cohort study. PloS One. (2021) 16:e0258421. doi: 10.1371/journal.pone.0258421
49. Fajnzylber J, Regan J, Coxen K, Corry H, Wong C, Rosenthal A, et al. SARS-CoV-2 viral load is associated with increased disease severity and mortality. Nat Commun. (2020) 11:5493. doi: 10.1038/s41467-020-19057-5
50. Nayak A, Dodagatta-Marri E, Tsolaki AG, Kishore U. An insight into the diverse roles of surfactant proteins, SP-A and SP-D in innate and adaptive immunity. Front Immunol. (2012) 3:131. doi: 10.3389/fimmu.2012.00131
51. Bizzotto J, Sanchis P, Abbate M, Lage-Vickers S, Lavignolle R, Toro A, et al. SARS-coV-2 infection boosts MX1 antiviral effector in COVID-19 patients. iScience. (2020) 23:101585. doi: 10.1016/j.isci.2020.101585
52. Furci F, Murdaca G, Allegra A, Gammeri L, Senna G, Gangemi S. IL-33 and the cytokine storm in COVID-19: from a potential immunological relationship towards precision medicine. Int J Mol Sci. (2022) 23(23):14532. doi: 10.3390/ijms232314532
53. Gao Y, Cai L, Li L, Zhang Y, Li J, Luo C, et al. Emerging effects of IL-33 on COVID-19. Int J Mol Sci. (2022) 23(21):13656. doi: 10.3390/ijms232113656
54. Markovic SS, Jovanovic M, Gajovic N, Jurisevic M, Arsenijevic N, Jovanovic M, et al. IL 33 correlates with COVID-19 severity, radiographic and clinical finding. Front Med (Lausanne). (2021) 8:749569. doi: 10.3389/fmed.2021.749569
55. Wang J, Yang X, Li Y, Huang JA, Jiang J, Su N. Specific cytokines in the inflammatory cytokine storm of patients with COVID-19-associated acute respiratory distress syndrome and extrapulmonary multiple-organ dysfunction. Virol J. (2021) 18:117. doi: 10.1186/s12985-021-01588-y
56. Luporini RL, Rodolpho JMA, Kubota LT, Martin A, Cominetti MR, Anibal FF, et al. IL-6 and IL-10 are associated with disease severity and higher comorbidity in adults with COVID-19. Cytokine. (2021) 143:155507. doi: 10.1016/j.cyto.2021.155507
57. Liu Y, Zhang C, Huang F, Yang Y, Wang F, Yuan J, et al. Elevated plasma levels of selective cytokines in COVID-19 patients reflect viral load and lung injury. Natl Sci Rev. (2020) 7:1003–11. doi: 10.1093/nsr/nwaa037
58. Chen G, Wu D, Guo W, Cao Y, Huang D, Wang H, et al. Clinical and immunological features of severe and moderate coronavirus disease 2019. J Clin Invest. (2020) 130:2620–9. doi: 10.1172/JCI137244
59. Liu J, Li S, Liu J, Liang B, Wang X, Wang H, et al. Longitudinal characteristics of lymphocyte responses and cytokine profiles in the peripheral blood of SARS-CoV-2 infected patients. EBioMedicine. (2020) 55:102763. doi: 10.1016/j.ebiom.2020.102763
60. Hilligan KL, Namasivayam S, Clancy CS, Baker PJ, Old SI, Peluf V, et al. Bacterial-induced or passively administered interferon gamma conditions the lung for early control of SARS-CoV-2. Nat Commun. (2023) 14:8229. doi: 10.1038/s41467-023-43447-0
61. Islam MM, Islam S, Ahmed R, Majumder M, Sarkar B, Himu MER, et al. Reduced IFN-γ levels along with changes in hematologic and immunologic parameters are key to COVID-19 severity in Bangladeshi patients. Exp Hematol. (2023) 118:53–64.e51. doi: 10.1016/j.exphem.2022.11.006
62. van der Made CI, Netea MG, van der Veerdonk FL, Hoischen A. Clinical implications of host genetic variation and susceptibility to severe or critical COVID-19. Genome Med. (2022) 14:96. doi: 10.1186/s13073-022-01100-3
63. Lowe KE, Zein J, Hatipoglu U, Attaway A. Association of smoking and cumulative pack-year exposure with COVID-19 outcomes in the cleveland clinic COVID-19 registry. JAMA Intern Med. (2021) 181:709–11. doi: 10.1001/jamainternmed.2020.8360
64. Bucholc M, Bradley D, Bennett D, Patterson L, Spiers R, Gibson D, et al. Identifying pre-existing conditions and multimorbidity patterns associated with in-hospital mortality in patients with COVID-19. Sci Rep. (2022) 12:17313. doi: 10.1038/s41598-022-20176-w
65. Dhindsa S, Zhang N, McPhaul MJ, Wu Z, Ghoshal AK, Erlich EC, et al. Association of circulating sex hormones with inflammation and disease severity in patients with COVID-19. JAMA Netw Open. (2021) 4:e2111398. doi: 10.1001/jamanetworkopen.2021.11398
66. Velavan TP, Pallerla SR, Ruter J, Augustin Y, Kremsner PG, Krishna S, et al. Host genetic factors determining COVID-19 susceptibility and severity. EBioMedicine. (2021) 72:103629. doi: 10.1016/j.ebiom.2021.103629
67. Barnkob MB, Pottegard A, Stovring H, Haunstrup TM, Homburg K, Larsen R, et al. Reduced prevalence of SARS-CoV-2 infection in ABO blood group O. Blood Adv. (2020) 4:4990–3. doi: 10.1182/bloodadvances.2020002657
68. Ellinghaus D. COVID-19 host genetics and ABO blood group susceptibility. Camb Prism Precis Med. (2023) 1:e10. doi: 10.1017/pcm.2022.12
69. Severe Covid GG, Ellinghaus D, Degenhardt F, Bujanda L, Buti M, Albillos A, et al. Genomewide association study of severe covid-19 with respiratory failure. N Engl J Med. (2020) 383:1522–34. doi: 10.1056/NEJMoa2020283
70. Choudhary S, Sreenivasulu K, Mitra P, Misra S, Sharma P. Role of genetic variants and gene expression in the susceptibility and severity of COVID-19. Ann Lab Med. (2021) 41:129–38. doi: 10.3343/alm.2021.41.2.129
71. Du J, Abdel-Razek O, Shi Q, Hu F, Ding G, Cooney RN, et al. Surfactant protein D attenuates acute lung and kidney injuries in pneumonia-induced sepsis through modulating apoptosis, inflammation and NF-κB signaling. Sci Rep. (2018) 8:15393. doi: 10.1038/s41598-018-33828-7
72. LeVine AM, Hartshorn K, Elliott J, Whitsett J, Korfhagen T. Absence of SP-A modulates innate and adaptive defense responses to pulmonary influenza infection. Am J Physiol Lung Cell Mol Physiol. (2002) 282:L563–572. doi: 10.1152/ajplung.00280.2001
73. Al-Qahtani AA, Murugaiah V, Bashir HA, Pathan AA, Abozaid SM, Makarov E, et al. Full-length human surfactant protein A inhibits influenza A virus infection of A549 lung epithelial cells: A recombinant form containing neck and lectin domains promotes infectivity. Immunobiology. (2019) 224:408–18. doi: 10.1016/j.imbio.2019.02.006
74. Hsieh IN, De Luna X, White MR, Hartshorn KL. The role and molecular mechanism of action of surfactant protein D in innate host defense against influenza A virus. Front Immunol. (2018) 9:1368. doi: 10.3389/fimmu.2018.01368
75. Stravalaci M, Pagani I, Paraboschi EM, Pedotti M, Doni A, Scavello F, et al. Recognition and inhibition of SARS-CoV-2 by humoral innate immunity pattern recognition molecules. Nat Immunol. (2022) 23:275–86. doi: 10.1038/s41590-021-01114-w
76. Li G, Siddiqui J, Hendry M, Akiyama J, Edmondson J, Brown C, et al. Surfactant protein-A–deficient mice display an exaggerated early inflammatory response to a beta-resistant strain of influenza A virus. Am J Respir Cell Mol Biol. (2002) 26:277–82. doi: 10.1165/ajrcmb.26.3.4584
77. Floros J, DiAngelo S, Koptides M, Karinch AM, Rogan PK, Nielsen H, et al. Human SP-A locus: allele frequencies and linkage disequilibrium between the two surfactant protein A genes. Am J Respir Cell Mol Biol. (1996) 15:489–98. doi: 10.1165/ajrcmb.15.4.8879183
78. DiAngelo S, Lin Z, Wang G, Phillips S, Ramet M, Luo J, et al. Novel, non-radioactive, simple and multiplex PCR-cRFLP methods for genotyping human SP-A and SP-D marker alleles. Dis Markers. (1999) 15:269–81. doi: 10.1155/1999/961430
79. Wang G, Umstead TM, Hu S, Mikerov AN, Phelps DS, Floros J. Differential effects of human SP-A1 and SP-A2 on the BAL proteome and signaling pathways in response to klebsiella pneumoniae and ozone exposure. Front Immunol. (2019) 10:561. doi: 10.3389/fimmu.2019.00561
80. Casals C, Campanero-Rhodes MA, Garcia-Fojeda B, Solis D. The role of collectins and galectins in lung innate immune defense. Front Immunol. (2018) 9:1998. doi: 10.3389/fimmu.2018.01998
81. Hadjadj J, Yatim N, Barnabei L, Corneau A, Boussier J, Smith N, et al. Impaired type I interferon activity and inflammatory responses in severe COVID-19 patients. Science. (2020) 369:718–24. doi: 10.1126/science.abc6027
82. Sefik E, Qu R, Junqueira C, Kaffe E, Mirza H, Zhao J, et al. Inflammasome activation in infected macrophages drives COVID-19 pathology. Nature. (2022) 606:585–93. doi: 10.1038/s41586-022-04802-1
83. Zheng M, Karki R, Williams EP, Yang D, Fitzpatrick E, Vogel P, et al. TLR2 senses the SARS-CoV-2 envelope protein to produce inflammatory cytokines. Nat Immunol. (2021) 22:829–38. doi: 10.1038/s41590-021-00937-x
84. Thorenoor N, Kawasawa YI, Gandhi CK, Zhang X, Floros J. Differential impact of co-expressed SP-A1/SP-A2 protein on AM miRNome; sex differences. Front Immunol. (2019) 10:1960. doi: 10.3389/fimmu.2019.01960
85. Guillot L, Balloy V, McCormack FX, Golenbock DT, Chignard M, Si-Tahar M. Cutting edge: the immunostimulatory activity of the lung surfactant protein-A involves Toll-like receptor 4. J Immunol. (2002) 168:5989–92. doi: 10.4049/jimmunol.168.12.5989
86. Koptides M, Umstead TM, Floros J, Phelps DS. Surfactant protein A activates NF-kappa B in the THP-1 monocytic cell line. Am J Physiol. (1997) 273:L382–388. doi: 10.1152/ajplung.1997.273.2.L382
87. Yamada C, Sano H, Shimizu T, Mitsuzawa H, Nishitani C, Himi T, et al. Surfactant protein A directly interacts with TLR4 and MD-2 and regulates inflammatory cellular response. Importance of supratrimeric oligomerization. J Biol Chem. (2006) 281:21771–80. doi: 10.1074/jbc.M513041200
88. Henning LN, Azad AK, Parsa KV, Crowther JE, Tridandapani S, Schlesinger LS. Pulmonary surfactant protein A regulates TLR expression and activity in human macrophages. J Immunol. (2008) 180:7847–58. doi: 10.4049/jimmunol.180.12.7847
89. Mohamadzadeh M, Coberley SS, Olinger GG, Kalina WV, Ruthel G, Fuller CL, et al. Activation of triggering receptor expressed on myeloid cells-1 on human neutrophils by marburg and ebola viruses. J Virol. (2006) 80:7235–44. doi: 10.1128/JVI.00543-06
90. Zhong WJ, Duan JX, Liu T, Yang HH, Guan XX, Zhang CY, et al. Activation of NLRP3 inflammasome up-regulates TREM-1 expression in murine macrophages via HMGB1 and IL-18. Int Immunopharmacol. (2020) 89:107045. doi: 10.1016/j.intimp.2020.107045
91. Sharapova TN, Romanova EA, Ivanova OK, Yashin DV, Sashchenko LP. Hsp70 interacts with the TREM-1 receptor expressed on monocytes and thereby stimulates generation of cytotoxic lymphocytes active against MHC-negative tumor cells. Int J Mol Sci. (2021) 22(13):6889. doi: 10.3390/ijms22136889
92. Minutti CM, García-Fojeda B, Sáenz A, de Las Casas-Engel M, Guillamat-Prats R, de Lorenzo A, et al. Surfactant protein A prevents IFN-γ/IFN-γ Receptor interaction and attenuates classical activation of human alveolar macrophages. J Immunol. (2016) 197:590–8. doi: 10.4049/jimmunol.1501032
93. Lopes-Pacheco M, Silva PL, Cruz FF, Battaglini D, Robba C, Pelosi P, et al. Pathogenesis of multiple organ injury in COVID-19 and potential therapeutic strategies. Front Physiol. (2021) 12:593223. doi: 10.3389/fphys.2021.593223
94. Akamatsu MA, de Castro JT, Takano CY, Ho PL. Off balance: Interferons in COVID-19 lung infections. EBioMedicine. (2021) 73:103642. doi: 10.1016/j.ebiom.2021.103642
95. Choi CY, Gadhave K, Villano J, Pekosz A, Mao X, Jia H. Generation and characterization of a humanized ACE2 mouse model to study long-term impacts of SARS-CoV-2 infection. J Med Virol. (2024) 96:e29349. doi: 10.1002/jmv.29349
Keywords: humanized transgenic mice, innate immunity, SARS-CoV-2 infection, signaling pathway, surfactant protein A, susceptibility
Citation: Jacob IB, Lawal AO, Mahmoud SS, Kopsack EM, Reynolds ES, Meng Q, Fan H, Massa PT, Thangamani S, Jia H and Wang G (2025) Differential immunoregulation by human surfactant protein A variants determines severity of SARS-CoV-2-induced lung disease. Front. Immunol. 16:1462278. doi: 10.3389/fimmu.2025.1462278
Received: 09 July 2024; Accepted: 10 March 2025;
Published: 02 April 2025.
Edited by:
Antigoni Kotsaki, University General Hospital Attikon, GreeceReviewed by:
Enrique Roberto Argañaraz, University of Brasilia, BrazilSung-Hee Kim, Yonsei University, Republic of Korea
Copyright © 2025 Jacob, Lawal, Mahmoud, Kopsack, Reynolds, Meng, Fan, Massa, Thangamani, Jia and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Guirong Wang, wangg@upstate.edu