APP下载

Blockage of ETS homologous factor inhibits the proliferation and invasion of gastric cancer cells through the c-Met pathway

2021-01-15MengLiGuXinXinZhouMengTingRenKeDaShiMoSangYuWenRuiJiaoYaMeiWangWeiXiangZhongFengJi

World Journal of Gastroenterology 2020年47期

Meng-Li Gu, Xin-Xin Zhou, Meng-Ting Ren, Ke-Da Shi, Mo-Sang Yu, Wen-Rui Jiao, Ya-Mei Wang, Wei-Xiang Zhong, Feng Ji

Abstract

Key Words: Gastric cancer; ETS homologous factor; c-Met; Antineoplastic target; Signaling pathway

INTRODUCTION

Gastric cancer (GC) is a common malignant tumor in the gastrointestinal tract and the third leading cause of cancer-related mortality worldwide[1]. Risk factors such as Helicobacter pylori infection, atrophic gastritis, a high-salt diet and heavy alcohol drinking contribute to the tumorigenesis and progression of GC, but its exact pathogenesis has not been fully elucidated[2,3]. Currently, surgical resection represents the best curative option for GC. On the other hand, advanced novel targets for GC diagnosis and treatment have resulted in decreased incidence and mortality rates in recent years[4]. However, the 5-year overall survival rate of patients with advanced GC is less than 30% due to delayed diagnosis and recurrence[5]. Therefore, revealing the oncogenic mechanisms of GC and further exploring the strategies for GC diagnosis and treatment are urgent.

The ETS family consists of various transcriptional factors that contain a conservative ETS structure domain, which specifically identifies and binds to the GGAA/T sequence within the promoter or enhancer region of target genes to regulate their transcription[6]. ETS homologous factor (EHF) is a member of the ETS family of transcription factors expressed in multiple epithelial cell types, which participates in the oncogenesis of various malignant tumors via regulating cell proliferation and differentiation[7,8]. The aberrant expression of EHF is involved in the epithelial-tomesenchymal transition (EMT) and is a risk factor for the increased recurrence and reduced overall survival of prostate cancer[9]. A recent study demonstrated that upregulated EHF expression is correlated with the poor prognosis of ovarian cancer patients and that the inhibition of EHF exerts an antineoplastic effect on ovarian cancer cells[10]. Additionally, EHF is highly expressed in colon cancers, and the downregulation of EHF leads to apoptosis of colon cancer cells[11,12]. However, the roles of EHF in the tumorigenesis and progression of GC need further investigation.

A previous study indicated that EHF regulates the expression of c-Met by directly binding to the promoter region of the c-Met coding sequence[7]. c-Met is a receptor tyrosine kinase (RTK), and its dysregulation is closely related to the pathogenesis and development of GC[13]. In advanced GC patients, c-Met-positive status is an independent prognostic factor for poor outcomes[14]. To date, whether EHF participates in the development of GC through the c-Met pathway remains unclear. Herein, we investigated EHF expression in GC and elucidated its potential function and regulatory mechanisms involving the oncogenesis of GC, providing important evidence for the exploration of novel therapeutic targets for GC.

MATERIALS AND METHODS

Human samples

A total of 10 paired GC tissues and adjacent normal tissues were obtained from The First Affiliated Hospital of Zhejiang University (Zhejiang, China). The patients did not receive radiotherapy or chemotherapy prior to collection. The tissues were immediately stored in liquid nitrogen postoperatively. This study was performed according to the Declaration of Helsinki and approved by The Research Ethics Committee of The First Affiliated Hospital, College of Medicine, Zhejiang University, No. 20191218. Informed consent was obtained from all subjects involved in the study.

Cell culture

The human normal gastric epithelial cell line (GES-1) and human GC cell lines (BGC-823, KATO III and SGC-7901) were purchased from the cell bank of Chinese Academy of Sciences (Shanghai, China). The cell lines were cultured in Roswell Park Memorial Institute (RPMI) 1640 medium (Gibco, Carlsbad, CA, United States) supplemented with 10% fetal bovine serum (FBS; Gibco) at 37 °C in a humidified atmosphere of 5% CO2and 95% air.

Transfection of small interfering RNAs

Small interfering RNA (siRNA) transfections were performed using Lipofectamine™ 3000 reagent (Invitrogen, San Diego, CA, United States) in antibiotic-free RPMI 1640 according to the manufacturer’s instructions. The GC cell lines were transfected with 10 nmol/L siRNA. The control siRNA and siRNA targeting EHF were purchased from Invitrogen.

Quantitative PCR

The total RNA of each cell line was extracted using RNAiso Plus reagent (Takara, Dalian, China) following the manufacturer's protocol. The quantitative polymerase chain reaction (qPCR) assay was performed as previously described[15]. The primers were obtained from Sangon Biotech (Shanghai, China). Primer sequences were listed as follows: EHF, forward 5'-GAAGAACAACAGCAGCATGACC-3' and reverse 5'-TCAGGTTTCGGTGTATGAGTTG-3'; GAPDH, forward 5'-AGAAGGCT GGGGCTCATTTG-3' and reverse 5'-AGGGGCCATCCACAGTCTT-3'. The relative gene expression levels of the tested genes were calculated using the 2-ΔΔCtmethod[16].

Cell proliferation

Cell proliferation was determined using the Cell Counting Kit-8 assay according to our previous study[15].

Colony formation assay

After transfection for 48 h, BGC-823, KATO III, and SGC-7901 cells were seeded into 6-well plates at a density of 2 × 103cells per well and incubated for 10 d at 37 °C. Then, the culture medium was removed, and cells were washed with phosphate-buffered saline (referred to as PBS). The colonies were fixed in 4% paraformaldehyde and stained with 0.1% crystal violet. The colony numbers were ultimately counted.

Cell cycle assay and cell apoptosis assay

Cells were inoculated into 6-well plates at a density of 4 × 105cells per well and incubated for 24 h. Then the cells were transfected with negative control and EHFtargeted siRNA for 72 h. Finally, cell cycle assays and cell apoptosis assays were performed as previously described[15].

Transwell assay

Following the transfection of EHF-targeted siRNA for 24 h, 2.5 × 104cells were plated into 24-well chambers (Corning, NY, United States) with medium containing 1% FBS. Medium containing 10% FBS was added into the lower wells. For the invasion assays, Matrigel was added to each transwell chamber and incubated for 5 h before the cells were seeded. After 48 h of incubation, the cells remaining in the upper chamber were carefully removed with cotton swabs. The membrane-penetrating cells were stained with crystal violet for 15 min. The cells that migrated and invaded were counted from at least five randomly selected fields.

Immunohistochemical analysis

The expression of EHF was evaluated by immunohistochemical staining. Paraffin blocks that contained formalin-fixed GC and adjacent normal tissues were sectioned at 4 μm. The sections were deparaffinized in 100% xylene and re-hydrated at graded ethanol concentrations. Endogenous peroxidase activity was blocked using fresh 3% hydrogen peroxide for 10 min at room temperature. Antigen retrieval was done by boiling the slides in sodium citrate buffer (pH 6.0). After blocking in 3% bovine serum albumin, sections were incubated with anti-EHF antibody (1:100 dilution in PBS; Abcam, Cambridge, MA, United States) overnight at 4 °C. Subsequently, the sections were incubated for 30 min at room temperature with anti-rabbit horseradish peroxidase-conjugated secondary antibody. The immunostaining was performed using diaminobenzidine solution. The sections were then counterstained with hematoxylin. dehydrated and mounted. Blinded evaluations of EHF immunostaining were carried out independently by two experienced pathologists.

Western blot analysis

According to the manufacturer's instructions, the total protein was extracted from cells by using cell lysis buffer (Cell Signaling, Beverly, MA, United States). Western blotting analysis was performed as previously described[15]. Rabbit antibodies corresponding to EHF were purchased from Abcam. Antibodies against c-Met, Ras, total c-Raf, phospho-c-Raf (p-c-Raf), total Erk1/2, p-Erk1/2, phosphatase and tensin homolog (PTEN), total glycogen synthase kinase 3 beta (GSK-3β), phospho-GSK-3β, Bcl-2-associated X protein (Bax), survivin, total signal transducer and activator of transcription 3 (STAT3), p-STAT3, Twist, Snail, E-cadherin, and GAPDH were purchased from Cell Signaling Technology, Inc.

Statistical analyses

Statistical analyses were performed using SPSS version 20 (IBM Corp., Armonk, NY, United States). Results are presented as the mean ± standard deviation. The statistical differences between two groups was assessed by the Student’s t-test. One-way analysis of variance followed by Tukey's multiple comparison test was used to determine the statistical difference among multiple groups. The significance level was set at P < 0.05. All experiments were performed in triplicate.

RESULTS

EHF is overexpressed in GC specimens and cell lines in which c-Met expression is also upregulated

The expression of EHF mRNA in paired adjacent normal tissues and GC tissues were first measured via qPCR and Western blotting. The mRNA expression level of EHF was significantly upregulated in GC tissues compared with adjacent control tissues (Figure 1A). The protein expression of EHF was increased in GC tissues overexpressing c-Met (Figure 1B). As shown in Figure 1C, increased expression of EHF was also detected in GC tissues by immunohistochemistry. Second, the gene and protein expression of EHF in human GES-1 and GC cell lines (BGC-823, KATO III and SGC-7901) was detected. The expression of EHF mRNA was increased in these GC cell lines compared with GES-1 (Figure 1D). In GC cell lines with upregulated EHF, the protein expression of c-Met was also increased in a similar pattern (Figure 1E).

Figure 1 ETS homologous factor and c-Met expression in gastric cancer tissues and cell lines. A: Quantitative PCR analysis of mRNAs revealed the upregulation of ETS homologous factor (EHF) mRNA in gastric cancer (GC) tissues compared with their matched adjacent tissues; B: Western blotting results showed that the expression of EHF and c-Met were both upregulated in GC tissues; C: Immunohistochemical staining of EHF indicated that increased expression of EHF was higher in GC tissues than that in matched adjacent tissues (400 ×); D: The mRNA expression level of EHF was enhanced in GC cell lines compared with that in normal gastric epithelial cells; E: EHF protein expression was increased in GC cell lines in which c-Met was also upregulated. aP < 0.05 vs matched adjacent tissues; eP < 0.001 vs normal gastric epithelial cells. T: GC tissues; N: Matched adjacent tissues.

Silencing EHF by siRNA leads to c-Met downregulation

To clarify the potential function of EHF in the tumorigenesis of GC, EHF silencing was performed by transfecting EHF-targeted siRNA into human GC cell lines (BGC-823, KATO III and SGC-7901). As shown in Figure 2, the mRNA and protein expression levels of EHF significantly decreased in the GC cell lines compared with the negative control group after transfection for 72 h (P < 0.05). In addition, the inhibition of EHF also led to the downregulation of c-Met in these GC cell lines.

EHF silencing suppresses GC cell viability and colony formation

The aforementioned results demonstrated that EHF expression was higher in GC cell lines and GC tissues. The potential functional roles of EHF in GC cells were further explored. Knockdown of EHF significantly reduced the cell viability of GC cell lines after transfection for 5 d (Figure 3A). In addition, EHF silencing also significantly inhibited the colony formation of GC cell lines compared with the negative control group (Figure 3B and C). Thus, EHF plays a crucial role in the regulation of proliferation in GC.

Depletion of EHF promotes apoptosis and induces cell cycle arrest

To investigate whether EHF can affect apoptotic cell death in GC, Annexin VFITC/propidium iodide (PI) staining and flow cytometry were performed. As shown in Figure 4A and B, EHF downregulation induced significant apoptosis in the GC cell lines (P < 0.05), suggesting that EHF exerts an antiapoptotic effect in these GC cell lines. To further assess the cell cycle status following EHF silencing, PI staining, and flow cytometry were performed. As shown in Figure 4C and D, the inhibition of EHF caused an increase in GC cells in the G1 phase, but a significant decrease in the S phase (P < 0.05). These data indicate that EHF participates in regulating cell cycle progression in GC cell lines by promoting the transition from the G1 phase to the S phase.

Downregulation of EHF inhibits the migration and invasion of GC cells

Transwell assays were used to detect cell migration after transfection for 24 h. As shown in Figure 5A and B, the migration of GC cell lines significantly decreased following EHF silencing (P < 0.01). The effect of EHF on GC cell invasion ability was further investigated. As shown in Figure 5C and D, the number of cells that invaded through the filters significantly decreased in the EHF silencing group compared with the negative control group in these GC cell lines (P < 0.01). These results suggested that EHF promotes the migration and invasion of BGC-823, KATO III and SGC-7901 cells, contributing to the high metastatic potential of GC.

EHF may regulate GC cell malignant behaviors through the c-Met pathway

Previous findings have demonstrated that c-Met has a pivotal role in the development and progression of GC[17]. Aberrantly activated c-Met triggers tumor growth, survival, metastasis, and EMT via mitogen-activated protein kinases (MAPK), PTEN, and STAT3 signaling pathways[18-20]. Consequently, the expression of c-Met protein and its downstream molecules in the GC cell lines following EHF downregulation were detected. As shown in Figure 6A, the results suggested that EHF silencing led to the considerable inhibition of c-Met and Ras expression in BGC-823, KATO III and SGC-7901 cells. Moreover, the inactivation of c-Raf and Erk1/2 were both observed after EHF downregulation. In addition to the MAPK pathway, expression of the tumor suppressor molecule PTEN in these GC cells was significantly upregulated (Figure 6B). The decreased phosphorylation of GSK-3β, the downstream target of PTEN, was observed in BGC-823 and KATO III cells. Moreover, the expression of antiapoptotic survivin was remarkably inhibited in these GC cells, whereas the expression of proapoptotic Bax was significantly upregulated in BGC-823 and KATO III cells. Furthermore, the data indicated that STAT3, a signal transducer of the c-Met pathway, was inactivated in BGC-823, KATO III, and SGC-7901 cells following the downregulation of EHF (Figure 6C). The protein expression of the downstream transcriptional factors of STAT3 was also examined. Interestingly, the expression of Twist decreased in BGC-823 and SGC-7901 cells, while the inhibitory effect of Snail expression was considerable in KATO III and SGC-7901 cells. E-cadherin is a target of Twist and Snail, while it is also known as a biomarker for the EMT. Herein, we found that E-cadherin was significantly increased in the EHF silencing group, indicating an inhibitory effect on EMT in GC cells.

Figure 2 Transfection of ETS homologous factor-targeted small interfering RNA reduces c-Met expression in gastric cancer cell lines. A: Quantitative PCR was performed to evaluate the mRNA levels of ETS homologous factor (EHF) in gastric cancer (GC) cells following transfection for 72 h; B: Western blotting was used to detect the protein expression of EHF and c-Met after small interfering RNA transfection. bP < 0.01 vs negative control (NC) group; eP < 0.001 vs NC group.

DISCUSSION

The ETS family is a transcription factor family consisting of various members, which can be divided into 11 subfamilies according to their structural characteristics[21]. EHF belongs to the epithelium-specific ETS (ESE) subfamily, and consequently is also known as ESE-3. According to previous studies, EHF is mainly expressed in multiple types of epithelial cells; its dysregulation is often detected in multiple malignant tumors, and it is associated with poor clinical outcomes[9,22-25]. However, the biological roles and potential regulatory mechanisms of EHF in GC is not completely clear.

In this study, qPCR was performed to detect the expression of EHF mRNA in GC cell lines and tissues compared with that in the gastric mucosal epithelial cell line and the adjacent tissues. The data demonstrated that EHF is overexpressed in GC tissues and three GC cell lines (BGC-823, KATO III, and SGC-7901), which is consistent with the alterations in protein expression. As a potential downstream target of EHF, the expression of c-Met was also increased in GC tissues and cell lines with high EHF expression. c-Met is an RTK encoded by the proto-oncogene MET, and hepatocyte growth factor is its only known ligand, whose aberrant expression and functional changes are closely related to the occurrence and development of GC[26]. The expression of c-Met in GC tissues is positively associated with distant metastasis and poor patient prognosis[14]. Several in vitro experiments have indicated that activated c-Met triggers its downstream MAPK, Akt/mTOR and other signaling pathways, thereby inhibiting the apoptosis of tumor cells and promoting their proliferation and migration[27]. Additionally, the inhibition of c-Met leads to the apoptosis of GC cells and promotes their sensitivity to chemotherapeutic agents[28].

Figure 3 The effects of ETS homologous factor silencing on gastric cancer cell proliferation. A: Cell viability was evaluated by the Cell Counting Kit-8 after transfection for 1, 3, 5, and 7 d; B and C: Cell colony formation was measured to verify the changes in gastric cancer (GC) cell proliferation following ETS homologous factor downregulation. eP < 0.001 vs negative control (NC) group.

The results demonstrated that EHF silencing can significantly reduce cell proliferation, promote apoptosis and induce cell cycle arrest in GC cell lines (BGC-823, KATO III and SGC-7901). On the other hand, both the migration and invasion abilities of these GC cell lines were significantly inhibited following the interference of EHF. These data indicate that EHF positively regulates the malignant biological behaviors of GC. Moreover, knockdown of EHF by specific siRNA led to the downregulation of c-Met in both GC tissues and cell lines. Currently, whether EHF participates in the tumorigenesis and progression of GC via the c-Met pathway remains unknown. Previous investigations have shown that the inhibition of c-Met in GC results in the suppression of cell proliferation both in vitro and in vivo[29]. Selective c-Met inhibitors block cell DNA synthesis and induce cell cycle arrest at the G1-S phase transition in GC cells[30]. The inactivation of c-Met can also reduce migration and invasion in vitro and ultimately block the malignant progression of GC[31]. EHF directly regulates c-Met expression via binding to its promoter region[7]. A recent study showed that EHF contributes to the metastasis and chemoresistance of oral squamous cell carcinoma through the positive regulation of the c-Met pathway[32]. Thus, EHF may be involved in the regulation of the biological behaviors of GC via promoting c-Met expression.

Recently, it has been verified that the abnormally upregulated expression of at least one RTK (human epidermal growth factor receptor 2, epidermal growth factor receptor, c-Met or fibroblast growth factor receptor 2) was found in approximately two-thirds of GC patients, indicating that the detection of RTK expression levels or their gene mutations can provide important evidence for the diagnosis and treatment of GC[33]. However, the identification of one single type of RTK expression or mutation is inadequate to accurately predict chemotherapy sensitivity or overall prognosis[34]. In recent years, the detection of multiple RTKs has improved the efficiency of predicting the prognosis to some extent, whereas the accuracy of RTK pathway analysis is still limited by the variable mutation sites and complex mutation types of these coding genes. It is possible to determine the abnormal expression or activation of RTKs by examining the expression level of the upstream regulators that affect RTK transcription regardless of their mutation status. Consequently, the combined detection of EHF and c-Met may serve as potential markers for the diagnosis and prognosis of GC.

Figure 4 Downregulation of ETS homologous factor induces apoptosis and cell cycle arrest in gastric cancer cells. A and B: The analysis of apoptosis was performed by flow cytometry to detect the induction of apoptosis caused by ETS homologous factor (EHF) silencing following V-FITC/propidium iodide (PI) staining; C and D: The cell cycle arrest induced by EHF downregulation was assessed by PI staining and further analyzed by flow cytometry. aP < 0.05 vs negative control (NC) group; bP < 0.01 vs NC group; eP < 0.001 vs NC group. GC: Gastric cancer.

Activation of c-Met and its key downstream MAPK cascades are responsible for maintaining the abnormal biological function of tumor cells[35]. Erk1/2 is a subtype of MAPK that is closely associated with cell proliferation, differentiation and migration[36]. The current study demonstrated that the inhibition of EHF leads to the impaired biological function of GC cell lines via the blockage of Ras-Erk1/2 cascades.

PTEN is a tumor suppressor, and its expression level is negatively correlated with the tumor-node-metastasis stage of GC. In vivo experiments have confirmed that the absence of PTEN promotes the tumorigenesis of GC[37]. It has been reported that PTEN interacts with the c-Met signaling pathway[38], and a high expression level of c-Met is associated with the downregulation of PTEN[39]. The inhibition of PTEN expression leads to the increased proliferation activity of GC cells via GSK-3β inactivation[40]. The present study demonstrated that EHF may affect the apoptosis and cell cycle progression of GC via the c-Met/PTEN/GSK-3β cascade. Interestingly, the upregulation of PTEN was detected in these GC cell lines, whereas considerable inactivation of GSK-3β was not observed in SGC-7901 cells. Previous findings have suggested that survivin is a direct target of GSK-3β, which inhibits survivin expression, leading to apoptosis in GC cells[41]. Herein, the expression of survivin increased in all the GC cells, indicating that survivin may not be completely controlled by GSK-3β in SGC-7901 cells.

In this study, STAT3 was deactivated following EHF silencing. In addition, as a biomarker for EMT, E-cadherin expression was upregulated in the GC cell lines. The activation of STAT3 promotes the expression of Twist and Snail that reduces Ecadherin expression, resulting in the EMT of malignant tumor cells[42-44]. EMT causes significant epithelial phenotypic changes and is associated with the progression, metastasis, and drug resistance of GC[45]. A previous study indicated that abnormally activated STAT3 plays an important role in c-Met-mediated EMT in GC[46]. Accordingly, the results suggested that EHF may regulate EMT in GC cell lines via the c-Met/STAT3 pathway. However, the downregulation of Twist was not observed in KATO III cells, and Snail expression was not significantly affected in BGC-823 cells, indicating that STAT3 may affect the EMT of BGC-823 and KATO III cells through Twist and Snail, respectively.

This study had some limitations. Further investigations with large sample sizes are needed to clarify the relationship between EHF and c-Met expression in GC. In the present study, the potential role of EHF in GC tumorigenesis was demonstrated and the possible mechanisms underlying the antineoplastic effects induced by EHF downregulation were further sought. Nevertheless, the advantage of EHF as an antineoplastic target for GC diagnosis and treatment has not yet been fully elucidated and remains to be validated both in vitro and in vivo.

CONCLUSION

In summary, EHF may promote cell proliferation and cell cycle progression and inhibit the apoptosis of GC cells via regulating the c-Met pathway. Furthermore, EHF may contribute to the migration and invasion of GC cells by inducing EMT via the STAT3 pathway. Therefore, by targeting the c-Met pathway, EHF can be regarded as a promising candidate target for the development of a novel antineoplastic strategy for GC treatment.

Figure 5 ETS homologous factor regulates the cell migration and invasion abilities of gastric cancer cells. A and B: The transwell assay was performed to assess the cell migration ability of gastric cancer (GC) cells following ETS homologous factor (EHF) downregulation; C and D: Invasion assays were performed by adding Matrigel into the Transwell chamber to determine the effects of EHF silencing on the invasion ability of GC cells. eP < 0.001 vs negative control (NC) group.

Figure 6 ETS homologous factor promotes the malignant biological behaviors of gastric cancer cells through the c-Met pathway. Western blotting was conducted to investigate the effects of ETS homologous factor (EHF) silencing on the expression and activities of signaling molecules in the c-Met pathway. A: Alterations in the Ras- extracellular signal-related kinase 1/2 (Erk1/2) cascade after EHF inhibition; B: The molecular changes in phosphatase and tensin homolog (PTEN) and glycogen synthase kinase-3β (GSK3β) and their potential targets following EHF silencing; C: The altered activities of signal transducer and activator of transcription 3 (STAT3) and the changed expression of its downstream effectors after EHF downregulation. GC: Gastric cancer; NC: Negative control group.

ARTICLE HIGHLIGHTS

Research methods

The expression of EHF mRNA in GC tissues and cell lines was measured by qPCR. Western blotting was performed to determine the protein expression of EHF, c-Met and its downstream signal molecules. The EHF expression in GC tissues was further checked by immunohistochemical staining. To investigate the role of EHF in GC oncogenesis, small interfering RNAs (siRNAs) against EHF were transfected into GC cells. The cell proliferation of GC cells was determined by cell counting kit-8 and colony formation assays. Flow cytometry was performed following annexin V/PI and PI staining to determine apoptosis and the cell cycle, respectively. Cell migration and invasion were assessed by transwell assays.

Research results

The data indicated that EHF was upregulated in GC tissues and cell lines in which increased expression of c-Met was also observed. Silencing of EHF by siRNA reduced cell proliferation in GC cells. Inhibition of EHF induced significant apoptosis and cell cycle arrest in GC cells. Cell migration and invasion were significantly inhibited. EHF silencing led to the c-Met downregulation and further blocked Ras/c-Raf/Erk1/2. Additionally, PTEN was upregulated and GSK-3β was deactivated. Moreover, inactivation of STAT3 was detected following EHF inhibition, leading to the blockage of the epithelial-to-mesenchymal transition (EMT).

Research conclusions

These results suggested that EHF plays a key role in cell proliferation, invasion, apoptosis, the cell cycle and EMT via the c-Met pathway.

Research perspectives

EHF may serve as an antineoplastic target for the diagnosis and treatment of GC.

ACKNOWLEDGEMENTS

We are very grateful to The State Key Laboratory for Diagnosis and Treatment of Infectious Diseases of The First Affiliated Hospital of Zhejiang University for providing excellent technical assistance.