- Primary research
- Open Access
CBFA2T2 is associated with a cancer stem cell state in renal cell carcinoma
Cancer Cell Internationalvolume 17, Article number: 103 (2017)
The Correction to this article has been published in Cancer Cell International 2018 18:73
Renal cell carcinoma (RCC) is the most common kidney cancer, accounting for approximately 80–90% of all primary kidney cancer. Treatment for patients with advanced RCC remains unsatisfactory. Rare cancer stem cells (CSCs) are proposed to be responsible for failure of current treatment.
OncoLnc was used as a tool for interactively exploring survival correlations. Gene manipulation and expression analysis were carried out using siRNA, RT-PCR and Western blotting. Wound healing and invasion assays were used for phenotypical characterization. Aldefluor assay and FACS sorting Sphere culture were used to determine the “stemness” of CSCs. Co-Immunoprecipitation (Co-IP) was used to examine the interaction between OCT4 and CBFA2T2. Student’s t-test and Chi square test was used to analyze statistical significance.
CBFA2T2 expression can significantly predict the survival of RCC patients. Knocking-down of CBFA2T2 can inhibit cell migration and invasion in RCC cells in vitro, and reduce ALDHhigh CSCs populations. CBFA2T2 expression is necessary for sphere-forming ability and cancer stem cells marker expression in RCC cell lines.
Our data suggest that CBFA2T2 expression correlates with aggressive characteristics of RCC and CBFA2T2 is required for maintenance of “stemness” through regulation of stem cells factors, thereby highlighting CBFA2T2 as a potential therapeutic target for RCC treatment.
There are approximately 62,700 new cases of kidney cancer and an estimated 14,240 cancer associated deaths in the United States annually . Renal cell carcinoma (RCC) represents one of the most common type of kidney cancer, accounting for 90% of adult renal malignancies . Over the past five decades, the incidence of RCC has increased at a rate of 2% per year . Surgery still serve as the standard way for localized and locally advanced RCC. Despite the surgical resection, the prognosis of advanced RCC is very poor, with a 5-year survival rate of 5–10% . In addition, around 20–40% of RCC patients developed recurrence after treatment . The mechanisms of how these metastasis and recurrence occur are largely uncharacterized. Recent studies have shown the presence of cancer stem cells (CSCs) in various solid tumor tissues [5,6,7]. CSCs are undifferentiated multipotential tumor cells that have high tumorigenic activity and the capacity to self-renew [5, 6, 8]. Moreover, CSCs promote cancer resistance to treatment and recurrence, leading to high mortality . Therefore, CSC-targeted therapy can be an essential part of cancer therapy.
Previously, CSCs in RCC have been identified by using either stem cell markers [CD133, CD44, CXCR4, CD105 and Spalt-Like Transcription Factor 4 (SALL4)] [9,10,11,12,13] or functional assays (sphere-forming ability, side population and ALDH activity) [14,15,16]. Especially, in RCC cell lines, ALDHhigh cells showed CSC properties in vitro, such as clonogenic and self-renewal ability and increased expression of OCT4 and NANOG . Importantly, in the metastatic RCC cell lines, ALDHhigh cells formed about 15% of the total number of cells and had higher sphere-forming capacity, self-renewal ability, and tumorigenicity than ALDHlow cells. While anther study showed that ALDH gene expression was correlated with tumor grade but not with tumor stage in patients with RCC . Overall pooled analysis suggested that high CSCs markers, including CD133, CD44, CXCR4 and CD105, expression predicted poor overall survival, cancer-specific survival, disease-free survival and progression-free survival . However, the molecular mechanism of how CSCs in RCC are maintained has not been fully explored.
Recently, corepressor CBFA2T2 has been shown to be involved in Prdm14-mediated germ cell formation through oligomerization to form a scaffold, which allows OCT4 to be stabilized on chromatin , indicating a role of CBFA2T2 in stem cell pluripotency regulation. CBFA2T2 has been known as a fusion partner of RUNX1 in acute myeloid leukemia , suggesting a role of CBFA2T2 in cancer development. Indeed, CBFA2T2 plays a key role in tumorigenesis in AOM/DSS colitis-associated carcinoma in mouse . However, the role of CBFA2T2 in CSCs regulation remain unknown.
In this study, we sought to identify the role of CBFA2T2 in human RCC and to determine its characteristics in RCC tumorigenesis and RCC CSCs maintenance. Here, we showed that the survival rate in the RCC patients is negatively correlated CBFA2T2 expression. In addition, we showed the CBFA2T2 silencing impaired wound-healing and invasion capabilities of 786-O and A-498 cells, two human RCC cell lines. Our ALDH staining demonstrated that CBFA2T2 depletion in RCC cell lines led to reduced CSCs proportions. Importantly, CBFA2T2 is required for sphere-forming ability and cancer stem cells marker expression in RCC cell lines. We also showed an interaction between OCT4 and CBFA2T2 by pulldown assay. Collectively, our findings may provide new insights for future CSC study and clinical anti-cancer therapy in RCC.
Association between CBFA2T2 and stem cell factors expression with the prognosis of RCC patients
To examine the association of the expression levels of CBFA2T2, as well as several important cancer stem cells factors, including SOX2, OCT4, CD133, NANOG and ALDH1A3, with the and prognosis of patients with RCC, Kaplan–Meier survival analyses were performed on patients in the high and low expression groups between these gene expression and clinical outcomes on the homepage of oncolnc . As shown in Fig. 1, the survival rate in the RCC patients who exhibited high CBFA2T2 expression was significantly reduced compared with patients with low CBFA2T2 expression by log-rank test (p = 0.000942). In the meantime, higher levels of OCT-4, ALDH1A3, NANOG and SALL4 expression were associated with worse survival (p = 0.048, p = 0.00228, p = 2.46e−05 and p = 3.98e−07, respectively) of RCC patients (Fig. 1b–e). However, the expression levels SOX2 and CD133 showed no significant correlation between survival rate of the RCC patients (data not shown). Indeed, our immunohistochemistry staining showed elevated expression of CBFA2T2 in RCC cancer compared to normal kidney tissue (Additional file 1: Figure S1). Our analysis using The Cancer Genome Atlas (TCGA) showed showing 0.4% of CBFA2T2—altered in RCC samples (Additional file 2: Figure S2). To investigate the relationship between CBFA2T2 expression and RCC further, we analyzed a cohort of 66 RCC patients. RCC were classified into 2 groups based on CBFA2T2 staining score: Low (0–1) and High [2, 3]. Our results indicated that increased expression of CBFA2T2 was associated with high grade RCC significantly (p = 0.023). Moreover, increased expression of CBFA2T2 was associated high stage disease, with a trend toward statistical significance (p = 0.104). However, our data showed that there were no significant differences for gender, lymph node, systematic metastasis and pathological types between the groups (Table 1).
CBFA2T2 is required for migration and invasion of RCC cells in vitro
To investigate whether CBFA2T2 play a role in RCC cancer development, we used siRNA to inhibit CBFA2T2 expression. Knockdown of CBFA2T2 in both 786-O and A-498 cell were confirmed by Western blot analysis (Fig. 2a). Our wound healing assay showed that in silencing of CBFA2T2 in RCC cell lines significantly impaired their migratory abilities (Fig. 2b, c). Control 786-O cells almost completely closed the wound (78.5%) in 20 h, whereas 786-O cells treated with CBFA2T2 siRNA 20 h closed 66.7% of the wound. Similarly, control A-498 cells after 20 h almost completely closed the wound (85.0%), whereas A-498 cells treated with CBFA2T2 siRNA closed only 71% of the wound (Fig. 2b, c). In the invasion assay, we found that CBFA2T2 siRNA groups demonstrated a 31.51 ± 10.86% decrease in cell migration in 786-O cells and a 65.13 ± 9.97% decrease in that of A-498 cells, compared with the controls. These results indicated that CBFA2T2 reduced the migration and invasion in Renal cancer 786-O and A-498 cells.
CBFA2T2 is required for maintenance of cancer stem cell property of RCC
In RCC, ALDHhigh subpopulation as assessed by the Aldefluor assay has been shown to be the cancer stem cells population, which has upregulated levels of OCT-4 and NANOG . Hence, we examined whether knockdown of CBFA2T2 could lead to decrease of ALDHhigh cells in vitro. As shown in our flow cytometry results (Fig. 3a), silencing of CBFA2T2 dramatically decreased the ALDHhigh population of 786-O cells and A-498 cells. Next, we evaluated the effect CBFA2T2 on sphere-forming ability of sorted ALDHhigh 786-O and A-498 cells in low-attachment plates. Depletion of CBFA2T2 decreased the number of spheres in sorted ALDHhigh 786-O cells relative to the control group (21.92 ± 6.02 vs 5.17 ± 3.59, p < 0.01), while depletion of CBFA2T2 decreased the number of spheres in sorted ALDHhigh A-498 cells relative to the control group (72.67 ± 20.95 vs 26.92 ± 11.24, p < 0.01) (Fig. 3b, c). This further supports that CBFA2T2 is required for maintenance of cancer stem cell property in RCC.
CBFA2T2 is required for ALDH1A3, OCT-4 and NANOG expression in CSCs of RCC
Previous reports have indicated that ALDHhigh cells showed elevated expression of OCT-4 and NANOG. To investigate whether CBFA2T2 was involved in stem cell marker genes expression of RCC CSCs, we first isolated ALDHhigh 786-O and A-498 cells. We maintained these cells as sphere in stem cells medium and treated with siRNA. We then examined the levels of these markers using quantitative RT-PCR approach. The qRT-PCR results confirmed that depletion of CBFA2T2 downregulated ALDH1A3, OCT-4 and NANOG expression in CSCs of 786-O and A-498 cells (Fig. 4a, b). Although, we didn’t detect any expression difference between two known CSC markers: CD133 and SOX2 (Fig. 4a, b). These findings suggest that CBFA2T2 maintained CSCs property by regulation of ALDH1A3, OCT4 and NANOG in RCC cell lines. Since CBFA2T2 normally functions as a transcription regulator, we reasoned that it may act as cofactor of stem cell transcriptional factor. To investigate the mechanism on how CBFA2T2 regulates CSCs in RCC, we examined the interaction between CBFA2T2 and OCT4. We found that CBFA2T2 was able to be pulled down by immunoprecipitation with OCT4 (Fig. 4c), suggesting it might serve as cofactor of OCT4 in regulating the stemness of CSCs population in RCC.
Cancer stem cells research is expanding rapidly based on the knowledge of CSCs from various cancer types. Understanding of how CSCs are regulated holds great implication in therapeutic resistance and relapse after initial therapy. In current study, we demonstrated that CBFA2T2 is an indication of the survival of RCC patients. In addition, our functional data supported a role of CBFA2T2 in ALDHhigh CSCs population isolated from RCC cell lines. Also, we showed that CBFA2T2 is required for expression of stem cells factors in ALDHhigh CSCs population.
In this report, we demonstrate that CBFA2T2 is required for CSCs maintenance in RCC cell lines. CBFA2T2 has been shown as an important factor for efficient inflammatory-associated colonic carcinogenesis . In acute myeloid leukemia, CBFA2T2 has been identified as a fusion partner of RUNX1 . Interestingly, fusion RUNX1/CBFA2T2 transcript is highly upregulated compared with wild-type CBFA2T2 gene , suggesting that CBFA2T2 alone might not be sufficient for the development of acute myeloid leukemia.
In a metastatic RCC cell line, the expression of the OCT-4 and NANOG gene was significantly overexpressed compared to all other RCC cell lines . Our analysis also showed that higher levels of OCT-4 and NANOG expression were associated with worse survival of RCC patients, indicating these stem cell pluripotency factors might be a common pathway for inhibition of CSCs.
Our data that CBFA2T2 could be pull down by immunoprecipitation with OCT4, suggesting it might serve as cofactor of OCT4 in regulating the stemness of CSCs population in RCC. Since OCT4 can regulate the expression of NANOG , it is highly possible that the interaction of OCT4 and CBFA2T2 is required for NANOG expression.
In summary, our results established a link between expression of CBFA2T2 and clinical prognosis in RCC. Our data also showed that the CBFA2T2 plays a critical role in the maintenance of CSC properties in human RCC cells.
We showed in our study that CBFA2T2 is required for maintenance of “stemness” through regulation of stem cells factors and CBFA2T2 could be a potential therapeutic target for RCC treatment.
Materials and methods
Cell culture and siRNA silencing
786-O and A-498 cell lines (ATCC, USA) were cultured in DMEM supplemented with 10% fetal bovine serum, penicillin (5000 U/ml) and streptomycin (5000 μg/ml) (Thermo Fisher Scientific, China) at 37 °C with 5% CO2. For siRNA assay, pool siRNAs consist of three to five target-specific 19- to 25-nucleotide siRNAs designed to knock down CBFA2T2 gene expression or scrambled siRNA were purchased from Santa Cruz Biotechnology. Monolayer or sphere 786-O and A-498 cells were transfected into cell line using Lipofectamine™ 2000 (Invitrogen) according to the manufacturer’s instructions.
Wound healing and invasion assays
For wound healing assay, 786-O and A498 cells transfected with the siRNAs against CBFA2T2 or control siRNA were plated on 6-well plates. When cells become confluent, cells were scraped with a plastic pipette tip . The wound closure was photographed at different time points (0 and 20 h after scraping). For invasion assay, 786-O and A-498 cells were transfected with the siRNAs for 48 h, and then digested and resuspended at 105 cells/ml. 200 μl of cell suspension was transferred into insert of the Transwell covered with Matrigel (BD Bioscience). After 48 h, the cells that migrated through the permeable membrane were fixed with paraformaldehyde, stained with haematoxylin and eosin, and finally counted. Aphidicolin (Sigma) was used to inhibit proliferation of cells during wound healing and invasion assay. Experiments were performed in triplicate.
Aldefluor assay and FACS sorting
The ALDH activity of 786-O and A-498 cells was determined by using the Aldefluor assay kit (Stem Cell Technologies) following the manufacturer’s guidance. Briefly, cells were dissociated into single cells by trypsin/EDTA digestion, washed twice in PBS and suspended in 1 ml Aldefluor assay buffer containing 5 μl ALDH substrate and incubated for 30 min at 37 °C in the dark. Diethylaminobenzaldehyde (DEAB) was added to a separate aliquot of the sample as negative control. Cells were then washed twice with buffer and maintained in ALDH buffer on ice. For FACS sorting, cells were sorted using an Aria II cell sorter (BD Biosciences). The sorting gates were established, by negative control cells which were treated with the ALDH inhibitor DEAB. Cells were sorted into DMEM medium containing 20% FBS, penicillin (5000 U/ml) and streptomycin (5000 μg/ml).
After FACS sorting, cells were then inoculated into 1:1 DMEM/F12 medium containing N2 supplement (100x, invitrogen), B27 supplement (50x, invitrogen), human recombinant epidermal growth factor (EGF; 10 ng/ml) and basic fibroblast growth factor (bFGF; 10 ng/ml) without serum at a density of 10,000 cells/well in ultra-low-attachment six-well plates (Corning, Inc., Corning, NY, USA). Fresh aliquots of stem cell medium were added every other day.
Real-time PCR and Western blot analysis
All sphere samples after siRNA treatment were centrifuged for 10 min at 300×g and the supernatant was removed. Total RNA was extracted with Trizol (Life Technologies, China) and reverse-transcribed with a first-strand cDNA synthesis kit (Roche Diagnostics). TaqMan PCR was performed as described [27, 28]. Stem cell marker mRNA levels were standardized by GAPDH levels and expressed as relative ratio to those of sphere cells treated with control siRNA. Primer sequences are:
ALDH1A3 forward: 5′-TGAATGGCACGAATCCAAGAG-3′
ALDH1A3 reverse: 5′-CACGTCGGGCTTATCTCCT-3′
PROM1 forward: 5′-AGTCGGAAACTGGCAGATAGC-3′
PROM1 reverse: 5′-GGTAGTGTTGTACTGGGCCAAT-3′
SOX2 forward: 5′-GCCGAGTGGAAACTTTTGTCG-3′
SOX2 reverse: 5′-GGCAGCGTGTACTTATCCTTCT-3′
OCT4 forward: 5′-CTGGGTTGATCCTCGGACCT-3′
OCT4 reverse: 5′-CCATCGGAGTTGCTCTCCA-3′
NANOG forward: 5′-TTTGTGGGCCTGAAGAAAACT-3′
NANOG reverse: 5′-AGGGCTGTCCTGAATAAGCAG-3′
GAPDH forward: 5′-GGAGCGAGATCCCTCCAAAAT-3′
GAPDH reverse: 5′-GGCTGTTGTCATACTTCTCATGG-3′
For Western blot assay, 786-O and A498 cells transfected with the siRNAs against CBFA2T2 or control siRNA for 72 h were washed two times with ice cold phosphate-buffered saline (PBS) and lysed in RIPA buffer (50 mM Tris pH 7.4, 250 mM NaCl, 5 mM EDTA, 1% NP-40, 0.1% SDS, 0.5% sodium deoxycholate, 1 mM phenylmethylsulphonyl fluoride) containing 1% protease inhibitor cocktail (Roche) . Cell lysates were centrifuged at 12,000×g for 10 min at 4 °C. Supernatant were collected for protein concentration measure using the BCA protein assay kit (Pierce). Total protein of 15 μg was subjected to SDS-PAGE, transferred to polyvinylidene fluoride (PVDF) membrane, and incubated with antibodies, followed by HRP-conjugated secondary antibodies. Specific proteins were detected by ECL Western blotting Detection Reagents (GE Healthcare Biosciences). Antibody against CBFA2T2 was purchased from Abcam (ab128164); antibody against α-tubulin was the products of Sigma-Aldrich (clone B-5-1-2).
Kaplan–Meier survival curves analysis
In this study, OncoLnc (http://www.oncolnc.org) was used as a tool for interactively exploring survival correlations . OncoLnc dataset contains survival data for 522 patients from kidney renal clear cell carcinoma (KIRC) cancer studies performed by The Cancer Genome Atlas (TCGA). The multivariate cox regressions were performed followed by a Kaplan–Meier analysis for CBFA2T2, OCT-4, ALDH1A3 and NANOG.
For statistical analysis, GraphPad Prism (version 7) was used. Student’s t-test was used to analyze statistical significance of the data. For Kaplan–Meier Survival, p-value represents the results of log-rank test. Chi square test was used for analyzing the correlation between clinicopathologic categories and CBFA2T2 expression. A p-value of less than 0.05 was considered to be statistically significant.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2016. CA Cancer J Clin. 2016;66:7–30.
Motzer RJ, Agarwal N, Beard C, et al. Kidney cancer. J Natl Compr Canc Netw. 2011;9:960–77.
Garcia JA, Rini BI. Recent progress in the management of advanced renal cell carcinoma. CA Cancer J Clin. 2007;57:112–25.
Richey SL, Culp SH, Jonasch E, et al. Outcome of patients with metastatic renal cell carcinoma treated with targeted therapy without cytoreductive nephrectomy. Ann Oncol. 2011;22:1048–53.
Fiori ME, Villanova L, De Maria R. Cancer stem cells: at the forefront of personalized medicine and immunotherapy. Curr Opin Pharmacol. 2017;35:1–11.
Yuan ZX, Mo J, Zhao G, Shu G, Fu HL, Zhao W. Targeting strategies for renal cell carcinoma: from renal cancer cells to renal cancer stem cells. Front Pharmacol. 2016;7:423.
Chen D, Wu M, Li Y, et al. Targeting BMI1+ cancer stem cells overcomes chemoresistance and inhibits metastases in squamous cell carcinoma. Cell stem cell. 2017;20:621–34.
Li J, Yu B, Deng P, et al. KDM3 epigenetically controls tumorigenic potentials of human colorectal cancer stem cells through Wnt/β-catenin signalling. Nat Commun. 2017;8:15146.
Bruno S, Bussolati B, Grange C, et al. CD133 + renal progenitor cells contribute to tumor angiogenesis. Am J Pathol. 2006;169:2223–35.
Gassenmaier M, Chen D, Buchner A, et al. CXC chemokine receptor 4 is essential for maintenance of renal cell carcinoma-initiating cells and predicts metastasis. Stem Cells. 2013;31:1467–76.
Debeb BG, Zhang X, Krishnamurthy S, et al. Characterizing cancer cells with cancer stem cell-like features in 293T human embryonic kidney cells. Mol Cancer. 2010;9:180.
Bussolati B, Bruno S, Grange C, Ferrando U, Camussi G. Identification of a tumor-initiating stem cell population in human renal carcinomas. FASEB J. 2008;22:3696–705.
Nicole L, Sanavia T, Veronese N, et al. Oncofetal gene SALL4 and prognosis in cancer: a systematic review with meta-analysis. Oncotarget. 2017;8:22968–79.
Ji Q, Luo YQ, Wang WH, Liu X, Li Q, Su SB. Research advances in traditional Chinese medicine syndromes in cancer patients. J Integr Med. 2016;14:12–21.
Zhong Y, Guan K, Guo S, et al. Spheres derived from the human SK-RC-42 renal cell carcinoma cell line are enriched in cancer stem cells. Cancer Lett. 2010;299:150–60.
Addla SK, Brown MD, Hart CA, Ramani VA, Clarke NW. Characterization of the Hoechst 33342 side population from normal and malignant human renal epithelial cells. Am J Physiol Renal Physiol. 2008;295:F680–7.
Wang LPP, Zhang H, La Marca F, Lin CY. Characterization of renal tumor-initiating cells in human renal carcinoma cell lines. Can Res. 2011;71(supplement 8):293.
Ueda K, Ogasawara S, Akiba J, et al. Aldehyde dehydrogenase 1 identifies cells with cancer stem cell-like properties in a human renal cell carcinoma cell line. PLoS ONE. 2013;8:e75463.
Cheng B, Yang G, Jiang R, et al. Cancer stem cell markers predict a poor prognosis in renal cell carcinoma: a meta-analysis. Oncotarget. 2016;7:65862–75.
Tu S, Narendra V, Yamaji M, et al. Co-repressor CBFA2T2 regulates pluripotency and germline development. Nature. 2016;534:387–90.
Guastadisegni MC, Lonoce A, Impera L, et al. CBFA2T2 and C20orf112: two novel fusion partners of RUNX1 in acute myeloid leukemia. Leukemia. 2010;24:1516–9.
Barrett CW, Fingleton B, Williams A, et al. MTGR1 is required for tumorigenesis in the murine AOM/DSS colitis-associated carcinoma model. Cancer Res. 2011;71:1302–12.
Anaya J. OncoLnc: linking TCGA survival data to mRNAs, miRNAs, and lncRNAs. PeerJ Comp Sci. 2016;2:e67.
Khan MI, Czarnecka AM, Lewicki S, et al. Comparative gene expression profiling of primary and metastatic renal cell carcinoma stem cell-like cancer cells. PLoS ONE. 2016;11:e0165718.
Rodda DJ, Chew JL, Lim LH, et al. Transcriptional regulation of nanog by OCT4 and SOX2. J Biol Chem. 2005;280:24731–7.
Li J, Guo Y, Feng X, et al. Receptor for activated C kinase 1 (RACK1): a regulator for migration and invasion in oral squamous cell carcinoma cells. J Cancer Res Clin Oncol. 2012;138:563–71.
Zhang YH, Wang Y, Yusufali AH, et al. Cytotoxic genes from traditional Chinese medicine inhibit tumor growth both in vitro and in vivo. J Integr Med. 2014;12:483–94.
Li Y, Deng H, Lv L, et al. The miR-193a-3p-regulated ING5 gene activates the DNA damage response pathway and inhibits multi-chemoresistance in bladder cancer. Oncotarget. 2015;6:10195–206.
Liu W, Zhou L, Zhou C, et al. GDF11 decreases bone mass by stimulating osteoclastogenesis and inhibiting osteoblast differentiation. Nat Commun. 2016;7:12794.
DCC, DC, YDL, YZW, CZA, XXZ, ML designed and performed the experiments. DCC, DC, YWY, YS, BY, JH, YG, JL, YZJ analyzed the datas and wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets generated and/or analysed during the current study are available from the corresponding author on reasonable request.
Consent for publication
Written informed consent was obtained from all patients.
Ethics approval and consent to participate
The study was approved by the institutional research ethics board.
This work was supported by National Natural Science Foundation of China (NSFC, 31501096 to M.L.; 81361120386, 31570751, 31270809 and 30930046 to R.C; 81500354 to Y.Z.J.); Shenzhen Science Foundation (JCYJ20160308104109234 to Y.Z.J); China Postdoctoral Science Foundation Grant (2016M602526 to Y.W.Y; 2016M600665 to X.X.Z.); Fundamental Research Funds for the Central Universities (20720150053 to M.L.); the National Basic Research Program of China (973 Programs 2013CB917802 to R.C.); the NSFC for Fostering Talents in Basic Research (J1310027 to J.L., Y.G. and X.C.); and XMU Training Program of Innovation and Entrepreneurship for Undergraduates (2015X0189 to X.W.; 2016Y0646 to Y.G., 103842017155 to J.L.).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.