BioMed Research International
Volume 2016 (2016), Article ID 7823429, 7 pages
http://dx.doi.org/10.1155/2016/7823429

Research Article

VGSC: A Web-Based Vector Graph Toolkit of Genome Synteny and Collinearity

Yiqing Xu,1,2,3 Changwei Bi,2,3 Guoxin Wu,¹ Suyun Wei,2,3 Xiaogang Dai,2,4 Tongming Yin,2,4 and Ning Ye2,3

¹School of Computer Science and Engineering, Southeast University, Nanjing, Jiangsu 211189, China
²The Southern Modern Forestry Collaborative Innovation Center, Nanjing Forestry University, Nanjing, Jiangsu 210037, China
³College of Information Science and Technology, Nanjing Forestry University, Nanjing, Jiangsu 210037, China
⁴College of Forest Resources and Environment, Nanjing Forestry University, Nanjing, Jiangsu 210037, China

Received 6 November 2015; Revised 4 January 2016; Accepted 3 February 2016

Academic Editor: Hesham H. Ali

Copyright © 2016 Yiqing Xu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Background. In order to understand the colocalization of genetic loci amongst species, synteny and collinearity analysis is a frequent task in comparative genomics research. However many analysis software packages are not effective in visualizing results. Problems include lack of graphic visualization, simple representation, or inextensible format of outputs. Moreover, higher throughput sequencing technology requires higher resolution image output. Implementation. To fill this gap, this paper publishes VGSC, the Vector Graph toolkit of genome Synteny and Collinearity, and its online service, to visualize the synteny and collinearity in the common graphical format, including both raster (JPEG, Bitmap, and PNG) and vector graphic (SVG, EPS, and PDF). Result. Users can upload sequence alignments from blast and collinearity relationship from the synteny analysis tools. The website can generate the vector or raster graphical results automatically. We also provide a java-based bytecode binary to enable the command-line execution.

1. Introduction

Synteny is the collection of contiguous genes located on the chromosome of different species. Collinearity is a particular kind of synteny in which the genes are conserved in the same order [1]. Understanding this colocalization of genetic loci amongst species is a frequent task in comparative genomics research, and it often relies on the accuracy of homology identification within or across genomes. During the evolution, eukaryotic genomes between different species reveal this synteny and collinearity in various levels [2]. There are many reasons for the structural variation of genes all over the long evolutionary history, such as whole-genome duplication (WGD), segmental duplication, inversions, and translocations [3, 4]. Genomes have been shaped and restructured dynamically. Related application includes [5] annotation of newly sequenced genomes [6], identification of conserved noncoding sequences [7], estimation of whole genome duplication events [1], prediction of chromosomal rearrangements, and the structure of ancestral genomes [8]. As a result, the procedure of synteny and collinearity analysis has become a hot topic in evolutionary biology as a standard step for elucidating the evolutionary histories of both genomes and gene families.

To meet the requirement of synteny and collinearity analysis, the majority of softwares focus on the detection and alignment of the original sequencing data. Using the traditional clustering of neighboring match of gene pairs, various softwares have been developed to match gene pairs, including ADHoRe [9], the Max-gap Clusters by Multiple Sequence Comparison (MCMuSeC) [10], and OrthoCluster [11, 12]. More recent methods apply dynamic algorithm to pairwise collinear genes chains, in which a matching system scores the adjacent collinear gene pairs, known as anchor genes, and penalizes the distance between anchor genes. This method has been implemented in software tools such as ColinearScan [13], MCScan [1], SyMAP [6], FISH [14], and CYNTENATOR [15]. Besides the pairwise collinear relationships among chromosomal regions, the multialignment (alignment of three or more regions) of collinear chromosomal regions (referred to as collinear blocks) is more important as it can reveal ancient WGD events [1] and complex chromosomal duplication/rearrangement relationships [16]. One of the early software packages providing analysis of collinearity within gene families is MicroSyn [17]. MCScan [1], Multiple Collinearity Scan, is another very popular algorithm in synteny and collinearity detection. It scans multiple genomes or subgenomes, identifies putative homologous chromosomal regions, and marks these gene regions with alignment anchors. The latest i-ADHoRe 3.0 [18] combines pairwise comparison with an iterative profile search, and it uses rigorous statistical tests to ensure that regions found are significant. All these software packages have focused on the process of data rather than downstream analysis. Many of them do not even provide visual graphic outputs.

Another class of synteny and collinearity tools works with the general-purpose genome browsers, which are softwares that allow the user to view genome annotations in the context of a reference sequence. Most of them use vector graphics to enable the scrolling and zooming through arbitrary regions of a genome. GBrowse-syn [19] is the plugin of GBrowse 2.0 [20, 21], one of the most powerful web-based applications to visualize genomic data. It allows the comparison of collinear regions of multiple genomes using the GBrowse-styled web page, in which the synteny and collinearity are displayed as traditional connection diagram. This kind of general-purpose software packages however only provides very basic drawings, as they are not designed to meet the advance visualization requirement of the synteny and collinearity representation.

As synteny and collinearity visualization becomes increasingly important, many specific software programs have been developed lately. Most of these software programs, such as SynChro [5], GSV [22], and Easyfig [23], however inherit the linear tradition in this area, which plots the synteny and collinearity relationship into lines and bars. A typical output style uses two bars for the chromosomes and lines for the colocational relationship. While it is easier and more convenient to use web based interface to generate the linear plot, it is difficult for research reporting, especially for paper pipelines. The extension package of MCScan named MCScanX [24] implements 15 utility programs for display and analyses. However, MCScanX provides a command-line based plotter with PNG output only. Another case in point is that the i-ADHoRe 3.0 [18] extends ADHoRe [9] and provides a package to draw dot plot in SVG vector graphics and PNG raster images.

Circos [25] is a well-known visualization tool using circular ideogram layout to facilitate the identification and analysis of similarities and differences found in comparisons of genomes. Raster or vector images can be created from GFF-style data inputs and hierarchical configuration files, which are popular in bioinformatics researches, making Circos suitable for rapid reporting pipelines. A typical case is C-Sibelia [26], which focuses on the synteny and collinearity analysis and outputs the Circos-formated file to plot. Many recent genetic research reports in Nature and Science have applied Circos-styled figures, but still it only provides circular plot.

There are many online platforms for genome evolution that are dedicated to synteny and collinearity analysis. Meanwhile, more and more researchers use their visualization services in their research procedures. Since the cost of calculation grows exponentially with the amount of data, particularly in the process of analysis, most of these platforms provide dotted or linear plot because it is much simpler and faster to accomplish. Examples of such platforms include Plant Genome Duplication DataBase [27], MIPS CrowsNest [28], and Yeast Gene Order Browser [29]. Only very few platforms can generate complex plots, such as circular plot and multialignment plot, for example, the famous Ensembl [30, 31]. In plant comparative genomics, PLAZA 3.0 is one of the most powerful all-in-one solutions in this area. It has collected a large quantity of data and developed the full utility sets to support research from analyses to visualizations [32]. And yet none of them provides full support of vector graphic outputs. The gap for multistyled vector-based plots in synteny and collinearity remains to be filled.

Generally, synteny and collinearity analysis is a frequent task in comparative genomics research. Many analysis software packages are available, but not effective in visualizing the result, shown in Table 1. The problems include lack of graphic visualization, simple representation, or inextensible output format. On the other hand, general-purpose visualization tools are powerful, but not specific for synteny and collinearity display. This requirement grows rapidly while higher throughput of datasets generates higher resolution outputs.

Table 1: Software list for synteny and collinearity visualization.

Software name
Publishing year
Graphical synteny
Visualization types
Vector graphics

MCMuSeC
2009
✕
✕
✕

OrthoCluster
2009
✓
Dual bar, linear
✕

i-ADHoRe 3
2011
✓
Dotted hierarchy

FISH
2003
✕
✕
✕

ColinearScan
2006
✕
✕
✕

MCScan
2008
✕
✕
✕

CYNTENATOR
2010
✕
✕
✕

SyMAP
2011
✓
Dual bar
✕

MCScanX
2012
✓
Dotted linear circular
✕

SynChro
2014
✓
Dual bar
✕

GSV
2011
✓
Dual bar
✕

EasyFig
2011
✓
Dual bar

C-Sibelia
2013
✕
Circos-format
✕

Gbrowse-syn
2010
✓
Dual bar
✕

Kegg
2000
✓
Network
✕

Circos
2009
✓
Circular
✓

WebLogo
2004
✕
Textual
✓

VGSC
2015

Dotted linear circular

In this paper, we introduce VGSC, a purpose-built toolkit in visualizing the synteny and collinearity into general graphical format, including both raster (JPEG, Bitmap, and PNG) and vector graphics (SVG, EPS, and PDF).

Vector graphics are a computational representation of graphical objects using vectors, a geometric object with a magnitude and a direction. In this way, vector graphics are normally combinations of geometrical primitives, such as points, lines, curves, shapes, and polygons. In contrast, raster images use dot matrix data to represent a generally rectangular grid of pixels or points of color. The advantages of vectors are scale-invariance, rotate-invariance, and transform-invariance. They enable the antialiasing feature, which means graphics can be magnified infinitely without loss of quality. Therefore, vector graphics are widely used in scientific research, especially in the bioinformatics research where a massive amount of data from the sequencing process generates various types of high-resolution graphs. A good case in point is WebLogo [33], which is a software package to generate sequence logos, the graphical representations of the patterns within a multiple sequence alignment. WebLogo is so popular that in some areas it becomes the gold standard. This tool is very effective and efficient because it provides both command line interface and web interface, as well as both raster and vector graphics as outputs.

2. Implementation and Result

2.1. Software Architecture

Vector Graphic toolkit of genome Synteny and Collinearity (VGSC) is a new web-based interface for synteny and collinearity representation. Its software architecture is shown in Figure 1, in which the command-line toolkits and web-based service are both illustrated. The workflow of plotting remains as simple as most visualization tools: the end user prepares the required datasets and configures the basic parameters; the software then plots accordingly. Many of these features have simplified the process of drawing, so that researches can focus more on the analysis and interpretation of the data.

Figure 1: System architecture of VGSC.

2.2. Data Input and Configuration

In Figure 2, three inputs from end users are required: (1) synteny and collinearity file, (2) gene annotation file, and (3) control file. And they are explained as follows:

(1)Synteny and collinearity file: VGSC operates on the preprocessed synteny and collinearity data. It is easy to convert results from all the common synteny and collinearity analysis software packages into the required format. The detailed requirement is available in the software manual.(2)Gene annotation file: this GFF3 annotation file (http://www.gmod.org/wiki/GFF3) provides the fundamental map for the plotting, which is widely used in gene assembling software and gene databases.(3)Control file: in this file, the detailed configuration sets the width, length, color, and so forth for the plot.

If end users run VGSC in the command line, these settings serve as inputs as textual parameters. A Java Runtime Environment 1.8 is mandatory, as the software is packaged as a Java executable. For users, synteny and collinearity file and annotation file should be uploaded, and the parameters in the control file can be configured directly in the web form. In addition, we have listed a set of data samples with preconfigured parameters in the “Example” section of the website to help end users carry out tests.

(a)

(b)

(c)

(d)

Figure 2: Four types of synteny and collinearity plot: (a) Circle Plot, (b) Bar Plot, (c) Dot Plot, and (d) Dual Synteny Plot. Chromosomes are labeled in species abbreviation plus chromosome ID. os, Oryza sativa; sb, Sorghum bicolor.

2.3. Output and Result

VGSC provides four different types of plots in six different file formats, with which the synteny and collinearity information can be drawn into circle, bars, dots, and dual synteny. Figure 2 demonstrates the four plots generated by a sample data set of the synteny and collinearity across Rice (Oryza sativa) and Sorghum (Sorghum bicolor) from MCScanX website (http://chibba.pgml.uga.edu/mcscan2).

In the command-line executable, we have implemented a plot manager to integrate all types of plots into one command, which has made the selection much easier. We have introduced a multiple file format adaptor, which enables both raster and vector graphics, so that the output file formats expend to SVG, EPS, PDF, JPEG, and BMP, in addition to the popular PNG format. This automatic configuration mechanism is also applied to all the parameter settings, and the detailed settings list is in the software manual, available at http://bio.njfu.edu.cn:8080/vgsc-web/static/downloads/vgsc-manual.pdf.

One of the most important features of VGSC is its ability to produce vector graphics. As Figure 3 demonstrates, compared with raster graphics (right), vector graphics (left) provide higher compatibility when the image is magnified. This is particularly noticeable when high-throughput datasets are concerned. High-quality images are often a requirement for scientific research reports and papers.

Figure 3: Resolution comparison between vector graphics and raster graphics.

For web users, there is a list of options, where end users can specify the type of plot. A dropdown menu is also available, where end users can choose the output file format. Once the settings are confirmed, results can be downloaded as a separate file when the “Download” button is clicked. In the online service, both vector graphics and raster images are provided.

2.4. Online System

Parallel with command-line toolkit, we have published a web-based system, VGSC online, to provide the plotting service and to improve the experience in plotting. It is available at http://bio.njfu.edu.cn:8080/vgsc-web. The VGSC online uses Java Web Technology and is compatible with most of web containers including Tomcat and jetty. Figure 4 shows the screenshot from the example pages in VGSC online. It lists all types of plots with sample data, providing end users with a visual scaffold. We have also published the command line executable for downloading, along with some sample data and relevant documentation. All these resources are provided free.

Figure 4: Screenshot of VGSC online.

3. Conclusion

While many synteny and collinearity tools have become available in recent years, their visual presentation has not been developed accordingly. For this reason, users often have to write additional programs or redraw the synteny and collinearity output files in order to plot a representative high-quality image. This incompleteness of visualization has reduced the efficiency of existing synteny and collinearity detection pipeline. VGSC has been created to fill this gap. A distinguishing feature of VGSC and its online service is that diverse tools for vector graphics of synteny and collinearity are incorporated, which enables rapid and convenient conversion of synteny and collinearity information into graphical insights. Additional plots for downstream analysis, such as plots for gene family, will be implemented in the coming version of VGSC. VGSC therefore will also be an effective tool for structural changes and evolution analysis, annotation for new genomes, and gene family history research.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Authors’ Contribution

Yiqing Xu and Changwei Bi contributed equally to this work.

Acknowledgments

The authors thank National Basic Research Program of China (973 Program) (2012CB114505), China National Funds for Distinguished Young Scientists (31125008), and National Natural Science Foundation of China (Grant no. 31570662 and Grant no. 31500533). They also acknowledge supports from Key Projects in the National Science & Technology Pillar Program during the Twelfth Five-year Plan Period (no. 2012BAD01B07) and the Fundamental Research Funds for the Central Non-profit Research Institution of CAF (CAFYBB2014QB015). This work is also enabled by the Innovative Research Team Program of the Educational Department of China, the Innovative Research Team Program in Universities of Jiangsu Province, and the PAPD (Priority Academic Program Development) program at Nanjing Forestry University.

References

H. Tang, J. E. Bowers, X. Wang, R. Ming, M. Alam, and A. H. Paterson, “Synteny and collinearity in plant genomes,” Science, vol. 320, no. 5875, pp. 486–488, 2008.

A. Coghlan, E. E. Eichler, S. G. Oliver, A. H. Paterson, and L. Stein, “Chromosome evolution in eukaryotes: a multi-kingdom perspective,” Trends in Genetics, vol. 21, no. 12, pp. 673–682, 2005.

B. Dujon, D. Sherman, G. Fischer et al., “Genome evolution in yeasts,” Nature, vol. 430, no. 6995, pp. 35–44, 2004.

Y. Nakatani, H. Takeda, Y. Kohara, and S. Morishita, “Reconstruction of the vertebrate ancestral genome reveals dynamic genome reorganization in early vertebrates,” Genome Research, vol. 17, no. 9, pp. 1254–1265, 2007.

G. Drillon, A. Carbone, and G. Fischer, “SynChro: a fast and easy tool to reconstruct and visualize synteny blocks along eukaryotic chromosomes,” PLoS ONE, vol. 9, no. 3, Article ID e92621, 2014.

C. Soderlund, M. Bomhoff, and W. M. Nelson, “SyMAP v3.4: a turnkey synteny system with application to plant genomes,” Nucleic Acids Research, vol. 39, no. 10, article e68, 2011.

E. Lyons, B. Pedersen, J. Kane et al., “Finding and comparing syntenic regions among Arabidopsis and the outgroups papaya, poplar, and grape: CoGe with rosids,” Plant Physiology, vol. 148, no. 4, pp. 1772–1781, 2008.

J. Ma, L. Zhang, B. B. Suh et al., “Reconstructing contiguous regions of an ancestral genome,” Genome Research, vol. 16, no. 12, pp. 1557–1565, 2006.

K. Vandepoele, Y. Saeys, C. Simillion, J. Raes, and Y. Van de Peer, “The automatic detection of homologous regions (ADHoRe) and its application to microcolinearity between Arabidopsis and rice,” Genome Research, vol. 12, no. 11, pp. 1792–1801, 2002.

X. Ling, X. He, and D. Xin, “Detecting gene clusters under evolutionary constraint in a large number of genomes,” Bioinformatics, vol. 25, no. 5, pp. 571–577, 2009.

I. A. Vergara and N. Chen, “Using OrthoCluster for the detection of synteny blocks among multiple genomes,” Current Protocols in Bioinformatics, chapter 6: unit 6.10, 2009.

X. Zeng, M. J. Nesbitt, J. Pei, K. Wang, I. A. Vergara, and N. Chen, “OrthoCluster: a new tool for mining synteny blocks and applications in comparative genomics,” in Proceedings of the 11th International Conference on Extending Database Technology: Advances in Database Technology (EDBT '08), pp. 656–667, ACM, Nantes, France, March 2008.

X. Wang, X. Shi, Z. Li et al., “Statistical inference of chromosomal homology based on gene colinearity and applications to Arabidopsis and rice,” BMC Bioinformatics, vol. 7, article 447, 2006.

P. P. Calabrese, S. Chakravarty, and T. J. Vision, “Fast identification and statistical evaluation of segmental homologies in comparative maps,” Bioinformatics, vol. 19, no. 1, pp. i74–i80, 2003.

C. Rödelsperger and C. Dieterich, “CYNTENATOR: progressive gene order alignment of 17 vertebrate genomes,” PLoS ONE, vol. 5, no. 1, Article ID e8861, 2010.

M. Abrouk, F. Murat, C. Pont et al., “Palaeogenomics of plants: synteny-based modelling of extinct ancestors,” Trends in Plant Science, vol. 15, no. 9, pp. 479–487, 2010.

B. Cai, X. Yang, G. A. Tuskan, and Z.-M. Cheng, “MicroSyn: a user friendly tool for detection of microsynteny in a gene family,” BMC Bioinformatics, vol. 12, article 79, 2011.

S. Proost, J. Fostier, D. De Witte et al., “i-ADHoRe 3.0—fast and sensitive detection of genomic homology in extremely large data sets,” Nucleic Acids Research, vol. 40, no. 2, article e11, 2012.

S. J. McKay, I. A. Vergara, and J. E. Stajich, “UNIT 9.12 using the generic synteny browser (GBrowse_syn),” in Current Protocols in Bioinformatics, chapter 9, John Wiley & Sons, Hoboken, NJ, USA, 2010.

L. D. Stein, C. Mungall, S. Shu et al., “The generic genome browser: a building block for a model organism system database,” Genome Research, vol. 12, no. 10, pp. 1599–1610, 2002.

L. D. Stein, “Using GBrowse 2.0 to visualize and share next-generation sequence data,” Briefings in Bioinformatics, vol. 14, no. 2, pp. 162–171, 2013.

K. V. Revanna, C.-C. Chiu, E. Bierschank, and Q. Dong, “GSV: a web-based genome synteny viewer for customized data,” BMC Bioinformatics, vol. 12, no. 1, article 316, 2011.

M. J. Sullivan, N. K. Petty, and S. A. Beatson, “Easyfig: a genome comparison visualizer,” Bioinformatics, vol. 27, no. 7, pp. 1009–1010, 2011.

Y. Wang, H. Tang, J. D. Debarry et al., “MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity,” Nucleic Acids Research, vol. 40, no. 7, article e49, 2012.

M. Krzywinski, J. Schein, I. Birol et al., “Circos: an information aesthetic for comparative genomics,” Genome Research, vol. 19, no. 9, pp. 1639–1645, 2009.

I. Minkin, H. Pham, E. Starostina, N. Vyahhi, and S. Pham, “C-Sibelia: an easy-to-use and highly accurate tool for bacterial genome comparison,” F1000Research, November 2013, http://f1000research.com/articles/2-258/v1

T.-H. Lee, H. Tang, X. Wang, and A. H. Paterson, “PGDD: a database of gene and genome duplication in plants,” Nucleic Acids Research, vol. 41, no. 1, pp. D1152–D1158, 2013.

T. Nussbaumer, M. M. Martis, S. K. Roessner et al., “MIPS PlantsDB: a database framework for comparative plant genome research,” Nucleic Acids Research, vol. 41, no. 1, pp. D1144–D1151, 2013.

K. P. Byrne and K. H. Wolfe, “The Yeast Gene Order Browser: combining curated homology and syntenic context reveals gene fate in polyploid species,” Genome Research, vol. 15, no. 10, pp. 1456–1461, 2005.

P. Flicek, M. R. Amode, D. Barrell et al., “Ensembl 2014,” Nucleic Acids Research, vol. 42, no. 1, pp. D749–D755, 2014.

F. Cunningham, M. R. Amode, D. Barrell et al., “Ensembl 2015,” Nucleic Acids Research, vol. 43, no. 1, pp. D662–D669, 2015.

S. Proost, M. Van Bel, D. Vaneechoutte et al., “PLAZA 3.0: an access point for plant comparative genomics,” Nucleic Acids Research, vol. 43, no. 1, pp. D974–D981, 2015.

G. E. Crooks, G. Hon, J.-M. Chandonia, and S. E. Brenner, “WebLogo: a sequence logo generator,” Genome Research, vol. 14, no. 6, pp. 1188–1190, 2004.