Endothelial cell infection and dysfunction, immune activation in severe COVID-19.

Rationale: Pulmonary vascular endotheliitis, perivascular inflammation, and immune activation are observed in COVID-19 patients. While the initial SARS-CoV-2 infection mainly infects lung epithelial cells, whether it also infects endothelial cells (ECs) and to what extent SARS-CoV-2-mediated pulmonary vascular endotheliitis is associated with immune activation remain to be determined.
Methods: To address these questions, we studied SARS-CoV-2-infected K18-hACE2 (K18) mice, a severe COVID-19 mouse model, as well as lung samples from SARS-CoV-2-infected nonhuman primates (NHP) and patient deceased from COVID-19. We used immunostaining, RNAscope, and electron microscopy to analyze the organs collected from animals and patient. We conducted bulk and single cell (sc) RNA-seq analyses, and cytokine profiling of lungs or serum of the severe COVID-19 mice.
Results: We show that SARS-CoV-2-infected K18 mice develop severe COVID-19, including progressive body weight loss and fatality at 7 days, severe lung interstitial inflammation, edema, hemorrhage, perivascular inflammation, systemic lymphocytopenia, and eosinopenia. Body weight loss in K18 mice correlated with the severity of pneumonia, but not with brain infection. We also observed endothelial activation and dysfunction in pulmonary vessels evidenced by the up-regulation of VCAM1 and ICAM1 and the downregulation of VE-cadherin. We detected SARS-CoV-2 in capillary ECs, activation and adhesion of platelets and immune cells to the vascular wall of the alveolar septa, and increased complement deposition in the lungs, in both COVID-19-murine and NHP models. We also revealed that pathways of coagulation, complement, K-ras signaling, and genes of ICAM1 and VCAM1 related to EC dysfunction and injury were upregulated, and were associated with massive immune activation in the lung and circulation.
Conclusion: Together, our results indicate that SARS-CoV-2 causes endotheliitis via both infection and infection-mediated immune activation, which may contribute to the pathogenesis of severe COVID-19 disease. © The author(s).

Introduction

SARS-CoV-2 has killed more than one million people globally. However, we have not yet developed highly effective therapeutics for treatment of the COVID-19 disease 1. This is partially due to a limited understanding of COVID-19 pathogenesis. COVID-19 patients display a wide range of symptoms including severe illness associated with disseminated platelet activation and coagulation, cytokine storm, multiple organ failure, and death. The presence of widespread microthrombi and endothelial cell (EC) damage throughout the pulmonary vasculature suggests that vasculopathy is important in COVID-19 pathogenesis 2, 3. A better understanding of these pathogenic mechanisms offers important opportunities to identify novel agents or repurpose existing drugs to attenuate viral attachment and replication, inflammation, EC damage, blood clotting and systemic microangiopathy 4, 5.

SARS-CoV-2 primarily infects pneumocytes expressing the SARS-CoV-2 receptor ACE2, which facilitates viral entry. Little is known about SARS-CoV-2-mediated direct or indirect vasculopathy and its possible role in disease outcomes 6, 7. Studies indicate that in a variety of tissues, and especially in the lung, ECs may be infected by the virus 2, 8-11, although they express lower ACE2 levels compared to pneumocytes 1. Direct EC infection may explain many clinical features of advanced disease, including disseminated activation of coagulation pathways and platelets and propensity to thrombosis, especially in the lung. It has been hypothesized that COVID-19-related acute respiratory distress syndrome (ARDS) is a vascular endotype or even an EC disease 12, 13. However, whether SARS-CoV-2 directly infects microvascular EC in organs, including the lung, remains uncertain 14. In addition, the identification of SARS-CoV-2 in human EC by electron microscopy (EM) has been demonstrated but remains a subject of scientific debate 11, 15-17. The direct infection of EC by SARS-CoV-2 in the lung and other organs have not been investigated experimentally 1. Further, using the COVID-19 rhesus macaques' model, previous study indicates that SARS-CoV-2-mediated indirect endothelial disruption and vascular thrombosis via interactions between immune and inflammatory pathways contribute to SARS-CoV-2-induced EC dysfunction 18. However, rhesus macaques do not develop severe COVID-19, and thus cannot provide information on the pathogenesis of severe disease in humans 18. Also, most of clinical studies showing EC infection and dysfunction were conducted on autopsied human samples, which makes it difficult to identify the pathogenic role of EC infection and dysfunction and the underlying mechanisms of alterations observed. Together, the relative contribution of SARS-CoV-2-mediated direct EC infection and indirect immune-mediated EC dysfunction in COVID-19 still requires investigation 5.

To address these questions, we used a mouse model of severe COVID-19, SARS-CoV-2-infected K18-hACE2 (K18) mice, as well as tissue from SARS-CoV-2-infected NHPs and samples from autopsy of COVID-19 patient. We show that, in addition to infecting pneumocytes, SARS-CoV-2 also infects EC. We observe that SARS-CoV-2-mediated EC dysfunction is associated with immune, complement and platelet activation/thrombosis in the lung of the K18 mice and NHPs. Our results suggest that SARS-CoV-2-mediated direct EC infection and indirect EC dysfunction may contribute to the pathogenesis of COVID-19. These results also indicate that pharmacologic suppression of EC infection and dysfunction may have a beneficial effect in treating or reducing the severity of COVID-19 and possibly diseases associated with emerging coronaviruses.

Results

SARS-CoV-2-infected K18 mice developed severe COVID-19

To investigate the potential roles of immune activation, EC infection, and dysfunction underlying severe COVID-19, we established and characterized SARS-CoV-2-infected K18 mice. K18 transgenic mice express hACE2 in airway epithelial cells under the control of the human cytokeratin 18 (K18) promoter 19. We infected K18 mice with a moderate dose of SARS-CoV-2 (2 X 105 TCID50) and euthanized the mice at 4-, 5-, 6-, and 7-days post infection (DPI). Consistent with previous findings 20-24, we found that SARS-CoV-2 infection of K18 mice induced a progressive body weight loss from 4 to 7 DPI and fatality at 7 DPI, compared with infected non-transgenic wild type (WT) mice (Figure ). The mice developed mild to severe interstitial pneumonia with progressive perivascular inflammation (Figure ), alveolar edema (Figure ), and hemorrhage (Figure ) at 4 to 7 DPI. Of note, these mice developed pathological changes with various degrees of edema, hemorrhage, or inflammation at different time intervals (Figure ), and these histological changes were similar to those seen in severe COVID-19 patients 25, 26. After 5 DPI, mice exhibited at least one of these histological changes, all of which contribute to the severe COVID-19 phenotype. Further, in addition to the lung histopathology described above, K18 mice euthanized at 6 DPI also exhibited some degree of the pulmonary thrombosis (Figure ). We then stained the lung sections for P-selectin and CD41, well established biomarkers of platelet activation and aggregation 27-29. Severe to mild pulmonary platelet activation was observed by P-selectin (Figure and ) and CD41 staining of the lungs of the infected K18 mice ().

Additionally, infected K18 mice had fewer white blood cells, including lymphocytes and eosinopenia compared to infected WT (or K18) mice at either 4 or 6 DPI (). Reduction in these cell counts has been linked to severe COVID-19 in humans 30, 31. Serum chemistry analysis showed increased levels of LDH, a routinely used clinical biomarker for organ damage 32, as well as AST and ALT, indicative of liver injury (). Such changes in serum chemistry analytes suggest involvement of other organs, especially the liver. We interpreted the decrease in blood glucose and triglyceride seen in infected-K18 mice compared to controls as the consequence of body weight loss in these mice ().

To examine SARS-CoV-2 infection in tissues of K18 mice, we used qRT-PCR to quantify the nucleocapsid (N) sub-genomic viral RNA, a correlate of active virus replication 33, 34, in organs daily until 6 DPI (). We found that the lungs of SARS-CoV-2-infected K18 mice had, on average, 105 (at 1 DPI), 109 (at 2 and 3 DPI), and 106 (at 4- to 6- DPI) sub-genomic N RNA copies per 100ng of tissue RNA. This dynamic growth of replicating virus was further confirmed by direct viral staining, a reliable and sensitive method for detecting replicating virus in tissues 34 (). Further studies revealed that the brains of infected K18 mice had 103 (at 1, 2, 3 DPI), and 108 (at 4-6 DPI) subgenomic viral RNA load (), suggesting that the lung had been infected earlier than the brain. In addition, SARS-CoV-2-infected K18 mice exhibited several logs higher genomic viral load in the liver but not in the spleen, compared to the infected nontransgenic littermate (or WT) mice (). It is notable that viral loads of infected K18 mice were much higher in lung and brain than in spleen and liver ().

Together, these findings indicate that K18 mice infected with a moderate dose of SARS-CoV-2 primarily results in infection of the lung and brain. The mice developed severe COVID-19 associated with severe lung histological changes.

Lung disease is the primary determinant of COVID-19 severity in K18 mice

The infected K18 mice also developed progressive body weight loss. Whether this body weight loss was related to lung or brain infection remains unclear 35, and this question is relevant to the selection of COVID-19 model for pathogenesis studies. To address this question, we analyzed a cohort of five infected K18 mice including two mice that died at 7 DPI (Figure ). Histologic analysis revealed moderate postmortem autolysis with pulmonary perivascular infiltration by low to moderate numbers of lymphocytes, widespread edema, multifocal neutrophilic inflammation (Figure left panel), and congested blood vessel (Figure right panel). We did not observe obvious brain lesions either in this mouse (Figure ) or the other three mice euthanized at 7 DPI (Figure ). We excluded another mouse that had died from the analysis due to severe postmortem autolysis. Direct staining of the brain for viral antigens was used to assess the severity of brain infections among these three mice. While the three mice had lost 16.6, 21.6, and 25.6% of their body weight respectively, body weight loss did not correlate with the severity of the infection within the brain (Figure ). These results suggest that lung infection is the primary determinant of the body weight loss in this mouse model.

We further explored SARS-CoV-2-mediated damage in other organs of the three mice euthanized at 7 DPI and one mouse that died at 7 DPI. We observed organ damage in 3 of 4 mice. One mouse exhibited multiple organ damage including cerebral infarction, minimal myocarditis, minimal splenic congestion and hepatocellular vacuolation (). Another mouse exhibited minimal myocarditis ( and the third mouse exhibited minimal hepatitis (). Our results demonstrate that severe COVID-19 in K18 mice results in damage to organs other than the lung as seen in severe COVID-19 patients, but the lung injury remains the primary determinant of the clinical manifestations of severe COVID-19 including pneumonia, body weight loss, and fatality.

SARS-CoV-2 infection causes EC activation and dysfunction

Our experimental results described above suggest that SARS-CoV-2 infection in K18 mice causes severe COVID-19 phenotype. In addition to the typical lung histological changes described above, perivascular pulmonary inflammation has been reported in many COVID-19 animal models, including nonhuman primates 36, mice 37, hamsters 38, and ferrets 35. We also observed the same perivascular inflammation in our severe COVID-19 mouse model (Figure ). Endotheliitis has been previously shown in the endothelium and perivascular space of SARS-CoV-2-infected humans and NHPs, and is previously characterized by hyperplastic endothelium, intimal proliferation and adherence of macrophages and lymphocytes 3, 18, 39. Accordingly, we found that SARS-CoV-2 infected K18 mice also exhibited prominent endotheliitis confirmed by histological studies (Figure ). The detection of increased intercellular adhesion molecule 1 (ICAM-1) and vascular cell adhesion molecule 1 (VCAM1) by immunostaining is widely used to assess EC activation and dysfunction 40. We found that infected K18 but not WT mice had dramatically increased co-localization of VCAM1 and ICAM-1 with the EC marker CD31 in the pulmonary vessels (Figure and ), and increased immune cells adherence to pulmonary EC, observed by ultra-structural analysis (Figure ). Additionally, we examined vascular barrier integrity by staining with VE-cadherin, which plays a central role in permeability changes by controlling the opening and closing of the EC barrier 41, 42. We found decreased expression of VE-cadherin in the pulmonary vasculature (Figure and ) and associated increased gap and decreased length of the tight junction in the areas of inter-EC contacts (Figure ) of K18 mice, but not that of non-transgenic littermates infected with SARS-CoV-2, which is indicative of EC barrier disruption and thus increased permeability. These results suggest the existence of extensive EC activation, thrombosis and disrupted EC integrity in K18 mice. Thus, we conclude that SARS-CoV-2 infection in K18 mice induces EC activation and dysfunction, leading to increased vascular permeability.

SARS-CoV-2 infects ECs

Previous studies show that SARS-CoV-2 infects pneumocytes in K18 mice 22, 23 and SARS-CoV-2 N protein was occasionally observed in the EC lining of the pulmonary vessels 23. Whether EC infection occurs in this model remains to be determined. Using immunostaining, we found SARS-CoV-2 co-localized in alveolar septa with CD31, a biomarker for ECs (Figure ) and 15% of SARS-CoV-2-stained cells co-stained with CD31 in the infected K18 mice at 3 DPI (Figure ). Using RNAscope, we also found that SARS-CoV-2 spike RNA co-localized with the CD31 RNA in the alveolar septa (Figure ). EM also revealed large numbers of SARS-CoV-2 viral particles overlying ECs in alveolar capillaries (Figure ). These results indicate that ECs can be infected by SARS-CoV-2. We further investigated the kinetics of the EC infection by performing the co-staining of the lungs of the infected K18 mice at 1 (early course of the disease) and 6 DPI (late course of the disease). We detected the co-staining of SARS-CoV-2 with CD31 at both time intervals although the mice at both times had less viral staining in the lungs compared to the mice at 3 DPI (). This suggests that SARS-CoV-2 infects ECs in an entire course of the disease and implies that direct EC infection may contribute to the EC dysfunction found in COVID-19. Additionally, we also demonstrated the co-localization of SARS-CoV-2 with hACE2 and serine transmembrane protein 2 (Tmprss2) (), an established receptor and recently known co-receptor for the viral infection, respectively 1, 43-46. These results further support the notion that transgenic expression of hACE2 in K18 mice sensitizes the mice to the infection and mouse Tmprss2 may also contribute to the infection 1, 43-46.

SARS-CoV-2-mediated EC activation is also seen in the lungs of SARS-CoV-2-infected NHPs and a human patient

Previous histological analysis revealed the existence of EC barrier disruption and vascular thrombosis in SARS-CoV-2-infected rhesus macaques 18.

We further evaluated EC dysfunction in using our published and archived SARS-CoV-2 infected NHP lungs 47. We also documented extensive VCAM1 expression in the pulmonary vessel of a SARS-CoV-2 infected African green monkey (AGM), NC34, that was euthanized due to ARDS 47 (Figure red arrow in left panel). Other infected AGM (NC33 and NC40) also exhibited significantly higher VCAM1 staining in the lung than naïve monkeys (Figure right panel). Both NC34 and NC33 developed ARDS characterized by fulminant respiratory distress, radiographic lung opacity, and diffuse alveolar damage at 22 DPI 47. NC40 survived to the end of the study (26 DPI) and showed mild interstitial pneumonia 47. Consistently, we observed EC infection evidenced by co-localization of SARS-CoV-2 with CD31 within the alveolar septa (Figure ). We also found EC dysfunction in an autopsied lung collected from a patient who died of severe COVID-19 as evidenced by VCAM1 staining (Figure ), which is associated with extensive edema (Figure ), and viral infection evidenced by extensive viral staining in the lung (Figure ). Together, we also demonstrate EC infection by SARS-CoV-2 and dysfunction mediated by SARS-CoV-2 infection in the lungs of COVID-19 NHP and human samples.

RNA seq, scRNA seq, lung tissue array, and serum cytokine analyses

Bulk RNAseq analyses of SARS-CoV-2 challenged K18 vs WT mice showed up-regulation of 1078 genes (DESeq) and downregulation of 689 genes in K18 mice compared to the infected littermate controls. Lung tissue mRNA array analysis at 3 and 6 PDI showed induction of inflammatory cytokines and chemokines including Cxcl10, Cxcl9, Cxcl13, IL-6, and CCL2 with prominent increases at 6 DPI (). Consistently, we only detected induction of IL-18, CxcL10 and CCL2 in sera collected from the infected K18 mice compared to the infected WT mice at 6 but not at 4 DPI (). Pathway analysis showed induction of both type I and type II interferons, inflammatory response, and TNF-alpha signaling via NF-κB pathways, the IL-6/STAT3 pathway, and the K-ras pathway, as well as complement and coagulation at 4 and 6 DPI in SARS-CoV-2 infected K18 mice compared to the infected WT mice (Figure ). We performed over-enrichment analysis of the genes Icam1, Vcam1 and Cdh5, which are respectively encoded for Icam1, Vcam1, and VE-cadherin proteins, all of which are the biomarker used for detecting EC dysfunction. The analyses revealed increased expression of Icam-1 and Vcam1 and reduced expression of Cdh5 in K18 mice compared to the WT mice after challenge with SARS-CoV-2 (). These results further support the notion that SARS-CoV-2 infection mediates EC dysfunction. We previously observed induction of the interferon-dependent chemokines Cxcl9 and Cxcl10 in the Ad5-ACE2 model 34, as was also observed in the K18 model. Single cell RNA-seq (scRNA-seq) analysis was performed at 4 DPI, which identified 12 distinct clusters of cells (Figure ).

scRNA-seq comparison between SARS-CoV-2 infected K18 and WT mice showed an increase in the Cx3cr1/Ccr2+ macrophage (Mϕ) cluster (), fibroblast cluster, and the neutrophil cluster, as well as IgDlo B cells and a trend towards increased type II pneumocytes (Figure ). We next analyzed data with a dual alignment to the mouse genome and the SARS-CoV-2 genome (WA1/2020 isolate, accession MN985325). TSNE analysis identified the same cell clusters in this algorithm. When mapped to the SARS-CoV-2 genome and analyzed by reclusterring tool in Loupe Browser 5.0, we found viral RNA in the following clusters: epithelial cells, endothelial cells, dendritic cells, Cd68+ macrophages and eosinophils (Figure ). Observation of viral RNA in ECs further supports our findings that the ECs can be infected by SARS-CoV-2. Our data also demonstrated greater expression of Cxcl9 and Cxcl10 in Twist2+ fibroblasts as well as in the EC cluster in SARS-CoV-2 infected K18 mice compared to infected WT mice () indicative of an upregulated interferon pathway. Upregulation of the K-ras pathway was also observed (Figure ), as we detected the expression of several serpin genes in the scRNAseq data (). These analyses revealed that the infected K18 expressed higher Serpina3f and 3n in the fibroblast and EC cluster () and higher Serpina3g was expressed in these cell clusters in addition to the Ccr2+ Mϕ cluster, the B cell cluster, and the T cell cluster () compared to the infected WT mice. Given that there was evidence of interferon dependent chemokine expression in the EC cluster as well as evidence of EC activation (Figures and 5), we analyzed expression of Vcam1 and Icam1 in the scRNA-seq data. Both molecules were expressed and up regulated in the EC and fibroblast clusters as well as the Mϕ/MC cluster in the infected K18 mice (). These results indicate that Icam1 is also expressed in B cells, T cells, neutrophils and Cd11c+ Mϕ (). To further investigate whether specific pathways were regulated in distinct specific cell clusters by the presence of viral RNA, we analyzed viral positive versus negative cells in the endothelial and monocyte cluster. In this analysis, we found upregulation of genes in the lysosome/apoptosis and complement pathways in viral RNA positive cells (). Together, these results shed lights on the importance of SARS-CoV-2-mediated immune activation for the endothelial dysfunction and injury.

Consistent with the bulk RNA-Seq data showing upregulation of the complement pathway, we observed that Ctss, C1qa and C1qb were highly expressed in the alveolar Mϕ cluster as well as the Cxc3r1/Ccr2+ Mϕ/MC cluster in SARS-CoV-2 infected K18 mice compared to infected WT mice (). Cfb was also expressed in the same cells, as well as in the fibroblast cluster (). Consistently, we also found increased C3 and a component of the membrane attack complex (MAC or C5b9) in the lung and ECs of SARS-CoV-2-infected K18 mice and AGMs compared to the infected WT mice and naïve AGMs ().

Discussion

Most current COVID-19 animal models including hamsters, other human ACE2 transgenic mice and non-human primates at young ages, show only a mild to moderate COVID-19 phenotype. We recently reported a severe COVID-19 phenotype in aged AGMs 48. In contrast, SARS-CoV-2-infected K18 mice at various ages are a newly established mouse model and have been successfully used for preclinical assessment of therapeutic and vaccine efficacy 49-51. In agreement with previous reports, we demonstrate that SARS-CoV-2 infected K18 mice at two to three months of age exhibit many features of severe COVID-19 seen in patients associated with severe lung histological changes including severe interstitial pneumonia with progressive perivascular inflammation, edema, hemorrhage, and platelet activation and aggregation. Mice exhibited at least one or more of these pulmonary histopathologic alterations. Thus, we recommend incorporating the entirety of histopathologic changes for scoring or quantitative analysis when applying this severe model for evaluating therapeutic and vaccine efficacy at a histological level. Additionally, in a previous study, flow analysis of peripheral blood mononuclear cells (PBMCs) showed marked decrease in the lymphocytes, including CD4+, CD8+ T and B cells in this severe model 22. In the current study, by CBC count, we also observed lymphocytopenia in this model and further document that the infected mice develop eosinopenia as well. Both lymphocytopenia and eosinopenia are common findings in COVID-19 patients upon admission to the intensive care unit (ICU) 52, 53. Further, infection in this model also triggered mild damage to extrapulmonary organs including cerebral infarction, myocarditis, and hepatitis, which are consistent with previous findings 21, 22. Here, we also provide evidence that SARS-CoV-2 infection of the lung, not the brain, is the major cause of death in these mice. Previous studies have also suggested that the infection of the lungs in this severe COVID-19 model is a main cause of death. In contrast, two recent studies indicated that neuroinvasion and brain inflammation following SARS-CoV-2 infection is associated with fatality in K18 mice 54, 55. Our results suggest that lung infection is the major cause of death in this severe COVID-19 model, reinforcing the notion that the K18 mice infected with a moderate dose of SARS-CoV-2 represent an optimal model of severe COVID-19 for preclinical and mechanistic studies. Although COVID-19 is considered primarily as a respiratory disease, increasing clinical evidence also shows a large number of COVID-19 patients experience neurological symptoms such as anosmia and ageusia 56. Autopsies of deceased COVID-19 patients suggest that SARS-CoV-2 produces pathological changes in CNS by either direct infection of the brain or indirect cytokine storm 57, 58. The severe lung disease accompanied with CNS infection in K18 mice is also consistent with clinical observations of severe cases of COVID-19, suggesting this model of a lower, sub-lethal dose of SARS-CoV-2 may be used for COVID-19-related CNS malfunction studies.

In this severe COVID-19 model, we identified previously recognized and un-recognized massive immune cell activation, infiltration of the lung and an array of cytokine and chemokine changes in the lung and circulation. We identified the overexpression of genes encoding for Ccl2, Ccl3, Ccl4, Ccl7, Ccl17, Cxcl9, Cxcl10, Cxcl13 in the lung and increased serum levels of Ccl2, Cxcl10 and IL-18 related to the infiltration and activation of immune cells specifically for inflammatory Mϕ/MC. These results are comparable with data reported in autopsy samples 59, PBMCs 60, 61, and bronchoalveolar lavage 62 from severe COVID-19 patients. In addition, we observed upregulation of type I and type II interferons, inflammatory response, and TNF-alpha signaling via NF-κB and the IL-6/STAT3 pathways. These findings are consistent with the observation obtained from bulk RNA-seq analysis and RT-PCR analysis of the lungs of infected K18 mice 22 and serum cytokine and chemokine arrays 20, 21, 63. Further, our scRNA-seq analysis showed a Cx3cr1/Ccr2+ Mϕ cluster, a neutrophil cluster, an IgDlo B cluster, and a trend towards increased type II pneumocytes in infected K18 mice. These data suggest that myeloid cell infiltration and activation contribute to the pathogenesis of severe COVID-19 similar to severe influenza pneumonia. scRNA-seq analysis also showed an increase in the fibroblast cluster, which has not been reported in COVID-19 patients or animal models before. It has recently been recognized that exuberant fibroblast activity compromises lung function via ADAMTS4 64. The increased activity of damage-responsive lung fibroblasts has been documented to drive lethal immunopathology during severe influenza virus infection 64. It remains to be determined if this pathway overlaps with COVID-19. Lack of the recognition of this finding in COVID-19 patients can be attributed to non-accessibility to their lungs in the early course of the disease and difficulty in conducting scRNA analysis of COVID-19 animal models under BSL3 restrictions. Whether the massive immune activation, including myeloid cell infiltration and activation, is linked to the increase in fibroblast cluster and how the increase in fibroblast cluster contributes to severe COVID-19 disease, remains unclear and requires further investigation.

We show that SARS-CoV-2 directly infects ECs in entire course of the COVID-19 disease in this severe murine and NHP COVID-19 models. Previous studies using these mice have mainly focused on SARS-CoV-2-induced lung histology changes and host immune responses rather than ECs 20-22, 65. Although early autopsy studies showed the presence of virus in ECs by EM 3, 8, concerns have been raised due to the difficulty in distinguishing virus particles from cellular structures in these studies 17, 66. In a COVID-19 patient autopsy cohort study, there was no evidence of viral RNA detectable in ECs 67. However, a recent study using in situ hybridization (ISH) to detect SARS-CoV-2 RNA found EC infection in 2 out of 32 cases 68. Other autopsy studies also showed the presence of virus in ECs of the stomach, ileum, liver 69 and heart 70. Since autopsies only reveal infection characteristics at late stages of COVID-19, during which infected cells are sporadically present 69, the infection of endothelium may be underestimated. A recent in vitro study suggests that SARS-CoV-2 can infect and induce a pro-inflammatory response in primary human ECs, even without productive viral infection 71. The strength of our study is the demonstration of SARS-CoV-2-infection of ECs using multiple complementary approaches, at both early and late stages of the disease and in a severe COVID-19 model. The functional consequence of the infection and the molecular mechanisms underlying the infection warrant further investigation. Further, we demonstrate that vast majority of endothelial cells infected by SARS-CoV-2 are alveolar capillary endothelial cells. Endothelial cells of K-18 mice express very low-level of hACE2 as compared with lung epithelial cells such as club and alveolar type (AT2) cells (), which are similar to human endothelial cells 1, 72. Whether hACE2 or mouse ACE2 in K-18 mice facilitate the alveolar capillary endothelial cell infection remains unknown and requires further investigation. As lung epithelial cells shed virus, the local SARS-CoV-2 levels increase dramatically. Alveolar epithelial cells and endothelial cells of the alveolar capillaries are in direct contact with one another which, may facilitate infection of endothelial cells by SARS-CoV-2 viral particles released by lung epithelial cells. There is also possible that similar to other intracellular pathogens such as listeria monocytogenes 73, 74, SARS-CoV-2 virus may be directly transmitted from epithelial cells to alveolar capillary endothelial cells through cell-cell contacts. This intercellular pathogen trafficking may enable infection of cells that express low level receptor such as hACE-2 for infectious virus.

Using the lungs of two severe COVID-19 models (K18 mice and AGMs 48) together with samples from autopsied patient, we demonstrate extensive pulmonary EC activation and dysfunction and increased vascular permeability, which explains the severe lung histopathology seen in these models including edema, hemorrhage and inflammatory and platelet activation. We used well-established cellular markers including increased ICAM and VCAM and decreased EV-cadherin, as markers of EC dysfunction. These biomarker changes in COVID-19 mouse models were further observed by heat map analysis of our bulk RNA-seq. Our findings showing lung EC dysfunction and damage are consistent with another study of SARS-CoV-2-infected rhesus macaques, which develop a mild COVID-19 disease 18. Moreover, we identify multiple pathways that may be associated with infection-mediated EC dysfunction and activation. First, we found that complement pathway activation may contribute to the EC dysfunction and activation in severe COVID-19. We observed increased complement C3 and MAC deposition in the EC layer of pulmonary vasculature. We also found increased transcript levels of C1qa, C1qb and Cfb in immune cells (mainly in Mϕ/MC), which are critical complement pathway components for complement activations. These results directly support the notion of massive acute and ongoing complement activation in the infected lung, which may contribute to EC dysfunction and activation. This prediction is further supported by previous studies by our group and others demonstrating that by products of complement activation contribute to pathogenesis of a variety of vascular diseases including atherosclerosis 75-77, ischemia and reperfusion injury 78, hyperacute graft rejection, vasculitis, hemolytic anemia 79-81, and the vascular complications of human diabetes 82. Second, we found evidence of extensive platelet and coagulation pathway activation in the lung. Previously, we reported the critical roles of platelet activation and coagulation pathway activation in fatal pulmonary arterial hypertension in a murine hemolysis model 27, 29. Together, these findings indicate that platelet activation and coagulation pathways could be direct consequences of EC infection and dysfunction and could also further exacerbate EC injury. Finally, we show the upregulation of the K-ras pathway in the fibroblast and EC clusters, a phenomenon previously unrecognized in COVID-19 patients and animal models. The effectors of K-ras signaling activated by SARS-CoV-2 infection all appear to mediate EC activation and permeability which play a critical role in cancer metastasis rather than cell proliferation, the other well-established effector of Kras signaling 83-85. The endothelial dysregulation of serpina3n has been found in atherosclerosis and aneurysm disease 86 and serpina3n attenuates rupture in a murine model of aortic aneurysm 87. Serpina3f, Serpina3g, and serpina3n are up regulated in blood brain barrier (BBB) in the mouse multiple sclerosis model 88. These findings indicate that SARS-CoV-2 infection directly or indirectly activates molecular pathways that may lead to increased vascular permeability and loss of EC barrier. Such a breach of EC permeability would have dramatic consequences for COVID-19 patients as the absence of the EC barrier virus may allow the virus to infect smooth muscle cells and cardiomyocytes which express ACE-2. The causative effect of the immune cell activation and the critical role of these cytokine and specific pathways on SARS-CoV-2 infection-mediated EC infection and dysfunction require further investigation.

In summary, we document SARS-CoV-2 infection-mediated massive immune activation, EC infection and dysfunction in a severe COVID-19 mouse. Our study sheds light on host viral interaction and highlights the pathogenic role of EC infection and dysfunction in the pathogenesis of severe COVID-19.

Methods

Mice and ethics statement

All animal experiments were approved by the Institutional Animal Care and Use Committees at Tulane University. K18(034860) and C57BL/6J (000664) mice were purchased from the Jackson Laboratory and housed in the animal facility of Tulane University.

SARS-CoV-2 strain and infection

We used SARS-CoV-2, Isolate USA-WA1/2020, NR-52281 deposited by the Centers for Disease Control and Prevention and obtained through BEI Resources, NIAID, NIH. We passaged the virus in VeroE6 cells in DMEM media with 2% FBS and sequenced the virus for verification as described previously 34. The harvested stock was determined to be 1.00E+6 TCID50/mL. Mice were intranasally infected by SARS-CoV-2(2 X 105 TCID or 1 X 105 TCID50) in ABSL3.

Histological analysis and quantification of histopathologic lesions. Sections of lung were processed routinely, stained with hematoxylin and eosin (H&E), and digitally scanned with a Zeiss Axio Scan.Z1. Inflammation and edema were quantified as previously described 34 and visually confirmed by a veterinary pathologist. Briefly, the stained lung sections were digitally scanned by a Zeiss Axio Scan.Z1 creating whole-slide images, and analyzed by a board certified pathologist with computer software (HALO v3.1, Indica Labs) using two algorithms (Multiplex IHC v2.3.4 and Spatial Analysis). Annotation regions were drawn around small arterioles, then all nucleated cells were counted using Multiplex IHC v2.3.4. Spatial analysis was used to quantify the number of nucleated cells within a 100 μm radius of the annotated vessels. Perivascular inflammation was reported as the density of nucleated cells within 100 μm radius of the tunica adventitia of small arterioles (nucleated cells/mm2).

RNA isolation

Tissues were collected in 1mL Trizol reagent (15596026; Invitrogen), and extracted with RNeasy Mini Kit (Cat. No.74104; QIAGEN) following the manufacturer's protocol. The concentration of RNA was determined by NanoDrop 2000.

Subgenomic N viral copy number detection

100ng total RNAs were mixed in Taqpath 1-Step Multiplex Master Mix (Cat. No.A15299; Thermo Fisher) and FAM-labeled primers (sgm-N-FOR: 5'-CGATCTCTTGTAGATCTGTTCTC-3', sgm-N-Probe: 5'-FAM-TAACCAGAATGGAGAACGCAGTGGG-TAMRA-3', sgm-N-reverse: 5'-GGTGAACCAAGACGCAGTAT-3'), following the manufacturer's instructions. Subgenomic N viral copy number was calculated by standard Cq values. The assay was performed under ABI QuantStudio 6 system.

Tissue mRNA array

RT² Profiler™ PCR Array Mouse Cytokines & Chemokines (PAMM-150Z, Qiagen, Hilden, Germany) was used to determine the expression of 84 cytokines and chemokines as previously decribed 89. The assay was performed under ABI QuantStudio 6 system.

Bulk RNA sequencing

Lung tissue RNA was used to perform RNA sequencing. RNA concentration was determined with a Qubit 3.0 Fluorometer (ThermoFisher Scientific). RNA integrity number and fragment sizes (DV200 metrics) were obtained with either the Agilent 2100 Bioanalyzer or the Agilent 4150TapeStation. Illumina TruSeq Stranded mRNA sample prep kit was used for library preparation. The cDNA libraries were pooled at a final concentration of 1.8 pM for cluster generation and sequenced using NextSeq 500/550 High Output Kit (150 Cycles) on an Illumina NextSeq 550 using a minimum of 4 ng of RNA according to SMART-Seq Stranded Kit User Manual (Takara Bio USA, Inc). Raw reads were processed and mapped, then gene expression and nucleotide variation were evaluated as previously described 90, 91. Raw read counts were normalized across all samples and then used for differential expression analysis using DESeq, EdgeR and Cuffdiff (Slug Genomics, UV Santa Cruz). The EdgeR output was used to generate volcano plots in R. Data was deposited in the Sequence Read Archive BioProject number: GSE175996.

Over-enrichment analysis. Enrichment analysis was performed with the R package 'ClusterProfiler' at default parameters. Enrichment scores were calculated against the following four pathway databases containing a priori-defined gene sets: Gene Ontology (GO) database, Hallmark (H) and Curated (C2) gene sets of the Molecular Signature database (MsigDB), and WikiPathways. Gene sets significantly enriched in the datasets (p < 0.05) were subsequently curated for those relevant to EC biological function. Enrichment plot was generated with the R software package ggplot2 and heatmaps with the ComplexHeatmap package.

Single cell RNA-seq analysis. Lungs were isolated and minced with forceps and scissors. They were further enzymatically digested with 2 mg/mL collagenase (Sigma) and 80 U/mL DNase I (Sigma) in 2 mL serum-free medium for 60 min at 37°C. Digested tissues were passed through a sterile 70 mm filter (Fisher Scientific) and 40 mm filter (Miltenyi Biotec) to generate a single cell suspension. After washing, cell suspensions were incubated with 2 ml of ACK Lysing Buffer (Cat #A1049201, Gibco) for 2 min and washed twice. 1x106 cells per condition were collected as whole lung single cell population. Cell numbers and viability were validated by Countess II (Thermo Fisher) prior to preparation of single cell RNAseq library.

10x single-cell RNA-seq assay. 5000 live cells per sample were targeted by using 10x single-cell RNA-seq technology provided by 10x Genomics (10X Genomics Inc, CA). Briefly, viable single cell suspensions were partitioned into nanoliter-scale Gel Beads-In-EMulsion (GEMs). Full-length barcoded cDNAs were then generated and amplified by PCR to obtain sufficient mass for library construction. Following enzymatic fragmentation, end-repair, A-tailing, and adaptor ligation, single cell 3' libraries comprising standard Illumina P5 and P7 paired-end constructs were generated. Library quality controls were performed using Agilent High Sensitivity DNA kit with Agilent 2100 Bioanalyzer and quantified by Qubit 2.0 fluorometer. Pooled libraries at a final concentration of 650 pM were sequenced with paired end dual index configuration by Illumina NextSeq 2000. Cell Ranger version 4.0.0 (10x Genomics) was used to process raw sequencing data and Loupe Cell Browser (10x Genomics) to obtain differentially expressed genes between specified cell clusters. In addition, Seurat suite version 2.2.1 90, 91 was used for further quality control and downstream analysis. Filtering was performed to remove multiplets and broken cells. Variable genes were determined by iterative selection based on the dispersion vs. average expression of the gene.

Immunohistochemistry assay

5 um sections of paraffin-embedded lung were baked overnight at 60 ºC, then dewaxed and rehydrated with xylene, graded ethanol and dd water. Slides were retrieved using a microwave oven heating method in Tris-EDTA buffer (ab93684, Abcam) for 20 min and cooled to room temperature. Sections were incubated with endogenous blocking solution (SP-6000-100, vector lab) for 10 min and 2.5% normal horse serum for 20 min at room temperature. Primary antibodies were incubated overnight at 4 ºC, then the slides were incubated with HRP horse anti-Rabbit IgG reagent for 30 min at room temperature. Immunoreactivity was detected using DAB system. Slides were digitally scanned by Zeiss Axio Scan. Z1.

RNA-scope

Z-fix-embedded lung sections were stained according to RNAscope® 2.5 HD Duplex detection kit guide. Z-fix-embedded lung sections were deparaffinized in fresh xylene and fresh 100% ethanol and then air dried. Target retrieval was performed (RNAscope® Target Retrieval Reagents, ACD Cat# 322000) after hydrogen peroxide treatment and was followed by protease treatment (RNAscope® H2O2 & Protease Plus, ACD Cat# 322330). Probes for RNAs of platelet endothelial cell adhesion molecule (RNAscope® Probe Mm-pcam1, ACD Cat# 316721) and spike protein (RNAscope® Probe V-nCoV2019-S-C2, ACD Cat# 848561-C2) were mixed and incubated on slides for in situ hybridization. Finally, slides were stained with 50% hematoxylin (Hematoxiylin Solution, Gill No.1, Sigma-Aldrich, Cat# GHS132) and mounted with VectaMount Permanent Mounting Media (ACD, Cat# 321584). Slides were digitally scanned by Zeiss Axio Scan.Z1.

Immunostaining

Immunostaining was performed as described previously 34. Zinc formalin fixed (Z-fixed), paraffin embedded lung was sectioned, deparaffinized, and subjected to heat induced epitope retrieval using both high-pH (Vector Labs H-3301), and low-pH (Vector Labs H-3300) solutions. Sections were blocked with 10% BSA or 1% donkey serum for 40 min, incubated with the primary antibodies for overnight at 4 0C and secondary antibodies for 40 min at room temperature (). Slides were digitally scanned by Zeiss Axio scan. Z1.

Scanning EM

Lung tissues were inflated and fixed with 2.5% glutaraldehyde buffer, washed 3 times with fresh cacodylate buffer (10 min each time), and postfixed by 1% OsO4 for 1 h at room temperature. Specimens were dehydrated in graded ethanol (30, 50, 70, 90% and 100%). Images were taken with Hitachi 4800 SEM.

Transmission EM

1mm3 lung cubes were fixed with 2.5% glutaraldehyde buffer followed by 1% OsO4. After dehydration, 90 nm thin sections were stained with UranyLess and lead citrate. Images were taken with JEM-1400.

Statistics

Data are shown as mean ± SEM. To compare values obtained from multiple groups over time, two-way analysis of variance (ANOVA) was used, followed by Bonferroni post hoc test. To compare values obtained from two groups, the unpaired Student's t-test was performed. Statistical significance was taken at the P < 0.05 level.

Supplementary figures and tables.

Click here for additional data file.