Skip to content

Advertisement

You're viewing the new version of our site. Please leave us feedback.

Learn more

BMC Bioinformatics

Open Access

IGDD: a database of intronless genes in dicots

BMC BioinformaticsBMC series – open, inclusive and trusted201617:289

https://doi.org/10.1186/s12859-016-1148-9

Received: 28 November 2015

Accepted: 13 July 2016

Published: 27 July 2016

Abstract

Background

Intronless genes are a significant characteristic of prokaryotes. Systematic identification and annotation are primary and crucial steps for determining the functions of intronless genes and understanding their occurrence in eukaryotes.

Description

In this paper, we describe the construction of the Intronless Genes Database in Dicots (IGDD; available at http://bio.njfu.edu.cn/igdd/), which contains data for five well-annotated plants including Arabidopsis thaliana, Carica papaya, Populus trichocarpa, Salix suchowensis and Vitis vinifera. Using highly visual settings, IGDD displays the structural and functional annotations, the homolog groups, the syntenic relationships, the expression patterns, and the statistical characteristics of intronless genes. In addition, useful tools such as an advanced search and local BLAST are available through a user-friendly and intuitive web interface.

Conclusion

In conclusion, the IGDD provides a comprehensive and up-to-date platform for researchers to assist the exploration of intronless genes in dicot plants.

Keywords

Intronless genesDatabaseDicots

Background

Genes that are not interrupted by introns are called intronless genes or single-exon genes. Depending on their structure, genes in eukaryotes fall into two categories: intronless genes and intron-containing genes. Intronless genes can serve as beacons in analyses of gene function and evolution. For example, intronless genes as a model, compared with intron-containing homologs, can enable an inverse approach to studying the considerable roles of introns, which are only found in eukaryotes [1]. Furthermore, studies on intronless genes should help to unravel some evolutionary issues including: (i) the major factors that have contributed to the emergence of intronless genes (gene duplications, inheritance from ancient prokaryotes, retroposition or other mechanisms); (ii) the evolutionary significance of retroposition (retrogenes are considered to be intronless); and (iii) the biological origins of introns (i.e., which hypothesis is more correct: the introns-early hypothesis or introns-late hypothesis?) [2]. Therefore, related studies and the construction of databases for intronless genes have received great attention in recent years.

IGD [3] and PIGD [4] are two existing intronless gene databases that contain information about intronless genes in human and Poaceae, respectively. The identification methods and knowledge of intronless genes in these two databases are worthy of reference. However, limitation still persists as they still are not the centralized platforms for intronless genes. Furthermore, IGD is derived from GenBank, which does not contain all the sequences available from genome projects; thus, it does not provide representative intronless gene sets. Compared with IGD, PIGD has an improved methodology, but contains less evolutionary information in its gene annotations.

To build a centralized platform, we developed IGDD—a comprehensive and integrated database for intronless genes that covers five dicots, including Arabidopsis thaliana (Arabidopsis), Carica papaya (papaya), Populus trichocarpa (poplar), Salix suchowensis (willow) and Vitis vinifera (grape). These plants were chosen as representatives not only because of their burgeoning sets of genome sequencing, but also because of their obvious patterns of genome duplication [5]. Specifically, the Arabidopsis thaliana genome has undergone two whole genome duplication events (α and β) within the crucifer lineage and one more ancient genome triplication event (γ) shared with most dicots [6], Carica papaya and Vitis vinifera have each experienced only the γ triplication and no subsequent polyploidies; whereas Populus trichocarpa and Salix suchowensis have undergone a salicoid-specific duplication (ρ) besides the γ triplication [7]. The large number of paralogous genes present in these species offers a unique opportunity to help us elucidate the evolutionary consequences of intronless genes.

IGDD not only contains comprehensive knowledge and links to authoritative external datasets, but also includes the important features within and among species. In addition to basic information on intronless genes, the annotations in IGDD include predicted protein domains, KEGG pathways, GO items, functional descriptions and subcellular localizations. Moreover, IGDD also lists paralogs, orthologs, syntenic blocks, and expression patterns across different developmental stages and tissues. By integrating massive functional and evolutionary information and developing valuable tools, we aim to provide a useful resource and versatile platform that will benefit the related research community.

Construction and content

Data source

The genome sequences of Arabidopsis, papaya, poplar, and grape were acquired from the publicly available database Phytozome [811]. In our database, the latest versions of data in use are TAIR10 for Arabidopsis, ASGPBv0.4 for papaya, v3.0 for poplar, and Genoscope.12X for grape separately. The genome of Salix suchowensis is a most updated in-house database developed by our research group [7].

Identification of intronless genes

Intronless genes were identified by Perl scripts. First, the scripts extracted the information for mRNAs that were not interrupted by “introns” from a GFF file, and genes containing only one continuous exon were considered putative intronless genes. Subsequently, the longest transcript was retained if two or more transcripts represented a gene that met the criteria in the first step. Considering their biological functions, only proteins encoded by intronless genes with lengths ≥ 30aa were selected.

Annotation of intronless genes

To provide informative clues for further functional analysis, we systematically annotated each intronless gene. Initially, we obtained the standard gene information from Phytozome [12], including PFAM [13], GO annotation [14] and PANTHER data [15]. Functional description was then added using the InterproScan program (Version 5.44.0) [16]. Pathway annotations from the KEGG database [17] were also included in IGDD. In addition, we added subcellular localization prediction information from WoLF PSORT [18].

Putative homolog annotation

To predict paralogs of the intronless genes, the method described by Blanc and Wolfe was used [19]. We performed all-against-all nucleotide sequence similarity searches among the transcribed sequences in each genome using the BLASTN software [20]. Sequences with strict cutoffs of alignment length ≥ 300 bp and identity ≥ 40 % were defined as paralogs. We also downloaded the paralogs of intronless genes annotated in Ensembl using Biomart [21]. The results obtained from both methods were included for accurate prediction of paralogs. Orthologs of intronless genes across the five dicots were identified by OrthoMCL with the default parameters [22].

Syntenic analysis

To detect intronless genes that had arisen from whole-genome duplication (WGD), syntenic blocks within and between species were identified using BLAST, OrthoMCL and the MCscanX software [20, 22, 23]. We then retrieved the syntenic blocks containing intronless genes. Finally, to show an overview of all blocks, the intra-genome and inter-genome syntenic relationships of the intronless genes were visualized with Circos [24]. The nonsynonymous (Ka) and synonymous (Ks) substitution rates of gene pairs in each block were determined by KaKs_Calculator 2.0 [25].

Expression patterns of intronless genes

An expression-based resource can provide an important bridge between genotype and phenotype through transcript profiles; therefore, we integrated a wide variety of expression data to explore the spatial and temporal expression patterns of intronless genes during development. Comprehensive microarray data for A. thaliana was acquired from TAIR across 63 tissues [26], for P.trichocarpa from NCBI-GEO (GSE13990) covering nine tissues [27], and for V.vinifera from NCBI-GEO (GSE36128), which included 54 samples representing green and woody tissues [28]. Expression data for S. suchowensis was obtained using RNA-Seq from five tissues [7]. The expression data of the intronless genes was imported into R and Bioconductor, and then the pheatmap package was used to make heatmaps.

Utility and discussion

IGDD collects sequence information from multiple resources to enrich the data for intronless genes. At present, IGDD contains 28,016 intronless genes from five dicots (Table 1). In IGDD, different categories of datasets are classified and compared within and among species with user-friendly graphics. Additionally, analytical tools such as BLAST are embedded in IGDD to help predict the putative orthologous groups of the intronless genes. Thus, IGDD is a solid web platform for searching, browsing, visualizing and downloading intronless genes (Fig. 1a). Below, we discuss the four main functional units in IGDD: (1) the comprehensive individual gene information, (2) the integration of data sources, (3) the BLAST sequence search engine, and (4) the interactive platform.
Table 1

List of five dicots currently served by IGDD

Species

Release version

Identified Intronless Genes

Percentage of Intronless Genes

Arabidopsis thaliana

TAIR 10

5733

20.9 %

Carica papaya

ASGPB V0.4

7791

28.1 %

Populus trichocarpa

V3.0

7209

17.4 %

Salix suchowensis

V 1.0

5657

21.3 %

Vitis vinifera

Genoscope.12X

1626

6.2 %

Fig. 1

An overview of the IGDD website. a The home page. b Comprehensive individual gene information (e.g., AT2G21910). c The browse page. d The comparison page. e The BLAST page

  1. (1)

    Comprehensive individual gene information

    IGDD provides detailed annotations for every putative intronless gene (Fig. 1b). Users can access the webpage in multiple ways, such as by clicking “Browse” on the main navigation bar or importing a gene name into the “Search” section directly. The annotation information for individual genes includes: (i) basic information, (ii) protein sequence features, (iii) homolog groups including paralogs and orthologs, (iv) syntenic blocks, and (v) expression profiles. In detail, basic information consists of the gene identifier, location, strand and annotations for the corresponding coding protein (such as length, isoelectric point (PI), and molecular weight (Mw)). The protein sequence features display protein identities (IDs) for PFAM, GO, PANTHER and KEGG, and functional descriptions. These IDs are hyperlinked to external databases to access more details. Additionally, subcellular localization prediction is shown on the page because it gives an important clue to the protein’s role(s). Gene duplication followed by rapid sequence divergence between paralogous pairs is considered to be the major mechanism for the emergence of new genes [29]. To determine the contribution of gene duplication to the amplification of intronless genes, we provide information on paralogous pairs among intronless genes, and between intronless genes and intron-containing genes. To explore cross-species evolutionary study, ortholog groups between intronless genes were identified. Among the syntenic blocks shown on the webpage, the corresponding intronless gene is highlighted in blue font to facilitate searching. The Ka and Ks values of paralogous pairs are also provided to determine their evolutionary distances. Knowing the expression patterns of intronless genes in different developmental stages and tissues is essential to illustrate whether they have corresponding biological functions. In IGDD, we extracted microarray and RNA-seq data to assess the temporal and spatial expression of intronless genes. When an intronless gene is found to be expressed, the expression data is visualized by a line diagram that allows users to view the expression changes across different tissues and developmental stages.

     
  2. (2)

    Integration of data sources

    Based on the collection of available data as aforementioned, individual genes were further analyzed and compared within and among species. We dedicated two main functions of the interface to this: the “Browse” and “Comparison” sections. Specifically, the comparative data for intronless genes within species was analyzed from a variety of perspectives, and visualized on the “Browse” webpage (Fig. 1c). From the interface, basic information on the classification under different attributes is displayed, such as the distribution of intronless genes on chromosomes, subcellular localization, distribution of protein length, PI and Mw. In particular, users can browse and download detailed data by clicking the corresponding section of the figure to view the differences. Apart from basic information, we also collected functional annotations through GO associations, PFAM domain information, and KEGG data. The “Browse” page shows the proportion of genes and gene product attributes associated with cellular components, biological processes, molecular functions, gene families, and pathways. To study genome organization and evolution, colinearity information within species can be applied to analyze segmental and WGD events. The Circos software [24] was employed to enable browsing of the syntenic relationships between intronless genes and other genes using different colored curves. In addition, to better understand whether intronless genes are associated with plant phenotypes, the gene expression in different tissues is clearly visible from the heat map. The heat map provides not only an overview of the global gene expression trends, but also conclusive evidence that these intronless genes are truly transcribed to mRNA.

    The “Comparison” section compares gene characteristics across species and provides detailed information by clicking on the image (Fig. 1d). We found that proportions of intronless genes against the total number of genes differed significantly across the five investigated species, with 20.9, 28.1, 17.4, 21.3, and 6.2 % in Arabidopsis, papaya, poplar, willow, and grape, respectively (Chi-square test, P < 0.05). The average protein sizes were 279, 182, 229, 333, and 178aa in Arabidopsis, papaya, poplar, willow, and grape, respectively. In pairwise comparisons, the differences in the average protein sizes were remarkable (μ-test; P < 0.05) except between papaya and grape. Conversely, we noticed that intronless genes had similar characteristics across the five dicots. For the protein features, the peak Mw value ranged from 8 to 13 kDa, and the PI distribution showed two main peaks at 6 and 9.

     
  3. (3)

    The BLAST sequence search engine

    The BLAST (Basic Local Alignment Search Tool) interface provides a flexible way to search for homologous genes of every intronless gene stored in IGDD. Users can submit query sequence(s) and adjust the BLAST parameters including database and search program. To make viewing the result page easier, we set the following two parameters: e-value cutoff of 1e-5 and retaining the top 100 hits based on e-values.

     
  4. (4)

    The interactive platform

    To expand the annotation scope of IGDD, reliable functional and evolutionary information was retrieved from the well-known databases PFAM, PANTHER, KEGG, GO and Expasy. Once plant genomics websites release the latest versions of the four dicots genomes, the datasets in IGDD would be updated synchronously. Compared with comparative genomics websites, IGDD focuses on the relationship among intronless genes, which is helpful to explore functional and evolutionary features on intronless genes across dicots. All of the information in IGDD is available for researchers at an FTP (File Transfer Protocol) site. In addition, IGDD encourages users to submit fully annotated intronless genes from other dicots to the database. When the quality of the submitted information meets our requirements, the IGDD curator will import the data into the database. Overall, the interactive platforms, including links to external databases, and download and submission components, make IGDD a comprehensive and systematic platform for the research community.

     

Conclusions

Our main goal with IGDD is to construct a comprehensive platform for the exploration of intronless genes in dicots. IGDD integrates various types of data and links with multiple external databases to provide rich annotation information, which can be browsed and retrieved through user-friendly web interfaces. Biological tools such as BLAST and a comparison platform are also provided to facilitate investigations into the functional and evolutionary consequences of intronless genes. In future, the authors will regularly review authoritative databases for new dicot genomes, and expand the content of IGDD. Moreover, we intend to add new functions and integrate multiple data sources to enhance the IGDD database. For example, expression profiles will focus not only on different tissues and developmental stages, but also experimental treatments, such as biotic and abiotic stress. Our ultimate goal is to construct a co-expression network that provides informative clues to the regulation of intronless gene expression. In conclusion, we hope that IGDD will serve as a useful resource and fundamental platform for studying intronless genes, especially their occurrence and evolution.

Availability and requirements

Database: IGDD

Database homepage: http://bio.njfu.edu.cn/igdd

Operating system(s): Unix

Programming language: C+, JavaScript, Perl

Other requirements: MySQL, Apache, PHP

These data are freely available without restrictions for use by academics. Inquiries concerning the database may be directed to IGDD@163.com.

Abbreviations

IGDD, intronless genes database in dicots

WGD, whole-genome duplication

PI, isoelectric point

Mw, molecular weight

BLAST, basic local alignment search tool

FTP, file transfer protocol

Declarations

Acknowledgments

Special thanks go for the anonymous reviewers for their enlightening comments in reformulating this paper.

Funding

This work was supported by the National Basic Research Project (2012CB114505), the Natural Science Foundation of China (31570662), and Anhui provincial Natural Science Foundation (1608085QC65). It was also enabled by China postdoctoral science foundation (2015 M581806), the Innovative Research Team of the Educational Department of China and the PAPD (Priority Academic Program Development) program at Nanjing Forestry University.

Availability of data and materials

The database is available at http://bio.njfu.edu.cn/igdd/, and no restrictions to its use by non-academics.

Authors’ contribution

HY performed integration, constructed the database platform, and wrote the manuscript. XD, KF and QM helped with the design of database platform and update of the database, and provided scientific suggestions and criticisms for improving the manuscript and website. TY participated in the design, helped write the manuscript and supervised the whole project. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Open AccessThis 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.

Authors’ Affiliations

(1)
Key Laboratory of Forest Genetics & Biotechnology, Nanjing Forestry University
(2)
Laboratory of Modern Biotechnology, Anhui Agricultural University

References

  1. Tine M, Kuhl H, Beck A, Bargelloni L, Reinhardt R. Comparative analysis of intronless genes in teleost fish genomes: insights into their evolution and molecular function. Mar Genomics. 2011;4(2):109–19.View ArticlePubMedGoogle Scholar
  2. Sakharkar M, Kangueane P. Genome SEGE: a database for‘intronless’ genes in eukaryotic genomes. BMC Bioinformatics. 2004;5(1):1.View ArticleGoogle Scholar
  3. Louhichi A, Fourati A, Rebaï A. IGD: a resource for intronless genes in the human genome. Gene. 2011;488(1):35–40.View ArticlePubMedGoogle Scholar
  4. Yan H, Jiang C, Li X, Sheng L, Dong Q, Peng X, Li Q, Zhao Y, Jiang H, Cheng B. PIGD: a database for intronless genes in the Poaceae. BMC Genomics. 2014;15(1):832.View ArticlePubMedPubMed CentralGoogle Scholar
  5. Tang H, Bowers JE, Wang X, Ming R, Alam M, Paterson AH. Synteny and collinearity in plant genomes. Science. 2008;320(5875):486–8.View ArticlePubMedGoogle Scholar
  6. Bowers JE, Chapman BA, Rong J, Paterson AH. Unravelling angiosperm genome evolution by phylogenetic analysis of chromosomal duplication events. Nature. 2003;422(6930):433–8.View ArticlePubMedGoogle Scholar
  7. Dai X, Hu Q, Cai Q, Feng K, Ye N, Tuskan GA, Milne R, Chen Y, Wan Z, Wang Z. The willow genome and divergent evolution from poplar after the common genome duplication. Cell Res. 2014;24:1274–7.View ArticlePubMedPubMed CentralGoogle Scholar
  8. Lamesch P, Berardini TZ, Li D, Swarbreck D, Wilks C, Sasidharan R, Muller R, Dreher K, Alexander DL, Garcia-Hernandez M. The Arabidopsis Information Resource (TAIR): improved gene annotation and new tools. Nucleic Acids Res. 2012;40(D1):D1202–10.View ArticlePubMedGoogle Scholar
  9. Ming R, Hou S, Feng Y, Yu Q, Dionne-Laporte A, Saw JH, Senin P, Wang W, Ly BV, Lewis KL. The draft genome of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus). Nature. 2008;452(7190):991–6.View ArticlePubMedPubMed CentralGoogle Scholar
  10. Tuskan GA, Difazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, Putnam N, Ralph S, Rombauts S, Salamov A. The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science. 2006;313(5793):1596–604.View ArticlePubMedGoogle Scholar
  11. Jaillon O, Aury J-M, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007;449(7161):463–7.View ArticlePubMedGoogle Scholar
  12. Goodstein DM, Shu S, Howson R, Neupane R, Hayes RD, Fazo J, Mitros T, Dirks W, Hellsten U, Putnam N. Phytozome: a comparative platform for green plant genomics. Nucleic Acids Res. 2012;40(D1):D1178–86.View ArticlePubMedGoogle Scholar
  13. Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, Pang N, Forslund K, Ceric G, Clements J. The Pfam protein families database. Nucleic Acids Res. 2012;40(D1):D290–301.View ArticlePubMedGoogle Scholar
  14. Consortium GO. The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res. 2004;32 suppl 1:D258–61.View ArticleGoogle Scholar
  15. Mi H, Muruganujan A, Thomas PD. PANTHER in 2013: modeling the evolution of gene function, and other gene attributes, in the context of phylogenetic trees. Nucleic Acids Res. 2013;41(D1):D377–86.View ArticlePubMedGoogle Scholar
  16. Jones P, Binns D, Chang H-Y, Fraser M, Li W, McAnulla C, McWilliam H, Maslen J, Mitchell A, Nuka G. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30(9):1236–40.View ArticlePubMedPubMed CentralGoogle Scholar
  17. Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M. The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004;32 suppl 1:D277–80.View ArticlePubMedPubMed CentralGoogle Scholar
  18. Horton P, Park K-J, Obayashi T, Fujita N, Harada H, Adams-Collier C, Nakai K. WoLF PSORT: protein localization predictor. Nucleic Acids Res. 2007;35 suppl 2:W585–7.View ArticlePubMedPubMed CentralGoogle Scholar
  19. Blanc G, Wolfe KH. Widespread paleopolyploidy in model plant species inferred from age distributions of duplicate genes. Plant Cell. 2004;16(7):1667–78.View ArticlePubMedPubMed CentralGoogle Scholar
  20. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.View ArticlePubMedPubMed CentralGoogle Scholar
  21. Herrero J, Muffato M, Beal K, Fitzgerald S, Gordon L, Pignatelli M, Vilella AJ, Searle SM, Amode R, Brent S. Ensembl comparative genomics resources. Database. 2016;2016:bav096.View ArticlePubMedPubMed CentralGoogle Scholar
  22. Li L, Stoeckert CJ, Roos DS. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003;13(9):2178–89.View ArticlePubMedPubMed CentralGoogle Scholar
  23. Wang Y, Tang H, DeBarry JD, Tan X, Li J, Wang X, Lee T-H, Jin H, Marler B, Guo H. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012;40(7):e49–9.View ArticlePubMedPubMed CentralGoogle Scholar
  24. Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19(9):1639–45.View ArticlePubMedPubMed CentralGoogle Scholar
  25. Wang D, Zhang Y, Zhang Z, Zhu J, Yu J. KaKs_Calculator 2.0: a toolkit incorporating gamma-series methods and sliding window strategies. Genomics Proteomics Bioinformatics. 2010;8(1):77–80.View ArticlePubMedGoogle Scholar
  26. Schmid M, Davison TS, Henz SR, Pape UJ, Demar M, Vingron M, Schölkopf B, Weigel D, Lohmann JU. A gene expression map of Arabidopsis thaliana development. Nat Genet. 2005;37(5):501–6.View ArticlePubMedGoogle Scholar
  27. Wilkins O, Nahal H, Foong J, Provart NJ, Campbell MM. Expansion and diversification of the Populus R2R3-MYB family of transcription factors. Plant Physiol. 2009;149(2):981–93.View ArticlePubMedPubMed CentralGoogle Scholar
  28. Fasoli M, Dal Santo S, Zenoni S, Tornielli GB, Farina L, Zamboni A, Porceddu A, Venturini L, Bicego M, Murino V. The grapevine expression atlas reveals a deep transcriptome shift driving the entire plant into a maturation program. Plant Cell. 2012;24(9):3489–505.View ArticlePubMedPubMed CentralGoogle Scholar
  29. Zhang J. Evolution by gene duplication: an update. Trends Ecol Evol. 2003;18(6):292–8.View ArticleGoogle Scholar

Copyright

© The Author(s). 2016

Advertisement