Volume 16 Supplement 1
Exome Sequencing of Normal and Isogenic Transformed Human Colonic Epithelial Cells (HCECs) Reveals Novel Genes Potentially Involved in the Early Stages of Colorectal Tumorigenesis
© Zhang et al; licensee BioMed Central Ltd. 2015
Published: 15 January 2015
We have generated a series of isogenically derived immortalized human colonic epithelial cell (HCEC 1CT and HCEC 2CT) lines, including parental un-immortalized normal cell strains. The CDK4 and hTERT immortalized colonic epithelial cell line (HCEC 1CT) is initially karyotypically normal diploid and expresses a series of epithelial cell markers including stem cell markers. Under stressful tissue culture conditions, a spontaneous aneuploidy event occurred in the HCEC 1CT line, resulting in a single chromosomal change leading to a stable trisomy 7 cell line (1CT7). Trisomy 7 occurs in about 40% of all benign human adenomas (polyps) and thus this specific chromosomal change in diploid HCEC 1CT cells appears to be non random. In addition, we have partially transformed the HCEC 1CT line by introducing stable knockdown of wild type APC and TP53, and ectopically introducing a mutant Krasv12 and a mutant version of APC (A1309), all commonly found mutations in colorectal cancer (CRC).
Whole exome sequencing and bioinformatic analyses were performed to comprehensively examine the genetic background of these isogenic cell lines.
Exome sequencing of these experimentally progressed cell lines recapitulates a list of genes previously reported to be involved in CRC tumorigenesis. In addition, sequencing revealed a collection of novel genes specifically detected in 1CT7 and A1309 cells but not normal diploid 1CT cells.
This study demonstrates the utility of using isogenic experimentally derived HCEC lines as a model to recapitulate CRC initiation and progression. Exome sequencing reveals a collection of novel genes that may play important roles in CRC tumorigenesis.
Colorectal cancer (CRC) is the third most commonly diagnosed cancer and third leading cause of cancer related mortality in the United States. It is well established that sporadic colorectal cancer (CRCs) arises through the acquisition of a series of sequential genetic mutations in both tumor suppressor genes and oncogenes . Mutational activation of oncogenes together with inactivation of tumor suppressor genes (TSG) contributes to colorectal tumor formation. It has been proposed that a minimum of four sequential genetic alterations are required for colorectal cancer evolution, including one oncogene (KRAS) and three TSGs (APC, SMAD4, TP53) as the main targets . The dominant or recessive nature of these genes predict that at least seven mutations (KRAS and six additional ones) are required for complete inactivation of important TSG function . The TSG mutations occur in most tumors, whereas KRAS mutations are found in approximately 50% of sporadic adenomas and carcinomas [3, 4]. However, additional changes are required to convert a normal colonic epithelial cell into a malignant carcinoma. While most CRCs have ~100 or more genomic changes, many of these are believed to be incidental or “passenger” alterations, and it is estimated that up to 15 “driver” oncogenic changes are required for transforming into full malignancy . Many of these changes are not frequently observed in CRC and thus it remains to be determined which less frequently mutated genes are involved in CRC initiation and development.
Recent advances in next generation sequencing (NGS) technology have allowed for rapid and efficient analysis of causative mutations in rare Mendelian disorders . Several studies have demonstrated the utility of exome sequencing in identifying novel driver mutations in various cancer types [7–10]. In particular, the whole exome and even the whole genome sequencing of colorectal tumors have delineated a comprehensive mutational landscape of genetic alterations in CRC [5, 11, 12]. However, the mutational events that contribute to CRC initiation are less well-studied, partly due to the lack of appropriate cellular reagents for validating important changes. We reasoned that examination of the landscape of genomic changes as early events in CRC initiation could be determined by introduction of specific alterations in the background of normal diploid HCECs. In the present study, we applied exome sequencing on a series of isogenically-derived immortalized human colonic epithelial cell (HCEC) lines generated from the same individual with defined genetic manipulations. Analysis of the mutation spectrum of these cell lines reveal expected changes and a list of novel candidate genes that may be involved in early stage of CRC tumorigenesis.
The culture conditions of HCECs and their isogenic series have been reported elsewhere . Briefly, HCECs were maintained under 2% oxygen and 5% carbon dioxide on Primaria® (BD Biosciences, San Jose, CA) plates in 4:1 high-glucose Dulbecco modified Eagle medium/medium 199 with 2% cosmic calf serum (Hyclone, Logan, UT) plus growth supplements: epidermal growth factor (EGF; 20 ng/ml; Peprotech, Rocky Hill, NJ), hydrocortisone (1 mg/ml), insulin (10 mg/ml), transferrin (2 mg/ml), and sodium selenite (5 nM) (all Sigma, St Louis, MO).
DNA and RNA Extraction
DNA and RNA were extracted from cell lysates using a DNeasy Blood & Tissue Kit or RNeasy Plus Mini Kit (Qiagen, Valencia, CA) according to the manufacturer’s instructions. Genomic DNA was used for exome capture.
Total RNA was isolated from cells using RNeasyMinikit (Qiagen, Chatsworth, CA) according to the manufacturer's protocol. Then 1 µg RNA was converted to cDNA using a First Strand cDNA Synthesis Kit (Roche, Indianapolis, IN). Real-time quantitative PCR reactions were set up in triplicate with Ssofast Master Mix (Biorad, Hercules, CA) and run on a LightCycler® 480 (Roche, Indianapolis, Indiana).
PCR was performed on cDNA from each cell line and purified PCR products were directly sequenced. Each read was aligned with reference sequence at Nucleotide BLAST website (http://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastn&PAGE_TYPE=BlastSearch&LINK_LOC=blasthome ). The forward and reverse primers for INCENP are: 5’-tctgcagggcagcaagag-3’ and 5’-tcctccttcatctgctccac-3’.
Exome capture using 3 µg of genomic DNA from each cell line was performed using the TargetSeq (TM) Exome Enrichment system (A14061) from Life Technologies according to the manufacturer's protocol. Sequencing was performed on the SOLiD(TM) 5500XL platform. Mapping to the hg19 version of the human genome and single nucleotide variations as well as small indels identification was performed using default settings of the LifeScope software (Life Technologies, Carlsbad, CA). High quality variants (with coverage >=10x and MQV>=20) were annotated and filtered using the SNP and Variation Suite (SVS) version 7 from Golden Helix. Novel and rare variants (with MAF <1%) were filtered against the NHLBI exome project database. SNVs were predicted damaging using the SIFT, Poly-Phen or the Mutation Taster software within the SVS7 pipeline.
Characteristics of the sequenced HCEC lines
The HCEC 1CT line used in these experiments was derived from non-malignant colonic tissue from a patient with a previous history of CRC who was undergoing routine colonoscopy screening. The cells derived from explants were immortalized with ectopic expression of CDK4 and hTERT as previously described . This cell line maintains a stable normal karyotype (46, XY) when continuously propagated in 2% oxygen and medium containing 2% serum . 1CT7 cells were spontaneously generated from 1CT cells after prolonged passage under serum-free condition . Trisomy in chromosome 7 is one of the earliest events occurring in up to ~ 40% of colonic benign adenomas [15–17]. 1CT7 cells have enhanced cell migration (in a scratch-wound assay) compared to 1CT cells when cultured under low (2%) oxygen conditions (data not shown). Additionally, 1CT7 cells have significant up-regulation of EGFR and c-Met, which are two chromosome 7-located receptor tyrosine kinase compared to 1CT cells . 1CTRPA A1309 (abbreviated as A1309) is a partially transformed cell line harboring TP53 and APC knockdowns (>90%), as well as ectopic expression of oncogenic KRASV12 and truncated APC1309, all of which are common mutations detected in CRC tumors [1, 2]. This cell line exhibits enhancement in cellular proliferation, anchorage independent growth as well as invasion through Matrigel® compared with the 1CT line which does not have any detectable tumorigenic characteristics. Both 1CT7 and A1309 cell line are not fully transformed because they lack the ability to form tumors in immunocompromised mice (data not shown).
Exome capture and sequencing results
Summary of exome sequencing results of three isogenic HCEC cell lines.
Total yield (bp)
Mappable yield (bp)
On-target yield (bp)
Coverage of target region (more than 10x)
Mean read depth of targeted region
Mean read depth of called variants
Number of total variants
Number of coding variants
Number of missense, nonsense, splice, and indel variants
After NHLBI ESP6500 Exomes filtering
Mutation spectrum of the isogenic 1CT cell lines
Previously, exome sequencing of 24 randomly selected colorectal adenomas revealed mutations involved in multiple known CRC related pathways, such as Wnt signaling, cadherin signaling, integrin signaling, inflammation, and angiogenesis . Comparison of cell autonomous specific mutations found in the present study shows the overlap of a subset of the genes implicated in these pathways (Additional file 2). Additionally, the distribution of all the high confidence hits in the present study exhibited similar pattern of biological processes to that of mutations detected in those adenomas (Additional file 3). Among these activities, metabolic, cell communication and transport are the most highly represented processes. Interestingly, a subset of mutations detected in 1CT7 or A1309 cells clustered on chromosome 11p15 (Additional file 4), consistent with a previous report of numerous aberrations detected on chromosome 11 in colorectal adenocarcinoma . Taken together, these results can be interpreted to suggest that existence of trisomy 7 and the other introduced genetic alterations lead to the acquisition of additional mutations that may drive CRC initiation and progression.1CT7 and partially transformed A1309 cells may harbor the genetic background mimicking early stages CRC. In addition, since many of these mutations are detected sequentially in an experimental in vitro manipulated setting, it suggests these mutations occur in a cell autonomous manner and are not dependent on the extracellular microenvironment that occurs in vivo.
Identification of novel candidate genes involved in CRC tumorigenesis
Identification of INCENP polymorphism in 1CT isogenic series
Characteristics of INCENP variants detected in this study.
Amino acid change
Single AA change
Single AA change
We performed whole exome sequencing of a series of isogenically derived human colonic epithelial cell lines (HCECs), including the non-cancerous diploid parental 1CT cells, the 1CT7 cells with spontaneously occurring trisomy 7 which is frequently observed as an early event in CRC, and partially transformed A1309 cells harboring commonly found mutations in CRC. On average, ~60 % of the bases were covered to 10X within the targeted bases with over 10,000 variants detected in these samples. The reason we chose 1CT as the control cell line but not 2CT or another known independent cell line is because 1CT, 1CT7 and A1309 are isogenically derived from the original pre-immortalized patient cells. Thus, the genes specifically mutated in 1CT7 and A1309 cells are more likely to be candidate “driver” genes instead of “passengers” involved in CRC tumorigenesis.
Based on the TCGA datasets examined, the mutations unique to 1CT7 occurs in 30.4% of all CRC cases whereas A1309 specific mutated genes are altered in 73.6% of all CRC cases. 1CT7 is a premalignant cell line containing only one early molecular change that occurs in about 40% of CRC cases. Additionally, many of the mutations found in this cell line are likely to be incidental events. Therefore, it is very likely amplification of chromosome 7 is an important early event in a reasonable fraction of sporadic CRC. The top ranked genes, PTPRT and CSMD1, which are unique to 1CT7 and A1309 cells, respectively, have previously been reported to be mutated in colorectal tumors [7, 19]. Comparison of our sequencing data with a previous exome sequencing study for 24 colorectal adenomas reveals the overlap of a subset of genetic mutations involved in CRC related pathways. These results suggest that the existence of trisomy 7 and the introduction of other genetic manipulation can lead to acquisition of additional genetic mutations that may contribute to CRC progression. Interestingly, knockdown of TP53 and expression of K-RasV12 in 1CT7 cells results in the emergence of trisomy 20, another nonrandom aneuploidy observed in ~85% of CRC . Therefore, 1CT7 and partially transformed A1309 cells may harbor the genetic background mimicking susceptibility to early stage colon cancer initiation and progression. These cell lines represent an ideal cell autonomous model to delineate the molecular events that contribute to CRC tumorigenesis.
Utilizing the ToppGene portal, we prioritized the candidate gene list based on protein-protein interactions and other functional annotations. A total of 27 genes are putative CRC hits using more than three methods. Many of these genes are frequently mutated in CRC tumors and patients with alteration in these genes exhibit overall poorer survival. Network analysis of these genes reveals additional and perhaps novel interactors that are also altered in CRC tumors. Therefore, this set of genes and their interacting partners may play important role in CRC tumorigenesis.
Epigenetic regulation of gene silencing is another pathway by which tumor suppressor genes are inactivated . Aberrant DNA methylation has been reported to contribute to colon cancer progression through CpG Island Methylator Phenotype (CIMP) [29, 30]. We can speculate that 1CT7 and A1309 cell lines may harbor a higher presence of aberrantly methylated genes compared to their normal isogenic counterpart, 1CT cells. Future investigation of the epigenetic signatures of these cell lines compared to preimmortalized normal epithelial cells as well as the authentic CRC samples are warranted.
In conclusion, the present study revealed the comprehensive mutation spectrums of a series of isogenically-derived HCEC lines. This has led to the identification of known CRC genes as well as a collection of novel candidate CRC genes, demonstrating the potential of utilizing these isogenic HCEC lines to unravel the early cell autonomous events that contribute to CRC initiation and progression. These newly identified important CRC “driver” genes can be potentially utilized as biomarkers for the diagnostic and prognostic applications. A collection of these candidate genes may be further pursued as novel therapeutic targets for CRC prevention and intervention.
Availability of supporting data
Raw sequencing data can be retrieved from DOI: 10.6070/H44M92HV.
List of abbreviations used
Human colonic epithelial cell
CpG island methylator phenotype
Chromosomal passenger complex
CUB and shushi multiple domain 1
Next generation sequencing
Protein tyrosine phosphatase receptor type T
Single nucleotide polymorphisms
Tumor suppressor gene
This work was supported by CPRIT Training Grant #RP130189 to LZ, NASA grants NNX09AU95G, NNX11AC15G, and NNX11AC54G to JWS and support from Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, Saudi Arabia, project no (HiCi-143216-4). This work was performed in laboratories constructed with support from NIH grant C06 RR30414.
Publication fees of this article has been funded by Centre of Excellence in Genomic Medicine Research (CEGMR), King Abdulaziz University, Jeddah, Saudi Arabia.
This article has been published as part of BMC Genomics Volume 16 Supplement 1, 2015: Selected articles from the 2nd International Genomic Medical Conference (IGMC 2013): Genomics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/16/S1
- Kinzler KW, Vogelstein B: Lessons from hereditary colorectal cancer. Cell. 1996, 87 (2): 159-170. 10.1016/S0092-8674(00)81333-1.View ArticlePubMedGoogle Scholar
- Fodde R, Smits R, Clevers H: APC, signal transduction and genetic instability in colorectal cancer. Nature reviews Cancer. 2001, 1 (1): 55-67. 10.1038/35094067.View ArticlePubMedGoogle Scholar
- Forrester K, Almoguera C, Han K, Grizzle WE, Perucho M: Detection of high incidence of K-ras oncogenes during human colon tumorigenesis. Nature. 1987, 327 (6120): 298-303. 10.1038/327298a0.View ArticlePubMedGoogle Scholar
- Bos JL, Fearon ER, Hamilton SR, Verlaan-de Vries M, van Boom JH, van der Eb AJ, Vogelstein B: Prevalence of ras gene mutations in human colorectal cancers. Nature. 1987, 327 (6120): 293-297. 10.1038/327293a0.View ArticlePubMedGoogle Scholar
- Wood LD, Parsons DW, Jones S, Lin J, Sjoblom T, Leary RJ, Shen D, Boca SM, Barber T, Ptak J, et al: The genomic landscapes of human breast and colorectal cancers. Science. 2007, 318 (5853): 1108-1113. 10.1126/science.1145720.View ArticlePubMedGoogle Scholar
- Ng SB, Buckingham KJ, Lee C, Bigham AW, Tabor HK, Dent KM, Huff CD, Shannon PT, Jabs EW, Nickerson DA, et al: Exome sequencing identifies the cause of a mendelian disorder. Nature genetics. 2010, 42 (1): 30-35. 10.1038/ng.499.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang K, Kan J, Yuen ST, Shi ST, Chu KM, Law S, Chan TL, Kan Z, Chan AS, Tsui WY, et al: Exome sequencing identifies frequent mutation of ARID1A in molecular subtypes of gastric cancer. Nature genetics. 2011, 43 (12): 1219-1223. 10.1038/ng.982.View ArticlePubMedGoogle Scholar
- Yan XJ, Xu J, Gu ZH, Pan CM, Lu G, Shen Y, Shi JY, Zhu YM, Tang L, Zhang XW, et al: Exome sequencing identifies somatic mutations of DNA methyltransferase gene DNMT3A in acute monocytic leukemia. Nature genetics. 2011, 43 (4): 309-315. 10.1038/ng.788.View ArticlePubMedGoogle Scholar
- Wei X, Walia V, Lin JC, Teer JK, Prickett TD, Gartner J, Davis S, Program NCS, Stemke-Hale K, Davies MA, et al: Exome sequencing identifies GRIN2A as frequently mutated in melanoma. Nature genetics. 2011, 43 (5): 442-446. 10.1038/ng.810.PubMed CentralView ArticlePubMedGoogle Scholar
- Varela I, Tarpey P, Raine K, Huang D, Ong CK, Stephens P, Davies H, Jones D, Lin ML, Teague J, et al: Exome sequencing identifies frequent mutation of the SWI/SNF complex gene PBRM1 in renal carcinoma. Nature. 2011, 469 (7331): 539-542. 10.1038/nature09639.PubMed CentralView ArticlePubMedGoogle Scholar
- Sjoblom T, Jones S, Wood LD, Parsons DW, Lin J, Barber TD, Mandelker D, Leary RJ, Ptak J, Silliman N, et al: The consensus coding sequences of human breast and colorectal cancers. Science. 2006, 314 (5797): 268-274. 10.1126/science.1133427.View ArticlePubMedGoogle Scholar
- Bass AJ, Lawrence MS, Brace LE, Ramos AH, Drier Y, Cibulskis K, Sougnez C, Voet D, Saksena G, Sivachenko A, et al: Genomic sequencing of colorectal adenocarcinomas identifies a recurrent VTI1A-TCF7L2 fusion. Nature genetics. 2011, 43 (10): 964-968. 10.1038/ng.936.PubMed CentralView ArticlePubMedGoogle Scholar
- Roig AI, Eskiocak U, Hight SK, Kim SB, Delgado O, Souza RF, Spechler SJ, Wright WE, Shay JW: Immortalized epithelial cells derived from human colon biopsies express stem cell markers and differentiate in vitro. Gastroenterology. 2010, 138 (3): 1012-1021. 10.1053/j.gastro.2009.11.052. e1011-1015View ArticlePubMedGoogle Scholar
- Ly P, Eskiocak U, Kim SB, Roig AI, Hight SK, Lulla DR, Zou YS, Batten K, Wright WE, Shay JW: Characterization of aneuploid populations with trisomy 7 and 20 derived from diploid human colonic epithelial cells. Neoplasia. 2011, 13 (4): 348-357. 10.1593/neo.101580.PubMed CentralView ArticlePubMedGoogle Scholar
- Habermann JK, Paulsen U, Roblick UJ, Upender MB, McShane LM, Korn EL, Wangsa D, Kruger S, Duchrow M, Bruch HP, et al: Stage-specific alterations of the genome, transcriptome, and proteome during colorectal carcinogenesis. Genes, chromosomes & cancer. 2007, 46 (1): 10-26. 10.1002/gcc.20382.View ArticleGoogle Scholar
- Ried T, Knutzen R, Steinbeck R, Blegen H, Schrock E, Heselmeyer K, du Manoir S, Auer G: Comparative genomic hybridization reveals a specific pattern of chromosomal gains and losses during the genesis of colorectal tumors. Genes, chromosomes & cancer. 1996, 15 (4): 234-245. 10.1002/(SICI)1098-2264(199604)15:4<234::AID-GCC5>3.0.CO;2-2.View ArticleGoogle Scholar
- Bomme L, Lothe RA, Bardi G, Fenger C, Kronborg O, Heim S: Assessments of clonal composition of colorectal adenomas by FISH analysis of chromosomes 1, 7, 13 and 20. International journal of cancer Journal international du cancer. 2001, 92 (6): 816-823. 10.1002/ijc.1275.View ArticlePubMedGoogle Scholar
- Wang Z, Shen D, Parsons DW, Bardelli A, Sager J, Szabo S, Ptak J, Silliman N, Peters BA, van der Heijden MS, et al: Mutational analysis of the tyrosine phosphatome in colorectal cancers. Science. 2004, 304 (5674): 1164-1166. 10.1126/science.1096096.View ArticlePubMedGoogle Scholar
- Farrell C, Crimm H, Meeh P, Croshaw R, Barbar T, Vandersteenhoven JJ, Butler W, Buckhaults P: Somatic mutations to CSMD1 in colorectal adenocarcinomas. Cancer biology & therapy. 2008, 7 (4): 609-613. 10.4161/cbt.7.4.5623.View ArticleGoogle Scholar
- Masica DL, Karchin R: Correlation of somatic mutation and expression identifies genes important in human glioblastoma progression and survival. Cancer research. 2011, 71 (13): 4550-4561. 10.1158/0008-5472.CAN-11-0180.PubMed CentralView ArticlePubMedGoogle Scholar
- Lee W, Jiang Z, Liu J, Haverty PM, Guan Y, Stinson J, Yue P, Zhang Y, Pant KP, Bhatt D, et al: The mutation spectrum revealed by paired genome sequences from a lung cancer patient. Nature. 2010, 465 (7297): 473-477. 10.1038/nature09004.View ArticlePubMedGoogle Scholar
- Nikolaev SI, Sotiriou SK, Pateras IS, Santoni F, Sougioultzis S, Edgren H, Almusa H, Robyr D, Guipponi M, Saarela J, et al: A single-nucleotide substitution mutator phenotype revealed by exome sequencing of human colon adenomas. Cancer research. 2012, 72 (23): 6279-6289. 10.1158/0008-5472.CAN-12-3869.View ArticlePubMedGoogle Scholar
- Tagawa Y, Sawai T, Nakagoe T, Morinaga M, Yasutake T, Ayabe H, Tomita M: Numerical aberrations of chromosomes 11 and 17 in colorectal adenocarcinomas. Surgery today. 1996, 26 (11): 869-874. 10.1007/BF00311787.View ArticlePubMedGoogle Scholar
- Chen J, Bardes EE, Aronow BJ, Jegga AG: ToppGene Suite for gene list enrichment analysis and candidate gene prioritization. Nucleic acids research. 2009, 37 (Web Server issue): W305-311.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen J, Xu H, Aronow BJ, Jegga AG: Improved human disease candidate gene prioritization using mouse phenotype. BMC bioinformatics. 2007, 8: 392-10.1186/1471-2105-8-392.PubMed CentralView ArticlePubMedGoogle Scholar
- Carmena M, Wheelock M, Funabiki H, Earnshaw WC: The chromosomal passenger complex (CPC): from easy rider to the godfather of mitosis. Nature reviews Molecular cell biology. 2012, 13 (12): 789-803. 10.1038/nrm3474.PubMed CentralView ArticlePubMedGoogle Scholar
- Adams RR, Eckley DM, Vagnarelli P, Wheatley SP, Gerloff DL, Mackay AM, Svingen PA, Kaufmann SH, Earnshaw WC: Human INCENP colocalizes with the Aurora-B/AIRK2 kinase on chromosomes and is overexpressed in tumour cells. Chromosoma. 2001, 110 (2): 65-74. 10.1007/s004120100130.View ArticlePubMedGoogle Scholar
- Wheatley SP, Carvalho A, Vagnarelli P, Earnshaw WC: INCENP is required for proper targeting of Survivin to the centromeres and the anaphase spindle during mitosis. Current biology : CB. 2001, 11 (11): 886-890. 10.1016/S0960-9822(01)00238-X.View ArticlePubMedGoogle Scholar
- Lao VV, Grady WM: Epigenetics and colorectal cancer. Nature reviews Gastroenterology & hepatology. 2011, 8 (12): 686-700. 10.1038/nrgastro.2011.173.View ArticleGoogle Scholar
- Goel A, Boland CR: Epigenetics of colorectal cancer. Gastroenterology. 2012, 143 (6): 1442-1460. 10.1053/j.gastro.2012.09.032. e1441PubMed CentralView ArticlePubMedGoogle 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 cited. 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.