- Poster presentation
- Open Access
An integrated database of Eucalyptusspp. genome project
BMC Proceedings volume 5, Article number: P170 (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).
To integrate all data described above, we developed a web site (Fig. 1) hosted in a Fedora Linux machine with MySQL database server. The web interface is based on a combination of CGI scripts using PERL language (including BioPerl module) and the Apache Web Server. The site contains many bioinformatics tools allowing the user perform keyword or local BLAST search in assembled unigenes. Also it is possible to connect these results with gene expression analysis. Moreover, the Gbrowse software (Generic Genome Browser) (Fig. 2) was used to visualization the data in a genomic context, integrating the different information by clickable tracks. The top track is the reference genome assembly and the other tracks correspond to assembled unigenes and RNA-seq data mapped into reference.
Baudet C, Dias Z: New EST Trimming Strategy. Brazilian Symposium on Bioinformatics. 2005, Lecture Notes in Bioinformatics – Berlin – Germany: Springer – Verlag, 3594: 206-209.
Huang X, Madan A: CAP3: A DNA Sequence Assembly Program. Genome Res. 1999, 9: 868-877. 10.1101/gr.9.9.868.
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.
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.
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.
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.
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.
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.
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.
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.
Trapnell C, Pachter L, Salzberg SL: TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009, 25 (9): 1105-1111. 10.1093/bioinformatics/btp120.
Audic S, Claverie JM: The significance of Digital Gene Expression Profiles. Genome Res. 1997, 7: 986-995.
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.
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.