Elasticity Theory and Shape Transitions of Viral Shells

Recently, continuum elasticity theory has been applied to explain the shape transition of icosahedral viral capsids - single-protein-thick crystalline shells - from spherical to buckled/faceted as their radius increases through a critical value determined by the competition between stretching and bending energies of a closed 2D elastic network. In the present work we generalize this approach to capsids with non-icosahedral symmetries, e.g., spherocylindrical and conical shells. One key new physical ingredient is the role played by nonzero spontaneous curvature. Another is associated with the special way in which the energy of the twelve topologically-required five-fold sites depends on the background local curvature of the shell in which they are embedded. Systematic evaluation of these contributions leads to a shape phase diagram in which transitions are observed from icosahedral to spherocylindrical capsids as a function of the ratio of stretching to bending energies and of the spontaneous curvature of the 2D protein network. We find that the transition from icosahedral to spherocylindrical symmetry is continuous or weakly first-order near the onset of buckling, leading to extensive shape degeneracy. These results are discussed in the context of experimentally observed variations in the shapes of a variety of viral capsids.Comment: 53 pages, 17 figure

Electrostatics and the Assembly of an RNA Virus

Electrostatic interactions play a central role in the assembly of single-stranded RNA viruses. Under physiological conditions of salinity and acidity, virus capsid assembly requires the presence of genomic material that is oppositely charged to the core proteins. In this paper we apply basic polymer physics and statistical mechanics methods to the self-assembly of a synthetic virus encapsidating generic polyelectrolyte molecules. We find that (i) the mean concentration of the encapsidated polyelectrolyte material depends on the surface charge density, the radius of the capsid, and the linear charge density of the polymer but neither on the salt concentration or the Kuhn length, (ii) the total charge of the capsid interior is equal but opposite to that of the empty capsid, a form of charge reversal. Unlike natural viruses, synthetic viruses are predicted not to be under an osmotic swelling pressure. The design condition that self-assembly only produces filled capsids is shown to coincide with the condition that the capsid surface charge exceeds the desorption threshold of polymer surface adsorption. We compare our results with studies on the self-assembly of both synthetic and natural viruses.Comment: 41 pages, 4 figure

PHAST: A Fast Phage Search Tool

PHAge Search Tool (PHAST) is a web server designed to rapidly and accurately identify, annotate and graphically display prophage sequences within bacterial genomes or plasmids. It accepts either raw DNA sequence data or partially annotated GenBank formatted data and rapidly performs a number of database comparisons as well as phage ‘cornerstone’ feature identification steps to locate, annotate and display prophage sequences and prophage features. Relative to other prophage identification tools, PHAST is up to 40 times faster and up to 15% more sensitive. It is also able to process and annotate both raw DNA sequence data and Genbank files, provide richly annotated tables on prophage features and prophage ‘quality’ and distinguish between intact and incomplete prophage. PHAST also generates downloadable, high quality, interactive graphics that display all identified prophage components in both circular and linear genomic views. PHAST is available at (http://phast.wishartlab.com)

Genome Stability of Lyme Disease Spirochetes: Comparative Genomics of Borrelia burgdorferi Plasmids

Lyme disease is the most common tick-borne human illness in North America. In order to understand the molecular pathogenesis, natural diversity, population structure and epizootic spread of the North American Lyme agent, Borrelia burgdorferi sensu stricto, a much better understanding of the natural diversity of its genome will be required. Towards this end we present a comparative analysis of the nucleotide sequences of the numerous plasmids of B. burgdorferi isolates B31, N40, JD1 and 297. These strains were chosen because they include the three most commonly studied laboratory strains, and because they represent different major genetic lineages and so are informative regarding the genetic diversity and evolution of this organism. A unique feature of Borrelia genomes is that they carry a large number of linear and circular plasmids, and this work shows that strains N40, JD1, 297 and B31 carry related but non-identical sets of 16, 20, 19 and 21 plasmids, respectively, that comprise 33–40% of their genomes. We deduce that there are at least 28 plasmid compatibility types among the four strains. The B. burgdorferi ∼900 Kbp linear chromosomes are evolutionarily exceptionally stable, except for a short ≤20 Kbp plasmid-like section at the right end. A few of the plasmids, including the linear lp54 and circular cp26, are also very stable. We show here that the other plasmids, especially the linear ones, are considerably more variable. Nearly all of the linear plasmids have undergone one or more substantial inter-plasmid rearrangements since their last common ancestor. In spite of these rearrangements and differences in plasmid contents, the overall gene complement of the different isolates has remained relatively constant

Borrelia burgdorferi EbfC defines a newly-identified, widespread family of bacterial DNA-binding proteins

The Lyme disease spirochete, Borrelia burgdorferi, encodes a novel type of DNA-binding protein named EbfC. Orthologs of EbfC are encoded by a wide range of bacterial species, so characterization of the borrelial protein has implications that span the eubacterial kingdom. The present work defines the DNA sequence required for high-affinity binding by EbfC to be the 4 bp broken palindrome GTnAC, where ‘n’ can be any nucleotide. Two high-affinity EbfC-binding sites are located immediately 5′ of B. burgdorferi erp transcriptional promoters, and binding of EbfC was found to alter the conformation of erp promoter DNA. Consensus EbfC-binding sites are abundantly distributed throughout the B. burgdorferi genome, occurring approximately once every 1 kb. These and other features of EbfC suggest that this small protein and its orthologs may represent a distinctive type of bacterial nucleoid-associated protein. EbfC was shown to bind DNA as a homodimer, and site-directed mutagenesis studies indicated that EbfC and its orthologs appear to bind DNA via a novel α-helical ‘tweezer’-like structure

Population Genomics and Phylogeography of an Australian Dairy Factory Derived Lytic Bacteriophage

In this study, we present the full genomic sequences and evolutionary analyses of a serially sampled population of 28 Lactococcus lactis–infecting phage belonging to the 936-like group in Australia. Genome sizes were consistent with previously available genomes ranging in length from 30.9 to 32.1 Kbp and consisted of 55–65 open reading frames. We analyzed their genetic diversity and found that regions of high diversity are correlated with high recombination rate regions (P value = 0.01). Phylogenetic inference showed two major clades that correlate well with known host range. Using the extended Bayesian Skyline model, we found that population size has remained mostly constant through time. Moreover, the dispersion pattern of these genomes is in agreement with human-driven dispersion as suggested by phylogeographic analysis. In addition, selection analysis found evidence of positive selection on codon positions of the Receptor Binding Protein (RBP). Likewise, positively selected sites in the RBP were located within the neck and head region in the crystal structure, both known determinants of host range. Our study demonstrates the utility of phylogenetic methods applied to whole genome data collected from populations of phage for providing insights into applied microbiology

P-value based visualization of codon usage data

Two important and not yet solved problems in bacterial genome research are the identification of horizontally transferred genes and the prediction of gene expression levels. Both problems can be addressed by multivariate analysis of codon usage data. In particular dimensionality reduction methods for visualization of multivariate data have shown to be effective tools for codon usage analysis. We here propose a multidimensional scaling approach using a novel similarity measure for codon usage tables. Our probabilistic similarity measure is based on P-values derived from the well-known chi-square test for comparison of two distributions. Experimental results on four microbial genomes indicate that the new method is well-suited for the analysis of horizontal gene transfer and translational selection. As compared with the widely-used correspondence analysis, our method did not suffer from outlier sensitivity and showed a better clustering of putative alien genes in most cases

Data for Millennia of genomic stability within the invasive Para C Lineage of Salmonella enterica: date estimation 1

Salmonella enterica serovar Paratyphi C is the causative agent of enteric (paratyphoid) fever. While today a potentially lethal infection of humans that occurs in Africa and Asia, early 20th century observations in Eastern Europe suggest it may once have had a wider-ranging impact on human societies. We recovered a draft Paratyphi C genome from the 800-year-old skeleton of a young woman in Trondheim, Norway, who likely died of enteric fever. Analysis of this genome against a new, significantly expanded database of related modern genomes demonstrated that Paratyphi C is descended from the ancestors of swine pathogens, serovars Choleraesuis and Typhisuis, together forming the Para C Lineage. Our results indicate that Paratyphi C has been a pathogen of humans for at least 1,000 years, and may have evolved after zoonotic transfer from swine during the Neolithic period