A Method for Improving the Accuracy and Efficiency of Bacteriophage Genome Annotation

Bacteriophages are the most numerous entities on Earth. The number of sequenced phage genomes is approximately 8000 and increasing rapidly. Sequencing of a genome is followed by annotation, where genes, start codons, and functions are putatively identified. The mainstays of phage genome annotation are auto-annotation programs such as Glimmer and GeneMark. Due to the relatively small size of phage genomes, many groups choose to manually curate auto-annotation results to increase accuracy. An additional benefit of manual curation of auto-annotated phage genomes is that the process is amenable to be performed by students, and has been shown to improve student recruitment to the sciences. However, despite its greater accuracy and pedagogical value, manual curation suffers from high labor cost, lack of standardization and a degree of subjectivity in decision making, and susceptibility to mistakes. Here, we present a method developed in our lab that is designed to produce accurate annotations while reducing subjectivity and providing a degree of standardization in decision-making. We show that our method produces genome annotations more accurate than auto-annotation programs while retaining the pedagogical benefits of manual genome curation

Comparative Genomics of 9 Novel Paenibacillus Larvae Bacteriophages

American Foulbrood Disease, caused by the bacterium Paenibacillus larvae, is one of the most destructive diseases of the honeybee, Apis mellifera. Our group recently published the sequences of 9 new phages with the ability to infect and lyse P. larvae. Here, we characterize the genomes of these P. larvae phages, compare them to each other and to other sequenced P. larvae phages, and putatively identify protein function. The phage genomes are 38–45 kb in size and contain 68–86 genes, most of which appear to be unique to P. larvae phages. We classify P. larvae phages into 2 main clusters and one singleton based on nucleotide sequence identity. Three of the new phages show sequence similarity to other sequenced P. larvae phages, while the remaining 6 do not. We identified functions for roughly half of the P. larvae phage proteins, including structural, assembly, host lysis, DNA replication/metabolism, regulatory, and host-related functions. Structural and assembly proteins are highly conserved among our phages and are located at the start of the genome. DNA replication/metabolism, regulatory, and host-related proteins are located in the middle and end of the genome, and are not conserved, with many of these genes found in some of our phages but not others. All nine phages code for a conserved N-acetylmuramoyl-L-alanine amidase. Comparative analysis showed the phages use the “cohesive ends with 30 overhang” DNA packaging strategy. This work is the first in-depth study of P. larvae phage genomics, and serves as a marker for future work in this area

Mechanisms of B cell Synapse Formation Predicted by Stochastic Simulation

The clustering of B cell receptor (BCR) molecules and the formation of the protein segregation structure known as the immunological synapse appears to precede antigen (Ag) uptake by B cells. The mature B cell synapse is characterized by a central cluster of BCR/Ag molecular complexes surrounded by a ring of LFA-1/ICAM-1 complexes. Recent experimental evidence shows receptor clustering in B cells can occur via mechanical or signaling-driven processes. An alternative mechanism of diffusion and affinity-dependent binding has been proposed to explain synapse formation in the absence of signaling-driven processes. In this work, we investigated the biophysical mechanisms that drive immunological synapse formation in B cells across the physiological range of BCR affinity (KA~10^6-10^10 M-1) through computational modeling. Our computational approach is based on stochastic simulation of diffusion and reaction events with a clearly defined mapping between probabilistic parameters of our model and their physical equivalents. We show that a diffusion-and-binding mechanism is sufficient to drive synapse formation only at low BCR affinity and for a relatively stiff B cell membrane that undergoes little deformation. We thus predict the need for alternative mechanisms: a difference in the mechanical properties of BCR/Ag and LFA-1/ICAM-1 bonds and/or signaling driven processes.Comment: 35 pages, 11 figures; Supplemental Materials adde

Characterization of CRISPR Spacer and Protospacer Sequences in Paenibacillus larvae and Its Bacteriophages

The bacterium Paenibacillus larvae is the causative agent of American foulbrood, the most devastating bacterial disease of honeybees. Because P. larvae is antibiotic resistant, phages that infect it are currently used as alternative treatments. However, the acquisition by P. larvae of CRISPR spacer sequences from the phages could be an obstacle to treatment efforts. We searched nine complete genomes of P. larvae strains and identified 714 CRISPR spacer sequences, of which 384 are unique. Of the four epidemiologically important P. larvae strains, three of these have fewer than 20 spacers, while one strain has over 150 spacers. Of the 384 unique spacers, 18 are found as protospacers in the genomes of 49 currently sequenced P. larvae phages. One P. larvae strain does not have any protospacers found in phages, while another has eight. Protospacer distribution in the phages is uneven, with two phages having up to four protospacers, while a third of phages have none. Some phages lack protospacers found in closely related phages due to point mutations, indicating a possible escape mechanism. This study serve a point of reference for future studies on the CRISPR-Cas system in P. larvae as well as for comparative studies of other phage-host systems

Complete Genome Sequences of Paenibacillus Larvae Phages BN12, Dragolir, Kiel007, Leyra, Likha, Pagassa, PBL1c, and Tadhana

We present here the complete genomes of eight phages that infect Paenibacillus larvae, the causative agent of American foulbrood in honeybees. Phage PBL1c was originally isolated in 1984 from a P. larvae lysogen, while the remaining phages were isolated in 2014 from bee debris, honeycomb, and lysogens from three states in the USA

Complete Genome Sequences of Mycobacterium smegmatis Phages Chewbacca, Reptar3000, and Riparian, Isolated in Las Vegas, Nevada

Here, we present the complete genome sequences of Mycobacterium smegmatis phages Chewbacca, Reptar3000, and Riparian, isolated from soil in Las Vegas, NV. The phages were isolated and annotated by undergraduate students enrolled in the Phage Discovery course offered by the School of Life Sciences at the University of Nevada, Las Vega

Monte Carlo Investigation of Diffusion of Receptors and Ligands that Bind Across Opposing Surfaces

Studies of receptor diffusion on a cell surface show a variety of behaviors, such as diffusive, sub-diffusive, or super-diffusive motion. However, most studies to date focus on receptor molecules diffusing on a single cell surface. We have previously studied receptor diffusion to probe the molecular mechanism of receptor clustering at the cell–cell junction between two opposing cell surfaces. Here, we characterize the diffusion of receptors and ligands that bind to each other across two opposing cell surfaces, as in cell–cell and cell–bilayer interactions. We use a Monte Carlo method, where receptors and ligands are simulated as independent agents that bind and diffuse probabilistically. We vary receptor–ligand binding affinity and plot the molecule-averaged mean square displacement (MSD) of ligand molecules as a function of time. Our results show that MSD plots are qualitatively different for flat and curved interfaces, as well as between the cases of presence and absence of directed transport of receptor–ligand complexes toward a specific location on the interface. Receptor–ligand binding across two opposing surfaces leads to transient sub-diffusive motion at early times provided the interface is flat. This effect is entirely absent if the interface is curved, however, in this instance we observe sub-diffusive motion. In addition, a decrease in the equilibrium value of the MSD occurs as affinity increases, something which is absent for a flat interface. In the presence of directed transport of receptor–ligand complexes, we observe super-diffusive motion at early times for a flat interface. Super-diffusive motion is absent for a curved interface, however, in this case we observe a transient decrease in MSD with time prior to equilibration for high-affinity values

Modeling of B cell Synapse Formation by Monte Carlo Simulation Shows That Directed Transport of Receptor Molecules Is a Potential Formation Mechanism

The formation of the protein segregation structure known as the “immunological synapse” in the contact region between B cells and antigen presenting cells appears to precede antigen (Ag) uptake by B cells. The mature B cell synapse consists of a central cluster of B cell receptor/Antigen (BCR/Ag) complexes surrounded by a ring of LFA-1/ICAM-1 complexes. In this study, we used an in silico model to investigate whether cytoskeletally driven transport of molecules toward the center of the contact zone is a potential mechanism of immunological synapse formation in B cells. We modeled directed transport by the cytoskeleton in an effective manner, by biasing the diffusion of molecules toward the center of the contact zone. Our results clearly show that biased diffusion of BCR/Ag complexes on the B cell surface is sufficient to produce patterns similar to experimentally observed immunological synapses. This is true even in the presence of significant membrane deformation as a result of receptor–ligand binding, which in previous work we showed had a detrimental effect on synapse formation at high antigen affinity values. Comparison of our model’s results to those of experiments shows that our model produces synapses for realistic length, time, and affinity scales. Our results also show that strong biased diffusion of free molecules has a negative effect on synapse formation by excluding BCR/Ag complexes from the center of the contact zone. However, synapses may still form provided the bias in diffusion of free molecules is an order-of-magnitude weaker than that of BCR/Ag complexes. We also show how diffusion trajectories obtained from single-molecule tracking experiments can generate insight into the mechanism of synapse formation