key: cord-1055234-9xpbwf9c
authors: Li, Jingjiao; Guo, Mingquan; Tian, Xiaoxu; Wang, Xin; Yang, Xing; Wu, Ping; Liu, Chengrong; Xiao, Zixuan; Qu, Yafei; Yin, Yue; Wang, Chunxia; Zhang, Yucai; Zhu, Zhaoqin; Liu, Zhenshan; Peng, Chao; Zhu, Tongyu; Liang, Qiming
title: Virus-host interactome and proteomic survey of PBMCs from COVID-19 patients reveal potential virulence factors influencing SARS-CoV-2 pathogenesis
date: 2020-07-21
journal: Med (N Y)
DOI: 10.1016/j.medj.2020.07.002
sha: c9947018541b5ab2410bed8c5221aeb7467e434e
doc_id: 1055234
cord_uid: 9xpbwf9c

Summary Background The ongoing coronavirus disease (COVID-19) pandemic caused by Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) is a global public health concern due to relatively easy person-to-person transmission and the current lack of effective antiviral therapy. However, the exact molecular mechanisms of SARS-CoV-2 pathogenesis remain largely unknown. Methods Genome wide screening was used to establish intra-viral and viral-host interactomes. Quantitative proteomics was used to investigate peripheral blood mononuclear cell (PBMC) proteome signature in COVID-19. Findings We elucidated 286 host proteins targeted by SARS-CoV-2 and more than 350 host proteins that are significantly perturbed in COVID-19 derived PBMCs. This signature in severe COVID-19 PBMCs reveals significant upregulation of cellular proteins related to neutrophil activation and blood coagulation, as well as downregulation of proteins mediating T cell receptor signaling. From the interactome, we further identified that non-structural protein 10 interacts with NF-kappa-B-repressing factor (NKRF) to facilitate interleukin-8 (IL-8) induction, which potentially contributes to IL-8-mediated chemotaxis of neutrophils and the overexuberant host inflammatory response observed in COVID-19 patients. Conclusions Our study not only presents a systematic examination of SARS-CoV-2-induced perturbation of host targets and cellular networks but also reveals insights into the mechanisms by which SARS-CoV-2 triggers cytokine storms, representing a powerful resource in the pursuit for therapeutic intervention.

In December 2019, Severe Acute Respiratory Syndrome Coronavirus-2 [SARS-CoV-2, previously known by the provisional name "2019 novel coronavirus" (2019-nCoV) given by the World Health Organization (WHO)], was first isolated from a cluster of patients with pneumonia. This novel strain of coronavirus has been identified as the causative agent for the epidemic of acute respiratory distress syndrome [also referred to as coronavirus disease 2019 (COVID- 19) ] in Wuhan, China. [1] [2] [3] Since then, the COVID-19 outbreak has rapidly swept all the provinces of China and has disseminated globally to 216 countries and territories through traveling. 4 On 31 January 2020, the WHO declared SARS-CoV-2 a Public Health Emergency of International Concern. 5 As of 25 June 2020, the current pandemic of SARS-CoV-2 has caused over 9 million confirmed infection cases including 479,133 deaths in 216 countries and territories, calling for extended efforts to understand the molecular drivers of its immune evasion and pathogenesis.

COVID-19 is a potential zoonotic disease with a low to moderate mortality rate (2-5%) and is primarily transmitted through droplet and direct contact with infected persons or incubation carriers. 6 COVID-19 patients develop dynamic innate and adaptive immune responses to SARS-CoV-2. 7-9 SARS-CoV-2 has a preferential tropism for airway epithelial cells mediated through the accessibility of the ACE2 receptor. [10] [11] [12] [13] The common clinical manifestations of fever and cough in conjunction with laboratory results of progressively increasing neutrophil counts and leukocytopenia in serum samples taken from COVID-19 patients at various stages of the illness indicates uncontrolled neutrophil extravasation and activation. [14] [15] [16] Neutrophils are the most abundant type of white blood cells in human blood and serve as the first responders in the line of defense against acute infection. 17 Chemical signals such as interleukin-8 (IL-8) and IL-6 stimulate neutrophil chemotaxis and recruitment to the site of infection where they congregate and cause the release of more pro-inflammatory mediators that in turn activate additional inflammatory signals, leading to cytokine storms that may progress to vital organ failure. 18 Indeed, a large number of macrophages and neutrophils were found in alveolar cavities in COVID-19 fatal patients, 19 suggesting their roles for cytokine storms. While uncontrolled neutrophil activity and aberrant cytokine response contribute to SARS-CoV-2 pathogenesis in COVID-19, how SARS-CoV-2 causes these immunological changes, especially the enhancement of neutrophil chemotaxis and activation remains largely undetermined.

SARS-CoV-2 is a single-stranded positive-sense RNA virus and belongs to the β-coronavirus family consisting of many animal species and human pathogens including SARS-CoV and Middle East Respiratory Syndrome (MERS)-CoV. The ~30 kb non-segmented genome of SARS-CoV-2 encodes 14 major open reading frames (ORFs), which are further processed into 4 structural proteins [glycoprotein spike (S), membrane (M), envelope (E), and nucleocapsid (N)], 16 non-structural proteins (nsp1 to nsp16), and at least 8 accessory proteins (ORF3a, ORF6, ORF7a, ORF7b, ORF8, ORF9b, ORF9c, and ORF10). 20, 21 Similar to SARS-CoV, spike protein of SARS-CoV-2 directly interacts with human ACE2 expressed on human airway epithelia to facilitate virus entry. [10] [11] [12] [13] However, host targets for other viral proteins remain largely unknown, which limits our understanding of SARS-CoV-2 immune evasion and pathogenesis, further impeding anti-viral drug development and the design of therapeutic approaches.

Limited by their relatively small genomes, viral proteins form different combinations of protein complexes to fulfill a variety of critical functions during the viral life cycle, making the investigation of intra-viral protein-protein interactions (PPIs) a priority. To systematically characterize the intra-viral PPI networks, both genome-wide yeast-two hybrid (Y2H) screens and co-immunoprecipitations (co-IP) experiments with SARS-CoV-2 genes were performed. All the predicted viral genes were synthesized based on the SARS-CoV-2 Wuhan-Hu-1 isolate (GenBank: MN908947). The entire library was cloned into the Y2H DNA-binding domain (pACT2) and activation domain (pGBTK7) vectors and then transformed into the Y2HGold yeast strain (Clontech). Nsp3 was truncated into two smaller fragments in this experiment since large genes tend to yield fewer interactions in the Y2H assay. All the possible pair-wise combinations were examined independently for the activation of the reporter genes HIS3 and ADE2 in triplicates. We captured 19 PPIs from the Y2H screens, including nsp15-nsp2, ORF7b-nsp1, N-nsp4, ORF3a-nsp4, ORF7b-nsp7, nsp7-nsp8, nsp12-nsp8, nsp7-nsp9, nsp14-nsp10, nsp16-nsp10, ORF6-nsp14, E-E, N-E, ORF3a-E, ORF9b-E, N-N, nsp7-ORF7a, nsp8-ORF8, and ORF9b-ORF9b ( Figure S1 ). Most intra-viral PPIs characterized by Y2H were confirmed by co-IP ( Figure S2A ). The Co-IP assay resulted in the identification of an additional 39 PPIs (Figures S1 and S2B) . Some viral proteins, such as M, N, E, nsp2, nsp5, nsp8, ORF6, ORF7a, ORF7b, ORF9b, and ORF10 displayed selfassociation ( Figure 1A) , suggesting their dimer or oligomer forms may carry out important functions during infection. In total, we characterized 58 distinct intra-viral PPIs amongst 28 SARS-CoV-2 genes that are potentially involved in virus replication and can be expected to be important to immune evasion and viral pathogenesis ( Figure 1A) .

Two independent intra-viral PPI networks of closely related SARS-CoV have been previously established with Y2H or mammalian two-hybrid systems, identifying 72 and 40 intra-viral PPIs, respectively. 22, 23 Due to the amino acid sequence similarity, 20 intra-viral PPIs are overlapped between SARS-CoV-2 and SARS-CoV, including nsp2-nsp15, nsp7-nsp8, nsp7-nsp9, nsp8-nsp12, nsp8-nsp13, nsp10-nsp14, nsp10-nsp16, nsp10-nsp16, nsp16-N, nsp16-ORF7a, M-N, M-ORF6, E-ORF9b, and ORF6-ORF7a, as well as selfassociations of nsp2, nsp5, nsp8, E, N, and ORF9b, suggesting these intra-viral PPIs play critical roles for maintaining fundamental functions in SARS coronavirus family.

Among the 4 classical structural proteins of SARS-CoV-2, 3 PPIs were detected ( Figure 1B) ; these interactions of structural proteins may be critical in the organization of viral particles during viral replication and assembly. We also found 8 PPIs between structural proteins and accessory proteins (ORF3a, ORF6, ORF7a, ORF7b, ORF9b, and ORF10) and 6 PPIs between structural proteins and non-structural proteins (nsp2, nsp4, nsp5, nsp8, and nsp16) (Figure 1B) , which are potentially involved in viral structural protein processing, modification, and trafficking, supporting their importance for viral particle assembly and egress.

While structural proteins are components for the capsid and the envelope of the virus, several nonstructural proteins (nsps), such as nsp3, nsp5, nsp12, nsp14, and nsp16 have associated enzymatic activities essential for coronavirus RNA replication and viral protein translation. In addition, nsps can function as virulence factors that inhibit the host's immune system defenses. We speculate that the relatively large number of nsp interactions detected (12) might imply a crucial role of nsps in the formation of multimeric complexes during the viral life cycle ( Figure 1C) . Indeed, nsp7 and nsp8 support the RNA-dependent RNA polymerase function of nsp12. [24] [25] [26] Similar to SARS-CoV, 27,28 SARS-CoV-2 nsp10 may facilitate the methyltransferase activities of nsp14 and nsp16 by nsp10-nsp14 and nsp10-nsp16 associations, respectively. While the functional importance of most intra-viral PPIs remains less understood, our network provides a powerful resource to dissect the role of these viral proteins in mediating the synthesis, assembly, and processing of viral proteins for SARS-CoV-2. Disruption of these intra-viral PPIs by repurposing small molecule or peptide drugs that target specific epitopes can be a treatment strategy against COVID- 19. 29 Genome-wide proteomic screen for SARS-CoV-2 viral-host protein-protein interactions SARS-CoV-2 infection can lead to imbalanced host immune responses such as diminished type I and type III interferon (IFN) signaling as well as strong chemotactic and inflammatory response in infected cell lines, animal models, and COVID-19 patients, 7, 8, [30] [31] [32] [33] suggesting this virus must encode components to modulate cellular signaling pathways and shape its unique immune signature. To better understand the roles of these viral proteins during the virus life cycle, it is pivotal to characterize the host targets of SARS-CoV-2 in human cells. We overexpressed the plasmids encoding each SARS-CoV-2 encoded gene with N-terminal 3xFlag-epitope in HEK293 cells and purified each SARS-CoV-2 protein complex by affinity purification. Copurified cellular proteins were subsequently analyzed by liquid chromatography and tandem mass spectrometry (AP-LC-MS). With statistical analysis by MiST scoring algorithm (MiST score > 0.8, spectral count > 8, Fold change > 5), 34 286 cellular proteins were identified to interact with SARS-CoV-2 proteins, resulting in 295 high-confidence interactions (Figures 2, S3A and Table S1 ). Host cellular pathways that are most likely affected by SARS-CoV-2 viral proteins are ATP biosynthesis and metabolic processes (M), mRNA transport (N), MDA5 and RIG-I RNA sensing signaling (nsp1), nucleotide-excision repair (nsp4), protein methylation and alkylation (nsp5), translation initiation (nsp9), cellular amino acid metabolic process (nsp10), ROS metabolic process (nsp14), Golgi to plasma membrane transport (nsp16), endoplasmic reticulum stress (ORF3a), and mRNA transport (ORF6) (Figures S3B-C) , which is consistent with the proteome analysis of SARS-CoV-2 infected cells. 30 It is important to note that our AP-LC-MS analysis shares ~16% overlap (45 targets) with the SARS-CoV-2-host protein interactions recently identified by the Krogan Laboratory using 2xStrep-tagged viral proteins isolated from HEK293T/17 cells. 35 Some cellular targets characterized by both interactomes include G3BP1/G3BP2 (anti-viral stress granule proteins), 36 RAE1 (an interferon-inducible mRNA nuclear export protein), 37 TBK1 (essential kinase for type I IFN pathway), 38 TRIM59 (ubiquitin E3 involved in antiviral response), 39 and SigmaR1 (involved in lipid remodeling and ER stress response). 35 The identification of common SARS-CoV-2 protein and human protein interactions suggest these cellular processes are vital molecular targets of SARS-CoV-2 infection. Indeed, PB28, an agonist of Sigma receptor, has potent anti-SARS-CoV-2 activity in cell culture models. 35 Notably, our AP-LC-MS analysis uncovered additional host targets pertinent to inflammatory and innate immune responses (Figure 2 ) and may help contextualize the unique immune signature identified in SARS-CoV-2-infected cells. 7, 9 Key viral-cellular PPIs that may influence the pathogenesis of SARS-CoV-2 include interactions between NKRF and nsp10, TBK1 and ORF6, TANK and ORF6, TRAF2 and ORF6, and KIT and ORF3a. NKRF (NF-kappa-B-repressing factor) is a transcriptional repressor and controls the induction of IL-8 to stimulate neutrophil chemotaxis and recruitment to the site of infection. 40, 41 Nsp10 may regulate IL-8 levels by targeting NKRF and shapes the unique immune signature in COVID-19. TBK1 plays critical roles for type I IFN and NF-κB pathways upon viral infection: it phosphorylates the key transcriptional factors, IRF3 and IRF7, to activate IFN promoters; 38 it also forms a ternary complex with TANK and TRAF2 to facilitate NF-κB activation. 42 ORF6 may interact with the TBK1-TANK-TRAF2 ternary complex to modulate both type I IFN and NF-κB pathways upon infection, facilitating SARS-CoV-2 to escape from host innate immune surveillance. The tyrosine-protein kinase KIT has been shown to facilitate the activation of STATs, 43 such as STAT1, STAT3, and STAT5, which are involved in signaling transduction upon IFN or IL-6 stimulation. ORF3a might regulate IFN or IL-6 signaling pathways by association with KIT.

Several viral proteins that form intra-viral PPIs ( Figure 1A ) are associated with host ubiquitin pathway components, such as interaction partners ORF7a-UBE2J1 (E2 ligase), ORF3a-TRIM59 (E3 ligase), and ORF6-RNF20, RNF40, RFWD3, and MYCBP2 (E3 ligase), suggesting the potential modulation of the host ubiquitin system by these viral protein complexes (Figure 2) . Our observations provide an overview of host proteins and signaling pathways that are modulated by SARS-CoV-2 viral proteins and these PPIs might serve as potential druggable targets towards the development of new antiviral drugs or repurposing of existing ones.

To elucidate the mechanism and effect of SARS-CoV-2 infection in patients, we determined the global proteome profile changes in PBMCs samples taken from healthy donors compared with PBMCs collected from COVID-19 patients with mild or severe symptoms. The PBMC samples were obtained from the Shanghai Public Health Clinical Center and include 6 healthy donors, 22 COVID-19 patients with mild symptoms (COVID-19-MS), and 13 COVID-19 patients with severe symptoms (COVID-19-SS) ( Table S2) . PBMCs were harvested from the leftover blood samples and all the patients were confirmed negative for influenza A virus. Following the Guideline of Novel Coronavirus Laboratory Safety (Edition 2), PBMCs were collected, lysed in denaturing buffer (STAR Methods) at the Shanghai Public Health Clinical Center, and sent for proteome analysis at the National Facility for Protein Science in Shanghai.

A total of 6,739 proteins were identified in our study, of which 4,274 proteins could be quantified in all PMBCs samples. Proteins with less than three quantitation values in each category did not pass the filtering threshold and were excluded from the quantification and statistical analysis. 5,521 proteins were quantified according to criteria and included in global proteomic profiling to reveal a list of 726 proteins found to be differentially expressed [|log2 (fold change)| ≥ 0.58, adjusted p value ≤0.01, n≥3] in the PBMCs of COVID-19-MS compared with that of healthy donors: 333 proteins showed a statistically significant increase while 393 proteins showed a reduction in expression when comparing COVID-19-MS with control ( Figures 3A, S4A -B and Table S3 ). Consistent with animal infection model, 7,9 while some interferon-stimulated genes (ISGs), such as OAS3, MX2, and IFIT5, were upregulated, there was no obvious type I IFN response or overall strong induction of ISGs, suggesting SARS-CoV-2 encoded components may evade cellular IFN response to facilitate virus replication (Figures 3A-B and S4C) . Functional analysis by Cluster Profiler revealed that the underexpressed proteins were enriched in specific cellular biological processes, including negative regulation of megakaryocyte differentiation, blood coagulation, hemostasis, and negative regulation of immune system processes (Figures 3C and S4C) . Analysis of the global proteomic changes reveals significantly reduced levels of PF4, GP1BA, SERPINA5, HBB, HBD, C1QBP, HRG, and RAC1 in the COVID-19-MS PBMC cohort ( Figure 3B ). PF4 (CXCL4) could neutralize the anticoagulant effect during platelet aggregation and serve as a chemotactic signal for neutrophils and monocytes. 44, 45 GP1BA is a surface membrane protein on platelets and is involved in the formation of platelet plugs. 46 SERPINA5, a serine protease inhibitor, acts as a procoagulant and inflammatory factor and plays hemostatic roles in blood plasma. 47, 48 C1QBP, a mitochondrial complement component binding protein, is involved in the inflammatory response to infections. 49-51 RAC1 is a plasma membrane-associated small GTPase that plays a general role in the platelet procoagulant response. 52 RAC1 also directly interacts with SigmaR1, a potential therapeutic target for COVID-19. 35 Overall, our functional analysis reveals the complexity of host responses to SARS-CoV-2 infection in mildly symptomatic COVID-19 patients, highlighting the abnormal regulations in blood coagulation and hemostasis.

Next, we compared the proteome change between PBMCs isolated from COVID-19-SS and COVID-19-MS patients. Amongst 6,235 quantified proteins (n≥3), 366 proteins were found to be significantly differentially expressed [|log2 (fold change)| ≥ 0.58, adjusted p value ≤0.01] between COVID-19-SS and COVID-19-MS's PMBCs. 209 of these proteins were upregulated and 157 of these proteins were downregulated in COVID-19-SS as compared to COVID-19-MS (Figures 3D-E, S4D -E, and Table S3 ). Functional analysis by Cluster Profiler revealed that these upregulated genes are highly enriched in pathways associated with neutrophil activation, neutrophil degranulation, gas transport, bicarbonate transport, myeloid cell differentiation, and blood coagulation (Figures 3F and S4F) . The expression levels of neutrophil surface receptors (CXCR1, Siglec5, and CD177) and neutrophil granule contents (DEFA1) were dramatically upregulated in COVID-19-SS ( Figures 3D-E) , which is consistent with the elevated neutrophil population in COVID-19 patients. 14, 16 Consistent with RNA sequencing results form COVID-19 patients, 7, 9, 32 these differentially expressed proteins support a model by which SARS-CoV-2 triggers cytokine release, such as IL-6 and IL-8, resulting in uncontrolled infiltration and activation of neutrophils to the sites of infection followed by eventual progression to cytokine storm and respiratory distress. Indeed, laboratory results reveal increased neutrophils counts to ~90% of the white blood cell count as COVID-19 progresses to critical stages. [14] [15] [16] The elevation of IL-8 and IL-6 but not IL-17 in PMBCs collected from severe symptomatic COVID-19 patients indicates critical roles for IL-8/IL-6 in the severe inflammatory response induced by SARS-CoV-2 infection (Figures 3G-H and S4H) . In addition, a cluster of genes related to blood coagulation was dramatically upregulated in COVID-19-SS as compared to COVID-19-MS ( Figure S4G) , supporting the hypothesis that the abnormal blood clotting contributes to death in some severe COVID-19 patients. 53 As compared to COVID-19-MS, the downregulated proteins of COVID-19-SS are highly enriched in T cell activation, T cell receptor (TCR) signaling, viral transcription, viral gene expression, translation initiation, and lymphocyte co-stimulation (Figures 3F and S4F) . The expression of TCR subunits (CD3e, CD3g, and CD247), T cell surface molecules (CD4, CD8a, CD8b, CD2, CD5, and CD7), MHC class II molecules (HLA-DRA), TCR signaling kinases and adaptors (LCK, and FYN, SKAP1, CARD11, CAMK4, and BTN3A1), as well as T cell differentiation and survival (BCL11B) were dramatically diminished in COVID-19-SS compared to COVID-19-MS (Figures 3D-E) , suggesting a failure of T cell activation and function. These results demonstrate that severe COVID-19 patients undergo a functional decline in adaptive immunity, which correlates with their lower T cell population. 14-16

From the quantitative proteome analysis of PBMCs, it is clear that SARS-CoV-2 triggers IL-8 expression and increases neutrophil count which may mediate the manifestation of acute respiratory distress syndrome in COVID-19 patients. Therefore, we hypothesized that SARS-CoV-2 virulence factors may potentiate or regulate the induction of IL-8 during infection, leading to the uncontrolled infiltration and activation of neutrophils. From our virus-host interactome, we discovered that NKRF interacts with SARS-CoV-2 nsp10 (Figures 2 and S3A) . Co-IP further confirmed that nsp10 binds to endogenous NKRF in cells ( Figure 3I) . Notably, nsp10 did not bind to NKRF in Y2H assay (Figure S4I) , suggesting this interaction is either not direct a PPI, or additional post-translational modifications are required for this interaction. Individually expressed nsp10 facilitated IL-8 induction in lung epithelial A549 cells, while individually expressed nsp12, nsp13, or nsp15 had little or no effect (Figure 3J) , suggesting that nsp10 specifically promotes IL-8 induction. The levels of IL-1β and TNFα were not significantly affected by nsp10 (Figures S4J-K) , indicating nsp10 does not affect NK-κB signaling transduction overall. Since NKRF is a transcriptional repressor that inhibits IL-8 induction by competing with NF-κB for IL-8 promoter binding, 40, 41 we hypothesize nsp10 targets NKRF to modulate IL-8 expression. With specific siRNA transfection, the NKRF level was diminished in A549 cells (Figures 3K and  S4L) , results in elevated induction of IL-8 ( Figure 3L) , confirming the critical role of NKRF on IL-8 production. Importantly, nsp10 increased IL-8 expression level in cells treated with scrambled siRNA (Figure 3L , red bar), while it has almost no effect on IL-8 in cells treated with NKRF specific siRNAs, as compared to vector control ( Figure 3L , blue bar), indicating that nsp10 specifically targets NKRF to regulate IL-8 induction.

Similar to other coronaviruses, SARS-CoV-2 has undergone rapid adaptive changes against evolutionary pressures since its isolation. 16, 54 Phylogenetic and genetic analysis 55 of 160 public SARS-CoV-2 sequences reveals 137 amino acid substitutions in 23 viral proteins with two mutations (A61G and F68L) occurring on nsp10 with very low variation rate during evolution (1/160 for each site) (Figures S4M-N and  Table S4 ), which indicates nsp10 has remained relatively stable during evolution. Accordingly, these two nsp10 mutations did not affect its interaction with NKRF in cells ( Figure S4O) . Furthermore, although there are four amino acid variations between nsp10 SARS-CoV and nsp10 SARS-CoV-2 (Figure S4P) , nsp10 SARS-CoV remains tightly associated with NKRF ( Figure S4Q ), indicating these four sites are not essential for nsp10-NKRF binding. Therefore, our data suggest that nsp10 is a potential stable drug target to curb SARS-CoV-2-mediated IL-8 activation and subsequent recruitment of neutrophils.

Establishing viral-viral and virus-host PPI networks are of high priority to understand host response to SARS-CoV-2 infection and guide the design the potential drug candidates for COVID-19. So far, limited data has been obtained describing these PPI networks, preventing our understanding of virus immune evasion and pathogenesis. Using systematic approaches, we characterized 58 intra-viral PPIs and 295 virus-host PPIs for SARS-CoV-2. From our observation, most intra-viral PPIs identified by Y2H were confirmed by co-IP and more PPIs were detected by co-IP than Y2H in yeast. This is common for most viruses because cellular proteins or post-translational modifications might be required for the intra-viral complex formation. Due to the sequence similarity, 20 of 58 SARS-CoV-2 intra-viral PPIs are overlapped with SARS-CoV, 22 ,23 such as nsp7-nsp8, nsp8-nsp12, nsp10-nsp14, nsp10-nsp16, and nsp14-nsp16, indicating the general roles of these PPIs for coronavirus replication. The additional 38 intra-viral PPIs might provide evidence for the difference between SARS-CoV and SARS-CoV-2, not only for virus replication efficiency but also for the unique modulation of cellular signaling pathways.

An independent SARS-CoV-2-host interactome was recently generated by the Krogan group in HEK293T/17 cells with 2xStrep tag purification. 35 Although different tag purifications from different cell lines could generate variability on the binding partners, our SARS-CoV-2-host interactome shares 45 key host targets with the interactome from the Krogan group, such as TBK1, G3BP1, G3BP2, RAE1, TRIM59, and SigmaR1, which are involved in host innate immune responses to virus infection. Since lack the specific antibodies for viral protein immunoprecipitation, detection of endogenous virus-host PPIs in the context of virus remains currently unavailable. Some signal peptides and viral co-factors might also affect expression and localization of viral protein in cells, which potentially changed the binding partners of viral proteins in physiological conditions. Some true positive PPIs may have been discarded during analysis since we omitted ribosomal proteins and proteins bound to more than three viral targets during MiST analysis. With MiST statistical analysis, 34 we obtained 295 additional high-confident PPIs using parameters with MiST score > 0.8, Average peptide > 8 from two independent biological replicates, and Fold change > 5 as compared to GFP control. While nsp10-NKRF interaction was not captured by the Krogan group, 35 we have confirmed that nsp10 is associated with endogenous NKRF in HEK293T cells, suggesting our dataset could serve as a valuable additional resource for host targets that are potentially involved in SARS-CoV-2 immune evasion and pathogenesis. In addition, we evaluated the effect of nsp10 on IL-8 induction, highlighting the potential role of nsp10 in SARS-CoV-2 pathogenesis. Further functional analysis of nsp10 in the context of SARS-CoV-2 infection is in progress in our lab. Therefore, a combination list of intra-viral and SARS-CoV-2-host PPI networks could serve as a powerful resource to dissect the molecular mechanism for each SARS-CoV-2 viral protein, outlining the unique feature of SARS-CoV-2-mediated pathogenesis at the molecular level.

Although recent RNA sequencing analysis of infected cell lines, infected animal models and patient samples have improved our understanding of host response to SARS-CoV-2 infection, 7, 8, [30] [31] [32] [33] there is still lack of direct evidence at the protein level. Our quantitative proteomic analysis using PBMC samples derived from both mild and severe COVID-19 patients begins to address this gap in knowledge. Since our viral-host interactome analysis indicates SARS-CoV-2 viral proteins may modulate and evade host innate immune response by interacting with key kinases and adaptors of immune signaling pathway, it is not surprising that COVID-19 patients fail to mount robust anti-viral IFN responses. From our quantitative proteomic analysis, neutrophilia, leukocytopenia, and hypercoagulation were observed in PBMCs of severe COVID-19 patients as compared to mild cases, supporting the hypothesis that imbalanced inflammatory cytokine production, such as IL-8 and IL-6 dysregulation, contributes to the overactivation of neutrophils and regulation of the blood clotting cascade. 16, 56 These molecular signatures could serve as diagnostic markers to prognosticate the severity of SARS-CoV-2 infection. Further proteomic analysis with sorted cells will specifically address which cell types may play critical roles in shaping these clinical signatures in the future.

Overall, our study integrated virus-host interaction network and quantitative proteome approaches to identify host targets of SARS-CoV-2 viral proteins, demonstrating that SARS-CoV-2 evolved multiple mechanisms to modulate fundamental host cellular processes and evade the host immune system. To our knowledge, this is the first study to characterize the PBMC proteome profile in both mild and severe cases of COVID-19 infection, and contextualizes the overall clinical features of COVID-19 patients at the protein level. These systematic analyses not only generate groundwork for future SARS-CoV-2 research on the nature of this virus but also provide valuable targets for specific drug development or repurposing. With these approaches, we identified nsp10 as the potential virulence factor for SARS-CoV-2 through the interaction with NKRF to facilitate IL-8 induction and potentially lead to uncontrolled infiltration and activation of neutrophils (Figure 4) . Specific agents targeting IL-8/IL-6 or the downstream JAK kinases could be repurposed as suitable drug candidates for clinical evaluation. These small molecule inhibitors (i.e. IL-6 or IL-8 receptor antagonists) or agents blocking virulence factor PPIs would have enormous value as anti-virulence therapies in the treatment of SARS-CoV-2 infections (Figure 4) . Follow-up work is required to tease out the detailed mechanisms involved and translate the therapeutic relevance of these findings in targeting the cytokine storm in COVID-19.

The experiment design for viral-host PPIs in this study may not provide the exact physiological environment in terms of protein processing, accessory proteins etc. and hence additional validation in a more physiologically relevant system, such as during SARS-CoV-2 infection, would be needed. Our study provides a preliminary picture of the PPI network, but further validation is needed to confirm these findings and explore their functional significance. 

The authors declare no competing interests. Figure  S1 . Different types of viral proteins are labelled with the indicated colors. Self-associated viral proteins are circled in blue. (B) Interactions between SARS-CoV-2 structural proteins and other viral proteins. (C) Interaction network amongst SARS-CoV-2 non-structural proteins. Also see Figures S1 and S2.

Network representation of the high-confidence SARS-CoV-2-host interactome in HEK293 cells. SARS-CoV-2 bait proteins and interacting host protein are labelled with red squares and circles, respectively. Subcellular localization of host proteins is labelled with indicated colors. Proteins with known functions in inflammation and immune responses or ubiquitination pathway are circled in purple or red, respectively. Also see Figure S3 and Table S1 . 

Further information and requests for resources and reagents should be directed to and will be fulfilled by Lead Contact Qiming Liang (liangqiming@shsmu.edu.cn).

All unique/stable reagents generated in this study are available from the Lead Contact with completed Materials Transfer Agreement.

The datasets generated during this study were deposited to the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the iProX partner repository 59 with the dataset identifier IPX0002285000. Additional Supplemental Items are available at http://dx.doi.org/10.17632/kjs2wnncjs.1 from Mendeley Data.

This study was approved by the Ethics Committee of Shanghai Public Health Clinical Center (#YJ-2020-S052-02). The proteomic analysis of PBMC samples was performed on existing samples collected during standard diagnostic tests, posing no extra burden to patients.

HEK293 (ATCC, #CCL-1573), HEK293T (ATCC, #CRL-11268), and A549 cells (ATCC, #CCL-185) were maintained in Dulbecco's modified Eagle's medium (DMEM; Gibco-BRL) containing 4 mM glutamine and 10% FBS. Vero cells were cultured in DMEM with 5% FBS and 4 mM glutamine. Transient transfections were performed with Lipofectamine 3000 (Thermo Fisher Scientific, #3000015). SARS-CoV-2 genes were synthesized as DNA fragments (Sangon Biotech) and these expression constructs were amplified by PCR and cloned into pLVX-3xFlag-MCS-P2A-tagRFP (puro), pEF-MCS-3xHA, pGBTK7, or pACT2 vectors. Mutations of nsp10 were generated using the QuikChange site-directed mutagenesis kit (Stratagene). All constructs were sequenced using an ABI PRISM 377 automatic DNA sequencer to verify 100% correspondence with the original sequence.

HEK293 were transfected with 3xFlag-tagged plasmids encoding each SARS-CoV-2 protein by Lipofectamine 3000 (Thermo Fisher Scientific, #L3000015). 3x10 7 HEK293 cells were harvested in 10 ml lysis buffer (50 mM Tris-HCl [pH 7.5], 150 mM NaCl, 0.5% Nonidet P40, 10% glycerol, phosphatase inhibitors, and protease inhibitors) at 72 h post-transfection, followed by centrifugation and filtration (0.45 µm) to remove debris. Cell lysates were precleared with 100 µl protein A/G resin and then incubated with 25 µl anti-Flag M2 resin (Sigma, #F2426) for 12 h at 4°C on a rotator. The resin was washed 5 times and transferred to a spin column with 40 µl of 3 X Flag peptide (Sigma, #F4799) for 1 h at 4°C on a rotator. The purified complexes were sent for mass spectrometry analysis for protein-protein interaction networks at the National Facility for Protein Science in Shanghai.

The samples (purified complexes) were precipitated by TCA and resolved in 8 M urea, and then treated with 5 mM TCEP and 10 mM IAA to reduce the disulfide bonds and alkylate the resulting thiol groups, sequentially. The mixture was digested for 16 h at 37 °C by trypsin at an enzyme-to-substrate ratio of 1:50 (wt/wt). The trypsin-digested peptides were desalted by monospin C18 and loaded on a home-made 40 cmlong pulled-tip analytical column (75 µm ID x 360µm OD, ReproSil-Pur C18-AQ 1.9 µm resin, Dr. Maisch GmbH), connected to a Thermo Easy-nLC1000 HPLC system. The samples were analyzed with a 120 min-HPLC gradient from 0 to 100% of buffer B (0.1% formic acid in 80% acetonitrile) at 300 nL/min：0-1 min, 5%-10% B;1-96 min, 10-40% B; 97-104min, 40-60% B; 104-105min, 60-100% B; 105-120min, 100% B. The eluted peptides were ionized and directly introduced into a Q-Exactive mass spectrometer using a nano-spray source with a distal 2.0-kV spray voltage. Survey full-scan MS spectra (from m/z 300-1,800) was acquired in the Orbitrap analyzer with resolution r = 70,000 at m/z 200. One acquisition cycle included one full-scan MS spectrum followed by top 20 MS/MS events, sequentially generated on the first to the twentieth most intense ions selected from the full MS spectrum at a 28% normalized collision energy. The acquired MS/MS data were analyzed against a UniProtKB Human database (database released on Sept. 30, 2018) containing SARS-CoV-2 viral proteins using Maxquant V1.6.10.43. To accurately estimate peptide probabilities and false discovery rates (FDR), we used a decoy database containing the reversed sequences of all the proteins appended to the target database. FDR was set at 0.01. Mass tolerance for precursor ions was set at 20 ppm. Trypsin was defined as a cleavage enzyme and the maximal number of missed cleavage sites was set at 3. Carbamidomethylation (+57.02146) of cysteine was considered as a static modification. Methionine oxidation was set as variable modifications. For AP-MS analysis, proteins with spectral counts no more than 8 and fold change no more than 5 were excluded as background. MiST algorithm was used to calculate the statistical significance of the dataset. 34 To screening the high-confidence interactions, the MiST score setting at 0.8 was used to filter the interacting host targets. The ribosome proteins were also removed from the list due to their high background for protein purification. Finally, the interactome network was generated by Cytoscape software (https://cytoscape.org/). The information of 35 patients with confirmed SARS-CoV-2 infection (13 COVID-19 patients with severe symptoms and 22 COVID-19 patients with mild symptoms) and 6 health donors as control are listed in Table S2 , including age (47-73), sex, and symptoms. This study followed the Guideline of Novel Coronavirus Laboratory Safety (Edition 2) and was approved by the Shanghai Public Health Clinical Center (#YJ-2020-S052-02). The peripheral venous whole blood samples were collected from all patients and PBMCs were prepared using density-gradient centrifugation followed by washing with cold PBS twice. The cell pellets were lysed by 8M Urea in 20mM Tris-HCl (pH 8.0) with protein inhibitor cocktail and phosphor-STOP cocktail (Roche) to denature the proteins and sonicated twice on ice to break the protein-protein or DNA-protein interactions. All the procedures were carried out within the certified laboratory for studies of infectious materials. For all the sample sets, individual patient PBMC samples of each category (health, MS or SS) were pooled into several pools to achieve the minimal starting amount of total protein. 60 Basically, protein samples were combined to generate 3 control sample pools from 6 healthy donors (one sample pool was excluded from the list for the sample quality issue), 13 mild-case sample pools from 22 COVID-19-MS patients, and 10 severe-case sample pools from 13 COVID-19-SS patients according to the protein concentrations of all the samples due to the limited protein amounts extracted from the variable PBMCs of each case (Table S2 ). The resulted lysates were considered individual cases and subjected to protein digestion for the following proteomic experiments.

Protein Digestion: The protein concentration was measured, and equal amount of proteins were used for tryptic digestion. TCEP [tris(2-carboxyethyl) phosphine, final concentration is 5 mM] (Thermo Scientific) and iodoacetamide (final concentration is 10mM) (Sigma) for reduction and alkylation were added to the solution and incubated at room temperature for 30 minutes, respectively. The protein mixture was diluted four times and digested with LysC at 1:50 (w/w) for two hours and then with Trypsin at 1:40 (w/w) (Promega) at 37 °C overnight. The tryptic-digested peptides soups were desalted by MonoSpin TM C18 (GL Science) and the eluents were dried out in speed vacuum. The peptides mixtures were finally dissolved in HEPES buffer (100 mM pH 8.0).

The peptide concentrations were measured with Quantitative Colorimetric Peptide Assay (Thermo Scientific) before TMT-labeling to make sure the amount of peptides in each channel for TMT labeling is equal. TMT10plex amino reactive reagents (0.8 mg per vial, Thermo Scientific) were suspended in 41 µl of anhydrous acetonitrile respectively and suitable for the labeling of 100 µg peptide soups in each channel. According to the amount of the peptide mixtures, equal amounts of the peptide mixtures (~90 µg) were labeling with three sets of the TMT10 labeling kit. Meanwhile, a peptide mixture consisted of aliquots of each sample was also labeled with TMT10 and spiked into each set of quantitation experiments as internal standards for protein normalization. Reactions were allowed to proceed at room temperature for 1 h, and then quenched by an additional 8 µl of 5% hydroxylamine for 15 min. The TMT-labeled samples were pooled at a 1:1:1:1:1:1:1:1:1:1 ratio. The mixture was vacuum centrifuged to near dryness and desalted on a MonoSpin TM C18 (GL Science).

The desalted TMT-labeled peptides were resuspended in 300 ul 10% buffer C (150 mM ammonium hydroxide, pH 10). Off-line fractionations were performed using Agilent 1260 HPLC with a 4.6 mm × 250 mm BEH C18, 3.5 µm column (Waters) by a 58 min gradient from 5% to 80% buffer B (100% acetonitrile) with buffer A (100% water) at a flow rate of 0.5 mL/min while buffer C was running constantly at 10% of total flow rate during the process. Fractions were collected over 100 min at 2 min intervals after the start of the gradient in 1.5 ml Eppendorf tubes. The peptide mixture was fractionated into a total of 50 fractions which were consolidated into 18 vials. Samples were subsequently desalted and vacuum centrifuged to near dryness. The desalted peptides were reconstituted in 0.1% formic acid for downstream LC-MS/MS analysis.

For proteome study, all the 18 fractions for each TMT experiment were analyzed by a home-made 40 cm-long pulled-tip analytical column (75 µm ID x 360µm OD, ReproSil-Pur C18-AQ 1.9 µm resin, Dr. Maisch GmbH), the column was then placed in-line with an Easy-nLC 1200 nano HPLC (Thermo Scientific) for mass spectrometry analysis. The analytical column temperature was set at 55 ℃ during the experiments. The mobile phase and elution gradient used for peptide separation were as follows: 0.1% formic acid in water as buffer A and 0.1% formic acid in 80% acetonitrile as buffer B, 0-1 min, 5%-10% B; 1-96 min, 10-40% B; 96-104 min, 40%-60% B, 104-105 min, 60%-100% B, 105-120 min, 100% B. The flow rate was set at 300 nL/min. Data-dependent tandem mass spectrometry (MS/MS) analysis was performed with a Q Exactive Orbitrap mass spectrometer (Thermo Scientific). Peptides eluted from the LC column were directly electrosprayed into the mass spectrometer with the application of a distal 2.2-kV spray voltage. A cycle of one full-scan MS spectrum (m/z 300-1800) was acquired followed by top 20 MS/MS events (first mass fixed at 100 m/z scan range), sequentially generated on the first to the twentieth most intense ions selected from the full MS spectrum at a 35% normalized collision energy. Full scan resolution was set to 70,000 with automated gain control (AGC) target of 3e6. MS/MS scan resolution was set to 35,000 with an isolation window of 1.8 m/z and AGC target of 1e5. The number of micro scans was one for both MS and MS/MS scans and the maximum ion injection time was 50 and 100 ms, respectively. The dynamic exclusion settings used were as follows: charge exclusion, 1 and >8; exclude isotopes, on; and exclusion duration, 30 seconds.

The acquired MS/MS data were analyzed against a UniProtKB human (database released on Sept. 30, 2018) by Maxquant V1.6.10.43 using the default setting. Tolerance of precursor mass and fragment mass were set to ±20ppm. The main search peptides tolerance was setting at 4.5ppm according to the feature of the instrument in the study. Carbamidomethylation of cysteine (+57.021 Da) and acetylation of protein N-terminus were set as static modifications. The quantification search type chooses the reporter ion ms2 using TMT 10plex method. Oxidation of methionine (+15.995 Da) was set as variable modification. Trypsin was defined as cleavage enzyme, and the maximum missing cleavage was set at 2. All identified proteins had an FDR ≤1%, which was calculated at the peptide level. The unique peptides were selected for protein quantitation.

The interpretation of three sets of proteomic results (protein Groups.txt) from MaxQuant was done in the Perseus software environment (version 1.6.2.3) and R4.0.0. Briefly, the datasets were merged by matching the protein accession numbers of each dataset using in-house developed Python scripts (comparepg.py). The protein abundances were normalized by the TMTlabeled internal standard samples which consisted of an aliquot of each sample and were spiked into the 3 sets of TMT labeling experiments. Protein features with less than 3 quantitative values in each condition (H&MS or MS&SS) were removed from the data set to make sure each quantitative protein in the list have at least 3 values for H cases, MS cases, and SS cases during the proteomics analysis. The remaining missing quantitative values were replaced by values imputed using multivariate imputation algorithms from the MICE R package. 61 Data were further normalized across the TMT 10plex sets using the median normalization method. 62 All the normalized reporter ion abundances of biological repeats were averaged for fold-change calculation. For proteome analysis, P-value was calculated by empirical Bayes moderated t-tests and then adjusted by Benjamini and Hochberg method, using the limma package 63 built for R environment. Proteins were required to have a |log2 (fold change) |≥ 0.58 and adj. P-value cut-off of 0.01 to be considered differentially expressed. 64 Gene ontology (GO) analysis was performed using ClusterProfiler v3.14.3. 58

28 SAR-CoV-2 genes were cloned into both pGBTK7 (bait) and pACT2 (prey) vectors. The reporter Y2HGold yeast strain (Clontech) was co-transformed with different bait-prey combinations as indicated in Figure S1 in accordance with instructions for Matchmaker GAL4 Two-Hybrid System 3 (Clontech). The transformed yeast was then screened in 4-dropout plates for the protein-protein interaction.

Total RNA was isolated from cells with the RNeasy Mini Kit (Qiagen, #74106) and treated with RNasefree DNase according to the manufacturer's protocol. All gene transcripts were quantified by quantitative PCR using qScript TM One-Step qRT-PCR Kit (Quanta Biosciences, #95057-050) on CFX96 real-time PCR system (Bio-Rad). Primer sequences for qPCR were as follows: IL- 

160 SARS-CoV-2 isolates with complete genome sequences from 18 different countries collected from December 23, 2019 to February 29, 2020 were used for molecular evolutionary analysis in this study. Multiple alignments were performed by Clustal W and a phylogenetic tree was constructed by Neighbor-Joining statistical method based on the Maximum Composite Likelihood model with 500 replicates of bootstrap testing using MEGA. 55

For AP-MS analysis (Figure 2 ), proteins with spectral counts no more than 8 and fold change no more than 5 were excluded as background. MiST algorithm was used to calculate the statistical significance of the dataset. 34 To screening the high-confidence interactions, the MiST score setting at 0.8 was used to filter the interacting host targets. The ribosome proteins were also removed from the list due to their high background for protein purification. The interactome network was generated by Cytoscape software (https://cytoscape.org/). For statistics of proteome analysis ( Figures 3A-F and S4A-F) , The acquired MS/MS data were analyzed against a UniProtKB human (database released on Sept. 30, 2018) by Maxquant V1.6.10.43 using the default setting. The interpretation of three sets of proteomic results (protein Groups.txt) from MaxQuant was done in the Perseus software environment (version 1.6.2.3) and R4.0.0. P-value was calculated by empirical Bayes moderated t-tests and then adjusted by Benjamini and Hochberg method, using the limma package 63 built for R environment. Proteins were required to have a |log2 (fold change) |≥ 0.58 and adj. P-value cut-off of 0.01 to be considered differentially expressed. 64 All data were expressed as mean ± SEM when appropriate. For parametric analysis, the F test was used to determine the equality of variances between the groups compared; statistical significance across two groups was tested by Student's t-test ( Figures 3K and S4J-K) ; one-way analysis of variance (ANOVA) followed by Bonferroni's post hoc test was used to determine statistically significant differences between multiple groups ( Figures S3G-H, 3J, 3L, and S4H) . P-values of less than 0.05 were considered significant. Figure S1 . Summary of SARS-CoV-2 intra-viral PPIs, Related to Figure 1 Virus-host PPI network

A new coronavirus associated with human respiratory disease in China

A pneumonia outbreak associated with a new coronavirus of probable bat origin

A Novel Coronavirus from Patients with Pneumonia in China

Travelers Give Wings to Novel Coronavirus (2019-nCoV)

Emerging understandings of 2019-nCoV

Clinical characteristics of 2019 novel coronavirus infection in China

A Dynamic Immune Response Shapes COVID-19 Progression

Immunology of COVID-19: current state of the science

SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor

Structure, Function, and Antigenicity of the SARS-CoV-2 Spike Glycoprotein. Cell

Structural basis for the recognition of the SARS-CoV-2 by full-length human ACE2

Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved αketoamide inhibitors

Clinical features of patients infected with 2019 novel coronavirus in

Identification of a novel coronavirus causing severe pneumonia in human: a descriptive study

Viral and host factors related to the clinical outcome of COVID-19

Neutrophils as emerging therapeutic targets

Neutrophil recruitment and function in health and inflammation

Aveolar Macrophage Activation and Cytokine Storm in the Pathogenesis of Severe COVID-19

The Architecture of SARS-CoV-2 Transcriptome

Genome Composition and Divergence of the Novel Coronavirus (2019-nCoV) Originating in China

Analysis of intraviral protein-protein interactions of the SARS coronavirus ORFeome

Genome-wide analysis of protein-protein interactions and involvement of viral proteins in SARS-CoV replication

Structure of the RNA-dependent RNA polymerase from COVID-19 virus

Structure of replicating SARS-CoV-2 polymerase

Structural basis for RNA replication by the SARS-CoV-2 polymerase

Crystal structure and functional analysis of the SARS-coronavirus RNA cap 2'-O-methyltransferase nsp10/nsp16 complex

Structural basis and functional analysis of the SARS coronavirus nsp14-nsp10 complex

Coronavirus nsp10/nsp16 Methyltransferase Can Be Targeted by nsp10-Derived Peptide In Vitro and In Vivo To Reduce Replication and Pathogenesis

Proteomics of SARS-CoV-2-infected host cells reveals therapy targets

Detection of SARS-CoV-2-Specific Humoral and Cellular Immunity in COVID-19 Convalescent Individuals

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

Heightened Innate Immune Responses in the Respiratory Tract of COVID-19 Patients

Scoring Large-Scale Affinity Purification Mass Spectrometry Datasets with MiST

A SARS-CoV-2 protein interaction map reveals targets for drug repurposing

Stress granule components G3BP1 and G3BP2 play a proviral role early in Chikungunya virus replication

VSV disrupts the Rae1/mrnp41 mRNA nuclear export pathway

IKKepsilon and TBK1 are essential components of the IRF3 signaling pathway

TRIM59 interacts with ECSIT and negatively regulates NF-κB and IRF-3/7-mediated signal pathways

The NF-κB Repressing Factor Is Involved in Basal Repression and Interleukin (IL)-1-induced Activation of IL-8 Transcription by Binding to a Conserved NF-κB-flanking Sequence Element

NF-κB repressing factor downregulates basal expression and mycobacterium tuberculosis induced IP-10 and IL-8 synthesis via interference with NF-κB in monocytes

NF-kappaB activation by a signaling complex containing TRAF2, TANK and TBK1, a novel IKK-related kinase

Mechanisms of STAT protein activation by oncogenic KIT mutants in neoplastic mast cells

Platelet factor 4 (CXCL4) seals blood clots by altering the structure of fibrin

A potent inhibitor of endothelial cell proliferation is generated by proteolytic cleavage of the chemokine platelet factor 4

Genetics of platelet reactivity in normal, healthy individuals

Protein C inhibitor acts as a procoagulant by inhibiting the thrombomodulin-induced activation of protein C in human plasma

Protein C inhibitor secreted from activated platelets efficiently inhibits activated protein C on phosphatidylethanolamine of platelet membrane and microvesicles

Chemotaxis of human monocyte-derived dendritic cells to complement component C1q is mediated by the receptors gC1qR and cC1qR

HCV core protein interaction with gC1q receptor inhibits Th1 differentiation of CD4+ T cells via suppression of dendritic cell IL-12 production

Inhibition of RIG-I and MDA5-dependent antiviral response by gC1qR at mitochondria

Agonist-induced platelet procoagulant activity requires shear and a Rac1-dependent signaling mechanism

COVID-19 Coagulopathy in Caucasian patients

Evolution and variation of 2019-novel coronavirus

MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms

Effects of IL-1β, IL-6 and IL-8 on erythrocytes, platelets and clot viscoelasticity

Novel Role of vBcl2 in the Virion Assembly of Kaposi's Sarcoma-Associated Herpesvirus

clusterProfiler: an R package for comparing biological themes among gene clusters

iProX: an integrated proteome resource

Comparative top down proteomics of peripheral blood mononuclear cells from kidney transplant recipients with normal kidney biopsies or acute rejection

Multiple imputation by chained equations for systematically and sporadically missing multilevel data

Normalyzer: a tool for rapid evaluation of normalization methods for omics data sets

limma powers differential expression analyses for RNA-sequencing and microarray studies

Proteome-Wide Evaluation of Two Common Protein Quantification Methods

The COVID-19 pandemic is caused by SARS-CoV-2, but little is known about the functions of its viral proteins. The authors characterized the SARS-CoV-2 intra-viral and virus-host interaction networks in human cells and identified 286 potential host targets. Quantitative proteomic analysis revealed elevated levels of IL-6 and IL-8 in PBMCs collected from severe COVID-19 patients compared with mild ones, and functional annotation of differentially expressed proteins implicate pathways involved in neutrophil activation, T-cell receptor signaling, and the coagulation cascade. Combining virus-host interactome with COVID-19 proteomic analysis, the authors found that nsp10 interacts with NKRF to mediate IL-8 expression, providing a potential molecular mechanism for SARS-CoV-2-induced cytokine storm and marking it as a possible emerging therapeutic target.