Xylella fastidiosa comparative genomic database is an information resource to explore the annotation, genomic features, and biology of different strains

Xylella fastidiosa comparative genomic database is an information resource to explore the annotation, genomic features, and biology of different strains

Alessandro M. VaraniI,^*; Claudia B. Monteiro-VitorelloI,^*; Luiz G.P. de Almeida^II; Rangel C. Souza^II; Oberdan L. Cunha^I; Wanessa C. Lima^I; Edwin Civerolo^III; Marie-Anne Van Sluys^I; Ana T.R. Vasconcelos^II

^IGenome and Transposable Elements Laboratory, Departamento de Botânica, Instituto de Biociências, Universidade de São Paulo, São Paulo, SP, Brazil

^IILaboratório de Bioinformática, Laboratório Nacional de Computação Científica, Petrópolis, RJ, Brazil

^IIIUnited States Department of Agriculture, Agricultural Research Service, San Joaquin Valley Agricultural Sciences Center, Parlier, CA, USA

^{Send correspondence to}

ABSTRACT

The Xylella fastidiosa comparative genomic database is a scientific resource with the aim to provide a user-friendly interface for accessing high-quality manually curated genomic annotation and comparative sequence analysis, as well as for identifying and mapping prophage-like elements, a marked feature of Xylella genomes. Here we describe a database and tools for exploring the biology of this important plant pathogen. The hallmarks of this database are the high quality genomic annotation, the functional and comparative genomic analysis and the identification and mapping of prophage-like elements. It is available from web site http://www.xylella.lncc.br.

Key words: genome annotation and assembly, comparative genomics, mobile genetic elements. Received: September 15, 2011; Accepted: December 9, 2011.

The quality of bacterial-genome annotation varies. The lack of a direct link between annotation in public databases, and the functional information accumulated over recent years, highlights how the importance of maintaining this up-to-date is becoming a crucial task in the genomics era (Parkhill et al., 2010). Although considerable literature has accumulated on Xylella fastidiosa over the last decade, this information has not been transferred to the annotation files in public databases. Xylella is a phytopathogenic bacterium that causes economically devastating losses in the yields of such crops as grapes, citrus fruits, almonds and other plant species (Van Sluys et al. (, 2002). The 9a5c strain, the causal agent of citrus variegated chlorosis, was the first bacterial plant pathogen to have its genome completely sequenced (Simpson et al. (, 2000). Nowadays, besides the six different genomes published, additional strains are part of ongoing sequencing projects.

Genomic studies have indicated extensive lateral gene transfer (LGT) related to prophage-like regions, which in turn are related to intra-genomic deletions, insertions and rearrangements (Monteiro-Vitorello et al., 2005; da Silva et al. (, 2007). Moreover, the presence of phage particles has also been demonstrated by both electron microscopy (Chen and Civerolo, 2008), and plaque propagation (Summer et al. (, 2010), all of which implying that phages are capable of playing a major role in genomic shaping and differentiation in Xylella strains (de Mello Varani et al. (, 2008).

Analysis of the genomic differences between closely related strains provides, not only a starting point towards understand functional and evolutionary processes, but also clues towards defining what makes one strain more pathogenic and/or aggressive than others. This information would be useful in epidemiological studies, all of which can potentially lead to the development of novel disease management strategies by identifying potential gene targets for mitigating infection and/or disease development.

We hereby report the first comprehensive and specialized up-to-date database comprising all the sequenced genomes of the different Xylella fastidiosa strains. The web-accessible application was developed, by using the SABIA package (System for Automated Bacterial Integrated Annotation), a public-domain software for the automated identification of genome landmarks that uses a user-friendly interface for browsing and retrieving data and information (Almeida et al., 2004).

Xylella fastidiosa strains were recently grouped into subspecies (Schaad et al. (, 2004, 2009), although the current database version follows the original strain identification. The database includes four complete and finished public genomic sequences of strains that cause citrus variegated chlorosis (9a5c), Pierce's disease (Temecula1), almond leaf scorch and Pierce's diseases (M23), and almond leaf scorch disease (M12). In addition, the public draft genomes of the strains associated with oleander leaf scorch (Ann1) and almond leaf scorch (Dixon) were assembled (closed but not finished) into candidate molecules representing the main replicon and plasmids. Additional information for finishing and gap-closures can be found in the supplementary material.

Prophage-like element identification was carried out using the methodology implemented by de Mello Varani et al. ( (2008). Orthologous clusters were identified using the bidirectional best-hit method (Overbeek et al. (, 1999). This database provides access to the latest annotations that can be downloaded in raw datasets, such as flat file and GenBank file format.

The high-quality annotation process was a collaborative effort among annotator specialists. The database can be searched by gene or protein names, as well as other functional annotation terms. The search engine is capable of further refining queries using SQL rules defined by the user.

Nucleotide and amino acid sequences can be searched by BLAST (Altschul et al. (, 1997).

A genome viewer provides a graphical overview of the position of a given selected gene on the chromosome, as well as of neighboring genes. The annotation integrates information on putative gene products, transcription regulatory sequences and ribosome binding sites. InterPro protein signatures, UniProt (Universal Protein Resource) and the NCBI non-redundant protein database were used with the BLAST program for orthology and similarity assignment. Putative protein localization is assigned by PSORT (Nakai and Horton, 1999), and possible membrane transport capacity using the TCDB database. The Enzyme Commission number (EC Number), Gene Ontology terms and COG phylogenetic classification, were used for functional categorization of the putative gene products. KEGG metabolic pathways are also available through tables and in a graphical overview interface, thereby facilitating user visualization and comparison of the complete set of pathways available in each strain.

All identified prophage-like elements and prophage remnants were characterized and annotated as special features in each strain. They are indicated with a special tag after the gene name, i.e. [phage-related protein, xfp3], wh'xfp3' represents the prophage-like element number three of the 9a5c strain. For other strains, the notation is as follows: xpd 1 to 9 for Temecula1, xap 1 to 11 for Dixon, xop 1 to 10 for Ann1, xmp 1 to 9 for M23, and xp 1 to 7 for M12 strains (for details see "Genome sequence alignment and comparative map of prophage regions of the six strains" in Supplementary Materials of the Xylella fastidiosa comparative database, http://www.xylella.lncc.br/supplementary.html).

The comparative interface consists of a precalculated similarity analysis of Xylella predicted genes against thirteen completely sequenced Proteobacteria genomes, by using reciprocal BLAST searches for the computation of BBH clusters (Table1). The core and pan genome calculation was estimated, and can be accessed as tables and graphs (Figure 1). Proteins involved in a common structural complex or metabolic pathway are highlighted and this information is associated with the identification of strain-specific regions that might be related to host specificity.

The database attempts to provide a comprehensive view of all sequence elements and their related functions in Xylella genomes, providing a valuable online resource for Xylella community researchers. Expectedly, its use will contribute to understanding the biology of Xylella, and to the study of the mechanisms involved in its pathogenicity. New sequenced Xylella genomes can be included in future versions of the database, after the complete annotation and curation process.

Acknowledgments

We thank Anamaria Aranha Camargo, Alexandre Morais do Amaral, Sergio Verjovski-Almeida, Luis Eduardo Aranha Camargo, Carlos F.M. Menck, Marilis do Valle Marques, Eliana Macedo Lemos, Manoel Vitor Lemos, Ana Lucia Nascimento, Mariana C. de Oliveira and Marcelo Zerillo for assistance in the re-annotation process, and Roger Paixão for extensive bioinformatics support. This work was supported by FAPESP (São Paulo, SP, Brazil) and CAPES / CNPq (Brasília, DF, Brazil) and USDA-ARS.

Internet Resources

• Almeida LG, Paixão R, Souza RC, Costa GC, Barrientos FJ, Santos MT, Almeida DF and Vasconcelos AT (2004) A System for Automated Bacterial (genome) Integrated Annotation SABIA. Bioinformatics 20:2832-2833.
• Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W and Lipman DJ (1997) Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res 25:3389-3402.
• Chen J and Civerolo EL (2008) Morphological evidence for phages in Xylella fastidiosa. Virol J 6:75.
• da Silva VS, Shida CS, Rodrigues FB, Ribeiro DC, de Souza AA, Coletta-Filho HD, Machado MA, Nunes LR and de Oliveira RC (2007) Comparative genomic characterization of citrus-associated Xylella fastidiosa strains. BMC Genomics 8:e474.
• de Mello Varani A, Souza RC, Nakaya HI, de Lima WC, Paula de Almeida LG, Kitajima EW, Chen J, Civerolo E, Vasconcelos AT and Van Sluys MA (2008) Origins of the Xylella fastidiosa prophage-like regions and their impact in genome differentiation. PLoS One 3:e4059.
• Monteiro-Vitorello CB, de Oliveira MC, Zerillo MM, Varani AM, Civerolo E and Van Sluys MA (2005) Xylella and Xanthomonas Mobil'omics. OMICS 9:146-159.
• Nakai K and Horton P (1999) PSORT: A program for detecting sorting signals in proteins and predicting their subcellular localization. Trends Biochem Sci 24:34-36.
• Overbeek R, Fonstein M, D'Souza M, Pusch GD and Maltsev N (1999) The use of gene clusters to infer functional coupling. Proc Natl Acad Sci USA 96:2896-2901.
• Parkhill J, Birney E and Kersey P (2010) Genomic information infrastructure after the deluge. Genome Biol 11:e402.
• Schaad NW, Postnikova E, Lacy G, Fatmi M and Chang CJ (2004) Xylella fastidiosa subspecies: X. fastidiosa subsp piercei, subsp. nov., X. fastidiosa subsp. multiplex subsp. nov., and X. fastidiosa subsp. pauca subsp. nov. Syst Appl Microbiol 27:290-300. Erratum in Syst Appl Microbiol 27:763.
• Simpson AJ, Reinach FC, Arruda P, Abreu FA, Acencio M, Alvarenga R, Alves LM, Araya JE, Baia GS, Baptista CS, et al. (2000) The genome sequence of the plant pathogen Xylella fastidiosa. Nature 406:151-159.
• Summer EJ, Enderle CJ, Ahern SJ, Gill JJ, Torres CP, Appel DN, Black MC, Young R and Gonzalez CF (2010) Genomic and biological analysis of phage Xfas53 and related prophages of Xylella fastidiosa. J Bacteriol 192:179-190.
• Van Sluys MA, de Oliveira MC, Monteiro-Vitorello CB, Miyaki CY, Furlan LR, Camargo LE, da Silva AC, Moon DH, Takita MA, Lemos EG, et al. ( (2002) Comparative analyses of the complete genome sequences of Pierce's disease and citrus variegated chlorosis strains of Xylella fastidiosa. J Bacteriol 185:1018-1026.
• InterPro protein sequence analysis & classification, http://www.ebi.ac.uk/interpro (August 10, 2011).
• KEGG: Kyoto Encyclopedia of Genes and Genomes, http://www.genome.ad.jp/kegg (August 10, 2011).
• The Gene Ontology, http://www.geneontology.org (August 10, 2011).
• Clusters of Orthologous Groups (COGs), http://www.ncbi.nlm.nih.gov/COG (August 10, 2011).
• Prediction of Protein Sorting Signals and Localization Sites in Amino Acid Sequences (PSORT), http://psort.hgc.jp (August 10, 2011).
• Functional and Phylogenetic Classification of Membrane Transport Proteins (TCDB), http://www.tcdb.org (August 10, 2011).
• Universal Protein Resource (UNIPROT), http://www.uniprot.org (August 10, 2011).

*

Current address: Departamento de Genética, Escola Superior de Agricultura "Luiz de Queiroz", Universidade de São Paulo, 13418-900 Piracicaba,SP, Brazil.

Associate Editor: Luís Carlos de Souza Ferreira

License information: This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Publication Dates

• Publication in this collection
  16 Feb 2012
• Date of issue
  2012

History

• Received
  15 Sept 2011
• Accepted
  09 Dec 2011