Cdc7 overexpression is an independent prognostic marker and a potential therapeutic target in colorectal cancer
- Nathaniel Melling†2,
- Johanna Muth†1,
- Ronald Simon1,
- Carsten Bokemeyer3,
- Luigi Terracciano4,
- Guido Sauter1,
- Jakob Robert Izbicki2 and
- Andreas Holger Marx1Email author
© Melling et al. 2015
Received: 28 March 2015
Accepted: 9 July 2015
Published: 25 July 2015
Cdc7 is a widely expressed protein kinase implicated in cell division, cell cycle checkpoint mechanisms and cancer progression. Recently, it has been suggested as a target for anti-cancer therapy.
To determine the relationship of Cdc7 protein expression with tumor phenotype, molecular features and prognosis, 1800 colorectal carcinomas were analyzed by immunohistochemistry on a tissue microarray.
Cdc7 expression was considered negative in 33.6 %, weak in 57.2 % and strong in 9.2 % of 1711 interpretable CRCs. Loss of Cdc7 expression was significantly associated with high tumor stage (p < 0.0001) and high tumor grade (p = 0.0077), but was unrelated to the nodal status (p = 0.5957). Moreover, a link between Cdc7 expression and the tubular histological tumor type was seen (p < 0.0001). p53 and Cdc7 expression were significantly linked to each other (p = 0.0013). In a multivariate survival analysis, strong Cdc7 expression of CRC was an independent marker of improved patient survival (p = 0.0031).
Our data show that Cdc7 is highly expressed in CRC and a potential therapeutic target in a subset of cancers with high p53 expression. Moreover, our findings strongly argue for a clinical utility of Cdc7 immunostaining as an independent prognostic biomarker in colorectal cancer enabling to select patients for adjuvant treatment.
KeywordsBiological markers Cell cycle checkpoints Colorectal neoplasms Immunohistochemistry Prognosis
Colorectal cancer (CRC) is the fourth most common malignant disease with over 1 million novel cases and over 500.000 deaths each year worldwide .
Although recent advances in the management of the disease have improved outcomes, CRC remains the second leading cause of cancer-related death in Western countries . In advanced metastatic colorectal cancer (mCRC), surgery alone is not curative and therefore adjuvant chemotherapy is needed. There is much promise in targeted anti-cancer therapies. Encouraged by the success of targeted therapies directed against key molecules of cell growth like the HER2, EGFR, or KIT tyrosine kinase receptors, researchers have extended their survey for promising candidate targets to genes that participate in the control of cell division and replication. Such genes represent the most downstream effectors of growth signaling pathways and therefore it is believed that their inhibition will be effective in a broad range of different tumor types characterized by rapid cell proliferation.
The cell division cycle 7- (Cdc7) related protein kinase is an enzyme that is essential for DNA replication in human cells [2, 3]. Overexpression of Cdc7 and its protein regulator Dbf4 has been reported in many human tumors , including ovarian cancer , melanoma [6, 7], diffuse large B-cell lymphoma , oral squamous cell carcinoma  and breast cancer . Cdc7 and Dbf4 overexpression are reported to cause cell-cycle arrest in S phase and it has been hypothesized that increased Cdc7 activity may aid recovery or repair of stalled replication forks to enhance survival of tumor cells . Therefore, it can be assumed that alterations in Cdc7/Dbf4 protein activity during tumorigenesis may have important consequences for tumor cell survival, underlining the potential of Cdc7 as an anticancer target. High expression of Cdc7 protein correlates with poor prognosis in patients with diffuse large B-cell lymphoma and is a marker of resistance to DNA-damaging agents in oral squamous cell carcinoma [9, 8]. Preclinical data show that Cdc7 is a novel and promising target for tumor-cell killing, as has been shown with different inhibitors [12–15]. Little is known about the prevalence and significance of Cdc7 in colorectal cancer. To study the potential role of Cdc7 in CRC a tissue microarray containing 1.800 tumor samples with clinical follow-up data was analyzed. Our data suggest that Cdc7 overexpression may point out a small but significant subset of immunohistochemically p53-positive CRC that could benefit from anti-Cdc7 treatment.
Patients and tissue microarray (TMA) construction
Clinical and pathological features of colorectal cancers
Mean: 69 (29 – 96)
Total number of patients
The utilization of tissues and clinical data was according to the Hamburger Krankenhaus Gesetz (§12 HmbKHG) and approved by our local Ethical Committee.
Standard indirect immunoperoxidase procedures were used for the detection of Cdc7 (abcam, clone SPM171, dilution 1:150). Sections were heated in an autoclave at 121 °C for 10 min in citrate buffer (pH 9.0). Diaminobenzidine was used as a chromogen and sections were counterstained with Mayer’s hematoxylin. Immunostaining was typically nuclear. For tumor tissue the percentage of positive cells was estimated and the staining intensity was recorded as 1+, 2+ or 3+. For statistical analyses, the staining results were categorized into three groups as previously described . Tumors without any staining were considered “negative”. Tumors showing at least weak Cdc7 staining were considered “positive”. Tumors with 1+ or 2+ positivity in up to 50 % or 3+ positivity in up to 20 % of cells were considered “weakly positive”. Tumors with 2+ staining in >50 % or 3+ staining in >20 % of cells were considered “strongly positive”. The molecular database attached to this TMA contained results on p53 expression in 1800 cancers.
Statistical calculations were performed with JMP® 10.0.2 software (2012 SAS Institute Inc., NC, USA). Contingency tables and the chi2-test were performed to search for associations between molecular parameters and tumor phenotype. Survival curves were calculated according to Kaplan-Meier. The Log-Rank test was applied to detect significant survival differences between groups. Cox proportional hazards regression analysis was performed to test the statistical independence and significance between pathological and clinical variables.
A total of 1711 (95.1 %) of tumor samples were interpretable in our TMA analysis (Table 1). Reasons for non-informative cases (289 spots; 4.9 %) included lack of tissue samples or absence of unequivocal cancer tissue in the TMA spot.
Cdc7 expression in colorectal cancer
Cdc7 IHC, clinico-pathological and molecular parameters in CRC
Association with tumor localization
Cdc7 expression levels were not related to tumor localization (Table 2, p = 0.2093).
Association with p53 expression
Multivariate analysis of Cdc7 in colorectal cancer
The results of this study show that in CRC strong Cdc7 expression is an independent favorable prognostic marker and significantly linked to immunohistochemical p53 positivity. Cdc7 expression was found in 66.4 % of all interpretable CRCs, and it was typically described as weak in normal tissue. Only few reports on results of Cdc7 expression in CRC exist. Chen et al. reported a significantly higher mRNA and protein expression of Cdc7 in 39 colorectal cancers as compared to tumor-adjacent normal colorectal tissue  and colorectal cancer is among the cancer types that shows the strongest up regulation in Cdc7 expression as determined from data obtained by mining Oncomine and Gene Expression Omnibus databases. Bonte et al. examined Cdc7 expression in an array of human malignancies, including colon cancer where they found low Cdc7 expression in 2 and high Cdc7 expression in 8 of 10 cases .
The percentage of CRC cases with "strong" Cdc7 staining in our study is 9.2 %. This may be attributable to a mismatch between mRNA and protein levels, but it is more likely caused by immunohistochemistry conditions. It is well known that varying antibody dilutions lead to significant changes in the rate of positive cases . This is all the more expected in the case of ubiquitously expressed proteins, such as Cdc7.
Various functional studies had previously demonstrated that Cdc7 is essential for cell proliferation [20, 21]. The strong association between Cdc7 and p53 overexpression found in our study was therefore expected, although p-values may appear stronger than the absolute values imply, due to the high caseload (e.g. p53 positive vs p53 negative 11.7 % vs 9.8 %; p = 0.0013). A previous study has shown that a p53-dependent pathway in normal fibroblasts actively prevents progression through a lethal S phase in the absence of sufficient Cdc7 kinase. In contrast, Cdc7 and p53-depleted tumor cells do not elicit a robust checkpoint response. Therefore, it has been suggested that p53 is required for cell viability in Cdc7 deficient cells . Moreover, it has been reported that ATR-dependent activation of p38 MAP kinase is responsible for apoptotic cell death in Cdc7-depleted cells . In CRC, up to 42 % of the cases show p53-mutations with a strong inverse correlation between p53 alterations and MSI . Thus, we hypothesize that Cdc7 inhibition may serve as a new therapeutical approach in a subgroup of p53-deficient CRC, which are predominately localized in the left colon .
In our study, Cdc7 overexpression proved to be an independent marker for good prognosis in CRC. This finding is in contrast to studies on other cancer types where Cdc7 overexpression is linked with worse prognosis. For example, Cheng et al. reported in their study on oral squamous cell carcinoma (OSCC) that Cdc7 overexpression is an unfavorable prognostic marker and suggested Cdc7 overexpression contributing to the resistance to DNA-damaging agents .
In another recent study of our group on more than 2.100 breast cancers, Cdc7 expression was seen in 57 % of the cases. Comparable to what was found in CRC, in breast cancer high levels of Cdc7 expression were also significantly related to high tumor grade, high Ki67 expression and p53 overexpression. Although the IHC scoring system was slightly different in this breast cancer study (4 categories vs 3 categories as described earlier) , the percentages of Cdc7 negative, weakly and strongly positive breast cancers were similar to our findings in CRC .
However, when comparing Cdc7 expression and tumor grade and stage, results in both studies opposed each other. Whereas in CRC low levels of Cdc7 expression were significantly associated with high tumor grade and stage, in breast cancer high levels of Cdc7 expression were linked with higher tumor grade and stage . The reason for this discrepancy is unknown but description of variable biological roles of one biomarker in different tumor entities is not uncommon. For example, AQP5 expression is associated with poor prognosis in breast cancer [25–27], hepatocellular carcinoma , colorectal carcinoma [29, 30], chronic myelogenous leukemia , lung cancer [32, 33], ovarian cancer [34, 35], cervical cancer  and oral squamous cell carcinoma , but to better outcome in gallbladder  and biliary tract cancer .
Furthermore, some biomarkers, such as GSK3ß, are known to play both tumor suppressive and promoting roles depending on the microenvironment It has, thus, been suggested that GSK3ß function may switch between tumor promotion and suppression depending on cell properties and the stage of tumorigenesis .
Therefore, varying oncogenic properties of Cdc7 might also be attributable to differing molecular backgrounds in various cancer types.
Of note, Cdc7 has been shown to represent a potent and highly specific anticancer target in p53-mutant, Her2-overexpressing and triple-negative breast cancers . Based on these findings, a number of different Cdc7 kinase inhibitors which showed anti-tumor activity have been developed recently [12–14]. Our data show that Cdc7 is expressed in the majority of colorectal cancers and therefore could serve as a potential therapy target.
It is well known that immunostaining of p53 can serve as a surrogate marker for p53 mutations in cancer . Thus, CRC, which show p53 overexpression are likely to harbor p53 mutations and therefore may be susceptible to anti-Cdc7 therapy .
However, experimental and clinical studies examining the efficacy of Cdc7 kinase inhibitors in CRC are still needed as well as data on Dbf4 protein expression, since colorectal cancers have been shown to have strong levels of Dbf4 mRNA expression and Dbf4 expression levels were uniformly high among most colorectal tumors (Bonte et al.). This data could well be recruited by a TMA approach similar to ours.
In summary, our data show that Cdc7 is highly expressed in CRC and may represent a potential therapeutic target in a subset of cancers with high p53 expression.
Finally, our findings strongly argue for a clinical utility of Cdc7 immunostaining as an independent prognostic biomarker in colorectal cancer enabling to select patients for adjuvant treatment.
We thank Martina Mirlacher, Christina Koop and Inge Brandt for excellent technical assistance.
- Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global cancer statistics. CA Cancer J Clin. 2011;61(2):69–90. doi:https://doi.org/10.3322/caac.20107.PubMedView ArticleGoogle Scholar
- Jiang W, Hunter T. Identification and characterization of a human protein kinase related to budding yeast Cdc7p. Proc Natl Acad Sci U S A. 1997;94(26):14320–5.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim JM, Yamada M, Masai H. Functions of mammalian Cdc7 kinase in initiation/monitoring of DNA replication and development. Mutat Res. 2003;532(1–2):29–40.PubMedView ArticleGoogle Scholar
- Bonte D, Lindvall C, Liu H, Dykema K, Furge K, Weinreich M. Cdc7-Dbf4 kinase overexpression in multiple cancers and tumor cell lines is correlated with p53 inactivation. Neoplasia. 2008;10(9):920–31.PubMed CentralPubMedView ArticleGoogle Scholar
- Kulkarni AA, Kingsbury SR, Tudzarova S, Hong HK, Loddo M, Rashid M, et al. Cdc7 Kinase Is a Predictor of Survival and a Novel Therapeutic Target in Epithelial Ovarian Carcinoma. Clin Cancer Res. 2009;15:2417–25.PubMedView ArticleGoogle Scholar
- Nambiar S, Mirmohammadsadegh A, Hassan M, Mota R, Marini A, Alaoui A, et al. Identification and functional characterization of ASK/Dbf4, a novel cell survival gene in cutaneous melanoma with prognostic relevance. Carcinogenesis. 2007;28(12):2501–10. doi:https://doi.org/10.1093/carcin/bgm197.PubMedView ArticleGoogle Scholar
- Clarke LE, Fountaine TJ, Hennessy J, Bruggeman RD, Clarke JT, Mauger DT, et al. Cdc7 expression in melanomas, Spitz tumors and melanocytic nevi. J Cutan Pathol. 2009;36(4):433–8.PubMedView ArticleGoogle Scholar
- Hou Y, Wang HQ, Ba Y. High expression of cell division cycle 7 protein correlates with poor prognosis in patients with diffuse large B-cell lymphoma. Med Oncol. 2012;29(5):3498–503. doi:https://doi.org/10.1007/s12032-012-0223-y.PubMedView ArticleGoogle Scholar
- Cheng AN, Jiang SS, Fan CC, Lo YK, Kuo CY, Chen CH, et al. Increased Cdc7 expression is a marker of oral squamous cell carcinoma and overexpression of Cdc7 contributes to the resistance to DNA-damaging agents. Cancer Lett. 2013;337(2):218–25. doi:https://doi.org/10.1016/j.canlet.2013.05.008.PubMedView ArticleGoogle Scholar
- Choschzick M, Lebeau A, Marx AH, Tharun L, Terracciano L, Heilenkotter U, et al. Overexpression of cell division cycle 7 homolog is associated with gene amplification frequency in breast cancer. Hum Pathol. 2010;41(3):358–65. doi:https://doi.org/10.1016/j.humpath.2009.08.008.PubMedView ArticleGoogle Scholar
- Guo B, Romero J, Kim BJ, Lee H. High levels of Cdc7 and Dbf4 proteins can arrest cell-cycle progression. Eur J Cell Biol. 2005;84(12):927–38. doi:https://doi.org/10.1016/j.ejcb.2005.09.016.PubMedView ArticleGoogle Scholar
- Koltun ES, Tsuhako AL, Brown DS, Aay N, Arcalas A, Chan V, et al. Discovery of XL413, a potent and selective CDC7 inhibitor. Bioorg Med Chem Lett. 2012;22(11):3727–31. doi:https://doi.org/10.1016/j.bmcl.2012.04.024.PubMedView ArticleGoogle Scholar
- Robertson S. 339 POSTER Identification of XL413, a selective Cdc7 kinase inhibitor which induces cell cycle arrest and exhibits potent antitumor activity. Eur J Cancer Suppl. 2008;6(12):108.View ArticleGoogle Scholar
- Frattini M, Shum D, O'Dwyer K, Brentjens R, Yeh R, Maslak P et al., editors. Discovery and validation of a novel class of small molecule inhibitors of the CDC7 kinase: modulation of tumor growth in vitro and in vivo. Blood (ASH Annual Meeting Abstracts); 2008.Google Scholar
- Natoni A, Murillo LS, Kliszczak AE, Catherwood MA, Montagnoli A, Samali A, et al. Mechanisms of action of a dual Cdc7/Cdk9 kinase inhibitor against quiescent and proliferating CLL cells. Mol Cancer Ther. 2011;10(9):1624–34. doi:https://doi.org/10.1158/1535-7163.MCT-10-1119.PubMedView ArticleGoogle Scholar
- Kononen J, Bubendorf L, Kallioniemi A, Barlund M, Schraml P, Leighton S, et al. Tissue microarrays for high-throughput molecular profiling of tumor specimens. Nat Med. 1998;4(7):844–7.PubMedView ArticleGoogle Scholar
- Marx A, Simon P, Simon R, Mirlacher M, Izbicki JR, Yekebas E, et al. AMACR expression in colorectal cancer is associated with left-sided tumor localization. Virchows Arch. 2008;453(3):243–8. doi:https://doi.org/10.1007/s00428-008-0646-1.PubMedView ArticleGoogle Scholar
- Chen HJ, Zhu Z, Wang XL, Feng QL, Wu Q, Xu ZP, et al. Expression of huCdc7 in colorectal cancer. World J Gastroenterol. 2013;19(20):3130–3. doi:https://doi.org/10.3748/wjg.v19.i20.3130.PubMed CentralPubMedView ArticleGoogle Scholar
- Schlomm T, Iwers L, Kirstein P, Jessen B, Kollermann J, Minner S, et al. Clinical significance of p53 alterations in surgically treated prostate cancers. Mod Pathol. 2008;21(11):1371–8. doi:https://doi.org/10.1038/modpathol.2008.104.PubMedView ArticleGoogle Scholar
- Masai H, Arai K. Cdc7 kinase complex: a key regulator in the initiation of DNA replication. J Cell Physiol. 2002;190(3):287–96.PubMedView ArticleGoogle Scholar
- Sclafani RA. Cdc7p-Dbf4p becomes famous in the cell cycle. J Cell Sci. 2000;113(Pt 12):2111–7.PubMedGoogle Scholar
- Montagnoli A, Tenca P, Sola F, Carpani D, Brotherton D, Albanese C, et al. Cdc7 inhibition reveals a p53-dependent replication checkpoint that is defective in cancer cells. Cancer Res. 2004;64(19):7110–6.PubMedView ArticleGoogle Scholar
- Im JS, Lee JK. ATR-dependent activation of p38 MAP kinase is responsible for apoptotic cell death in cells depleted of Cdc7. J Biol Chem. 2008;283(37):25171–7. doi:https://doi.org/10.1074/jbc.M802851200.PubMedView ArticleGoogle Scholar
- Elsaleh H, Powell B, McCaul K, Grieu F, Grant R, Joseph D, et al. P53 alteration and microsatellite instability have predictive value for survival benefit from chemotherapy in stage III colorectal carcinoma. Clin Cancer Res. 2001;7(5):1343–9.PubMedGoogle Scholar
- Lee SJ, Chae YS, Kim JG, Kim WW, Jung JH, Park HY, et al. AQP5 expression predicts survival in patients with early breast cancer. Ann Surg Oncol. 2014;21(2):375–83. doi:https://doi.org/10.1245/s10434-013-3317-7.PubMedView ArticleGoogle Scholar
- Jung HJ, Park JY, Jeon HS, Kwon TH. Aquaporin-5: a marker protein for proliferation and migration of human breast cancer cells. PloS one. 2011;6(12):e28492. doi:https://doi.org/10.1371/journal.pone.0028492.PubMed CentralPubMedView ArticleGoogle Scholar
- Shi Z, Zhang T, Luo L, Zhao H, Cheng J, Xiang J, et al. Aquaporins in human breast cancer: identification and involvement in carcinogenesis of breast cancer. J Surg Oncol. 2012;106(3):267–72. doi:https://doi.org/10.1002/jso.22155.PubMedView ArticleGoogle Scholar
- Guo X, Sun T, Yang M, Li Z, Li Z, Gao Y. Prognostic value of combined aquaporin 3 and aquaporin 5 overexpression in hepatocellular carcinoma. Biomed Res Int. 2013;2013:206525. doi:https://doi.org/10.1155/2013/206525.PubMed CentralPubMedGoogle Scholar
- Kang SK, Chae YK, Woo J, Kim MS, Park JC, Lee J, et al. Role of human aquaporin 5 in colorectal carcinogenesis. Am J Pathol. 2008;173(2):518–25. doi:https://doi.org/10.2353/ajpath.2008.071198.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang W, Li Q, Yang T, Bai G, Li D, Li Q, et al. Expression of AQP5 and AQP8 in human colorectal carcinoma and their clinical significance. World J Surg Oncol. 2012;10:242. doi:https://doi.org/10.1186/1477-7819-10-242.PubMed CentralPubMedView ArticleGoogle Scholar
- Chae YK, Kang SK, Kim MS, Woo J, Lee J, Chang S, et al. Human AQP5 plays a role in the progression of chronic myelogenous leukemia (CML). PloS One. 2008;3(7):e2594. doi:https://doi.org/10.1371/journal.pone.0002594.PubMed CentralPubMedView ArticleGoogle Scholar
- Chae YK, Woo J, Kim MJ, Kang SK, Kim MS, Lee J, et al. Expression of aquaporin 5 (AQP5) promotes tumor invasion in human non small cell lung cancer. PloS One. 2008;3(5):e2162. doi:https://doi.org/10.1371/journal.pone.0002162.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhang Z, Chen Z, Song Y, Zhang P, Hu J, Bai C. Expression of aquaporin 5 increases proliferation and metastasis potential of lung cancer. J Pathol. 2010;221(2):210–20. doi:https://doi.org/10.1002/path.2702.PubMedView ArticleGoogle Scholar
- Yang JH, Shi YF, Cheng Q, Deng L. Expression and localization of aquaporin-5 in the epithelial ovarian tumors. Gynecol Oncol. 2006;100(2):294–9. doi:https://doi.org/10.1016/j.ygyno.2005.08.054.PubMedView ArticleGoogle Scholar
- Yan C, Yang J, Shen L, Chen X. Inhibitory effect of Epigallocatechin gallate on ovarian cancer cell proliferation associated with aquaporin 5 expression. Arch Gynecol Obstet. 2012;285(2):459–67. doi:https://doi.org/10.1007/s00404-011-1942-6.PubMedView ArticleGoogle Scholar
- Zhang T, Zhao C, Chen D, Zhou Z. Overexpression of AQP5 in cervical cancer: correlation with clinicopathological features and prognosis. Med Oncol. 2012;29(3):1998–2004. doi:https://doi.org/10.1007/s12032-011-0095-6.PubMedView ArticleGoogle Scholar
- Ishimoto S, Wada K, Usami Y, Tanaka N, Aikawa T, Okura M, et al. Differential expression of aquaporin 5 and aquaporin 3 in squamous cell carcinoma and adenoid cystic carcinoma. Int J Oncol. 2012;41(1):67–75. doi:https://doi.org/10.3892/ijo.2012.1445.PubMedGoogle Scholar
- Sekine S, Shimada Y, Nagata T, Sawada S, Yoshioka I, Matsui K, et al. Role of aquaporin-5 in gallbladder carcinoma. Eur Surg Res. 2013;51(3–4):108–17. doi:https://doi.org/10.1159/000355675.PubMedView ArticleGoogle Scholar
- Sekine S, Shimada Y, Nagata T, Moriyama M, Omura T, Watanabe T, et al. Prognostic significance of aquaporins in human biliary tract carcinoma. Oncol Rep. 2012;27(6):1741–7. doi:https://doi.org/10.3892/or.2012.1747.PubMedGoogle Scholar
- Zhang Q, Bhojani MS, Ben-Josef E, Spalding AC, Kuick R, Sun Y, et al. Glycogen Synthase Kinase 3beta in Pancreatic Cancer and its Implications in Chemotherapy and Radiation Therapy. J Carcinog Mutagen. 2013;4(3):147. doi:https://doi.org/10.4172/2157-2518.1000147.PubMed CentralPubMedGoogle Scholar
- Rodriguez-Acebes S, Proctor I, Loddo M, Wollenschlaeger A, Rashid M, Falzon M, et al. Targeting DNA replication before it starts: Cdc7 as a therapeutic target in p53-mutant breast cancers. Am J Pathol. 2010;177(4):2034–45. doi:https://doi.org/10.2353/ajpath.2010.100421.PubMed CentralPubMedView ArticleGoogle Scholar
- Yemelyanova A, Vang R, Kshirsagar M, Lu D, Marks MA, Shih Ie M, et al. Immunohistochemical staining patterns of p53 can serve as a surrogate marker for TP53 mutations in ovarian carcinoma: an immunohistochemical and nucleotide sequencing analysis. Mod Pathol. 2011;24(9):1248–53. doi:https://doi.org/10.1038/modpathol.2011.85.PubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.