Motivation: Completing the genome sequence of an organism is an important task in comparative, functional and structural genomics. However, this remains a challenging issue from both a computational and an experimental viewpoint. Genome scaffolding (i.e. the process of ordering and orientating contigs) of de novo assemblies usually represents the first step in most genome finishing pipelines. Results: In this article we present MeDuSa (Multi-Draft based Scaffolder), an algorithm for genome scaffolding. MeDuSa exploits information obtained from a set of (draft or closed) genomes from related organisms to determine the correct order and orientation of the contigs. MeDuSa formalizes the scaffolding problem by means of a combinatorial optimization formulation on graphs and implements an efficient constant factor approximation algorithm to solve it. In contrast to currently used scaffolders, it does not require either prior knowledge on the microrganisms dataset under analysis (e.g. their phylogenetic relationships) or the availability of paired end read libraries. This makes usability and running time two additional important features of our method. Moreover, benchmarks and tests on real bacterial datasets showed that MeDuSa is highly accurate and, in most cases, outperforms traditional scaffolders. The possibility to use MeDuSa on eukaryotic datasets has also been evaluated, leading to interesting results. Availability and implementation: MeDuSa web server: http://combo.dbe.unifi.it/medusa . A stand-alone version of the software can be downloaded from https://github.com/combogenomics/medusa/releases. All results presented in this work have been obtained with MeDuSa v. 1.3. Contact:marco.fondi@unifi.it Supplementary information: Supplementary data are available at Bioinformatics online.
10 June 2016 in Proceedings of the National Academy of Sciences
Significance: Comparative analysis of multiple strains within a species is a powerful way to uncover pathoadaptive genetic acquisitions. Hundreds of genome sequences are now available for the human pathogen Staphylococcus aureus, mostly known for its antibiotic-resistant variants that threaten the emergence of panresistant superbugs. In this study, genome-scale models of metabolism are used to analyze the shared and unique metabolic capabilities of this pathogen and its strain-specific variants. The models are used to distinguish S. aureus strains responsible for severe infections based solely on growth capabilities and presence of different virulence factors. The results identify metabolic similarities and differences between S. aureus strains that provide insights into the epidemiology of S. aureus and may help to combat its spread.
Background: In the bacterium Caulobacter crescentus, CtrA coordinates DNA replication, cell division, and polar morphogenesis and is considered the cell cycle master regulator. CtrA activity varies during cell cycle progression and is modulated by phosphorylation, proteolysis and transcriptional control. In a phosphorylated state, CtrA binds specific DNA sequences, regulates the expression of genes involved in cell cycle progression and silences the origin of replication. Although the circuitry regulating CtrA is known in molecular detail in Caulobacter, its conservation and functionality in the other alpha-proteobacteria are still poorly understood. Results: Orthologs of Caulobacter factors involved in the regulation of CtrA were systematically scanned in genomes of alpha-proteobacteria. In particular, orthologous genes of the divL-cckA-chpT-ctrA phosphorelay, the divJ-pleC-divK two-component system, the cpdR-rcdA-clpPX proteolysis system, the methyltransferase ccrM and transcriptional regulators dnaA and gcrA were identified in representative genomes of alpha-proteobacteria. CtrA, DnaA and GcrA binding sites and CcrM putative methylation sites were predicted in promoter regions of all these factors and functions controlled by CtrA in all alphas were predicted. Conclusions: The regulatory cell cycle architecture was identified in all representative alpha-proteobacteria, revealing a high diversification of circuits but also a conservation of logical features. An evolutionary model was proposed where ancient alphas already possessed all modules found in Caulobacter arranged in a variety of connections. Two schemes appeared to evolve: a complex circuit in Caulobacterales and Rhizobiales and a simpler one found in Rhodobacterales.
Background: The pioneering ancestor of land plants that conquered terrestrial habitats around 500 million years ago had to face dramatic stresses including UV radiation, desiccation, and microbial attack. This drove a number of adaptations, among which the emergence of the phenylpropanoid pathway was crucial, leading to essential compounds such as flavonoids and lignin. However, the origin of this specific land plant secondary metabolism has not been clarified.Results: We have performed an extensive analysis of the taxonomic distribution and phylogeny of Phenylalanine Ammonia Lyase (PAL), which catalyses the first and essential step of the general phenylpropanoid pathway, leading from phenylalanine to p-Coumaric acid and p-Coumaroyl-CoA, the entry points of the flavonoids and lignin routes. We obtained robust evidence that the ancestor of land plants acquired a PALviahorizontal gene transfer (HGT) during symbioses with soil bacteria and fungi that are known to have established very early during the first steps of land colonization. This horizontally acquired PAL represented then the basis for further development of the phenylpropanoid pathway and plant radiation on terrestrial environments.Conclusion: Our results highlight a possible crucial role of HGT from soil bacteria in the path leading to land colonization by plants and their subsequent evolution. The few functional characterizations of sediment/soil bacterial PAL (production of secondary metabolites with powerful antimicrobial activity or production of pigments) suggest that the initial advantage of this horizontally acquired PAL in the ancestor of land plants might have been either defense against an already developed microbial community and/or protection against UV.Reviewers: This article was reviewed by Purificación López-García, Janet Siefert, and Eugene Koonin.