{'fraction_non_alphanumeric': 0.047117732766935504, 'fraction_numerical': 0.02994110283404133, 'mean_word_length': 5.074545001246572, 'pattern_counts': {'":': 0, '<': 22, '<?xml version=': 0, '>': 10, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 1, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Multi-organ proteomic landscape of COVID-19 autopsies

The molecular pathology of multi-organ injuries in COVID-19 patients remains unclear, preventing effective therapeutics development. Here, we report a proteomic analysis of 144 autopsy samples from seven organs in 19 COVID-19 patients. We quantified 11,394 proteins in these samples, in which 5,336 were perturbed in the COVID-19 patients compared to controls. Our data showed that cathepsin L1, rather than ACE2, was significantly upregulated in the lung from the COVID-19 patients. Systemic hyperinflammation and dysregulation of glucose and fatty acid metabolism were detected in multiple organs. We also observed dysregulation of key factors involved in hypoxia, angiogenesis, blood coagulation, and fibrosis in multiple organs from the COVID-19 patients. Evidence for testicular injuries includes reduced Leydig cells, suppressed cholesterol biosynthesis, and sperm mobility. In summary, this study depicts a multi-organ proteomic landscape of COVID-19 autopsies that furthers our understanding of the biological basis of COVID-19 pathology.


In Brief
A proteomics analysis of 144 autopsy samples from seven organs in 19 COVID-19 patients shows elevated expression of cathepsin L1, rather than ACE2, in the lung tissue and highlights dysregulation of angiogenesis, coagulation, and fibrosis in multiple organs, in addition to systemic hyperinflammation.

INTRODUCTION
The ongoing COVID-19 pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), led to more than 55 million infected individuals and over 1.3 million deaths by the middle of November 2020. Morphological characterization of autopsies, mainly focused on the pulmonary lesions, has greatly advanced our understanding of COVID-19-caused deaths (Carsana et al., 2020;Su et al., 2020;Wichmann et al., 2020;Wu et al., 2020;Xu et al., 2020;Yao et al., 2020). Mechanistic studies of SARS-CoV-2-infected cell line models (Bojkova et al., 2020;Bouhaddou et al., 2020;Gordon et al., 2020) offered new insights into virus-perturbed biochemical processes of COVID-19 and suggested potential therapies. SARS-CoV-2-infected mouse models Hassan et al., 2020;Jiang et al., 2020) and rhesus macaque models (Chandrashekar et al., 2020;Deng et al., 2020) generated by adenovirus transduction of human ACE2 have been established for preclinical selection of antiviral therapeutic agents and vaccines as well as for investigating pathogenesis. Few studies have characterized host responses at the molecular level from clinical specimens. We and others have studied the host responses by proteomic and metabolomic analysis of patient sera (Messner et al., 2020;Shen et al., 2020), but molecular changes in infected tissues and consequentially affected organs remain elusive. To date, little knowledge has been obtained concerning how SARS-CoV-2 virus induces injuries in multiple organs (Bian, 2020;Tian et al., 2020;Wichmann et al., 2020) including lung, kidney (Kudose et al., 2020), liver, heart, spleen, thyroid, and testis , and how to prevent and revert them. Ackermann et al. (2020) analyzed the lung transcriptome of seven COVID-19 autopsies and reported intussusceptive angiogenesis that may be induced by hypoxemia. Latest advances of proteomics technologies allow effective and robust analysis of formalin-fixed tissue samples Zhu et al., 2019). Comparative analysis of mRNA and protein expression in tissue samples showed that proteins measured by mass spectrometry were much more stable than transcripts (Shao et al., 2019). Here, we report a multi-organ proteomic profiling of 144 autopsy tissue samples collected from the lung, spleen, liver, heart, kidney, thyroid, and testis of 19 patients died from COVID-19 and 74 control tissue samples from 56 non-COVID-19 patients. Using tandem mass tagging (TMT)-based shotgun proteomics, we quantified 11,394 proteins, of which 5,336 were significantly dysregulated in at least one organ in the COVID-19 patients. This data resource offers a unique channel for understanding multi-organ injuries in the COVID-19 patients and nominating potential therapeutics.

RESULTS AND DISCUSSION
Generation and characterization of proteomic landscape We first performed proteomic profiling of 144 autopsy tissue samples from seven organs, namely lung (N = 15, n = 30; N represents the number of patients, n represents the number of samples), spleen (N = 9, n = 9 from white pulp and red pulp and N = 8, n = 8 from red pulp), liver (N = 10, n = 24), kidney (N = 10, n = 18 from renal cortex and N = 10, n = 16 from renal medulla), heart (N = 9, n = 19), testis (N = 5, n = 5), and thyroid (N = 15, n = 15). The samples were from 19 COVID-19 cases, ten of which have been described previously , compared with 74 control samples from 56 non-COVID-19 cases with other diseases via surgeries (Figure 1A; Tables S1 and S2). All 19 COVID-19 patients died from SARS-CoV-2 pneumonia or respiratory failure, among which seven also developed terminal multiple organ dysfunction syndrome (MODS). Other confounders, including the effects of smoking, gender, hypertension, diabetes, and coronary heart disease, all showed no difference between the COVID-19 and non-COVID-19 groups by Fisher's exact tests (Table S1). Due to missing values of body mass index (BMI) information, this study could not be utilized to investigate the impact of BMI. Detailed information of patients including medication history during hospitalization, laboratory test data, pathological changes, and cause of death are summarized in Figures S1 and S2 and Tables S1 and S5.
Altogether, we quantified 11,394 proteins from above samples with a false discovery rate (FDR) less than 1% at both peptide and protein levels ( Figure 1A; Table S3). The number of identified proteins ranged from 5,828 (heart) to 9,544 (kidney) across seven types of organs. We analyzed 37 technical replicates of randomly selected tissue samples, as well as 18 pooled controls for each TMT batch ( Figure S3A; Table S2). Proteins quantification in these technical replicates and control samples showed relatively low median coefficient of variance (CV) of 6.88% and 2.47%, respectively ( Figures S3B and S3C). A total of 5,336 dysregulated proteins were characterized from the seven types of organs between COVID-19 and control groups (Benjamini-Hochberg [B-H] adjusted p value <0.05 and |log 2 (fold change) | >log 2 [1.2]) (Figure 1A; Table S4). Splenic red pulp samples were excluded from downstream analysis, because they did not show any statistically significant proteomic regulations ( Figure 1A). Hierarchical clustering of the differentially expressed proteins from each organ type (Figure S3D) showed that these proteins well separated COVID-19 samples and controls. The t-distributed stochastic neighbor embedding (t-SNE) analysis shows that dysregulated proteomes of each organ type clustered tightly apart from each other (Figure S3E), consolidating that these selected proteins well-resolved different organ types. Except for testis, the other six organ types shared only 27 dysregulated proteins, suggesting different organs responded to viral infection via diverse pathways. Among them, the acute inflammatory protein C-reactive protein (CRP) and the scavenger receptor cysteine-rich type 1 protein M130 (CD163), which is highly expressed marker in M2 type macrophage (Etzerodt and Moestrup, 2013), were the most upregulated proteins (Figure 1B), probably reflecting both the hyperinflammatory and repairing state in patients died from COVID-19.
In our study, 14 out of 19 COVID-19 patients exhibited sepsis (Table S1), which is associated with systemic hyperinflammatory response, multi-organ injuries, and MODS (Hotchkiss et al., 2016). We compared the multi-organ proteomes of the COVID-19 patients with MODS (N = 7) to those without MODS (N = 12). As shown in Table S4, three types of organs including liver, renal cortex, and testis showed significant dysregulation of protein expression between the two groups. In the liver of MODS group, most of elevated proteins in regulated exocytosis are involved in immune response, such as acute phase response, cytokine secretion, and neutrophil degranulation (Table S6), suggesting dysfunction of the first line of host defense. These inflammatory cells and mediators all contribute to the development of MODS (Lobo et al., 2003).

Six clusters of proteins relevant to SARS-CoV-2 infection
We then focused on six clusters of proteins including viral receptors and proteases, transcription factors (TFs), cytokines (and their receptors), coagulation system-related proteins, angiogenesis-associated proteins, and fibrosis markers due to their relevance to SARS-CoV-2 infection (Figures 2A and 2B). After cellular entry mediated by receptors and proteases, SARS-CoV-2 hijacks the host translation machinery and induces host inflammatory response via TFs, leading to hyper-inflammatory state, which might be associated with the clinically observed blood hypercoagulability as measured by blood tests ( Figure S2; Table S5), fibrosis and microthrombosis as examined by ll pathologists (Figure S1), and enhanced angiogenesis as reported by morphologic and gene expression examination in the lung of COVID-19 (Ackermann et al., 2020). More details of the regulated proteins are provided in the supplementary tables.
Our data showed substantial regulation of TFs in COVID-19 autopsies. 395 out of 1,117 quantified TFs were altered in at least one tissue type ( Figure S4A; Table S5), and they were signifi-cantly enriched in spliceosome and viral carcinogenesis, among others, as shown in Figure S4B. By matching the experimental fold-change with the predicted activation state in Ingenuity Pathway Analysis (IPA), ten of these dysregulated TFs showed the same regulatory trend. Our data showed that six of these proteins are involved in the inflammatory responses and all of them were upregulated in multiple COVID-19 organs. The six The quantified and dysregulated proteins across multiple organs. The outermost (first) ring represents the type of samples. The number of samples and patients (n/N) is labeled respectively. The second ring (in blue) refers to the missing/undetected proteins for each organ. The numbers in black represent the quantified proteins in the specific organ. The third ring (in light green) refers to unregulated proteins . The numbers in white represent the significantly dysregulated proteins in specific organ type (B-H adjusted p value <0.05; |log 2 [fold change of COVID-19 versus non-COVID-19]| >log 2 [1.2]). The innermost ring refers to the number of significantly dysregulated proteins for each organ (pink, upregulated; dark green, downregulated). (B) Protein expression of CRP and CD163 across six organs (except testis). The y axis stands for the protein expression ratio by TMT-based quantitative proteomics. Pairwise comparison of each protein between COVID-19 and non-COVID-19 groups was performed with Student's t test. The cutoff of dysregulated proteins has been set at B-H adjusted p value <0.05 and |log 2 (FC)| >log 2 (1.2). B-H adjusted p value: *p < 0.05; **p < 0.01; ***p < 0.001. See also Figures S1, S2, and S3 and Tables S1, S2, S3, and S4 proteins include nuclear factor kB (NF-kB) subunit 2 (NFKB2) (Lawrence, 2009), transcription factor p65 (RELA) (Rayet and Gé linas, 1999), CCAAT/enhancer-binding protein b (C/EBPB) (Greenbaum et al., 1998), signal transducer and activator of transcription 1 and 3 (STAT1/3) (Yu et al., 2009), and transcription factor jun-B (JUNB) (Hess et al., 2004) (Figure 2B; Table S5). We also detected dysregulation of TFs involved in tissue injuries and hypoxia state. The RB transcriptional corepressor 1 (RB1) was upregulated in the liver and renal cortex ( Figures 2B and  S4A), which may negatively modulate the cell cycle (Dyson, 2016) and induce a higher degree of mitochondrial membrane permeabilization and apoptosis (Hilgendorf et al., 2013). The La-related protein 1 (LARP1) is an RNA-binding protein regulating mRNA translation and modulated by mTOR signaling (Tcherkezian et al., 2014). Recently, LARP1 has been reported to interact with SARS-CoV-2 (Gordon et al., 2020). Interestingly, our data showed that it was upregulated in the lung, spleen, liver, renal medulla, and thyroid of these COVID-19 patients (Figure 2B), suggesting it might be a critical player in the SARS-CoV-2 replication. The hypoxia-inducible factor 1-alpha (HIF1A) was upregulated in renal cortex (Figures 2B and S4A) and was also predicted to be activated in the lung, liver, renal cortex, and medulla by upstream regulator analysis, although not detected in these samples, suggesting systemic hypoxia state (Lee et al., 2019b) in these COVID-19 patients. Hepatocyte nuclear factor 4-alpha (HNF4A) can modulate the transcription of genes regulating lipid and bile acid synthesis and gluconeogenesis (Babeu and Boudreau, 2014). The hypoxia state and upregulated STAT3 could reduce the expression of HNF4A (Cairo and Buendia, 2012;Mazure et al., 2001), leading to downregulation of microsomal triglyceride transfer protein large subunit (MTTP) and apolipoprotein B-100 (ApoB), and upregulation of scavenger receptor class B member 1 (SCARB1) (Hayhurst et al., 2001), as confirmed by our proteomic data (Figure 2B; Table S4). Reduced MTTP and ApoB and increased SCARB1 have been associated with hepatic steatosis (Hayhurst et al., 2001) ( Figure S1).
Among the 242 quantified cytokines and related proteins, 112 were significantly dysregulated and enriched in angiogenesis, response to growth factor, and other pathways ( Figures  S4C and S4D; Table S5). Nicotinamide phosphoribosyl transferase (NAMPT), glucocorticoid receptor (NR3C1), and interferon-gamma receptor 1 (IFNGR1) were dysregulated in most organs. NAMPT, which participates in multiple signaling pathways (e.g., interleukin [IL]6-STAT3 and NF-kB) (Garten et al., 2015), was upregulated in the six organs except for testis. NR3C1 is the receptor for glucocorticoid and a transcriptional factor promoting anti-inflammatory responses (Baschant and Tuckermann, 2010). It is classified as cytokine and associated proteins in the ImmPort database. The downregulation of NR3C1 was identified in the five organs except for testis and thyroid ( Figure S4C), suggesting reduced anti-inflammatory responses. IFNGR1, which triggers host immune responses on viral infection (Xia et al., 2018) (e.g., induction of phagocyte oxidase system, nitric oxide [NO] production, and lysosomal enzymes activation for microbe destruction in macrophages) (Schroder et al., 2004), was upregulated in the five organs except for testis and thyroid.
CTSL, rather than ACE2, was upregulated in lungs Our data identified six reported potential receptors or proteases for the virus entry ( Figure 2C), namely angiotensin-converting enzyme 2 (ACE2) (Hoffmann et al., 2020), C-type lectin domain family 4 member M (CLEC4M) (Jeffers et al., 2004) and member L (CD209) (Yang et al., 2004), Niemann-Pick C1(NPC1) (Cô té et al., 2011), carcinoembryonic antigen-related cell adhesion molecule 1 (CEACAM1) (Tsai et al., 2003), and cathepsin L1 (CTSL) (Liu et al., 2020b;Ou et al., 2020). ACE2, the known receptor mediating SARS-CoV-2 entry, did not show significant regulation in the lung, suggesting that ACE2 inhibitors might not be an effective therapy for severe and critical COVID-19 patients. However, it was downregulated in both kidney and heart ( Figure 2C), which may be associated with its modulatory roles on angiotensin II, including inflammation, vasoconstriction, and thrombosis (Liu et al., 2020a). The other two receptors for SARS (Chan et al., 2006;Jeffers et al., 2004), namely CD209 and CLEC4M, were not significantly dysregulated in the lung either ( Figure 2C). Interestingly, our results showed CTSL, the serine protease of SARS-CoV-2 in the endosomal pathway, was significantly upregulated in the lung ( Figure 2C), nominating it as a potential therapeutic target for COVID-19 (Bittmann et al., 2020). In the SARS-CoV-2 infection cell model, the inhibitor of CTSL has been proven to be effective for blocking the virus entry (Ou et al., 2020). Besides CTSL, our data showed upregulation of several other cathepsins including CTSB, CTSD, CTSE, CTSH, CTSK, CTSS, and CTSZ in the lung of the COVID-19 patients (Table S4). Cathepsins have been reported to facilitate the synthesis of pro-IL-1b and activate NLRP3 inflammasome, leading to pyroptosis (Orlowski et al., 2015;Yap et al., 2020). IL-1 signaling was also enriched and activated in the spleen, renal medulla, and thyroid of the COVID-19 patients (Table S6). In the lung of COVID-19 patients, myeloperoxidase (MPO) for neutrophil primary granules was upregulated (Table S4), indicating infiltration of neutrophils recruited by IL-1. Interestingly, Figure 2. Six functional clusters of dysregulated proteins from seven organs between COVID-19 and non-COVID-19 patients (A) Counts of dysregulated proteins in six clusters of molecules, including potential virus receptors and proteases, fibrosis markers, cytokines (and their receptors), transcription factors (TFs), coagulation system, and angiogenesis-associated proteins are shown in a bar chart. Each column along y axis represents a type of organ. The number of proteins is shown in x axis. (B) Landscape of 5336 significantly dysregulated proteins in seven organs. The dysregulated proteins in the six clusters are labeled as circles (solid, upregulated proteins; hollow, downregulated proteins). The size of circle indicates |log 2 (FC)|. (C) Protein expression of potential virus receptors across multiple organs. The y axis stands for the protein expression ratio by TMT-based quantitative proteomics. Pairwise comparison of each protein between COVID-19 and non-COVID-19 groups was performed using Student's t test. The cutoff of dysregulated proteins has been set at B-H adjusted p value <0.05 and |log 2 (FC)| >log 2 (1.2). B-H adjusted p value: *p < 0.05; **p < 0.01; ***p < 0.001. See also Figures S4 and S5 and Table S5. we also found upregulation of lactate dehydrogenase (LDH) in the lung (Table S4), which is an indicator of pyroptosis and tissue injuries. Thus, the inhibition of CTSL may also contribute to the excessive inflammatory activity in the COVID-19 patients. Of note, we also observed upregulation of CTSL in the spleen, renal medulla, and thyroid ( Figure 2C). . Coagulation, angiogenesis-associated proteins, and fibrosis markers regulated in multiple organs (A-C) Chord diagrams show dysregulated and multi-organ shared proteins in coagulation system (A), angiogenesis associated proteins (B), and potential fibrosis markers (C) between COVID-19 and non-COVID-19 patients across multiple organs. The cutoff of dysregulated proteins has been set at B-H adjusted p value <0.05 and | log 2 (FC)| >log 2 (1.2). The length of the brick for each protein corresponds to the sum of |log 2 (FC)| in multiple organs. The length of the brick for each organ corresponds to the sum of |log 2 (FC)| in one or more proteins. See also Figure S5.

Multi-organ coagulation, angiogenesis, and fibrosis
Our data suggested systematic dysregulation of coagulation, angiogenesis, and fibrosis in COVID-19 patients, as shown in Figures 3 and S5 and Table  S5. Microthrombi in the COVID-19 patients were observed in the lung, kidney (Figure S1), and lower extremity veins (Table S1) of the COVID-19 patients, in agreement with laboratory test results, such as elevated D-dimer ( Figure S2). The formation of microthrombi is due to imbalance among coagulation, anticoagulation, and fibrinolytic systems (Palta et al., 2014). Our data uncovered multiple dysregulated proteins participating in coagulation, anticoagulation, and fibrinolytic system that might have contributed to the coagulation disorders in COVID-19 ( Figure 3A; Table S5). Four coagulation factors, namely coagulation factors prothrombin (F2), XI, XII, and XIIIa (F11, F12, and F13A1), were dysregulated in the COVID-19 patients ( Figure 3A). F13A1 is activated in the last step of coagulation, which could induce hemostasis and stabilize the fibrin clot to avoid fibrinolysis (Muszbek et al., 2011). We detected increase of F13A1 in the renal cortex, which might be associated with the observed blood clotting in the renal cortex ( Figure S1). VWF is a glycoprotein that binds to the factor VIII (F8), protects F8 from degradation by vitamin K-dependent protein C (PROC), and could trigger the platelet aggregation following vascular injury (Peyvandi et al., 2011). We found that VWF was increased in the renal cortex in COVID-19 ( Figure 3A), indicating higher risk of thrombosis. Our data also showed increase of fibrinogen alpha chain, gamma chain, and beta chain (FGA, FGG, and FGB) in the COVID-19 lung ( Figure 3A). These proteins could be cleaved into fibrin that contributes to formation of blood clots (Gö bel et al., 2018). We detected dysregulation of multiple serine protease inhibitors in the COVID-19 organs (Figure 3A). Heparin cofactor 2 (SERPIND1) is a serine proteinase inhibitor, which acts as the inhibitor of thrombin (He et al., 2002) and the cofactor for heparin. Plasminogen activator inhibitor 1 (SERPINE1) is a major inhibitor of plasmin, which could break down the blood clotting (Chapin and Hajjar, 2015). The downregulation of SERPIND1 and the upregulation of SERPINE1 in the renal cortex might have contributed to the microthrombi observed in the COVID-19 kidney cortex ( Figure S1).
Abnormal angiogenesis resulted from aberrant coagulation, and tissue hypoxia has been reported in the COVID-19 lungs with specialized technologies using microvascular corrosion casting coupled with scanning electron microscopy (Ackermann et al., 2020). In our study, a total of 139 angiogenesis-related proteins from the nCounter PanCancer Progression Panel (NanoSring Technologies) were significantly dysregulated (Figure 3B; Table S5), suggesting abnormal angiogenesis.
Fibrosis has been observed in the lungs of COVID-19 patients ( Figure S1; Table S1). It is usually divided into four stages (i.e., initiation, inflammation, proliferation, and modification) (Wynn, 2003). A cascade of stress and immune responses are triggered at the initiation stage, followed by activation of multiple inflammatory signaling pathways, such as the chemokine signaling, complement system, macrophage activation, and NF-kB signaling, among others. At the proliferation stage, fibroblasts differentiate and proliferate. Finally, the extracellular matrix (ECM) composed of immune cells and fibroblast cells is re-structured in the modification stage. In the COVID-19 samples, we characterized 335 dysregulated proteins involved in the four stages of fibrosis according to nCounter Fibrosis Panel (NanoSring Technologies), 179 of which are stage-specific proteins (Table S5). These 179 proteins formed 447 interactions according to String  and Cytoscape (Shannon et al., 2003) in our data (Figure S5). The dysregulated proteins associated with the fibrosis process in the lung have probably contributed to the pathological changes observed in these patients ( Figure S1). Although the proteins involved in the modification stage were most dominant in the lung, the other organs exhibited relatively less perturbation in these proteins, consistent with the microscopic examination ( Figure S1). Nevertheless, our data suggested the fibrosis processes have been triggered in these organs, although microscopically morphological changes were not discernable. These molecular changes might be exploited to instruct tissue fibrosis treatment for the COVID-19 patients.
Thrombospondin 1 (THBS1) is an extracellular glycoprotein involved in the inflammation and proliferation stages, and it is highly expressed by stromal fibroblasts and activated platelets (Adams and Lawler, 2011). Studies have reported its diverse roles in regulating cell-matrix interactions, platelet aggregation, angiogenesis, and formation of collagen matrix during wound healing (Bornstein, 2001;Sweetwyne and Murphy-Ullrich, 2012). In the COVID-19 patients, THBS1 was upregulated in the liver, heart, and kidney (Table S4), suggesting potential fibrosis in the COVID-19 patients.
In the modification stage, multiple lysosomal cathepsins, including CTSL and CTSD, have been reported to play a pro-fi-brogenic role in the liver, kidney, heart, and lung, which regulates the ECM degradation and tissue remodeling (Fox et al., 2016;Manchanda et al., 2017). The upregulated CTSL and CTSD in the COVID-19 patients (Table S5) might contribute to fibrosis in the lung, spleen, thyroid, liver, and heart.
In addition, we detected 29 fibrosis-associated proteins in our dataset according to manual literature mining ( Figure 3C; Table S5). SERPINE1 and Chitinase 3 Like 1 (CHI3L1) were upregulated in most organs from the COVID-19 patients ( Figure 3C). SERPINE1 is a potent inhibitor of proteolytic urokinase and tissue plasminogen activator (uPA/tPA). Elevated SERPINE1 could inhibit degradation of ECM, contributing to multi-organ fibrogenesis (Ghosh and Vaughan, 2012). CHI3L1 is a secreted chitinaselike protein modulating fibroblasts proliferation, immune cell differentiation, ECM reorganization, and angiogenesis in response to cytokines and stresses such as hypoxia . The upregulation of SERPINE1 and CHI3L1 were detected not only in the lung, but also in the liver, heart, and kidney, suggesting a multi-organ profibrotic state in the COVID-19 patients.
Dysregulated protein translation, glucose, and fatty acid metabolism To obtain a systematic understanding of biological processes represented by 5,336 dysregulated proteins, we performed pathway enrichment analysis for each tissue type using IPA. Comparisons of the most enriched or dysregulated pathways (-log 10 [p value] >10 or ratio >0.35 or absolute [Z score] >5) among seven organ types are shown in Figure 4A and Table  S6. EIF2 signaling is involved in the regulation of mRNA translation (Roux and Topisirovic, 2018), which has been reported to be affected by virus infection (Bojkova et al., 2020). The lung, liver, and thyroid shared a similar pattern of mRNA translation (Figure 4A), although there was little evidence of SARS-CoV-2 infection in the liver or thyroid. Most dysregulated proteins specific to the lung belonged to L13a-mediated translational silencing of ceruloplasmin expression ( Figures 4B and 6A), which has been reported as an innate immune mechanism after the virus infection (Mazumder et al., 2014). We observed suppression of multiple metabolic processes including glycogenolysis, galactose degradation, and glycolysis ( Figure 4A; Table S6). In contrast, fatty acid b oxidation (FAO) and oxidative phosphorylation were activated in most organs, suggesting a switch to high-efficiency energy production mode to support virus replication in the lung and mRNA translation in the liver (Heaton and Randall, 2011). In addition, dysregulated FAO and oxidative phosphorylation also led to excessive generation of reactive oxygen species (ROS) and release of pro-apoptotic proteins, which induced liver necrosis ( Figure S1). In the kidney, FAO was inhibited, which is thought to be a contributor to acute kidney injury (AKI)-induced renal fibrogenesis (Kang et al., 2015). Indeed, AKI has been observed in most COVID-19 patients in our study (Table S1).

SARS-CoV-2-associated protein regulation in the lung
The lung is the major target attacked by SARS-CoV-2. In our data, COVID-19 lung proteome showed unique enrichment of pathways that are known to be associated with virus infection, including mRNA decay and translation shutoff (Tanaka et al., 2012) (Figures 4C, 5A, and 6). The ARE-mediated mRNA The top pathways dysregulated across multiple organs. Pathway analysis was performed using all dysregulated proteins in the specific organ using IPA. The size of circle represents the -log 10 (p value) and the color represents the Z score by IPA. (B) The pathways enriched by Metascape for translation initiation relating proteins that are differentially expressed only in lung or liver, respectively. (C) Translation-associated pathway comparison across multiple organs. The size of circle represents the -log 10 (P) and the color represents the Z score by IPA. (D) Heatmap of SARS-CoV-2 interacting proteins dysregulated in the lung. The significance (''Sig.'' as the short term in figures) was calculated using Student's t test. B-H adjusted p value: *p < 0.05; **p < 0.01; ***p < 0.001. The cutoff of dysregulated proteins has been set at B-H adjusted p value <0.05 and | log 2 (FC)| >log 2 (1.2). See also Figure S7 and Table S6. ll degradation pathway was inhibited in the lung, while activated in all the other organs ( Figure 4C). During coronavirus replication, the translation of nested subgenomic mRNAs is cap-dependent and required for the cap-binding protein, eukaryotic translation initiation factor 4E (EIF4E), which involves cap recognition in the translation initiation. In addition, EIF4E has been reported as a potential target to block human coronavirus 229E (HCoV-229E) infection (Cencic et al., 2011). Our data showed that EIF4E was upregulated only in the lung of the COVID-19 patients ( Figure 5A), which might be perturbed and hijacked by SARS-CoV-2. Further, double-stranded RNA (dsRNA) and uncapped mRNA of the virus act as viral pathogen-associated molecular patterns (PAMPs) that trigger the innate immune response through recognition by pattern recognition receptors (PRRs) in the cytoplasm (Lin and Cao, 2020). By comparing with the SARS-CoV-2 protein interaction map (Gordon et al., 2020), 12 virus-host interacting proteins were dysregulated only in the lung, including stress granule-related factor G3BP1, mitochondrial protein TIMM10, transcription regulator eIF4H, RING-type E3 ubiquitin ligase MIB1, pro-inflammatory cytokine receptor IL17RA, and member of Cullin RING E3 ligase 2 complex ZYG11B ( Figure 4D). These proteins have been reported to promote virus replication, inhibit host mRNA expression, mediate the delivery of viral DNA through viral nuclear pore complex, participate in pulmonary fibrosis, and degrade virus restriction factors (Gordon et al., 2020). G3BP1 was reported to interact with SARS-CoV-2 nucleocapsid (N) protein, which could be sequestered by viruses to promote their replication (Gordon et al., 2020). TIMM10 has been reported to be targeted by viruses to enhance their replication (Williamson et al., 2012). The interaction between SARS-CoV-2 Nsp9 and eIF4H may indicate the inhibition of host mRNA expression (Gordon et al., 2020). The interaction between Nsp9 and MIB1 may mediate the delivery of viral DNA through nuclear pore complex (Gordon et al., 2020). IL17RA is associated with elevation of collagen and pulmonary fibrosis, and its inhibitor has been reported to reduce fibrosis in SARS infection (Mi et al., 2011). The interaction between Orf10 and ZYG11B may be hijacked for degradation of virus restriction factors or be blocked to protect itself from degradation (Gordon et al., 2020).

Immune responses in multiple organs
We next investigated the multi-organ immune responses in the COVID-19 patients based on our proteomics data (Figures 5,  6, and S7). We found that spleen and lung exhibited similar immune response patterns ( Figure 5A). In the lung, we detected upregulation of two immune checkpoint proteins, namely carcinoembryonic antigen-related cell adhesion molecule 1 (CEACAM1) and CD276 ( Figure 5A). CEACAM1 is a cell-cell adhesion protein expressed in lymphocytes and suppresses cytokine production, proliferation, and cytotoxic activity of T cells in response to virus (Gray-Owen and Blumberg, 2006). Another immune checkpoint protein CD274 was upregulated in the spleen ( Figure 5A). The elevation of checkpoint proteins suggests suppression of adaptive immunity in the lung and spleen of the COVID-19 patients. Further, in the lung and spleen, we detected downregulation of lymphocyte-specific tyrosine-protein kinase (LCK), which is enriched in T cells. Thus, the downregulated LCK suggests suppression of T cell-mediated immune response in the lung and spleen. Next, we performed pathway analysis of dysregulated proteins in the spleen, which highlighted activation of PD-1 and PD-L1 pathway, and inhibition of B cell receptor signaling ( Figure 6A; Table S6), further consolidating decreased adaptive immune response in the COVID-19 patients.
We further checked multiple markers for various immune cells in the spleen samples using immunohistochemistry (IHC) and found that the numbers of T and B lymphocytes were significantly reduced by IHC staining of CD3 for the total T cells, CD4 for CD4+ T cells, CD8 for CD8+ T cells, and CD20 for B cells, especially in the white pulp of COVID-19 patients, whereas the number of macrophages, M2 macrophages increased in the spleen of COVID-19 patients by IHC staining of CD68 and CD163, respectively ( Figure S6). These findings agreed with our proteomic data ( Figures 5A and 6A).
T cell exhaustion and upregulation of monocytes biomarkers in the spleen (Figures 5A and S6) suggested hyperinflammation, which may have damaged the integrity of the gas exchange barrier and induced hypoxia ( Figure 6A). The hypoxia state would further stimulate the inflammatory responses (Eltzschig and Carmeliet, 2011). Interestingly, our pathway analysis for the immunological proteins from GSEA-immunologic gene sets highlighted activation of NF-kB signaling and acute phase response in the liver (Figures 5B, 6A, and 6B; Table S6), probably induced by increased cytokines in the circulation (Alonzi et al., 2001;Israë l, 2010). We detected upregulation of the kinase subunit of IKK complex and its modulator subunit, namely IKBKG and IKBKB, in the liver. This also supports the activation of canonical NF-kB cascade. We observed vasogenic edema in the heart ( Figure S1), which might be due to vascular hyperpermeability induced by circulating cytokines, acute phase proteins, and other molecules including histamine (Nagy et al., 2008). Our proteomic data did not measure metabolites, but we found downregulation of histamine N-methyltransferase (HNMT) (Figures 6B and S7A), which catabolizes histamine, supporting accumulation of histamine.
Sepsis and its common complication, acute kidney injury (AKI), characterized by systemic inflammatory cascade (Poston and Koyner, 2019), were observed in most of the COVID-19 patients in our study (Table S1). In the renal cortex of the COVID-19 patients, we detected activation of multiple pathways involving inflammatory response, including the LPS/IL-1 mediated inhibition of RXR function, acute phase response, Toll-like receptor signaling, IL-6 signaling, and NF-kB signaling (Table S6). Of B-H adjusted p value: *p < 0.05; **p < 0.01; ***p < 0.001. See also Figure S6. ll note, Toll-like receptor 2 (TLR2) and its co-receptor CD14, which were reported to recognize the pathogenic molecules and mediate host innate immune responses (Oliveira-Nascimento et al., 2012), were both upregulated in the renal cortex ( Figures  6B and S7B).

Testicular injuries
Testis is one of the very few organs with immune privilege, and usually remains intact and unaffected from the host response during introduction of antigens (Zhao et al., 2014). Compared with the other organ types, the number of differentially expressed proteins in testes was small. Only ten proteins were regulated and all of them were downregulated ( Figure 7A), however, these changes suggest unusual pathological processes in the COVID-19 testis tissues. Insulin-like factor 3 (INSL3), the most abundantly expressed proteins in Leydig cells (Uhlé n et al., 2015), was the most dramatically decreased protein in the COVID-19 testicular tissue ( Figures 7A and 7D), suggesting impaired Leydig cell functions or a reduced Leydig cell population. Indeed, the histological examination revealed a reduction of Leydig cells ( Figures 7B and 7C), consistent with our previous pathological report . We also found five downregulated proteins related to cholesterol biosynthesis ( Figures 7A and 7D). All steroid hormones, including testosterone, are derived from cholesterol. E3 ubiquitin-protein ligase (RNF216) is essential for spermatogenesis and male fertility (Melnick et al., 2019). The dynein regulatory complex subunit 7 (DRC7) is a sperm motility factor and its deletion leads to aberrant tail formation in mouse spermatozoa that phenocopies patients with multiple morphological abnormalities of the sperm flagella (MMAF) (Morohoshi et al., 2020). These two proteins were reduced in the COVID-19 patients ( Figures 7A and 7D), suggesting impairment of spermatogenesis and sperm motility caused by SARS-CoV-2 infection.

Effect of coronary heart disease on SARS-CoV-2 infection
Coronary heart disease (CHD) has been reported as a high-risk factor for the COVID-19 associated mortality (Guan et al., 2020). We then compared the lung proteomes of COVID-19 patients with and without CHD, and identified 77 upregulated proteins and 309 downregulated ones (Table S4). The upregulated proteins are enriched in pathways including retinoic acid-inducible gene I (RIG-I) signaling pathway and ribosome biogenesis and its assembly ( Figure 7E). Multiple proteins involved in the RIG-I signaling pathway were upregulated, including probable ATP-dependent RNA helicase (DHX58), pumilio homolog 1 (PUM1), LSM14 homolog A (LSM14A), SLP adaptor, and CSKinteracting membrane protein (SCIMP) in the CHD group (Figure 7F). Interestingly, none of these proteins were significantly regulated in the lung of COVID-19 patients without CHD compared with the non-COVID-19 patients ( Figure 7F), suggesting that these changes might be CHD-specific in the COVID-19 patients. Both RIG-I-like receptor (RLR) and LSM14A belong to the PRR family, serving as sentinels for viral RNAs invasion and inducing the production of antiviral and proinflammatory cytokines, exerting both antiviral and tissue-damaging effects (Lee et al., 2019a;Li et al., 2012). On virus entry, the RLR family member protein DHX58 promotes the recognition of dsRNA by RLRs and enhances the production of IFN-beta and other antiviral genes including itself. Accumulation of DHX58 will then inhibit RLR signaling, which forms a negative feedback (Rehwinkel and Gack, 2020). The elevation of PUM1, a suppressor of DHX58, further exposes elegant modulation of the RIG-I pathway in the CHD patients (Liu et al., 2017). LSM14A has been reported to recognize the viral RNA and promote the expression of IFN-beta (Li et al., 2012). SCIMP is a transmembrane adaptor of TLR4, which triggers the production of proinflammatory cytokines, such as IL-6, in the macrophage (Luo et al., 2017). The upregulation of SCIMP and LSM14A suggests enhanced host defense and hyperinflammation in the lung of COVID-19 patients with CHD.

Limitations of study
This study is limited by the sample size. Due to the relatively small number of patients, the analyses of comorbidities await further investigation in larger independent cohorts. Despite the range of age between the COVID-19 patients and control individuals not being fully balanced, we did not find substantial proteomic difference between elder and younger patient groups (Table  S4). Although the non-COVID-19 tissue samples are from surgeries of individuals with certain diseases (Table S1), they were histologically healthy as examined by two independent senior pathologists ( Figure S1). In addition, future in-depth investigation of the perturbed pathways and the nominated therapeutics is needed.

Conclusions
In summary, we have quantified 11,394 proteins in seven types of organs from patients that died from COVID-19 and identified 5,336 significantly dysregulated proteins compared to non-COVID-19 patients. This proteomic atlas uncovered multiple biological and pathological processes regulated in COVID-19, (B) A hypothetical systems view of the multiple organs' responses to SARS-CoV-2 infection. In the lung, the virus and its released RNA could induce immune response and hijack the host translation mechanism. The innate and adaptive immune cells in the spleen and the cytokine induce acute phase proteins secreted by hepatic cells in response to antiviral defense. Such hyperinflammatory status across the whole body through circulatory system leads to multi-organ injuries. Red boxes, upregulated proteins/pathways; green boxes, downregulated proteins/pathways; blue boxes, pathological processes. See also Figures S6 and S7. ll which include, but not limited to, immune response, protein translation, coagulation disorder, angiogenesis, and profibrotic process. Crosstalk among multiple organs further linked the aforementioned processes by the hyperinflammatory environment with tissue hypoxia after SARS-CoV-2 infection. This systematic proteomic investigation provides a rich resource for improving our understanding of the molecular pathogenesis of SARS-CoV-2 infection and offers clues for therapeutics.

STAR+METHODS
Detailed methods are provided in the online version of this paper and include the following:  (2020).

AUTHOR CONTRIBUTIONS
T  (F) Dysregulated proteins in RIG-I signaling pathway in the lung. The y axis stands for the protein expression ratio by TMT-based quantitative proteomics. Pairwise comparison of each protein among the non-COVID-19, COVID-19 patients with CHD, and without CHD groups was performed using Student's t test. The cutoff of dysregulated proteins has been set at B-H adjusted p value <0.05 and |log2(FC)| > log2(1.2). *p < 0.05; **p < 0.01; ***p < 0.001. See also Table S6.

RESOURCE AVAILABILITY
Lead contact Further information should be directed to and will be fulfilled by the Lead Contact Tiannan Guo (guotiannan@westlake.edu.cn).

Materials availability
This study did not generate new unique reagents.
Data and code availability All data are available in the manuscript or the supplementary material. The accession number for the proteomics data reported in this paper is iProX: IPX0002393000.  (National Health Commission, 2020. All the cases were from Wuhan city and admitted to Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, according to their respiratory symptoms, chest imaging and SARS-CoV-2 nucleic acid test of oropharyngeal swabs. With the consent by the families of the late COVID-19 patients, ultrasound-guided percutaneous multi-point postmortem core biopsies were performed within two hours after death in a negative pressure isolation ward, from February to April 2020. Ultrasound examination was performed using an ultrasound scanner (Mindray M9, Shenzhen Mindray Bio-Medical Electronics Co., Ltd, Guangdong, China) equipped with a C5-1s/SP5-1s/L12-4s probe. All autopsies were performed using the large core (14-guage) needle by the multitask team consisting of pathologists, thoracic surgeons, and ultrasound doctors. Procured tissues were preserved in the 10% neutral formalin immediately after the procedure and fixed for over 24 hours, and then routinely processed for paraffin blocks. Histological slices were prepared using a microtome and stained with hematoxylin and eosin (H&E) as described previously . Immunohistochemical stains (IHC) were performed for CD3 (Dako, Copenhagen, Denmark), CD4 (Dako), CD8 (Dako), CD20 (Roche, Tucson, AZ, USA), CD68 (Dako) and CD163 (Cell Marque, California, USA) according to the manufacturers' protocols on a Dako Link 48 automated stainer (for CD3, CD4, CD8, CD68 and CD163) or a Roche Benchmark XT Ultra system (for CD20). All slides were examined by at least two senior pathologists independently. A total of 144 microscopy-guided dissection samples were analyzed from 73 specimens, including heart (9 specimens), lung (15 specimens), liver (10 specimens), kidney (10 specimens), spleen (9 specimens), testis (5 specimens), and thyroid (15 specimens), as shown in Tables S1 and S2.

EXPERIMENTAL MODEL AND SUBJECT DETAILS
We also collected 74 control samples from archived blocks of 56 non-COVID-19 patients (Tables S1 and S2). Control samples in our study were selected from patients who had undergone surgery in the Union Hospital, Tongji Medical College, Huazhong University of Science and Technology from November 2019 to March 2020. All control samples are histologically healthy tissue samples from the non-COVID-19 patients who were mostly victims of injuries. In case of lack of injured patients, benign lesions from cancer patients were procured. We only selected the tissue samples without bleeding, degeneration or any other pathological lesions. The appropriateness of these control samples has been confirmed by at least two senior pathologists independently, and their representative microscopic pictures are shown in Figure S1. More details of the patients including age, gender and other potential confounding clinical characteristics are provided in Tables S1 and S5.

METHOD DETAILS
Proteomics data acquisition Around 1~1.5 mg FFPE tissue samples were processed to generate peptide samples using accelerated Pressure Cycling Technology (PCT) assisted sample preparation method as described previously Zhu et al., 2019). Briefly, the samples were dewaxed by heptane and rehydrated by ethanol solution of different concentrations; then undergone acidic hydrolysis by 0.1% formic acid, and basic hydrolysis by 0.1 M Tris-HCl (pH 10.0); followed by lysis with 6 M urea/2M thiourea (Sigma) buffer, reduction by Tris(2 carboxyethyl)phosphine (TCEP, Sigma) and alkylation by iodoacetamide (IAA, Sigma) in PCT. Then the lysates were digested using PCT by a mix of Lys-C and trypsin (Hualishi Tech. Ltd, Beijing, China) and the digestion was quenched by trifluoroacetic acid. After C18 cleaning, 7 mg of peptides from each of the 218 tissue samples, and 37 technical replicates, 34 pooling samples including common pooled samples and tissue specific pooled samples (Table S2) were taken for the following TMTpro 16plex labeling based proteomics analysis. The common pooled samples were prepared by mixing equal amount of peptides from six different human organs including testis, lung, kidney, spleen, heart and liver. The tissue specific pooled samples were prepared by mixing equal amount of peptides from all samples of a specific organ. And then all these peptide samples were labeled with TMTpro 16plex reagent Figure S1. Comparison of histopathological features of organs between COVID-19 and non-COVID-19 patients (H&E 3200), related to Figure 1 The left column shows the results of the non-COVID-19 patients and the right part shows the results of the COVID-19 patients. A. Pathological features of lungs from the COVID-19 patients. (i) The lung showed diffuse alveolar damage, with the alveolar epithelia replaced by hyperplastic type II alveolar epithelia and falling of type II alveolar epithelia (a). Alveolar septa were thickened with the proliferation of fibroblasts/myofibroblasts and fibrosis (b). Fibrinous exudation and hyaline membrane formation (c) were observed, with neutrophils aggregation in the alveolar cavity (d). (ii) There were some microthrombus in small vessels of the lung (a), and megakaryocytes in the alveolar septal capillaries (b). B. Pathological findings of the spleen from COVID-19 patients. (i) The red pulp of the spleen expanded, and splenic sinus extended with hyperemia accompanying with macrophage proliferation. (ii) The white pulp of spleen exhibited atrophy with significantly reduced lymphocytes. C. Pathological findings of the liver from COVID-19 patients. (i) Coagulative necrosis of hepatocytes was observed in zone III. (ii) The hepatocytes exhibited prominent steatosis. D. Pathological findings of the heart from COVID-19 patients. (i) Atrophic myocardia and scant lymphocytes (a) were present in the edematous cardiac interstitium. (ii) The myocardium showed hydropic degenerative change. E. Pathological findings of the kidney from the COVID-19 patients. (i) In the renal cortex, the proximal tubules showed prominent acute tubule injury manifested as the loss of brush border and epithelial cells necrosis in the tubular lumens (a). Microthrombi in peritubular capillaries and glomeruli (b) were frequently observed. (ii) In the renal medulla, the collecting ducts showed occasional cellular swelling and atrophy (a), and edema without significant inflammation (b). F. Pathological findings of the thyroid from the COVID-19 patients. Lymphoid infiltration was found in some of interfollicular region. Neither neutrophilic infiltration nor necrosis was present.