Rosa et al. find that genes associated with Interferon (IFN) signaling, neutrophil degranulation and innate immune pathways are induced in the lungs of macaques in response to SARS-CoV-2 infection, whereas genes associated with collagen formation and regulation are reduced. Authors also note an overlap of genes associated with SARS-CoV-2 infection and tuberculosis, altogether providing new insights into the immunopathogenesis of COVID-19 disease.
SARS-CoV-2 virus has infected more than 92 million people worldwide resulting in the Coronavirus disease 2019 (COVID-19). Using a rhesus macaque model of SARS-CoV-2 infection, we have characterized the transcriptional signatures induced in the lungs of juvenile and old macaques following infection. Genes associated with Interferon (IFN) signaling, neutrophil degranulation and innate immune pathways are significantly induced in macaque infected lungs, while pathways associated with collagen formation are downregulated, as also seen in lungs of macaques with tuberculosis. In COVID-19, increasing age is a significant risk factor for poor prognosis and increased mortality. Type I IFN and Notch signaling pathways are significantly upregulated in lungs of juvenile infected macaques when compared with old infected macaques. These results are corroborated with increased peripheral neutrophil counts and neutrophil lymphocyte ratio in older individuals with COVID-19 disease. Together, our transcriptomic studies have delineated disease pathways that improve our understanding of the immunopathogenesis of COVID-19.
IntroductionCoronavirus disease 2019 (COVID-19), caused by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), emerged as a pandemic disease during the end of 2019 and the beginning of 2020. Infected individuals develop a cytokine storm1, which can initiate viral sepsis and inflammation-induced lung injury leading to complications including pneumonitis, acute respiratory distress syndrome (ARDS), respiratory failure, shock, organ failure, and potentially death1,2.By combining established principles of anti-viral immunity with analysis of immune responses in COVID-19 patients, a picture of the host defense response against SARS-CoV-2 is beginning to emerge3,4. Upon infection of the mucosal epithelium, SARS-CoV-2 is detected by intracellular pattern recognition receptors (PRRs) that bind viral RNA and DNA. PRR signaling triggers activation of transcription factors and induces Interferon (IFN) signaling, which in turn activates resident macrophages. Infected macrophages induce cytokine secretion that consequently triggers the recruitment of myeloid cells, likely resulting in a feed-back loop that aggravates immunopathogenesis and promotes disease progression.Analyses of the transcriptomic response of host cells upon virus infection have the potential to identify the host immune response dynamics and gene activated regulatory networks5,6. Recent studies have reported transcriptional changes in cells in the broncho-alveolar lavage (BAL) and peripheral blood mononuclear cells (PBMCs) of COVID-19 patients7–11. Single-cell RNA-seq has recently identified initial cellular targets of SARS-CoV-2 infection in model organisms12 and patients13 characterizing peripheral and local immune responses in severe COVID-1914, thereby associating the severe disease with a cytokine storm and increased neutrophil accumulation. However, most studies have been performed in a limited number of moderate or severe COVID-19 patients within limited age ranges14. To get more in-depth understanding of the transcriptional changes during COVID-19, we have developed a SARS-CoV-2 macaque model, where both juvenile and old macaques were infected and exhibited clinical symptoms that reflect human COVID-19 disease that is self-limiting. In the current study, we have characterized the transcriptional signatures induced in the lungs of juvenile and old rhesus macaques following SARS-CoV-2 infection. We show that genes associated with IFN signaling, neutrophil degranulation, and innate immune pathways are significantly induced in the lungs in response to SARS-CoV-2 infection. This is associated with a downregulation of genes associated with collagen formation and regulation of collagen pathways. In COVID-19, increasing age is a significant risk factor for poor prognosis of infection15. We demonstrate that specific immune pathways, namely Type I IFN and Notch signaling, are significantly upregulated in juvenile macaques when compared with old macaques infected with SARS-CoV-2. These results are corroborated with increased peripheral neutrophil counts and neutrophil-lymphocyte ratio in older individuals with COVID-19 disease. In contrast, the VEGF pathway is downregulated in old infected macaques. Incidently, levels of VEGF protein are increased in plasma of older COVID-19 patients, emphasizing the importance of studying both local and peripheral responses. Finally, we report that neutrophil degranulation, innate immune system, and IFN gamma (IFN-γ) signaling pathways are upregulated in both tuberculosis (TB) and COVID-19, two pulmonary infectious diseases where neutrophils accumulation is associated with increased severity. Together, our study has delineated disease pathways that can serve as a valuable tool in understanding the immunopathogenesis of SARS-CoV-2 infection and progressive COVID-19.ResultsGenes upregulated in COVID-19 infected macaques represent pathways characteristic of neutrophil degranulation and IFN signalingWe recently assessed the ability of SARS-CoV-2 to infect rhesus macaques during a longitudinal two-week infection study16. Indian-origin, SPF-rhesus macaques (Macaca mulatta) were infected by multiple routes (ocular, intratracheal, and intranasal) with the sixth-passage virus at a target dose of 1.05 × 106 PFU/per animal and studied for two weeks. The macaques were grouped as naïve (uninfected), and infected (juvenile or old) macaques All infected animals developed clinical signs of viral infection16. Both juvenile and old macaques exhibited comparable clinical disease, and equivalent longitudinal viral loads in the BAL, nasopharyngeal and buccopharyngeal swabs, as well as lungs at endpoint. This was followed by comparable viral clearance. In the current study in order to fully understand the immune pathways regulated upon SARS-CoV-2 infection, RNA was extracted and RNA sequencing was carried out from a lung biopsy from juvenile macaques (n = 3, 1 male and 2 females) and old macaques infected with SARS-CoV-2 (n = 5, 1 male and 4 females) and naive uninfected macaques (n = 4, 2 males and 2 females) (Supplementary Data 1). An average of 68.6 million reads was generated, with an average of 20.3 million fragments (read pairs or orphaned reads) mapping to macaque coding sequences, following analytical processing and mapping (Supplementary Data 1). Principal components analysis (PCA) based on whole-transcriptome gene expression levels17 showed that despite within-group variability for the COVID-19 infected samples, the naive samples grouped separately, suggesting substantial overall transcriptomic differences resulting from the infection (Fig. 1a). Differential gene expression analysis (DESeq217) with the juvenile and old COVID-19 samples grouped together identified 1,026 genes significantly (P ≤ 0.01) upregulated in response to infection, while 1109 genes were significantly downregulated (Fig. 1b). Expression, annotation, and differential expression data for all genes is available in Supplementary Data 2. Complete lists of differentially expressed genes for each comparison of interest (described below) ranked by P value, with Z-scores for expression visualization are available in Supplementary Data 3, and significant pathway enrichment (Reactome18, KEGG19, and Gene Ontology20) for all comparisons is shown in Supplementary Data 4.Fig. 1Gene expression profiles for lung tissue samples from macaques infected with SARS-CoV-2 are distinct from naive uninfected samples.a PCA plot showing the clustering of samples based on overall transcriptomic profiles. b Differential gene expression plot showing the relative normalized gene expression levels (FPKM) for each gene, with genes significantly differentially regulated by COVID-19 are indicated.Evaluation of the top 30 most significantly upregulated genes in the lungs of SARS-CoV-2-infected macaques revealed significantly higher expression of CTSG (Cathepsin G), ATP6AP2 (ATPase H+ transporting accessory protein 2), IFNγR1 (Interferon Gamma Receptor), CD36, and CD58, in comparison to expression in uninfected macaque lungs (Fig. 2a). Cathepsin G is a serine protease prominently found in neutrophilic granules. IFNγR1 associates with IFNγR2 to form a receptor for the cytokine interferon gamma (IFNγ)21–24, and required for activation of antiviral responses, such as IRF3 (IFN regulatory factor-3), nuclear factor ΚB (NF-ΚB) and JAK (Janus kinase)/STAT (signal transducer and activator of transcription) signaling pathways25. Reactome pathway analysis on upregulated and downregulated genes in the lungs of SARS-CoV-2-infected rhesus macaques showed that genes significantly upregulated by infection, included pathway enrichment for genes involved in “Neutrophil degranulation”, “Innate Immune system”, “Immune system” and “IFN signaling” (Supplementary Table 1 and Supplementary Data 4A). The upregulation of CD36 during COVID-19 in lungs is in conformity with these enriched pathways, since CD36, a scavenger receptor expressed in multiple cell types, mediates lipid uptake, immunological recognition, inflammation, molecular adhesion, and apoptosis26, and is a Matrix Metalloproteinase-9 substrate that induces neutrophil apoptosis. CD58 molecule (lymphocyte function-associated antigen-3) is expressed on human hematopoietic and non-hematopoietic cells, including dendritic cells, macrophages, and endothelial cells27–30, and interacts with its receptor CD2 molecule31,32 on CD8+ cytotoxic T lymphocytes and NK cells to mediate cytotoxic reactions33–35. The complete ranked list of the 1026 genes upregulated during COVID-19 is shown in Supplementary Data 3A.Fig. 2Genes differentially regulated in SARS-CoV-2-infected macaques represent pathways characteristic of neutrophil degranulation, IFN signaling, collagen degradation and TFG-β signaling.The top 20 most significantly a upregulated genes and b downregulated genes in SARS-CoV-2-infected macaque lungs. Expression values are visualized by Z scores of normalized expression data (FPKM) per sample, and Log2 Fold Change and -Log P values are from the DESeq2 output. Genes are sorted by P value. c Gene expression plot showing the relative normalized gene expression levels (FPKM) for each gene, with “neutrophil degranulation” (R-HSA-6798695) pathway genes highlighted, including those that were significantly upregulated during COVID-19 (blue) and those that were not (orange). d Multilabel confocal immunofluorescence microscopy of FFPE lung sections from SARS-CoV-2-infected rhesus macaques with SARS-CoV-2 Spike specific antibody (green), neutrophil marker CD66abce (red), and DAPI (blue) at 63X magnification.ATP6AP2 was the most significant of the 65 genes upregulated within the enriched “neutrophil degranulation” (R-HSA-6798695) pathway (Supplementary Data 3B). It interacts with renin or prorenin to cause activation of intracellular signaling pathways, resulting in secretion of inflammatory and fibrotic factors36. CEACAM8 (Carcinoembryonic Antigen-Related Cell Adhesion Molecule 8) is the gene that encodes for CD66b, a well-characterized marker of degranulation37. Significant upregulation of neutrophil degranulation genes was observed in SARS-CoV-2-infected macaques (Fig. 2c). Indeed, CD66b+ neutrophils accumulate in the lungs of macaques infected with SARS-CoV-2 (Fig. 2d). We have also previously demonstrated that neutrophils are heavily recruited early to the alveolar space following SARS-CoV-2 infection of macaques16. Additional genes strongly upregulated during COVID-19 in the neutrophil degranulation pathway are IDH-1 (Isocitrate Dehydrogenase (NADP(+) 1) which regulates neutrophil chemotaxis, and FPR2 (Formyl Peptide Receptor 2), a G-coupled surface receptor which has a deleterious role to play in viral infection including influenza38. LTA4H (Leukotriene A4 hydrolase) is an enzyme that generates a neutrophil chemoattractant, leukotriene B4, a marker for ARDS39. Expression of 162 genes belonging to the “immune system” (R-HSA-168256) pathway was upregulated in SARS-CoV-2 infected macaques (Supplementary Data 3C). These included LAMP-2 (Lysosomal Associated Membrane Protein 2), and ATG7 (Autophagy Related 7), key genes involved in autophagy. LAMP-2 is known to influence phagosomal maturation in neutrophil40. The IFN response constitutes the major first line of defense against viruses. Consistent with this, we found upregulation of genes associated with the IFN signaling pathways, specifically Interferon Induced Protein with Tetratricopeptide Repeats 1 (IFIT3), IFN alpha receptor 1 (IFNAR1), IFN gamma receptor 1 (IFNGR1), and OAS 1 protein (2’-5’-Oligoadenylate Synthetase 1). Together, these results suggest that upregulation of neutrophil degranulation, Type I IFN signaling, and the innate immune system is a characteristic feature of host responses to SARS-CoV-2 infection.Genes downregulated following SARS-CoV-2 infection in macaques represent pathways characteristic of collagen degradation and TFG-β signalingUp to 40% of patients with COVID-19 develop ARDS, and 20% of ARDS cases are severe41. A well-documented sequela of ARDS is the development of fibrotic disease42,43. We found that the 1109 genes downregulated in SARS-CoV-2-infected macaques were significantly enriched for collagen degradation, regulation, and formation (Fig. 2b, Supplementary Table 2; Supplementary Data 3D and Supplementary Data 4B). For example, among the “collagen degradation” (R-HSA-1442490) enriched pathway genes (Supplementary Data 3E), COLA1 (collagen type I chain), other members of the collagen gene family (COL4A2 COL16A1 COL4A4 COL6A2 COL6A1 COL5A1 COL9A1 COL13A1 COL12A1 COL1A2) were all significantly downregulated in COVID-19 diseased lungs when compared with expression in lungs of uninfected controls. Although matrix metallopeptidase (MMPs) are proteolytic enzymes responsible for extracellular matrix protein degradation, in our study matrix metalloproteases such as MMP23B, MMP15 and MMP14 were downregulated in COVID-19 diseased lungs compared with expression in lungs of uninfected controls. In addition, Reactome pathway enrichment prominently featured pathways downregulated in COVID-19 disease in macaques comprised of “collagen degradation”, “collagen chain trimerization”, “degradation of extracellular matrix” and “collagen formation” (Supplementary Table 2). Increased collagen degradation is essential for the prevention of fibrosis, a sequelae of COVID-19 and ARDS. Therefore, regulation of collagen degradation and extracellular matrix modeling suggests that this may be a feature of SARS-CoV-2 infection of rhesus macaques being a self-limiting model with early and robust anamnestic responses. TGFβ (Transforming Growth Factor Beta 1) is involved in normal tissue repair following lung injury, and in mediating fibrotic tissue remodeling by increasing the production and decreasing the degradation of connective tissue44. Our results indicate downregulation of genes associated with TGFβ signaling (Supplementary Table 2), including the genes PARD3 (par-3 family cell polarity regulator) and PARD6A (par-6 family cell polarity regulator alpha), which are involved in regulating epithelial cell apico-basolateral polarization, SMURF (SMAD specific E3 ubiquitin-protein ligase 1), a negative regulator of TGFβ pathway, and FURIN, which is a TGFβ converting enzyme (Supplementary Data 3F). While the interaction of the genes within these pathways is complex, our results project a broad downregulation of mechanisms that contribute to lung repair and remodeling in animals with anamnestic control of SARS-CoV-2 infection.Type I IFN signaling and Notch signaling pathways are upregulated in young macaques but not old macaques with COVID-19 diseaseAge is a significant risk factor for increased morbidity and mortality in COVID-19 disease15. In order to identify the differential immune responses associated with SARS-CoV-2 infection in old macaques, we carried out differential expression analysis between juvenile (n = 3) vs naive (n = 4), and old (n = 5) vs naive (n = 4) groups. In order for a gene to be considered to be differentially expressed only in the juvenile macaques, we required a stringent P value for significance ≤0.01 in the juvenile COVID-19 vs naive, and a P value for significance ≥0.1 in the old COVID-19 vs naive comparison. This approach identified 86 genes significantly upregulated (Fig. 3a and Supplementary Data 3G) and 96 genes significantly downregulated (Fig. 3b and Supplementary Data 3H) with COVID-19 disease only in juveniles. Note that no genes were significantly upregulated in juveniles and significantly downregulated in old, and vice-versa. Of these genes, the top 30 most significantly differential between juvenile and old are shown for upregulated genes in Fig. 4a and for downregulated genes in Fig. 4b. No pathways were found to be significantly enriched among the 96 genes significantly downregulated only in juveniles, but the Reactome and KEGG pathways significantly enriched among the 86 genes upregulated only in juveniles are shown in Supplementary Table 3. Complete gene lists per pathway, and all significant pathways enrichment results including for Gene Ontology (GO) are available in Supplementary Data 4C.Fig. 386 genes were significantly upregulated and 96 genes were significantly downregulated with COVID-19 only in juvenile macaques.Heatmaps visualizing the significance values of COVID-19 a upregulated and b downregulated genes, in juvenile and old macaques. Green shaded areas contain genes significant only in juveniles, and red shaded areas contain genes significant only in old macaques.Fig. 4Genes related to Type I interferon signaling are upregulated in juvenile macaques compared to old macaques during COVID-19-infection.The top 20 most significantly a upregulated genes and b downregulated genes in SARS-CoV-2 infected juvenile macaque lungs but not in old macaques. Expression values are visualized by Z scores of normalized expression data (FPKM) per sample, and Log2 Fold Change and -Log P values are from the DESeq2 output. Genes are sorted by P value. c The relative gene expression of ACE2 and ADAM17 among naive, juvenile and old SARS-CoV-2-infected macaques.The genes with significantly upregulated expression in SARS-CoV-2 infected juvenile but not old macaques included MX1 (MX Dynamin Like GTPase 1), MX2 (MX Dynamin Like GTPase 2), and USP18 (Ubiquitin Specific Peptidase 18) (Fig. 5a, b). This is consistent with and highlights the role of the Reactome pathway “interferon alpha/beta signaling” being enriched in juvenile macaques during SARS-CoV-2 infection (Supplementary Table 3 and Supplementary Data 4C). Other genes in this pathway that exhibited increased expression included IFIT1 and IFIT2 (Fig. 5a, b). In addition, by KEGG analysis, the Notch signaling pathway was observed to be significantly upregulated in juvenile infected macaques when compared with old infected macaques. ADAM17 (ADAM Metallopeptidase Domain 17), a key component of the Notch signaling pathways is known to be involved in the shedding of the surface protein ACE2 (Angiotensin-converting enzyme 2)45. Therefore, it is interesting that a linear correlation in the expression of ACE2 and ADAM17 exists in infected macaques (Fig. 4c). Note that we also see a significant upregulation of ACE2 across SARS-CoV-2-infected samples and a substantially larger upregulation among the juvenile samples (7.1-fold, P = 3.4 × 10-4). In addition, the induction of DLL4, a Notch ligand, was increased in the infected juvenile macaques. Finally, the differential induction of DTX3L (Deltex E3 Ubiquitin Ligase 3L) in juvenile infected macaques compared to old infected macaques is important because Deltex stabilizes the receptor in the endocytic compartment allowing signal transduction to proceed in Notch signaling46. Of the Hepatitis-induced pathway genes that are upregulated in juvenile COVID-19 diseased lungs, CXCL-10 (C-X-C Motif Chemokine Ligand 10) is a chemokine associated with severe disease in COVID-19 in humans47, but can also be involved in the recruitment of CXCR3 (C-X-C Motif Chemokine Receptor 3) expressing immune cells. 14-3-3 (otherwise called YWHAG) interacts with MDA5 (melanoma differentiation-associated protein 5), which belongs to the RIG-I-like receptor family and drive anti-viral immunity. Together, these results suggest that specific pathways including Type I IFN and Notch signaling are highly induced in juvenile macaques during SARS-CoV-2 infection, when compared to similarly infected old macaques.Fig. 5Interferon alpha signaling genes are significantly upregulated in juvenile SARS-CoV-2-infected macaques but not old SARS-CoV-2-infected macaques.a Differential expression significance values for genes from the “Interferon-alpha signaling” (R-HSA-909733) pathway, in both the juvenile vs naive and old vs naive comparisons. Included are significance values for all genes (gray), pathway genes significantly upregulated in juvenile but not old (orange) and other pathway genes (green). P value cutoffs (0.01 and 0.1) are indicated with dashed lines. b The relative expression levels (FPKM) for the five “interferon-alpha signaling” genes belonging to this gene set are shown. P values represent FDR-corrected significance values from DESeq2.Genes related to VEGF signaling are downregulated in old macaques but not juvenile macaques during COVID-19-diseaseUsing the same approach as for the juvenile macaque-specific differentially regulated genes, we identified 97 genes significantly upregulated (Fig. 3a and Supplementary Data 3I) and 160 genes significantly downregulated (Fig. 3b and Supplementary Data 3J) with COVID-19 disease only in infected old macaques, and not infected juveniles. Pathway enrichment analysis only identified significant functional enrichment among the downregulated gene set (Supplementary Table 4 and Supplementary Data 4D). Our results show that in the lungs of old macaques, the only Reactome pathways enriched among genes downregulated during COVID-19 included genes involved in the “VEGF-VEGFR2 Pathway” and “Signaling by VEGF” (Figs. 6a, b and 7a, b). Vascular endothelial growth factor (VEGF) is a signaling protein that promotes angiogenesis, and is a key factor that promotes ARDS. Previous research depicts that ACE2 antagonizes and down-regulates VEGFA48, improving lung function following acute lung injury49. Here, we observe both a significant increase in ACE2 in response to COVID-19 and a significant decrease in VEGF pathways in old macaques, which may be due to this antagonistic relationship. VEGFA, p21-activated kinase (PAK2), cytoplasmic tyrosine kinase (SRC), RhoA/ROCK signaling [ROCK1(Rho Associated Coiled-Coil Containing Protein Kinase 1) and WASF2 (WASP Family Member 2) are all essential for multiple aspects of VEGF-mediated angiogenesis and are all significantly downregulated in old macaques with COVID-19 (Fig. 7a, b). Overall, despite juvenile and old macaques having a comparable clinical course with resolution, our data suggest that there are significant differences in signaling pathways, especially those related to VEGF signaling that may ultimately result in differences in long-term outcomes. Thus, our results suggest that downregulation of VEGF pathways is associated with increasing age, in a macaque model of the self-limiting disease.Fig. 6Genes related to VEGF signaling are downregulated in old macaques compared to juvenile macaques during COVID-19.The top 20 most significantly a upregulated genes and b downregulated genes in infected old macaque lungs but not in juvenile macaques. Expression values are visualized by Z scores of normalized expression data (FPKM) per sample, and Log2 fold change and -Log P values are from the DESeq2 output. Genes are sorted by P value.Fig. 7VEGF pathway genes are significantly downregulated in old SARS-CoV-2- infected macaques but not juvenile SARS-CoV-2-infected macaques.a Differential expression significance values for genes from the “Signaling by VEGF” (R-HSA-194138) pathway, in both the juvenile vs naive and old vs naive comparisons. Included are significance values for all genes (gray), pathway genes significantly downregulated in old but not juvenile (orange) and other pathway genes (green). P value cutoffs (0.01 and 0.1) are indicated with dashed lines. b The relative expression levels (FPKM) for the seven “Signaling by VEGF” genes belonging to this gene set are shown. P values represent FDR-corrected significance values from DESeq2.Aged COVID-19 patients exhibit increased plasma VEGF protein levels and high peripheral neutrophil to lymphocyte ratioTo further address if our findings were relevant in the human setting of SARS-CoV-2 infection, we stratified COVID-19 patients into an aged group (>60 years) and a group of COVID-19 patients <60 years (Supplementary Data 5). We found that with increasing age, there was an increased association of disease parameters and comorbidities (Supplementary Data 5). We measured the levels of human plasma proteins levels for IFN-α, IFN-β, and IFN-γ. While levels of plasma IFN-α, and IFN-β were below the levels of reliable detection, we found that the COVID-19 patients who were <60 years expressed significantly higher plasma IFN-γ levels when compared to levels in the plasma of healthy controls (Fig. 8a). Although plasma levels of IFN-γ protein were also increased in the aged COVID-19 patient group, levels were not significantly different from healthy controls (Fig. 8a). This was in contrast to plasma protein levels of VEGF, which were significantly higher in aged individuals with COVID-19 disease when compared with levels in individuals with COVID-19 disease who were <60 years old (Fig. 8b). The increased levels of VEGF in aged COVID-19 patients coincided with significantly increased peripheral neutrophil counts as well as increased peripheral neutrophil to lymphocyte ratios, when compared with both healthy controls and COVID-19 group <60 years old (Fig. 8c, d). These results show that plasma protein levels of VEGF and accumulation of peripheral neutrophils is increased in aged individuals with COVID-19 disease, when compared to younger individuals with COVID-19 disease.Fig. 8VEGF and peripheral neutrophil counts are higher in old COVID-19 patients.Peripheral blood samples were obtained from a cohort of patients with laboratory-confirmed SARS-CoV-2 infection at hospital admission. Levels of different immune markers were determined by Luminex assay in plasma samples from COVID-19 and healthy volunteer controls. COVID-19 patients were stratified by age as younger than or older than 60 years. a Levels of IFN-γ and b levels of VEGF proteins were measured in plasma of COVID-19 and healthy controls. c Peripheral neutrophil counts and d neutrophil to lymphocyte ratio (NLR) values were retrieved from the medical records of COVID-19 patients and compared between age groups. Each point in the figure corresponds to an independent biological sample for each group. Mean ± SD for each group is shown in the figure. Statistical significance was verified by one-way anova: P = *(<0.05), **(<0.01) with Tukey’s correction (a and b) or by Student’s t-test: P = *(<0.05), **(<0.01) (c and d).Neutrophil degranulation and IFN pathways overlap between COVID-19 and TB diseaseTB is a pulmonary granulomatous disease caused by infection with Mycobacterium tuberculosis. TB disease in humans and macaques is associated with neutrophil and IFN signature50. Thus, we next compared and contrasted the transcriptional profile of genes and pathways that are shared by the two diseases, and those that are unique to COVID-19. There was not a substantial overlap between differentially expressed genes in response to COVID-19 and TB. However, of the 97 genes that were commonly upregulated in TB and COVID-19 (Fig. 9a and Supplementary Data 3K), the Reactome pathway enrichment was well featured in “Neutrophil degranulation”, “Innate immune response”, and “Interferon-gamma signaling” (Fig. 9b and Supplementary Data 4E). Nearly as many genes (76) had opposite differential expression patterns (upregulated in COVID-19, downregulated in TB), as genes upregulated in both diseases (Fig. 10a and Supplementary Data 3L). These genes were associated with blood vessel morphogenesis and angiogenesis including leptin receptor (LEPR) and TGFβ2 (Fig. 10b and, Supplementary Data 4F). These results suggest that both TB and COVID-19 share features of neutrophil accumulation of IFN signaling, but that COVID-19 disease immunopathogenesis uniquely features vascularization of the lung.Fig. 9Genes higher in expression during both COVID-19 and TB share common pathways.a The top 20 (of 97) most significantly upregulated genes in COVID-19 infected and TB infected macaques. Expression values are visualized by Z scores of normalized expression data (FPKM) per sample, and Log2 fold change and -Log P values are from the DESeq2 output. Genes are sorted by P value. b Significant Reactome pathway enrichment among the 97 genes.Fig. 10Genes higher in expression during COVID-19 than TB are related to blood morphogenesis pathways.a The top 20 (of 76) most significantly upregulated genes in SARS-CoV-2-infected compared to Mtb-infected macaques. Expression values are visualized by Z scores of normalized expression data (FPKM) per sample, and Log2 fold change and -Log P values are from the DESeq2 output. Genes are sorted by P value. b Significant Gene Ontology pathway enrichment among the 76 genes.DiscussionLack of understanding of the complexity of COVID-19 immunopathogenesis hampers the identification of therapeutic strategies for COVID-19. While immune profiling in COVID-19 patients has shed light on related immune mechanisms of this disease, this approach has primarily involved peripheral samples obtained from moderate to severe COVID-19 patients, who are generally also older. To overcome these limitations and serve as a hypothesis-generating platform, we have developed a nonhuman primate model (rhesus macaque) of SARS-CoV-2 infection. This model reflects several features of the immunopathogenesis of human COVID-19, and allows us to interrogate the immune pathways that mediate disease versus protection, especially in the context of young versus older hosts. In this study, we show that upregulation of pathways characteristic of neutrophil degranulation and IFN signaling are characteristic of COVID-19 disease in infected hosts. Importantly, the induction of genes associated with Type I IFN signaling pathway and Notch signaling in young macaques infected with SARS-CoV-2 is a key determinant that distinguishes them from infected old macaques. Lungs of old macaques infected with COVID-19 however, uniquely feature downregulation of VEGF signaling pathways. Importantly, in PBMCs of humans infected with SARS-CoV-2 we found increased levels of VEGF and peripheral neutrophil counts in individuals >60 years when compared to younger individuals. These results together provide novel insights into the immunopathogenesis of COVID-19 disease, especially from the unique perceptive of age as a contributing factor.As we learn more about the pathophysiology of COVID-19, it is becoming clear that disease severity is associated with hyperinflammation which in turn induces lung and multiorgan injury and mortality via a cytokine storm1,2,51. While therapeutic options that focus on immunomodulatory agents such as corticosteroids are being considered and used, immunomodulators may also inhibit protective pathways. Therefore, a thorough understanding of the host inflammatory responses during SARS-CoV-2 infection is needed before precise immunomodulators can be specifically designed to limit inflammation without regulating protective mechanisms of action. The distinct role of myeloid cells in COVID-19 lung injury and immunopathogenesis is beginning to be described in humans, where recruited macrophages, inflammatory monocytes and neutrophils are associated with more severe COVID-19 disease and hyper-inflmmatory cytokine environment8,9,11. Furthermore, patients with COVID-19 exhibited hyper-inflammatory signatures across all types of cells among PBMCs, particularly upregulation of the TNF/IL-1β-driven inflammatory response and type I IFN responses10. Indeed, we have clearly shown that neutrophils are recruited to the macaque lung after SARS-CoV-2 infection. Neutrophils play a protective role contributing to early antiviral defense52, but also can be pathological due to processes associated with degranulation and lysis, thereby promoting lung inflammation. Consistent with this notion, an increased peripheral neutrophil-to-lymphocyte ratio is observed in severe COVID-19 cases, and in some studies is also associated with unfavorable prognosis53. These results in human studies are consistent with our macaque studies that describe neutrophil degranulation as one of the transcriptional pathways upregulated in the lungs of COVID-19 macaques when compared to uninfected controls. In this regard, the expression of Cathepsin G is noteworthy since it is a prominent serine protease that amplifies inflammation by stimulating the production of cytokines and chemokines that drive immune cell recruitment to the lung54, and activates metalloproteases to cleave extracellular matrix proteins, thereby promoting neutrophil migration55. Cathepsin G also induces potent chemotactic recruitment of monocytes, neutrophils and antigen-presenting cells in addition to promoting endothelial and epithelial permeability56. The latter function of Cathepsin G could enhance viral invasion to extra-alveolar sites while increased epithelial permeability might also explain the gastrointestinal route of transmission16. In addition, ATP6AP2, causes the secretion of inflammatory and fibrotic factors36, CD36, that induces neutrophil apoptosis, and CECAM8 whose cross-linking induces IL-8 production, all of which are highly expressed in COVID-19 diseased lungs. In patients with severe COVID-19, a higher frequency of CD66b+ neutrophils exist57. Furthermore, our studies shed light on the importance of the membrane glycoprotein, CD36 in the response to SARS-CoV-2 infection. CD36 is expressed on platelets, macrophages, and even epithelial cells. In addition to its well-characterized apoptotic function, CD36 is also a receptor for thrombospondin-1 and related proteins and can function as a negative regulator of angiogenesis58. This is particularly important given that angiogenesis is an important feature in patients with COVID-19 and associated ARDS59. CD36 also binds long-chain fatty acids and facilitates their transport into cells, leading to muscle utilization, coupled with fat storage. This contributes to the pathogenesis of metabolic disorders, such as diabetes and obesity and atherothrombotic disease59. A recent single-cell analysis revealed significantly higher CD36 expression in association with ACE2-expressing human lung epithelial cells60. Increased CD36 expression may therefore provide a protective role from an extreme lung injury during COVID-19, which is observed in the macaques. Our novel findings that CD36 (as well as other prominent signaling pathways) may be involved in the pathogenesis of COVID-19 has implications for host-directed therapy for SARS-CoV-2 infection. In contrast, neutrophils are recruited into the lung very early following macaque infection with SARS-CoV-216. In addition, in the lungs of deceased individuals with severe COVID-19 disease, neutrophil infiltration occurred in pulmonary capillaries and was accompanied by extravasation of neutrophils into the alveolar space, and neutrophilic mucositis61. In the case of COVID-19, neutrophils could also be a source of excess neutrophil extracellular traps62. These results together suggest a scenario in the lung where induction of the cytokine storm drives the recruitment of neutrophils, thereby contributing to inflammation. Thus, degranulation of neutrophils and formation of NETs may further promote cytokine responses and inflammation and disease immunopathogenesis.The IFN response constitutes the major first line of defense against viruses. Recognition of viral infections by innate immune sensors activates both the type I and type III IFN response. While some studies have shown that serum of COVID-19 patients contains increased expression of proinflammatory cytokines and chemokines, without detectable levels of type I and III IFNs63, other studies suggest that the IFN response may be delayed. Importantly, elevated IFNs correlate with more severe disease64,65. Thus, it is unclear if type I IFNs are protective or pathological in COVID-1966. It is possible that severe infection drives the higher expression of IFN pathway genes, but may not lead to viral containment, but instead drives pathological damage. On the other hand, increased induction of type I IFN pathways in SARS-CoV-2 infected macaques, as well as in juvenile macaques, could support a role for IFN signaling in protection rather than disease progression. Our studies provide data to support the recently proposed hypothesis that IFN induction may be compromised in older hosts66. When the early IFN response is not optimal to control viral infection, it is possible that delayed or inadequate IFN responses may lead to inflammation-mediated damage. Further testing the protective versus pathological roles of IFNs in the macaque model with the availability of IFNAR blocking reagents should further clarify the specific role of IFN pathways in COVID-19.ARDS in influenza, MERS, and SARS have been associated with fibrotic irreversible interstitial lung disease67,68. Pulmonary fibrosis is a recognized sequelae of ARDS42. Pulmonary fibrosis can develop either following chronic inflammation or as a consequence of genetically associated and age-related fibroproliferative process, as in idiopathic pulmonary fibrosis (IPF)69. Fibrosis is the hardening, and/or scarring of tissues due to excess deposition of extracellular matrix components including collagen. Since a significant proportion of COVID-19 patients develop severe ARDS, it is predicted that a similar outcome of fibrosis will be associated with COVID-19. Also, since the risk factors associated with COVID-19 including increasing age, male and associated co-morbidities coincide with risk factors, it is expected that COVID-19 patients will experience fibrotic lung disease. Despite these associations, there is no evidence currently that “scarring of the lung” experienced by COVID-19 patients is fibrotic or progressive and an outcome of COVID-19 disease post recovery. Therefore, our results provide unique insights into the role of fibrosis during SARS-CoV-2 infection. Most notably, we find downregulation of collagen degradation pathways, as well as pathways associated with collagen formation, collagen trimerization, and assembly. Furthermore, the role of TGF-β and ECM degradation is well documented in fibrosis. Indeed, the genes associated with these pathways are also significantly downregulated. These results for the first time provide important insights into the early pathological events occurring during COVID-19 in the lungs with relevance to underlying immune mechanisms associated with canonical fibrosis pathways. While long-term consequences of the pulmonary COVID-19 such as fibrosis remain to be determined, our results on downregulation of collagen degradation and TGF-β pathways may represent important early events on the lungs of SARS-CoV-2 infected individuals.Finally, we provide novel insights into the transcriptional regulation of immune pathways that are induced and regulated by age, an important risk factor for COVID-19 disease and outcome. This is a significant component of risk for disease and prognosis of COVID-19. We find higher induction of genes associated with Type I IFN signaling and Notch signaling in the young macaques. While we had sufficient macque samples to compare young and old macaques, we do not have sufficient representation of both sexes to study sex-associated modulation of disease transcriptional profiles. Upregulation of these significant Type I IFN signaling genes suggest that in a model of self-limited clinical disease in macaques, Type I IFN induction may be differentially regulated by age-associated factors. Age-specific regulation of this pathway has been demonstrated in the murine model of TB70. There is also a well-documented relationship between Notch signaling and viral infections. For example, Human Papilloma Virus and Simian Virus 40 can highjack the cellular machinery, including components of Notch signaling, and these events re-associated with cancer progression71. While innate and T cell responses are comparable between juvenile and old macaques following infection, SARS-CoV-2 specific antibody is generated at significantly higher levels in the plasma of juveniles, relative to old macaques16. Since Notch signaling regulates multiple stages of B-cell differentiation and shapes the antibody repertoire72, higher expression of many of the Notch pathway member genes in juvenile macaques may be responsible for the development of stronger antibody responses in these animals, impacting disease progression. Similarly, Type I IFN responses are critical for the downstream breadth of antibody production and recognition73–75. Thus, while T cell responses are comparable in juvenile and old macaques, differences in critical signaling pathways uncovered by our RNA-seq analysis potentially explain why juvenile macaques mount significantly stronger antibody responses, and consequently why younger subjects have reduced susceptibility to COVID-19. While this has not been recapitulated in the macaque model, older patients of COVID-19 are more susceptible to progression. This is consistent with increased disease progression when COVID-19 patients were stratified based on age. A previous study found that peripheral VEGF concentrations were significantly higher in COVID-19 patients than in healthy controls76. We also find this effect in our human samples (Fig. 8b) where people with COVID-19 that are older than 60 years of age have more VEGF protein in their peripheral blood. However, we also find significantly lower levels of VEGF pathway gene transcripts in the lungs of macaques with SARS-CoV-2 infection, especially older macaques (Fig. 6b and Fig. 7b). Our study further demonstrates that the changes in VEGF signaling may be associated with increasing age rather than just with disease severity. VEGF pathways promote angiogenesis and induce vascular leakiness and permeability. Our results therefore suggest that higher levels of VEGF in the periphery, while a biomarker for COVID-19, may be driven as a compensatory mechanism due to lower levels of VEGF signaling at the site of infection, i.e., the lung. These results further underscore the value of studying responses to SARS-CoV-2 infection in the lung compartment. By uncovering new aspects of the role of these signaling pathway in SARS-CoV-2 infection in the lung compared to the periphery using animal models and human samples, will shed further light on pathways that can be harnessed for therapeutics for COVID-19 disease.TB and COVID-19 both primarily affect lung function. TB was already one of the leading causes of death due to an infectious disease prior to the emergence of COVID-19. In the current scenario the clinical management of both TB and COVID together, particularly in the endemic regions is another rapidly emerging healthcare challenge needing immediate attention. Our results are the first to clearly demarcate the main differences in the manifestation of both the diseases in the alveolar niche. Neutrophil degranulation was one of the most significantly enriched pathways in both the disease conditions and therefore appears as a promising druggable target for efficient management of severe co-morbid TB COVID-19 condition. However, the selective enrichment of angiogenesis and vascular permeability in observed in the lungs of SARS-CoV-2 infected macaques is not seen in animal models, or patients of TB. These results have the potential to generate additional, specific druggable targets for COVID-19.Overall, we interrogated transcriptional profiles of lungs from juvenile and old macaques infected with SARS-CoV-2. This study has provided fundamentally new information on the host response in young and old macaques infected with SARS-CoV-2. This model provides relevant insights necessary for further vaccine and therapeutic development for COVID-19.MethodsMacaquesAll of the infected animals were housed in Animal Biosafety Level 3 (ABSL3) at the Southwest National Primate Research Center, Texas Biomedical Research Institute, where they were treated per the standards recommended by AAALAC International and the NIH Guide for the Care and Use of Laboratory Animals. Sham controls were housed in ABSL2. The animal studies in each of the species were approved by the Animal Care and Use Committee of the Texas Biomedical Research Institute and as an omnibus Biosafety Committee protocol.Animal studies, and tissue harvest for RNA sample preparationRhesus macaques (Macaca mulatta) animals enrolled in this study have been described in detail16, and the infection of these animals with 1.05 × 106 pfu SARS-CoV-2 isolate USA-WA1/2020 (BEI Resources, NR-52281, Manassas, VA) has also been described earlier16. Control (SARS-CoV-2 uninfected) samples were obtained from opportunistic necropsies conducted on rhesus macaques from the same colony in the past few months. Infected animals were euthanized for tissue collection at necropsy, including lung specimens. Lung tissue from three juvenile (3-yrs old) and five old (average 17-yrs old) rhesus macaques (Supplementary Data 1) were homogenized, snap-frozen in RLT buffer, and DNAse-treated total RNA was extracted using the Qiagen RNeasy Mini kit (Qiagen) for RNA-seq analysis as described earlier50.Viral RNA determinationSARS-CoV-2 RNA isolation and measurement of viral RNA in lung homogenates using RTqPCR have been described16.RNA-sequencing and analysiscDNA libraries were prepared from RNA samples using the Clontech SMARTer universal low input RNA kit to maximize yield, and samples were sequenced on Illumina NovaSeq S4 XP (paired 150 bp reads). After adapter trimming using Trimmomatic v0.3977, sequenced RNA-seq reads were aligned to the Macaca mulatta genome (version 10, Ensembl release 10078) using the STAR aligner v2.7.3a79 (2-pass mode, basic). All raw RNA-Seq fastq files were uploaded to the NCBI Sequence Read Archive (SRA80), and complete sample metadata and accession information are provided in Supplementary Data 1. Read fragments (read pairs or single reads) were quantified per gene per sample using featureCounts v1.5.181. Significantly differentially expressed genes between naive, controller, and progressor sample sets were identified using DESeq2 v1.4.517 with default settings, and a minimum P value significance threshold of 0.01 (after False Discovery Rate [FDR82] correction for the number of tests). PCA also was calculated using DESeq2 output (default settings, using the top 500 most variable genes). FPKM (fragments per kilobase of gene length per million reads mapped) normalization was performed using DESeq2-normalized read counts. Pathway enrichment analysis among differentially expressed gene sets of interest was performed for (a) Reactome18 pathways, using the human orthologs as input into the WebGestalt83 web server (p ≤ 0.05 after FDR correction, minimum 3 genes per term) and (b) KEGG19 pathways and Gene Ontology20 terms, using the g:profiler web server84 which has a database of these annotations matched to macaque ENSEMBL gene IDs (p ≤ 0.05 after FDR correction, minimum 3 genes per term). Mapped fragment counts, relative gene expression levels, gene annotations, and differential expression data for every macaque gene are available in Supplementary Data 2, along with orthology matches to human genes retrieved from ENSEMBL78 and identifications of differentially expressed (DE) genes belonging to enriched pathways of interest, for genes of interest in Supplementary Data 3, and significant functional enrichment for Reactome, KEGG and Gene Ontology pathways, among differentially gene sets of interest in Supplementary Data 4. In addition, genes significantly differentially regulated during the progression of tuberculosis (in both the macaque gene and the corresponding mouse ortholog) were identified from a previous transcriptomic study of tuberculosis-infected lung tissue50, and the upregulated and downregulated gene sets were intersected with the COVID-19 results from the current study.ImmunohistochemistryFormalin-fixed paraffin-embedded non-human primate lung tissues were sectioned into 5-μm thick cuts. Baking of sections was done at 65 °C with subsequent de-paraffinization using Xylene and re-hydration with decreasing gradations of ethanol and dH2O. Heat-induced epitope retrieval method was used for antigen unmasking using Sodium citrate buffer (10 mM, pH 6.0) for 20 min at 95 °C. Sections were then blocked using 3% BSA in TBST at RT for 45 min. SARS CoV-2 Spike rabbit polyclonal antibody (ProSci, USA, 1:200, 2 h at 37 °C) and Mouse anti-human CD66abce-APC conjugated (Miltenyi Biotech, USA, 1:20, 2 h at 37 °C) were used for the detection of SARS CoV-2 virus and neutrophils, respectively. Goat anti-rabbit IgG (H + L), Alexa Fluor 488 conjugate secondary antibody (Thermo Fisher Scientific, USA, 1:400, 1 h at 37 °C) was used for labeling SARS CoV-2 Spike primary antibody. Further, specimens were incubated with DAPI for nuclear staining (Thermo Fisher Scientific, USA, 1:5000, 5 min at 37 °C) followed by mounting with Prolong Diamond Antifade mountant (Thermo Fisher Scientific, USA). The stained sections were then visualized using Ziess LSM 800 confocal microscope (×63 magnification).Human sample collectionPlasma samples were collected from COVID-19 patients that attended the emergency room of the Instituto Nacional de Ciencias Médicas y Nutrición Salvador Zubirán (INCMNSZ), and the Instituto Nacional de Enfermedades Respiratorias Ismael Cosío Villegas (INER) in Mexico City, from March to June of 2020. Detection of SARS-CoV-2 was performed by real-time polymerase chain reaction (RT-PCR) in swab samples, bronchial aspirates (BA), or bronchoalveolar lavage (BAL). For this purpose, viral RNA was extracted from clinical samples with the MagNA Pure 96 system (Roche, Penzberg, Germany). The RT-PCR reactions were performed in a total volume of 25 μL, containing 5 μL of RNA, 12.5 μL of 2× reaction buffer provided with the Superscript III one-step RT-PCR system with Platinum Taq Polymerase (Invitrogen, Darmstadt, Germany; containing 0.4 mM of each deoxyribose triphosphates (dNTP) and 3.2 mM magnesium sulfate), 1μL of reverse transcriptase/ Taq mixture from the kit, 0.4 μL of a 50 mM magnesium sulfate solution (Invitrogen), and 1 μg of nonacetylated bovine serum albumin (Roche). All oligonucleotides were synthesized and provided by Tib-Molbiol (Berlin, Germany). Thermal cycling was performed at 55 °C for 10 min for reverse transcription, followed by 95 °C for 3 min and then 45 cycles of 95 °C for 15 s, 58 °C for 30 s. Primer and probe sequences are as follows: RdRP gene [RdRp-SARSr-F:GTGARATGGTCATGTGTGGCGG,RdRp-SARSr-P2: FAMCAGGTGGAACCTCATCAGGAGATGCBBQ,RdRP_SARSrP1:FAMCCAGGTGGWACRTCATCMGGTGATGCBBQ,RdRp_SARSrR:CARATGTTAAASACACTATTAGCATA], E gene [E_Sarbeco_F:ACAGGTACGTTAATAGTTAATAGCGT,E_Sarbeco_P1:FAMACACTAGCCATCCTTACTGCGCTTCGBBQ,E_Sarbeco_R:ATATTGCAGCAGTACGCACACA], N gene [N_Sarbeco_F:CACATTGGCACCCGCAATC,N_Sarbeco_P1:FAMACTTCCTCAAGGAACAACATTGCCABBQ, N_Sarbeco_R:GAGGAACGAGAAGAGGCTTG]. Clinical and demographic data were retrieved from the medical records of all participants. These data included age, gender, anthropometrics, comorbidities, symptoms, triage vital signs, and initial laboratory test results. Initial laboratory tests were defined as the first test results available (typically within 24 h of admission) and included white blood cell counts (WBC), neutrophil, and lymphocyte counts (Supplementary Data 5). For human studies, written informed consent from participants under protocols approved by the Ethics Committee of the Instituto Nacional de Enfermedades Respiratorias Ismael Cosio Villegas, Mexico. The study was approved by the Institutional Research Committee of Instituto Nacional de Enfermedades Respiratorias Ismael Cosio Villegas, Mexico.Cytokine levels in human plasma samplesPeripheral blood samples were obtained from all participants at hospital admission. Plasma levels of interferon-gamma (IFN-γ) and VEGF, were determined by Luminex assays using the Luminex platform Bio-Plex Multiplex 200 (Bio-Rad Laboratories, Inc., Hercules, CA, USA). Plasma samples from four healthy volunteer donors were used as controls.Statistics and reproducibilityWhere applicable, comparisons between two groups were performed using the unpaired two-tailed Student t-test, using the GraphPad Prism 5 software. In all instances, P values ≤ 0.05 were considered as significant. DESeq2 (version 1.4.5)17 with default settings, and a minimum P value significance threshold of 0.01 (after False Discovery Rate [FDR82] correction for the number of tests was used for all RNA-Seq comparisons of differential expression.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article.Supplementary informationSupplementary InformationDescription of Additional Supplementary FilesSupplementary Data 1Supplementary Data 2Supplementary Data 3Supplementary Data 4Supplementary Data 5Supplementary Data 6Reporting Summary
Viral infection Interferons