- Open Access
GenomeCAT: a versatile tool for the analysis and integrative visualization of DNA copy number variants
BMC Bioinformatics volume 18, Article number: 19 (2017)
The analysis of DNA copy number variants (CNV) has increasing impact in the field of genetic diagnostics and research. However, the interpretation of CNV data derived from high resolution array CGH or NGS platforms is complicated by the considerable variability of the human genome. Therefore, tools for multidimensional data analysis and comparison of patient cohorts are needed to assist in the discrimination of clinically relevant CNVs from others.
We developed GenomeCAT, a standalone Java application for the analysis and integrative visualization of CNVs. GenomeCAT is composed of three modules dedicated to the inspection of single cases, comparative analysis of multidimensional data and group comparisons aiming at the identification of recurrent aberrations in patients sharing the same phenotype, respectively. Its flexible import options ease the comparative analysis of own results derived from microarray or NGS platforms with data from literature or public depositories. Multidimensional data obtained from different experiment types can be merged into a common data matrix to enable common visualization and analysis. All results are stored in the integrated MySQL database, but can also be exported as tab delimited files for further statistical calculations in external programs.
GenomeCAT offers a broad spectrum of visualization and analysis tools that assist in the evaluation of CNVs in the context of other experiment data and annotations. The use of GenomeCAT does not require any specialized computer skills. The various R packages implemented for data analysis are fully integrated into GenomeCATs graphical user interface and the installation process is supported by a wizard. The flexibility in terms of data import and export in combination with the ability to create a common data matrix makes the program also well suited as an interface between genomic data from heterogeneous sources and external software tools. Due to the modular architecture the functionality of GenomeCAT can be easily extended by further R packages or customized plug-ins to meet future requirements.
DNA copy number variants represent the greatest source of genetic variability in humans  and are the underlying cause of many human diseases. Array CGH is recognized as a first-tier test for DNA copy number variants (CNV)  and accordingly, many laboratories have already established their pipelines for pre-processing of array CGH data and CNV calling. In many cases these pipelines are based on software packages provided by the companies selling DNA microarrays or scanners such as BlueFuse , CytoSure  or CytoGenomics . Yet, the scope of these tools is focused on the identification of CNVs and their evaluation in the context of gene content and frequency of a given variant in the healthy population. Comparative analysis, which integrates data obtained from multiple patients, or other experiment types are hardly supported, in particular when they are based on different array platforms or NGS technology.
Such kind of meta-analysis needs the implementation of additional commercial or free software. Each of the currently existing software solutions have their particular strength and focus. Some are particularly useful for the identification of genomic regions significantly associated with a given phenotype [4, 6–12] or have implemented algorithms specifically designed to detect and query copy number changes in SNP data sets [8, 13]. Others provide a gene centered view on copy number aberrations [14, 15] or examine CNVs in a clinical context . Only a few free software packages offer a comprehensive spectrum of visualization and analysis tools for multidimensional array data operable via a graphical user interface [12–17]. What these tools have in common is that they have been designed with the intention to analyze microarray data. NGS data are usually displayed in alternative data browsers such as the Integrative Genome Viewer - IGV , or the Integrated Genome Browser – IGB . These browsers also support visualization of array data when present in the appropriate format. However, as in the case of the IGV, analysis of array data that goes beyond visualization requires the export to the GenePattern software , where several web-based features for DNA copy number analysis are provided.
In light of the increasing relevance of multi-dimensional data analysis several commercial softwares have been brought to market, including Partek , GenomicWorkbench , Genedata Expressionist for Genomic Profiling , Array Studio , GenomeStudio , CGH Fusion , Nexus Expression , CLC Workbench  and Subio . Yet, these programs are neither open source nor free in most instances. Thus, considerable licensing fees have to be paid and advancement of this software is solely dependent on the company.
Proceeding on the experiences with our previous analysis software CGHPRO , we aimed to create a versatile tool that facilitates the meta-analysis of array CGH results and corresponding data from other experiment types and platforms. We designed GenomeCAT under the premise that it is easy to install and use, and offers a broad spectrum of flexible visualization and analysis options without the need of specialized computer skills or the obligation to upload sensible patient data to web servers.
GenomeCAT is a desktop application developed in Java using the NetBeans Platform. It is an open source software and is provided as a free download. The program has a modular structure, which supports the program-aided updating and the implementation of new plug-ins. At the center of the program is a MySQL database, designed to maintain experiment data, metadata and annotation tracks. The current version refers to the human genome only, but the database is designed to be adaptable to any other genome when necessary (Fig. 1). An installation wizard guides through the installation process that comprises the set-up of the desktop application, the MySQL database and an R environment.
The execution of R packages is embedded in the desktop application. Users can enter and modify method-specific parameters via a dialog box. The R packages themselves run as a background process, the progress of which is reported on the screen (Additional file 1). Results produced by the R session are automatically stored into GenomeCAT. The design of this interface eases the future addition of R packages in order to further increase the functionality of GenomeCAT.
Integration of heterogeneous data
Multidimensional data are frequently produced by means of different experiment platforms. This implies that comparative analysis has to be preceded by the creation of a common data matrix. GenomeCAT is capable to address this issue in three ways: data binning based on annotation attributes (for example genes), genomic bins of variable size or user defined intervals. The resulting data matrix can not only be accessed within GenomeCAT, but can also be exported for further analysis in external programs (e.g., to visualize data as interactive heatmaps in Gitools  or as network attributes in Cytoscape ).
The search for features overlapping with a given genomic interval is a recurrent procedure in multidimensional data analysis. This applies also to the mapping to annotation attributes or genomic intervals as described above. In order to accelerate these queries in GenomeCAT we took advantage of the Spatial Index, an extension of the MySQL database. There, genomic locations or intervals are stored as geometric objects, which are indexed via R-Trees . Based on this indexing scheme we employ the MBRIntersect function [34, 35] for queries and filters, which speeds up the processing time by a factor of 4 on average compared to the use of a composite index. In contrast to other software packages  in GenomeCAT the number of cases that can be simultaneously analyzed is not confined. It is only limited by the heapsize of the Java application. Loading of data and computational intensive calculations are parallelized in order to optimally exploit the potential of the multi-core architecture of modern CPUs.
Results and Discussion
GenomeCAT supports different ways of data import. While users can choose the traditional way – to import an array platform and add sample data afterwards - the preferred route may be the direct import of data formatted in BED style (chromosome start stop score). This format is simple, platform independent and enables more flexible entry points into the analysis with GenomeCAT. Thus users can stick to their well-established pipeline for primary data analysis and start the analysis in our software with a set of already predefined CNVs. Moreover, the format is ideal for comparisons of own array CGH data with results from other experiment types or, for example, CNVs that have been reported in literature as a list of genomic intervals. For maximal flexibility, our software also offers the option to interactively compile the necessary data from more comprehensive tables.
In addition, GenomeCAT has import routines for experiment data available in GEOs SOFT file format and can also process BAM files for the import of NGS data. All data are stored in a MySQL database together with metadata such as phenotype information. Data are organized in a hierarchical structure that is searchable and can be filtered by various criteria.
Module 1: single view
Single case analysis can be accomplished with SingleView, the first of the three modules that make up GenomeCAT. In this module users can display array CGH profiles as familiar ratio plots along the chromosome ideograms. Optionally, annotation tracks such as GC content, CNVs from the Database of Genomic Variants  and segmental duplications  can be depicted beside (Fig. 2). The layout can be customized and the plots are interactive. For example, genomic coordinates are provided as mouseover event and regions can be zoomed-in or directly viewed in the UCSC Genome Browser . CNV calling can be performed by means of customizable fixed and dynamic threshold settings and in combination with CBS  or HMM . For each genomic region extracted in this way, GenomeCAT calculates a quality score by dividing the averaged ratio value within each extracted region by the median average deviation. If CBS has been employed before, the median average deviation is calculated based on the deviation from the averaged ratio value within the segments as defined by CBS. Moreover, Ringo  is implemented to facilitate peak finding in data derived from chromatin immunoprecipitation (ChIP) experiments. All processed data are stored as separate tracks and each processing step is recorded together with its parameter settings. Thus it is easy to recapitulate the analysis procedure and to go back to the original data if necessary, or to compare results obtained with different parameter settings. All tracks generated in the course of data analysis can be exported as tab delimited file in BED or BedGraph format, which is suitable for direct visualization in the UCSC Genome Browser .
A different type of single case visualization can be performed in the Region of Interest (ROI) viewer. Proceeding on a list of user defined genomic intervals this feature sorts these intervals according to their average experiment values and displays the value distribution of each defined interval in a heatmap with a resolution of 10 bins per interval. Applications of this feature can be the pre-screening of array CGH results based on a list of genomic intervals recurrently altered in genomic disorders (Fig. 3) or the identification of genes with highest scores in ChIP experiments to name a few.
Module 2: comparative view
The second module of GenomeCAT is dedicated to the simultaneous visualization of multiple tracks. These can be the results of the same sample processed with different parameters, data from the very same patient but different array or NGS-based experiment types or data from different patients. Plots produced by this module are interactive, including the options to re-sort, to zoom-in and re-scale, and to view the intervals of interest in the UCSC Genome Browser (Fig. 4).
One practical issue that complicates the integrative analysis of data derived from different experiment types is the fact, that they usually do not share a common coordinate system. Oligonucleotides of gene expression and CGH arrays hardly overlap and both platforms are not directly comparable to NGS data. This problem is addressed by our mapping feature, which allows the creation of a common data matrix for all experiments opened in this view either based on genomic bins of selectable size, genes or custom defined intervals. The resulting table is automatically stored in the internal database, but can also be exported as tab delimited file for down-stream analysis by statistical packages or visualization in other tools. For example, mapping on genes can be used to display array CGH ratios as attributes in a Cytoscape network.
Module 3: group explorer
The considerable variability of the human genome in health and disease complicates the interpretation of CNVs or patterns of copy number alterations. Recurrence of aberrations within a group of patients with similar phenotype or differences between patient groups has proven a valuable criterion to filter for biological meaningful alterations. The third module of GenomeCAT has been designed to facilitate such group comparisons. All experiment results stored in the GenomeCAT database can be filtered and selected for simultaneous visualization based on metadata such as phenotype. Separate colors can be assigned to each group or even to individual cases. The latter option can be used to highlight particular cases to ease their identification in the overview later on. CNVs can be displayed as colored bars along the chromosomes with the option to control opacity and color saturation of a given CNV by its ratio or quality score. In this way it is possible to discriminate homozygous deletions from heterozygous ones and moderate gains of DNA from high copy amplifications, respectively. Also this graphical user interface is interactive. Clicking on individual CNVs in the plot highlights case details in the adjacent table and vice versa. Regions of interest can be zoomed-in or checked in the UCSC Genome Browser.
While this mode is well suited to present the absolute numbers of CNVs and their genomic location, the relative frequency plot - also included in this module- can be employed to compare CNV frequencies independent of group size. As demonstrated in Fig. 5, this way of visualization facilitates the recognition of phenotype specific aberrations. However, the application of this tool is not restricted to CNV analysis, but it can also be used to depict the probability of epigenetic modifications in regions frequently affected by copy number changes in a specific tumor type and so forth. Relative frequencies as calculated by GenomeCAT can be exported as CSV files for further statistical analysis or visualization in external software packages.
GenomeCAT provides comprehensive tools for the analysis of DNA copy number variants and facilitates the evaluation of their biological relevance in the context of genome annotations and results obtained from different experiment types. Its flexible import options ease the comparative analysis of own results with data from literature or public depositories. Moreover, GenomeCAT can act as an interface to other software tools since results generated in GenomeCAT can be exported in standard file formats.
Availability and requirements
Project name: GenomeCAT
Project home page: http://genomecat.github.io/genomeCATSuite
Source code: https://github.com/genomeCAT/genomeCATSuite
Operating system: Linux and Windows 7
Programming language: Java, SQL, R
Other requirements: Java 1.8 or higher, MySQL database and R
License: GNU General Public License
Any restriction to use by non-academics: Contact authors
Browser extensible data
Circular binary segmentation
Comparative genomic hybridization
Copy number variant
Hidden Markov Model
Integrated genome browser
Integrative genome viewer
Minimum bounding rectangle
Region of interest
Conrad DF, Pinto D, Redon R, Feuk L, Gokcumen O, Zhang Y, Aerts J, Andrews TD, Barnes C, Campbell P, et al. Origins and functional impact of copy number variation in the human genome. Nature. 2010;464(7289):704–12.
Riggs ER, Wain KE, Riethmaier D, Smith-Packard B, Faucett WA, Hoppman N, Thorland EC, Patel VC, Miller DT. Chromosomal microarray impacts clinical management. Clin Genet. 2014;85(2):147–53.
BlueFuse Multi software. http://www.illumina.com/clinical/clinical_informatics/bluefuse.html.
CytoSure Interpret Software. http://www.ogt.com/products/246_cytosure_interpret_software.
Beroukhim R, Getz G, Nghiemphu L, Barretina J, Hsueh T, Linhart D, Vivanco I, Lee JC, Huang JH, Alexander S, et al. Assessing the significance of chromosomal aberrations in cancer: methodology and application to glioma. Proc Natl Acad Sci U S A. 2007;104(50):20007–12.
Diskin SJ, Eck T, Greshock J, Mosse YP, Naylor T, Stoeckert Jr CJ, Weber BL, Maris JM, Grant GR. STAC: a method for testing the significance of DNA copy number aberrations across multiple array-CGH experiments. Genome Res. 2006;16(9):1149–58.
Lingjaerde OC, Baumbusch LO, Liestol K, Glad IK, Borresen-Dale AL. CGH-Explorer: a program for analysis of array-CGH data. Bioinformatics. 2005;21(6):821–2.
Morganella S, Pagnotta SM, Ceccarelli M. Finding recurrent copy number alterations preserving within-sample homogeneity. Bioinformatics. 2011;27(21):2949–56.
Fridlyand J, Snijders AM, Pinkel D, Albertson DG, Jain AN. Hidden Markov models approach to the analysis of array CGH data. J Multivar Anal. 2004;90(1):132–53.
Zhang J, Feng B. cghMCR: Find chromosome regions showing common gains/losses. R package version 1.24.0.
Cai H, Kumar N, Baudis M. arrayMap: a reference resource for genomic copy number imbalances in human malignancies. PLoS One. 2012;7(5):e36944.
Korn JM, Kuruvilla FG, McCarroll SA, Wysoker A, Nemesh J, Cawley S, Hubbell E, Veitch J, Collins PJ, Darvishi K, et al. Integrated genotype calling and association analysis of SNPs, common copy number polymorphisms and rare CNVs. Nat Genet. 2008;40(10):1253–60.
Gribov A, Sill M, Luck S, Rucker F, Dohner K, Bullinger L, Benner A, Unwin A. SEURAT: visual analytics for the integrated analysis of microarray data. BMC Med Genet. 2010;3:21.
Al-Shahrour F, Minguez P, Tarraga J, Montaner D, Alloza E, Vaquerizas JM, Conde L, Blaschke C, Vera J, Dopazo J. BABELOMICS: a systems biology perspective in the functional annotation of genome-scale experiments. Nucleic Acids Res. 2006;34(Web Server issue):W472–6.
Knierim E, Schwarz JM, Schuelke M, Seelow D. CNVinspector: a web-based tool for the interactive evaluation of copy number variations in single patients and in cohorts. J Med Genet. 2013;50(8):529–33.
La Rosa P, Viara E, Hupe P, Pierron G, Liva S, Neuvial P, Brito I, Lair S, Servant N, Robine N, et al. VAMP: visualization and analysis of array-CGH, transcriptome and other molecular profiles. Bioinformatics. 2006;22(17):2066–73.
Thorvaldsdottir H, Robinson JT, Mesirov JP. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2013;14(2):178–92.
Nicol JW, Helt GA, Blanchard Jr SG, Raja A, Loraine AE. The Integrated Genome Browser: free software for distribution and exploration of genome-scale datasets. Bioinformatics. 2009;25(20):2730–1.
Reich M, Liefeld T, Gould J, Lerner J, Tamayo P, Mesirov JP. GenePattern 2.0. Nat Genet. 2006;38(5):500–1.
CLC Genomics Workbench. http://www.clcbio.com/products/clc-genomics-workbench/.
Genedata Expressionist for Genomic Profiling. https://www.genedata.com/lp/expressionist/genomic-profiling/.
GenomeStudio Software. http://www.illumina.com/techniques/microarrays/array-data-analysis-experimental-design/genomestudio.html.
CGH Fusion. http://www.infoquant.com/cghfusion.
Nexus Expression. http://www.biodiscovery.com/software/nexus-expression/.
CLC Workbench [http://www.clcbio.com/products/clc-main-workbench/].
Subio Platform [http://www.subio.jp/products/platform].
Chen W, Erdogan F, Ropers HH, Lenzner S, Ullmann R. CGHPRO -- a comprehensive data analysis tool for array CGH. BMC Bioinformatics. 2005;6:85.
Perez-Llamas C, Lopez-Bigas N. Gitools: analysis and visualisation of genomic data using interactive heat-maps. PLoS One. 2011;6(5):e19541.
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–504.
Guttman A. R-trees: a dynamic index structure for spatial searching. SIGMOD Rec. 1984;14(2):47–57.
OpenGIS Implementation Standard for Geographic information - Simple feature access - Part 2: SQL option [http://www.opengeospatial.org/standards/sfs].
MySQL-Specific Spatial Relation Functions That Use Minimum Bounding Rectangles (MBRs) [http://dev.mysql.com/doc/refman/5.1/en/spatial-relation-functions-mysql-specific.html#function_mbrintersects].
Chi B, de Leeuw RJ, Coe BP, Ng RT, MacAulay C, Lam WL. MD-SeeGH: a platform for integrative analysis of multi-dimensional genomic data. BMC bioinformatics. 2008;9:243.
MacDonald JR, Ziman R, Yuen RK, Feuk L, Scherer SW. The Database of Genomic Variants: a curated collection of structural variation in the human genome. Nucleic Acids Res. 2014;42(Database issue):D986–92.
Bailey JA, Yavor AM, Massa HF, Trask BJ, Eichler EE. Segmental duplications: organization and impact within the current human genome project assembly. Genome Res. 2001;11(6):1005–17.
Zhu J, Sanborn JZ, Benz S, Szeto C, Hsu F, Kuhn RM, Karolchik D, Archie J, Lenburg ME, Esserman LJ, et al. The UCSC cancer genomics browser. Nat Methods. 2009;6(4):239–40.
Venkatraman ES, Olshen AB. A faster circular binary segmentation algorithm for the analysis of array CGH data. Bioinformatics. 2007;23(6):657–63.
Toedling J, Skylar O, Krueger T, Fischer JJ, Sperling S, Huber W. Ringo--an R/Bioconductor package for analyzing ChIP-chip readouts. BMC bioinformatics. 2007;8:221.
Guelen L, Pagie L, Brasset E, Meuleman W, Faza MB, Talhout W, Eussen BH, de Klein A, Wessels L, de Laat W, van Steensel B. Domain organization of human chromosomes revealed by mapping of nuclear lamina interactions. Nature. 2008;453:948–51.
We thank Alischo Ahmed, Udo Georgi and Ines Mueller for technical assistance, software testing and helpful feedback and the reviewers for their helpful comments.
This work was supported by the Deutsche Forschungsgemeinschaft (UL342/2-1; UL342/2-2) and the German Ministry of Defense.
KT was the principal programmer of GenomeCAT. VB, AS, MP and GE contributed ideas and tested the software. KT and RU designed the software and wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Additional file 1:
Implementation of R Packages in GenomeCAT. Word document demonstrating how R packages are implemented in GenomeCAT. (DOCX 23 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
About this article
Cite this article
Tebel, K., Boldt, V., Steininger, A. et al. GenomeCAT: a versatile tool for the analysis and integrative visualization of DNA copy number variants. BMC Bioinformatics 18, 19 (2017). https://doi.org/10.1186/s12859-016-1430-x
- DNA copy number variants
- Integrative visualization