- Open Access
Transcriptional landscape of circulating platelets from patients with COVID-19 reveals key subnetworks and regulators underlying SARS-CoV-2 infection: implications for immunothrombosis
Cell & Bioscience volume 12, Article number: 15 (2022)
Thrombosis and coagulopathy are pervasive pathological features of coronavirus disease 2019 (COVID-19), and thrombotic complications are a sign of severe COVID-19 disease and are associated with multiple organ failure and increased mortality. Platelets are essential cells that regulate hemostasis, thrombus formation and inflammation; however, the mechanism underlying the interaction between platelets and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) remains unclear.
The present study performed RNA sequencing on the RNA isolated from platelets obtained from 10 COVID-19 patients and eight healthy donors, and discovered that SARS-CoV-2 not only significantly altered the coding and non-coding transcriptional landscape, but also altered the function of the platelets, promoted thrombus formation and affected energy metabolism of platelets. Integrative network biology analysis identified four key subnetworks and 16 risk regulators underlying SARS-CoV-2 infection, involved in coronavirus disease-COVID-19, platelet activation and immune response pathways. Furthermore, four risk genes (upstream binding transcription factor, RNA polymerase II, I and III subunit L, Y-box binding protein 1 and yippee like 2) were found to be associated with COVID-19 severity. Finally, a significant alteration in the von Willebrand factor/glycoprotein Ib-IX-V axis was revealed to be strongly associated with platelet aggregation and immunothrombosis.
The transcriptional landscape and the identification of critical subnetworks and risk genes of platelets provided novel insights into the molecular mechanisms of immunothrombosis in COVID-19 progression, which may pave the way for the development of novel therapeutic strategies for preventing COVID-19-associated thrombosis and improving the clinical outcome of COVID-19 patients.
Coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)  and represents one of the greatest public health challenges since the 1918 influenza pandemic over 100 years ago . The main epidemiological characteristics of COVID-19 are rapid spread, high infectiousness and an approximately 1–3% mortality rate [3,4,5,6]. Due to the lack of specific treatment options available [7, 8], symptomatic treatment and supportive care remain the primary treatment strategy.
The respiratory symptoms associated with COVID-19 can lead to acute respiratory distress syndrome, in severe cases . In addition, accumulating evidence has suggested that thrombosis and coagulopathy are pervasive pathological features of COVID-19 [10, 11]. Compared with other respiratory infectious diseases, COVID-19 shows a higher cumulative incidence of thrombotic complications [12,13,14], especially in COVID-19 patients admitted to intensive care units (ICUs) [15, 16]. It is noteworthy that hospitalized children with COVID-19 have been reported to be at an increased risk of developing thrombosis [17, 18]. In addition, the thrombotic complications of COVID-19 patients are associated with multiple organ failure and increased mortality [19, 20].
Autopsies of COVID-19 patients have shown that platelet-rich thrombi exist in the microcapillaries of multiple organs. The number of megakaryocytes in the heart and lungs is abnormally increased [21,22,23]. Previous studies have observed that the platelets of COVID-19 patients appear to be overreactive and may interact with SARS-CoV-2 to promote coagulation dysfunction during COVID-19 infection [24,25,26]. Platelets, the second most abundant type of cell in the peripheral blood , are best known as mediators of thrombus formation and hemostasis . In recent years, studies have reported that platelets were the key sentinel and effector cells in infectious diseases, such as Dengue virus (DENV) and malaria [29, 30]. During pathogen invasion, the transcriptome and proteome of platelets have been reported to be altered to augment host defense mechanisms [31, 32]; however, these changes may also result in adverse outcomes. Activated platelets release various cytokines, including chemokines CXCL1, PF (platelet factor)-4, to amplify thrombin production, enhance leukocyte recruitment, promote neutrophil extracellular trap formation, upregulate the endothelial expression of proinflammatory cytokines and, finally, induce immunothrombosis . Nevertheless, the mechanism underlying the interaction between platelets and SARS-CoV-2 remains unknown to the best of our knowledge.
The present study performed RNA sequencing (RNA-seq) on RNA isolated from platelets obtained from 10 COVID-19 patients and eight healthy donors and further examined the transcriptional dysregulation in circulating platelets from COVID-19 patients. Then, integrative network biology analysis was performed to identify critical subnetworks and risk regulators underlying SARS-CoV-2 infection. Our findings may provide novel insights into the molecular mechanisms of immunothrombosis in COVID-19 progression and pave the way for developing therapeutic strategies for preventing COVID-19-associated thrombosis and improving the clinical outcome of COVID-19 patients.
Clinical characteristics of healthy donors and hospitalized patients with COVID-19
As shown in Table 1, COVID-19 patients were matched with healthy donors by age and sex. Comorbidities, including diabetes, hypertension and cancer, were present in 40% of the COVID-19 cases, which is consistent with previous reports [33,34,35]. The symptoms of the patient’s initial diagnosis mainly include fever (80%), cough (80%) and sputum (50%), occasionally sore throat (30%), fatigue (10%), muscle aches (20%), and nausea (10%). And the respiratory rate and heart rate of the hospitalized patients were in the normal range. Additionally, in Additional file 1: Fig. S1, the imaging of all patients was consistent with the changes of viral pneumonia. The CRP levels were significantly increased in the COVID-19 patient group, whereas the lymphocyte count was notably decreased, which is consistent with previous reports that found that elevated CRP levels and lymphopenia were the main laboratory characteristics of COVID-19 patients [36,37,38,39]. As shown in Fig. 1, the platelet count, plateletcrit (PCT), mean platelet volume (MPV), platelet distribution (PDW), prothrombin time (PT), and activated partial thromboplastin time (APTT) of the COVID-19 patient group were either within the normal ranges or only slightly exceeded the normal ranges, but the platelet count was substantially lower, and the PT and APTT were significantly increased compared with the control group. Furthermore, the alanine aminotransferase and creatinine levels between the COVID-19 patients and the healthy donors were not significantly different and were all within the normal range, indicating that liver and kidney damage or failure was not present in the COVID-19 patient cohort. The above results support that all patients enrolled in the group are ordinary new coronavirus patients.
SARS-CoV-2 infection alters the coding and non-coding transcriptional landscape of human platelets
To determine whether SARS-CoV-2 infection altered the transcriptome of circulating human platelets, we performed RNA-seq on RNA isolated from platelets obtained from 10 COVID-19 patients and eight healthy donors. Hierarchical clustering of transcriptome-wide RNA expression showed differential grouping of COVID-19 patients and healthy donors (Additional file 2: Figure S2), suggesting that SARS-CoV-2 infection altered the transcriptional landscape of platelets. Differential expression analysis identified 1,223 differentially expressed RNAs between COVID-19 patients and healthy donors (adjusted P-value < 0.05 and absolute log2FC > 1.0), including 191 and 13 significantly upregulated mRNAs and lncRNAs, respectively, and 883 and 136 significantly downregulated mRNAs and lncRNAs, respectively (Fig. 2 and Additional file 5: Table S1).
We performed GO functional term and KEGG signaling pathway enrichment analyses to determine the involvement of differentially expressed RNAs in biological processes. GO functional enrichment analysis revealed that differentially expressed mRNAs and lncRNAs were significantly enriched in biological processes involved in ‘hemostasis’, ‘platelet activation’, ‘immune response’ and ‘metabolic process and energy’ (Fig. 3A, B). The top 20 enriched KEGG signaling pathways are shown in Fig. 3C, in which differentially expressed mRNAs and lncRNAs were highly clustered in several signaling pathways associated with ‘coronavirus disease-COVID-19’, ‘platelet activation’ and ‘immune’. In addition, differentially expressed mRNAs and lncRNAs were observed to be preferentially expressed in the blood, spleen and bone marrow, as determined using tissue and cell-specific enrichment analysis (Fig. 3D). These data suggested that SARS-CoV-2 infection may alter the expression of platelet RNAs, change the function of platelets, promote thrombosis and affect the biological processes of energy metabolism.
Integrative network biology analysis of circulating platelets from COVID-19 patients
To gain further insight into the biological network of transcriptional changes associated with SARS-CoV-2 infection and platelet activation, we constructed a dysGCN based on differentially expressed mRNAs and lncRNAs by calculating the PCC. The dysGCN was built and included 960 nodes (118 lncRNAs and 842 mRNAs) and 30,716 connections between them (Fig. 4A). After module mining, we identified a total of 17 co-expression modules ranging in size from 6 to 130 RNAs. Functional characterization using GO functional term and KEGG signaling pathway enrichment analyses showed that modules-2, -8 and -16 were enriched in ‘coronavirus disease-COVID-19’ and ‘platelet activation pathways’. In addition, module-2 was also associated with ‘immune response pathways’ and module-15 was enriched in the ‘coronavirus disease-COVID-19 pathway’ only (Fig. 4B, C). These results suggested that modules-2, -8, -15 and -16 are key subnetworks involved in circulating platelets during SARS-CoV-2 infection.
Systematic identification of key regulators in circulating platelets during SARS-CoV-2 infection
The four modules formed a key subnetwork, consisting of 196 nodes (17 lncRNAs and 179 mRNAs) and 8,822 edges between them, and were mainly involved in the biological process and RNA interaction of platelet activation caused by SARS-CoV-2 infection. To discover critical regulators of this subnetwork during SARS-CoV-2 infection, we mapped these mRNAs to ‘coronavirus disease-COVID-19’ and ‘platelet activation pathways’, and identified 16 key candidate regulators (Fig. 5A). These 16 key candidate regulators were predominantly expressed in the blood and spleen, as shown in tissue and cell-specific enrichment analysis (Additional file 3: Figure S3A). Further analysis using disease and viral perturbation datasets from the GEO database also showed that these key candidate regulators were associated with blood coagulation-related diseases (Additional file 3: Figure S3B) and were differentially expressed in SARS-CoV-2- or SARS virus-associated GEO datasets (Additional file 3: Figure S3C).
Furthermore, we calculated the expression correlation of these 16 key candidate regulators in healthy donors and disease samples. As shown in Fig. 5B, the expression levels of these 16 key candidate regulators were significantly correlated in COVID-19 patients, which were not observed in healthy donors. These observations indicated that these regulators were differentially expressed and exhibited a highly synergistic pattern during SARS-CoV-2 infection. We observed a significant overlap between the co-expression relationship and protein–protein interactions (PPIs) among the 16 key candidate regulators when considering the physical interaction.
Finally, we combined the results obtained from the enrichment analyses, identification of key subnetworks and regulators, PPI network and previous studies to infer the key active pathway involved in platelet aggregation during SARS-CoV-2 infection. As shown in Fig. 5D, von Willebrand factor (VWF) bound to the platelet membrane glycoprotein (GP) Ib-IX-V complex and initiated a signaling cascade, which subsequently activated LYN proto-oncogene, Src family tyrosine kinase (LYN). The downstream pathway was initiated by LYN and propagated through PI3K, phospholipase C γ2 and protein kinase C, which promoted ATP secretion and finally caused αIIbβ3-mediated platelet activation and aggregation.
We measured the levels of four DEmRNAs (F13A1, ITGB3, ITGA2B and VWF) involved in platelet aggregation during SARS-CoV-2 infection in 5 patients with pneumonia and 6 normal donors by RT-qPCR. As shown in Additional file 4: Figure S4, we found that the ITGA2B and VWF were significantly decreased in patients with pneumonia compared with control donors. In addition, the F13A1 and ITGB3 also showed lower expression levels in patients with pneumonia compared with control donors although the difference did not reach statistical significance which might be due to the limitation of sample size. These RT-qPCR results are in line with our observation form RNA-seq.
Identification of risk regulators associated with COVID-19 severity
COVID-19 patients are often accompanied by many complications, and the development of these complications often contributes to the progression of the disease and even death [40, 41]. We used RNAs identified in the key subnetwork to construct an RNA-disease association network to investigate how COVID-19-associated transcriptional characteristics of circulating platelets were shared with other diseases. Information related to disorders, genes, and possible association was obtained from the Online Mendelian Inheritance in Man (OMIM) and CTD databases. We classified the associated disorders into 19 categories. As shown in Fig. 6A, the number of interconnected genes with diseases suggested that key subnetworks in circulating platelets underlying SARS-CoV-2 infection were closely associated with immunological, cardiovascular, liver and hematological diseases, and cancer.
To further examine whether RNAs in the key subnetwork were risk regulators associated with COVID-19 severity, we compared the expression levels of RNAs in the key subnetwork between non-ICU COVID-19 patients and healthy individuals [mild (M)—health (H)], ICU COVID-19 patients and healthy individuals [severe (S)—H], and ICU COVID-19 patients and non-ICU COVID-19 patients (S—M) in both our and Bhanu’s cohorts. According to their expression levels, we identified four RNAs as candidate risk regulators associated with COVID-19 severity (Fig. 6B). As shown in Fig. 6C, RNA polymerase II, I and III subunit L (POLR2L), Y-box binding protein 1 (YBX1) and yippee like 2 (YPEL2) were revealed to be upregulated, while upstream binding transcription factor (UBTF) was downregulated.
COVID-19 patient platelet-specific transcriptional characteristics can be used to identify potential drugs for therapeutic repurposing
We used core mRNAs to identify potential drugs and antiviral agents to analyze their interactions with different drugs and constructed a gene-drug network (Fig. 7A). We identified 54 drugs that acted against the 16 key regulators, implying that these key regulators have been widely used as drug targets to treat different diseases, and genes that have no connection with drugs in the network may be potential drug targets. In addition, we replaced these drugs with their ATC classification and constructed a gene-drug ATC coding network (Fig. 7B) to gain further information about the drugs. Most of the drugs were used for blood and blood-forming organs and acted as antineoplastic and immunomodulating agents. Since these genes were identified as potential drug targets, we inputted these genes into the CMAP database to predict the candidate chemicals. The top 10 drugs associated with the downregulated expression of mRNAs in COVID-19 patients are shown in Table 2, amongst which gabexate, valproic acid, estradiol, ketorolac and pronetalol have already been reported in previous studies. As these signature drugs were detected to target the key regulators, these drugs may represent potential drugs for COVID-19 and stimulate platelet activation.
Since the outbreak of COVID-19, an increasing number of studies have reported the hyperactivity of activated platelets in COVID-19 patients [26, 42, 43], and their role in immunothrombosis induced by SARS-CoV-2 [20, 44]. Emerging evidence has shown that biomarkers of platelet activation were significantly increased in COVID-19 patients and were associated with thrombosis and increased mortality in hospitalized COVID-19 patients, including CD40 ligand (sCD40L), P-selectin, the metabolite of thromboxane A2, thromboxane B2 (TxB2) [45, 46]. Thus, it remains an important research priority to clarify the mechanisms underlying the interaction between platelets and SARS-CoV-2 to provide potential treatment strategies for immunothrombosis during COVID-19.
We performed RNA-seq on RNA isolated from platelets obtained from 10 COVID-19 patients and eight healthy donors in the present study. We found that SARS-CoV-2 significantly altered the coding and non-coding transcriptional landscape (Fig. 2 and Additional file 1: Figure S1). However, the platelets and coagulation indicators of COVID-19 patients were within the normal range or only slightly exceeded the normal range (Fig. 1 and Table 1). We further performed functional and tissue enrichment analysis and discovered that SARS-CoV-2 changed the function of the platelets, promoted thrombosis and affected their energy metabolism. However, the activation and hyperresponsiveness of platelets in patients with infectious diseases are known to be caused by numerous factors. For example, during excessive inflammation, the accumulation of proinflammatory cytokines and chemokines has been found to activate platelets indirectly . Viruses, such as DENV, influenza virus  and human immunodeficiency virus , have also been identified to activate platelets directly. Therefore, the current study conducted further analysis to explore how SARS-CoV-2 may activate and regulate platelets.
We constructed a dysregulated gene co-expression network based on differentially expressed mRNAs and lncRNAs and identified four key subnetworks in the circulating platelets during SARS-CoV-2 infection. The four key subnetworks were found to be involved in ‘coronavirus disease-COVID-19’, ‘platelet activation’ and ‘immune response pathways’. These results indicated that the four key subnetworks might contain the key factor that SARS-CoV-2 infection regulates to promote platelet activation. We subsequently mapped these mRNAs to ‘coronavirus disease-COVID-19’ and ‘platelet activation pathways’ and identified 16 key candidate regulators (Fig. 5A). In previous studies, Palma et al. found that PI3Kδ inhibition could be used as a potential therapeutic target for COVID-19 to reduce inflammation and patient death . A single-center cross-sectional study report showed that VWF activity and antigen levels continued to increase as the disease progressed , which was associated with the poor prognosis of COVID-19 patients . Another comparative study discovered that integrin and integrin signaling genes, including integrin subunit β 3 and integrin subunit α 2b, were highly expressed in lung samples and were suggested to provide SARS-CoV-2 with a more effective competitive advantage for invading lung cells . The previous published experimental work provides further evidence for the core role of the 16 key candidate regulators in the regulation of platelets following SARS-CoV-2 infection.
Notably, we observed a significant overlap between the co-expression relationship and PPI among the 16 key candidate regulators. Thus, we combined the results obtained from the enrichment analyses, identification of key subnetworks and regulators, the PPI network and reported studies to hypothesize that a significant alteration in the VWF/GP Ib-IX-V axis may be strongly associated with platelet aggregation and immunothrombosis. It is well-established that the entry of SARS-CoV-2 into cells relies on angiotensin-converting enzyme 2 (ACE2) and transmembrane protease serine 2 (TMPRSS2) for viral spike protein initiation, and this receptor is highly expressed by nasopharyngeal airway and alveolar epithelial cells, vascular endothelial cells and lung macrophages . However, an increased number of studies have suggested that platelets do not express ACE2 and TMPRSS2 and the activation and regulation of platelets by SARS-CoV-2 seems to occur independently of ACE and TMPRSS2 [26, 42, 54]. VWF is a multimeric glycoprotein in the plasma, which plays a vital role in hemostasis and thrombosis by mediating the adhesion of platelets to damaged and activated blood vessels . During infection, inflammation may regulate the formation of VWF-platelet thrombi in large arteries and small blood vessels by increasing VWF levels, enhancing VWF responsiveness and regulating the level and activity of regulatory molecules in the circulation . This explains the significant changes in platelet VWF/GP Ib-IX-V levels in patients with COVID-19 in the present study, and further indicates that the cascade initiated by the VWF/GP Ib-IX-V complex may play a pivotal role in the activation and regulation of platelets by SARS-CoV-2 infection.
Furthermore, the course of COVID-19 develops rapidly and can lead to severe and fatal complications, such as acute myocardial injury and chronic damage to the cardiovascular system [40, 41]. To facilitate the stratification of high-risk patients, reliable biomarkers related to disease progression are urgently required in the clinic. In the present study, four risk regulators (UBTF, POLR2L, YBX1 and YPEL2) were identified to be associated with COVID-19 severity. Therefore, these four regulators may represent potential biomarkers for the clinical management of COVID-19 patients.
Immunothrombosis is strongly associated with the progression of COVID-19 disease. Timely and accurate anticoagulation and vigilant monitoring of thrombotic complications are key to COVID-19 patient management . The current study identified 54 drugs targeting the 16 key regulators and five out of the top 10 drugs had been previously reported in the literature. Thus, they may represent potential drugs for the treatment of COVID-19-associated thrombosis.
The findings of the present study provided novel evidence to suggest that SARS-CoV-2 may significantly alter the coding and non-coding transcriptional landscape of platelets and subsequently alter their function. The altered transcriptional landscape and the identification of the critical subnetwork and risk regulators of platelets provided novel insights into the molecular mechanisms of immunothrombosis in COVID-19 progression, which paves the way for developing new therapeutic strategies for preventing COVID-19-associated thrombosis and improving the clinical outcome of COVID-19 patients.
Subjects and specimen collection
The current study was approved by the Institutional Review Board (IRB) of The Second Affiliated Hospital of Wenzhou Medical University (IRB# LCKY2020-09 and LCKY2020-222), and written informed consent was obtained from each participant prior to the commencement of the study. A total of 10 patients with COVID-19 were recruited from The Second Affiliated Hospital of Wenzhou Medical University between 6 and 16 February 2020, and the control group consisting of eight donors were recruited between 4 May and 15 August 2020. COVID-19 cases were confirmed by detecting SARS-CoV-2 RNA using reverse transcription-quantitative PCR. All donors of the control group underwent preliminary screening and the following exclusion criteria were applied: (i) Use of antiplatelet drugs; (ii) recent infections; and (iii) diagnosis of myelosuppressive diseases. K2-EDTA-anticoagulated whole blood was collected from hospitalized COVID-19 patients before antiplatelet or anticoagulation therapy. In addition, the blood count and C reactive protein (CRP) were detected by Sysmex XE 5000 hematology analyzer (Sysmex Corporation, Kobe, Japan), the coagulation parameters were detected using an automation coagulator STA-R Evolution (Stago, NJ, USA), and biochemical Indicators were measured using Cobas 6000 c501 (Roche Diagnostics, Rotkreuz, Switzerland). The clinical data of all participants, including age, sex, dietary history, medication use and laboratory parameters, were collected and are shown in Table 1.
Another platelet RNA-seq cohort of 10 COVID-19 patients (including six non-ICU and four ICU COVID-19 patients) was retrieved from NCBI short-read archives under PRJNA634489 .
Platelet isolation and sequencing
The samples of all hospitalized patients were taken from the first morning after the patients were admitted to the hospital. The same standardized operating procedure was used to collect the whole blood from COVID-19 patients and healthy donors. Platelet parameters were measured using a Sysmex XE 5000 hematology analyzer (Sysmex Corporation, Kobe, Japan). Platelets were isolated from the whole blood within 12 h and leukocyte-removed platelets were obtained using the previously described method [59,60,61]. Briefly, whole blood collected in the 5 ml EDTA-coated vacuum tube was centrifuged at 120×g for 20 min to separate the platelet-rich plasma (PRP) from the nucleated blood cells. The PRP was then centrifuged using a 20-min 360×g step to pellet the platelets. Subsequently, 9/10th of the upper PRP was carefully removed and the platelet pellet was collected and resuspended in RNAlater, then frozen at − 80 °C. The counts and the purity of the separated platelets were confirmed using flow cytometry.
For next-generation RNA-seq, total RNA was isolated from 1 × 109 separated platelets using TRIzol® reagent (Invitrogen; Thermo Fisher Scientific, Inc., Waltham, MA, USA), followed by DNase treatment. The quantity and quality of total RNA were determined and used for RNA library construction.
Quality control of raw RNA-seq data was conducted using Trimmomatic (v0.39) , including read trimming and quality filtering. Clean RNA-seq reads were mapped to the human reference genome (hg19) that included GENCODE v29 gene annotation using the HISAT2 program (v2.2.1, with default parameters). Alignment files were processed and the final BAM files were obtained using SAMtools (V1.7) . The number of reads mapping each human gene (as annotated in the GENCODE v29 gene annotation) was counted using HTSeq v0.11.2 . Differential expression analysis was conducted using the DESeq2 package , and statistical significances in expression between COVID-19 patients and healthy individuals were determined using a false discovery rate-adjusted P-value < 0.05 and a log2fold change (FC) > 1.0.
Gene enrichment analyses
Functional enrichment analyses of differentially expressed RNAs, including Gene Ontology (GO) functional term and Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling pathway enrichment analyses, and tissue and cell-specific enrichment analysis , were performed and visualized using the Metascape gene enrichment analysis tool (https://metascape.org)  and R package ‘clusterProfiler’ . Enrichment analysis of differentially expressed RNAs in virus perturbations datasets obtained from the Gene Expression Omnibus (GEO) database was conducted using Enrichr (https://amp.pharm.mssm.edu/Enrichr) .
Co-expression network analysis
The Pearson’s correlation coefficient (PCC) was calculated to measure the linear expression correlation between each pair of differentially expressed RNAs. A dysregulated gene co-expression network (dysGCN), which was composed of differentially expressed protein-coding genes (PCGs) and long non-coding RNAs (lncRNAs) as points and a PCC of > 0.8 as edges, was then constructed and visualized using Cytoscape (http://www.cytoscape.org). The densely connected and bipartite network modules were identified using Molecular Complex Detection (http://apps.cytoscape.org/apps/mcodeplugin) of Cytoscape software .
Drug prediction analysis
Proteins are the direct targets of drugs and the bridge connecting genes and targeted drugs; thus, we mapped the RNA to protein using the UniProt database (https://www.uniprot.org). The UniProt ID was input into the DrugBank database (https://go.drugbank.com) to identify drugs related to the proteins. To obtain more information about the drugs, we also recorded the Anatomical Therapeutic Chemical (ATC) classification of the drugs and constructed a gene-drug ATC coding network. Since these genes can be used as disease targets, some candidate chemicals have been identified to predict the chemicals that target these genes and are related to the treatment of COVID-19 and bleeding. The Connectivity Map (CMAP) database  uses gene-expression signatures to predict small molecular compounds for a specific disease. Upregulated and downregulated genes were subsequently uploaded to the CMAP database and small molecule drugs that may treat COVID-19 and bleeding were searched for.
Real-Time Quantitative PCR (RT-qPCR)
A total of 5 patients with pneumonia and 6 control donors were recruited from The Second Affiliated Hospital of Wenzhou Medical University on 17 December 2021. The samples of patients were taken from the first morning after the patients were admitted to the hospital and the samples of control donors were collected on the morning of the physical examination. Platelet total RNA from all samples was extracted according to the method described above within 12 h. Total RNA was reversely transcribed into cDNA using the reverse transcriptase reaction kit (vazyme). DNA removal was conducted at 42 °C for 2 min with the addition of 4ul 4XgDNA wiper mix and 12ul RNA, and then reverse transcription reaction was conducted at 50 °C for 15 min and 85 °C for 5 s with the addition of 4ul 5X hiscript qRT supermix. RT-qPCR was carried out on a Bio-Rad C1000 quantitative PCR instrument with the SYBR Select mastermix (Applied Biosystems) kits. The amplification conditions were: 2 min at 50 °C, and 2 min at 95 °C; 15 s at 95 °C, and 1 min at 60 °C (40cycles). Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as the internal reference for the mRNAs. The ΔCq was used to calculate the expression of the target gene in the patient's group and the control group. The primer sequences were shown in Additional file 6: Table S2. The primers were synthesized by Sangon Biotech (Shanghai, China).
Availability of data and materials
The datasets used and/or analyzed during the present study are available from the corresponding author on reasonable request.
Coronavirus disease 2019
Severe acute respiratory syndrome coronavirus 2
Intensive care units
C reactive protein
Gene Expression Omnibus
Pearson’s correlation coefficient
Activated partial thromboplastin time
Von Willebrand factor
Online Mendelian Inheritance in Man
Angiotensin-converting enzyme 2
Transmembrane protease serine 2
Wan Y, Shang J, Graham R, Baric RS, Li F. Receptor Recognition by the Novel Coronavirus from Wuhan: an analysis based on decade-long structural studies of sars coronavirus. J Virol. 2020;94:7.
Tncpere T. The Epidemiological characteristics of an outbreak of 2019 novel Coronavirus diseases (COVID-19)—China, 2020. CCDCW. 2020;2019:113–22.
Zhou P, Yang XL, Wang XG, Hu B, Zhang L, Zhang W, et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature. 2020;579(7798):270–3.
Wu ZMJ. Characteristics of and Important Lessons From the Coronavirus Disease 2019 (COVID-19) Outbreak in China: Summary of a Report of 72 314 Cases From the Chinese Center for Disease Control and Prevention. JAMA. 2020;323(13):1239–42.
Ruan Q, Yang K, Wang W, Jiang L, Song J. Clinical predictors of mortality due to COVID-19 based on an analysis of data of 150 patients from Wuhan. China Intensive Care Med. 2020;46(5):846–8.
Tian S, Hu N, Lou J, Chen K, Kang X, Xiang Z, et al. Characteristics of COVID-19 infection in Beijing. J Infect. 2020;80(4):401–6.
Scavone CBS, Bertini M, Sportiello L, Rafaniello C, Zoccoli A, Berrino L, Racagni G, Rossi F, Capuano A. Current pharmacological treatments for COVID-19: What’s next? Br J Pharmacol. 2020;177(21):4813–24.
Zhao Z, Bourne PE. Structural Insights into the Binding Modes of Viral RNA-Dependent RNA polymerases using a function-site interaction fingerprint method for RNA Virus Drug Discovery. J Proteome Res. 2020;19(11):4698–705.
Shi R, Lai C, Teboul JL, Dres M, Moretto F, De Vita N, et al. COVID-19 ARDS is characterized by higher extravascular lung water than non-COVID-19 ARDS: the PiCCOVID study. Crit Care. 2021;25(1):186.
Matsuyama T, Kubli SP, Yoshinaga SK, Pfeffer K, Mak TW. An aberrant STAT pathway is central to COVID-19. Cell Death Differ. 2020;27(12):3209–25.
Sriram K, Insel PA. Inflammation and thrombosis in COVID-19 pathophysiology: proteinase-activated and purinergic receptors as drivers and candidate therapeutic targets. Physiol Rev. 2021;101(2):545–67.
Ackermann M, Verleden SE, Kuehnel M, Haverich A, Welte T, Laenger F, et al. Pulmonary vascular endothelialitis, thrombosis, and angiogenesis in covid-19. N Engl J Med. 2020;383(2):120–8.
Iba T, Levy JH, Levi M, Connors JM, Thachil J. Coagulopathy of coronavirus disease 2019. Crit Care Med. 2020;48(9):1358–64.
Chiumello D, Busana M, Coppola S, Romitti F, Formenti P, Bonifazi M, et al. Physiological and quantitative CT-scan characterization of COVID-19 and typical ARDS: a matched cohort study. Intensive Care Med. 2020;46(12):2187–96.
Klok FA, Kruip M, van der Meer NJM, Arbous MS, Gommers D, Kant KM, et al. Confirmation of the high cumulative incidence of thrombotic complications in critically ill ICU patients with COVID-19: An updated analysis. Thromb Res. 2020;191:148–50.
Helms J, Tacquard C, Severac F, Leonard-Lorant I, Ohana M, Delabranche X, et al. High risk of thrombosis in patients with severe SARS-CoV-2 infection: a multicenter prospective cohort study. Intensive Care Med. 2020;46(6):1089–98.
Goldenberg NA, Sochet A, Albisetti M, Biss T, Bonduel M, Jaffray J, et al. Consensus-based clinical recommendations and research priorities for anticoagulant thromboprophylaxis in children hospitalized for COVID-19-related illness. J Thromb Haemost. 2020;18(11):3099–105.
Feldstein LR, Rose EB, Horwitz SM, Collins JP, Newhams MM, Son MBF, et al. Multisystem Inflammatory Syndrome in US Children and Adolescents. N Engl J Med. 2020;383(4):334–46.
Tang N, Li D, Wang X, Sun Z. Abnormal coagulation parameters are associated with poor prognosis in patients with novel coronavirus pneumonia. J Thromb Haemost. 2020;18(4):844–7.
McFadyen JD, Stevens H, Peter K. The Emerging Threat of (Micro)Thrombosis in COVID-19 and Its Therapeutic Implications. Circ Res. 2020;127(4):571–87.
Lax SF, Skok K, Zechner P, Kessler HH, Kaufmann N, Koelblinger C, et al. Pulmonary Arterial Thrombosis in COVID-19 with fatal outcome : results from a prospective, single-center. Clinicopathologic Case Series Ann Intern Med. 2020;173(5):350–61.
Wichmann D, Sperhake JP, Lutgehetmann M, Steurer S, Edler C, Heinemann A, et al. Autopsy findings and venous thromboembolism in patients with COVID-19: a prospective cohort study. Ann Intern Med. 2020;173(4):268–77.
Rapkiewicz AV, Mai X, Carsons SE, Pittaluga S, Kleiner DE, Berger JS, et al. Megakaryocytes and platelet-fibrin thrombi characterize multi-organ thrombosis at autopsy in COVID-19: A case series. EClinicalMedicine. 2020;24:100434.
Shen S, Zhang J, Fang Y, Lu S, Wu J, Zheng X, et al. SARS-CoV-2 interacts with platelets and megakaryocytes via ACE2-independent mechanism. J Hematol Oncol. 2021;14(1):72.
Aid M, Busman-Sahay K, Vidal SJ, Maliga Z, Bondoc S, Starke C, et al. Vascular disease and thrombosis in SARS-CoV-2-infected rhesus macaques. Cell. 2020;183(5):1354–66.
Manne BKDF, Middleton EA, Portier I, Rowley JW, Stubben C, Petrey AC, Tolley ND, Guo L, Cody M, Weyrich AS, Yost CC, Rondina MT, Campbell RA. Platelet gene expression and function in patients with COVID-19. Blood. 2020;136(11):1317–29.
Junt TSH, Chen Z, Massberg S, Goerge T, Krueger A, Wagner DD, Graf T, Italiano JE Jr, Shivdasani RA, von Andrian UH. Dynamic visualization of thrombopoiesis within bone marrow. Science. 2007;317(5845):1767–70.
Zm R. Platelets in atherothrombosis. Nat Med. 2002;8(11):1227–34.
Petrey AC, Obery DR, Kessler SP, Zawerton A, Flamion B, de la Motte CA. Platelet hyaluronidase-2 regulates the early stages of inflammatory disease in colitis. Blood. 2019;134(9):765–75.
Li Guo MTR. The era of thromboinflammation: platelets are dynamic sensors and effector cells during infectious diseases. Front Immunol. 2019;10:89.
Rondina MT, Weyrich AS. Regulation of the genetic code in megakaryocytes and platelets. J Thromb Haemost. 2015;13(Suppl 1):S26-32.
Rowley JW, Oler AJ, Tolley ND, Hunter BN, Low EN, Nix DA, et al. Genome-wide RNA-seq analysis of human and mouse platelet transcriptomes. Blood. 2011;118(14):e101–11.
Guan WJ, Ni ZY, Hu Y, Liang WH, Ou CQ, He JX, et al. Clinical characteristics of coronavirus disease 2019 in China. N Engl J Med. 2020;382(18):1708–20.
Dennis JM, Mateen BA, Sonabend R, Thomas NJ, Patel KA, Hattersley AT, et al. Type 2 Diabetes and COVID-19-related mortality in the critical care setting: a national cohort study in England, March-July 2020. Diabetes Care. 2021;44(1):50–7.
Richardson S, Hirsch JS, Narasimhan M, Crawford JM, McGinn T, Davidson KW, et al. Presenting characteristics, comorbidities, and outcomes among 5700 patients hospitalized With COVID-19 in the New York City Area. JAMA. 2020;323(20):2052–9.
Li S, Jiang L, Li X, Lin F, Wang Y, Li B, et al. Clinical and pathological investigation of patients with severe COVID-19. JCI Insight. 2020;5:12.
Wu C, Chen X, Cai Y, Xia J, Zhou X, Xu S, et al. Risk factors associated with acute respiratory distress syndrome and death in patients with coronavirus disease 2019 Pneumonia in Wuhan. China JAMA Intern Med. 2020;180(7):934–43.
Chen G, Wu D, Guo W, Cao Y, Huang D, Wang H, et al. Clinical and immunological features of severe and moderate coronavirus disease 2019. J Clin Invest. 2020;130(5):2620–9.
Iadecola C, Anrather J, Kamel H. Effects of COVID-19 on the Nervous System. Cell. 2020;183(1):16–27.
Zheng YY, Ma YT, Zhang JY, Xie X. COVID-19 and the cardiovascular system. Nat Rev Cardiol. 2020;17(5):259–60.
Siddiqi HK, Libby P, Ridker PM. COVID-19 - A vascular disease. Trends Cardiovasc Med. 2021;31(1):1–5.
Zaid Y, Puhm F, Allaeys I, Naya A, Oudghiri M, Khalki L, et al. Platelets Can Associate with SARS-Cov-2 RNA and Are Hyperactivated in COVID-19. Circ Res. 2020;9:78.
Koupenova M, Freedman JE. Platelets and COVID-19: Inflammation, Hyperactivation and Additional Questions. Circ Res. 2020;127(11):1419–21.
Lippi G, Sanchis-Gomar F, Favaloro EJ, Lavie CJ, Henry BM. Coronavirus Disease 2019-Associated Coagulopathy. Mayo Clin Proc. 2021;96(1):203–17.
Barrett TJ, Lee AH, Xia Y, Lin LH, Black M, Cotzia P, et al. Platelet and vascular biomarkers associate with thrombosis and death in coronavirus disease. Circ Res. 2020;127(7):945–7.
Hottz ED A-QI, Palhinha L, Teixeira L, Barreto EA, Pão CRR, Righy C, Franco S, Souza TML, Kurtz P, Bozza FA, Bozza PT. Platelet activation and platelet-monocyte aggregates formation trigger tissue factor expression in severe COVID-19 patients. Blood. 2020;136(11):1330–41.
Campbell RA, Schwertz H, Hottz ED, Rowley JW, Manne BK, Washington AV, et al. Human megakaryocytes possess intrinsic antiviral immunity through regulated induction of IFITM3. Blood. 2019;133(19):2013–26.
Youssefian T, Drouin A, Masse JM, Guichard J, Cramer EM. Host defense role of platelets: engulfment of HIV and Staphylococcus aureus occurs in a specific subcellular compartment and is enhanced by platelet activation. Blood. 2002;99(11):4021–9.
Palma G, Pasqua T, Silvestri G, Rocca C, Gualtieri P, Barbieri A, et al. PI3Kdelta Inhibition as a Potential Therapeutic Target in COVID-19. Front Immunol. 2020;11:2094.
Goshua G, Pine AB, Meizlish ML, Chang CH, Zhang H, Bahel P, et al. Endotheliopathy in COVID-19-associated coagulopathy: evidence from a single-centre, cross-sectional study. The Lancet Haematology. 2020;7(8):e575–82.
Wibowo A, Pranata R, Lim MA, Akbar MR, Martha JW. Endotheliopathy Marked by High von Willebrand Factor (vWF) Antigen in COVID-19 is Associated with Poor Outcome: A Systematic Review and Meta-analysis. Int J Infect Dis. 2021;9:8.
Islam A, Khan MA, Ahmed R, Hossain MS, Kabir SMT, Islam MS, et al. Transcriptome of nasopharyngeal samples from COVID-19 patients and a comparative analysis with other SARS-CoV-2 infection models reveal disparate host responses against SARS-CoV-2. J Transl Med. 2021;19(1):32.
Xu H, Zhong L, Deng J, Peng J, Dan H, Zeng X, et al. High expression of ACE2 receptor of 2019-nCoV on the epithelial cells of oral mucosa. Int J Oral Sci. 2020;12(1):8.
Middleton EA, Rowley JW, Campbell RA, Grissom CK, Brown SM, Beesley SJ, et al. Sepsis alters the transcriptional and translational landscape of human and murine platelets. Blood. 2019;134(12):911–23.
Je S. Biochemistry and genetics of von Willebrand factor. Annu Rev Biochem. 1998;67:395–424.
Chen JCD. Inflammation, von Willebrand factor, and ADAMTS13. Blood. 2018;132:141–7.
Shaw RJ, Bradbury C, Abrams ST, Wang G, Toh CH. COVID-19 and immunothrombosis: emerging understanding and clinical management. Br J Haematol. 2021;194(3):518–29.
Manne BK, Denorme F, Middleton EA, Portier I, Rowley JW, Stubben C, et al. Platelet gene expression and function in patients with COVID-19. Blood. 2020;136(11):1317–29.
Best MG, Veld S, Sol N, Wurdinger T. RNA sequencing and swarm intelligence-enhanced classification algorithm development for blood-based disease diagnostics using spliced blood platelet RNA. Nat Protoc. 2019;14(4):1206–34.
Best MG, Sol N, Veld S, Vancura A, Muller M, Niemeijer AN, et al. Swarm intelligence-enhanced detection of non-small-cell lung cancer using tumor-educated platelets. Cancer Cell. 2017;32(2):238–52.
Campbell RA, Franks Z, Bhatnagar A, Rowley JW, Manne BK, Supiano MA, et al. Granzyme A in human platelets regulates the synthesis of proinflammatory cytokines by monocytes in aging. J Immunol. 2018;200(1):295–304.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Anders S, Pyl PT, Huber W. HTSeq–a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10(1):1523.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.
Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, et al. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016;44(W1):W90–7.
Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003;4:2.
This study was supported by the Emergency scientific research project (ZY202001-1) of Wenzhou City, Natural Science Foundation of Zhejiang Province (LY18H160058). The funders had no roles in study design, data collection and analysis, publication decision, or manuscript preparation.
Ethics approval and consent to participate
The current study was approved by the Institutional Review Board (IRB) of The Second Affiliated Hospital of Wenzhou Medical University (IRB# LCKY2020-09 and LCKY2020-222), and written informed consent was obtained from each participant prior to the commencement of the study.
Consent for publication
The authors have declared no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The imaging of all patients.
Hierarchical clustering of transcriptome-wide RNA expression. Hierarchical clustering of samples according to global gene expression segregates SARS-CoV-2-infected patients (red numbers) from healthy donors (gray numbers). SARS-CoV-2, severe acute respiratory syndrome coronavirus 2.
Functional characterization of key regulators in circulating platelets during SARS-CoV-2 infection. (A) Tissue and cell-specific enrichment analysis. (B) Disease enrichment analysis. (C) Enrichment analysis of differentially expressed RNAs in viral perturbations datasets from the Gene Expression Omnibus database. SARS-CoV-2, severe acute respiratory syndrome coronavirus 2.
Boxplots showing expression levels of four key regulators using RT-qPCR.
Lists of differentially expressed RNAs between COVID-19 patients and healthy donors.
Primers sequences for four candidate mRNAs and reference mRNA.
About this article
Cite this article
Ji, W., Chen, L., Yang, W. et al. Transcriptional landscape of circulating platelets from patients with COVID-19 reveals key subnetworks and regulators underlying SARS-CoV-2 infection: implications for immunothrombosis. Cell Biosci 12, 15 (2022). https://doi.org/10.1186/s13578-022-00750-5