An integrated database of Eucalyptusspp. genome project
© Nascimento et al; licensee BioMed Central Ltd. 2011
Published: 13 September 2011
The species of the genus Eucalyptus are the most planted for the fiber crop in the world. They are mainly utilized for timber, pulp and paper production. Brazil, helped by the favorable weather conditions, appears as a big producer and exporter of eucalyptus derivates. In 2002, the Brazilian network research of the Eucalyptus Genome (Genolyptus) was established with the goal of integrating several academic and private institutions currently working with eucalyptus genomics in Brazil. This project generated around 200.000 ESTs from several tissues and conditions. Consequently, several individual projects have been implemented generating other transcriptome databases, in special, using RNA-Seq technology. In 2010, a draft genome (http://eucalyptusdb.bi.up.ac.za) of the specie E. grandis was produced by researches of the Joint Genome Institute (DOE-JGI) and the Eucalyptus Genome Network (EUCAGEN). The main goal of this work is to develop an Eucalyptusdatabase (http://www.lge.ibi.unicamp.br/genolyptus) integrating public and private data in a friendly and secure web interface with bioinformatics tools that allowing the users perform complex searches.
Results and discussion
First, the public and private ESTs (130,290 from Genolyptus and 36,981 from NCBI) were assembled producing 48,760 unigenes (17,795 contigs and 30,765 singlets). Basically, the bdtrimmer  and CAP3  programs were used to perform sequence trimming (exclude vector, ribosomal, low quality and too short reads) and sequence assembly, respectively.
The autofact pipeline  was used to perform an automatic annotation of the assembled unigenes based on BLAST  searches, e-value cutoff of 1e-5, against some protein databases, including: non-redundant (NR) database of NCBI, uniref90 and uniref100 – databases containing only curated proteins , pfam – database of proteins families , kegg – database of metabolic pathways  and Gene Ontology (GO) – database of functional annotation .
The Genomic and Expression Laboratory at State University of Campinas (http://www.lge.ibi.unicamp.br) sequenced ten RNA-Seq libraries from four species (E. Urograndis, E. globulus, E. grandis and E. urophylla) using the Illumina/Solexa technology. Additionally, three RNA-seq libraries  were downloaded from NCBI (SRA – sequence read archive). All RNA-seq reads were aligned against the assembled unigenes and genome assembly using the SOAP2  and TopHat  aligners, configured to allow up two mismatches, discard sequences with “N”s and return all optimal alignments.
In order to perform a differential expression analysis between ESTs or RNA-seq libraries some normalization pipelines and statistical tests have been implemented.From ESTs, the differentially expressed genes between libraries were performed applying AC test  in assembled unigenes. The results are available to the users by a web interface (called Electronic Northern) that allows searches by gene or library name. Furthermore, it is possible to compare the gene expression between two or more libraries. From the RNA-seq libraries, the DEG-seq software  was used to perform normalization and statistical analysis considering 99% of confidence rate (cut-off of 0.01).
The authors would like to acknowledge all researches of the Joint Genome Institute (DOE-JGI) and the Eucalyptus Genome Network (EUCAGEN), responsible to produce the draft genome of the E. grandis. Moreover, we thank the CNPQ (Conselho Nacional de Desenvolvimento Científico e Tecnológico - Brazil) for the financial support of this work.
- Baudet C, Dias Z: New EST Trimming Strategy. Brazilian Symposium on Bioinformatics. 2005, Lecture Notes in Bioinformatics – Berlin – Germany: Springer – Verlag, 3594: 206-209.Google Scholar
- Huang X, Madan A: CAP3: A DNA Sequence Assembly Program. Genome Res. 1999, 9: 868-877. 10.1101/gr.9.9.868.PubMed CentralView ArticlePubMedGoogle Scholar
- Koski LB, Gray LW, Lang BF, Burger G: AutoFACT: An Automatic Functional Annotation and Classification Tool. BMC Bioinformatics. 2005, 6: 151-10.1186/1471-2105-6-151.PubMed CentralView ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucl. Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- Suzek BE, Huang H, McGarvey P, Mazumber R, Wu CH: Uniref: comprehensive and non-redundant UniProt reference clusters. Bioinformatics. 2007, 23 (10): 1282-1288. 10.1093/bioinformatics/btm098.View ArticlePubMedGoogle Scholar
- Bateman A, Birney E, Cerruti L, Durbin R, Etwiller L, Griffiths-Jones S, Howe KL, Marshall M, Sonnhammer ELL: The Pfam Protein Families Database. Nucl. Acids Res. 2002, 30 (1): 276-280. 10.1093/nar/30.1.276.PubMed CentralView ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S: KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucl. Acids Res. 2000, 28 (1): 27-30. 10.1093/nar/28.1.27.PubMed CentralView ArticlePubMedGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene Ontology: tool for the unification of biology. Nature Genetics. 2000, 2525-29.Google Scholar
- Mizrachi E, Hefer CA, Ranik M, Joubert F, Myburg AA: De novoassembled expressed gene catalogue f a fast-growing Eucalyptus tree produced by Illumina mRNA-Seq. BMC Genomics. 2010, 11 (681): 1471-2164.Google Scholar
- Li R, Yu C, Li Y, Lam T-W, Yiu S-M, Kristiansen K, Wang J: SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009, 25 (15): 1966-1967. 10.1093/bioinformatics/btp336.View ArticlePubMedGoogle Scholar
- Trapnell C, Pachter L, Salzberg SL: TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009, 25 (9): 1105-1111. 10.1093/bioinformatics/btp120.PubMed CentralView ArticlePubMedGoogle Scholar
- Audic S, Claverie JM: The significance of Digital Gene Expression Profiles. Genome Res. 1997, 7: 986-995.PubMedGoogle Scholar
- Wang L, Feng Z, Wang X, Wang X, Zhang X: DEGSeq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010, 26 (1): 136-138. 10.1093/bioinformatics/btp612.View 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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.