Li Wen1†
Liping Guo2†Wen Zhang3†Yajian Li1Weixing Jiang1Xuebing Di2JianHui Ma1Lin Feng2*Kaitai Zhang2*Jianzhong Shou1*- 1Department of Urinary Surgery, National Cancer Center, National Clinical Research Center for Cancer, Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
- 2State Key Laboratory of Molecular Oncology, Department of Etiology and Carcinogenesis, National Cancer Center, National Clinical Research Center for Cancer, Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
- 3Department of Immunology, National Cancer Center, National Clinical Research Center for Cancer, Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
Most renal cell carcinoma (RCC) patients die from metastasis or recurrence after the spread of cancer to another organ, but the mechanisms underlying the intravascular survival of circulating tumor cells (CTCs) have not been completely deciphered. Additionally, although elevated plasma C-reactive protein (CRP) levels and thrombocytosis are strongly correlated and both indicate a poor prognosis for RCC patients, the bridge connecting inflammation and coagulation remains poorly understood. To explore the complicated relationship among inflammation, the coagulation system and CTC survival, we obtained viable CTC counts and clinical information from 106 treatment-naïve patients. In addition, we performed RNA sequencing on peripheral blood leukocytes from 21 of these patients. Patients with elevated CRP and fibrinogen (FIB) levels had higher CTC counts than patients with normal levels of these indexes. Each pair of the three variables (CTC count, CRP level and FIB level) was positively correlated. According to transcriptomic analysis of blood leukocytes, the functions of the 257 genes identified as being positively correlated with the CTC count indicated neutrophil extracellular trap (NET) formation. Indeed, gene set enrichment analysis (GSEA) suggested that NET formation or increased levels of NET markers would promote CTC viability. Additionally, the calculated NET score was positively correlated with the plasma FIB concentration, and both of these values were increased in patients with elevated CRP levels. Moreover, immunofluorescence staining showed that NETs were entangled with viable renal cancer cells and that the NET frameworks were decorated with NET-derived tissue factor (TF). Finally, analysis of 533 RCC samples from The Cancer Genome Atlas (TCGA) indicated that the NET score and TF value are independent prognostic factors for RCC patients. Collectively, NETs formed by intravascular neutrophils further activate the coagulation system. Both the DNA scaffold sprouted and fibrin net triggered by NETs anchor and shield CTCs from attack. Thus, degrading this framework maybe could destroy the double shelter of CTCs, the pioneers of metastasis.
Introduction
Although the global incidence of renal cancers ranked 16th among 35 common cancers in 2018 (1), renal cancers carry a poor prognosis. The prognosis of renal cell carcinoma (RCC), the most common kidney tumor, is particularly poor; the 5-year survival rate decreases sharply to 8% in metastatic patients (2). In addition, the recurrence rate is approximately 40% in patients treated with localized disease (3). Therefore, most RCC patients die from the metastasis or recurrence of the cancer after it spreads to another organ. However, little is known about the mechanisms by which circulating tumor cells (CTCs) survive intravascularly, which, in turn, reduces the opportunities to inhibit metastasis.
After shedding from the primary tumor into the circulatory system, a single CTC encounters countless survival challenges, such as attacks by T cells and natural killer (NK) cells, shear forces and oxidative stress (4). However, a few surviving cells can be detected and overtly metastasize, and studies on the mechanisms underlying CTC survival have thus far focused on two mutually dependent aspects. The first aspect involves CTC self-changes, such as the antiphagocytic signal CD47 being overexpressed in circulating colorectal cancer cells (5). The other aspect is the self-arming and self-protection of CTCs with suppressive cellular components; for example, CTCs recruit platelets to aggregate on their surface (6), and neutrophils adhere to circulating liver tumor cells expressing Mac-1 (7).
In RCC patients, elevated plasma C-reactive protein (CRP) levels and thrombocytosis are strongly correlated with poor prognosis and recurrence (8), indicating that inflammation and hypercoagulability may be critical environmental components required for CTC survival in circulation. Regarding inflammation, certain subtypes of neutrophils have been increasingly acknowledged as important components of protumor-associated inflammatory cells (9, 10). Moreover, neutrophils not only mediate inflammation but also initiate the coagulation system by forming a specific structure—neutrophil extracellular traps (NETs) (11).
NETs were first reported as a distinct antimicrobial mechanism in 2004 (12). Specifically, neutrophils are activated with lipopolysaccharide (LPS), phorbol myristate acetate (PMA), interleukin-8 (IL-8), bacteria (13), or granulocyte colony stimulating factor (G-CSF) (14), and histone H3 in neutrophils is citrullinated by the protein arginine deiminase 4 (PAD4) enzyme, which is essential for NET formation. The neutrophil nucleus then undergoes chromatin decondensation and extrudes a meshwork of chromatin fibers. Furthermore, the extruded DNA scaffold is decorated with a variety of enzymes, such as neutrophil elastase (NE), myeloperoxidase (MPO), cathepsin G, matrix metalloproteinase 9 (MMP9) and tissue factor (TF) (13, 15, 16). NET-derived TF, a trigger factor of the extrinsic coagulation pathway, is the link between inflammation and coagulation. Although the initial observation of NET formation was in the context of an innate immune response in an antimicrobial defense process, NETs also participate in various pathologies, such as autoimmune diseases and cancer, and function as a double-edged sword. Most studies suggest that NETs formed in the cancer environment promote tumor progression by enhancing the abilities of cancer cells to invade (14) and proliferate (17), and NETs awaken dormant tumor cells (18). Additionally, NETs formed in the circulatory system can bind CTCs and thus enhance liver cancer metastasis (19, 20), and those formed in the peritoneal cavity can promote colon metastatic colonization (21). Although NETs form and function as tumor promoters in breast liver and colon cancer, their functions in RCC are unclear.
Since RCC is strongly associated with inflammation (22) and frequently accompanied by tumor thrombus, we hypothesize that a bridge links inflammation and thrombosis and that this bridge also shelters CTCs, promoting their survival. The bridge may be the DNA scaffold decorated with neutrophil-derived TF extruded by circulating neutrophils. Thus, we investigated the relationship between the viable CTC counts and peripheral blood leukocyte statuses of RCC patients.
Materials and Methods
Collection of Patient Specimens and Clinical Information
In this double-blinded study,106 treatment-naïve patients with renal masses who were going to be treated with surgery at the Cancer Institute and Hospital of the Chinese Academy of Medical Sciences (CAMS) between December 2015 and April 2017 were enrolled. Of the 106 enrolled patients, 95 patients were eventually diagnosed with RCC, 9 patients were diagnosed with benign renal conditions according to the final pathology results, and 2 patients who did not undergo surgery were diagnosed with RCC based on the follow-up results. The RCC stage was designated according to the American Joint Committee on Cancer staging system (23). The clinical and histopathological characteristics of the enrolled patients are shown in Table 1. Additionally, hematological parameters, including neutrophil counts and plasma fibrinogen (FIB) and CRP concentrations, were extracted from laboratory records prior to the operation. The cutoff values for defining “high” and “normal” plasma FIB and CRP levels were determined according to the diagnostic criteria reported by the clinical laboratory at our hospital. 4 ml peripheral blood in K2-EDTA vacuum tubes (BD Biosciences, Franklin Lakes, NJ, USA; cat. no. 367844) for CTC count determination, and RNA were extracted from the remaining peripheral blood for RNA sequencing (RNAseq). This study was approved by the Ethics Committee of the Cancer Institute and Hospital of the CAMS.
Table 1. Patient demography.
Counting and Validation of CTCs With the oHSV1-hTERT-GFP Method
The preoperative CTC counts in 106 patients were measured by laboratory personnel blinded to the diagnosis by the oHSV1-hTERT-GFP method; the detailed protocol for this method is available in publications by Zhang and Wang (24, 25). To verify the feasibility of this CTC detection method in RCC patients, we imaged CTCs using an ImageStream® X Mark II imaging flow cytometer (Merck Millipore, Darmstadt, Germany) with a second marker, a human carbonic anhydrase IX/CA9 APC-conjugated antibody (R&D Systems, Minneapolis, MN, USA; cat. no. FAB2188A) (25).
Next-Generation Sequencing and Data Processing
Peripheral blood leukocytes from 21 treatment-naïve patients were obtained after platelets were removed by centrifugation, and red blood cells (RBCs) were lysed with erythrocyte lysis buffer (Qiagen, Hilden, Germany; cat. no. 79217). Subsequently, total RNA was extracted from the peripheral blood leukocytes using TRIzol reagent (Invitrogen, Carlsbad, CA, USA; cat. no. 15596018), and the RNA quality was verified by capillary gel electrophoresis on a Bioanalyzer 2,100 instrument using an Agilent RNA 6,000 Nano Kit (Agilent, CA, USA; cat. no. 5067–1511). Next-generation sequencing (RNAseq) libraries were prepared from qualified samples using the NEBNext® Ultra™ RNA Library Prep Kit for Illumina (New England Biolabs, Ipswich, MA, UK; cat. no. E7530L) according to the manufacturer's instructions. Then, eligible libraries were clustered and sequenced on an Illumina HiSeq X 10 platform (150-bp paired-end reads) using HiSeq X Reagent Kits (Illumina, San Diego, CA, USA; cat. no. FC-501-2521).
Adapter sequences and low-quality reads were removed from the fastq files using Cutadapt (26) and Sickle (http://github.com/najoshi/sickle/) software. Clean reads in fastq files were quantified against an Ensembl catalog (GRCh37) at the transcript level using Salmon (27) and aggregated to the gene level using tximport (28); transcripts per million reads (TPM) values were obtained with Salmon gene expression quantification software, and unexpressed genes were removed. The Pearson correlation coefficients between the log2 transformed TPM values and CTC counts were calculated. Genes with a correlation coefficient > 0.4 and a p-value < 0.05 were defined as being positively correlated with the CTC count.
Another dataset including the mRNA expression data (RNAseqV2, RSEM) and corresponding clinical information for patients with clear cell renal cell carcinoma (ccRCC) (n = 533) was downloaded from The Cancer Genome Atlas (TCGA) (https://tcga-data.nci.nih.gov/tcga/). Subsequent analysis was performed with the log2 transformed mRNA expression values.
Immune-Related Gene Sets
Gene sets, including gene signatures used for decomposition of the 24 immune cell types, angiogenesis marker genes and antigen presentation signatures, were obtained from a publication by Senbabaoglu et al. (29).Combining known key molecules in the formation of NETs and the genes positively correlated with CTC counts, we defined a gene expression signature for NETs (Table S1) (13, 15, 16, 30–32). In addition, the signature for the coagulation cascade is shown in Table S1. The c7 immunological gene sets were downloaded from the Molecular Signatures Database v6.0 (MSigDB).
Decomposition of Immune Cell Types (Immunophenoscore)
We used three methods to identify immune cell types and immune phenotypes in peripheral blood. The relative counts of 24 immune cells and their subtypes, the angiogenesis and coagulation activity levels, the antigen presentation conditions, and NET formation were evaluated using the single-sample gene set enrichment analysis (ssGSEA) method and the gene set variation analysis (GSVA) method implemented in the R package “GSVA” (33). The counts and activities of NETs calculated with the ssGSEA method are hereafter referred to as the NET score. The above methods produced scores for each immune cell type and phenotype in every sample, hereafter referred to as the immunophenoscore. In addition, we analyzed the immune phenotypes of the 21 peripheral blood transcriptomes by the GSEA method considering their associated CTC counts as continuous labels with the c7 immunological gene sets (34).
RCC Cell Line
ACHN cells, which were purchased from the National Infrastructure of Cell Line Resource in China, were cultivated in RPMI-1640 medium (HyClone Laboratories; SH30809.01B) supplemented with 10% fetal bovine serum (FBS; Gibco; 10099–141) and 100 units/mL penicillin and streptomycin.
Detection of NETs in RCC Patients
Peripheral venous blood (4 mL) from RCC patients was collected into K2-EDTA vacuum tubes (BD Biosciences; 367844). Polymorphonuclear cells (PMNs), mainly neutrophils, were isolated using Polymorphprep (Axis-Shield), and the purity of neutrophils was evaluated as previously described (35). Isolated neutrophils were then resuspended in serum-free RPMI-1640 medium and plated on poly-L-lysine-coated chamber coverslips (Thermo Fisher Scientific; 88–8824) with or without LPS (100 ng/mL) stimulation for 4 h. Next, cells were fixed, permeabilized, blocked and stained with an anti-MPO antibody [2C7] (1:400 dilution, Abcam; ab25989), anti-histone H3 (citrulline R2 + R8 + R17) antibody (1:400 dilution, Abcam; ab5103), anti-TF antibody (1:100 dilution, Abcam; ab48647), Hoechst 33342 (1:1,000 dilution, Invitrogen; H3570), and the appropriate secondary antibody, either goat anti-mouse IgG (Alexa 488) (1:2,000 dilution, Abcam; ab150113) or goat anti-rabbit IgG (Alexa 568) (1:2,000 dilution, Abcam; ab175471).
ACHN cells were stained with PKH26 (Sigma; MIDI26), a general cell membrane-labeling probe that binds lipids, according to the manufacturer's instructions. Subsequently, stained ACHN cells (3,000 cells) were spiked into neutrophils isolated from 4 mL of peripheral venous blood. ACHN cells and neutrophils were cocultured and stained with anti-MPO antibody and Hoechst 33342 according to the abovementioned protocols. Stained chamber coverslips were observed with a confocal laser scanning microscope (Leica TCS SP8).
Additional Bioinformatic and Statistical Analyses
To distinguish the differentially expressed immune cells between the two groups (stratified by the concentration of plasma CRP, details in Table 1), linear models and empirical Bayes methods were applied with the R package “limma.” Immune cells with an absolute fold change > 1.2 and a p-value < 0.05 were defined as differentially active immune cells. Gene Ontology (GO) analysis and hierarchy relation analyses were conducted with the R package “topGO” (36). Gene networks were constructed in Cytoscape 3.6.1 with the STRING plugin (37, 38). The Wilcoxon rank test and Pearson and Spearman correlations were implemented by the R package “stat.” Moreover, survival curves were generated by the Kaplan-Meier method and evaluated with the log-rank test (two-tailed). The hazard ratio was determined by Cox multivariate analysis with the R package “survival.” All tests were two-sided, and p < 0.05 indicated statistical significance.
Results
Hematological Parameters Suggested That the Viable CTC Counts Were Related to the Inflammation and Coagulation Statuses of the Patients
We investigated the preoperative viable CTC counts and other clinical characteristics of the 106 enrolled patients (the detailed demographic characteristics of the patients are shown in Table 1) and confirmed that the oHSV1-hTERT-GFP method was feasible in RCC patients (Figure S1). We observed that the CTC counts correlated with the plasma concentrations of CRP and FIB (Figure 1). Specifically, the CTC counts in the group with elevated CRP levels were higher than those in the group with normal CRP levels (Wilcoxon test, p = 0.0104), and both variables displayed a positive correlation (Spearman correlation analysis, R = 0.29, p = 0.0072) (Figure 1A). FIB, a terminal coagulation factor in the coagulation cascade, exhibited a similar pattern; the CTC counts in the group with elevated FIB levels were higher than those in the group with normal FIB levels (Wilcoxon test, p = 0.0124), and higher FIB levels were correlated with higher viable CTC counts (Pearson correlation analysis, R = 0.32, p = 0.0042) (Figure 1B). In addition, CRP and FIB exhibited a positive correlation (Spearman correlation analysis, R = 0.73, p = 6.1 × 10−14) (Figure 1C). These results indicated a probable link between inflammation and coagulation. Moreover, the plasma concentrations of CRP and FIB were positively correlated with not only the viable CTC count but also with the percentage of circulating neutrophils (Figure 1D). Thus, neutrophils may play a crucial role in these intricate relationships. The relationships between viable CTC counts and other clinical characteristics, such as the pathological subtype, stage, presence of microscopic tumor embolus, and other hematological parameters, are shown in Figure S2.
Figure 1. Correlation of hematological parameters (CRP and FIB) with the viable CTC count. CTC counts were increased in the group with elevated circulating CRP levels (A) and the group with elevated circulating FIB levels (B), and both had a positive correlation. The plasma CRP and FIB levels showed a strong positive relationship (C). The circulating neutrophil count was positively correlated with the CTC count and with the plasma CRP and FIB levels. However, the circulating lymphocyte count was negatively correlated with the CTC count (D).
Collectively, these results lead to the hypothesis that inflammation and hypercoagulability are favorable for CTC survival and that neutrophils act as a critical node linking inflammation and coagulation.
Neutrophil Activation Was Related to CTC Survival
To investigate whether neutrophils acted as an assistant and linker in this process, transcriptomic analysis was performed on peripheral blood leukocytes from 21 treatment-naïve patients. At the molecular level, we explored the relationship between the leukocyte transcriptome and the viable CTC count and identified 257 genes positively correlated with the CTC count (see “Next-generation sequencing and data processing” within the Materials and Methods section for details; Figure 2A and Table S2). In addition, the 21 patients were divided into two groups by unsupervised clustering of the 257 genes (Figure 2B). Specifically, members of the group with higher expression levels of these 257 genes tended to have higher CTC counts and plasma CRP levels, but this group included patients in both early and late clinical stages of disease (Figure 2B). To annotate the functions enriched in the 257 genes of interest, GO analysis was performed. The 20 most significantly enriched biological process GO terms were associated with the innate immune response (particularly neutrophil degranulation and phagocytosis), extracellular matrix organization and the integrin-mediated signaling pathway (Figure 2C). Additionally, the GO hierarchy chart of the 10 most significantly enriched biological process terms converged into one term: neutrophil degranulation (GO: 0043312, Figure 2D). Furthermore, we extracted the maximal subnetwork constructed from the 257 genes of interest and highlighted the nearest neighbors of the ELANE gene, such as MPO and CTSG. The ELANE, MPO and CTSG genes encode elastase, MPO, and cathepsin G, respectively. These enzymes are not only released from primary azurophils during degranulation (39) but also decorate the meshwork of chromatin fibers in NETs (40) (Figure 2E). Thus, activated neutrophils are, to some extent, correlated with the survival of CTCs.
Figure 2. Gene profiles positively correlated with the CTC count. (A) Schematic of the discovery pipeline for genes of interest. (B) Heat map of the genes of interest. The heat map displays the log2 transformed TPM values of the 257 genes of interest. The x-axis shows the samples ordered by hierarchical clustering with Ward's linkage. The top bar on the x-axis represents the clinical stage of the patients, the second bar indicates the measured viable CTC count in 4 mL of preoperative peripheral blood, the third bar indicates the group stratified by the CRP concentration (cutoff of 0.6 mg/dL), and the fourth bar indicates the log2 transformed CRP concentration. (C) The 20 most significant GO terms enriched in the genes of interest. The vertical axis shows the number of positively correlated significantly enriched genes for each term. The horizontal axis displays the names of the GO terms ordered by the degree of significance. The color indicates the–log10 value (weighted Fisher's p-value). (D) Hierarchy chart of the 10 most significantly enriched GO biological process terms. The nodes are displayed in a circle or box shape and include information on the GO ID, GO name, weighted Fisher's p-values, and counts of genes of interest/total count of genes in the term. The box represents the 10 most significantly enriched terms. The colors of the box indicate the significance and range from dark red (most significant) to light yellow (least significant). The black arrows indicate is-a relationships. (E) The maximal subnetwork constructed with the genes of interest. The green nodes represent the nearest neighbors of ELANE. The thickness and color of the edges indicate the coexpression values; lower values have lower thicknesses and darker colors.
NET Formation Cooperated With Coagulation to Promote CTC Survival
Because neutrophil-mediated inflammation was shown to promote CTC survival, we sought to determine whether NETs induced by circulating neutrophils were conducive to CTC metastasis. Indeed, the CTC count and expression value of each key molecule in NET production were positively correlated (Figure S3); particularly PAD4, an essential enzyme mediating histone hypercitrullination during NET formation (30). To further confirm this finding, we explored the data with another approach. When the CTC count was treated as a continuous variable, GSEA showed that the leukocytes in patients with high CTC counts were similar to S. pneumoniae rather than E. coli infection or to G-CSF-treated peripheral blood mononuclear cells (PBMCs) (Figure 3A). The S. pneumoniae infection signature suggested that low-density neutrophils (LDNs) copurified with mononuclear cells during density gradient centrifugation were activated (41). It has been acknowledged that unstimulated LDNs can produce NETs spontaneously (17, 42), and G-CSF can prime neutrophils to generate NETs (13). Additionally, we used gene set-based scoring systems (immunophenoscore obtained using the ssGSEA method; details are shown in the Materials and Methods section) to evaluate the functions and counts of physiological processes and immune cells. Interestingly, among the whole set of immune cells and their specific physiological processes, NET formation exhibited the strongest positive correlation with the viable CTC count (Pearson correlation analysis, R = 0.69, p = 6 × 10−4, Figures 3B,C). Increased neutrophil activity, particularly the NET formation mechanism, was unanimously associated with an increased viable CTC count.
Figure 3. Relationship among the NET score, FIB level, and CTC count. (A) GSEA of the peripheral blood leukocyte RNAseq dataset. GSEA was performed with continuous-label CTC counts. The GSEA mountain plots exhibited the 2 most divergent signatures (with a minimum enrichment score of 0.6 and a false discovery rate (FDR) cutoff of 0.05). The gene set names were refined to fit the figure. (B) Relationship between the CTC count and the immunophenoscore as assessed by ssGSEA for immune cell subgroups and selected physiological processes. The vertical axis displays the correlation coefficients of the CTC counts and each subgroup of immune cells. In addition, the colors of the bars correspond to the correlation coefficient. (C) Relationship between the CTC count and the NET score. (D) Relationship between the NET score and the FIB concentration. (E) Plot of the enriched (yellow) and depleted (blue) immune phenotypes between the group with elevated CRP levels and the group with normal CRP levels based on the immunophenoscore (GSVA).
Moreover, the NET score was positively correlated with the plasma concentration of FIB (Pearson correlation analysis, R = 0.49, p = 0.033, Figure 3D). Not only the transcriptomic data but also the image of immunofluorescence staining showed the similar result (Figure 4B) that NET-derived TF adheres to the DNA framework of NETs. Because TF is a factor that triggers an extrinsic coagulation cascade, eventually leading to crosslinking with fibrin, the defined list of genes involved in the coagulation system was composed of the key coagulation factors in the extrinsic coagulation pathway (details in the Materials and Methods section). The volcano plot of the immunophenoscores (GSVA scores) showed that the coagulation system and NET formation were stimulated consistently in patients with elevated CRP levels (Figure 3E). Thus, NETs may be the link between inflammation and coagulation.
Figure 4. Formation of NETs in the peripheral blood of patients with ccRCC. (A) DNA, MPO, and H3Cit were colocalized with NETs formed in unstimulated and LPS-treated blood. (B) DNA, MPO, and TF were colocalized with NETs formed in LPS-treated blood. Neutrophil-derived TF was produced during NET formation. (C) NETs were entangled with ACHN cells in the blood of ccRCC patients.
NETs Were Entangled With and Sheltered Spiked CTCs in the Circulatory System of RCC Patients
Next, we confirmed that NETs were morphologically detectable in unstimulated and LPS-treated peripheral blood leukocytes of RCC patients (Figure 4A). NETs were colocalized with released MPO and citrullinated histone 3 (H3Cit), which adhered to the extruded DNA scaffold. We also confirmed via immunofluorescence staining that the framework of NETs was decorated with NET-derived TF (Figure 4B). Thus, NETs function as a bridge to link inflammation and coagulation. Moreover, we observed NETs entangled with the spiked ACHN cells, which were labeled with PKH26, in both untreated and LPS-stimulated peripheral blood leukocytes of ccRCC patients (Figure 4C). Collectively, these results indicate that NETs induced by circulating neutrophils could form a protective shelter and vehicle conducive to the metastasis of CTCs.
Inflammation and Coagulation Indexes in the Tumor Microenvironment Are Independent Prognostic Factors in RCC Patients
Although persistent inflammation and hypercoagulability increased the counts of viable metastatic seeds (CTCs) in the circulatory system, we sought to determine whether these factors contribute to the generation of CTCs or the microenvironment of primary tumors. To this end, we investigated RNAseq data in the TCGA database from 533 patients with ccRCC. A marked separation of the survival curves was confirmed regardless of whether the samples were stratified by quartiles or median NET scores (p = 0.00056 and 0.0359, respectively). Furthermore, the NET score was an independent prognostic factor when the model was adjusted for clinical stage (Figure 5A). Interestingly, the triggering factor of the extrinsic coagulation pathway, TF, could separate the survival curves and was an independent prognostic factor (Figure 5B). In addition, the NET scores and TF expression values were positively correlated (Figure 5C). Via immunofluorescence, neutrophils were also observed to release NETs decorated with F3 (Figure 4B). In a sense, both NET-derived TF-mediated fibrin crosslinking and the DNA scaffold released from NETs create a secure scaffold for cancer cell invasion.
Figure 5. NETs and F3 in tissues are independent prognostic factors in ccRCC patients. (A) Kaplan-Meier plots of the overall survival of ccRCC patients (n = 533) from the TCGA database stratified by the NET score quartile (left, p = 0.000566) and median NET score (right, p = 0.0359). The bottom list shows the adjusted HR for NETs by clinical stage (HR = 16.18). (B) Kaplan-Meier plots of the overall survival of ccRCC patients (n = 533) from the TCGA database stratified by the TF expression quartile (left, p = 0.49 × 10−5) and median TF expression value (right, p = 0.0134). The bottom list shows the adjusted HR for TF expression by clinical stage (HR = 1.12). (C) The NET score and TF expression value were positively correlated.
Discussion
Most RCC patients die from cancer that relapses after CTCs spread to other organs rather than from early localized carcinoma. However, to date, little is known about the mechanisms by which CTCs survive in the circulatory system to seed other inhospitable tissues. Based on our research results, we propose a potential process by which circulating renal cancer cells metastasize to other organs from the perspective of dysregulated inflammation and coagulation.
NETs are the bridge connecting inflammation and coagulation, and the NET DNA scaffold not only provides a location for platelet and red blood cell adhesion (43) but is also decorated with TF (11). In a functioning host, the immune system cooperates with the coagulation system for infection defense or trauma repair to maintain homeostasis. Among the diverse defense mechanisms, immunothrombosis generated by immune cells and specific thrombosis-related molecules generate an intravascular scaffold that helps to defend against exogenous pathogens (44). However, cancer is known as the wound that does not heal (45). Thus, uncontrolled immunothrombosis in a cancer-burdened host might become a major pathological process fostering foreign CTCs.
Our data indicated that peripheral blood leukocytes of RCC patients tend to be activated by a specific process, similar to the spontaneous activation of LDNs or G-CSF priming in vitro—namely, the formation of NETs—particularly in samples with a higher viable CTC count. In addition, our results indicated that increased NET formation was accompanied by increased plasma levels of FIB, which have been reportedly associated with poor prognosis due to NK cell inhibition (46, 47), and several studies have reported that NET formation promotes fibrin and thrombin generation (48–51). Thus, intravascular neutrophils form TF-decorated tethers and slings to activate the coagulation system, simultaneously anchoring CTCs to these NETs as a protective vehicle.
Moreover, intravascular NETs not only protect CTCs from attack but also promote their extravasation, as dysregulated NET formation contributes to inflammatory vascular injury, endothelial cell shrinkage and tissue damage (50, 52), and NET-remodeled laminin activates integrin α3β1 signaling (18), which was a significantly enriched GO term associated with the 257 genes of interest herein (Figure 2C). NETs and their decorated molecules (e.g., TF), which were herein demonstrated to be indicators of poor prognosis in the TCGA ccRCC RNAseq data (Figure 5), may cooperate to facilitate cancer cell penetration into the vasculature. In addition, NETs have been shown to promote angiogenesis in vitro (49), but this phenomenon has not yet been observed in vivo.
Additionally, patients with microscopically visible microtumor thrombosis exhibited higher numbers of viable single CTCs than patients without visible thrombosis (Figure S2C). This finding, combined with the abovementioned information, revealed that large renal venous tumor thromboses originate from a single CTC captured by NETs and accumulate gradually. However, the hypothesis of venous tumor-immunothrombosis formation needs further validation. Except for venous tumor-immunothrombosis, research has suggested that the appearance of CTCs is a risk factor for venous thromboembolism (VTE) in early and metastatic breast cancer (53, 54). Similar to our result, viable CTC counts and plasma FIB displayed a positive correlation (Figure 1B).
RCC is not only complicated by thrombus but is also strongly associated with inflammation, and elevated levels of plasma CRP, an indicator of inflammation, predict poor survival in patients with localized RCC (55) and metastatic RCC (56). However, the mechanisms driving the association of CRP, an acute phase reactant protein produced by the liver, with cancer-specific survival are poorly understood. CRP levels have been reportedly correlated with the production of interleukin-6 (IL-6) (57), which mediates ongoing inflammatory processes. We showed that both persistent inflammation and the coagulation system were activated in patients with elevated CRP levels (Figure 3E). Thus, the presence of CRP, a non-specific indicator of inflammation, may represent a microenvironment suitable for NET formation and coagulation system activation, which is conducive to the intravascular survival of CTCs, metastatic seeds. Conversely, fewer alive CTCs were detected when peripheral blood dendritic cells (DCs) and antigen-presenting machinery were activated (Figure 3B). Consistent with conventional cognition, both the activated protumor immune response and the compromised function of the antitumor immune response (such as DCs) support the survival of CTCs (58).
In this study, we did not identify a direct relationship between the CTC count and prognosis due to the lack of complete follow-up information for the enrolled patients. We observed only spiked cancer cells entangled in NET tethers in untreated ccRCC patients (Figure 4B). We will attempt to elucidate the process by which immunothrombosis (NET entanglement with the fibrin network) is connected to CTCs in vivo in the future.
In conclusion, we proposed a survival mechanism of intravascular CTCs—intravascular neutrophils form TF-decorated tethers and slings, namely, NETs, to activate the coagulation system and further anchor CTCs. Thus, degrading this framework (the double shelter of CTCs—the DNA scaffold and the fibrin network initiated by NET-derived TF) may be a potential way to inhibit metastasis.
Data Availability
The data for the 21 raw RNA sequences generated for this study can be found in the Genome Sequence Archive (59) at the BIG Data Center (60), Beijing Institute of Genomics (BIG), Chines Academy of Sciences, under accession number HRA000042. [http://bigd.big.ac.cn/gsa-human]. Data download needs to be in agreement with the corresponding author.
Ethics Statement
This study was carried out in accordance with the recommendations of the Ethics Committee of the Cancer Institute and Hospital of the Chinese Academy of Medical Sciences with written informed consent from all subjects. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by the Ethics Committee of the Cancer Institute and Hospital of the Chinese Academy of Medical Sciences.
Author Contributions
JS, LF, and KZ contributed to the conception and design of the study. LW, YL and WJ collected the samples and clinical information. LG performed the experiments. XD provided instruction for the experiments. LF and LG performed the statistical analysis. LG wrote the first draft of the manuscript. LF, WZ, JM, KZ, and JS wrote sections of the manuscript. All authors contributed to the manuscript revision and read and approved the submitted version.
Funding
This work was supported by the Chinese Academy of Medical Sciences (CAMS) Initiative for Innovative Medicine under grant number 2016-I2M-1-007, the Chinese Academy of Medical Sciences (CAMS) Initiative for Innovative Medicine under grant number 2017-I2M-2-003, and the PUMC Fund of the Funds for the Central Universities (3332016027).
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.
Acknowledgments
We would like to thank the patients who were enrolled in this study. We also thank Dr. Qingya Yang of Qilu Hospital of Shandong University for helping with this study. We thank the staff of the National Center for Protein Sciences Beijing (Peking University) and Ms. Fei Wang for her technical assistance with the flow cytometry imaging.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2019.00504/full#supplementary-material
References
1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2018) 68:394–424. doi: 10.3322/caac.21492
2. Choueiri TK, Motzer RJ. Systemic therapy for metastatic renal-cell carcinoma. N Engl J Med. (2017) 376:354–66. doi: 10.1056/NEJMra1601333
3. Robert J, Motzer MD, Thomas E, Hutson DO, Pharm D, Piotr Tomczak MD, et al. Sunitinib versus interferon alfa in metastatic renal-cell carcinoma. N. Engl. J. Med. (2007) 356:115–24. doi: 10.1056/NEJMoa065044
4. Massagué J, Obenauf AC. Metastatic colonization by circulating tumour cells. Nature. (2016) 529:298–306. doi: 10.1038/nature17038
5. Steinert G, Schölch S, Niemietz T, Iwata N, García SA, Behrens B, et al. Immune escape and survival mechanisms in circulating tumor cells of colorectal cancer. Cancer Res. (2014) 74:1694–704. doi: 10.1158/0008-5472.CAN-13-1885
6. Labelle M, Begum S, Hynes RO. Direct signaling between platelets and cancer cells induces an epithelial-mesenchymal-like transition and promotes metastasis. Cancer Cell. (2011) 20:576–90. doi: 10.1016/j.ccr.2011.09.009
7. Spicer JD, McDonald B, Cools-Lartigue JJ, Chow SC, Giannias B, Kubes P, et al. Neutrophils promote liver metastasis via Mac-1–mediated interactions with circulating tumor cells. Cancer Res. (2012) 72:3919–27. doi: 10.1158/0008-5472.CAN-11-2393
8. Ito K, Asano T, Yoshii H, Satoh A, Sumitomo M, Hayakawa M. Impact of thrombocytosis and C-reactive protein elevation on the prognosis for patients with renal cell carcinoma. Int J Urol. (2006) 13:1365–70. doi: 10.1111/j.1442-2042.2006.01563.x
9. Cools-Lartigue J, Spicer J, Najmeh S, Ferri L. Neutrophil extracellular traps in cancer progression. Cell Mol Life Sci. (2014) 71:4179–94. doi: 10.1007/s00018-014-1683-3
10. Wu L, Saxena S, Awaji M, Singh RK. Tumor-associated neutrophils in cancer: going pro. Cancers. (2019) 11:e564. doi: 10.3390/cancers11040564
11. Kambas K, Mitroulis I, Ritis K. The emerging role of neutrophils in thrombosis-the journey of TF through NETs. Front Immunol. (2012) 3:385. doi: 10.3389/fimmu.2012.00385
12. Volker B, Ulrike R, Christian G, Beatrix F, Yvonne U, Weiss DS, et al. Neutrophil extracellular traps kill bacteria. Science. (2004) 303:1532–5. doi: 10.1126/science.1092385
13. Demers M, Wagner DD. Neutrophil extracellular traps: a new link to cancer-associated thrombosis and potential implications for tumor progression. Oncoimmunology. (2013) 2:e22946. doi: 10.4161/onci.22946
14. Park J, Wysocki RW, Amoozgar Z, Maiorino L, Fein MR, Jorns J, et al. Cancer cells induce metastasis-supporting neutrophil extracellular DNA traps. Sci Transl Med. (2016) 8:361ra138. doi: 10.1126/scitranslmed.aag1711
15. Papayannopoulos V, Metzler KD, Hakkim A, Zychlinsky A. Neutrophil elastase and myeloperoxidase regulate the formation of neutrophil extracellular traps. J Cell Biol. (2010) 191:677–91. doi: 10.1083/jcb.201006052
16. Carmona-Rivera C, Zhao W, Yalavarthi S, Kaplan MJ. Neutrophil extracellular traps induce endothelial dysfunction in systemic lupus erythematosus through the activation of matrix metalloproteinase-2. Ann Rheum Dis. (2015) 74:1417–24. doi: 10.1136/annrheumdis-2013-204837
17. Kanamaru R, Ohzawa H, Miyato H, Matsumoto S, Haruta H, Kurashina K, et al. Low density neutrophils (LDN) in postoperative abdominal cavity assist the peritoneal recurrence through the production of neutrophil extracellular traps (NETs). Sci. Rep. (2018) 8:632. doi: 10.1038/s41598-017-19091-2
18. Albrengues J, Shields MA, Ng D, Park CG, Ambrico A, Poindexter ME, et al. Neutrophil extracellular traps produced during inflammation awaken dormant cancer cells in mice. Science. (2018) 361:eaao4227. doi: 10.1126/science.aao4227
19. Tohme S, Yazadani HO, Alkhafaji AB, Chidi AP, Lougharn P, Mowen K, et al. Neutrophil extracellular traps promote the development and progression of liver metastases after surgical stress. Cancer Res. (2016) 76:1367. doi: 10.1158/0008-5472.CAN-15-1591
20. Najmeh S, Cools-Lartigue J, Rayes RF, Gowing S, Vourtzoumis P, Bourdeau F, et al. Neutrophil extracellular traps sequester circulating tumor cells via it-integrin mediated interactions. Int J Cancer. (2017) 140:2321–30. doi: 10.1002/ijc.30635
21. Al-Haidari AA, Algethami N, Lepsenyi M, Rahman M, Syk I, Thorlacius H. Neutrophil extracellular traps promote peritoneal metastasis of colon cancer cells. Oncotarget. (2019) 10:1238–49. doi: 10.18632/oncotarget.26664
22. Blay JY, Rossi JF, Wijdenes J, Menetrier-Caux C, Schemann S, Négrier S, et al. Role of interleukin-6 in the paraneoplastic inflammatory syndrome associated with renal-cell carcinoma. (1997) 72:424–30. doi: 10.1002/(SICI)1097-0215(19970729)72:3<424::AID-IJC9>3.0.CO;2-R
23. Stephen E, David RB, Carolyn C, April GF, Frederick LG, Andy T, et al. AJCC Cancer Staging Manual. 7 ed. New York, NY: Springer (2010).
24. Zhang W, Bao L, Yang S, Qian Z, Dong M, Yin L, et al. Tumor-selective replication herpes simplex virus-based technology significantly improves clinical detection and prognostication of viable circulating tumor cells. Oncotarget. (2016) 7:39768. doi: 10.18632/oncotarget.9465
25. Wang Y, Guo L, Feng L, Zhang W, Xiao T, Di X, et al. Single nucleotide variant profiles of viable single circulating tumour cells reveal CTC behaviours in breast cancer. Oncol Rep. (2018) 39:2147–59. doi: 10.3892/or.2018.6325
26. Martin M. Cutadapt removes adapte sequences from high-throughput sequencing reads. Embnet J. (2011) 17:10–12. doi: 10.14806/ej.17.1.200
27. Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. (2017) 14:417–9. doi: 10.1038/nmeth.4197
28. Soneson C, Love MI, Robinson MD. Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences. F1000Research. (2015) 4:1521. doi: 10.12688/f1000research.7563.1
29. Senbabaoglu Y, Gejman RS, Winer AG, Liu M, Van Allen EM, de Velasco G, et al. Tumor immune microenvironment characterization in clear cell renal cell carcinoma identifies prognostic and immunotherapeutically relevant messenger RNA signatures. Genome Biol. (2016) 17:231. doi: 10.1186/s13059-016-1092-z
30. Leshner M, Shu W, Lewis C, Han Z, Chen XA, Santy L, et al. PAD4 mediated histone hypercitrullination induces heterochromatin decondensation and chromatin unfolding to form neutrophil extracellular trap-like structures. Front Immunol. (2012) 3:307. doi: 10.3389/fimmu.2012.00307
31. Kuchler R, Schroeder BO, Jaeger SU, Stange EF, Wehkamp J. Antimicrobial activity of high-mobility-group box 2: a new function to a well-known protein. Antimicrob Agents Chemother. (2013) 57:4782–93. doi: 10.1128/AAC.00805-13
32. Tadie J-M, Bae H-B, Jiang S, Park DW, Bell CP, Yang H, et al. HMGB1 promotes neutrophil extracellular trap formation through interactions with Toll-like receptor 4. Am J Physiol Lung Cell Mol Physiol. (2013) 304:L342–9. doi: 10.1152/ajplung.00151.2012
33. Hänzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-Seq data. BMC Bioinform. (2013) 14:7. doi: 10.1186/1471-2105-14-7
34. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. (2005) 102:15545–50. doi: 10.1073/pnas.0506580102
35. Guo L, Chen G, Zhang W, Zhou L, Xiao T, Di X, et al. A high-risk luminal A dominant breast cancer subtype with increased mobility. Breast Cancer Res Treat. (2019) 175:459–72. doi: 10.1007/s10549-019-05135-w
36. Alexa A, Rahnenfuhrer J, Lengauer T. Improved scoring of functional groups from gene expression data by decorrelating GO graph structure. Bioinformatics. (2006) 22:1600–7. doi: 10.1093/bioinformatics/btl140
37. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. (2003) 13:2498. doi: 10.1101/gr.1239303
38. Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein–protein interaction networks, integrated over the tree of life. Nucleic Acids Res. (2015) 43:D447–52. doi: 10.1093/nar/gku1003
39. Coffelt SB, Wellenstein MD, de Visser KE. Neutrophils in cancer: neutral no more. Nat Rev Cancer. (2016) 16:431–46. doi: 10.1038/nrc.2016.52
40. Kaplan MJ, Radic M. Neutrophil extracellular traps: double-edged swords of innate immunity. J Immunol. (2012) 189:2689–95. doi: 10.4049/jimmunol.1201719
41. Ramilo O, Allman W, Chung W, Mejias A, Ardura M, Glaser C, et al. Gene expression patterns in blood leukocytes discriminate patients with acute infections. Blood. (2007) 109:2066–77. doi: 10.1182/blood-2006-02-002477
42. Lood C, Blanco LP, Purmalek MM, Carmonarivera C, Ravin SSD, Smith CK, et al. Neutrophil extracellular traps enriched in oxidized mitochondrial DNA are interferogenic and contribute to lupus-like disease. Nat Med. (2016) 22:146–53. doi: 10.1038/nm.4027
43. Fuchs TA, Brill A, Duerschmied D, Schatzberg D, Monestier M, Myers DD, et al. Extracellular DNA traps promote thrombosis. Proc Natl Acad Sci USA. (2010) 107:15880–5. doi: 10.1073/pnas.1005743107
44. Engelmann B, Massberg S. Thrombosis as an intravascular effector of innate immunity. Nat Rev Immunol. (2013) 13:34–45. doi: 10.1038/nri3345
45. Dvorak HF. Tumors: wounds that do not heal-redux. Cancer Immunol Res. (2015) 3:1–11. doi: 10.1158/2326-6066.CIR-14-0209
46. Palumbo JS, Talmage KE, Massari JV, Jeunesse CML, Flick MJ, Kombrinck KW, et al. Platelets and fibrin(ogen) increase metastatic potential by impeding natural killer cell. Blood. (2005) 105:178–85. doi: 10.1182/blood-2004-06-2272
47. Pichler M, Hutterer GC, Stojakovic T, Mannweiler S, Pummer K, Zigeuner R. High plasma fibrinogen level represents an independent negative prognostic factor regarding cancer-specific, metastasis-free, as well as overall survival in a European cohort of non-metastatic renal cell carcinoma patients. Br J Cancer. (2013) 109:1123–9. doi: 10.1038/bjc.2013.443
48. Hisada Y, Grover SP, Maqsood A, Houston R, Ay C, Noubouossie DF, et al. Neutrophils and neutrophil extracellular traps enhance venous thrombosis in mice bearing human pancreatic tumors. Haematologica. (2019). doi: 10.3324/haematol.2019.217083. [Epub ahead of print].
49. Jung HS, Gu J, Kim JE, Nam Y, Song JW, Kim HK. Cancer cell-induced neutrophil extracellular traps promote both hypercoagulability and cancer progression. PLoS ONE. (2019) 14:e0216055. doi: 10.1371/journal.pone.0216055
50. Li B, Liu Y, Hu T, Zhang Y, Zhang C, Li T, et al. Neutrophil extracellular traps enhance procoagulant activity in patients with oral squamous cell carcinoma. J Cancer Res Clin Oncol. (2019) 132:2399. doi: 10.1007/s00432-019-02922-2
51. Seo JD, Gu JY, Jung HS, Kim YJ, Kim HK. Contact system activation and neutrophil extracellular trap markers: risk factors for portal vein thrombosis in patients with hepatocellular carcinoma. Clin Appl Thromb Hemost. (2019) 25:1076029618825310. doi: 10.1177/1076029618825310
52. McInturff AM, Cody MJ, Elliott EA, Glenn JW, Rowley JW, Rondina MT, et al. Mammalian target of rapamycin regulates neutrophil extracellular trap formation via induction of hypoxia-inducible factor 1 alpha. Blood. (2012) 120:3118–25. doi: 10.1182/blood-2012-01-405993
53. Mego M, Karaba M, Minarik G, Benca J, Sedlackova T, Tothova L, et al. Relationship between circulating tumor cells, blood coagulation, and urokinase-plasminogen-activator system in early breast cancer patients. Breast J. (2015) 21:155–60. doi: 10.1111/tbj.12388
54. Mego M, Zuo Z, Gao H, Cohen EN, Giordano A, Tin S, et al. Circulating tumour cells are linked to plasma D-dimer levels in patients with metastatic breast cancer. Thromb Haemost. (2015) 113:593–8. doi: 10.1160/TH14-07-0597
55. Komai Y, Saito K, Sakai K, Morimoto S. Increased preoperative serum C-reactive protein level predicts a poor prognosis in patients with localized renal cell carcinoma. BJU Int. (2007) 99:77–80. doi: 10.1111/j.1464-410X.2006.06497.x
56. Fujikawa K, Matsui Y, Oka H, Fukuzawa S, Takeuchi H. Serum C-reactive protein level and the impact of cytoreductive surgery in patients with metastatic renal cell carcinoma. J Urol. (1999) 162:1934–7. doi: 10.1016/S0022-5347(05)68072-X
57. Blay JY, Negrier S, Combaret V, Attali S, Goillot E, Merrouche Y, et al. Serum level of interleukin 6 as a prognosis factor in metastatic renal cell carcinoma. Cancer Res. (1992) 52:3317–22.
58. Mego M, Gao H, Cohen EN, Anfossi S, Giordano A, Tin S, et al. Circulating tumor cells (CTCs) are associated with abnormalities in peripheral blood dendritic cells in patients with inflammatory breast cancer. Oncotarget. (2017) 8:35656–68. doi: 10.18632/oncotarget.10290
59. Wang Y, Song F, Zhu J, Zhang S, Yang Y, Chen T, et al. GSA: genome sequence archive < sup/>. Genomics Proteomics Bioinform. (2017) 15:14–8. doi: 10.1016/j.gpb.2017.01.001
Keywords: renal cell carcinoma, circulating tumor cells, neutrophil extracellular trap, hypercoagulability, C-reactive protein
Citation: Wen L, Guo L, Zhang W, Li Y, Jiang W, Di X, Ma J, Feng L, Zhang K and Shou J (2019) Cooperation Between the Inflammation and Coagulation Systems Promotes the Survival of Circulating Tumor Cells in Renal Cell Carcinoma Patients. Front. Oncol. 9:504. doi: 10.3389/fonc.2019.00504
Received: 23 February 2019; Accepted: 28 May 2019;
Published: 17 June 2019.
Edited by:
Janice P. Dutcher, Cancer Research Foundation, United StatesReviewed by:
Michal Mego, Campus Bio-Medico University, ItalyElena Ranieri, University of Foggia, Italy
Copyright © 2019 Wen, Guo, Zhang, Li, Jiang, Di, Ma, Feng, Zhang and Shou. 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: Lin Feng, ZmVuZ2xpbiYjeDAwMDQwO2NpY2Ftcy5hYy5jbg==; Kaitai Zhang, emhhbmdrdCYjeDAwMDQwO2NpY2Ftcy5hYy5jbg==; Jianzhong Shou, c2hvdWp6aCYjeDAwMDQwOzEyNi5jb20=
†These authors have contributed equally to this work