Search

Metagenome-assembled genomes (MAGs) are microbial genomes reconstructed from metagenome data. In the last few years, many thousands of MAGs have been reported in the literature, for a variety of environments and host-associated microbiota, including humans. MAGs have helped us better understand microbial populations and their interactions with the environment where they live; moreover most MAGs belong to novel species, therefore helping to decrease the so-called microbial dark matter. However, questions about the biological reality of MAGs have not, in general, been properly addressed. In this review, I define the notions of hypothetical MAGs and conserved hypothetical MAGs. These notions should help with the understanding of the biological reality of MAGs, their worldwide occurrence, and the efforts to improve MAG recovery processes.

The experimental determination of a bacteriophage host is a laborious procedure. For this reason, there is a pressing need for reliable computational predictions of bacteriophage hosts in phage research in general and in phage therapy in particular. Here, we present a new program called vHULK for phage host prediction based on 9,504 phage genome features. These features take into account alignment significance scores between predicted-protein sequences in the phage genomes and a curated database of viral protein families. The features were fed to a deep neural network, and four distinct models were trained to predict 61 different host genera and 52 host species. In random controlled test sets, the program obtained 99% and 98% accuracy values at the genus and species levels, respectively. On a validation dataset with 2,178 phage genomes, mean accuracies were 82% and 52% at the genus and species levels, respectively. When compared against other phage host prediction programs on the same validation dataset, vHULK achieved substantially better performance, therefore demonstrating that the program is an advance on the state-of-art in phage host prediction. vHULK is freely available at https://github.com/LaboratorioBioinformatica/vHULK.

The Gram-negative bacterium Xylella fastidiosa colonizes plant xylem vessels and is obligately vectored by xylem sap-feeding hemipteran insects. X. fastidiosa causes diseases in many plant species but in a variety of its plant hosts this bacterium behaves as a commensal endophyte. Originally confined to the Americas, infecting mainly grapevine, citrus and coffee plants, X. fastidiosa has spread to several plant species in Europe, causing devastating crop diseases. Although many pathogenicity and virulence factors have been identified in X. fastidiosa which enable the bacterium to successfully establish in the xylem tissue, the mechanisms by which distinct X. fastidiosa strains colonize and cause disease in specific plant hosts have not been fully elucidated. Here we present comparative analyses of 94 publicly available whole-genome sequences of X. fastidiosa strains with the goal of providing insights into plant host specificity determinants for this phytopathogen as well as of expanding the knowledge of its mobile genetic elements (MGE) content, mainly prophages. Our results revealed a pangenome of 4,549 protein coding sequences (CDSs) which is still open. The core- and accessory genomes comprise 954 and 2,219 CDSs, respectively. Phylogenetic tree construction using all core genome CDSs grouped the strains in three major clades of subspecies fastidiosa, multiplex and pauca, with subclades related to the strains’ sequence type (ST) obtained from multi-locus sequence typing (MLST). The geographic region where the strains were collected showed stronger association with the clades of X. fastidiosa strains rather than the plant species from which they were isolated. Among the CDS related to virulence and pathogenicity found in the core genome, those related to lipopolysaccharide (LPS) synthesis and trimeric autotransporter adhesins (TAA) are somewhat related with the plant host of a given strain according to phylogenetic inference. The X. fastidiosa accessory genome is represented by an abundant and heterogeneous mobilome, which includes a diversity of prophage regions. In summary, the genome comparisons reported here will enable a better understanding of the diversity of phylogenetically close genomes and warrant further investigation of LPS and TAAs as potential X. fastidiosa host-specificity determinants.Impact statementThe bacterium Xylella fastidiosa is a pathogen that infects many plant species and has caused devastating diseases in grapevine, citrus, coffee, and olive plants. This phytopathogen X. fastidiosa is original from the Americas and has emerged in Europe where it is causing severe economic losses for olive producers, mainly in Italy. Although many pathogenicity and virulence factors have been identified in X. fastidiosa, which enable this bacterium to successfully establish in the xylem vessels network, the mechanisms by which distinct X. fastidiosa strains colonize and cause disease in the different plant host species have not been fully elucidated. The comparative analyses of 94 whole-genome sequences from X. fastidiosa strains from diverse hosts and geographic regions provide insights into host specificity determinants for this phytopathogen as well as expand the knowledge of its mobile genetic elements (MGE) content, mainly prophages. Our results contribute for a better understanding of the diversity of phylogenetically close genomes and warrant further experimental investigation of lipopolysaccharide and trimeric autotransporter adhesins as potential host-specificity determinants for X. fastidiosa.Data summaryAll genomic sequences were accessed from publicly available GenBank RefSeq database at NCBI (National Center for Biotechnology Information). A full listing of NCBI accession numbers for X. fastidiosa strains described in this paper is available in Table S1 (available in the online version of this article).

Background Composting is an important technique for environment-friendly degradation of organic material, and is a microbe-driven process. Previous metagenomic studies of composting have presented a general description of the taxonomic and functional diversity of its microbial populations, but they have lacked more specific information on the key organisms that are active during the process. Results Here we present and analyze 60 mostly high-quality metagenome-assembled genomes (MAGs) recovered from time-series samples of two thermophilic composting cells, of which 47 are potentially new bacterial species; 24 of those did not have any hits in two public MAG datasets at the 95% average nucleotide identity level. Analyses of gene content and expressed functions based on metatranscriptome data for one of the cells grouped the MAGs in three clusters along the 99-day composting process. By applying metabolic modeling methods, we were able to predict metabolic dependencies between MAGs. These models indicate the importance of coadjuvant bacteria that do not carry out lignocellulose degradation but may contribute to the management of reactive oxygen species and with enzymes that increase bioenergetic efficiency in composting, such as hydrogenases and N2O reductase. Strong metabolic dependencies predicted between MAGs revealed key interactions relying on exchange of H+, NH3, O2 and CO2, as well as glucose, glutamate, succinate, fumarate and others, highlighting the importance of functional stratification and syntrophic interactions during biomass conversion. Our model includes 22 out of 49 MAGs recovered from one composting cell data. Based on this model we highlight that Rhodothermus marinus, Thermobispora bispora and a novel Gammaproteobacterium are dominant players in chemolithotrophic metabolism and cross-feeding interactions. Conclusions The results obtained expand our knowledge of the taxonomic and functional diversity of composting bacteria and provide a model of their dynamic metabolic interactions.

Serratia liquefaciens strain FG3 (SlFG3), isolated from the flower of Stachytarpheta glabra in the Brazilian ferruginous fields, has distinctive genomic, adaptive, and biotechnological potential. Herein, using a combination of genomics and molecular approaches, we unlocked the evolution of the adaptive traits acquired by S1FG3, which exhibits the second largest chromosome containing the largest conjugative plasmids described for Serratia. Comparative analysis revealed the presence of 18 genomic islands and 311 unique protein families involved in distinct adaptive features. S1FG3 has a diversified repertoire of genes associated with Nonribosomal peptides (NRPs/PKS), a complete and functional cluster related to cellulose synthesis, and an extensive and functional repertoire of oxidative metabolism genes. In addition, S1FG3 possesses a complete pathway related to protocatecuate and chloroaromatic degradation, and a complete repertoire of genes related to DNA repair and protection that includes mechanisms related to UV light tolerance, redox process resistance, and a laterally acquired capacity to protect DNA using phosphorothioation. These findings summarize that SlFG3 is well-adapted to different biotic and abiotic stress situations imposed by extreme conditions associated with ferruginous fields, unlocking the impact of the lateral gene transfer to adjust the genome for extreme environments, and providing insight into the evolution of prokaryotes. Coordenacao de Aperfeicoamento de Pessoal de Nivel Superior -Brazil (CAPES)CAPES [001, CFP 51/2013, 3385/2013]; National Council of Technological and Scientific Development (CNPq)National Council for Scientific and Technological Development (CNPq) [481226/2013-3]; Foundation of Protection to Research of the State of Minas Gerais - FAPEMIG [APQ-02387-14, APQ-02357-17]; CNPqNational Council for Scientific and Technological Development (CNPq); Fundect-MS [TO 141/2016, TO 007/2015, 007/2015 SIAFEM 025139]; UFOP grants; Laboratory of Biochemistry and Molecular Biology of the Faculty of Agrarian and Veterinary Sciences UNESP, Jaboticabal campus Thanks to all members of the Laboratory of Biochemistry and Molecular Biology (LBBM, Federal University of Ouro Preto, UFOP) and to the Laboratory of Biochemistry and Molecular Biology of the Faculty of Agrarian and Veterinary Sciences UNESP, Jaboticabal campus, for their support. This study was financed in part by the Coordenacao de Aperfeicoamento de Pessoal de Nivel Superior -Brazil (CAPES) - Finance Code 001 - (the BIGA Project, CFP 51/2013, process 3385/2013), National Council of Technological and Scientific Development (CNPq Process 481226/2013-3), and Foundation of Protection to Research of the State of Minas Gerais - FAPEMIG (process APQ-02387-14 and process APQ-02357-17). LMM, JCS, AMV, NFA and JAF have a research fellowship from CNPq. NFA has grants from Fundect-MS (TO 141/2016 and TO 007/2015) and Fundect-MS (007/2015 SIAFEM 025139). LAGM and LMM have UFOP grants. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Assessing the compatibility between gene family phylogenies is a crucial and often computationally demanding step in many phylogenomic analyses. Here, we describe the Evolutionary Similarity Index (IES), a means to assess shared evolution between gene families using a weighted orthogonal distance regression model applied to sequence distances. The utilization of pairwise distance matrices circumvents comparisons between gene tree topologies, which are inherently uncertain and sensitive to evolutionary model choice, phylogenetic reconstruction artifacts, and other sources of error. Furthermore, IES enables the many-to-many pairing of multiple copies between similarly evolving gene families. This is done by selecting non-overlapping pairs of copies, one from each assessed family, and yielding the least sum of squared residuals. Analyses of simulated gene family data sets show that IES’s accuracy is on par with popular tree-based methods while also less susceptible to noise introduced by sequence alignment and evolutionary model fitting. Applying IES to an empirical data set of 1,322 genes from 42 archaeal genomes identified eight major clusters of gene families with compatible evolutionary trends. The most cohesive cluster consisted of 62 genes with compatible evolutionary signal, which occur as both single-copy and multiple homologs per genome; phylogenetic analysis of concatenated alignments from this cluster produced a tree closely matching previously published species trees for Archaea. Four other clusters are mainly composed of accessory genes with limited distribution among Archaea and enriched toward specific metabolic functions. Pairwise evolutionary distances obtained from these accessory gene clusters suggest patterns of interphyla horizontal gene transfer. An IES implementation is available at https://github.com/lthiberiol/evolSimIndex.

Summary Large-scale mapping of antigens and epitopes is pivotal for developing immunotherapies but challenging, especially for eukaryotic pathogens, owing to their large genomes. Here, we developed an integrated platform for genome phage display (gPhage) to show that unbiased libraries of the eukaryotic parasite Trypanosoma cruzi enable the identification of thousands of antigens recognized by serum samples from patients with Chagas disease. Because most of these antigens are hypothetical proteins, gPhage provides evidence of their expression during infection. We built and validated a comprehensive map of Chagas disease antibody response to show how linear and putative conformation epitopes, many rich in repetitive elements, allow the parasite to evade a buildup of neutralizing antibodies directed against protein domains that mediate infection pathogenesis. Thus, the gPhage platform is a reproducible and effective tool for rapid simultaneous identification of epitopes and antigens, not only in Chagas disease but perhaps also in globally emerging/reemerging acute pathogens., Graphical abstract, Highlights • Genomic shotgun phage display (gPhage) of eukaryotes is feasible and promising. • gPhage allows rapid antigen ID and epitope mapping, including 3D structures. • Conformation epitopes can be identified and validated by using the gPhage platform. • Most Chagas disease antigens are hypothetical proteins rich in repetitive elements., Parasitology; Sequence analysis; Systems biology

Assessing the phylogenetic compatibility between individual gene families is a crucial and often computationally demanding step in many phylogenomics analyses. Here we describe the Evolutionary Similarity Index (IES) to assess shared evolution between gene families using a weighted Orthogonal Distance Regression applied to sequence distances. This approach allows for straightforward pairing of paralogs between co-evolving gene families without resorting to multiple tests, or a priori assumptions of molecular interactions between protein products from assessed genes. The utilization of pairwise distance matrices, while less informative than phylogenetic trees, circumvents error-prone comparisons between trees whose topologies are inherently uncertain. Analyses of simulated gene family evolution datasets showed that IES was more accurate and less susceptible to noise than popular tree-based methods (Robinson-Foulds and geodesic distance) for assessing evolutionary signal compatibility, since it bypasses phylogenetic reconstruction and its inherent uncertainty. Applying IES to a real dataset of 1,322 genes from 42 archaeal genomes identified eight major clusters of gene families with compatible evolutionary trends. Four of these clusters included genes with a taxonomic distribution across all archaeal phyla, while other clusters included a subset of taxa that do not map to generally accepted archaeal clades, indicating possible shared horizontal transfers by clustered gene families. We identify one strongly connected set of 62 genes from the same cluster, occurring as both single-copy and multiple homologs per genome, with compatible phylogenetic reconstructions closely matching previously published species trees for Archaea. An IES implementation is available at https://github.com/lthiberiol/evolSimIndex.

Background: Thermophilic composting is a semi-engineered process carried out by diverse microbial communities. Composting is an environment friendly way of degrading biomass; its study may help uncover important biomass-degrading organisms and key enzymes. DNA sequence-based previous studies have presented a general description of the microbial-molecular features of composting, but they have lacked more specific information on the key organisms that are active during the process and their genomes. Methods: We present an analysis of metagenome-assembled genomes (MAGs) obtained from time-series samples of a thermophilic composting process in the São Paulo Zoological Park (Brazil). Our results are based on a careful analysis of MAG gene content and on metabolic modeling of their interactions. Results: We recovered 60 MAGs from sequencing datasets from two separate composting cells. Phylogenetic analysis shows that 47 of these MAGs represent novel taxa at the genus or higher levels. We have analyzed the gene repertoire of these MAGs in terms of lignocellulose degradation, secondary metabolite production, antibiotic resistance genes, denitrification genes, sulfur metabolism, hydrogen metabolism, and oxygen metabolism. For one of the composting cells we also had metatranscriptome data, which allowed a deeper analysis of 49 MAGs. This analysis showed the presence of three distinct clusters of MAGs with varying activity during the 99-day composting process. The interaction model pointed to Sphaerobacter thermophilus and Thermobispora bispora as key players in the process, as well as other bacteria that are novel. Our results also show the importance of coadjuvant bacteria and of microbial functions related to efficient bioenergetic processes during biomass conversion, such as N2O reduction and hydrogenases. A novel acidobacteria MAG encodes N2O reductase hallmark genes (nosZD). Strong metabolic dependencies predicted between MAGs revealed that cross-feeding in composting can be determined by complementary functions found in the genomes of producers and consumers, supporting the Black Queen hypothesis for co-evolutionary interactions. Conclusions: This study reveals for the first time the key bacterial players in thermophilic composting and provides a model of their dynamic metabolic interactions. These findings pave the way for more rational composting procedures and provide information that could help the development of novel biomass-degrading technologies.

Summary This protocol describes the genomic phage (gPhage) display platform, a large-scale antigen and epitope mapping technique. We constructed a gPhage display peptide library of a eukaryotic organism, Trypanosoma cruzi (causative agent of Chagas disease), to map the antibody response landscape against the parasite. Here, we used an organism with a relatively large but intronless genome, although future applications could include other prevalent or (re)emerging infectious organisms carrying genomes with a limited number of introns. For complete details on the use and execution of this protocol, please refer to Teixeira et al. (2021)., Graphical abstract, Highlights • A streamlined protocol to build genome shotgun phage display libraries • Scalable to different organisms, eukaryotes or prokaryotes • Optimizes key steps for high throughput identification of antibody antigens/epitopes • gPhage is also a platform to study conformational and difficult-to-validate antigens, This protocol describes the genomic phage (gPhage) display platform, a large-scale antigen and epitope mapping technique. We constructed a gPhage display peptide library of a eukaryotic organism, Trypanosoma cruzi (causative agent of Chagas disease), to map the antibody response landscape against the parasite. Here, we used an organism with a relatively large but intronless genome, although future applications could include other prevalent or (re)emerging infectious organisms carrying genomes with a limited number of introns.

Made available in DSpace on 2021-06-25T10:44:05Z (GMT). No. of bitstreams: 0 Previous issue date: 2020-11-01 Here we present and analyze the complete genome of Alcaligenes faecalis strain Mc250 (Mc250), a bacterium isolated from the roots of Mimosa calodendron, an endemic plant growing in ferruginous rupestrian grasslands in Minas Gerais State, Brazil. The genome has 4,159,911 bp and 3,719 predicted protein-coding genes, in a single chromosome. Comparison of the Mc250 genome with 36 other Alcaligenes faecalis genomes revealed that there is considerable gene content variation among these strains, with the core genome representing only 39% of the protein-coding gene repertoire of Mc250. Mc250 encodes a complete denitrification pathway, a network of pathways associated with phenolic compounds degradation, and genes associated with HCN and siderophores synthesis; we also found a repertoire of genes associated with metal internalization and metabolism, sulfate/sulfonate and cysteine metabolism, oxidative stress and DNA repair. These findings reveal the genomic basis for the adaptation of this bacterium to the harsh environmental conditions from where it was isolated. Gene clusters associated with ectoine, terpene, resorcinol, and emulsan biosynthesis that can confer some competitive advantage were also found. Experimental results showed that Mc250 was able to reduce (∼60%) the virulence phenotype of the plant pathogen Xanthomonas citri subsp. citri when co-inoculated in Citrus sinensis, and was able to eradicate 98% of juveniles and stabilize the hatching rate of eggs to 4% in two species of agricultural nematodes. These results reveal biotechnological potential for the Mc250 strain and warrant its further investigation as a biocontrol and plant growth-promoting bacterium. Núcleo de Pesquisas em Ciências Biológicas (NUPEB) Universidade Federal de Ouro Preto Departamento de Ciências Biológicas (DECBI) Instituto de Ciências Exatas e Biológicas (ICEB) Universidade Federal de Ouro Preto (UFOP) Instituto Pristino Departamento de Bioquimica (DBQ) Instituto de Quimica (IQ) Universidade de São Paulo (USP) Instituto de Biotecnologia Aplicada a Agropecuária (BIOAGRO) Universidade Federal de Vicosa (UFV) Departamento de Tecnologia Faculdade de Ciências Agrárias e Veterinárias de Jaboticabal (FCAV) Universidade Estadual Paulista (UNESP) Faculdade de Computação (FACOM) Universidade Federal de Mato Grosso do Sul Departamento de Tecnologia Faculdade de Ciências Agrárias e Veterinárias de Jaboticabal (FCAV) Universidade Estadual Paulista (UNESP)

Mycobacterium bovis is an important animal pathogen worldwide, parasitizing wild and domesticated vertebrate live stocks as well as humans. We investigated evolutionary patterns and processes in 38 sequenced genomes of M. bovis from the Americas, Europe, Asia, and Africa. We found that US strains are evolving distinctly (though with low magnitude) from most South American lineages regarding core-coding size, GC-content, and codon usage bias, and we further employ agreement to this pattern as proxy for phylogenetic reliability. Three data types were used for tree inference: 1. Single copy core-coding genes (core-coding); 2. Indel blocks coded as numeric characters (indels); 3. Gene presence/absence similarly coded (pangenome). Core-coding´s performance was below expected regarded we used 2,000+ genes, but indels and pangenome individually estimated reasonable trees, with HGT rates as high as 200x the DNA substitution rate. In addressing why core-coding performed poorly, we excluded phylogenetic method bias, homologous recombination, positive, and purifying selection as possible causes, whilst transient retention of random mildly deleterious mutations due to recent diversification, and also ongoing pseudogeneization, may be plausible. Although total evidence´s retention index was as good as presence/absence data individually in identifying US samples, it additionaly found a South American clade with more lineages excluding only two Brazilian samples. We further inferred divergence times using indels + pangenome due to lack of informativeness of core-coding, recovering relatively similar date estimates to previous studies on the origin of M. bovis (< 10,000 y) regardless of the dating prior (strict vs. uncorrelated clock, uniform vs. exponential clock, different age/rate standard deviation priors), evincing the suitability and unsuspected robustness of using strictly genomic-wide presence/absence data to dating, while also decisively rejecting a strict clock for our dataset. All trees supported a relatively basal position of African strains, which were isolated from Homo sapiens, suggesting it was an important region for early diversification and confirming we were one of the earliest hosts of the species. Evolution of US strains was affected by lateral acquisition of 57 exclusive genes from, including genes important for virulence. A clade of three US localities, distant from one another, is under even faster evolutionary rates supported by core-coding, indels, and pangenome. Fil: Patané, José S.L.. Universidade de Sao Paulo; Brasil Fil: Martins, Joaquim. Universidade de Sao Paulo; Brasil Fil: Castelão, Ana Beatriz. Universidade Federal do Mato Grosso do Sul; Brasil Fil: Nishibe, Christiane. Universidade Federal do Mato Grosso do Sul; Brasil Fil: Montera, Luciana. Universidade Federal do Mato Grosso do Sul; Brasil Fil: Bigi, Fabiana. Instituto Nacional de Tecnología Agropecuaria. Centro Nacional de Investigaciones Agropecuarias; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina Fil: Zumárraga, Martín José. Instituto Nacional de Tecnología Agropecuaria. Centro Nacional de Investigaciones Agropecuarias; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina Fil: Cataldi, Ángel Adrián. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina. Instituto Nacional de Tecnología Agropecuaria. Centro Nacional de Investigaciones Agropecuarias; Argentina Fil: Junior, Antônio Fonseca. Ministerio da Agricultura Pecuaria e Abastecimento de Brasil. Empresa Brasileira de Pesquisa Agropecuaria; Brasil Fil: Roxo, Eliana. Instituto Biológico de São Pablo; Brasil Fil: Osório, Ana Luiza A. R.. Universidade Federal do Mato Grosso do Sul; Brasil Fil: Jorge, Klaudia S.. Universidade Federal do Mato Grosso do Sul; Brasil Fil: Thacker, Tyler C.. United States Department of Agriculture; Estados Unidos Fil: Almeida, Nalvo F.. Universidade Federal do Mato Grosso do Sul; Brasil Fil: Araújo, Flabio R.. Ministerio da Agricultura Pecuaria e Abastecimento de Brasil. Empresa Brasileira de Pesquisa Agropecuaria; Brasil Fil: Setubal, João C.. Universidade de Sao Paulo; Brasil. Biocomplexity Institute of Virginia Tech; Estados Unidos

Despite some notable progress in data sharing policies and practices, restrictions are still often placed on the open and unconditional use of various genomic data after they have received official approval for release to the public domain or to public databases. These restrictions, which often conflict with the terms and conditions of the funding bodies who supported the release of those data for the benefit of the scientific community and society, are perpetuated by the lack of clear guiding rules for data usage. Existing guidelines for data released to the public domain recognize but fail to resolve tensions between the importance of free and unconditional use of these data and the “right” of the data producers to the first publication. This self-contradiction has resulted in a loophole that allows different interpretations and a continuous debate between data producers and data users on the use of public data. We argue that the publicly available data should be treated as open data, a shared resource with unrestricted use for analysis, interpretation, and publication.

Made available in DSpace on 2020-12-12T01:43:24Z (GMT). No. of bitstreams: 0 Previous issue date: 2019-10-11 Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) Xanthomonas citri pv. aurantifolii pathotype B (XauB) and pathotype C (XauC) are the causative agents respectively of citrus canker B and C, diseases of citrus plants related to the better-known citrus canker A, caused by Xanthomonas citri pv. citri. The study of the genomes of strains of these related bacterial species has the potential to bring new understanding to the molecular basis of citrus canker as well as their evolutionary history. Up to now only one genome sequence of XauB and only one genome sequence of XauC have been available, both in draft status. Here we present two new genome sequences of XauB (both complete) and five new genome sequences of XauC (two complete). A phylogenomic analysis of these seven genome sequences along with 24 other related Xanthomonas genomes showed that there are two distinct and well-supported major clades, the XauB and XauC clade and the Xanthomonas citri pv. citri clade. An analysis of 62 Type III Secretion System effector genes showed that there are 42 effectors with variable presence/absence or pseudogene status among the 31 genomes analyzed. A comparative analysis of secretion-system and surface-structure genes showed that the XauB and XauC genomes lack several key genes in pathogenicity-related subsystems. These subsystems, the Types I and IV Secretion Systems, and the Type IV pilus, therefore emerge as important ones in helping explain the aggressiveness of the A type of citrus canker and the apparent dominance in the field of the corresponding strain over the B and C strains. Programa de Pós-graduação em Biotecnologia Núcleo de Pesquisas em Ciências Biológicas Universidade Federal de Ouro Preto Laboratório Especial de Ciclo Celular Instituto Butantan Departamento de Tecnologia Universidade Estadual Paulista UNESP, Campus de Jaboticabal Departamento de Fitopatologia e Nematologia Escola Superior de Agricultura Luiz de Queiroz Universidade de São Paulo Departamento de Bioquímica Instituto de Química Universidade de São Paulo Fundo de Defesa da Citricultura (FUNDECITRUS) Faculdade de Computação Universidade Federal de Mato Grosso do Sul Departamento de Ciências Biológicas Instituto de Ciências Exatas e Biológicas Universidade Federal de Ouro Preto Departamento de Tecnologia Universidade Estadual Paulista UNESP, Campus de Jaboticabal

Made available in DSpace on 2020-12-12T01:43:17Z (GMT). No. of bitstreams: 0 Previous issue date: 2019-01-01 Fundo de Defesa da Citricultura Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) Fundação de Amparo à Pesquisa do Estado de Minas Gerais (FAPEMIG) Background. In Citrus cultures, three species of Xanthomonas are known to cause distinct diseases. X. citri subsp. citri patothype A, X. fuscans subsp. aurantifolii pathotypes B and C, and X. alfalfae subsp. citrumelonis, are the causative agents of cancrosis A, B, C, and citrus bacterial spots, respectively. Although these species exhibit different levels of virulence and aggressiveness, only limited alternatives are currently available for proper and early detection of these diseases in the fields. The present study aimed to develop a new molecular diagnostic method based on genomic sequences derived from the four species of Xanthomonas. Results. Using comparative genomics approaches, primers were synthesized for the identification of the four causative agents of citrus diseases. These primers were validated for their specificity to their target DNA by both conventional and multiplex PCR. Upon evaluation, their sensitivity was found to be 0.02 ng/µl in vitro and 1.5 × 104 CFU ml−1 in infected leaves. Additionally, none of the primers were able to generate amplicons in 19 other genomes of Xanthomonas not associated with Citrus and one species of Xylella, the causal agent of citrus variegated chlorosis (CVC). This denotes strong specificity of the primers for the different species of Xanthomonas investigated in this study. Conclusions. We demonstrated that these markers can be used as potential candidates for performing in vivo molecular diagnosis exclusively for citrus-associated Xanthomonas. The bioinformatics pipeline developed in this study to design specific genomic regions is capable of generating specific primers. It is freely available and can be utilized for any other model organism. Núcleo de Pesquisas em Ciências Biológicas Universidade Federal de Ouro Preto Departamento de Tecnologia Faculdade de Ciências Agrárias e Veterinárias de Jaboticabal Universidade Estadual Paulista-Unesp Departamento de Fitopatologia e Nematologia Escola Superior de Agricultura ‘‘Luiz de Queiroz’’ Departamento de Bioquímica Instituto de Química Universidade de São Paulo Instituto de Computação Universidade Estadual de Campinas Faculdade de Computação Universidade Federal de Mato Grosso do Sul Departamento de Ciências Biológicas Universidade Federal de Ouro Preto Departamento de Tecnologia Faculdade de Ciências Agrárias e Veterinárias de Jaboticabal Universidade Estadual Paulista-Unesp Fundo de Defesa da Citricultura: 007/2015 SIAFEM 025139 CAPES: 3385/2013 CNPq: 481226/2013-3 FAPEMIG: APQ-02387-14 CAPES: CFP 51/2013

The genetic determinants of bacterial pathogenicity are highly variable between species and strains. However, a factor that is commonly associated with virulent Gram-negative bacteria, including manyAeromonasspp., is the type 3 secretion system (T3SS), which is used to inject effector proteins into target eukaryotic cells. In this study, we developed a bioinformatics pipeline to identify T3SS effector proteins, applied this approach to the genomes of 105Aeromonasstrains isolated from environmental, mutualistic, or pathogenic contexts and evaluated the cytotoxicity of the identified effectors through their heterologous expression in yeast. The developed pipeline uses a two-step approach, where candidate families are initially selected using HMM profiles with minimal similarity scores against the Virulence Factors DataBase (VFDB), followed by strict comparisons against positive and negative control datasets, greatly reducing the number of false positives. Using our approach, we identified 21AeromonasT3SS likely effector groups, of which 8 represented known or characterized effectors, while the remaining 13 had not previously been described inAeromonas. We experimentally validated ourin silicofindings by assessing the cytotoxicity of representative effectors inSaccharomyces cerevisiaeBY4741, with 15 out of 21 assayed proteins eliciting a cytotoxic effect in yeast. The results of this study demonstrate the utility of our approach, combining a novelin silicosearch method within vivoexperimental validation, and will be useful in future research aimed at identifying and authenticating bacterial effector proteins from other genera.

Motivation: With the recent advances in DNA sequencing technologies, the study of the genetic composition of living organisms has become more accessible for researchers. Several advances have been achieved because of it, especially in the health sciences. However, many challenges which emerge from the complexity of sequencing projects remain unsolved. Among them is the task of assembling DNA fragments from previously unsequenced organisms, which is classified as an NP-hard (nondeterministic polynomial time hard) problem, for which no efficient computational solution with reasonable execution time exists. However, several tools that produce approximate solutions have been used with results that have facilitated scientific discoveries, although there is ample room for improvement. As with other NP-hard problems, machine learning algorithms have been one of the approaches used in recent years in an attempt to find better solutions to the DNA fragment assembly problem, although still at a low scale. Results: This paper presents a broad review of pioneering literature comprising artificial intelligence-based DNA assemblers—particularly the ones that use machine learning—to provide an overview of state-of-the-art approaches and to serve as a starting point for further study in this field.

Here we present MARVEL, a tool for prediction of double-stranded DNA bacteriophage sequences in metagenomic bins. MARVEL uses a random forest machine learning approach. We trained the program on a dataset with 1,247 phage and 1,029 bacterial genomes, and tested it on a dataset with 335 bacterial and 177 phage genomes. We show that three simple genomic features extracted from contig sequences were sufficient to achieve a good performance in separating bacterial from phage sequences: gene density, strand shifts, and fraction of significant hits to a viral protein database. We compared the performance of MARVEL to that of VirSorter and VirFinder, two popular programs for predicting viral sequences. Our results show that all three programs have comparable specificity, but MARVEL achieves much better performance on the recall (sensitivity) measure. This means that MARVEL should be able to identify many more phage sequences in metagenomic bins than heretofore has been possible. In a simple test with real data, containing mostly bacterial sequences, MARVEL classified 58 out of 209 bins as phage genomes; other evidence suggests that 57 of these 58 bins are novel phage sequences. Sao Paulo State Research Foundation (FAPESP)Fundacao de Amparo a Pesquisa do Estado de Sao Paulo (FAPESP) [2014/16450-8]; Brazilian Federal Agency CAPESCAPES; CNPqNational Council for Scientific and Technological Development (CNPq); FAPESPFundacao de Amparo a Pesquisa do Estado de Sao Paulo (FAPESP) [2011/50870-6]; CAPESCAPES [3385/2013] DA was supported in part by fellowship Grant No. 2014/16450-8 from the Sao Paulo State Research Foundation (FAPESP). DA and LB were supported by a fellowship from Brazilian Federal Agency CAPES. JS and AdS wish to acknowledge their respective research fellowships from CNPq. This work was supported in part by FAPESP Grant No. 2011/50870-6 and by CAPES Grant No. 3385/2013.

This chapter covers the theory and practice of ortholog gene set computation. In the theoretical part we give detailed and formal descriptions of the relevant concepts. We also cover the topic of graph-based clustering as a tool to compute ortholog gene sets. In the second part we provide an overview of practical considerations intended for researchers who need to determine orthologous genes from a collection of annotated genomes, briefly describing some of the most popular programs and resources currently available for this task.

This chapter covers the theory and practice of ortholog gene set computation. In the theoretical part we give detailed and formal descriptions of the relevant concepts. We also cover the topic of graph-based clustering as a tool to compute ortholog gene sets. In the second part we provide an overview of practical considerations intended for researchers who need to determine orthologous genes from a collection of annotated genomes, briefly describing some of the most popular programs and resources currently available for this task.

MotivationWhole-genome alignment and pan-genome analysis are useful tools in understanding the similarities and differences of many genomes in an evolutionary context. Here we introduce the concept of pan-synteny graphs, an analysis method that combines elements of both to represent conservation and change of multiple prokaryotic genomes at an architectural level. Pan-synteny graphs represent a reference free approach for the comparison of many genomes and allows for the identification of synteny, insertion, deletion, replacement, inversion, recombination, missed assembly joins, evolutionary hotspots, and reference based scaffolding.ResultsWe present an algorithm for creating whole genome multiple sequence comparisons and a model for representing the similarities and differences among sequences as a graph of syntenic gene families. As part of the pan-synteny graph creation, we first create a de Bruijn graph. Instead of the alphabet of nucleotides commonly used in genome assembly, we use an alphabet of gene families. This de Bruijn graph is then processed to create the pan-synteny graph. Our approach is novel in that it explicitly controls how regions from the same sequence and genome are aligned and generates a graph in which all sequences are fully represented as paths. This method harnesses previous computation involved in protein family calculation to speed up the creation of whole genome alignment for many genomes. We provide the software suite Panaconda, for the calculation of pan-synteny graphs given annotation input, and an implementation of methods for their layout and visualization.AvailabilityPanaconda is available at https://github.com/aswarren/pangenome_graphs and datasets used in examples are available at https://github.com/aswarren/pangenome_examplesContactAndrew Warren anwarren@vt.edu

Microbial consortia selected from complex lignocellulolytic microbial communities are promising alternatives to deconstruct plant waste, since synergistic action of different enzymes is required for full degradation of plant biomass in biorefining applications. Culture enrichment also facilitates the study of interactions among consortium members, and can be a good source of novel microbial species. Here, we used a sample from a plant waste composting operation in the Sao Paulo Zoo (Brazil) as inoculum to obtain a thermophilic aerobic consortium enriched through multiple passages at 60°C in carboxymethylcellulose as sole carbon source. The microbial community composition of this consortium was investigated by shotgun metagenomics and genome-centric analysis. Six near-complete (over 90%) genomes were reconstructed. Similarity and phylogenetic analyses show that four of these six genomes are novel, with the following hypothesized identifications: a new Thermobacillus species; the first Bacillus thermozeamaize genome (for which currently only 16S sequences are available) or else the first representative of a new family in the Bacillales order; the first representative of a new genus in the Paenibacillaceae family; and the first representative of a new deep-branching family in the Clostridia class. The reconstructed genomes from known species were identified as Geobacillus thermoglucosidasius and Caldibacillus debilis. The metabolic potential of these recovered genomes based on COG and CAZy analyses show that these genomes encode several glycoside hydrolases (GHs) as well as other genes related to lignocellulose breakdown. The new Thermobacillus species stands out for being the richest in diversity and abundance of GHs, possessing the greatest potential for biomass degradation among the six recovered genomes. We also investigated the presence and activity of the organisms corresponding to these genomes in the composting operation from which the consortium was built, using compost metagenome and metatranscriptome datasets generated in a previous study. We obtained strong evidence that five of the six recovered genomes are indeed present and active in that composting process. We have thus discovered three (perhaps four) new thermophillic bacterial species that add to the increasing repertoire of known lignocellulose degraders, whose biotechnological potential can now be investigated in further studies.

Ralstonia solanacearum is an important soil-borne plant pathogen with broad geographical distribution and the ability to cause wilt disease in many agriculturally important crops. Genome sequencing of multiple R. solanacearum strains has identified both unique and shared genetic traits influencing their evolution and ability to colonize plant hosts. Previous research has shown that DNA methylation can drive speciation and modulate virulence in bacteria, but the impact of epigenetic modifications on the diversification and pathogenesis of R. solanacearum is unknown. Sequencing of R. solanacearum strains GMI1000 and UY031 using Single Molecule Real-Time technology allowed us to perform a comparative analysis of R. solanacearum methylomes. Our analysis identified a novel methylation motif associated with a DNA methylase that is conserved in all complete Ralstonia spp. genomes and across the Burkholderiaceae, as well as a methylation motif associated to a phage-borne methylase unique to R. solanacearum UY031. Comparative analysis of the conserved methylation motif revealed that it is most prevalent in gene promoter regions, where it displays a high degree of conservation detectable through phylogenetic footprinting. Analysis of hyper- and hypo-methylated loci identified several genes involved in global and virulence regulatory functions whose expression may be modulated by DNA methylation. Analysis of genome-wide modification patterns identified a significant correlation between DNA modification and transposase genes in R. solanacearum UY031, driven by the presence of a high copy number of ISrso3 insertion sequences in this genome and pointing to a novel mechanism for regulation of transposition. These results set a firm foundation for experimental investigations into the role of DNA methylation in R. solanacearum evolution and its adaptation to different plants.

Background Among viruses, bacteriophages are a group of special interest due to their capacity of infecting bacteria that are important for biotechnology and human health. Composting is a microbial-driven process in which complex organic matter is converted into humus-like substances. In thermophilic composting, the degradation activity is carried out primarily by bacteria and little is known about the presence and role of bacteriophages in this process. Results Using Pseudomonas aeruginosa as host, we isolated three new phages from a composting operation at the Sao Paulo Zoo Park (Brazil). One of the isolated phages is similar to Pseudomonas phage Ab18 and belongs to the Siphoviridae YuA-like viral genus. The other two isolated phages are similar to each other and present genomes sharing low similarity with phage genomes in public databases; we therefore hypothesize that they belong to a new genus in the Podoviridae family. Detailed genomic descriptions and comparisons of the three phages are presented, as well as two new clusters of phage genomes in the Viral Orthologous Clusters database of large DNA viruses. We found sequences encoding homing endonucleases that disrupt a putative ribonucleotide reductase gene and an RNA polymerase subunit 2 gene in two of the phages. These findings provide insights about the evolution of two-subunits RNA polymerases and the possible role of homing endonucleases in this process. Infection tests on 30 different strains of bacteria reveal a narrow host range for the three phages, restricted to P. aeruginosa PA14 and three other P. aeruginosa clinical isolates. Biofilm dissolution assays suggest that these phages could be promising antimicrobial agents against P. aeruginosa PA14 infections. Analyses on composting metagenomic and metatranscriptomic data indicate association between abundance variations in both phage and host populations in the environment. Conclusion The results about the newly discovered and described phages contribute to the understanding of tailed bacteriophage diversity, evolution, and role in the complex composting environment. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3729-z) contains supplementary material, which is available to authorized users.

Motivation Information about metabolic pathways in a comparative context is one of the most powerful tool to help the understanding of genome-based differences in phenotypes among organisms. Although several platforms exist that provide a wealth of information on metabolic pathways of diverse organisms, the comparison among organisms using metabolic pathways is still a difficult task. Results We present TabPath (Tables for Metabolic Pathway), a web-based tool to facilitate comparison of metabolic pathways in genomes based on KEGG. From a selection of pathways and genomes of interest on the menu, TabPath generates user-friendly tables that facilitate analysis of variations in metabolism among the selected organisms. Availability and implementation TabPath is available at http://200.239.132.160:8686.

The Xanthomonadaceae family consists of species of non-pathogenic and pathogenic γ-proteobacteria that infect different hosts, including humans and plants. In this study, we performed a comparative analysis using 69 fully sequenced genomes belonging to this family, with a focus on identifying proteins enriched in phytopathogens that could explain the lifestyle and the ability to infect plants. Using a computational approach, we identified seven phytopathogen-enriched protein families putatively secreted by type II secretory system: PheA (CM-sec), LipA/LesA, VirK, and four families involved in N-glycan degradation, NixE, NixF, NixL, and FucA1. In silico and phylogenetic analyses of these protein families revealed they all have orthologs in other phytopathogenic or symbiotic bacteria, and are involved in the modulation and evasion of the immune system. As a proof of concept, we performed a biochemical characterization of LipA from Xac306 and verified that the mutant strain lost most of its lipase and esterase activities and displayed reduced virulence in citrus. Since this study includes closely related organisms with distinct lifestyles and highlights proteins directly related to adaptation inside plant tissues, novel approaches might use these proteins as biotechnological targets for disease control, and contribute to our understanding of the coevolution of plant-associated bacteria.

Viruses represent the most abundant life forms on the planet. Recent experimental and computational improvements have led to a dramatic increase in the number of viral genome sequences identified primarily from metagenomic samples. As a result of the expanding catalog of metagenomic viral sequences, there exists a need for a comprehensive computational platform integrating all these sequences with associated metadata and analytical tools. Here we present IMG/VR (https://img.jgi.doe.gov/vr/), the largest publicly available database of 3908 isolate reference DNA viruses with 264 413 computationally identified viral contigs from >6000 ecologically diverse metagenomic samples. Approximately half of the viral contigs are grouped into genetically distinct quasi-species clusters. Microbial hosts are predicted for 20 000 viral sequences, revealing nine microbial phyla previously unreported to be infected by viruses. Viral sequences can be queried using a variety of associated metadata, including habitat type and geographic location of the samples, or taxonomic classification according to hallmark viral genes. IMG/VR has a user-friendly interface that allows users to interrogate all integrated data and interact by comparing with external sequences, thus serving as an essential resource in the viral genomics community.

Ancestral genome reconstruction can be understood as a phylogenetic study with more details than a traditional phylogenetic tree reconstruction. We present a new computational system called REGEN for ancestral bacterial genome reconstruction at both the gene and replicon levels. REGEN reconstructs gene content, contiguous gene runs, and replicon structure for each ancestral genome. Along each branch of the phylogenetic tree, REGEN infers evolutionary events, including gene creation and deletion and replicon fission and fusion. The reconstruction can be performed by either a maximum parsimony or a maximum likelihood method. Gene content reconstruction is based on the concept of neighboring gene pairs. REGEN was designed to be used with any set of genomes that are sufficiently related, which will usually be the case for bacteria within the same taxonomic order. We evaluated REGEN using simulated genomes and genomes in the Rhizobiales order. Published version

Diverse gene products including phytotoxins, pathogen-associated molecular patterns, and type III secreted effectors influence interactions between Pseudomonas syringae strains and plants, with additional yet uncharacterized factors likely contributing as well. Of particular interest are those interactions governing pathogen-host specificity. Comparative genomics of closely related pathogens with different host specificity represents an excellent approach for identification of genes contributing to host-range determination. A draft genome sequence of Pseudomonas syringae pv. tomato T1, which is pathogenic on tomato but nonpathogenic on Arabidopsis thaliana, was obtained for this purpose and compared with the genome of the closely related A. thaliana and tomato model pathogen P. syringae pv. tomato DC3000. Although the overall genetic content of each of the two genomes appears to be highly similar, the repertoire of effectors was found to diverge significantly. Several P. syringae pv. tomato T1 effectors absent from strain DC3000 were confirmed to be translocated into plants, with the well-studied effector AvrRpt2 representing a likely candidate for host-range determination. However, the presence of avrRpt2 was not found sufficient to explain A. thaliana resistance to P. syringae pv. tomato T1, suggesting that other effectors and possibly type III secretion system–independent factors also play a role in this interaction.

Ralstonia solanacearum is the causative agent of bacterial wilt of potato. Ralstonia solanacearum strain UY031 belongs to the American phylotype IIB, sequevar 1, also classified as race 3 biovar 2. Here we report the completely sequenced genome of this strain, the first complete genome for phylotype IIB, sequevar 1, and the fourth for the R. solanacearum species complex. In addition to standard genome annotation, we have carried out a curated annotation of type III effector genes, an important pathogenicity-related class of genes for this organism. We identified 60 effector genes, and observed that this effector repertoire is distinct when compared to those from other phylotype IIB strains. Eleven of the effectors appear to be nonfunctional due to disruptive mutations. We also report a methylome analysis of this genome, the first for a R. solanacearum strain. This analysis helped us note the presence of a toxin gene within a region of probable phage origin, raising the hypothesis that this gene may play a role in this strain’s virulence., We also thank COST action Sustain from the European Union for funding and Nemo Peeters and Stéphane Genin for hosting MP for a short stay to carry out UY031 effector annotation. RGS has a Ph.D. fellowship from FAPESP, Brazil. JCS has an investigator fellowship from the Conselho Nacional de Desenvolvimento Cientifico e Tecnologico, Brazil.

We present a bacterial genome computational analysis pipeline, called GenVar. The pipeline, based on the program GeneWise, is designed to analyze an annotated genome and automatically identify missed gene calls and sequence variants such as genes with disrupted reading frames (split genes) and those with insertions and deletions (indels). For a given genome to be analyzed, GenVar relies on a database containing closely related genomes (such as other species or strains) as well as a few additional reference genomes. GenVar also helps identify gene disruptions probably caused by sequencing errors. We exemplify GenVar's capabilities by presenting results from the analysis of four Brucella genomes. Brucella is an important human pathogen and zoonotic agent. The analysis revealed hundreds of missed gene calls, new split genes and indels, several of which are species specific and hence provide valuable clues to the understanding of the genome basis of Brucella pathogenicity and host specificity. PHS HHS [HHSN266200400035C]

Background Schistosomiasis is one of the most prevalent parasitic diseases worldwide and is a public health problem. Schistosoma mansoni is the most widespread species responsible for schistosomiasis in the Americas, Middle East and Africa. Adult female worms (mated to males) release eggs in the hepatic portal vasculature and are the principal cause of morbidity. Comparative separate transcriptomes of female and male adult worms were previously assessed with using microarrays and Serial Analysis of Gene Expression (SAGE), thus limiting the possibility of finding novel genes. Moreover, the egg transcriptome was analyzed only once with limited bacterially cloned cDNA libraries. Methodology/Principal findings To compare the gene expression of S. mansoni eggs, females, and males, we performed RNA-Seq on these three parasite forms using 454/Roche technology and reconstructed the transcriptome using Trinity de novo assembly. The resulting contigs were mapped to the genome and were cross-referenced with predicted Smp genes and H3K4me3 ChIP-Seq public data. For the first time, we obtained separate, unbiased gene expression profiles for S. mansoni eggs and female and male adult worms, identifying enriched biological processes and specific enriched functions for each of the three parasite forms. Transcripts with no match to predicted genes were analyzed for their protein-coding potential and the presence of an encoded conserved protein domain. A set of 232 novel protein-coding genes with putative functions related to reproduction, metabolism, and cell biogenesis was detected, which contributes to the understanding of parasite biology. Conclusions/Significance Large-scale RNA-Seq analysis using de novo assembly associated with genome-wide information for histone marks in the vicinity of gene models constitutes a new approach to transcriptome analysis that has not yet been explored in schistosomes. Importantly, all data have been consolidated into a UCSC Genome Browser search- and download-tool (http://schistosoma.usp.br/). This database provides new ways to explore the schistosome genome and transcriptome and will facilitate molecular research on this important parasite., Author Summary Schistosomiasis is a public health problem caused by parasites of the genus Schistosoma, of which S. mansoni is the primary causative agent. The parasite has a complex life cycle; their sexual reproductive stage is dependent on female and male adult worms mating inside the mesenteric circulation of the human host, with the female releasing hundreds of eggs daily. This phase of the life cycle is responsible for the development of pathology, which is proportional to the number of eggs accumulating in the liver and intestine of the human host. Genome and transcriptome sequencing of this parasite represent important advances in schistosome research, but there is still a need for integrated analyses to better understand the biology of this parasite. In this study, we describe the first large-scale transcriptomes of eggs, and female and male adult worms, the parasite forms that are mainly responsible for the pathology of schistosomiasis. We were able to cross-reference the gene transcription regions with promoter regions, thus improving the gene annotations. Moreover, we identified the expression of novel protein-coding genes not yet described in the current genome annotation, advancing the biological knowledge regarding this parasite.

Computational analyses of four bacterial genomes of the Xanthomonadaceae family reveal new unique genes that may be involved in adaptation, pathogenicity, and host specificity. The Xanthomonas genus presents 3636 unique genes distributed in 1470 families, while Xylella genus presents 1026 unique genes distributed in 375 families. Among Xanthomonas-specific genes, we highlight a large number of cell wall degrading enzymes, proteases, and iron receptors, a set of energy metabolism genes, second copy of the type II secretion system, type III secretion system, flagella and chemotactic machinery, and the xanthomonadin synthesis gene cluster. Important genes unique to the Xylella genus are an additional copy of a type IV pili gene cluster and the complete machinery of colicin V synthesis and secretion. Intersections of gene sets from both genera reveal a cluster of genes homologous to Salmonella's SPI-7 island in Xanthomonas axonopodis pv citri and Xylella fastidiosa 9a5c, which might be involved in host specificity. Each genome also presents important unique genes, such as an HMS cluster, the kdgT gene, and O-antigen in Xanthomonas axonopodis pv citri; a number of avrBS genes and a distinct O-antigen in Xanthomonas campestris pv campestris, a type I restriction-modification system and a nickase gene in Xylella fastidiosa 9a5c, and a type II restriction-modification system and two genes related to peptidoglycan biosynthesis in Xylella fastidiosa temecula 1. All these differences imply a considerable number of gene gains and losses during the divergence of the four lineages, and are associated with structural genome modifications that may have a direct relation with the mode of transmission, adaptation to specific environments and pathogenicity of each organism.

Submitted by Guilherme Lemeszenski (guilherme@nead.unesp.br) on 2014-02-26T17:15:29Z No. of bitstreams: 1 WOS000180834300035.pdf: 1288581 bytes, checksum: c80b0ab1625eeba1cd15e9da6e465b26 (MD5) Made available in DSpace on 2014-02-26T17:15:29Z (GMT). No. of bitstreams: 1 WOS000180834300035.pdf: 1288581 bytes, checksum: c80b0ab1625eeba1cd15e9da6e465b26 (MD5) Previous issue date: 2003-02-01 Submitted by Vitor Silverio Rodrigues (vitorsrodrigues@reitoria.unesp.br) on 2014-05-20T15:25:38Z No. of bitstreams: 1 WOS000180834300035.pdf: 1288581 bytes, checksum: c80b0ab1625eeba1cd15e9da6e465b26 (MD5) Made available in DSpace on 2014-05-20T15:25:38Z (GMT). No. of bitstreams: 1 WOS000180834300035.pdf: 1288581 bytes, checksum: c80b0ab1625eeba1cd15e9da6e465b26 (MD5) Previous issue date: 2003-02-01 Xylella fastidiosa is a xylem-dwelling, insect-transmitted, gamma-proteobacterium that causes diseases in many plants, including grapevine, citrus, periwinkle, almond, oleander, and coffee. X. fastidiosa has an unusually broad host range, has an extensive geographical distribution throughout the American continent, and induces diverse disease phenotypes. Previous molecular analyses indicated three distinct groups of X.fastidiosa isolates that were expected to be genetically divergent. Here we report the genome sequence of X. fastidiosa (Temecula strain), isolated from a naturally infected grapevine with Pierce's disease (PD) in a wine-grape-growing region of California. Comparative analyses with a previously sequenced X.fastidiosa strain responsible for citrus variegated chlorosis (CVC) revealed that 98% of the PD X.fastidiosa Temecula genes are shared with the CVC X. fastidiosa strain 9a5c genes. Furthermore, the average amino acid identity of the open reading frames in the strains is 95.7%. Genomic differences are limited to phage-associated chromosomal rearrangements and deletions that also account for the strain-specific genes present in each genome. Genomic islands, one in each genome, were identified, and their presence in other X.fastidiosa strains was analyzed. We conclude that these two organisms have identical metabolic functions and are likely to use a common set of genes in plant colonization and pathogenesis, permitting convergence of functional genomic strategies. Univ São Paulo, Inst Biociencias, Dept Bot, BR-05508900 São Paulo, Brazil Univ São Paulo, Inst Quim, BR-05508900 São Paulo, Brazil Univ São Paulo, Inst Ciências Biomed, BR-05508900 São Paulo, Brazil Univ São Paulo, Escola Super Agr Luiz de Queiroz, BR-13418900 Piracicaba, Brazil Univ São Paulo, Ctr Energia Nucl Agr, BR-13400970 Piracicaba, Brazil Inst Agron Campinas, Ctr Citricultura Sylio Moreira, BR-13490970 Cordeiropolis, Brazil Univ Estadual Paulista, Fac Ciências Agr & Vet, BR-14884900 Jaboticabal, Brazil Univ Estadual Campinas, Ctr Biol Mol & Engn Genet, BR-13083970 Campinas, SP, Brazil Univ Estadual Campinas, Inst Biol, BR-13083970 Campinas, SP, Brazil Univ Estadual Campinas, Inst Comp, BR-13083970 Campinas, SP, Brazil Univ São Paulo, Fac Ciências Farmaceut, BR-14040901 Ribeirao Preto, Brazil Univ São Paulo, Fac Filosofia Ciências & Letras, BR-14040901 Ribeirao Preto, Brazil Inst Ludwig Pesquisa Canc, BR-01509010 São Paulo, Brazil Univ Mogi das Cruzes, Nucleo Integrado Biotecnol, BR-08780911 Mogi das Cruzes, Brazil Inst Agron Campinas, BR-13001970 Campinas, SP, Brazil Inst Biol, BR-04014002 São Paulo, Brazil Univ Estadual Paulista, Fac Med Vet & Zootecnia, BR-18603970 Botucatu, SP, Brazil Univ Estadual Paulista, Fac Ciências Agron, BR-18603970 Botucatu, SP, Brazil Univ Estadual Paulista, Inst Biociencias, BR-18618000 Botucatu, SP, Brazil Univ Fed Sao Carlos, BR-13600970 São Paulo, Brazil Univ Fed Mato Grosso do Sul, Dept Comp & Estatist, BR-79070900 Campo Grande, MS, Brazil USDA ARS, Parlier, CA 93648 USA Univ Estadual Paulista, Fac Ciências Agr & Vet, BR-14884900 Jaboticabal, Brazil Univ Estadual Paulista, Fac Med Vet & Zootecnia, BR-18603970 Botucatu, SP, Brazil Univ Estadual Paulista, Fac Ciências Agron, BR-18603970 Botucatu, SP, Brazil Univ Estadual Paulista, Inst Biociencias, BR-18618000 Botucatu, SP, Brazil

Understanding the evolutionary history and potential of bacterial pathogens is critical to prevent the emergence of new infectious bacterial diseases. Xanthomonas axonopodis subsp. citri (Xac) (synonym X. citri subsp. citri), which causes citrus canker, is one of the hardest-fought plant bacterial pathogens in US history. Here, we sequenced 21 Xac strains (14 XacA, 3 XacA* and 4 XacA(w)) with different host ranges from North America and Asia and conducted comparative genomic and evolutionary analyses. Our analyses suggest that acquisition of beneficial genes and loss of detrimental genes most likely allowed XacA to infect a broader range of hosts as compared with XacA(w) and XacA*. Recombination was found to have occurred frequently on the relative ancient branches, but rarely on the young branches of the clonal genealogy. The ratio of recombination/mutation ρ/θ was 0.0790±0.0005, implying that the Xac population was clonal in structure. Positive selection has affected 14% (395 out of 2822) of core genes of the citrus canker-causing Xanthomonas. The genes affected are enriched in 'carbohydrate transport and metabolism' and 'DNA replication, recombination and repair' genes (P0.05). Many genes related to virulence, especially genes involved in the type III secretion system and effectors, are affected by positive selection, further highlighting the contribution of positive selection to the evolution of citrus canker-causing Xanthomonas. Our results suggest that both metabolism and virulence genes provide advantages to endow XacA with higher virulence and a wider host range. Our analysis advances our understanding of the genomic basis of specialization by positive selection in bacterial evolution.

▪ Abstract This review deals with a comparative analysis of seven genome sequences from plant-associated bacteria. These are the genomes of Agrobacterium tumefaciens, Mesorhizobium loti, Sinorhizobium meliloti, Xanthomonas campestris pv campestris, Xanthomonas axonopodis pv citri, Xylella fastidiosa, and Ralstonia solanacearum. Genome structure and the metabolism pathways available highlight the compromise between the genome size and lifestyle. Despite the recognized importance of the type III secretion system in controlling host compatibility, its presence is not universal in all necrogenic pathogens. Hemolysins, hemagglutinins, and some adhesins, previously reported only for mammalian pathogens, are present in most organisms discussed. Different numbers and combinations of cell wall degrading enzymes and genes to overcome the oxidative burst generally induced by the plant host are characterized in these genomes. A total of 19 genes not involved in housekeeping functions were found common to all these bacteria.

We used PFGE followed by CsCl bisbenzymide centrifugation to obtain sufficient DNA for pyrosequencing of the Ca. Liberibacter americanus (Lam) strain “Sao Paulo” genome. The complete circular genomic DNA sequence of Lam is 1,195,201 bp, with an average GC content of 31.12%, somewhat lower than other Liberibacters. There are 1,056 predicted Lam genes, with 1,002 encoding proteins, 9 encoding rRNA genes and 45 encoding tRNAs. The overall gene organization and structure of the Lam genome is more similar to Lso than to Las. There are 951 genes common to Lam, Lso and Las, 27 genes found in Lam and Lso but not Las, and only 8 genes common to Lam and Las but not found in Lso. Many pseudogenes or truncated genes were found among the unique genes of all 3 species. As with Las, two prophage were confirmed in Lam, with SP2 being 39,941 bp and SP1 being 16,398 bp in size; as in Las, the one that appears to replicate as an excision plasmid prophage carries putative lysogenic conversion genes, specifically peroxidases and a Type Vc secreted adhesin. These predicted peroxidases and adhesin were found in both Las and Lam, but appeared fragmentary or degenerated in Lso, indicating their potential for citrus host range determination. Liberibacter genomes seem to be under selective pressure to reduce GC% content and to lose unimportant genes. Although Lam have an outer membrane, most of the genes required for biosynthesis of lipopolysaccharide, which can trigger ROS production, are missing from the Lam genome.

The 5.67-megabase genome of the plant pathogenAgrobacterium tumefaciensC58 consists of a circular chromosome, a linear chromosome, and two plasmids. Extensive orthology and nucleotide colinearity between the genomes ofA. tumefaciensand the plant symbiontSinorhizobium melilotisuggest a recent evolutionary divergence. Their similarities include metabolic, transport, and regulatory systems that promote survival in the highly competitive rhizosphere; differences are apparent in their genome structure and virulence gene complement. Availability of theA. tumefacienssequence will facilitate investigations into the molecular basis of pathogenesis and the evolutionary divergence of pathogenic and symbiotic lifestyles.

Instituto Ludwig de Pesquisa sobre o Câncer, Rua Prof. Antonio Prudente, 109-4 andar, 01509-010, Sao Paulo-SP

Background Citrus bacterial canker is a disease that has severe economic impact on citrus industries worldwide and is caused by a few species and pathotypes of Xanthomonas. X. citri subsp. citri strain 306 (XccA306) is a type A (Asiatic) strain with a wide host range, whereas its variant X. citri subsp. citri strain Aw12879 (Xcaw12879, Wellington strain) is restricted to Mexican lime. Results To characterize the mechanism for the differences in host range of XccA and Xcaw, the genome of Xcaw12879 that was completed recently was compared with XccA306 genome. Effectors xopAF and avrGf1 are present in Xcaw12879, but were absent in XccA306. AvrGf1 was shown previously for Xcaw to cause hypersensitive response in Duncan grapefruit. Mutation analysis of xopAF indicates that the gene contributes to Xcaw growth in Mexican lime but does not contribute to the limited host range of Xcaw. RNA-Seq analysis was conducted to compare the expression profiles of Xcaw12879 and XccA306 in Nutrient Broth (NB) medium and XVM2 medium, which induces hrp gene expression. Two hundred ninety two and 281 genes showed differential expression in XVM2 compared to in NB for XccA306 and Xcaw12879, respectively. Twenty-five type 3 secretion system genes were up-regulated in XVM2 for both XccA and Xcaw. Among the 4,370 common genes of Xcaw12879 compared to XccA306, 603 genes in NB and 450 genes in XVM2 conditions were differentially regulated. Xcaw12879 showed higher protease activity than XccA306 whereas Xcaw12879 showed lower pectate lyase activity in comparison to XccA306. Conclusions Comparative genomic analysis of XccA306 and Xcaw12879 identified strain specific genes. Our study indicated that AvrGf1 contributes to the host range limitation of Xcaw12879 whereas XopAF contributes to virulence. Transcriptome analyses of XccA306 and Xcaw12879 presented insights into the expression of the two closely related strains of X. citri subsp. citri. Virulence genes including genes encoding T3SS components and effectors are induced in XVM2 medium. Numerous genes with differential expression in Xcaw12879 and XccA306 were identified. This study provided the foundation to further characterize the mechanisms for virulence and host range of pathotypes of X. citri subsp. citri.

Two groups independently sequenced the Agrobacterium tumefaciens C58 genome in 2001. We report here consolidation of these sequences, updated annotation, and additional analysis of the evolutionary history of the linear chromosome, which is apparently limited to the biovar I group of Agrobacterium .

Brucella species are Gram-negative bacteria that infect mammals. Recently, two unusual strains (Brucella inopinata BO1T and B. inopinata-like BO2) have been isolated from human patients, and their similarity to some atypical brucellae isolated from Australian native rodent species was noted. Here we present a phylogenomic analysis of the draft genome sequences of BO1T and BO2 and of the Australian rodent strains 83-13 and NF2653 that shows that they form two groups well separated from the other sequenced Brucella spp. Several important differences were noted. Both BO1T and BO2 did not agglutinate significantly when live or inactivated cells were exposed to monospecific A and M antisera against O-side chain sugars composed of N-formyl-perosamine. While BO1T maintained the genes required to synthesize a typical Brucella O-antigen, BO2 lacked many of these genes but still produced a smooth LPS (lipopolysaccharide). Most missing genes were found in the wbk region involved in O-antigen synthesis in classic smooth Brucella spp. In their place, BO2 carries four genes that other bacteria use for making a rhamnose-based O-antigen. Electrophoretic, immunoblot, and chemical analyses showed that BO2 carries an antigenically different O-antigen made of repeating hexose-rich oligosaccharide units that made the LPS water-soluble, which contrasts with the homopolymeric O-antigen of other smooth brucellae that have a phenol-soluble LPS. The results demonstrate the existence of a group of early-diverging brucellae with traits that depart significantly from those of the Brucella species described thus far., IMPORTANCE This report examines differences between genomes from four new Brucella strains and those from the classic Brucella spp. Our results show that the four new strains are outliers with respect to the previously known Brucella strains and yet are part of the genus, forming two new clades. The analysis revealed important information about the evolution and survival mechanisms of Brucella species, helping reshape our knowledge of this important zoonotic pathogen. One discovery of special importance is that one of the strains, BO2, produces an O-antigen distinct from any that has been seen in any other Brucella isolates to date.

Inversions are one of the most frequent large-scale rearrangements observed in actual genomes. While a large body of literature exists on mathematical problems related to the computation of the inversion distance between abstract genomes, these works generally do not take into account that most inversions in bacterial chromosomes are symmetric or roughly symmetric with respect to the origin of replication. We define a new problem: how to sort genomes (or permutations) using almost-symmetric inversions. We show an algorithm that can sort any permutation using only almost-symmetric inversions. Two variants of this algorithm are presented that have better performance in practice. We explore the question of determining the minimum number of almost-symmetric inversions needed to sort a genome by presenting lower and upper bounds and results for special permutation families. The results obtained are the first steps in exploring this interesting new problem.

Background Decreasing costs of DNA sequencing have made prokaryotic draft genome sequences increasingly common. A contig scaffold is an ordering of contigs in the correct orientation. A scaffold can help genome comparisons and guide gap closure efforts. One popular technique for obtaining contig scaffolds is to map contigs onto a reference genome. However, rearrangements that may exist between the query and reference genomes may result in incorrect scaffolds, if these rearrangements are not taken into account. Large-scale inversions are common rearrangement events in prokaryotic genomes. Even in draft genomes it is possible to detect the presence of inversions given sufficient sequencing coverage and a sufficiently close reference genome. Results We present a linear-time algorithm that can generate a set of contig scaffolds for a draft genome sequence represented in contigs given a reference genome. The algorithm is aimed at prokaryotic genomes and relies on the presence of matching sequence patterns between the query and reference genomes that can be interpreted as the result of large-scale inversions; we call these patterns inversion signatures. Our algorithm is capable of correctly generating a scaffold if at least one member of every inversion signature pair is present in contigs and no inversion signatures have been overwritten in evolution. The algorithm is also capable of generating scaffolds in the presence of any kind of inversion, even though in this general case there is no guarantee that all scaffolds in the scaffold set will be correct. We compare the performance of sis, the program that implements the algorithm, to seven other scaffold-generating programs. The results of our tests show that sis has overall better performance. Conclusions sis is a new easy-to-use tool to generate contig scaffolds, available both as stand-alone and as a web server. The good performance of sis in our tests adds evidence that large-scale inversions are widespread in prokaryotic genomes.

Xanthomonas axonopodis pv. citrumelo is a citrus pathogen causing citrus bacterial spot disease that is geographically restricted within the state of Florida. Illumina, 454 sequencing, and optical mapping were used to obtain a complete genome sequence of X. axonopodis pv. citrumelo strain F1, 4.9 Mb in size. The strain lacks plasmids, in contrast to other citrus Xanthomonas pathogens. Phylogenetic analysis revealed that this pathogen is very close to the tomato bacterial spot pathogen X. campestris pv. vesicatoria 85-10, with a completely different host range. We also compared X. axonopodis pv. citrumelo to the genome of citrus canker pathogen X. axonopodis pv. citri 306. Comparative genomic analysis showed differences in several gene clusters, like those for type III effectors, the type IV secretion system, lipopolysaccharide synthesis, and others. In addition to pthA , effectors such as xopE3, xopAI , and hrpW were absent from X. axonopodis pv. citrumelo while present in X. axonopodis pv. citri. These effectors might be responsible for survival and the low virulence of this pathogen on citrus compared to that of X. axonopodis pv. citri. We also identified unique effectors in X. axonopodis pv. citrumelo that may be related to the different host range as compared to that of X. axonopodis pv. citri. X. axonopodis pv. citrumelo also lacks various genes, such as syrE1, syrE2 , and RTX toxin family genes, which were present in X. axonopodis pv. citri. These may be associated with the distinct virulences of X. axonopodis pv. citrumelo and X. axonopodis pv. citri. Comparison of the complete genome sequence of X. axonopodis pv. citrumelo to those of X. axonopodis pv. citri and X. campestris pv. vesicatoria provides valuable insights into the mechanism of bacterial virulence and host specificity.

We present a linear-time algorithm that can generate a contig scaffold for a draft genome sequence represented in contigs given a reference genome. The algorithm is aimed at prokaryotic genomes and relies on the presence of matching sequence patterns between the query and reference genomes that can be interpreted as the result of large-scale inversions; we call these patterns inversion signatures. Large-scale inversions are common rearrangement events in prokaryotic genomes. Even in draft genomes it is possible to detect the presence of inversions given sufficient sequencing coverage and a sufficiently close reference genome. Our algorithm is capable of correctly generating a scaffold if at least one member of every inversion signature pair is present in contigs and no inversion signatures have been overwritten in evolution. The algorithm is also capable of generating scaffolds in the presence of any kind of inversion, although in this general case there is no guarantee that the scaffold will be completely correct. We compare the performance of SIS, the program that implements the algorithm, to five other scaffold-generating programs. The results from two batches of tests using real genomes and artificial contig boundaries show that SIS has significantly better performance.

With the increase of the number of available prokaryotic genomes, studies focusing on the higher level evolutionary events, including the evolution of whole genomes, are now possible. Simulated datasets are essential components in both the development and evaluation stages for all evolutionary studies. In this paper, we implemented the first random events-based high level bacterial genome evolution simulator with models for a large number of possible events at both the gene- and replicon-scale. We call this simulator Prokaryote Evolutionary Genomics Simulator (PEGsim). PEGsim also implements the concept of conserved blocks, which can be biologically related to both Operons and/or regulons. We present results showing genome sizes and replicon numbers in simulated data and in real data. These results demonstrate that PEGsim is capable of generating evolutionary-related bacterial genome datasets with diverse realistic properties.

Background Bacterial spot of tomato and pepper is caused by four Xanthomonas species and is a major plant disease in warm humid climates. The four species are distinct from each other based on physiological and molecular characteristics. The genome sequence of strain 85-10, a member of one of the species, Xanthomonas euvesicatoria (Xcv) has been previously reported. To determine the relationship of the four species at the genome level and to investigate the molecular basis of their virulence and differing host ranges, draft genomic sequences of members of the other three species were determined and compared to strain 85-10. Results We sequenced the genomes of X. vesicatoria (Xv) strain 1111 (ATCC 35937), X. perforans (Xp) strain 91-118 and X. gardneri (Xg) strain 101 (ATCC 19865). The genomes were compared with each other and with the previously sequenced Xcv strain 85-10. In addition, the molecular features were predicted that may be required for pathogenicity including the type III secretion apparatus, type III effectors, other secretion systems, quorum sensing systems, adhesins, extracellular polysaccharide, and lipopolysaccharide determinants. Several novel type III effectors from Xg strain 101 and Xv strain 1111 genomes were computationally identified and their translocation was validated using a reporter gene assay. A homolog to Ax21, the elicitor of XA21-mediated resistance in rice, and a functional Ax21 sulfation system were identified in Xcv. Genes encoding proteins with functions mediated by type II and type IV secretion systems have also been compared, including enzymes involved in cell wall deconstruction, as contributors to pathogenicity. Conclusions Comparative genomic analyses revealed considerable diversity among bacterial spot pathogens, providing new insights into differences and similarities that may explain the diverse nature of these strains. Genes specific to pepper pathogens, such as the O-antigen of the lipopolysaccharide cluster, and genes unique to individual strains, such as novel type III effectors and bacteriocin genes, have been identified providing new clues for our understanding of pathogen virulence, aggressiveness, and host preference. These analyses will aid in efforts towards breeding for broad and durable resistance in economically important tomato and pepper cultivars.

Agrobacterium sp. H13-3, formerly known as Rhizobium lupini H13-3, is a soil bacterium that was isolated from the rhizosphere of Lupinus luteus. The isolate has been established as a model system for studying novel features of flagellum structure, motility and chemotaxis within the family Rhizobiaceae. The complete genome sequence of Agrobacterium sp. H13-3 has been established and the genome structure and phylogenetic assignment of the organism was analysed. For de novo sequencing of the Agrobacterium sp. H13-3 genome, a combined strategy comprising 454-pyrosequencing on the Genome Sequencer FLX platform and PCR-based amplicon sequencing for gap closure was applied. The finished genome consists of three replicons and comprises 5,573,770 bases. Based on phylogenetic analyses, the isolate could be assigned to the genus Agrobacterium biovar I and represents a genomic species G1 strain within this biovariety. The highly conserved circular chromosome (2.82Mb) of Agrobacterium sp. H13-3 mainly encodes housekeeping functions characteristic for an aerobic, heterotrophic bacterium. Agrobacterium sp. H13-3 is a motile bacterium driven by the rotation of several complex flagella. Its behaviour towards external stimuli is regulated by a large chemotaxis regulon and a total of 17 chemoreceptors. Comparable to the genome of Agrobacterium tumefaciens C58, Agrobacterium sp. H13-3 possesses a linear chromosome (2.15Mb) that is related to its reference replicon and features chromosomal and plasmid-like properties. The accessory plasmid pAspH13-3a (0.6Mb) is only distantly related to the plasmid pAtC58 of A. tumefaciens C58 and shows a mosaic structure. A tumor-inducing Ti-plasmid is missing in the sequenced strain H13-3 indicating that it is a non-virulent isolate. Copyright 2011 Elsevier B.V. All rights reserved.