key: cord-0714912-zhpvg9t6
authors: Hasan, Md Zobaer; Islam, Syful; Matsumoto, Kenichi; Kawai, Taro
title: SARS-CoV-2 infection initiates interleukin-17-enriched transcriptional response in different cells from multiple organs
date: 2021-08-19
journal: Sci Rep
DOI: 10.1038/s41598-021-96110-3
sha: d74cac898e1b073add1da15b358a1db2828b104d
doc_id: 714912
cord_uid: zhpvg9t6

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection has emerged as a pandemic. Paucity of information concerning the virus and therapeutic interventions have made SARS-CoV-2 infection a genuine threat to global public health. Therefore, there is a growing need for understanding the molecular mechanism of SARS-CoV-2 infection at cellular level. To address this, we undertook a systems biology approach by analyzing publicly available RNA-seq datasets of SARS-CoV-2 infection of different cells and compared with other lung pathogenic infections. Our study identified several key genes and pathways uniquely associated with SARS-CoV-2 infection. Genes such as interleukin (IL)-6, CXCL8, CCL20, CXCL1 and CXCL3 were upregulated, which in particular regulate the cytokine storm and IL-17 signaling pathway. Of note, SARS-CoV-2 infection strongly activated IL-17 signaling pathway compared with other respiratory viruses. Additionally, this transcriptomic signature was also analyzed to predict potential drug repurposing and small molecule inhibitors. In conclusion, our comprehensive data analysis identifies key molecular pathways to reveal underlying pathological etiology and potential therapeutic targets in SARS-CoV-2 infection.

www.nature.com/scientificreports/ Understanding the hyperinflammation and type of immune response is of utmost importance to design an effective therapy. The innate immune system mounts the first line of defense against viral infection at an early stage. Innate sensing of viral material initiates antiviral response by producing type I interferon (IFN) through interferon regulatory factors (IRFs) and elicits a pro-inflammatory cytokine response via NF-κB-dependent pathways 16 . Reportedly, monocytes from elderly subjects showed diminished IFNα/β level while producing pro-inflammatory cytokines upon stimulation 17 . This suggests that aging specifically impairs the IFN activation but not pro-inflammatory cytokine production. Indeed, COVID-19 and SARS-CoV patients exhibited hypercytokinemia 18, 19 with an aberrant IFN response. Hence, a more comprehensive and concerted study involving different in vitro or in vivo models are required to understand the interplay between virus and host innate immune response.

Here, we sought to characterize spatial distribution of host transcriptome in different cells upon SARS-CoV-2 infection and whether the infection generates unique and distinctive transcriptomic signature. In this regard, we have analyzed and integrated publicly available gene expression datasets. First, we have analyzed RNA-seq data from SARS-CoV-2-infected normal human bronchial epithelial (NHBE) cells, A549 cells (human lung carcinoma), primary human airway epithelial (pHAE) cultures, cardiomyocytes and liver organoids. To further understand, we have compared these data with other respiratory viruses such as SARS-CoV, MERS, Influenza A virus (IAV), and Respiratory Syncytial Virus (RSV). SARS-CoV-2 infection elicited a clearly distinctive host response in multiple cells compared with those observed in other viral infections. A computational analysis was carried out to predict potential drug candidates based on their ability to reverse SARS-CoV-2-mediated host transcriptional signature.

Dataset and sample information. The NCBI Gene Expression Omnibus (GEO) database (http:// www. ncbi. nlm. nih. gov/ geo/) has been queried using keywords "SARS-CoV-2" and "COVID-19". The datasets were selected with stringent selection of data generated on "Homo-sapiens". Finally, following datasets were included for this study: GSE147507 20 , GSE150392 21 , GSE153970 22 and GSE151803 23 . Briefly, three independent biological replicates of primary human lung/airway epithelial cells or A549 cells that were mock treated or infected with SARS-CoV-2 at a Multiplicity Of Infection (MOI) of 2 (GSE147507) or 0.25 (GSE153970), for 24 h and 48 h respectively. A low MOI of 0.1 was used for 72 h in GSE150392 using iPSC derived cardiomyocytes. Liver organoid infections were performed at a MOI of 0.1 for 24 h in GSE 151,803. SARS-CoV-2 (USA-WA1/2020) was used in all these studies. In addition to SARS-CoV-2, GSE47960 and GSE100504 have been used to analyze host response to SARS-CoV and MERS viruses respectively. Primary human lung/airway epithelial cells were infected at an MOI of 2 and 5 for SARS-CoV and MERS, respectively in these studies. GSE147507 dataset was used for RSV-and IAV-infected host transcriptome. A complete list containing cell types, MOI, infection period has been prepared as supplementary table 1 (see supplementary information file).

Data processing and analysis. Raw sequencing reads were quality controlled and aligned to the human genome (hg19) using the RNA-Seq Alignment App on Basespace (Illumina, CA, USA), followed by differential expression analysis using DESeq2 24 . Differentially expressed genes (DEGs) were characterized for each sample (p adjusted-value < 0.05 and log2FC > 1). Volcano plots were constructed using custom scripts in R. A large volume of DEGs were identified with low count and highly dispersed genes, hence shrinkage estimator 'apeglm' 25 was used in order to remove noise and visualize the DEGs with significant differences in mock treated vs SARS-CoV-2 treated group. Particularly this statistical approach preserves true, large differences in log fold change (LFC) across conditions and is superior to common methods in ranking of genes by LFC in presence of low counts. In brief, apeglm utilizes heavy-tailed Cauchy distribution instead of normal distribution on the effect sizes, with fixed shape and scale adapted to the distribution of observed maximum likelihood estimates (MLE) for all genes. Then a LaPlace approximation was used to provide shrinkage estimates and corresponding standard deviation. This eliminated further needs for filtering rules or pseudo counts and maximizes the power of the current data to estimate the effect size for each gene.

A web-based tool, Metascape 26 was used for functional enrichment of DEGs. Metascape queries publicly available databases, e.g. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes (KEGG) 27 and assigns DEGs to their respective enriched pathways by calculating the pairwise similarity between any two terms. Hypergeometric test and Benjamini-Hochberg p-value correction algorithm have been used to identify statistically significant enriched ontology terms. The Bioconductor cluster Profiler package v3.14.3 28 and STRING v11 29 have also been used to develop gene to pathway and protein-protein interaction (PPI) networks of upregulated DEGs and Cytoscape application v3.7.1 30 was used for plotting genes based on the count and p-value.

Disease similarities and drug predictive analysis. DisGeNET database 31 was employed to identify disease similarities with enriched DEGs. A minimum count of 3 genes, with a p-value < 0.01 and an enrichment factor > 1.5 (the enrichment factor is the ratio between the observed counts and the counts expected by chance) were collected and grouped into clusters based on their membership similarities. Accumulative hypergeometric distribution was employed for calculating p-value. The L1000FDW web tool 32 was utilized for potential drug candidate search for the treatment of SARS-CoV-2 infection. L1000FWD computes similarity between a given set of genes with the Library of Integrated Network-based Cellular Signatures (LINCS)-L1000 data and predicts compounds or drugs that reverse the input transcriptomic signature. 

DEGs of SARS-CoV-2 infection. In order to determine and compare host cell response against SARS-CoV-2 infection, we first analyzed various publicly available RNA-sequencing datasets. We identified numerous DEGs in lung, liver and heart cells (Fig. 1A ). Due to large number of DEGs, we employed more stringent parameter and used 'apeglm' algorithm in Deseq2 program to identify highly differential genes. This significantly narrowed down the gene list and interestingly most of the DEGs belong to an upregulated category (Fig. 1A) . We could not find any significant association among those few downregulated genes. For subsequent analysis, mostly upregulated genes have been used. Notably, chemokine genes (CXCL1, CXCL3, CXCL5, CXCL8, CCL20) and ISGs (IFITM1, IFI44L, IFI27) dominated this list (See supplementary file 1 for the gene lists). An MA plot showing A549 cells had very few DEGs compared with other cells used in the study (See Figure  S1 ). However, supplementation with putative SARS-CoV-2 receptor ACE2 in A549 cells resulted in higher number of modulated transcriptome (see Figure S2A ). A circos plot showing association of DEGs among various datasets. Particularly, NHBE cells share many DEGs with pHAE cells (Fig. 1B) . Then, the DEGs were screened for their approximate association in disease-gene networks using DisGeNET database. Strikingly, enrichment of our upregulated genes could link most of the clinical features such as pneumonia, influenza, myocardial ischemia, hemorrhagic shock, etc. observed in COVID-19 patients. (Fig. 1C ).

Functional annotation and pathway enrichment of DEGs. Next, we performed gene set functional enrichment analysis using KEGG pathway. Circos plot showing many distinct genes in one dataset are ontologically connected to the ontological features in another dataset ( Fig. 2A ). Key pathways shared by genes from most of the datasets were cytokine-cytokine receptor signaling, IL-17 signaling pathway, NOD-like receptor-mediated signaling, and Measles (Fig. 2B ). Interestingly, very few genes were upregulated in A549 cells and were not enriched into IL-17 signaling pathway. We then analyzed transcriptome of A549 cells supplemented with putative SARS-CoV-2 receptor, ACE2. ACE2 receptor supplementation induced expression of transcripts that were shared with other cells (See Figure S2A ) and enriched in IL-17-mediated signaling pathways (See Figure S2B ).

DEGs were then ranked based on their frequency of distribution in datasets and number of associated pathways. We selected top 50 upregulated genes that were present in at least two different cells. Pathway enrichment analysis based on gene count assigned in KEGG pathway for upregulated DEGs revealed IL-17 signaling, NOD-like receptor-mediated signaling, and the TNF signaling (Fig. 2C ) amongst the most regulated pathways www.nature.com/scientificreports/ by these genes. Among other statistically significant enriched terms, pathways related to NF-κ B and rheumatoid arthritis were also evident (Fig. 2C ). In addition, analyzed data from BALF collected from COVID-19 patients in China (Xiong et al. 2020) followed a similar trend of enrichment identified from our analysis (See Figure S3 ). However, we did not observe strong similarities with our host transcriptome data when PBMC transcripts were used for analysis from the same study group in China. (Fig. 3A,B ). Next, we compared genes that were upregulated in SARS-CoV, SARS-CoV-2-and MERS-infected cells (Fig. 3C ). Enrichment of 40 shared genes showed 'Measles' , 'NOD-like receptor signaling' , 'TLR signaling' , 'RLR signaling' , and 'TNF signaling' (Fig. 3D ). IL-17 signaling was also among the enriched pathway, however it was not as strongly activated as we observed in all the datasets of SARS-CoV-2 infection. Next, we performed a gene set enrichment of all uniquely upregulated DEGs for each of the infections. We have observed during time course infection of SARS-CoV and MERS, grouping all DEGs from different time points enriched into TLR and RLR signaling pathways, along with the NF-κ B and MAPK pathways. Of note, SARS-CoV-2-infected transcripts did not trigger RLR signaling while induced strong cytokine responses mediated by NF-κ B signaling (Fig. 4A) . On the other hand, SARS-CoV-2-infected unique transcriptomes predominantly enriched into 'IL-17 signaling pathway' while both SARS-CoV and MERS-infected unique transcripts did not. This data indicates that SARS-CoV-2 infection induces a plethora of chemokines that mediates a strong inflammatory response driven by IL-17 signaling (Fig. 4B) .

We then sought to relate these unique IL-17 specific transcripts to the biological processes. Notably, 'leukocyte aggregation' , 'positive regulation of granulocytes/neutrophils' , and 'fever generation' were the top biological processes modulated by these transcripts (Fig. 4C ). This demonstrated striking similarities between phenotypes observed in COVID-19 patients and biological processes controlled by the IL-17 specific transcripts. (Fig. 5A) . Further analysis showed that SARS-CoV-2 transcriptional regulation was mainly controlled by NF-κ B (RELA) while IRF mediated signaling was not activated. On the contrary, both SARS-CoV and MERS mediated transcriptional responses triggered IRF mediated signaling cascade (Fig. 5B ). Next, we developed protein-protein interaction (PPI) networks with the upregulated DEGs in respiratory infections using cytoscape. KEGG categories "IL-17 signaling", "Influenza A signaling" and "TNF signaling" formed most densely connected subnetworks, which are mainly composed of ribosomal proteins and chemokines, respectively. Most of the genes in IL-17 signaling were representative of SARS-CoV-2-infected transcriptome as denoted by color (Fig. 5C ). On the other hand, Influenza A signaling pathway was shared by most of the respiratory viruses and possibly produce similar phenotypes observed in those respiratory illness mediated by viral infection.

Potential therapeutic intervention based on transcriptome signature. Reverse signature perturbation analysis was performed using the upregulated DEGs identified upon SARS-CoV-2 infection (Fig. 6 ). Potential drugs have been clustered according to their ability to reverse the upregulated DEGs using the webbased tool L1000CDS2. Among the identified drugs; Saracatinib showed in vitro inhibition of MERS replication 33 . Besides, Dasatinib and Imatinib from our list have also been reported to inhibit SARS-CoV and MERS infection in Vero E6 cells 34 . This suggests that the repurposing of drugs based on host transcriptome may be an effective strategy until a safe and effective vaccine is developed.

Viruses utilize host cell machinery and modulate host cell transcriptome either for their replication or for evading host immune responses. Hence, transcriptional profiling of host cell would help us realize changes in the genetic landscape during viral infection. In this regard, we have analyzed differentially expressed genes in SARS-CoV-2-mediated infection using various in vitro cellular systems. Datasets used in our study are mostly obtained using various cells from lungs along with cardiomyocytes and liver organoids. Despite being a respiratory virus, SARS-CoV-2 has been reported to impair liver function 35 and has active transcription site in cardiomyocytes www.nature.com/scientificreports/ of deceased patients 36 . Notably, ACE2 expression is high in cardiomyocytes 37 and in general the expression is relatively higher in liver and heart when compare with lungs 38 . Thus, it is important to carry out a comparative transcriptional response in cells from different organs. We observed a consistent enrichment of cytokine response in all the datasets used in this meta-analysis. Notably, monocyte or neutrophil recruiting chemokines such as CCL20, CXCL1, CXCL3 and CXCL8 have been upregulated. This data is in sync with the elevated level of circulating neutrophils observed in COVID-19 patients along with a decrease in lymphocytes [39] [40] [41] .

Pattern recognition receptors such as Toll-like receptors (TLRs) and RIG-I like receptors (RLRs) sense viral single or double-stranded RNA to produce inflammatory cytokines and type I IFN 42 . We observed that IRF7 and IRF9 were upregulated in SARS-CoV-2 infection. This in turn may stimulate production of interferon stimulated genes (ISGs) to produce an antiviral state. Indeed, several ISGs such as IFI1-3, IFI27, IFITM3, OAS1 and OAS3 have been upregulated. This possibly indicates initial host cellular response in order to inhibit viral replication. However, severe viral load may overwhelm type I IFN response and dictate outcome of infection. Of note, SARS-CoV and MERS-CoV infection exhibited similar outcome depending on the timing of type I IFN production. A delayed type I IFN production results in higher viral replication, cellular damage and production of cytokine storm, whereas an earlier production of type I IFN could protect from lethal effect 43, 44 .

COVID-19 patients with severe respiratory symptoms had higher levels of Tumor necrosis factor-α (TNF-α) and IL-6 45 . It is worth mentioning that IL-6 was not only expressed in SARS-CoV-2-infected cells but was also involved with most of the enriched pathways in our analysis. A recent study 45, 46 showed effective inhibition of cytokine storm by blocking IL-6 with monoclonal antibody Tocilizumab, supporting the notion that IL-6-mediated inflammation contributes to the disease severity. Both IL-6 and IL-17 have been reported with protecting virally infected cells from apoptosis and thereby promoting viral persistence 47 .

All the SARS-CoV-2 datasets used in this study showed particularly strong inflammatory response triggered with IL-17 activation. Previous studies observed elevated IL-17 expression was conjoined with impaired IRF7 upregulation and IFNα responses in herpes simplex virus 2 infection of aged mice. Indeed, IL-17 inhibition could prevent death of aged mice during viral infection 48, 49 . Our data analysis indicated the importance of IL-17 in SARS-CoV-2 infection. Meanwhile IL-17 has also been reported to be increased in critically ill patients compared with non-severe patients 3 . However, whether this upregulation is the direct effect of SARS-CoV-2 infection or not remains unclear. A very recent study published in April 2021 showed that SARS-CoV-2 open reading frame 8 (ORF8) binds to IL17 receptor and activates IL17 signaling pathway. Blocking IL17RA with an antibody reduced IL17 mediated inflammation in lung and liver in SARS-CoV2 ORF8 pseudovirus infected It is still obscure how viruses induce a proinflammatory response whilst keeping immune homeostasis at bay. SARS-CoV-2 or H1N1 2009 infections showed higher IL-1RA and IL-6 levels along with very low viral load in the lung in severely ill people 54 . This raises an interesting possibility as to whether viruses blindside the innate immune system by maintaining a low amplification or self-replication state while causing local inflammation or tissue damage. Our comparative transcriptional analysis of SARS-CoV-2 and other viruses supports this notion. As evident from our observation, NF-κB is the most significant transcription factor that modulates host transcriptional responses upon SARS-CoV-2 infection. Despite producing ISGs, an absence of IRF3 or IRF7 regulated transcriptional processes was notable in SARS-CoV-2 infection. On the other hand, SARS-CoV and MERS transcriptomics were at least partly regulated by IRF1, IRF3 and IRF7. Of note, a recent work showed that SARS-CoV-2 proteases could cleave IRF3 directly and resulted in blunted IFN production 55 . This accumulating evidence suggests that SARS-CoV-2 launches a pro-inflammatory response while specifically blocking antiviral responses.

Despite showing a variable clinical spectrum, severe COVID-19 patients commonly exhibited shortness of breath, and production of pro-inflammatory cytokines such as IL-6, IL-8, IL-1β, IL-1RA 3,56 . Classic clinical symptoms predominantly indicate lung associated phenotypes, however recent evidence indicated that mortality rate is higher in cardiac patients 57, 58 . Our gene set enrichment also showed that SARS-CoV-2 infected cardiomyocytes have higher expression of genes that may induce myocardial ischemia and cerebral arterial infarction.

A consistent observation in each of our SARS-CoV-2 datasets was the presence of chemokines and inflammatory cytokines. Hence, we used a reverse transcriptome signature approach for predicting drugs that are approved by FDA or in clinical trials. This analysis predicted that Dasatinib, CGP-60474, Canertinib, Alvocidib, Saracatinib etc. can perturb and reverse host transcriptional signature of SARS-CoV-2 infection. Reportedly, CGP-60474 could inhibit IL-6 and TNFα to alleviate LPS-induced sepsis in mice 59 .

We have observed distinct transcriptional responses mediated by SARS-CoV-2 in cell culture-based studies. However, MERS and IAV datasets used in our study have different MOI than that of SARS-CoV and SARS-CoV-2. MERS have higher MOI but similar infection period of 24-48 h than two other coronaviruses used in the analysis. On the other hand, IAV infection was done at slightly higher MOI than SARS-CoV or SARS-CoV-2. Interestingly, MERS infection showed late viral peak compare with SARS-CoV and SARS-CoV-2 60-62 , whilst IAV showed faster inflammatory and antiviral responses than SARS-CoV-2 63 . This clearly explains the reasoning of using slightly higher MOI for MERS and shorter infection period with IAV. www.nature.com/scientificreports/ Despite comprehensive analysis of the publicly available datasets, there may be some limitations as well. As the percentage of lethality is higher in aged population, future studies using aged or senescent cells would clarify the discrepancies between young and aged in terms of infection and mortality rate. In addition, all the studies used a single cell type in the absence of immune cells, which may not reflect true intercellular communications take place in vivo. As activation of immune system varies in age groups, further research should be done using a multicellular system. As we have seen that these cellular models of infection secret numerous chemokines and particularly able to trigger immune cell activation and chemotaxis, future studies should be directed towards validating our data in moderate and severe COVID-19 patients.

Taken together, our comprehensive data analysis showed that SARS-CoV-2 initiated a distinct IL-17-driven inflammatory response irrespective of the cells used in various studies (Fig. 7) . Overall, our data analysis showed unique and distinctive SARS-CoV-2-induced transcriptional responses compared with other respiratory viruses.

The datasets and code used in this study are publicly available and comprehensively described in the methods section. Any other necessary information would be provided upon reasonable request. www.nature.com/scientificreports/

Neurologic Manifestations of Hospitalized Patients With Coronavirus Disease

COVID-19 and the cardiovascular system

Clinical features of patients infected with 2019 novel coronavirus in Wuhan

Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus-infected pneumonia in Wuhan

Expression of the SARS-CoV-2 cell receptor gene ACE2 in a wide variety of human tissues

Hospital outbreak of Middle East respiratory syndrome coronavirus

Clinical features and short-term outcomes of 144 patients with SARS in the greater Toronto area

Transcriptomic characteristics of bronchoalveolar lavage fluid and peripheral blood mononuclear cells in COVID-19 patients

Overly exuberant innate immune response SARS-CoV-2 infect

Association of a Province-Wide Intervention With Salt Intake and Hypertension in Shandong Province

SARS-associated viral hepatitis caused by a novel coronavirus: report of three cases

Histopathology of Middle East respiratory syndrome coronovirus (MERS-CoV) infection-Clinicopathological and ultrastructural study

Liver injury in COVID-19: management and challenges

Myocardial injury characterized by elevated cardiac troponin and in-hospital mortality of COVID-19: an insight from a meta-analysis

Prevalence and impact of cardiovascular metabolic diseases on COVID-19 in China

Innate immune recognition of viral infection

Mx1 reveals innate pathways to antiviral resistance and lethal influenza disease

Impaired type I interferon activity and inflammatory responses in severe COVID-19 patients

Expression profile of immune response genes in patients with Severe Acute Respiratory Syndrome

SARS-CoV-2 launches a unique transcriptional signature from in vitro, ex vivo, and in vivo systems

Human iPSC-derived cardiomyocytes are susceptible to SARS-CoV-2 infection

Type I and Type III IFN Restrict SARS-CoV-2 infection of human airway epithelial cultures

Identification of candidate COVID-19 therapeutics using hPSC-derived lung organoids

Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2

Heavy-tailed prior distributions for sequence count data: Removing the noise and preserving large differences

Metascape provides a biologist-oriented resource for the analysis of systems-level datasets

Kyoto Encyclopedia of Genes and Genomes

STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets

Cytoscape: a software environment for integrated models of biomolecular interaction networks

The DisGeNET knowledge platform for disease genomics: 2019 update

L1000FWD: Fireworks visualization of drug-induced transcriptomic signatures

Saracatinib inhibits middle east respiratory syndrome-coronavirus replication in vitro

Repurposing of clinically developed drugs for treatment of Middle East respiratory syndrome coronavirus infection

SARS-CoV-2 infection of the liver directly contributes to hepatic impairment in patients with COVID-19

Evidence of SARS-CoV-2 transcriptional activity in cardiomyocytes of COVID-19 patients without clinical signs of cardiac involvement

Mediators of SARS-CoV-2 entry are preferentially enriched in cardiomyocytes

Expression of the SARS-CoV-2 cell receptor gene ACE2 in a wide variety of human tissues

Characterization of an asymptomatic cohort of SARS-COV-2 infected individuals outside of Wuhan, China

Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: A descriptive study

Functional exhaustion of antiviral lymphocytes in COVID-19 patients

The role of pattern-recognition receptors in innate immunity: update on Toll-like receptors

Dysregulated type I interferon and inflammatory monocyte-macrophage responses cause lethal pneumonia in SARS-CoV-infected mice

IFN-I response timing relative to virus replication determines MERS coronavirus infection outcomes

Complex immune dysregulation in COVID-19 patients with severe respiratory failure

Tocilizumab treatment in COVID-19: a single center experience

Interleukin-6 (IL-6) and IL-17 synergistically promote viral persistence by inhibiting cellular apoptosis and cytotoxic T cell function

Aging promotes neutrophil-induced mortality by augmenting IL-17 production during viral infection

Aging impairs IFN regulatory factor 7 upregulation in plasmacytoid dendritic cells during TLR9 activation

ORF8 contributes to cytokine storm during SARS-CoV-2 infection by activating IL-17 pathway

SARS-CoV-2 genomic surveillance in Taiwan revealed novel ORF8-deletion mutant and clade possibly associated with infections in Middle East

Effects of a major deletion in the SARS-CoV-2 genome on the severity of infection and the inflammatory response: an observational cohort study

Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: A descriptive study

Role of host immune response and viral load in the differential outcome of pandemic H1N1 (2009) influenza virus infection in Indian patients

SARS-CoV-2 proteases cleave IRF3 and critical modulators of inflammatory pathways (NLRP12 and TAB1): implications for disease presentation across species and the search for reservoir hosts

Potential Effects of Coronaviruses on the Cardiovascular System: A

The Variety of Cardiovascular Presentations of COVID-19

Association of Cardiac Injury With Mortality in Hospitalized Patients With COVID-19 in Wuhan

LINCS L1000 dataset-based repositioning of CGP-60474 as a highly potent anti-endotoxemic agent

Clinical progression and viral load in a community outbreak of coronavirus-associated SARS pneumonia: a prospective study

Viral load kinetics of MERS coronavirus infection

Temporal profiles of viral load in posterior oropharyngeal saliva samples and serum antibody responses during infection by SARS-CoV-2: an observational cohort study

Infection of human Nasal Epithelial Cells with SARS-CoV-2 and a 382-nt deletion isolate lacking ORF8 reveals similar viral kinetics and host transcriptional profiles

/scientificreports/ Reprints and permissions information is available at www.nature.com/reprints

This work was supported by KAKENHI Grants-in-Aid for Research Activity 20H03468 (T. Kawai). We sincerely thank Amy Leigh Harding, Ph.D. of University of Sheffield for her critical advises on the manuscript. 

The authors declare no competing interests.

The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598-021-96110-3.Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.