KIFCI, a novel putative prognostic biomarker for ovarian adenocarcinomas: delineating protein interaction networks and signaling circuitries
- Shrikant Pawar†1,
- Shashikiran Donthamsetty†1,
- Vaishali Pannu1,
- Padmashree Rida1,
- Angela Ogden1,
- Nathan Bowen2,
- Remus Osan3,
- Guilherme Cantuaria4Email author and
- Ritu Aneja1Email author
© Pawar et al.; licensee BioMed Central Ltd. 2014
Received: 14 April 2014
Accepted: 24 April 2014
Published: 12 May 2014
Amplified centrosomes in cancers are recently garnering a lot of attention as an emerging hub of diagnostic, prognostic and therapeutic targets. Ovarian adenocarcinomas commonly harbor supernumerary centrosomes that drive chromosomal instability. A centrosome clustering molecule, KIFC1, is indispensable for the viability of extra centrosome-bearing cancer cells, and may underlie progression of ovarian cancers.
Centrosome amplification in low- and high- grade serous ovarian adenocarcinomas was quantitated employing confocal imaging. KIFC1 expression was analyzed in ovarian tumors using publically-available databases. Associated grade, stage and clinical information from these databases were plotted for KIFC1 gene expression values. Furthermore, interactions and functional annotation of KIFC1 and its highly correlated genes were studied using DAVID and STRING 9.1.
Clinical specimens of ovarian cancers display robust centrosome amplification and deploy centrosome clustering to execute an error-prone mitosis to enable karyotypic heterogeneity that fosters tumor progression and aggressiveness. Our in silico analyses showed KIFC1 overexpression in human ovarian tumors (n = 1090) and its upregulation associated with tumor aggressiveness utilizing publically-available gene expression databases. KIFC1 expression correlated with advanced tumor grade and stage. Dichotomization of KIFC1 levels revealed a significantly lower overall survival time for patients in high KIFC1 group. Intriguingly, in a matched-cohort of primary (n = 7) and metastatic (n = 7) ovarian samples, no significant differences in KIFC1 expression were detectable, suggesting that high KIFC1 expression may serve as a marker of metastases onset. Nonetheless, KIFC1 levels in both primary and matched metastatic sites were significantly higher compared to normal tissue . Ingenuity based network prediction algorithms combined with pre-established protein interaction networks uncovered several novel cell-cycle related partner genes on the basis of interconnectivity, illuminating the centrosome clustering independent agenda of KIFC1 in ovarian tumor progression.
Ovarian cancers display amplified centrosomes, a feature of aggressive tumors. To cope up with the abnormal centrosomal load, ovarian cancer cells upregulate genes like KIFC1 that are known to induce centrosome clustering. Our data underscore KIFC1 as a putative biomarker that predicts worse prognosis, poor overall survival and may serve as a potential marker of onset of metastatic dissemination in ovarian cancer patients.
Centrosome amplification (CA) is a hallmark of cancers [1–3]. Recent evidence suggests that amplified centrosomes can drive malignant transformation  and perhaps fuel metastatic dissemination . Logically, excess centrosomes would orchestrate a multipolar spindle which might result in inviable progeny and jeopardize the survival of extra centrosome bearing cancer cells. Cancer cells however, overcome this paradoxical situation in their favor through multiple mechanisms including centrosome clustering, a cancer cell-specific trait [5–7]. It is being revealedthat cancer cells have indeed evolved quite a sophisticated and extensive arsenal of ‘clever tactics’ to hijack cellular mechanisms and deploy them to cluster supernumerary centrosomes into two polar groups to allow formation of a pseudo-bipolar mitotic spindle [6–9]. A centrosome clustering molecule, KIFC1 (also known as HSET), a minus end-directed motor protein of the kinesin-14 family, is essential for the viability of extra centrosome-bearing cancer cells .
Recently, supernumerary centrosomes and high KIFC1 expression have been associated with chromosome missegregation that results in low-grade aneuploidy, a landmark of cancers [6–9]. Given that normal cells most often have two centrosomes, they do not rely on centrosome clustering mechanisms; as a result, targeting KIFC1 is an attractive chemotherapeutic strategy . Proof-of-concept comes from the recent discovery and preclinical development of two novel KIFC1 small molecule inhibitors AZ82  and CW069 , that cause centrosome declustering exclusively in cancer cells with amplified centrosomes . Given the excellent promise of KIFC1 as a therapeutic target, its role as a negative prognosticator merits investigation in cancers with amplified centrosomes. Since KIFC1 has been shown to predict non-small cell lung cancer metastasis to the brain , we were inquisitive to examine its usefulness as a risk predictor and/or negative prognosticator in other aggressive cancer types which harbor amplified centrosomes. Although centrosomal aberrations in epithelial ovarian cancers (EOC) have been a relatively understudied area, abnormalities in centrosomes have been reported in ovarian tumors . Centrosomal aberrations may be an early event in ovarian carcinogenesis and implicated in ovarian tumor progression. The features of many epithelial tumors, including EOC, are the presence of aneuploidy, a consequence of chromosomal instability (CIN) that arises due to aberrant CA . In several solid malignancies, amplified centrosomes are a potential indicator of cancer aggressiveness .
Herein we examined the severity and extent of centrosome amplification in low- and high- grade serous ovarian adenocarcinomas employing multicolor immunofluorescence confocal imaging. In these clinical specimens, we also visualized various cell cycle stages and spindle architecture of cells in mitosis to gain insights into the propensity of ovarian tumors to undergo aberrant cell divisions that foster CIN, karyotypic heterogeneity and generation of aneuploid clones. Further, we quantitated the extent of spindle polarity, in particular, cells with multipolar spindle configurations as well as cells with pseudo-bipolar spindles with centrosomes corralled at the two poles. Given the link between presence of excess centrosomes and upregulation of KIFC1, an important member of the centrosome clustering arsenal, we evaluated KIFC1’s potential as a negative prognostic indicator in EOC. Using independent gene expression datasets, we identified that KIFC1 gene in ovarian cancer is expressed at least 2-fold (logarithm to base 2 scale) higher than in normal ovaries. Our in silico data also suggest a correlation between KIFC1 and grade, stage, and clinical outcomes in EOC. The gene expression profiling-based identification of KIFC1 as a negative prognosticator in EOC may improve evaluation of disease course. In addition, an in silico-guided mechanistic understanding of KIFC1 gene interactions have delineated pathways and protein interactions which illuminate previously unrecognized partners of this centrosome clustering molecule to unravel the biological behavior of ovarian tumors.
Epithelial ovarian cancers harbor amplified centrosomes
Epithelial ovarian cancers exhibit mitotic aberrations and abnormal spindle architecture
Given the notable differences in centrosomal aberrations between low and high grade ovarian cancers, we next asked if differences in centrosomal aberrations between tumor grades translated into differences in mitotic and spindle aberrations. We found that low- and high- grade tumors significantly differed in the proportion of cells harboring aberrant mitotic spindles. Over 50% mitotic cells in low-grade ovarian tumors exhibited multipolar spindles in stark contrast to a mere 20% multipolar mitotic cells in high-grade tumors (Figure 1C, D). Intriguingly, multipolar spindles observed in both low- and high- grade tumors were predominantly tripolar, suggesting the enabling role of spindle multipolarity in promoting chromosome missegregation that underlies low-grade aneuploidy as opposed to mitotic catastrophe.
Epithelial ovarian cancers display clustering of amplified centrosomes
The discordance observed in the extent of CA and multipolar mitosis in high-grade tumors naturally led us to hypothesize that high-grade ovarian tumors manage their “excess baggage” of supernumerary centrosomes by clustering them into “pseudobipolar” spindles. Thus we evaluated the extent of centrosome clustering by counting the number of mitotic cells harboring “pseudo-bipolar” spindles. Intriguingly, our observations revealed that 80% of the mitotic cells in high-grade tumors display centrosome clustering compared to only 20% in low-grade tumors (Figure 1C, E). Taken together, our data suggest that ovarian cancers display a high degree of CA which is grade-dependent. Supernumerary centrosomes drive the assembly of multipolar spindle during mitosis to facilitate higher degree of chromosomal heterogeneity, mostly numerical in low-grade tumors. However, high-grade tumors tend to cluster their extra centrosomes and form a “pseudobipolar” spindle to enable low-level of chromosome missegregation as compared to low-grade tumors.
KIFC1 gene expression level is high in ovarian cancers compared to uninvolved normal ovarian tissue
KIFC1 gene expression increases with grade in ovarian cancers
Having found that KIFC1 expression was higher in EOC as compared to normal ovarian tissue, we next assessed if KIFC1 expression changes with grade and stage within EOC. Ovarian carcinoma is categorized by four stages, ranging from stage I in which the cancer is localized to the ovary/ovaries to stage IV in that the cancer has spread outside of the peritoneal cavity . Grades for ovarian carcinoma are classified based on their histological appearance. Grade 1 (well differentiated, grade 2 moderately differentiated, and grade 3 poorly differentiated) [18, 19]. Based on 468 EOC samples, average KIFC1 expression levels were ~6.53 for grade 1, ~7.85 for grade 2, and ~8.25 for grade 3, which were significantly different (p < 0.01) amongst these sub-grades (Figure 2C). As a result, higher KIFC1 is associated with increased grade. Average KIFC1 expression levels were approximately 7.62 for stage I, 8.00 for stage II, 8.10 for stage III, and 7.82 for stage IV, and these differences were not statistically significant amongst these sub-stages (Figure 2B). Although we found a significant correlation between KIFC1 expression levels within sub-grades, such differences were not observed within sub-stages. Given our finding that KIFC1 expression increases with grade, we were next interested in determining whether it is expressed at higher levels in metastatic versus primary ovarian carcinomas (Figure 2D). Although both primary and metastatic tumors showed significantly higher expression of KIFC1 as compared to normal tissue, there was no significant difference in the expression value of KIFC1 between matched primary and metastatic ovarian carcinomas (Figure 2D). These fourteen matched sets of primary and metastatic (omental) samples were collected from seven advanced staged (III/IV) ovarian cancer (serous adenocarcinoma) patients and were a part of a recently published paper from John McDonald’s group .
Increased KIFC1 expression is associated with poorer overall survival in age-specific ovarian cancer patients
Comparing survival in days for patients with ovarian cancer categorized in different age groups and with high and low KIFC1 expression levels
Age groups (Years)
Survival (days) for high KIFC1 expression
Survival (days) for low KIFC1 expression
831.7 (N = 11)
1262.2 (N = 12)
916.2 (N = 57)
1100.2 (N = 54)
1025.1 (N = 89)
1147.3 (N = 103)
1239.4 (N = 61)
975.9 (N = 64)
Interactions and functional annotation of KIFC1 and its highly correlated genes using DAVID and STRING 9.1
Pathways and genes involved with KIFC1 and its correlated genes
No of genes
% of genes involved
Cell cycle checkpoint
Mitotic cell cycle
Confidence values for KIFC1 and respective interactions
Over the recent years, the role of KIFC1 in centrosome clustering in cancer cells with supernumerary centrosomes has been well recognized. Our present study shows that both low- and high- grade ovarian tumors display CA. This finding is consistent with a previous report showing presence of CA in ovarian tumors in a stage-dependent manner . Interestingly, we found that while low-grade tumors display a higher proportion of mitotic cells displaying multipolar mitosis compared to high-grade ones, centrosome clustering was almost an exclusive feature of high-grade tumors compared to low-grade ones. As part of the tumor evolution agenda, we rationalize that when tumors are low-grade, multipolar mitosis reflects their proclivity to undergo aberrant mitoses and subsequent cell divisions. Perhaps multipolar mitosis enables them to maintain low-grade aneuploidy which fosters tumor growth and progression. On the other hand, maintaining enhanced centrosomal clustering as in the case of high-grade tumors may serve cancer cells by helping them attain a more aggressive phenotype . We believe that centrosome clustering confers cancer cells with cytoskeletal advantages that may enhance cell polarization, Golgi-dependent vesicular trafficking, stromal invasion, and other aspects of metastatic progression .
Since high-grade ovarian cancers display a high degree of centrosomal amplification coupled with a configuration that keeps them in a “bundled” or clustered state, it is reasonable to expect that they overexpress proteins that will aid in centrosome clustering. To this end, we examined the expression of KIFC1, a known centrosome clustering molecule, and evaluated its prognostic power in ovarian cancer. We found KIFC1 levels to be significantly higher in ovarian cancer compared to normal ovarian epithelia. It is noteworthy that increasing KIFC1 levels were associated with increasing grades. This is especially important from a clinical viewpoint as the likely course of disease could be predicted accurately by measuring KIFC1 expression in different grades. Even though KIFC1 levels in both primary and metastatic tissue were significantly higher compared to normal tissue, KIFC1 expression between matched primary and metastatic ovarian carcinoma were not different. We speculate that the primary tumor expressed KIFC1 to a high enough level to cause epithelial to mesenchymal transition (EMT) that marks the beginning of the metastatic journey. Thus, it is likely that KIFC1 may serve to be a marker for metastatic onset. The overall survival was also lower in patients with high KIFC1 expression, implying its value as a prognostic biomarker in ovarian cancer. We further delineated the proteins that might be interacting with KIFC1 gene and found that many of these potentially interacting proteins were cell-cycle related genes. KIFC1 is a kinesin involved in various cellular processes such as mitotic spindle assembly , centrosome clustering [25, 26], and vesicle transport  in cancer cells. Thus, its interaction with an array of cell cycle-specific proteins is ostensible. Nevertheless, these data implicate the role of KIFC1 in the regulation of cell-cycle and interaction with these proteins. It is likely that KIFC1 has clustering-independent role in cancer cells which require further investigation.
Taken together, we demonstrate that ovarian cancers display amplified centrosomes, a feature of aggressive tumors. To cope up with the abnormal centrosomal load and at the same time circumvent mitotic catastrophe, ovarian cancer cells upregulate genes like KIFC1 that are known to induce centrosome clustering, a “tactic” that tumor cells have evolved to execute mitosis in a pseudo-bipolar state. Our data compellingly underscores that KIFC1 can be a prognostic biomarker in ovarian cancers. Our interactome data have discovered some “purely” novel potential binding partners based upon pathway connectivity, which merit further screening and investigation to shed more light into the possibly clustering independent roles of KIFC1.
Material and methods
A) In silico analysis of KIFC1 gene expression
A.I. Data collection
List of Gene ID
Normal samples GEO Series ID
Cancer samples GEO Series ID
GSE14407, GSE18520 (N = 38)
GSE20565, GSE14764, GSE12418, GSE41498, GSE9890, GSE9891 (N = 494)
The gene expression data for the fourteen matched sets of primary and metastatic (omental) samples were procured from a recently published paper from John McDonald’s group .
A.II. Data pre-processing
One channel micro array data was Mas5.0 normalized , and was further taken for processing.
A.III. Identification of KIFC1 gene expression
Logarithm to the base 2 transformed KIFC1 expression levels of ovarian cancer patients were extracted from the TCGA and GEO patients and compared to their normal pairs. Identification of KIFC1 expression levels for 14 primary and secondary ovarian cancer samples was done following AI, AII, and AIII protocols.
A.IV. KIFC1 gene expression and associated clinical outcome information
Associated grade, stage and clinical information for 468 patients from GEO database were plotted for KIFC1 gene expression values.
B) Interactions of KIFC1 protein
KIFC1 gene was correlated with 22277 affymetrix probe id’s, genes with a correlation value of >0.5 and between -0.1 and -0.5 were fed in the IPA tool to identify pathways associated. Proteins and their interactions amongst them with confidence values were extracted from STRING 9.1, and inputted in Cytoscape 3.0 for building interactomes.
C) tissue specimens
In this study, a total of 14 human serous ovarian cancer specimens and their normal adjacent were procured from Northside Hospital (Atlanta, GA). There were 7 each of low- and high- grade cancer tissues.
Slides were first deparaffinized by baking in oven at 60°C for 2 h followed by 3 xylene baths. Rehydration was then performed in a series of ethanol baths (100%, 90%, 75% and 50%). Antigen retrieval was achieved by citrate buffer (pH 6.0) in a pressure-cooker (15 psi) for 3 min. Primary antibodies (1:2000 dilution) were incubated with the slides for 45 min at 37°C. The cells were washed 10× with PBS at room temperature before incubating at 37°C with a 1:2000 dilution of conjugated secondary antibodies. Cells were washed 5× with PBS and then mounted with Prolong-Gold antifade reagent that contained DAPI (Invitrogen).
E) Statistical analysis
Statistical analysis was performed using Student’s t-test and the criteria for statistical significance was p < 0.05.
This study was supported by grants to RA from the National Cancer Institute at the National Institutes of Health (R00CA131489, R01 CA169127).
- Ogden A, Rida PC, Aneja R: Let’s huddle to prevent a muddle: centrosome declustering as an attractive anticancer strategy. Cell Death Differ 2012,19(8):1255–1267. 10.1038/cdd.2012.61PubMed CentralView ArticlePubMedGoogle Scholar
- Ogden A, Rida PC, Aneja R: Heading off with the herd: how cancer cells might maneuver supernumerary centrosomes for directional migration. Cancer Metastasis Rev 2013,32(1–2):269–287.PubMed CentralView ArticlePubMedGoogle Scholar
- Pannu V, Rida PC, Ogden A, Clewley R, Cheng A, Karna P, Lopus M, Mishra RC, Zhou J, Aneja R: Induction of robust de novo centrosome amplification, high-grade spindle multipolarity and metaphase catastrophe: a novel chemotherapeutic approach. Cell Death Dis 2012, 3: e346. 10.1038/cddis.2012.82PubMed CentralView ArticlePubMedGoogle Scholar
- Mahjoub MR: The importance of a single primary cilium. Organogenesis 2013,9(2):61–69. 10.4161/org.25144PubMed CentralView ArticlePubMedGoogle Scholar
- Fukasawa K: Oncogenes and tumour suppressors take on centrosomes. Nat Rev Cancer 2007,7(12):911–924. 10.1038/nrc2249View ArticlePubMedGoogle Scholar
- Acilan C, Saunders WS: A tale of too many centrosomes. Cell 2008,134(4):572–575. 10.1016/j.cell.2008.08.007View ArticlePubMedGoogle Scholar
- Quintyne NJ, Reing JE, Hoffelder DR, Gollin SM, Saunders WS: Spindle multipolarity is prevented by centrosomal clustering. Science 2005,307(5706):127–129. 10.1126/science.1104905View ArticlePubMedGoogle Scholar
- Fukasawa K: Centrosome amplification, chromosome instability and cancer development. Cancer Lett 2005,230(1):6–19. 10.1016/j.canlet.2004.12.028View ArticlePubMedGoogle Scholar
- Kramer A, Maier B, Bartek J: Centrosome clustering and chromosomal (in) stability: a matter of life and death. Mol Oncol 2011,5(4):324–335. 10.1016/j.molonc.2011.05.003View ArticlePubMedGoogle Scholar
- Wu J, Mikule K, Wang W, Su N, Petteruti P, Gharahdaghi F, Code E, Zhu X, Jacques K, Lai Z, Yang B, Lamb ML, Chuaqui C, Keen N, Chen H: Discovery and mechanistic study of a small molecule inhibitor for motor protein KIFC1. ACS Chem Biol 2013,8(10):2201–2208. 10.1021/cb400186wView ArticlePubMedGoogle Scholar
- Watts CA, Richards FM, Bender A, Bond PJ, Korb O, Kern O, Riddick M, Owen P, Myers RM, Raff J, Gergely F, Jodrell DI, Ley SV: Design, synthesis, and biological evaluation of an allosteric inhibitor of HSET that targets cancer cells with supernumerary centrosomes. Chem Biol 2013,20(11):1399–1410. 10.1016/j.chembiol.2013.09.012PubMed CentralView ArticlePubMedGoogle Scholar
- Grinberg-Rashi H, Ofek E, Perelman M, Skarda J, Yaron P, Hajduch M, Jacob-Hirsch J, Amariglio N, Krupsky M, Simansky DA, Ram Z, Pfeffer R, Galernter I, Steinberg DM, Ben-Dov I, Rechavi G, Izraeli S: The expression of three genes in primary non-small cell lung cancer is associated with metastatic spread to the brain. Clin Cancer Res 2009,15(5):1755–1761. 10.1158/1078-0432.CCR-08-2124View ArticlePubMedGoogle Scholar
- Hsu LC, Kapali M, DeLoia JA, Gallion HH: Centrosome abnormalities in ovarian cancer. Int J Cancer 2005,113(5):746–751. 10.1002/ijc.20633View ArticlePubMedGoogle Scholar
- Bayani J, Paderova J, Murphy J, Rosen B, Zielenska M, Squire JA: Distinct patterns of structural and numerical chromosomal instability characterize sporadic ovarian cancer. Neoplasia 2008,10(10):1057–1065.PubMed CentralView ArticlePubMedGoogle Scholar
- Chan JY: A clinical overview of centrosome amplification in human cancers. Int J Biol Sci 2011,7(8):1122–1144.PubMed CentralView ArticlePubMedGoogle Scholar
- Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res 2002,30(1):207–210. 10.1093/nar/30.1.207PubMed CentralView ArticlePubMedGoogle Scholar
- Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, Haussler D: The human genome browser at UCSC. Genome Res 2002,12(6):996–1006. 10.1101/gr.229102. Article published online before print in May 2002PubMed CentralView ArticlePubMedGoogle Scholar
- Chu CS, Rubin SC: Screening for ovarian cancer in the general population. Best Pract Res Clin Obstet Gynaecol 2006,20(2):307–320. 10.1016/j.bpobgyn.2005.10.012View ArticlePubMedGoogle Scholar
- Singh N, Gilks CB, Wilkinson N, McCluggage WG: Assignment of primary site in high-grade serous tubal, ovarian and peritoneal carcinoma: a proposal. Histopathology 2014. epub ahead of printGoogle Scholar
- Lili LN, Matyunina LV, Walker LD, Wells SL, Benigno BB, McDonald JF: Molecular profiling supports the role of epithelial-to-mesenchymal transition (EMT) in ovarian cancer metastasis. J Ovarian Res 2013,6(1):49. 10.1186/1757-2215-6-49PubMed CentralView ArticlePubMedGoogle Scholar
- Badgwell D, Bast RC Jr: Early detection of ovarian cancer. Dis Markers 2007,23(5–6):397–410.PubMed CentralView ArticlePubMedGoogle Scholar
- Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, Lin J, Minguez P, Bork P, von Mering C, Jensen LJ: STRING v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res 2013,41(Database issue):D808-D815.PubMed CentralView ArticlePubMedGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 2003,13(11):2498–2504. 10.1101/gr.1239303PubMed CentralView ArticlePubMedGoogle Scholar
- Cai S, Weaver LN, Ems-McClung SC, Walczak CE: Proper organization of microtubule minus ends is needed for midzone stability and cytokinesis. Current biology : CB 2010,20(9):880–885. 10.1016/j.cub.2010.03.067PubMed CentralView ArticlePubMedGoogle Scholar
- Kwon M, Godinho SA, Chandhok NS, Ganem NJ, Azioune A, Thery M, Pellman D: Mechanisms to suppress multipolar divisions in cancer cells with extra centrosomes. Genes Dev 2008,22(16):2189–2203. 10.1101/gad.1700908PubMed CentralView ArticlePubMedGoogle Scholar
- Kleylein-Sohn J, Pollinger B, Ohmer M, Hofmann F, Nigg EA, Hemmings BA, Wartmann M: Acentrosomal spindle organization renders cancer cells dependent on the kinesin HSET. J Cell Sci 2012,125(Pt 22):5391–5402.View ArticlePubMedGoogle Scholar
- Mukhopadhyay A, Quiroz JA, Wolkoff AW: Rab1a regulates sorting of early endocytic vesicles. Am J Physiol Gastrointest Liver Physiol 2014,306(5):G412-G424. 10.1152/ajpgi.00118.2013PubMed CentralView ArticlePubMedGoogle Scholar
- Pepper SD, Saunders EK, Edwards LE, Wilson CL, Miller CJ: The utility of MAS5 expression summary and detection call algorithms. BMC Bioinforma 2007, 8: 273. 10.1186/1471-2105-8-273View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. 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.