Comprehensive Comparison of RNA-Seq Data of SARS-CoV-2, SARS-CoV and MERS-CoV Infections: Alternative Entry Routes and Innate Immune Responses

Background The pathogenesis of COVID-19 emerges as complex, with multiple factors leading to injury of different organs. Some of the studies on aspects of SARS-CoV-2 cell entry and innate immunity have produced seemingly contradictory claims. In this situation, a comprehensive comparative analysis of a large number of related datasets from several studies could bring more clarity, which is imperative for therapy development. Methods We therefore performed a comprehensive comparative study, analyzing RNA-Seq data of infections with SARS-CoV-2, SARS-CoV and MERS-CoV, including data from different types of cells as well as COVID-19 patients. Using these data, we investigated viral entry routes and innate immune responses. Results and Conclusion First, our analyses support the existence of cell entry mechanisms for SARS and SARS-CoV-2 other than the ACE2 route with evidence of inefficient infection of cells without expression of ACE2; expression of TMPRSS2/TPMRSS4 is unnecessary for efficient SARS-CoV-2 infection with evidence of efficient infection of A549 cells transduced with a vector expressing human ACE2. Second, we find that innate immune responses in terms of interferons and interferon simulated genes are strong in relevant cells, for example Calu3 cells, but vary markedly with cell type, virus dose, and virus type.


INTRODUCTION
Coronaviruses are non-segmented positive-sense RNA viruses with a genome of around 30 kilobases (1,2). The genome has a 5' cap structure along with a 3' poly (A) tail, which acts as mRNA for translation of the replicase polyproteins (1,2). The replicase gene occupies approximately two thirds of the entire genome and encodes 16 non-structural proteins (nsps) (1,2). The remaining third of the genome contains open reading frames (orfs) that encode accessory proteins and four structural proteins, including spike (S), envelope (E), membrane (M), and nucleocapsid (N) (1,3).
Over the past 20 years, three epidemics or pandemics of lifethreatening diseases have been caused by three closely related coronavirusessevere acute respiratory syndrome coronavirus (SARS-CoV), which emerged with nearly 10% mortality (4) 1 in 2002-2003 and spread to 26 countries before being contained; Middle East respiratory syndrome coronavirus (MERS-CoV), with mortality around 34% (5) 2 starting in 2012 and since then spreading to 27 countries; SARS-CoV-2, emerging in late 2019 (6), which has caused many millions of confirmed cases and > 2 million deaths by March 17, 2021 3 . SARS-CoV, MERS-CoV or SARS-CoV-2 can cause mild to severe forms of respiratory illness with symptoms including fever, cough, and shortness of breath (7). 4 SARS-CoV-2 is a new coronavirus, but its similarity to SARS-CoV [amino acid sequences about 76% identical (8)] and MERS-CoV suggests comparisons to these earlier epidemics. Despite the difference in the total number of cases caused by SARS-CoV and SARS-CoV-2 due to different transmission rates, the outbreak caused by SARS-CoV-2 resembles the outbreak of SARS: both emerged in winter and were linked to exposure to wild animals sold at markets. Although MERS-CoV has high morbidity and mortality rates, lack of autopsies from MERS-CoV cases has hindered our understanding of MERS-CoV pathogenesis in humans.
In response to COVID-19, the public health measures such as face mask wearing, physical distancing, contact tracing, case isolation and quarantine are generally adopted to halt transmission (9)(10)(11)(12). Recently, SARS-CoV-2 vaccines have been developed to control COVID-19 (13,14). However, specific direct anti-SARS-CoV-2, anti-SARS-CoV or anti-MERS-CoV therapeutics are still highly needed for effective treatment of moderate and severe cases to reduce overall mortality rate (15)(16)(17). There are several points of attack for potential anti-SARS-CoV-2/SARS-CoV/MERS-CoV therapies, e.g. intervention on cell entry mechanisms to prevent virus invasion, or acting on the host immune system to kill the infected cells and thus prevent replication of the invading viruses. A better understanding of virus entry mechanisms and the immune responses can therefore guide the development of novel therapeutics.
Virus entry into host cells is the first step of the viral life cycle. It is an essential component of cross-species transmission and an important determinant of virus pathogenesis and infectivity (18,19), and also constitutes an antiviral target for treatment and prevention (20). It seems that SARS-CoV and SARS-CoV-2 use similar virus entry mechanisms (21). The infection of SARS-CoV or SARS-CoV-2 in target cells was initially identified to occur by cell-surface membrane fusion (22,23). Some later studies have shown that SARS-CoV can infect cells through receptor mediated endocytosis (24,25) as well. Both mechanisms require the S protein of SARS-CoV or SARS-CoV-2 to bind to angiotensin converting enzyme 2 (ACE2), and S protein of MERS-CoV to dipeptidyl peptidase 4 (DPP4) (26), respectively, through their receptor-binding domain (RBD) (27). In addition to ACE2, some recent studies suggest that there are possible other coronavirus-associated receptors and factors that facilitate the infection of SARS-CoV-2 (28,29), including the cell surface proteins Basignin (BSG or CD147) (30), CD209 (31) and the tyrosine-proteinkinase receptor UFO (AXL) (32). Recently, clinical data have revealed that SARS-CoV-2 can infect several organs where ACE2 expression could not be detected in healthy individuals (33,34), which highlights the need of closer inspection of virus entry mechanisms.
The binding of S protein to a cell-surface receptor is not sufficient for infection of host cell (35). In the cell-surface membrane fusion mechanism, after binding to the receptor, the S protein requires proteolytic activation by cell surface proteases like TMPRSS2, TPMRSS4, or other members of the TMPRSS family (23,36,37), followed by the fusion of virus and target cell membranes (27,38). In the alternative receptor mediated endocytosis mechanism, the endocytosed virion is subjected to an activation step in the endosome, resulting in the fusion of virus and endosome membranes and the release of the viral genome into the cytoplasm. The endosomal cysteine proteases cathepsin B (CTSB) and cathepsin L (CTSL) (39) might be involved in the fusion of virus and endosome membranes. Availability of these proteases in target cells largely determines whether viruses infect the cells through cellsurface membrane fusion or receptor mediated endocytosis. How the presence of these proteases impacts efficiency of infection with SARS-CoV-2, SARS-CoV and MERS-CoV, still remains elusive.
When the virus enters a cell, it may trigger an innate immune response, a crucial component of the defense against viral invasion. Compounds that regulate innate immune responses can be introduced as antiviral agents (19). The innate immune system is initialized as pattern recognition receptors (PRRs) such as Toll-like receptors (TLRs) and cytoplasmic retinoic acidinducible gene I (RIG-I) like receptors (RLRs) recognize molecular structures of the invading virus (40,41). This pattern recognition activates several signaling pathways and then downstream transcription factors such as interferon regulator factors (IRFs) and nuclear factor kB (NF-kB) (42)(43)(44) expression of type I (a or b) and type III (l) interferons (IFNs) (42,45,46). IFN-a (IFNA1, IFNA2, etc), IFN-b (IFNB1) and IFN-l (IFNL1-4) are important cytokines of the innate immune responses. IFNs bind and induce signaling through their corresponding receptors (IFNAR for IFN-a/b and IFNLR for IFN-l), and subsequently induceexpression of IFN-simulated genes (ISGs) (e.g. MX1, ISG15 and OASL) and proinflammatory chemokines (e.g. CXCL8 and CCL2) to suppress viral replication and dissemination (47,48). Dysregulated inflammatory host response results in acute respiratory distress syndrome (ARDS), a leading cause of COVID-19 mortality (49).
One attractive therapy option to combat COVID-19 is to harness the IFN-mediated innate immune responses. Clinical trials with type I and type III IFNs for treatment of COVID-19 have been conducted and many more are still ongoing (50,51). In this regard, the kinetics of the secretion of IFNs in the course of SARS-CoV-2 infection needs to be defined. Unfortunately, some results on the host innate immune responses to SARS-CoV-2 are apparently at odds with each other (52-56), e.g. it is unclear whether SARS-CoV-2 infection induces low IFNs and moderate ISGs (52), or robust IFN responses and markedly elevated expression of ISGs (53)(54)(55)(56). This has to be clarified. The use of IFNs as a treatment in COVID-19 is now a subject of debate as well (57). Thus, the kinetics of IFN secretion relative to the kinetics of virus replication need to be thoroughly examined to better understand the biology of IFNs and ISGs in the course of SARS-CoV-2 infection and thus provide guidance to identify the temporal window of therapeutic opportunity.

Data Collection
After the successful release of the virus genome into the cytoplasm, a negative-strand genomic-length RNA is synthesized as the template for replication. Negative-strand subgenome-length mRNAs are formed as well from the virus genome as discontinuous RNAs, and used as the templates for transcription. In the public data we collected for the analysis, there are two main library preparation methods to remove the highly abundant ribosomal RNAs (rRNA) from total RNA before sequencing. One is polyA+ selection, the other is rRNAdepletion (63). It is known that coronavirus genomic and subgenomic mRNAs carry a polyA tail at their 3' ends, so in the polyA+ RNA-Seq, we have (1) virus genomic sequence from virus replication, i.e. replicated genomic RNAs from negativestrand as template, and (2) subgenomic mRNAs from virus transcription; in the rRNA-depletion RNA-Seq we have (1) virus genomic sequence from virus replication: both replicated genomic RNAs from negative-strand as template and the negative-strand templates themselves, and (2) subgenomic mRNAs from virus transcription. PolyA+ selection was used if not specifically stated in this study, "total RNA" is used to specify that the rRNA-depletion method was used to prepare the sequencing libraries.
The raw FASTQ data of different cell types infected with SARS-CoV-2, SARS-CoV and MERS-CoV, and lung samples of COVID-19 patients and healthy controls were retrieved from NCBI (64) (https://www.ncbi.nlm.nih.gov/) and ENA (65) (https://www.ebi.ac.uk/ena) (accession numbers GSE147507 (52), GSE56189, GSE148729 (58) and GSE153940 (66)). The raw FASTQ data of PBMC and BALF samples of COVID-19 patients and corresponding controls were downloaded from BIG Data Center (67) (https://bigd.big.ac.cn/) (accession number CRA002390) (59), and the raw FASTQ data for BALF healthy control samples were downloaded from NCBI (accession numbers SRR10571724, SRR10571730, and SRR10571732 under project PRJNA434133 (60)). The preprocessed single cell RNA-Seq data of BALF samples from 6 severe COVID-19 patients and 3 moderate COVID-19 patients were downloaded from NCBI with accession number GSE14−5926 (61). The preprocessed single cell RNA-Seq data of BALF sample from a healthy control was retrieved from NCBI (accession number GSM3660650 under project PRJNA526088 (62)). Detailed information about these public datasets is available in the with parameters "-b 30 -single -l 180 -s 20" for single-end sequencing data and with parameter "-b 30" for paired-end sequencing data. Expression levels were calculated and summarized as transcripts per million (TPM) on gene levels with Sleuth (72), and logFC was then calculated for each condition. The single cell RNA-Seq data were summarized across all cells to obtain "pseudo-bulk" samples. R packages EDASeq (73) and org.Hs.eg.db (74) were used to obtain gene length, and TPM was calculated with the "calculateTPM" function of R package scater (75). logFC was then calculated for each patient. No corresponding 24h mock control samples for H1299 cells, 36h mock control samples were used instead. *There are three replicates, but when the manuscript was in preparation only two of them are available for downloading.
The clean RNA-Seq data were also aligned to the virus genome with Bowtie 2 (76) (version 2.2.6) and the aligned BAM files were created, and the mapping rates to the virus genomes were obtained as well. SAMtools (77) (version 1.5) was then used for sorting and indexing the aligned BAM files. The "SAMtools depth" command was used to produce the number of aligned reads per site along the virus genome.
The heatmap in Figure 3I was made by pheatmap R package (78), "complete" clustering method was used for clustering the rows and "euclidean" distance was used to measure the cluster distance. The heatmap in Figure 4A was made by ComplexHeatmap R package (79). "ward.D" clustering method was used for clustering the rows and "euclidean" distance was used to measure the cluster distance. Our analysis shows that the infection efficiency of viruses can be both cell type dependent and virus dose dependent (Figure 1). MERS-CoV can efficiently infect MRC5 and Vero E6 cells. However, the infection efficiency is influenced strongly by MOI in the same type of cells. Cells infected with low MOI, say 0.1, have significantly lower mapping rates than those with high MOI, say 3 ( Figure 1). For SARS-CoV and SARS-CoV-2, the infection efficiency is influenced strongly by cell type. For SARS-CoV-2, there is efficient virus infection in A549-ACE2, Calu3, Caco2, and Vero E6 cells, but not in A549, H1299, or NHBE cells (Figure 1 and Supplementary Table S1). The mapping rates in A549, H1299, and NHBE cells are low even at high MOIs ( Figure 1 and Supplementary Table S1). Similar to SARS-CoV-2, the infection by SARS-CoV is also cell type dependent, Vero E6 cells and Calu3 cells show high mapping rates to SARS-CoV genome, but the mapping rates of SARS-CoV in MRC5 and H1299 cells are close to zero even at the high MOI of 3 ( Figure 1 and Supplementary Table  S1). Since "total RNA" (see Methods/Data collection) includes additional negative-strand templates of virus, the mapping rates are usually much higher than those that used the PolyA+ selection method in the same condition ( Figure 1, Supplementary Table S1).

Evidence for Multiple Entry Mechanisms for SARS-CoV-2 and SARS-CoV
To examine the detailed replication and transcription status of these viruses in the cells, we calculated the number of reads (depth) mapped to each site of the corresponding virus genome ( Figure 2). For better comparison, these read numbers were log 10 transformed. The replication and transcription of MERS-CoV, SARS-CoV-2 and SARS-CoV share an uneven pattern of expression along the genome, typically with a minimum depth in the first half of the viral genome, and the maximum towards the end (Figures 2A-E). Among the parts with very high levels, there are especially coding regions for structural proteins, including S, E, M, and N proteins, as well as the first coding regions with nsp1 and nsp2. Interestingly, there is an exception for BALF samples in COVID-19 patients, which show a more irregular, fluctuating behavior along the genome ( Figure 2B). The deviation from the cellular expression pattern is not surprising because BALF is not a well-organized tissue but a mixture of many components, some of which will probably digest viral RNA.
Interestingly, the mentioned uneven transcription pattern of efficient infections with SARS-CoV-2, SARS-CoV, and MERS-CoV, is also visible for inefficient infection with SARS-CoV-2 in A549, NHBE, and H1299 cells, and SARS-CoV in H1299 and MRC5 cells ( Figures 2C, D), although there the total mapping rates to their corresponding virus genomes are much lower (Figure 1).
To further elucidate the corresponding entry mechanisms for different types of cells, we examined the expression levels of those receptors and proteases that have already been described as facilitating target cell infection (Figure 3).
Our analysis shows that MERS-CoV can efficiently infect MRC5 and Vero E6 cells (Figure 1 and Figure 2E) that both express DPP4 ( Figure 3A), though compared to Vero E6 cells, MRC5 cells infected with MERS-CoV have higher expression levels of DPP4 ( Figure 3A), but lower mapping rates to the virus genome ( Figure 1). These observations show that higher expression levels of the receptor (DPP4) do not guarantee higher MERS-CoV infection efficiency in cells. This is also true for SARS-CoV-2 receptor ACE2, which is expressed three orders of magnitudes higher in A549-ACE2 cells than in Vero E6 cells ( Figure 3B), while both cells produce about the same amount of virus ( Figure 1).
Although SARS-CoV-2 can efficiently infect A549-ACE2 cells (Figure 1 and Figure 2), there is no expression of TMPRSS2 or TMPRSS4 ( Figures 3C, D), needed for the canonical cell-surface membrane fusion mechanism ( Figure 3J). However, there are considerable expression levels of CTSB and CTSL ( Figures 3E, F), which are involved in endocytosis ( Figure 3J).
In A549, H1299, and MRC5 cells, which do express small amounts of SARS-CoV-2 and SARS-CoV virus (Figure 1,  Figures 2C, D), there is no ACE2 expression at all ( Figure 3B). This could point to an alternative ACE2-independent entry mechanism for SARS-CoV-2 and SARS-CoV ( Figure 3J). Since there were already reports about alternative SARS-CoV-2 receptors such as BSG/CD147 and CD209 (30, 31), we examined their expressions in these cells as well (Figures 3G, H). For all cells, the expression of BSG is at the same level of 2-3 ( Figure 3G), and the expression of CD209 is very low. Certainly, CD209 and BSG alone cannot explain the differences in virus expression (Figure 1), nor can we exclude other low efficiency entry mechanisms. It could e.g. be that relatively inefficient alternative entry paths are often present but in some cells masked by more efficient entry via ACE2/TPMRSS.
To gain a comprehensive overview we clustered cells with respect to gene expression levels of coronavirus-associated receptors and factors ( Figure 3I), and summarized conceivable mechanisms accordingly ( Figure 3J). Since all cells show high expression levels of CTSB and CTSL, the major differences between these cells lie in the expression levels of ACE2, TMPRSS2 and TPMRSS4.
Cell-surface membrane fusion ( Figure 3J, 1a) might be mainly used in SARS-CoV-2 infection of Calu3, Caco2, and NHBE cells where there are low to moderate expression of ACE2 and moderate expression of TMPRSS2 and TMPRSS4. Endocytosis ( Figure 3J, 1b) might be mainly used in SARS-CoV-2 infection of A549-ACE2 cells where ACE2 is expressed at high levels but there is no expression of TMPRSS2 or TMPRSS4. An alternative ACE2-independent way ( Figure 3J, 1c) in absence of ACE2, TMPRSS2, or TMPRSS4 could

Strength of IFN/ISG Response Varies Between Cell Lines and Viruses, With Strong Response to SARS-CoV-2 in Relevant Cells
As a virus enters a cell, it may trigger an innate immune response, i.e. the cell may start expression of various types of innate immunity molecules at different strengths. There is currently an intense debate about which of these molecules, especially IFNs and ISGs, are expressed how strongly (52)(53)(54)(55)(56). We therefore focused in our analysis on innate immunity molecules such as IFNs and ISGs. To broaden the basis for conclusions, we analyzed, apart from cell lines, bulk RNA-Seq data of lung, PBMC, and BALF samples of COVID-19 patients, and single-cell RNA-Seq data of BALF samples from moderate and severe COVID-19 patients; for each type of patient data, we also included healthy controls. Gene expressions were compared quantitatively in terms of TPM (transcripts per million), as well as log fold changes (logFC) with respect to healthy controls (human samples) or mock-infected cultures (cell lines) (Supplementary Figure S1, S2). The heatmap and clustering dendrogram of the logFC of IFNs and ISGs in Figure 4A reveal broadly two groups of samples with fundamentally different expression of PRRs, downstream transcription factors, IFNs and ISGs.
The top cluster in Figure 4A Figure S8), which explains why there are almost no virus reads in these tissues.
One of the two lung samples (accession number: S A M N 1 4 5 6 3 3 8 7 ) h a s s l i g h t l y u p r e g u l a t e d I F N L 1 (Supplementary Figure S6), which had been ignored in the original publication (52), although the total mapping rates to virus genome are both 0.00% for these two lung samples. We then checked the detailed coverage along the virus genome. There were a small number of virus reads aligned to SARS-CoV-2 genome in this sample (Supplementary Figure S7). Different from the other lung sample that did not express ACE2, this lung sample expressed ACE2 at a considerable level (5.45 TPM, Supplementary Table S2). This result implies when the lung tissue chosen for sequencing are successfully infected by SARS-CoV-2, IFNs (at least IFNL1) can be upregulated.
Calu3 cells infected with SARS-CoV and SARS-CoV-2, and A549-ACE2 cells infected with SARS-CoV-2 at a high MOI of 2 have upregulated IFNB1, IFNL1, IFNL2 and IFNL3 ( Figures 4B-E), as well as other interferons ( Figure 4A), although not as significantly as IFNB1 and INFL1-3. A549, H1299, NHBE ( Figures 4B-E), and MRC5 cells (Supplementary Figure S3), which do not support efficient virus infection, show no upregulation of IFNs. Low expression levels of IFNs are also observed in Caco2 cells, which are efficiently infected with SARS-CoV and SARS-CoV-2. The same is true for A549-ACE2 cells infected with SARS-CoV-2 at low MOI of 0.2. In Vero E6 cells IFNL1 is upregulated as well in infected with SARS-CoV and SARS-CoV-2, but not with MERS-CoV ( Figure 4F). In BALF samples of moderate and severe COVID-19 patients, upregulation of IFNs was not as obvious as in Calu3 cells, but is still present in some patients ( Figure 4G). These observations demonstrate that the innate immune response depends in complex ways on cell type, viral dose, and virus.

DISCUSSION
One attractive potential anti-SARS-CoV-2 therapy is intervention in the cell entry mechanisms (21). However, the entry mechanisms of SARS-CoV-2 into human cells are partly unknown. During the last few months scientists have confirmed that SARS-CoV-2 and SARS-CoV both use human ACE2 as entry receptor, and human proteases like TMPRSS2 and TMPRSS4 (8,23,36), and lysosomal proteases like CTSB and CTSL (39) as entry activators. Since ACE2 is beneficial in cardiovascular diseases such as hypertension or heart failure (80), treatments targeting ACE2 could have a negative effect. Inhibitors of CTSL (81) or TMPRSS2 (23) are seen as potential treatment options for SARS-CoV and SARS-CoV-2. However, recently alternate coronavirus-associated receptors and factors including BSG/CD147 (30) and CD209 (31) have been proposed to facilitate virus invasion. Additionally, clinical data of SARS-CoV-2 infection have shown that SARS-CoV-2 can infect several organs where ACE2 expression could not be detected (33,34), urging us to explore other potential entry routes.
First, our analyses here have shown that even without expression of TMPRRS2 or TMPRSS4, high SARS-CoV-2 infection efficiency in cells is possible (Figures 1A, C) with considerable expression levels of CTSB and CTSL ( Figures 2E, F). This suggests receptor mediated endocytosis (24,25,39) as an alternative major entry mechanism. Given this TMPRSS-independent route, TMPRSS inhibitors will likely not provide complete protection. The studies designed to predict the tropism of SARS-CoV-2 by profiling the expression levels of ACE2 and TMPRSS2 across healthy tissues (82,83) may need to be reconsidered as well.
Second, the evidence presented in our study suggests further, possibly undiscovered entry mechanism for SARS-CoV-2 and SARS-CoV ( Figure 2). Although BSG/CD147 has been recently proposed as an alternate receptor (30), later experiments  reported there was no evidence supporting the role of BSG/ CD147 as a putative spike-binding receptor (84). The expression patterns of BSG/CD147 in different types of cells observed in our study could not explain the difference in virus loads observed in these cells either. CD209 and CD209L were recently reported as attachment factors to contribute to SARS-CoV-2 infection in human cells as well (31). However, CD209 expression in the cell lines included here is low. Another reasonable hypothesis could be that the inefficient ACE2-independent entry mechanism we observed could be macropinocytosis, one endocytic pathway that does not require receptors (85). Until now there is still no direct evidence for macropinocytosis involvement in SARS-CoV-2 and SARS-CoV entry mechanism. To confirm such an involvement, specific experiments are needed. Moreover, this ACE2independent entry mechanism, only enables inefficient infection by SARS-CoV and SARS-CoV-2 ( Figure 2) and therefore cannot be a major entry mechanism. Figure 3J summarizes the outcomes of our study with respect to entry mechanisms. The observations with the broad range of transcriptome data can only be explained if there are several entry routes. This is certainly a challenge to be reckoned with in the development of antiviral therapeutics (86).
Another attractive potential anti-SARS-CoV-2 point of attack is supporting the human innate immune system to kill the infected cells and, thus disrupt viral replication. Not surprisingly, research in this area is flourishing but sometimes generates conflicting results, especially on the involvement of type I and III IFNs and ISGs (52)(53)(54)(55)(56). The results of our analyses could help to dissolve the confusion on the involvement of IFNs and ISGs.
We found that immune responses in Calu3 cells infected with SARS-CoV and SARS-CoV-2 resemble those of BALF samples of moderate and severe COVID-19 patients, with elevated levels of type I and III IFNs, robust ISG induction as well as markedly elevated pro-inflammatory cytokines, in agreement with recent studies (53)(54)(55)(56). However this picture differs from the one reported by (52) with low levels of IFNs and moderate ISGs. This latter study was partially based on A549 cells and NHBE cells with nearly no ACE2 expression and very low mapping rate to the viral genome, and lung samples of two patients (both show 0.00% mapping rate to virus genome). Hence, given that there was no efficient virus infection in these cells, the low levels of IFNs and ISGs were to be expected. However, in one of the lung samples sequenced by (52) (accession number: SAMN14563387), we observed a slight upregulation of IFNL1 (Supplementary Figure S6), which was ignored in the original publication, together with considerable ACE2 expression (Supplementary Table S2) (5.45 TPM), and a few virus reads aligned to SARS-CoV-2 genome (Supplementary Figure S7). This result suggests that levels of IFNs are ISGs are associated with viral load and severity of virus infection.
We found low induction of IFNs and moderate expression of ISGs in PBMC samples and BALF samples of COVID-19 patients (Figure 4, Supplementary Figure S5). In these PBMC samples, there are no (0.00%) virus reads mapping to the SARS-CoV-2 genome. The failure to detect virus reads in these three PBMC samples can be explained by the absence of efficient entry routes (e.g. no expression of ACE2 in PBMC samples of healthy individuals, Supplementary Figure S8), or with the cell types being otherwise incompatible with viral replication. This observation is consistent with the studies on SARS-CoV (87)(88)(89) with abortive infections of macrophages, monocytes, and dendritic cells; moreover, replication of SARS-CoV in PBMC samples is also self-limiting. However, due to the limited number of PBMC, BALF and lung samples included in this study, and the lack of the information of sample collection time points, infection stage and infection severity of these COVID-19 patients, the assessment of IFNs and ISGs as well as the infection of SARS-CoV-2 in these samples may not be representative of host response against SASR-CoV-2. Future studies that include also other affected organs of more patients with different infection stages and severity are necessary for a better understanding of the immune responses.
Several unexpected observations need further investigations. First, A549-ACE2 and Caco2 cells are efficiently infected with low MOI of 0.2 and 0.3 respectively (Figure 1), but fail to upregulate INF expression ( Figures 4B-E). Their cellular immune responses are more similar to those of cells that cannot support efficient virus infection ( Figure 4A). These results suggest that in Caco2 and A549-ACE2 cells the invasion of SARS-CoV-2 or SARS-CoV at low MOI shuts down or fails to activate the innate immune system.
Based on the results observed above, multiple factors including disease severity, different organs, cell types and virus dose contribute to the variability in the innate immune responses. For a better characterization of the innate immune responses, a more comprehensive profiling is necessary, including patients with infections in different stages, different levels of severity, and different clinical outcomes of the infection. Further, a larger array of cell types should be profiled over time after infection with different virus doses. In this way we would be better able to understand the kinetics of IFNs and ISGs in response to SARS-CoV-2 infection. It also needs be noted that our research is based on transcriptomic data, so the expression levels of coronavirusassociated receptors and factors, IFNs and ISGs, are measured as mRNA levels. For further confirmation of our conclusion, more proteomic studies data are required.
In summary, our study has comparatively analyzed an extensive data collection from different cell types infected with SARS-CoV-2, SARS-CoV and MERS-CoV, and from COVID-19 patients. We have presented evidence for multiple SARS-CoV-2 entry mechanisms. We could also dissolve apparent conflicts on innate immune responses in SARS-CoV-2 infection (52)(53)(54)(55)(56), by drawing upon a larger set of cell types and infection severity. The results emphasize the complexity of interactions between host and SARS-CoV-2, offer new insights into pathogenesis of SARS-CoV-2, and can inform development of antiviral drugs.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

AUTHOR CONTRIBUTIONS
PH and YC conceived the research. DH, PH, and YC designed the analyses. YC and XX conducted the analyses. All authors contributed to the article and approved the submitted version.

FUNDING
This work was partially funded by grant 01Kl20185B (SECOVIT) of the German Federal Ministry of Education and Research.