Skip to main content


Front. Oncol., 26 April 2019
Sec. Cancer Epidemiology and Prevention
Volume 9 - 2019 |

Is Epstein-Barr Virus Infection Associated With Thyroid Tumorigenesis?—A Southern China Cohort Study

Shi-Tong Yu1 Jun-Na Ge1 Rui-Chen Li2 Zhi-Gang Wei1 Bai-Hui Sun1 Yu-Ming Jiang1 Jing-Yi Luo3 Hao Liu1 Shang-Tong Lei1*
  • 1Department of General Surgery, Nanfang Hospital, Southern Medical University, Guangzhou, China
  • 2Department of Radiation Oncology, Eye & ENT Hospital, Fudan University, Shanghai, China
  • 3State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou, China

Background: Epstein-Barr virus (EBV) is associated with many epithelial malignancies. A few reports on the association between EBV and thyroid tumorigenesis have been investigated. However, the conclusion is highly contradictory. We aimed to explore the role of EBV in thyroid nodule development and its clinical significance in a cohort from southern China.

Method: We conducted a retrospective data abstraction study of patients who underwent thyroidectomy between December 2017 and June 2018. We retrospectively analyzed the clinicopathological parameters and EBV infection status (serological antibodies and in situ hybridization).

Result: The cohort comprised 384 patients with newly diagnosed thyroid diseases, including 261 papillary thyroid carcinomas, 87 nodular goiters, 21 follicular adenomas, 12 follicular thyroid carcinomas, and 3 medullary thyroid carcinomas. Forty-two (10.9%) patients were identified as being serological antibody positive. However, there was no association between the clinicopathological parameters and serological antibody positivity. Additionally, none of the patients showed EBER expression in thyroid normal/cancer cell nuclei in in situ hybridization.

Conclusion: In this study, no correlation between EBV and thyroid diseases was found in a cohort from southern China.


Epstein-Barr virus (EBV) is a well-known human tumor virus with a very high prevalence in the population, especially for children and youth. EBV serum (IgG) is positive in an estimated 95% of the world's population (1). EBV infection is associated with epithelial and lymphoid malignancies, including nasopharyngeal carcinoma (NPC), gastric cancer, Hodgkin's lymphoma and Burkitt's lymphoma (2, 3). Although EBV has B-lymphocyte tropism, it can also infect T lymphocytes, myocytes, and epithelial cells in the oropharynx and stomach (4). Once EBV infects a host cell, it starts to induce a lytic or latent infection with diverse genes expressed. EBV nuclear antigens (EBNA 1, 2, 3A, 3B, 3C, and LP), the latent membrane proteins (LMP 1, 2A, and 2B) and two small noncoding RNAs (EBV-coded small RNA, EBER-1, and EBER-2) are expressed during the infection (3). These genes collaborate to induce tumorigenesis by causing systematic inflammation, suppressing the antitumoral immune system, and preventing anoikis.

Whether EBV infects the thyroid gland remains controversial. To date, only a handful of reports on the association between EBV and thyroid tumorigenesis have been investigated. Stamatiou et al. (5) summarized publications regarding the EBV detection rate in thyroid cancer specimens from 2001 to 2015. The conclusion was inconclusive because the results were highly contradictory, ranging from negative to 100% positive for EBV infection.

Therefore, in the current study, we explored the role of EBV in thyroid nodule development and its clinical significance in a cohort from southern China.


Patient Selection and Sample Collection

A total of 384 patients, including 261 with papillary thyroid carcinomas, 87 with nodular goiters, 21 with follicular adenomas, 12 with follicular thyroid carcinomas, and 3 with medullary thyroid carcinomas, who underwent thyroidectomy between December 2017 and June 2018 at the Department of General Surgery, Nanfang Hospital, Southern Medical University, were identified as being eligible for the study. All the data were extracted from the database of the Department of General Surgery, Nanfang Hospital, Southern Medical University. The patients included in the study met the following criteria: (1) primary thyroid neoplasms (including thyroid cancer, nodular goiter and follicular adenoma) confirmed by post-operative pathological results, (2) no history of thyroid/neck surgery, (3) no previous diagnosis of nasopharyngeal carcinoma (NPC) or other EBV-related disease, (4) no history of neck radiotherapy, (5) exclusion of cervical metastatic cancer, parathyroid neoplasms and Graves' disease, and (6) sufficient medical history. All patients were invited to donate a 5-mL blood sample for storage when they received preoperative examinations in our department. Blood samples were allowed a maximum of 6 h at room temperature before serologic analysis. The serum samples were then separated and divided into 4 tubes. All the serum samples were stored at −80°C at the Department of General Surgery, Nanfang Hospital, Southern Medical University. All the tissues samples were formalin fixed and paraffin embedded (FFPE) and then were cut into 4-μm-thick sections. All the clinicopathological parameters were recorded and evaluated according to the criteria of the American Joint Committee on Cancer (AJCC, 8th edition). This study was approved by the ethical committee of Nanfang Hospital, Southern Medical University. Informed consents were obtained from all involved patients when they were admitted to the hospital to provide their test results/specimens for future medical research.

In situ Hybridization

The EBER-ISH (in situ hybridization) test is the most commonly employed method and is regarded as the gold standard to diagnose EBV-infected diseases. We evaluated the presence of EBV in thyroid cells by in situ hybridization (ISH) analysis using EBV-encoded small RNA (EBER) probes, PNA probe/FITC (code Y5200), and the PNA ISH detection kit (code K5201) (Dako, Denmark) on the FFPE samples. The protocol was performed according to the manufacturer instructions. Briefly, the slides were baked at 60°C for 2 h and then were deparaffinized using a standard protocol. Deparaffinized slides were boiled for 20 min in pH 6.0 citric buffer for antigen retrieval. The slides were further permeabilized by Protease III treatment at room temperature for up to 10 min. The slides were then hybridized at 40°C for 2 h, followed by amplification and detection by adding Amp 1–4. The primary antibody incubation time at 4°C varied from 12 to 16 h. The substrate was incubated for 60 min, followed by counterstaining with eosin and mounting using Aquamount (Dako). EBER expression was graded from negative (–) to slight (+), moderate or intensive (+++), where most cells express EBER-RNA. The procedure was conducted by an experienced pathologist. Additionally, ISH was performed individually using 3 slides in each patient. Previously known cases of EBV-positive Hodgkin's lymphoma and NPC were used as positive controls. Because there is no established cutoff for EBER in solid malignancies other than NPC or gastric cancer (6), we referred this cutoff (>5%) of dark-brown staining of the tumor nucleus as positive for EBER transcript expression.

Serologic Antibody Analysis

Serologic antibody analysis was conducted in all patients using EBV-specific capsid antigen (VCA/IgA) antibodies and EBV-specific early antigen (EA/IgA) antibodies. Serological analysis was performed on July 2018 at the Department of Laboratory Medicine, Nanfang hospital, Southern Medical University because VCA-IgA and EA-IgA were routinely tested for screening NPC in Southern China institutions (79). Validation of the serological data was existing results from patients who underwent these two tests at our hospital between December 2017 and June 2018. Positive groups were patients with NPC, and healthy controls were the normal population from the medical examination center. Their ages and sex were matched with those in the thyroid disease cohort (1:1). EBV-specific VCA/IgA antibodies and EA/IgA antibodies were measured by ELISA (Euroimmun, Lubeck, Germany). The levels of these seromarkers were measured photometrically, according to the manufacturers' instructions. Additionally, these results were standardized by calculating the ratio of the optical density (OD) of the sample over that of the reference control (rOD). If the rOD value was >1, the sample was regarded as positive (10, 11). Either VCA/IgA or EA/IgA was positive, and the patients were regarded as serological antibody positive.

Statistical Analysis

All continuous variables were expressed as medians (Percentile 25, Percentile 75). Statistical analysis was conducted using SPSS 22.0 (SPSS, Inc, Chicago, IL, USA). To explore the relationship between EBV and the clinical pathological features in PTC, such as gender, age, and tumor size, chi-squared test and Fisher's exact test were used as appropriate; p < 0.05 considered as statistically significant.


The clinicopathological features are presented in Table 1. The median age of the included patients was 45 years (40, 61). The median number of resected lymph nodes in Level VI (for those with central neck dissection) and Level II-Vb (for those with lateral neck dissection) were 11.5 (7, 15) and 40 (25.75, 58.75), respectively. Two hundred eighty eight patients came from Guangdong Province (Southern China), and 96 patients came from other provinces of China (Figure 1).


Table 1. Baseline characteristics of all thyroid neoplams patients.


Figure 1. Overall distribution of the study population. In total, 384 patients with thyroid disease who underwent thyroidectomy were recruited. A total of 187 patients were from the Pearl River Delta (Green Area), 24 (12.8%) patients were positive for serological antibodies test, 42 patients were from the Eastern Guangdong (Blue Area), 5 (11.9%) patients were serological positive, 36 patients and 23 patients were from the Western (Red Area) and the Northern (Yellow Area) Guangdong, with a positive rate of 11.1% (4 patients) and 8.7% (2 patients), respectively. Ninety-six patients were from other provinces of China, and 7 patients (7.3%) were serologically positive.

Serological Antibodies

The VCA/IgA and EA/IgA antibodies were tested in all 384 patients with thyroid nodules as well as the population of positive/normal controls. There were 42 (10.9%) patients who were either VCA/IgA antibody or EA/IgA antibody positive or positive for both antibodies; therefore, these patients were designated as the serological antibody positive group. Comparing all clinicopathological features between the serological antibody groups, none of these parameters with statistical significance were identified (Supplementary Tables 1, 2). The number of patients with available serological results was 857 and 6,923 for NPC (positive control) and normal population (healthy control). After matched ages and sex, 384 NPC patients (positive control) and 384 healthy controls (normal population) were compared with those of the thyroid cohort (Supplementary Table 3). There were 311 (81.0%) patients who were serological positive in positive controls, and 40 (10.4%) patients with positive serological tests in the normal population.

In situ hybridization

None of the samples analyzed by ISH showed EBER expression in thyroid normal/cancer cell nuclei, even those with positive results in serological antibody analysis (Figure 2).


Figure 2. Results of hematoxylin-eosin (HE) staining of Hodgkin's lymphoma: (A) 10×; (B) 40×; Results of EBER-in situ hybridization (ISH) of Hodgkin's lymphoma: (C) 10×; (D) 40×; Results of HE staining of nasopharyngeal carcinoma (NPC): (E) 10×; (F): 40×; Results of EBER-ISH of NPC: (G) 10×; (H) 40×; Results of EBER-ISH of papillary thyroid carcinoma: (I) 10×; (J) 40×.


EBV has been revealed to be associated with the development of many cancers, such as gastric cancer, NPC, and Hodgkin's lymphoma (3). Chronic inflammation induced by EBV infection may play a significant role in the progression of cancer (4). However, the relationship between thyroid tumorigeneses and EBV has not been fully elucidated with conflicting results. The preliminary investigation of EBV in thyroid lymphoma was inspired by EBV persistently infecting B lymphocytes, contributing to lymphoma formation. In 2003, Shimakage et al. first reported EBV infection in other types of thyroid malignancies with a Japanese cohort (12). Shimakage et al. (12) explored the potential involvement of EBV expression in the progression of thyroid cancer by examining different thyroid neoplasm specimens, ranging from PTC to anaplastic thyroid carcinoma (ATC). The specimens were subjected to PCR, RT-PCR, ISH and indirect immunofluorescence staining. The results showed that mRNA and protein were positive for all carcinoma specimens and their expression was prominent in ATCs. For benign nodules, they showed no signal or very few signals during ISH. A similar result was conducted by Moghoofei et al. (13). The authors determined EBV infection by qRT-PCR and revealed that the EBV detection rate in PTC was similar to that of the healthy control. However, EBV positivity was associated with the tumor stage. Additionally, based on the PCR result, several inflammatory factors, such as Survivin, CD44, NF-kappaB, and IL-6, were higher expressed in the EBV-positive groups, and the mRNA expression of EBER1 and EBER2 was higher in thyroid tumor group. The reports from Almeida et al. (14) and Homayouni et al. (15) conducted similar results based on PCR and the ISH test in Brazilian and Iranian populations, respectively.

However, negative results for the association between thyroid tumors and EBV infection have been obtained. Despite a few infiltrating lymphoid cells in 1 (2.2%) of the specimens, none of the 45 PTCs were positive for the EBER-ISH test in a Japanese cohort (16). Bychkov et al. (17) reported that 1 of 20 thyroid cancer tissues contained single EBER-positive inflammatory cells. Cancer cells and normal thyroid tissues were consistently negative for ISH. Additionally, Tsai et al. (18) reported a negative association between benign tumors and EBV infection using ISH or PCR or Southern hybridization in a Taiwanese population.

To the best of our knowledge, this is the first study regarding EBV and thyroid neoplasms based on serological antibodies and ISH analyses in a cohort from the southern part of China. Similar to other southeast Asia neighbors, EBV is highly prevalent in southern China. Additionally, NPC, which has a closer relationship with EBV infection, is more endemic than in any part of the world, especially in Guangdong Province and Hong Kong (19). VCA/IgA and EA/IgA antibodies could reflect the status of recent viral infection and, therefore, are widely used biomarkers for screening NPC in the southern China population (8). In the current study, the positivity rate of serological antibody analysis was 10.9%, which is similar to that in previous national population-based studies conducted in the 1970s (79). However, there is no meaningful results between the clinicopathological parameters and VCA/IgA or EA/IgA antibodies. Furthermore, we failed to detect EBER signals in thyroid cancer cells or normal thyroid cells based on the ISH test. Interestingly, two PTC patients who were serological test positive with concurrent NPC were enrolled in this cohort. These two patients had pathologically confirmed NPC by preoperative fiber-laryngoscopy accidentally. They all received radiotherapy for NPC after thyroidectomy and were disease-free at the last follow up. However, we failed to detect the EBER signal in thyroid cancer cells, although the NPC specimens showed nasopharyngeal carcinoma cells that were positive for ISH. Our data indicated that the EBV detection rate in thyroid tumor samples from the southern China population is extremely low, consistent with previous EBER-based reports but differed from PCR-based tests. Several reasons may explain our findings. First, there is difference in populations and geography among different studies. Second, diversity in EBV detection techniques and their interpretation may cause conflicting results. Bychokov A believed that not only viral load assays (PCR or qPCR) but histochemical assays (ISH, immunofluorescence and immunohistochemistry) should be implied to ensure the precise tissue detection of EBV (17). A false-positive result may be received because any tissue containing B lymphocytes may have traces of EBV DNA (20). Additionally, there is no standard for the interpretation of EBV positivity for the ISH EBER test in thyroid tumors. Previous studies with positive results, which set a cutoff of few (<5%) EBER-positive cancer cells and even a single EBER-positive cancer cell are questionable because, in hematological malignancy, <1% of the EBER-positive background lymphocytes are usually regarded as negative during evaluation (20).

Indeed, in the current study, there are many limitations. First, the retrospective study nature may cause inevitable bias. Second, serological antibody analysis only involved VCA/IgA and EA/IgA antibodies and may be insufficient for EBV viral load evaluation; more comprehensive assays should be included in further studies. Third, only the EBER-based ISH test was included in the current study. Other viral detection methodologies, such as the LMP-1 immunohistochemical test, could be implied to better evaluate EBV infection in a future study.


We found no correlation between EBV and thyroid diseases in our study based on the current evidence. Future studies are warranted to reveal the significance of EBV in thyroid tumor developments.

Ethics Statement

This study was approved by the ethics committee of Nanfang Hospital, Southern Medical University. Informed consents were obtained from all enrolled patients.

Author Contributions

S-TL and HL designed the study. S-TY, J-NG, and R-CL conducted the analysis. S-TY and R-CL drafted the manuscript. Z-GW, B-HS, Y-MJ, and J-YL participated in result interpretation. All authors have read the manuscript and approved for publication.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


S-TY received a grant from President Foundation of Nanfang Hospital, Southern Medical University (No. 2018C024). The funders played no role in study design, collection, analysis, interpretation of data, writing of the report, or in the decision to submit the paper of publication. They accept no responsibility for the contents. We thank Professor Shuang Wang (Department of Pathology, Nanfang Hospital, Southern Medical University) and Dr. Yan Li (Department of Pathology, Sun Yat-sen University Cancer Center) for pathological guidance. We also thank Dr. Huasong Zhang (Department of Otolaryngology, Sun Yat-sen Memorial Hospital, Sun Yat-sen University) for plotting the figures in this paper.

Supplementary Material

The Supplementary Material for this article can be found online at:


1. Morrison BJ, Labo N, Miley WJ, Whitby D. Serodiagnosis for tumor viruses. Semin Oncol. (2015) 42:191–206. doi: 10.1053/j.seminoncol.2014.12.024

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Tsai MH, Raykova A, Klinke O, Bernhardt K, Gartner K, Leung CS, et al. Spontaneous lytic replication and epitheliotropism define an Epstein-Barr virus strain found in carcinomas. Cell Rep. (2013) 5:458–70. doi: 10.1016/j.celrep.2013.09.012

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Tsao SW, Tsang CM, To KF, Lo KW. The role of Epstein-Barr virus in epithelial malignancies. J Pathol. (2015) 235:323–33. doi: 10.1002/path.4448

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Thorley-Lawson DA, Hawkins JB, Tracy SI, Shapiro M. The pathogenesis of Epstein-Barr virus persistent infection. Curr Opin Virol. (2013) 3:227–32. doi: 10.1016/j.coviro.2013.04.005

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Stamatiou DP, Derdas SP, Zoras OL, Spandidos DA. Herpes and polyoma family viruses in thyroid cancer. Oncol Lett. (2016) 11:1635–44. doi: 10.3892/ol.2016.4144

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Delecluse HJ, Feederle R, O'Sullivan B, Taniere P. Epstein Barr virus-associated tumours: an update for the attention of the working pathologist. J Clin Pathol. (2007) 60:1358–64. doi: 10.1136/jcp.2006.044586

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Yi Z, Yuxi L, Chunren L, Sanwen C, Jihneng W, Jisong Z, et al. Application of an immunoenzymatic method and an immunoautoradiographic method for a mass survey of nasopharyngeal carcinoma. Intervirology. (1980) 13:162–8. doi: 10.1159/000149121

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Cao SM, Liu Z, Jia WH, Huang QH, Liu Q, Guo X, et al. Fluctuations of Epstein-Barr virus serological antibodies and risk for nasopharyngeal carcinoma: a prospective screening study with a 20-year follow-up. PLoS ONE. (2011) 6:e19100. doi: 10.1371/journal.pone.0019100

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Liu Y, Huang Q, Liu W, Liu Q, Jia W, Chang E, et al. Establishment of VCA and EBNA1 IgA-based combination by enzyme-linked immunosorbent assay as preferred screening method for nasopharyngeal carcinoma: a two-stage design with a preliminary performance study and a mass screening in southern China. Int J Cancer. (2012) 131:406–16. doi: 10.1002/ijc.26380

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Li RC, Du Y, Zeng QY, Tang LQ, Zhang H, Li Y, et al. Epstein-Barr virus glycoprotein gH/gL antibodies complement IgA-viral capsid antigen for diagnosis of nasopharyngeal carcinoma. Oncotarget. (2016) 7:16372–83. doi: 10.18632/oncotarget.7688

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Ji MF, Huang QH, Yu X, Liu Z, Li X, Zhang LF, et al. Evaluation of plasma Epstein-Barr virus DNA load to distinguish nasopharyngeal carcinoma patients from healthy high-risk populations in Southern China. Cancer. (2014) 120:1353–60. doi: 10.1002/cncr.28564

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Shimakage M, Kawahara K, Sasagawa T, Inoue H, Yutsudo M, Yoshida A, et al. Expression of Epstein-Barr virus in thyroid carcinoma correlates with tumor progression. Hum Pathol. (2003) 34:1170–7. doi: 10.1053/j.humpath.2003.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Moghoofei M, Mostafaei S, Nesaei A, Etemadi A, Sadri Nahand J, Mirzaei H, et al. Epstein-Barr virus and thyroid cancer: the role of viral expressed proteins. J Cell Physiol. (2018) 234:3790–9. doi: 10.1002/jcp.27144

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Almeida JFM, Campos AH, Marcello MA, Bufalo NE, Rossi CL, Amaral LHP, et al. Investigation on the association between thyroid tumorigeneses and herpesviruses. J Endocrinol Invest. (2017) 40:823–9. doi: 10.1007/s40618-017-0609-y

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Homayouni M, Mohammad Arabzadeh SA, Nili F, Razi F, Amoli MM. Evaluation of the presence of Epstein-Barr virus (EBV) in Iranian patients with thyroid papillary carcinoma. Pathol Res Pract. (2017) 213:854–6. doi: 10.1016/j.prp.2017.01.020

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Kijima Y, Hokita S, Takao S, Baba M, Natsugoe S, Yoshinaka H, et al. Epstein-Barr virus involvement is mainly restricted to lymphoepithelial type of gastric carcinoma among various epithelial neoplasms. J Med Virol. (2001) 64:513–8. doi: 10.1002/jmv.1079

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Bychkov A, Keelawat S. Epstein-Barr virus and thyroid cancer: the controversy remains. J Endocrinol Invest. (2017) 40:891–2. doi: 10.1007/s40618-017-0703-1

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Tsai JH, Tsai CH, Cheng MH, Lin SJ, Xu FL, Yang CC. Association of viral factors with non-familial breast cancer in Taiwan by comparison with non-cancerous, fibroadenoma, and thyroid tumor tissues. J Med Virol. (2005) 75:276–81. doi: 10.1002/jmv.20267

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Chang ET, Adami HO. The enigmatic epidemiology of nasopharyngeal carcinoma. Cancer Epidemiol Biomarkers Prev. (2006) 15:1765–77. doi: 10.1158/1055-9965.EPI-06-0353

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Gulley ML, Tang W. Laboratory assays for Epstein-Barr virus-related disease. J Mol Diagn. (2008) 10:279–92. doi: 10.2353/jmoldx.2008.080023

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: thyroid cancer (TC), Epstein-Barr virus (EBV) infection, EBER = EBV-encoded small RNA, EBV-specific capsid antigen (VCA/IgA), EBV-specific early antigen (EA/IgA), Southern China

Citation: Yu S-T, Ge J-N, Li R-C, Wei Z-G, Sun B-H, Jiang Y-M, Luo J-Y, Liu H and Lei S-T (2019) Is Epstein-Barr Virus Infection Associated With Thyroid Tumorigenesis?—A Southern China Cohort Study. Front. Oncol. 9:312. doi: 10.3389/fonc.2019.00312

Received: 22 January 2019; Accepted: 05 April 2019;
Published: 26 April 2019.

Edited by:

Hajo Zeeb, Leibniz Institute for Prevention Research and Epidemiology (LG), Germany

Reviewed by:

Sylvia L. Asa, University of Toronto, Canada
Liting Liu, Sun Yat-Sen University Cancer Center (SYSUCC), China

Copyright © 2019 Yu, Ge, Li, Wei, Sun, Jiang, Luo, Liu and Lei. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Shang-Tong Lei,

These authors have contributed equally to this work